Modelo Hidrologico Swat
Modelo Hidrologico Swat
Modelo Hidrologico Swat
t
=
0
+(
du
-
su]
-
u
-
sccp
-
qw
)
t
=1
donde SWt es el contenido final de agua en el suelo (mmH2O); SW0 es el contenido de agua
inicial en el da i(mmH2O); t es el tiempo (da); Rday es la cantidad de precipitacin en el da i
(mmH2O); wseep es la cantidad de agua acumulada en la zona no saturada del perfil del suelo
en da i (mmH2O); Qgw es la cantidad de flujo de retorno en el da i (mmH2O); (Arnold
etal.,1991).
Los principales componentes del modelo incluyen: el clima, hidrologa, sedimentacin,
temperatura del suelo, crecimiento de cultivos, nutrientes, pesticidas y manejo de cultivos.
En SWAT, una cuenca se divide en subcuencas mltiples, que estn divididos en unidades de
respuesta hidrolgica (URHs) que consisten en reas homogneas de uso de la tierra y
caractersticas del suelo, las cuales son identificadas dentro de un espacio de simulacin de
SWAT.
2
Para ms informacin visitar la pgina web: http://swatmodel.tamu.edu/
10
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
I. FUENTES DE INFORMACIN
A partir de la informacin que fue suministra por el Ministerio de Agricultura, el Servicio
Nacional de Meteorologa e Hidrologa del Per (SENAMHI), la reserva Piasajistica Nor-
Yauyos-Cocha y la Junta de Regantes, inicialmente se hizo una detallada y minuciosa
bsqueda con el fin de extraer informacin valiosa para las entradas de los modelos y
establecer un primer panorama de la cuenca de inters. En la siguiente tabla se presenta las
fuentes de la informacin cartogrfica seleccionada:
TABLA 1. FUENTES DE LA INFORMACIN SELECCIONADA INICIALMENTE
*Proveido por la Reserva Paisajistica Nor Yauyos-Cochas
Uso del Suelo/Cobertura: Se utilizaron los datos de las Cartas Nacionales para las
provincias de Junn y Lima, realizado por el IGN (Instituto Geogrfico Nacional) en los aos
1970 y 1971 a escala 1:100,000. En este estudio se emplearon mtodos
estereofotogrametricos de fotografas areas tomadas en 1962 con la elaboracin de un
control horizontal y vertical realizado por el IGM y una clasificacin de campo en el ao
1969.
En la Tabla 2 se presentan las Coberturas Vegetales de este estudio y a su vez en la Figura 5
se puede observar la distribucin de estas coberturas dentro de la cuenca.
De acuerdo al estudio Inventario y evaluacin de los recursos naturales de la reserva
paisajstica NorYauyos Cochas en la cuenca alta del Rio Caete, se describen unas
coberturas que fue objeto del estudio.
CUENCA RIO CAETE
Cobertura
Coberturas Generales a Escala 1:100.000, Cartas Nacionales del Instituto Geografico
Nacional (IGN) - 1970
Suelos
Suelos Generales a Escala 1:100.000, Oficina Nacional de Evaluacion de Recursos
Natulales (ONERN) y Corporacion Departamental de Desarrollo de Lima (CORLIMA) - 1989
DEM DEM SRTM con Resolucin de 90 m, CIAT
Cartografa
Base
Drenajes, Centros poblados, Divisin Poltico Administrativa, Nevados, etc.
*
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
TABLA 2. COBERTURA VEGETAL
Fuente: Cartas Nacionales, IGN 1971.
FIGURA 2. MAPA DE COBERTURA
COBERTURA CAETE
Matorrales
Planicies costeras y estribaciones andinas sin vegetacion
Tierras altoandinas sin vegetacin
Herbazal de tundra
Lagos y Lagunas
Pajonal / Cesped de puna
Bofedal
Queoal
Pajonal de puna
Cultivos agropecuarios
Nevados
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
. COBERTURA VEGETAL DE LA CUENCA DEL RIO CAETE
. MAPA DE COBERTURA VEGETAL DE LA CUENCA DEL RIO CAETE.
DESCRIPCION
Comunidades arbustivas de ambientes secos
Planicies costeras y estribaciones andinas sin vegetacion Desierto costero sin vegetacin
Areas altoandinas sin vegetacin (influencia de glaciares)
Herbaceas de poco desarrollo, sobre los 4500 msnm (Pastoreo)
Cuerpos de agua
Agrupa los tipos de asociaciones Pajonal de puna + Cesped de puna
Herbceas de piso, pegadas a ras del suelo.
Bosque altoandino, sobre los 3800 msnm.
Herbceas altoandinas hasta 1 m de alto con hojas punzo-cortantes
Areas con cultivos agrcolas + Areas con pastizales cultivados
Nevados
11
36626
.cgiar.org
SIMBOLO AREA (Ha) AREA (%)
Ma 188600 31.3
Pl ce Sv 78730 13.1
Al Sv 46290 7.7
Ht 157100 26.0
Lag 3225 0.5
Pj/Cp 109100 18.1
Bo 1271 0.2
Q 7095 1.2
Pj pu 174 0.0
Cuap 6785 1.1
Nv 4758 0.8
TOTAL 603128 100
12
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
Bosque Queoa (Bq)
Este bosque est representado por el gnero Polylepis queua, oriundo de Sudamrica; es
un rbol achaparrado, con tronco retorcido y abundancia de ritidomas exfoliables de
consistencia papirceo de color pardo rojizo; para los andes peruanos se han reportado ms
de diez especies, y mayormente se les utiliza para la formacin de barreras vivas por encima
de los 3500 msnm. La copa es generalmente difusa, presenta una madera dura, pesada,
resistente a la humedad y laborable. Frecuentemente utilizada para lea, es una de las
principales fuentes de energa para las comunidades altoandinas. Dentro de sus
propiedades fisicoqumicas resalta su alta densidad, mediana resistencia a la compresin y
alta resistencia al cizallamiento. Las hojas y ramas trituradas y hervidas son utilizadas como
colorante marrn claro. (Gil Mora et al, 2007).
Pajonal de Puna / Csped de puna (Pj/Cp)
Esta gran unidad de vegetacin conformada por comunidades de hierbas se encuentra
ampliamente distribuida, extendindose en una superficie de 142 409 ha, que representa el
64.28% del rea estudiada. Se encuentra ubicada en la porcin superior de la cuenca del ro
Caete, desde los 3800 hasta aproximadamente 3900 msnm, ocupando las laderas y cimas
de paisajes colinosos y montaosos. Comprende las zonas de vida pramo muy hmedo
Subalpino Tropical y tundra pluvial Altino Tropical.
Bofedal (Bo)
Esta unidad de vegetacin se encuentra ubicada en reas depresionadas con mal drenaje,
tanto en planicies como en superficies inclinadas. Se encuentra a menudo en los alrededores
de las lagunas y cochas. Est conformada por comunidades de hierbas de porte
almohadillado o en cojn, conocida como turbera, la cual permanece siempre verdes
durante el ao gracias a la humedad permanente.
Informacin de Suelos: Esta informacin fue suministrada por el Ministerio de Agricultura
y corresponde al Mapa de Suelos y Capacidad de Uso Mayor realizado por la Oficina
Nacional de Evaluacin de Recursos Naturales (ONERN) y la Corporacin Departamental de
Desarrollo de Lima (CORLIMA) en el ao 1989 a escala 1:100.000.
Para obtener el mapa final de Suelos de la cuenca del Rio Caete fue necesario realizar un
trabajo previo de georeferenciacion y digitalizacin de dos planchas cartogrficas en
formato anlogo que conforman la regin de estudio; actividades adicionales que fueron
necesarias para obtener el mapa de suelos en la mejor escala de trabajo posible.
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
FIGURA 3. DIGITALIZACIN MAP
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
. DIGITALIZACIN MAPA DE SUELOS DE LA CUENCA DEL RIO CAETE
13
36626
.cgiar.org
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
FIGURA 4. MAPA DE SUELOS DE LA CUENCA
Modelo Digital de Elevacin
digital del terreno se utilizo DEM SRTM (Shutt
espacial 90 aportado por CIAT
han realizado los ajustes necesarios en el llenado de vacos de informacin. La cuenca del
Ro Caete presenta una gran v
se puede apreciar en la figura
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
LA CUENCA DEL RIO CAETE
de Elevacin: Para calcular la red hdrica de la cuenca, se utilizo el modelo
digital del terreno se utilizo DEM SRTM (Shuttle Radar Topography Mission) de resolucin
espacial 90 aportado por CIAT (Centro Internacional de Agricultura Tropical), al cual se le
han realizado los ajustes necesarios en el llenado de vacos de informacin. La cuenca del
presenta una gran variedad en altitudes que van desde los 3 a 5861
se puede apreciar en la figura 5.
14
36626
.cgiar.org
Para calcular la red hdrica de la cuenca, se utilizo el modelo
le Radar Topography Mission) de resolucin
, al cual se le
han realizado los ajustes necesarios en el llenado de vacos de informacin. La cuenca del
msnm como
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
FIGURA 5. MODELO DIGITAL DE
Informacin Hidrometeorolgica
como Humedad Relativa, Precipitacin, Temperatura (Mximas, Mnima
Radiacin Solar y Aforos de Caudal
Hidrologa (SENAMHI). La informacin se analiz con el fin de est
podran tenerse en cuenta para la modelacin, cumpliendo con registros
para el Perodo Simulado (1991
inters.
En la siguiente tabla se presenta las es
Modelo:
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
IGITAL DE ELEVACIN DE LA CUENCA DEL RIO CAETE
lgica: Se cont con datos de Estaciones Hidrometeorolgicas
como Humedad Relativa, Precipitacin, Temperatura (Mximas, Mnimas y Promedios),
Aforos de Caudal, suministrada por el Servicio Nacional de Meteorologa e
informacin se analiz con el fin de establecer las estaciones que
podran tenerse en cuenta para la modelacin, cumpliendo con registros diarios
1991-2009) y que se encontraran dentro o cerca de
En la siguiente tabla se presenta las estaciones que se establecieron como entradas al
15
36626
.cgiar.org
rometeorolgicas
s y Promedios),
, suministrada por el Servicio Nacional de Meteorologa e
ablecer las estaciones que
diarios y continuos
cerca de la cuenca de
e establecieron como entradas al
16
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
TABLA 3. ESTACIONES CLIMATOLGICAS UTILIZADA EN EL MODELO
NOMBRE ESTE (m) NORTE (m)
ALTITUD
(m)
PRECIPITACION
DIARIA
TEMPERATURA
DIARIA
CLIMATOLOGICA
MENSUAL
CAETE 355705.042 8554188.667 158 X X X
CARANIA 405202.688 8635139.789 3875 X
HUANGASCAR 409529.310 8573884.037 2533 X
PACARAN 385534.360 8577888.072 721 X X X
TANTA 389341.689 8660300.059 4323 X X
VILCA 410126.563 8660591.417 3864 X
YAURICOCHA 421394.860 8638284.486 4675 X
YAUYOS 401017.784 8618868.347 2327 X X X
HUAROCHIRI 365733.826 8658322.424 3154 X
NICOLAS 381970.587 8573780.842 650 X
El caudal del Rio caete es medido en la actualidad por cinco estaciones hidromtricas:
Socsi, Chavn, Tinco, Aguas Calientes y Tanta, siendo solo la estacin Socsi operada por
SENAMHI y las otra cuatro por ELECTROPERU. Para el presente proyecto se conto con
informacin de la estacin SOCSI y registros de caudales realizadas por las asociaciones de
Regantes del valle del Rio Caete en el mismo lugar.
TABLA 4. ESTACIONES HIDROMETEOROLGICAS UTILIZADA EN EL MODELO
NOMBRE VARIABLE ENTIDAD/FUENTE PERIODO LATITUD LONGITUD ELEVACION
SOCSI Caudales SENAMHI 1994 Y 2008 -76.1948 -13.0286 500
SOCSI Caudales REGANTES 1993-2009 -76.1948 -13.0286 500
Segn el estudio Evaluacin y Ordenamiento de los Recursos Hdricos de la cuenca del
Caete en la actualidad las mediciones en la estacin SOCSI son efectuadas controlando la
altura de la lamina de agua en un limnimetro ubicado al lado de la caseta limnigrafica
(margen derecha aguas arriba del puente) que no cuenta con un limnigrafo, aguas abajo, en
el puente Socsi, es aforado el rio con un correntmetro Universal en forma diaria en la
17
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
poca de estiaje, con lo cual se calibra continuamente la curva altura-caudal. En poca de
lluvias se realiza aforos solamente hasta que el caudal lo permita, puesto que la estacin no
cuenta con los respectivos equipos accesorio. El emplazamiento de la estacin no es el ms
adecuado, se ubica inmediatamente aguas debajo de una curva, adems que en la seccin de
control se produce acumulacin de material solido de arrastre (margen izquierdo). Otro
gran inconveniente es que dicha seccin no controla el total del ancho del cauce.
II. METODOLOGA PARA LA MODELACIN
HIDROLGICA
Toda la informacin espacial suministrada fue necesario realizar un proceso de re-
proyeccin al Sistema de Referencia Proyectado Transverse Mercator - UTM Zona 18S,
Coordenadas Falsas N: 1000.000 m y E: 500.000 m, Datum WGS (Elipsoide WGS 1984). Esto
con el fin de estandarizar toda la informacin utilizada y generada de las modelaciones,
debido a que la informacin inicial provena de diversas fuentes y se encontraban en
diferentes Sistemas de Referencia. El proceso de re proyeccin se hizo por medio del
Software ArcGIS teniendo en cuenta la Regionalizacin de los Parmetros de
Transformacin.
Con el fin de cortar la informacin espacial inicial (Cobertura Vegetal, Suelos y Modelo
Digital de Elevacin) dado que abarcaban una regin ms grande, se realiz una primera
delimitacin de la cuenca en el modelo SWAT, quedando definido la forma y el rea de sta.
Con el lmite definido y la informacin ya recortada se inicia el proceso de elaboracin y
edicin de la Base de Datos del modelo ArcSWAT para realizar la simulacin de la Lnea
Base.
A. ASIGNACIN DE INFORMACIN A LA COBERTURA VEGETAL
A las coberturas del suelo presentes en las cuencas (Figuras 5, 6 y 7) se les realiz la
relacin de los cdigos de las coberturas de la base de datos de SWAT (Tabla 5) y se
modificaron los valores del nmero de curva para cada uno de stas de acuerdo a la
experiencia obtenida en el comportamiento de las coberturas en la Regin de los Andes.
18
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
TABLA 5. ASIGNACIN CDIGO SWAT PARA LAS COBERTURAS DEL SUELO DE LA CUENCA CAETE.
CN: Numero de curva; n: Coeficiente de Manning
B. ESTRUCTURACIN DE LA INFORMACIN DE SUELOS
Esta informacin contaba con el estudio fsico qumico de las unidades de suelos en la zona
de estudio; pero careca de las variables de Conductividad, Densidad Aparente y
Disponibilidad de agua en el suelo. Por tal razn y teniendo en cuenta que estos datos son
indispensables para la modelacin hidrolgica en SWAT, se determinaron sus valores
utilizando el tringulo textural y la herramienta Soil Characteristics Tool (Saxton, 1986).
Se debe de tener presente que la herramienta de Saxton, 1986, funciona mejor para suelos
minerales con bajo contenido de materia orgnica (mximo 8%) y alta densidad aparente.
En la Tabla 8 se encuentran descritos los valores para las primeras 19 unidades de suelo a
emplear dentro de la cuenca y en el Anexo 1 se encuentran la descripcin desde la unidad
20 hasta la 46.
COBERTURA COD. SWAT CN n REA (Ha.)
Matorrales
MESQ 30-46-63-66
0.10 188600
Planicies costeras y estribaciones andinas sin vegetacion RNGB 39-61-74-80
0.15 78730
Tierras altoandinas sin vegetacin
RNGB 39-61-74-80
0.15 46290
Herbazal de tundra BLUG 31-59-72-79
0.10 157100
Lagos y Lagunas
WATR 92-92-92-92
0.01 3225
Pajonal / Cesped de puna FESC 31-59-72-79
0.10 109100
Bofedal
WETL 49-69-79-84
0.05 1271
Queoal OLIV 45-66-77-83
0.14 7095
Pajonal de puna
FESC 31-59-72-79
0.10 174
Cultivos agropecuarios AGRL 67-77-83-87
0.14 6785
Nevados
RNGE 49-69-79-84
0.15 4758
19
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
TABLA 6. CARACTERSTICAS FSICO-QUMICAS DE LOS SUELOS EN LA CUENCA CAETE
No. SNAM (K) USLE
SOL_ Z
profundidad
(mm)
SOL_BD
densidad
aparente
(g/cm3)
SOL_AWC
disponibilidad
de agua en el
suelo (mm/mm)
SOL_K
conductividad
hidraulica
(mm/hr)
SOL_CBN
carbono (%)
%
ARCILLA
(CLAY)
% LIMOS
(SILT)
% ARENA
(SAND)
% ROCK
0.13 225 1.34 0.26 212.59 2.42 13.5 29.0 57.5 20
0.13 500 1.41 0.28 114.94 1.99 21.0 28.5 50.5 20
0.14 1100 1.45 0.29 93.23 1.51 23.5 28.5 48.0 20
0.14 200 1.56 0.17 835.23 0.51 9.0 14.0 77.0 20
0.14 350 1.57 0.17 920.35 0.39 9.0 14.0 77.0 20
0.18 550 1.59 0.26 595.82 0.39 12.0 26.0 62.0 20
0.19 750 1.47 0.42 112.85 1.00 15.0 58.0 27.0 20
0.11 200 1.40 0.17 378.77 1.96 9.0 14.0 77.0 40
0.14 500 1.54 0.17 562.06 0.71 9.0 14.0 77.0 40
0.18 1000 1.59 0.24 402.09 0.42 12.0 26.0 62.0 40
0.20 1350 1.56 0.42 101.15 0.42 15.0 58.0 27.0 40
0.14 170 1.42 0.24 232.33 1.78 12.8 26.4 60.8 40
0.16 442 1.56 0.23 263.94 0.72 15.7 21.3 63.0 40
0.17 950 1.59 0.24 317.46 0.44 14.8 23.3 62.0 40
0.13 157 1.40 0.28 158.46 1.91 15.7 25.5 52.2 40
0.16 464 1.51 0.27 174.22 1.09 16.9 27.9 55.1 40
6 1.2d 0.16 188 1.52 0.23 357.03 0.97 11.0 23.8 65.3 40
7 1.2i 0.16 188 1.52 0.23 357.03 0.97 11.0 23.8 65.3 40
0.17 200 1.49 0.23 377.92 0.83 11.3 23.7 65.0 20
0.17 533 1.60 0.25 580.62 0.32 13.3 24.7 62.0 20
9 1.3b 0.16 167 1.52 0.22 433.63 0.97 12.0 21.0 67.0 20
0.11 200 1.05 0.19 299.41 4.84 8.0 15.0 77.0 20
0.15 350 1.59 0.22 323.98 0.49 20.0 13.0 67.0 20
0.15 900 1.60 0.25 239.76 0.23 23.0 16.0 61.0 20
0.13 150 1.44 0.22 278.22 1.77 14.0 19.0 67.0 20
0.15 400 1.55 0.24 277.19 0.80 18.0 19.0 63.0 20
0.12 800 1.49 0.23 165.49 1.67 21.0 13.0 66.0 20
0.13 175 1.38 0.30 74.23 2.16 25.0 31.0 44.0 20
0.15 400 1.48 0.30 61.46 1.05 27.5 28.0 44.5 20
0.15 700 1.50 0.30 40.47 0.71 31.0 26.5 42.5 20
0.14 150 1.30 0.27 173.33 2.62 13.7 34.3 52.0 40
0.16 350 1.50 0.25 315.42 1.14 10.7 28.3 61.0 40
0.16 717 1.48 0.27 212.95 1.17 14.0 30.3 55.7 40
0.18 1133 1.55 0.32 179.86 0.66 15.3 38.0 46.7 40
0.13 200 1.38 0.25 138.80 2.31 18.0 25.5 56.5 40
0.16 600 1.54 0.25 176.20 0.83 18.5 23.5 58.0 40
0.15 1100 1.54 0.26 101.84 0.75 23.5 22.0 54.5 40
0.15 250 1.49 0.23 255.76 1.27 13.8 22.8 63.4 40
0.15 450 1.55 0.19 528.73 0.66 9.5 17.0 73.5 40
0.16 690 1.56 0.23 284.46 0.73 15.0 21.0 64.0 40
0.12 175 1.37 0.28 57.54 2.39 27.0 26.3 46.8 40
0.16 463 1.50 0.30 52.52 0.82 27.0 30.0 43.0 40
17 2.1i 0.12 200 1.28 0.28 38.00 3.31 32.5 29.0 38.5 40
18 2.1e 0.13 150 1.21 0.24 209.28 3.50 12.0 26.7 61.3 40
0.16 200 1.48 0.34 115.87 1.04 19.5 40.5 40.0 20
0.17 425 1.51 0.34 120.95 0.79 20.0 39.0 41.0 20
0.17 650 1.54 0.31 117.49 0.62 22.5 31.5 46.0 20
15 2.1c
16 2.1d
19 2.4a
13 2.1a
14 2.1b
10 1.4a
11 1.4c
12 1.5d
5 1.2c
8 1.3a
3 1.2a
4 1.2b
1 1.1a
2 1.1b
20
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
C. PREPARACIN Y CORRECCIN DEL MDE
Al igual que la informacin de Cobertura Vegetal y de Suelos; el MDE tambin fue
recortado, adems se le practic una correccin por llenado de sumideros puesto que
es necesario encontrar las zonas donde no es posible que el agua drene hacia alguna
de las direcciones, pues la direccin del flujo determina el sentido y orientacin en que
fluye el agua, de esta forma se puede observar cmo se acumula el flujo. Todas estas
trayectorias de drenaje de agua se generan en funcin de las condiciones topogrficas.
Fuente: Baquero Olga Luca.
FIGURA 6. DIRECCIN DEL FLUJO DE AGUA
Como se puede apreciar en la siguiente figura el pxel 763 est rodeado de alturas mayores,
si existiese un flujo de agua que pasara por esa zona, sta se quedara estancada all y no
permitira la generacin de un drenaje continuo y coherente de acuerdo a la topografa de la
zona.
21
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
FIGURA 7. ESTANCAMIENTO DEL FLUJO DE AGUA
D. ESTRUCTURACIN DE LA INFORMACIN HIDROMETEOROLGICA
En primer medida se realiz la proyeccin al Sistema de Referencia seleccionado para el
proyecto. Posteriormente se realiz el ordenamiento y la estandarizacin de la informacin
de acuerdo al formato en que el modelo SWAT requiere los datos climatolgicos de entrada.
Para realizar la calibracin y validacin del modelo es necesario contar con registros de
caudal y sedimentos en lo posible a nivel diario, el caudal del Rio caete es medido en la
actualidad por cinco estaciones hidromtricas: Socsi, Chavn, Tinco, Aguas Calientes y Tanta,
siendo solo la estacin Socsi operada por SENAMHI y las otra cuatro por ELECTROPERU.
Para el presente proyecto se conto con informacin de la estacin SOCSI para los aos 1994
y 2008 y registros de caudales realizadas por las asociaciones de Regantes del valle del Rio
Caete en el mismo lugar para un periodo de 1993 y 2009. Es necesario aclarar que por el
momento no se cuenta con informacin de Sedimentos registrados que permita realizar la
calibracin y validacin de los resultados obtenidos en la modelacin.
Las Figuras 16 presentan la ubicacin geogrfica de las estaciones de precipitacin,
temperatura y caudales a emplear en la simulacin de las tres cuencas.
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
FIGURA 8. UBICACIN DE LAS ESTACIONES CLIMATOLG
E. GENERACIN DE LAS UN
Definido el lmite de la cuenca con un rea de
subdivididas en subcuencas (Figuras
clculos de caudales y sedimentos presentes
definir las URH (Unidades de Repuesta Hidrolgica).
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
STACIONES CLIMATOLGICAS EN LA CUENCA CAETE
GENERACIN DE LAS UNIDADES DE RESPUESTA HIDROLGICA
Definido el lmite de la cuenca con un rea de 5793.979 km2, es necesario realizar una
vididas en subcuencas (Figuras 9) con el objetivo de mejorar la exactitud en los
clculos de caudales y sedimentos presentes; para posteriormente en el modelo SWAT
definir las URH (Unidades de Repuesta Hidrolgica).
22
36626
.cgiar.org
CAETE
HIDROLGICA
es necesario realizar una
con el objetivo de mejorar la exactitud en los
para posteriormente en el modelo SWAT
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
FIGURA 9. MAPA DE SUBCUENCAS
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
. MAPA DE SUBCUENCAS CUENCA CAETE
23
36626
.cgiar.org
24
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
TABLA 7. PARMETROS GENERALES DE LAS SUBCUENCAS
Una vez establecidos el lmite de la cuenca y la divisin en subcuencas de la misma, se
ingres la informacin de suelos, uso de la tierra, clima, relieve y red de drenaje al modelo
hidrolgico SWAT (Soil and Water Assessment Tool) versin ArcSWAT bajo la plataforma
de ArcGIS (9.3). Este modelo permite analizar variables hidrolgicas de la cuenca como
produccin de agua al caudal, produccin de sedimentos, evapotranspiracin potencial y
real, contenido de agua en el suelo, flujo lateral, percolacin, recarga de acuferos, entre
otros, de manera integral. Es decir, los valores de estas variables son el resultado de la
interaccin agua-suelo-clima-uso-relieve en la cuenca para el periodo de tiempo simulado.
Este tiempo de simulacin esta dado por el periodo de tiempo de informacin climtica
disponible. Es decir, en este caso se realiz una modelacin hidrolgica para un periodo de
19 aos (1991-2009).
Es por esta razn, que el territorio de una cuenca no presenta un comportamiento
hidrolgico homogneo, ya que pueden darse mltiples interacciones de este tipo en la
medida que el suelo, el uso, el clima, y el relieve varan espacialmente. Estas diferentes
combinaciones dan resultado a porciones del territorio con un comportamiento hidrolgico
especfico llamadas Unidades de Respuesta Hidrolgica - URH.
Por medio de la integracin de la informacin mencionada anteriormente, se determinaron
3138 Unidades de Respuesta Hidrolgica (URH) en la cuenca de estudio, las cuales
corresponden a unidades del territorio que presentan condiciones de uso de tierra, suelo y
topografa homogneos y por lo tanto producen un impacto particular sobre la cantidad y
Subcuenca Area (Ha) Lat Long_ Elev ElevMin ElevMax
1 45889.74 -12.12 -75.98 4831 4022 5819
2 15536.61 -12.14 -75.70 4467 3897 4913
3 15628.95 -12.30 -75.66 4747 3773 5036
4 11042.73 -12.25 -75.89 4486 3276 5861
5 37260.81 -12.12 -75.85 4540 3105 5404
6 13535.10 -12.23 -75.72 4209 3122 5010
7 18403.20 -12.38 -75.73 4414 2929 5271
8 42042.24 -12.52 -75.73 4381 2633 5272
9 33211.62 -12.39 -75.88 3733 2310 5623
10 32460.75 -12.54 -76.05 4471 1741 4970
11 23195.16 -12.72 -75.79 4547 1178 5135
12 61682.31 -12.79 -75.71 4613 1095 5115
13 79835.22 -12.64 -75.93 3786 1082 5282
14 51748.47 -12.95 -75.78 3444 986 4860
15 97924.95 -12.87 -76.04 1183 308 4418
Area Cuenca 579397.9
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
calidad del agua de la cuenca (Figura
identificar qu reas en la cuenca son las que contribuyen ms al aporte de sedimentos y
agua al caudal.
FIGURA 10. MAPA DE UNIDADES DE
F. ANLISIS DE SENSIBIL
Para realizar un uso eficiente de los modelos, es necesario realizar un anlisis de
sensibilidad; el anlisis de sensibilidad utilizado en SWAT20
cual combina el mtodo de diseo OAT (One factor At a Time) y el mtodo de muestreo LH
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
(Figuras 10). Con la informacin de las URH, es posible
reas en la cuenca son las que contribuyen ms al aporte de sedimentos y
MAPA DE UNIDADES DE RESPUESTA HIDROLGICA CUENCA CAETE
ANLISIS DE SENSIBILIDAD, CALIBRACIN Y VALIDACIN
Para realizar un uso eficiente de los modelos, es necesario realizar un anlisis de
sensibilidad; el anlisis de sensibilidad utilizado en SWAT2005 es el mtodo LH
cual combina el mtodo de diseo OAT (One factor At a Time) y el mtodo de muestreo LH
25
36626
.cgiar.org
. Con la informacin de las URH, es posible
reas en la cuenca son las que contribuyen ms al aporte de sedimentos y
VALIDACIN
Para realizar un uso eficiente de los modelos, es necesario realizar un anlisis de
05 es el mtodo LH-OAT, el
cual combina el mtodo de diseo OAT (One factor At a Time) y el mtodo de muestreo LH
26
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
(Latin Hypercube), tomando el muestreo de Latin Hypercube como puntos iniciales para el
diseo OAT. El anlisis de sensibilidad identifica los parmetros que tienen una influencia
significativa en los resultados del modelo con respecto a las observaciones reales.
Aunque se han desarrollado mtodos de campo para medir parmetros del suelo (como la
conductividad hidrulica), virtualmente todos los modelos requieren una calibracin para
ajustar al menos algunos de sus parmetros (Sorooshian et., 1998). La calibracin de un
modelo como SWAT consiste en realizar un test de los parmetros de salida (caudal,
sedimentos, escorrenta, etc.) entre los caudales simulados y los observados o tomados en
campo; con la finalidad de poder determinar cules parmetros de entrada deben ser
ajustados en el modelo. En muchos casos el ajuste entre los caudales observados y
simulados obtenidos durante la validacin es menos satisfactorio que el obtenido en la
calibracin, pero es ms representativo de la exactitud de las predicciones que se hagan con
el modelo (Palacios, 1986).
El ajuste de parmetros se realiz de acuerdo con lo propuesto por el manual de SWAT, el
cual indica los parmetros de entrada que ms afectan las salidas del modelo y su orden de
importancia; en el diagrama adjunto se presenta la metodologa empleada para la
calibracin del modelo.
FIGURA 11.
ESQUEMA
METODOLGICO
DE CALIBRACIN
Y VALIDACIN
SI
NO
Entrada de datos periodo de calibracin
Proceso de Simulacin
Resultado Simulacin
Salidas = datos
de entrada
Ajuste de parmetros
Anlisis de la informacin
Informe
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
Los procesos de calibracin y va
simulacin de 19 aos (1991-2009
caudales diarios aforados con algunos huecos para ciertos periodos en los registros de las
asociaciones de Regantes y con solo dos aos (1994 y 2008) para la estacin SOCSI
Segn lo propuesto en el manual de
necesario dividir el rango de datos disponibles en dos grupos
parmetros ms sensitivos en
correspondiente a 1993-2000
se tom el rango 2001-2009. No se
porque dentro de la modelacin
se lleve a cabo.
Fuente: Modificado Manual SWAT.
FIGURA 12. PERIODOS CALIBRACI
Con base en los sedimentos trasportados en el caudal del
aproximacin para ver el aporte de volu
analizada, los valores obtenidos no pudieron ser validados, ya que
cont con informacin de aforos de sedimentos.
agua, se seleccionaron las URH con mayores aportes de
1991 1993 2000 2001
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
Los procesos de calibracin y validacin se realizaron en forma diaria para un periodo de
2009), contando tambin para este periodo con datos de
con algunos huecos para ciertos periodos en los registros de las
y con solo dos aos (1994 y 2008) para la estacin SOCSI
Segn lo propuesto en el manual de Calibracin/Validacin de SWAT (Figura 12
necesario dividir el rango de datos disponibles en dos grupos; la calibracin de los
sensitivos en los procesos modelados se realiz para el periodo
0 y para la validacin o evaluacin de la capacidad del modelo
No se tiene en cuenta los dos primeros ao (1991 y 1992
porque dentro de la modelacin se debe de tener un tiempo para que la recarga de acuferos
. PERIODOS CALIBRACIN/VALIDACIN.
trasportados en el caudal del ro se hizo una primera
el aporte de volumen de sedimentos por hectrea en la cuenca
, los valores obtenidos no pudieron ser validados, ya que por el momento
cont con informacin de aforos de sedimentos. Con el modelo ajustado en produccin de
las URH con mayores aportes de agua, tanto por hectrea como por
Real
2001 2009
27
36626
.cgiar.org
un periodo de
datos de
con algunos huecos para ciertos periodos en los registros de las
y con solo dos aos (1994 y 2008) para la estacin SOCSI.
12) es
; la calibracin de los
ra el periodo
y para la validacin o evaluacin de la capacidad del modelo
los dos primeros ao (1991 y 1992),
se debe de tener un tiempo para que la recarga de acuferos
una primera
cuenca
por el momento no se
produccin de
, tanto por hectrea como por
28
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
el total del rea de la unidad. As mismo se hizo una primera aproximacin a las reas que
deben tener una mayor produccin de sedimentos basndose nicamente en la simulacin
del proceso de prdida de suelo no calibrada, ya que no se cont con informacin de
produccin de sedimentos en la cuenca.
III. RESULTADOS MODELACIN HIDROLGICA
A. ANLISIS DE SENSIBILIDAD DE LOS PARMETROS DEL MODELO
El periodo de anlisis de sensibilidad vara segn la disponibilidad de registros reales, los
cuales en este caso corresponden a un periodo de 19 aos (1991-2009); en la Tabla 8 se
presenta los resultados obtenidos en el anlisis de sensibilidad para los 20 parmetros que
corresponden o tienen influencia en la variable del caudal (Stream Flow) del modelo; para
determinar el grado de sensibilidad de los parmetros evaluados, el modelo estadstico en
SWAT2005 los clasifica en un rango (Rank) de prioridad de 1-20 (rango que vara
dependiendo del nmero de parmetros evaluados); identificando cules parmetros tienen
mayor influencia sobre los resultados de caudal en el modelo. El rango de clasificacin de
parmetros est dado por:
1 muy importante
2 9 importante
10 20 relativamente importante
21 no importante
29
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
TABLA 8. RESULTADO DE ANLISIS DE SENSIBILIDAD PARA LOS PARMETROS DE CAUDAL
En los resultados del anlisis de sensibilidad se identific que de los 20 parmetros
analizados hay 1 parmetro como el ms sensible (sensibilidad global 1) que hace parte del
proceso de enrutamiento en el cauce principal. Adems, hay 8 parmetros importantes
(sensibilidad global >2 y sensibilidad global <= de 9) relacionados a los procesos de Flujo
Base, Escorrenta, Nieve y de movimiento de agua en el suelo. Por ltimo, hay 9 parmetros
"ligeramente importantes" (sensibilidad global >10 y sensibilidad global <=20), y 2
parmetros que no causan cambio a la salida del modelo (sensibilidad global de 21).
Este resultado permite identificar de manera ms clara cuales son los parmetros que
necesariamente deben ser tenidos en cuenta en el momento de la calibracin del modelo. En
el anexo 2, se observan cuales fueron los valores asignados a estas variables despus del
proceso de calibracin.
A manera de resumen se determina que el parmetro ms sensible en el modelo
corresponde al Conductividad Hidrulica Efectiva en el Canal (CH_K2) este parmetro se
calcula para la corriente o canal principal; donde las corrientes en el Modelo SWAT2005 se
pueden clasificar dependiendo de la relacin que existe con el sistema de aguas
subterrneas. Una corriente situada en un rea de descarga que recibe el flujo de agua
subterrnea es una corriente que gana mayor flujo (Figura 13-a). Este tipo de corriente se
caracteriza por un aumento en la descarga aguas abajo. Una corriente situada en un rea de
PARAMETRO RANK COMENTARIO PROCESO
Ch_K2 1 Channel effective hydraulic conductivity [mm/hr] Enrutamiento
Gwqmn 2 Threshold water depth in the shallow aquifer for flow [mm] Flujo Base
Canmx 3 Maximum canopy storage [mm] Escorrentia
Esco 4 Soil evaporation compensation factor Suelos
Timp 5 Snow pack temperature lag factor Nieve
Gw_Revap 6 Groundwater "revap" coefficient Flujo Base
Blai 7 Maximum potential leaf area index Cultivos
Cn2 8 Initial SCS CN II value Escorrentia
Alpha_Bf 9 Baseflow alpha factor [days] Flujo Base
Sol_Z 10 Soil depth [mm] Suelos
Sol_K 11 Saturated hydraulic conductivity [mm/hr] Suelos
Gw_Delay 12 Groundwater delay [days] Flujo Base
Ch_N2 13 Manning's nvalue for main channel Enrutamiento
Slope 14 Average slope steepness [m/m] Geomorfologia
Sol_Awc 15 Available water capacity [mm H20/mm soil] Suelos
Surlag 16 Surface runoff lag time [days] Escorrentia
Biomix 17 Biological mixing efficiency Suelos
Epco 18 Plant uptake compensation factor Evaporacion
Revapmn 21 Threshold water depth in the shallow aquifer for "revap" [mm] Flujo Base
Sol_Alb 21 Moist soil albedo Suelos
Telfono: Colombia: +57 2 4450000;
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali
recarga es una corriente de prdida
disminucin en la descarga aguas abajo. Una corrient
al flujo de aguas subterrneas de dos maneras
que a la vez recibe y pierde el
13-d). Para arroyos perennes con
conductividad efectiva ser igual a cero.
Fuente: Manual Terico SWAT2005.
FIGURA 13. RELACIN ENTRE LA
Dentro de los parmetros que
encontramos:
Umbral de profundidad de agua en el acufero superficial (
de la profundidad del agua requerida para que ocurra el flujo de retorno o aporte de aguas
subterrneas al flujo del cauce principal, cuanto ms alto sea este valor, los aportes de flujo
base hacia la corriente son ms restringidos, al contrario de si este valor es menor los
aportes sern mayores.
Capacidad Mxima de Almacenamiento del Dosel (CANMX)
cantidad de agua que puede ser atrapada en el dosel cuando la cubierta est completamente
desarrollada y que tiene influencia directa en el clculo de la infiltracin, escorrenta
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
recarga es una corriente de prdida de flujo. Este tipo de corriente se caracteriza por una
disminucin en la descarga aguas abajo. Una corriente de prdida de flujo se puede conectar
al flujo de aguas subterrneas de dos maneras (Figura 13-b) o (Figura 13-c). Una corriente
que a la vez recibe y pierde el flujo de agua subterrnea presenta un flujo continuo (Figura
Para arroyos perennes con la contribucin de las aguas subterrneas continua, la
conductividad efectiva ser igual a cero.
. RELACIN ENTRE LA CORRIENTE Y EL AGUA SUBTERRNEA
Dentro de los parmetros que tienen un nivel alto de importancia de sensibilidad,
bral de profundidad de agua en el acufero superficial (GWQMN): este es el umbral
de la profundidad del agua requerida para que ocurra el flujo de retorno o aporte de aguas
l flujo del cauce principal, cuanto ms alto sea este valor, los aportes de flujo
base hacia la corriente son ms restringidos, al contrario de si este valor es menor los
Capacidad Mxima de Almacenamiento del Dosel (CANMX) el cual determina
cantidad de agua que puede ser atrapada en el dosel cuando la cubierta est completamente
desarrollada y que tiene influencia directa en el clculo de la infiltracin, escorrenta
30
36626
.cgiar.org
se caracteriza por una
prdida de flujo se puede conectar
c). Una corriente
flujo continuo (Figura
la contribucin de las aguas subterrneas continua, la
tienen un nivel alto de importancia de sensibilidad,
es el umbral
de la profundidad del agua requerida para que ocurra el flujo de retorno o aporte de aguas
l flujo del cauce principal, cuanto ms alto sea este valor, los aportes de flujo
base hacia la corriente son ms restringidos, al contrario de si este valor es menor los
etermina la mxima
cantidad de agua que puede ser atrapada en el dosel cuando la cubierta est completamente
desarrollada y que tiene influencia directa en el clculo de la infiltracin, escorrenta
31
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
superficial y la evapotranspiracin; ya que la densidad de la cobertura del suelo y la
morfologa de las plantas ejercen una influencia sobre la intensidad de la precipitacin que
cae al suelo y por consiguiente en los cambios de cantidad de agua que alcanza el suelo.
Factor de compensacin para la evaporacin del suelo (ESCO), este parmetro permite al
usuario modificar la profundidad del perfil del suelo a la cual la demanda de evaporacin es
permitida para satisfacer la demanda de evaporacin del suelo en el efecto de la accin
capilar; este parmetro vara entre (0.01 y 1); cuando este valor es reducido, el modelo es
capaz de extraer ms de la demanda de evaporacin de los niveles ms bajos. SWAT2005
supone que el 50% de la demanda de evaporacin se extrae de los primeros 10 mm del
suelo y el 95% de la demanda de evaporacin se extrae de los primeros 100 mm del suelo.
Factor de retraso de la temperatura en el rea de nieve (TIMP), La influencia de la
temperatura del da anterior en el rea de nieve en la temperatura del da actual de esta
rea es controlado por un factor de retraso. El factor de retraso representa la cantidad y
densidad de la nieve, la profundidad de la capa de nieve, la exposicin y otros factores que
afectan la temperatura en el paquete de nieve. TIMP puede variar entre 0,01 y 1,0. Al
acercarse a 1.0, la temperatura media del aire en el da actual ejerce una influencia cada vez
mayor sobre la temperatura de la nieve y la temperatura de la capa de nieve del da anterior
ejerce una influencia cada vez menor. Cuando TIMP tiende a cero, la temperatura de la capa
de nieve ser menos influenciado por la temperatura del aire del da actual.
Factor Alfa para el flujo base (ALPHA_BF), este es el umbral de la profundidad del agua
requerida para que ocurra el flujo de retorno o aporte de aguas subterrneas al flujo del
cauce principal, cuanto ms alto sea este valor, los aportes de flujo base hacia la corriente
son ms restringidos, al contrario de si este valor es menor los aportes sern mayores.
B. CALIBRACIN Y VALIDACIN DEL MODELO
Tomando como base la experiencia obtenida en la modelacin de cuencas en la regin
Andina y de acuerdo a las recomendaciones sugeridas en estos casos por los autores de
SWAT, se modificaron algunos parmetros del modelo segn los resultados obtenidos en el
anlisis de sensibilidad y se realiz una nueva simulacin. A continuacin se presenta de
forma grfica el promedio mensual multianual del caudal simulado vs. el aforado.
Para el proceso de validacin se compararon los caudales diarios observados y simulados,
con los parmetros ajustados en el proceso de calibracin. En las Figuras 14 y 15 se
presentan los datos de produccin de agua diarios tanto para el periodo de calibracin
(1993-2000) como para el periodo de validacin (2001-2009).
32
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
FIGURA 14. CAUDAL SIMULADO VS. AFORADO PERIODO CALIBRACIN
FIGURA 15. CAUDAL SIMULADO VS. AFORADO PERIODO VALIDACIN
TABLA 9. ESTADSTICA DESCRIPTIVA DEL PROMEDIO DEL FLUJO DIARIO
0
100
200
300
400
500
600
1/1/1993 1/1/1994 1/1/1995 1/1/1996 1/1/1997 1/1/1998 1/1/1999 1/1/2000
C
a
u
d
a
l
(
m
3
/
s
)
Caudal Diario Simulado Vs. Registros Socsi (1993-2000)
Simualdo
Observados
0
100
200
300
400
500
1/1/2001 1/1/2002 1/1/2003 1/1/2004 1/1/2005 1/1/2006 1/1/2007 1/1/2008 1/1/2009
C
a
u
d
a
l
(
m
3
/
s
)
Caudal Diario Simulado Vs. Registros Socsi (2001-2009)
Simulado
Observados
1.32 0.04 0.50 0.10 0.09 51.50 0.60 3.72 0.65 0.00 43.27 55.84 promedio
454.67 550.00 max
1.34 2.14 min
0.62 correlacion
0.30 Nash-Sutcliffe
Precipitaci
n (mm)
Escorrentia
(mm)
Flujo
Lateral
(mm)
Aporte
acufero
(mm)
Percolacion
(mm)
Agua en
el suelo
(mm)
Evapotras.
(mm)
Evapotras.
Potencial
(mm)
Lamina de
Agua
(mm)
Sedimentos
(t/ha)
Caudal
Simulado
(m3/seg)
Caudal
Aforado
(m3/seg)
33
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
La Tabla 9 contiene el resumen del balance hidrolgico a nivel diario para todo el periodo
de simulacin, al igual que la correlacin entre los datos observados y los simulados que fue
de 0.62 y el Coeficiente de Nash-Sutcliffe que fue de 0.3. En general los caudales simulados
se ajustan bien a los caudales aforados, obteniendo ajustes razonables donde la mayor
diferencia se presenta en los caudales mximos y mnimos para ciertos eventos puntuales,
considerando que existe un nivel alto de incertidumbre en los datos observados en la
estacin SOCSI y que es necesario ajustar en el modelo en el aporte de desglaciacin con
datos de registros histricos para el periodo de simulacin.
FIGURA 16. CURVA FRECUENCIA DURACIN DE CAUDAL SIMULADO VS OBSERVADO.
A continuacin se presenta de forma grafica los resultados obtenidos a nivel mensual versus
el promedio de los registros de la estacin SOCSI.
0
100
200
300
400
500
600
0.0 0.1 0.3 0.4 0.6 0.7 0.8 1.0
C
a
u
d
a
l
(
m
3
/
s
)
CURVA FRECUENCIA DURACION - DIARIA
Simulado
Aforado
34
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
FIGURA 17. PROMEDIO MENSUAL CAUDAL SIMULADO VS. AFORADO
A partir de estos resultados mensuales se gener una tabla resumen en donde se realiza un
anlisis estadstico del promedio del flujo mensual en cada una de las cuencas. En esta se
evala el criterio de eficiencia de Nash-Sutcliffe, el cual se define como:
( )
( )
2
1
2
1
1
i i
i
n
sim obs
i
n
obs obs
i
Q Q
NS
Q Q
=
=
Dnde:
i
sim
Q : Caudal Simulado en el Periodo i.
i
obs
Q : Caudal Observado en el Periodo i.
obs
Q
: Promedio de Caudales Observados en el Periodo de Anlisis.
Este coeficiente se utiliza para evaluar el poder predictivo de los modelos hidrolgicos y
puede variar desde - hasta 1. Donde una eficiencia de 1 (NS=1) corresponde a un perfecto
ajuste entre los datos del modelo y los observados. Una eficiencia de 0 (NS=0) indica que
las predicciones del modelo son tan precisos como la media de los datos observados,
mientras que una eficiencia menor que cero (NS<0) ocurre cuando la varianza residual
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
Simulado 93.4 105.3 107.9 65.2 36.0 20.6 22.3 22.5 25.0 37.3 36.5 57.8
Observados 91.7 140.6 150.2 87.4 44.0 20.8 13.9 11.9 10.6 11.8 20.1 44.5
0
20
40
60
80
100
120
140
160
C
a
u
d
a
l
(
m
3
/
s
)
Caudal Simulado Vs. Registros Socsi
35
Telfono: Colombia: +57 2 4450000; USA: +1 650 8336625 / Fax: Colombia: +57 2 4450073; USA: +1 650 8336626
Direccin postal Colombia: CIAT, A.A. 6713, Km 17, Recta Cali-Palmira, Cali, Colombia [email protected] www.ciat.cgiar.org
(numerador de la expresin), es mayor que la varianza de los datos (denominador de la
ecuacin). Bsicamente, cuanto ms cercano a 1, ms exactos son los resultados.
TABLA 10. ESTADSTICA DESCRIPTIVA DEL PROMEDIO DEL FLUJO MENSUAL
Estadsticos
Estacin Socsi
Caudal
Simulado
(m
3
/s)
Caudal Aforado
(m
3
/s)
Media 52.5 54.0
Coeficiente de Nash 0.82
Correlacin 0.96
Para la simulacin mensual, en general, el flujo simulado se adapta muy bien con los valores
observados. La eficiencia de Nash y el coeficiente de correlacin son valores que
demuestran un buen ajuste para el perodo de anlisis, lo que demuestra la gran capacidad
que tiene el modelo en predecir los valores promedios mensuales de caudal.
C. LOCALIZACIN ESPACIAL DE URH Y SU PRODUCCIN DE SEDIMENTOS Y
AGUA