Estudio Hidráulico Cuenca El Morro
Estudio Hidráulico Cuenca El Morro
MODELACIÓN HIDROLÓGICA
INFORME FINAL
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Página 1 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Página 2 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
1 INTRODUCCIÓN
La cuenca hidrográfica El Morro ha sufrido un importante proceso de erosión hídrica en las últimas
décadas, manifestando un crecimiento exponencial de su red de drenaje superficial a partir de la
generación abrupta de nuevos cauces sobre el territorio y su posterior desarrollo (Contreras et al.
2012). El origen de los nuevos cursos se encontraría dado a partir de remociones en masa del terreno,
motivadas por procesos de auto-colapso espontáneo de depósitos sedimentarios loéssicos en estado
suelto, desencadenándose en sectores que manifiestan anegamientos locales asociados a niveles
freáticos cercanos a la superficie (Jobbágy et al. 2015).
Este fenómeno es seguido por eventos de erosión retrógrada causada a partir de mecanismos de
sapping, los cuales son catalizados por el aporte de flujos subsuperficiales y las fuertes pendientes del
terreno. A su vez, el proceso continúa avanzando hacia aguas abajo del sector inicial del colapso y
sobrepasando el explayamiento aluvional inicial, seguido por el desarrollo de fenómenos de erosión
fluvial que provocan la rápida profundización y crecimiento de los valles (Becker 2006, Barbeito et al.
2008 y 2009, Degiovanni et al. 2008). Estos episodios de expansión horizontal y vertical suelen ser
abruptos y coinciden con eventos de lluvias intensas, sin que sean necesariamente excepcionales en
el registro histórico (Jobbágy et al. 2015)
De esta manera, la formación de nuevos ríos trae consigo la remoción y transporte de grandes masas
de sedimento que luego se depositan aguas abajo en forma de aluvión, cubriendo grandes extensiones
de terreno y causando la obstrucción de vías de comunicación, fragmentación y pérdida de áreas
productivas y diversos daños sobre infraestructura existente.
Desde el punto de vista hidrológico se observa un ascenso sostenido de los niveles freáticos
generalizado en la región, el cual avanzó progresivamente en las últimas décadas (Barbeito et al. 2008
y 2009, Galván y Collado 2009, Contreras et al. 2012, Sánchez et al. 2013). La causante principal de
este proceso se encontraría dada a partir del cambio en el balance hidrológico que sufre el territorio,
motivado fundamentalmente por el cambio en el uso de suelo y la dinámica subsuperficial resultante
entre los excesos hídricos y el basamento hidrogeológico de la cuenca (Galván 2012, Contreras et al.
2012, Sánchez et al. 2013, Sáenz et al. 2016).
Página 3 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Vale resaltar que la configuración hidrogeológica de las áreas afectadas se caracteriza por la presencia
de un basamento rocoso cercano a la superficie, conformando un acuífero libre de baja potencia y en
general variable tanto en sentido norte-sur como en dirección este-oeste. Dicha variabilidad
condicionaría de manera significativa la respuesta de los niveles freáticos frente a la dinámica hídrica
existente (Barbeito et al. 2008 y 2009, Sánchez et al. 2013).
Se destaca como ejemplo el caso del denominado Río Nuevo en la cuenca El Morro, el cual inició su
formación alrededor de la década de 1980 y en la actualidad presenta valles de hasta 25 metros de
profundidad y 50 metros de ancho. Otro ejemplo reciente corresponde al caso del río La Paraguaya en
la provincia de Córdoba, el cual se desarrolló en forma abrupta a finales del año 2015 generando
obstrucción y daños sobre la Ruta Nacional N°8.
2 OBJETIVOS
Dada la evidencia existente respecto al aumento de los niveles freáticos y su relación con el inicio de
los nuevos cursos, se propone realizar un estudio basado en modelación hidrológica continua, con el
objeto de valorar la incidencia del cambio de uso de suelo y el aumento de las precipitaciones sobre
las recargas del acuífero freático en la cuenca del Morro.
Por su parte, dadas las escalas espaciales y temporales involucradas en ciertos procesos, la
evaluación de los sistemas hídricos únicamente mediante el análisis de periodos extendidos resulta
insuficiente para comprender todos los fenómenos involucrados en la alteración del medio físico. Los
eventos de crecida son un ejemplo claro de esto último, requiriendo el estudio de variables y
condicionantes a corto plazo para evaluar la relación lluvia-escorrentía con la dinámica hidráulica
resultante. De esta manera, el siguiente objetivo de este estudio consiste en evaluar el impacto del
cambio de uso de suelo, el aumento de áreas anegadas y la evolución de la red hídrica sobre la
magnitud de los eventos de crecida que afectan a la cuenca del Morro.
El presente análisis se justifica a los fines de diagnosticar la situación actual del medio físico y servirá
de base para la toma de decisiones respecto a la planificación del territorio y tendencias de crecimiento
en la infraestructura.
Página 4 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
3 METODOLOGÍA
La herramienta implementada corresponde al SWAT (Soil and Water Assessment Tool). La misma, de
dominio público, es un modelo hidrológico a escala de cuenca desarrollado por el Servicio de
Investigación en Agricultura (Agricultural Reaserch Service) del USDA (United States Department of
Agriculture, EE.UU.), con el objeto de predecir el impacto de las prácticas de manejo del suelo,
sedimentos y químicos agrícolas en áreas de aporte de gran extensión y complejidad, considerando
diversos tipos de suelo, usos y condiciones de manejo en periodos extendidos de tiempo.
Figura 4-1. Representación esquemática de los principales procesos modelados a partir de SWAT.
Página 5 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
El modelo parte del balance hidrológico para evaluar todos los fenómenos que ocurren en la cuenca.
La simulación hidrológica se puede separar en dos grandes grupos. Por un lado se caracterizará la
fase terrestre del ciclo hidrológico. Esta fase controla la cantidad de agua, sedimentos, nutrientes y
pesticidas que ingresaran al curso principal de cada subcuenca. La segunda fase la constituye el
traslado de caudales y sedimentos a lo largo de la red de ríos y arroyos de la cuenca hasta su punto
final o de cierre.
Donde:
Con el objeto de mejorar la discretización del modelo a nivel de subcuenca, SWAT permite discretizar
espacialmente el dominio de cálculo en Unidades de Respuesta Hidrológica (HRU´s), las cuales
comparten propiedades derivadas del mismo tipo de suelo, cobertura vegetal y pendiente topográfica.
Los volúmenes de escorrentía son simulados a escala diaria y a nivel de HRU, a partir de una
modificación sobre método del número de curva propuesto por el S.C.S, en el cual los parámetros
asignados a cada tipo de cobertura varían en forma no lineal con el contenido de humedad del suelo.
SWAT utiliza un único modelo de crecimiento basado en Heat Units (HU) para todos los tipos de
cobertura, diferenciando los cultivos anuales y perennes. El crecimiento de las raíces y la evolución
del Índice de Área Foliar son descriptos en forma explícita a partir de funciones adaptadas para cada
tipo de cobertura. Asimismo, el modelo es capaz de representar periodos de stress generados por
temperatura, falta de nutrientes y baja disponibilidad de agua, así como también simular distintas
consignas en el uso de la tierra (siembra, fertilización, cosecha, rotaciones, irrigación).
Las variables climáticas requeridas por SWAT consisten en la precipitación diaria, la temperatura
máxima y mínima daría del aire, la radiación solar, la velocidad del viento y la humedad relativa, y
Página 6 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
pueden provenir de registros de valores medidos o pueden generarse durante la simulación a través
de modelos climáticos. Los modelos climáticos WG (Weather Generator), se utilizan (de ser necesario)
para generar valores diarios de las variables de entrada del modelo a partir de valores medios
mensuales.
La interacción evapotranspirativa de la zona vadosa con la zona de raíces puede aproximarse a partir
del factor EPCO (factor de compensación de captación), el cual varía desde 0 a 1. De esta forma, en
aquellos periodos en los cuales la zona de raíces presenta un déficit hídrico, la demanda de cada
cobertura puede ser satisfecha en forma parcial captando parte del contenido de agua de la zona no
saturada (la compensación se incrementa a medida que el factor se acerca a 1). En forma análoga, el
factor ESCO (factor de compensación de la evaporación del suelo) se utiliza para representar la
evaporación directa del suelo hacia la atmósfera en función de la profundidad. A medida que el
parámetro se acerca a 1, la evaporación se extiende con mayor profundidad sobre el perfil del suelo.
Por su parte, el flujo desde la zona saturada hacia la zona de raíces es simulado a partir del factor
REVAP. De esta manera, el intercambio vertical ascendente se calcula en forma proporcional a la
evapotranspiración potencial de referencia afectada por el factor de compensación que varía de 0 a 1.
Finalmente, el flujo que abandona la zona radicular es computado como recarga. Estrictamente, este
volumen no alcanza la freática en forma directa sino que primero recarga los poros de la zona vadosa
para luego incorporarse a la zona saturada con un consecuente decaimiento y retardo. Si bien la
representación a partir de compartimientos discretos presenta una simplificación respecto a los
procesos mencionados, la misma permite efectuar una primera aproximación sobre el contraste en la
captación del recurso subterráneo para las distintas coberturas dominantes.
El modelo aplicado corresponde al HEC-HMS. El mismo, de dominio público, desarrollado por el Centro
de Ingeniería Hidrológica (Hydrologic Engineering Center) del cuerpo de ingenieros de la armada de
los EE.UU. (US Army Corps of Engineers), surge como evolución del conocido y ampliamente utilizado
HEC-1, con varias mejoras con respecto a éste, entre las que destaca la interfase gráfica de usuario
que facilita las labores de pre-proceso y pos-proceso, y la posibilidad mediante HEC GeoHMS de
conectar con el sistema de información geográfica QGIS.
Página 7 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
El programa es un sistema de modelación generalizada que permite simular varias cuencas y sub
cuencas. El modelo se construye separando el ciclo hidrológico en partes manejables y definiendo los
límites alrededor de la cuenca de interés. El software constituye un sistema integrado de trabajo,
incluyendo:
Se efectuó una recopilación y análisis estadístico detallado de todas las fuentes de datos
meteorológicos y pluviométricos disponibles del Servicio Meteorológico Nacional (SMN), Instituto
Nacional de Tecnología Agropecuaria (INTA) y la Red de Estaciones Meteorológicas de San Luis
(REM). En forma complementaria, se procesaron bases de datos satelitales a partir de los modelos
globales GPM (Global Precipitation Measurement) y TRMM (Tropical Rainfall Measuring Mission).
El Tropical Rainfall Measuring Mission (TRMM) es una misión espacial conjunta entre la Administración
Nacional de la Aeronáutica y del Espacio (NASA) y la Agencia de Exploración Aeroespacial Japonesa
(JAXA), diseñada para monitorear y estudiar las lluvias tropicales y sub-tropicales a escala global. Con
registros desde el año 1998, este recurso brinda información acerca de la evolución espacio-temporal
de las precipitaciones con una resolución temporal de 3 horas y una resolución espacial de 0.25º (≈25
Km). El producto implementado en el presente estudio corresponde al algoritmo 3B42.v7, el cual realiza
un ajuste de los registros utilizando información satelital de otras misiones, tales como GMS, GOES-
E, GOES-W, y ajustes a partir de mediciones realizadas en el terreno, con el objeto de mejorar la
confiabilidad de la información suministrada.
Por otra parte, el Global Precipitation Measurement (GPM) consiste en una red internacional de
satélites diseñada para estudiar las lluvias a escala global, extendiendo el dominio espacial respecto
a su sucesor (TRMM) y mejorando su resolución espacio-temporal. Con registros desde el año 2014,
este recurso brinda información pluviográfica distribuida con una resolución temporal de 30 minutos y
una resolución espacial de 0.1° (≈10 Km). El producto implementado en el presente estudio
corresponde al algoritmo Integrated Multi-satellite Retrievals (IMERG), el cual efectúa una inter-
calibración entre los distintos satélites que componen la constelación del GPM, en conjunto con
validaciones a partir de mediciones en el terreno.
Página 8 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
De esta manera se confeccionó una serie histórica de 50 años (1968-2018) de datos de precipitación
diaria a partir de los registros de la estación INTA Villa Mercedes (1968-2018). Asimismo, se
recopilaron diversos estudios antecedentes que involucran el desarrollo de curvas Intensidad-
Duración-Recurrencia y curvas de Decaimiento Areal, con el objeto de evaluar las características
espacio-temporales de los eventos extremos de precipitación en el área de estudio.
Figura 5-1. Curvas IDF antecedentes para la ciudad de Mendoza (INA CRA, 2008).
Figura 5-2. Evaluación espacio-temporal de la precipitación a través del GPM, cuenca del Morro.
Página 9 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
El modelo digital de elevaciones (MDE) fue construido a partir de topografía radar obtenida del SRTM
(Shuttle Radar Topographic Mission) en su resolución de 30 m. La Misión Topográfica Shuttle Radar
(SRTM) es un proyecto internacional entre la Agencia Nacional de Inteligencia-Geoespacial (NGA) y
la Administración Nacional de la Aeronáutica y del Espacio (NASA). Esta base cartográfica ha sido
ampliamente utilizada en diferentes campos del conocimiento relacionados con la geomática y puede
descargarse gratuitamente a través de Internet. Se adjunta a continuación una tabla de referencia con
los valores de precisión especificados en la documentación técnica oficial del SRTM.
Tabla 5-1. Especificaciones técnicas del Shuttle Radar Topographic Mission. (SRTM Technical Guide)
A su vez, se recopilaron cartas topográficas con curvas de nivel, complementadas con una serie de
puntos trigonométricos (Red de Alta Precisión, Precisión y Topográfica) pertenecientes a la base de
datos del Instituto Geográfico Nacional (I.G.N), con el objeto de validar la información topográfica
implementada. Finalmente se realizó una calibración de ambas fuentes de datos, obteniendo valores
satisfactorios en las correlaciones resultantes y definiendo así el modelo digital de elevaciones (MDE)
a implementar en los cálculos.
Figura 5-3. Calibración del M.D.E. en base a la Red de Alta Precisión del I.G.N.
Página 10 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 5-4. Modelo digital de elevaciones (SRTM), curvas de nivel y red de puntos altimétricos (I.G.N).
Se elaboraron detallados mapas de cobertura vegetal discretizados por periodos (1975, 1985, 1995,
2015) en base al procesamiento de imágenes satelitales, análisis del Índice de Vegetación de
Diferencia Normalizada (NDVI) y verificación de campo. Las clases fueron agrupadas en 3 coberturas
principales: cultivos anuales (maíz y soja), bosque nativo y pasturas/pastizales. La distribución espacial
por sub- cuenca fue procesada con ayuda de herramientas S.I.G, obteniendo así una evolución
espacio-temporal discreta de la cobertura vegetal en la zona de estudio.
Figura 5-5. Evolución temporal de la cobertura vegetal para la cuenca del Morro.
Página 11 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Página 12 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Página 13 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
A partir del análisis de fotografías aéreas en la región central de la cuenca del Morro durante el periodo
1903-1962 se estima que la cobertura vegetal nativa se encontraba conformada por monte (50%) y
pasturas/pastizales (50%) (Santoni, 2004). Se presenta a continuación un cuadro resumen con los
escenarios considerados:
Página 14 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
En el caso particular de la cuenca del Morro se identifican dos unidades principales en términos de los
depósitos cuaternarios superficiales, señaladas como planicies loéssicas arenosas y planicies
medanosas, las cuales cubren aproximadamente el 90% de la región intermontana. La primera de ellas
se caracteriza por la presencia de sedimentos eólicos loessoides y arenosos intercalados, apoyados
sobre basamento cristalino y depósitos piroclásticos terciarios. Los suelos se identifican como
Molisoles y Entisoles con textura franco arenosa, compuestos por arcilla (7%), limo (23%) y arenas
finas a muy finas (70%). El ambiente presenta un neto predominio de procesos de denudación y erosión
hídrica laminar y concentrada en cárcavas sobre un relieve ondulado con alternancia de lomas amplias,
bajos plano-cóncavos y encajonamientos sobre valles fluviales (Barbeito et al. 2008). La vegetación
predominante consiste en remanentes de bosque y pastizales nativos (algarrobo, quebrachos y caldén)
con un fuerte reemplazo por uso agrícola.
Por otra parte, la planicie medanosa presenta un relieve ondulado compuesto por paleomédanos de
baja altura fijados por vegetación. El paisaje incluye cubetas de deflación con afloramientos freáticos
en superficie, encadenamiento de lagunas por colmatación y acumulaciones aluvionales por
explayamientos de los cursos fluviales que derivan de la Planicie Loéssica (Barbeito et al. 2008). Se
observan procesos fluviales de fuerte actividad, motivados por el inicio y posterior crecimiento de la
red hidrográfica, determinando una progresiva acumulación y transporte de sedimentos por
explayamientos y posterior sobrepaso de los mismos. A su vez, se identifican amplias zonas afectadas
por salinización de suelos, asociadas a niveles freáticos aflorantes o cercanos a la superficie.
Página 15 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 5-11. Extracción de muestras alteradas de suelo sobre valle del Río Nuevo.
Con el objeto de profundizar el conocimiento de los suelos dominantes en la región desde un enfoque
geotécnico, se realizó una caracterización en detalle de los depósitos cuaternarios superficiales a partir
de ensayos de laboratorio sobre muestras alteradas de suelo extraídas en uno de los valles que
conforman la cabecera del Río Nuevo. A su vez, se efectuó una recopilación y análisis de estudios
antecedentes disponibles en la región, los cuales incluyen diversos sondeos en las localidades de Villa
Mercedes, Juan Jorba, Fraga y las represas El Chañar y Zelegua (Córdoba).
Los estudios de campo indican la presencia de un manto sedimentario superficial conformado por
depósitos cuaternarios de origen eólico (loessoides), compuestos por arenas muy finas con un alto
contenido de limo no plástico y algunas concreciones calcáreas de fuerte cementación con nódulos
compactos. Ensayos de picnometría indican pesos específicos de partículas de 2.5 gr/cm3 mientras
que su densidad según ensayos Proctor Normal se reduce en el entorno de 1.5-1.6 gr/cm3. La
categorización según el Sistema Unificado de Clasificación de Suelos corresponde a suelos
limoarenosos (ML) no plásticos. El sistema H.R.B los clasifica como suelos A-4 (5).
Por su parte, se observa una matriz de suelo de carácter macroporoso con una graduación muy pobre
de tamaños, con baja capacidad de retención de humedad y gran eficiencia al drenaje. Dichas
propiedades favorecen la infiltración y posterior percolación hacia el manto freático, potenciando la
recarga del acuífero subterráneo.
Figura 5-12. Preparación de muestras alteradas de suelo obtenidas del valle del Río Nuevo.
Página 16 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Vale resaltar que las características físicas del material revelan la presencia de suelos potencialmente
colapsables y dispersivos. El análisis efectuado se realizó en base a ensayos Pin Hole sobre muestras
alteradas de suelo compactadas al 95% del Proctor y estudios antecedentes de doble edometría sobre
muestras inalteradas de suelo. La colapsibilidad de estos mantos se explica a partir de la sensibilidad
de su comportamiento mecánico frente a los cambios de humedad, en los cuales se genera una
variación volumétrica brusca de su matriz a partir de la saturación, sin necesidad de un aumento en
las presiones efectivas sino bien debido a una atenuación de los esfuerzos capilares intergranulares.
En general, este tipo de materiales pertenecen a fracciones loessiales con estructuras macroporosas,
granulometrías predominantemente finas de pobre distribución y con fuertes contenidos de limo. Si
bien la fracción de arcilla es relativamente escasa, vale destacar que la misma presenta una influencia
significativa sobre el comportamiento mecánico y la estructura intergranular de estos suelos (Rocca et
al. 2006).
Página 17 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
materiales sin plasticidad y con bajos contenidos de arcilla. Sin embargo, dada la semejanza que
presentan los patrones morfológicos erosivos observados en campo, los antecedentes observados en
obra, la analogía de los resultados obtenidos a partir de tratamientos correctivos y el análisis de
distintos ensayos de laboratorio, es usual definir este tipo de suelos como dispersivos. Si bien esta
clasificación sigue siendo materia de discusión entre diversos autores, se puede establecer que la
dispersividad de estos suelos no se encuentra fundada en factores fisicoquímicos sino en las
características geomecánicas de la matriz del suelo, las cuales manifiestan de forma concluyente la
gran sensibilidad de estos materiales frente a la erosión.
Estudios geofísicos antecedentes a partir de sondeos eléctricos verticales en la cuenca del Morro
indican la presencia de un basamento rocoso cercano a la superficie, conformado por rocas
sedimentarias, metamórficas e ígneas efusivas, atribuidas al periodo Paleozoico y Terciario. La
estructura de este basamento configura un manto sedimentario de espesor variable tanto en sentido
norte-sur como en dirección este-oeste. Esta variabilidad condicionaría de manera significativa la
respuesta de los niveles freáticos frente a la dinámica hídrica existente en la región (Sánchez et al.
2013).
Desde un punto de vista estructural, se identifica un acuífero libre conformado por sedimentos sueltos
(no consolidados) atribuidos al Cuaternario con un marcado predominio de materiales finos (arenas,
limos y arcillas), los cuales conservan su porosidad primaria. El acuífero suprayace a rocas
sedimentarias terciarias, las cuales han perdido parte de su porosidad primaria por compactación y/o
cementación y cuyo límite conforma el basamento hidrogeológico de la cuenca.
Página 18 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
- Afloramiento de agua en superficie asociado a la dinámica existente entre las recargas del
acuífero y su interacción con el basamento hidrogeológico.
- Expansión, crecimiento y elongación de los vasos de agua.
- Remociones en masa abruptas asociadas a eventos extremos de escorrentía y sus
consecuentes explayamientos aluvionales.
- Interconexión de los cuerpos de agua y explayamientos a partir de la consolidación de las
incisiones del terreno por profundización, encajonamiento y expansión lateral de los cauces
formados.
Figura 6-1. Esquema conceptual del desarrollo evolutivo de los nuevos cauces.
Diversos autores (Becker 2006, Doffo 2007, Barbeito et al. 2008, Degiovanni et al. 2009) señalan este
mismo proceso en otras cuencas semiáridas de la región Pampeana Occidental. En este sentido, se
identifica la evolución de un régimen semipermanente-endorreico con predominio de inundación y
sedimentación, a permanente-exorreico donde prevalece la erosión y el transporte. Dichas
investigaciones vinculan la integración de las microcuencas al retroceso y coalescencia de cabeceras
de canales efímeros y discontinuos, asociado a ciclos húmedos (1953-1968, 1972-actualidad) y
acciones antrópicas.
Resulta enfático destacar que el proceso descripto es de carácter reciente en la región y no muestra
señales de estabilización, observándose que la densidad de drenaje ha mantenido un crecimiento
exponencial en las últimas tres décadas, y que hasta el presente no ha cambiado su tendencia
(Santoni, 2004, Contreras et al., 2012, Jobbágy et al., 2015).
Página 19 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 6-3. Explayamiento aluvional sobre el bajo La Salada, posterior a evento extremo de precipitación en el
mes de Febrero del año 2015.
Página 20 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Página 21 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Teniendo en consideración los severos daños registrados en la región a lo largo de las últimas dos
décadas producto de la apertura de cauces y explayamientos aluvionales, se efectuó un análisis
histórico de eventos con el objeto de caracterizar las condiciones climáticas asociadas al desarrollo de
estos fenómenos.
Vale señalar que los impactos documentados sobre la región son variados y de diversa magnitud. Los
mismos incluyen daños puntuales tales como la obstrucción de rutas Nacionales y Provinciales,
fragmentación de áreas productivas, destrucción de caminos vecinales y daños sobre diversa
infraestructura pública y privada. Por otro lado, se han logrado identificar impactos difusos a mayor
escala como el anegamiento de superficies y la salinización de suelos (Jobbágy et al., 2015).
A partir del procedimiento descripto, se procede a identificar los eventos más significativos en las
últimas dos décadas y caracterizar su gravedad. Asimismo, se analizaron las condiciones climáticas
teniendo en consideración la magnitud de la lluvia, su duración y la condición de humedad antecedente
a corto plazo.
En función del presente análisis se logra identificar una concordancia entre la apertura de
cauces/explayamientos aluvionales y la manifestación de eventos extremos de escorrentía. Vale
señalar que dichos procesos no solo corresponden a lluvias máximas históricas, sino que cubren un
espectro amplio definido por la condición de humedad antecedente, el estado de
conectividad/consolidación de la red hidrográfica y las características del evento de precipitación.
Teledetección y geoprocesamiento:
Se estima una superficie total de explayamiento de 200 Ha, concentrado principalmente en la cuenca
media del Río Nuevo. Se reconoce a este evento como uno de los desencadenantes en el nacimiento
del cauce, teniendo en consideración las últimas dos décadas. Posteriores eventos en los meses de
Octubre del año 2001 y Enero-Febrero del año 2002, propagan la extensión del explayamiento inicial.
Página 22 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
- La estación INTA Villa Mercedes (EMC) registra un valor de 93.60 mm para el día 30/09/2001.
- El resto de las estaciones pluviométricas registran fallas durante este evento.
- La condición de humedad antecedente resulta grave (alto potencial de escurrimiento, condición
III según S.C.S.).
Figura 6-7. Evolución de la precipitación diaria para el periodo Septiembre-Octubre del año 2001.
Página 23 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Observaciones:
Teledetección y geoprocesamiento:
Se estima una superficie total de explayamiento de 360 Ha, concentrado principalmente en la cuenca
baja del Río Nuevo. Se reconoce a este evento como un estadio previo al encuentro del cauce con la
Ruta Nacional N°8.
Figura 6-9. Explayamiento aluvional asociado a una sucesión de eventos en los meses de Enero-Febrero del
año 2008.
Estudio Hidráulico Cuenca El Morro - Informe Final.docx
Página 24 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
- La estación INTA Villa Mercedes (EMC) registra una serie de eventos entre el 15/01 y el 15/02
(4 eventos de 45 mm, 42 mm, 40 mm y 42 mm, respectivamente).
- Las estaciones La Esquina, La Punilla y Villa Reynolds presentaron fallas durante los eventos
de fines de enero.
- La condición de humedad antecedente es grave para el evento del 29-01-2008 (alto potencial
de escurrimiento, condición III según S.C.S.).
Figura 6-10. Evolución de la precipitación diaria para el periodo Enero-Febrero del año 2008.
Figura 6-11. Evolución de la precipitación cada 3 horas para los días 25 y 26-01-2008.
Página 25 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Observaciones:
- Los valores de precipitación mensual acumulada para Enero de 2008 son elevados. La estación
INTA Villa Mercedes (EMC) registró un total de 193 mm mensuales, lo cual responde a la serie
de máximos históricos (+70% de la media en Diciembre)
- En el trienio 2005-2006-2007 la precipitación anual acumulada fue de 595 mm, 765 mm y 920
mm, superior (+15%) a la media de las últimas 5 décadas, lo cual favoreció las recargas del
acuífero y contribuyó a aumentar el contenido de humedad de toda la zona vadosa. Esta
componente se ve agravada por la precipitación antecedente del mes.
Teledetección y geoprocesamiento:
Se estima una superficie total de explayamiento de 690 Ha, distribuido en las cuencas bajas del
Quebrachal, Río Nuevo y Zanjón del Cerro Negro. Se reconoce a este evento como el encuentro del
Río Nuevo con la Ruta Nacional N°8.
Página 26 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 6-13. Evolución de la precipitación diaria para el mes de Diciembre del año 2009.
Página 27 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Observaciones:
- Los valores de precipitación mensual acumulada para Diciembre de 2009 son elevados. La
estación INTA Villa Mercedes (EMC) registró un total de 206 mm mensuales, lo cual responde
a la serie de máximos históricos (+100% de la media en Diciembre).
- En el trienio 2006-2007-2008 la precipitación anual acumulada fue de 765 mm, 920 mm y
810 mm, superior (+25%) a la media de las últimas 5 décadas, lo cual favoreció
significativamente las recargas del acuífero y contribuyó a aumentar el contenido de humedad
de toda la zona vadosa.
- Vale señalar el contraste en la intensidad de precipitación registrada para la cuenca de Zanjón
del Cerro Negro, lo cual estaría relacionado con la distribución de los explayamientos
observados.
Página 28 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Teledetección y geoprocesamiento:
Se estima una superficie total de explayamiento de 870 Ha, distribuido principalmente en la cuenca
baja del Río Nuevo. Se reconoce a este evento como el de mayor magnitud registrada en las últimas
dos décadas.
Todas las estaciones cercanas presentan fallas durante la noche del 24 de Febrero.
Figura 6-17. Evolución de la precipitación diaria (GPM) para el mes de Febrero del año 2015.
Página 29 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 6-18. Evolución de la precipitación cada media hora para el día 24-02-2015.
Observaciones:
- Los valores de precipitación mensual acumulada para Febrero de 2015 son máximos históricos.
El satélite registra valores cercanos a los 400 mm en la cuenca del Morro (la Estación Villa
Mercedes EEA y La Esquina registraron 380 mm y 305 mm, respectivamente). La precipitación
mostró valores mensuales máximos históricos (+340% de la media en Febrero).
- En el trienio 2012-2013-2014 la precipitación anual acumulada fue de 987 mm, 945 mm y
862 mm, muy superior (+40%) a la media de las últimas 5 décadas, lo cual favoreció
significativamente las recargas del acuífero y contribuyó a aumentar el contenido de humedad
de toda la zona vadosa.
Estudio Hidráulico Cuenca El Morro - Informe Final.docx
Página 30 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
El análisis de datos históricos sobre pozos existentes muestra un ascenso sostenido de los niveles
freáticos generalizado en la región, el cual avanzó progresivamente en las últimas décadas. A su vez,
vale resaltar que la configuración hidrogeológica de las áreas afectadas se caracteriza por la presencia
de un basamento rocoso cercano a la superficie, conformando un acuífero libre de baja potencia y en
general variable tanto en sentido norte-sur como en dirección este-oeste. Dicha variabilidad
condicionaría de manera significativa la respuesta de los niveles freáticos frente a la dinámica hídrica
existente. Debido a la ausencia de series continuas de caudales de superficie, la validación del modelo
se realizó a partir de una valoración cuali-cuantitativa de las distintas salidas, comparando las mismas
con valores de referencia de estudios antecedentes y con el conocimiento local de la zona de estudio.
La definición de las subcuencas fue efectuada a partir del análisis conjunto del MDE y la definición de
una red digital de cauces (Digital Stream Network), utilizando rutinas de geoprocesamiento con el
software QGIS y ArcSWAT. Los límites resultantes fueron comparados con las subcuencas definidas
por Sánchez et.al 2013, obteniendo resultados aceptables. Vale destacar que debido a las bajas
pendientes, la configuración de la cuenca media/baja se encuentra altamente influenciada por la traza
de los nuevos cursos de agua, las canalizaciones y su respectiva vinculación con los sectores
topográficamente deprimidos, dificultando la definición de los límites propuestos.
Para caracterizar la cobertura vegetal actual, se elaboraron detallados mapas de uso del suelo para la
campaña 2015-2016 en base al procesamiento de imágenes satelitales, análisis del Índice de
Vegetación de Diferencia Normalizada (NDVI) y verificación de campo. Las clases fueron agrupadas
en 3 coberturas principales: cultivos anuales (maíz y soja), bosque nativo y pasturas/pastizales. La
distribución final obtenida corresponde a un 59% de cultivos anuales (60% soja y 40% maíz), 38% de
área cubierta por pastura/pastizales y un remanente del 3% de bosques nativos para la cuenca
hidrográfica el Morro. La profundidad máxima de raíces fue estimada en 6 metros para el bosque, 3
metros para pasturas/pastizales y de 2 metros para cultivos de soja/maíz. Con el objeto de mejorar la
discretización del modelo a nivel de subcuenca, se definieron 120 HRU´s, dando como resultado un
área de análisis de 2.200 Km2 para la cuenca del Morro.
Respecto a las forzantes climáticas, vale resaltar que la precipitación media anual en las inmediaciones
de la localidad de Villa Mercedes, correspondiente a mediciones efectuadas en un período de 109 años
es de 603 mm. Dicho promedio se reduce a 552,5 mm evaluando el periodo 1903-1957 y se incrementa
a 654,8 mm entre los años 1957-2011, dando como resultado un aumento promedio anual del orden
de 1,00 mm/año para el ciclo analizado.
Página 31 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 6-20. Evolución de la precipitación anual acumulada para la serie Villa Mercedes (1903-2011).
Figura 6-21. Evolución de la precipitación media móvil para la serie Villa Mercedes (1903-2011).
Figura 6-22. Temperatura máxima y mínima promedio mensual para la serie INTA Villa Mercedes (1968-2013).
Página 32 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 6-23. Humedad relativa y radiación solar promedio mensual para la serie INTA V.M. (1968-2013).
Figura 6-24. Evolución promedio mensual de las principales componentes del ciclo hidrológico.
Se observa que la evapotranspiración potencial media anual alcanza valores cercanos a los 1280 mm
para el periodo analizado, los cuales se acentúan en la época estival en concordancia con las mayores
temperaturas. En escala mensual, esta componente supera los valores medios de precipitación
registrados a lo largo de todo el año, lo cual señala la amplia capacidad del medio atmosférico para
absorber los excesos de agua que pudiesen existir en el sistema. Los resultados medios mensuales
simulados a partir del método de Penman-Monteith y su evolución a lo largo del periodo modelado
fueron comparados con los registros históricos del INTA y un modelo meteorológico global de NASA,
observando valores satisfactorios para toda la serie.
Por su parte, la evapotranspiración media anual entendida como un promedio espacial ponderado para
las distintas coberturas alcanza valores cercanos a los 600 mm para la cuenca del Morro. La evolución
media mensual de esta componente para cada cobertura fue comparada con el Índice de Vegetación
de Diferencia Normalizada (NDVI) de estudios antecedentes (Santoni, 2004), obteniendo una
valoración satisfactoria.
Página 33 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Vale señalar que en dos trienios de las últimas cinco décadas (1983/85 y 2006/08), las precipitaciones
anuales observadas en las inmediaciones de la cuenca del Morro superaron substancialmente los
valores de los módulos pluviométricos característicos, alcanzando valores máximos históricos de 1060
mm (1979), 1105 mm (1984) y 1047 mm (2007), coincidentes con posteriores eventos de remoción en
masa y aluviones de gran magnitud registrados en la región (declaración de Desastre Agropecuario
Provincial en el año 1985 y sobrepaso de corredores viales nacionales en el año 2008). Se observa
que dichos periodos también coinciden con los valores pico de recarga y escorrentía estimados a partir
del modelo.
Figura 6-25. Evolución de la escorrentía y las recargas anuales para la cuenca del Morro.
Página 34 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 6-27. Evolución de la escorrentía y las recargas anuales bajo cultivos de Soja.
Figura 6-29. Evolución de la escorrentía y las recargas anuales bajo cultivos de Maíz.
Página 35 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Página 36 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Se observa que las coberturas de soja y maíz presentan valores promedios anuales de recarga de 45
mm y 57 mm respectivamente, registrando picos que superan los 200 mm (años 2007 y 2014). Los
volúmenes de escorrentía anuales promedio corresponden a 67 mm y 87 mm, con valores máximos
del orden de los 150-200 mm. La evapotranspiración anual promedio para esta clase de coberturas se
ubica entre los 550-580 mm. Por otro lado, las pasturas/pastizales presentaron valores promedios
anuales de recarga de 36 mm, con picos que superan los 150 mm. El volumen de escorrentía anual
promedio corresponde a 22 mm, con valores máximos del orden de los 60-90 mm. La
evapotranspiración anual promedio corresponde a 635 mm. Finalmente, vale resaltar el caso de los
bosques nativos, los cuales mostraron valores de recarga nula para todo el periodo simulado. La
evapotranspiración anual promedio para este tipo de cobertura corresponde a 685 mm. Se observan
volúmenes de escorrentía anuales despreciables, con valores pico que llegan a alcanzar los 40 mm.
Con el objeto de valorar la incidencia del cambio de uso de suelo y el aumento de la precipitación sobre
el balance hidrológico de la cuenca, se compararon los resultados obtenidos con 2 escenarios de
simulación complementarios. El primero de ellos se define a partir de una cobertura nativa sobre la
cuenca El Morro, compuesta por 50% de bosque y 50% de pasturas/pastizales. Por otra parte,
manteniendo el escenario de cobertura original (59% de cultivos anuales, 38% de área cubierta por
pastura/pastizales y un remanente del 3% de bosque) se redujo la forzante pluviométrica en un 15%,
compatible con el escenario del periodo 1903-1957.
Estudio Hidráulico Cuenca El Morro - Informe Final.docx
Página 37 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 6-35. Evolución comparativa de las recargas para los escenarios simulados.
En relación al balance hídrico, se observa que el reemplazo de las coberturas nativas (bosque y
pasturas/pastizales) por cultivos anuales de secano (soja y maíz) ha provocado un incremento de los
volúmenes anuales promedio de recarga y escorrentía para la cuenca del Morro. Asimismo, se han
acentuado los valores picos históricos en aquellos años donde la precipitación superó
substancialmente los módulos pluviométricos característicos. Esto se explica dada la mayor capacidad
de las coberturas nativas para acceder y evapotranspirar el agua de lluvia gracias a una mayor
profundidad de raíces. Por su parte, el aumento del módulo pluviométrico en el último siglo también ha
Página 38 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
demostrado una incidencia sobre las recargas del mismo orden de magnitud que el cambio del uso del
suelo, aunque no así sobre los volúmenes de escorrentía.
Si bien el presente análisis se basa en la comparación volumétrica de las principales componentes del
ciclo hidrológico en una escala macro espacial y temporal extendida (a nivel de cuenca y a partir de
patrones inter-anuales/mensuales), los resultados obtenidos permiten apoyar las observaciones
efectuadas en campo y brinda los lineamientos base para el desarrollo de modelos matemáticos
distribuidos de mayor complejidad.
Vale destacar que los resultados obtenidos con SWAT respaldan las investigaciones realizadas por la
Dra. Santoni (2004) en el marco de los estudios experimentales realizados en la cuenca del Morro, en
los cuales se efectuaron intensas tareas de muestreo y análisis de laboratorio con el objeto de
caracterizar la dinámica vertical del agua y su relación con el conjunto suelo-vegetación. Asimismo,
dichos resultados fueron validados a partir de la aplicación de modelos matemáticos unidimensionales
con el software HYDRUS 1D, obteniendo valores congruentes con la evidencia de campo.
El rendimiento hídrico de las cuencas montañosas, definido como la fracción de la precipitación que
abandona las mismas en forma líquida, no sólo depende del clima y de la configuración topográfica y
geológica de las cuencas, sino también de la vegetación que las ocupa (Zhang et al. 2001; Brown et
al. 2005; Huxman et al. 2005; Jackson et al. 2009; Nosetto et al. 2012). Desde hace casi un siglo, se
han documentado respuestas del caudal de cuencas montañosas tras cambios en la vegetación,
incluyendo alteraciones en el caudal total anual, en su estacionalidad, y en su partición entre "picos de
caudal" o crecientes inmediatamente posteriores a las lluvias y "caudal base" asociado al aporte más
lento del sistema hidrológico subterráneo (Farley et al. 2005). Las transformaciones hidrológicas de
una cuenca no solo afectan la provisión de agua para consumo humano, industrial, y agrícola, sino
también la regulación de crecientes e inundaciones (Jackson et al. 2001; Bradshaw et al. 2007). Desde
el punto de vista de la provisión de agua, los cambios en el caudal base pueden ser más importantes
que los del caudal total, siendo el primero el que define la oferta segura o estable de agua para sus
distintos usos, mientras que desde el punto de vista de la regulación de crecientes, los picos de caudal
son el componente más relevante.
En regiones con balances hídricos negativos (evapotranspiración potencial > precipitación), las áreas
montañosas juegan un papel clave como proveedoras de agua líquida (Viviroli & Weingartner 2004).
Bajo este contexto climático, los territorios llanos y sedimentarios devuelven casi la totalidad de la
precipitación que reciben a la atmósfera por transpiración o evaporación directa (Scanlon et al. 2006;
Jobbágy et al. 2008), mientras que en los sectores montañosos o serranos las altas pendientes y la
rocosidad suelen impedir este consumo evaporativo exhaustivo, favoreciendo la alimentación de
cursos de agua tanto por escurrimiento (transporte superficial) como por drenaje profundo (transporte
subterráneo), aún bajo climas áridos (Jobbágy et al. 2011). En este contexto, la caracterización del
rendimiento hídrico base resulta de vital importancia dado que el agua subterránea es uno de los
recursos más utilizados para suministro de agua potable, ganadería y riego.
El objetivo principal de este análisis consiste en caracterizar la magnitud y variabilidad del rendimiento
hídrico base en regiones cercanas a la cuenca del Morro, a partir de mediciones estacionales
realizadas durante un periodo de 15 años. La zona de estudio abarca seis cuencas serranas ubicadas
sobre la vertiente oriental de la Sierra de Comechingones, con una precipitación media anual de 800
mm y una evapotranspiración potencial que alcanza los 1200 mm/año.
Los registros instantáneos de caudal base se efectuaron en forma estacional sobre las inmediaciones
de la Ruta Nacional Nº 8 y abarcan los arroyos Chaján, La Paraguaya, Las Cortaderas, Zelegua, Suco
y Sampacho.
Página 39 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Desde un punto de vista geomorfológico, las Sierras Pampeanas constituyen una macro unidad morfo
estructural compuesto por una sucesión de mega bloques integrados fundamentalmente por un
basamento plutónico-metamórfico (metamorfitas de edad precámbrica, intruidas en el Paleozoico
inferior por batolitos graníticos). Estas características implican una configuración que las posiciona
como una entidad hidrológica receptora y a la vez emisora de agua, con incidencia fundamental en el
sistema hidrológico subterráneo. Las rocas ígneas y metamórficas de las sierras están drenadas por
importantes ríos y numerosos arroyos menores que, en general, desarrollan una red de drenaje de alta
densidad, condicionada fuertemente por la presencia de fracturas y diaclasas. En conjunto, las sierras
constituyen una zona receptora de precipitaciones, de vital importancia para el funcionamiento
hidrogeológico ya que la mayoría de los ríos y arroyos que la drenan se infiltran al salir de la sierra y
alimentan los acuíferos de pedemonte y llanuras onduladas cercanas. Los arroyos Chaján, Las
Cortaderas, Zelegua y Sampacho son ejemplos de este comportamiento (Blarasin et al. 2014)
Figura 6-38. Ubicación del área de estudio y disposición de la red hidrográfica monitoreada.
Página 40 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Tabla 6-1. Caracterización morfométrica y valores de descarga base para las cuencas de estudio.
Figura 6-39. Evolución interanual del caudal base y precipitación obtenida por TRMM.
Figura 6-40. Evolución interanual del rendimiento hídrico base y precipitación obtenida por TRMM.
Página 41 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 6-41. Evolución temporal extendida del caudal base y precipitación obtenida por TRMM.
Las cuencas estudiadas mantienen una descarga hídrica base de ~100 mm/año, equivalente a una
octava parte de los aportes de la precipitación, mientras que el volumen restante se perdería por
evapotranspiración y escorrentía superficial. La distribución estacional de la descarga base es de 52%-
48% para el periodo Octubre-Marzo y Abril-Septiembre, mientras que la distribución de precipitaciones
corresponde al 80%-20% para los mismos intervalos.
Cabe aclarar que bajo condiciones semiáridas como las de la zona de estudio, se requiere solamente
un pequeño aumento de la evapotranspiración para producir un efecto significativo sobre la descarga
dado que naturalmente la misma representa una fracción minoritaria de los aportes de precipitación
(Jobbágy et al. 2008). De esta manera, resulta importante contemplar estos efectos al proyectar las
políticas de uso del suelo a nivel regional, en especial en áreas potenciales de aporte a los principales
focos de consumo de agua.
Página 42 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Cabe mencionar que este sistema involucra áreas muy activas, que han modificado significativamente
su red de drenaje, con muy pocos relevamientos y antecedentes topográficos e hidrográficos. Por lo
tanto, la evaluación comparativa incluye un análisis de sensibilidad de las diferentes variables
hidráulicas, así como la evaluación cuali-cuantitativa de las diferentes respuestas del sistema ante
escenarios históricos.
Tiene nacientes a partir del aporte de la vertiente occidental de la Sierra de Comechingones y la oriental
de la sierra del Portezuelo, en el primer caso mediante el Aº La Punilla y en el segundo, mediante el
Aº El Portezuelo, Ambos de caudal permanente se unen a la altura de La Cañada, para constituir un
único canal de carácter permanente y bien definido, en algunos tramos encajados entre 8 y 10 m de
profundidad, que luego de un recorrido de 38 km., recibe el Aº Cañada Del Zanjón Negro con nacientes
en la vertiente occidental de las Sierras Comechingones. A partir de este punto de encuentro en la
actualidad tiene lugar un único sistema con desagüe en el río Quinto luego de un recorrido de 73 km.
Tiene nacientes en la parte norte de la depresión central del Cerro El Morro a más de 1700 m.s.n.m.
mediante dos arroyos de caudal permanente que luego de su confluencia, abandonan la depresión
circular por la quebrada De la Cal. A lo largo del ambiente de pedemonte se define la traza de un cauce
fluvial encajonado con un recorrido de 34 km., para luego perder tal condición e infiltrarse a la altura
de la Ea la Guardia. El cauce vuelve a adoptar un carácter permanente por afloramiento de la napa
freática al norte del establecimiento Las Nenas, a partir de donde con sentido norte–sur, evoluciona
hacia un complejo sistema de lagunas y bañados ubicado al norte de la localidad de Juan Jorba. En
este sistema de bañados se integra otra cuenca con nacientes en el faldeo oriental del Cerro El Morro
y en el occidental de la Sierra del Portezuelo. Todo el conjunto originalmente con desagüe N-S en
dirección al río Quinto, en la actualidad se encuentra desviado mediante una canalización paralela a
la Ruta Nacional Nº 8 trasvasando sus excedentes al bajo La Salada, en el que desagua el sistema
del Río Nuevo.
Página 43 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 7-1. Vista aérea del valle del Río Nuevo, cuenca del Morro
Este sistema se organiza a partir de nacientes con desarrollo en la vertiente oriental de la Sierra del
Yulto, que al ingresar al ambiente de la Planicie Medanosa adopta una marcada anarquía,
concentrándose en áreas deprimidas (cubetas de deflación) con y sin afloramiento del nivel freático y
formación de lagunas con encadenamiento durante precipitaciones importantes. En su progresión
hacia el sur el sistema continúa con carácter netamente temporal regido por la interconexión de bajos
interdunas, hasta la depresión de La Salada.
El estudio de la red de drenaje y la delimitación de las subcuencas integrantes del sistema del Morro
fue realizado a partir de información de base del Instituto Nacional del Agua, la Subsecretaría de
Recursos Hídricos y el Instituto Geográfico Nacional, junto con apoyo de herramientas SIG para su
definición detallada.
Se utilizó la extensión ArcSWAT, desarrollada por Stone Environmental Inc., en colaboración con el
Laboratorio de Ciencias Espaciales de Texas para definir los cursos de agua. La extensión ArcSWAT
utiliza un Modelo Digital de Elevaciones (MDE) para calcular el comportamiento de los cursos de agua
de una zona, de acuerdo a la topografía y la pendiente. La red de drenaje obtenida a partir de la
información de base fue corregida empleando imágenes satelitales Landsat actuales, eliminando
algunos cauces menores y redefiniendo otros completamente según resultara necesario.
Una vez ajustada la información de base, se ejecutó la aplicación ArcSWAT. Allí se ingresó el MDE y
luego se incorporó el shape de drenaje corregido, de forma tal que la aplicación interprete los canales
por los que se evacúa el agua de la cuenca y a los que debe llegar cualquier gota de agua que cae
Página 44 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
sobre ella. Con toda esta información la aplicación calcula automáticamente la dirección de los flujos y
las zonas en las cuáles se acumulan. Una vez completado este paso, se configuraron los cierres de
cada subcuenca de forma tal de interceptar los flujos sobre la Ruta Nacional N° 8.
Finalmente se realizó una revisión manual de las subcuencas, pues es posible encontrar diferencias
tanto en las curvas de nivel como en la construcción del MDE. Para ello se compararon los contornos
obtenidos con las curvas de nivel y con imágenes satelitales, verificando que los límites fueran
respetando siempre la topografía divisoria de agua.
Figura 7-2. Configuración de la red hidrográfica y subcuencas sobre la cuenca del Morro.
Página 45 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
En base a las curvas de nivel, el modelo digital de elevaciones, la red de drenaje y el uso de suelo, se
determinaron por medio de herramientas S.I.G. las características hidrológicas de las subcuencas que
conforman el sistema del Morro.
La expresión propuesta por Kirpich (1940) fue desarrollada a partir de estudios del Soil Conservation
Service en siete cuencas rurales de Tennessee con canales bien definidos y pendientes empinadas (3
a 10%).
0.77
L
Tc 0.066
S
Donde:
Página 46 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
La ecuación propuesta por Dooge (1973) fue calibrada para cuencas rurales con áreas entre 140 y 930
km2.
A0.41
Tc 21.88
S 0.17
Donde:
La expresión de Johnstone y Cross (1949) fue derivada a partir del método de Clark aplicado a 19
cuencas con áreas entre 65 y 4200 km2.
0.5
L
Tc 5
S
Donde:
La ecuación propuesta por el Department of Public Works (1995) fue desarrollada a partir del estudio
de pequeñas cuencas en áreas de montaña en California.
0.385
11.9.L3
Tc 60
H
Donde:
La expresión propuesta por Giandotti (1934) fue desarrollada a partir de estudios sobre 12 cuencas
con áreas entre 170 y 70.000 km2.
4 A 1,5.L
Tc
0.8.H
Donde:
Página 47 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
La expresión propuesta por el S.C.S (1975) surge como una suma individual de tiempos de viaje
calculados a partir de velocidades medias de escurrimiento.
L
Tc
60.V
Donde:
La expresión propuesta por Témez (1991) fue desarrollada para cuencas con un tiempo de
concentración entre 0.25-24 hs y una superficie de hasta 3000 km2.
0.76
L
Tc 0,3. 0.25
S
Donde:
La expresión propuesta por Passini es usualmente implementada para cuencas de gran superficie.
L ( A.L)1 / 3
Tc .
A S 0.5
Donde:
La expresión propuesta por Bransby y Williams también es de frecuente uso en cuencas de gran
superficie.
Donde:
Página 48 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Se calcularon los valores de Número de Curva (CN) representativos para cada subcuenca y escenario,
teniendo en cuenta las características de la cobertura vegetal y la capacidad de infiltración del suelo.
Tabla 7-3. Números de curva recomendados por el S.C.S. para una condición de humedad antecedente tipo II
y una pendiente media del terreno del 5%. (SCS, 1972)
Estudio Hidráulico Cuenca El Morro - Informe Final.docx
Página 49 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
El Servicio de Conservación de Recursos Naturales de EEUU (NRCS), clasifica los suelos en cuatro
grupos hidrológicos de acuerdo con sus características de infiltración:
- A: (Bajo potencial de escorrentía). El suelo presenta una tasa de infiltración elevada incluso en
condiciones de humedad antecedente altas. En general se trata de mantos profundos
excesivamente drenados. Arenas, gravas y loess profundos.
- B: estos suelos presentan una tasa de infiltración moderada, en condiciones saturadas. Se trata
de profundos o bastante profundos de sedimentos con texturas moderadamente finas a
moderadamente gruesas. Presentan una tasa de infiltración moderada. Loess poco profundos,
marga arenosa.
- C: Suelos que tienen una tasa de infiltración lenta incluso en condiciones saturadas. En general
presentan una capa subyacente que impide la infiltración vertical del agua o tienen una textura
fina a muy fina. Marga arcillosa, marga arenosa poco profunda, suelos de bajo contenido
orgánico y suelos generalmente con alto contenido de arcilla.
- D: (Suelos con alto potencial de escorrentía). Estos suelos se caracterizan por presentar una
tasa de infiltración muy baja incluso en condiciones saturadas. Suelos que aumentan de
volumen cuando están mojados, arcillas plásticas pesadas y algunos suelos salinos.
Por su parte, el SCS considera tres condiciones posibles de humedad antecedente señaladas:
Tabla 7-4. Condición de humedad antecedente en función del estado de desarrollo fenológico de la cobertura
vegetal y la lluvia acumulada de cinco días. (S.C.S, 1972)
Teniendo en cuenta la alta capacidad de infiltración de los suelos existentes (grupo hidrológico A) y en
función de los valores recomendados por el Soil Conservation Service se adoptó un CN=65 para
cultivos anuales de secano, CN=45 para coberturas tipo pasturas/pastizales y CN=35 para monte. Vale
señalar que dichos valores son representativos para una condición de humedad antecedente tipo II y
una pendiente media de la cuenca del 5%, con lo cual requieren un ajuste posterior para ser aplicados
en forma directa en los cálculos.
En este sentido, se efectuó una corrección de los valores de CN ponderados para cada subcuenca a
través de la siguiente expresión propuesta por Williams (1995):
(CN 3 CN 2 )
CN 2 s 1 2. exp( 13.86.slp) CN 2
3
Donde:
CN2s es el número de curva ajustado por pendiente para el estado antecedente II.
CN3 es el número de curva sin ajuste para el estado antecedente III.
CN2 es el número de curva sin ajuste para el estado antecedente II.
Página 50 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
El ajuste del CN en función del estado de humedad antecedente fue efectuado a partir de las
expresiones recomendadas por Chow et al. 1988:
4.2CN II
CN I
10 0.058CN II
23CN II
CN III
10 0.13CN II
Finalmente, teniendo en cuenta la distribución espacial de cobertura vegetal por periodo, se definen
los CN para cada escenario y subcuenca:
Tabla 7-5. Números de curva (CN) estimados por subcuenca, periodo y condición de humedad antecedente.
Página 51 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 7-3. Evolución del número de curva (CN) por periodo, para un estado de humedad antecedente tipo II.
El análisis estadístico de las precipitaciones diarias máximas anuales fue aproximado mediante la
expresión empírica de Weibull. El procedimiento consiste en ordenar los datos en forma decreciente y
asociar a cada valor un período de retorno dado por la siguiente expresión:
( N 1)
Tr
m
Donde:
Figura 7-4. Serie anual de precipitación máxima diaria durante el periodo 1968-2017.
Cabe aclarar que por precipitación diaria máxima anual se entiende a la máxima lámina registrada en
veinticuatro horas en un determinado año. Estas se obtienen de los valores registrados por estaciones
automáticas o de las planillas mensuales confeccionadas por los observadores de los pluviómetros.
Dicha información fue sistematizada a los fines de proceder al tratamiento estadístico de valores
extremos mediante las fórmulas de distribución de probabilidad de Gumbel, Pearson III y el Método
Lognormal, con el objeto de determinar las precipitaciones máximas probables correspondientes a
distintos períodos de retorno.
El ajuste analítico de los datos mediante el modelo probabilístico de Gumbel permite obtener el valor
de precipitación máxima diaria (PMD) asociado a un período de retorno determinado, a partir de la
siguiente expresión deducida de la función de distribución acumulada:
Página 52 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
1
ln ln 1
Tr
PT
a
Donde:
PT= precipitación diaria máxima anual (en mm) correspondiente a un período de retorno (Tr).
Tr= período de retorno en años.
m; a= parámetros de la distribución de Gumbel que se calculan mediante las fórmulas siguientes:
a Sn / S x
X m (S x Yn / S n )
( X i X m )2
Sx
2
( N 1)
Donde:
Los parámetros de la variable reducida (Yn; Sn) son sólo función del tamaño de la muestra o extensión
del registro de datos (N). Se obtienen de un cuadro provisto por los manuales en los que figura la
distribución de Gumbel.
Página 53 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
1 x
x
f x
1
e
para x
( )
La función de variable aleatoria x puede expresarse como:
1 x
x
F x Px x1
1
e
.dx
( )
Donde:
2 .
x1
2
2
2
Cs
S
x S.
n.i 1 ( xi x) 3
n
CS
(n 1).(n 2).S 3
Página 54 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 7-6. Ajuste de bondad gráfico de PMD estimada a partir de la distribución tipo Pearson III.
El método Lognormal se basa en que dada una variable aleatoria x, si la función y= ln(x) tiene una
distribución normal, se dice que los valores de x presentan una distribución del tipo logarítmico-normal.
La expresión que define la distribución lognormal de los eventos de la muestra de datos es la siguiente:
1 2
x
f masa x
1
e 2
2
2
A su vez, la expresión que determina la función de probabilidad acumulada lognormal se obtiene de la
integración de la anterior en el rango de interés:
1 2
x
f acum xi f masa x dx
xi xi 1
e 2
dx
2
0 0 2
Donde:
ƒmasa (x) = función masa o de densidad de la distribución lognormal de los eventos de la muestra.
ƒacum. (xi) = función de distribución de probabilidad acumulada lognormal para el evento xi.
x = ln (Pmáx)= logaritmo natural de la precipitación diaria máxima anual.
μ= promedio de los logaritmos naturales de la serie de precipitaciones diarias máximas anuales.
σ= desviación estándar de los logaritmos naturales de la base de datos.
e = base de los logaritmos naturales o neperianos (e = 2,718).
Página 55 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 7-7. Ajuste de bondad gráfico de la PMD estimada a partir de la distribución LogNormal.
Al construir la curva de frecuencias de una serie de eventos máximos anuales se observa, en algunos
casos, que los mayores valores de la variable (principalmente el mayor) se apartan de la tendencia
indicada por el resto de la serie. Ello se atribuye a que durante el período de observaciones
probablemente se hayan producido eventos cuya magnitud corresponda a intervalos de recurrencia
diferentes (generalmente mayores) al del número de años de registro.
Este apartamiento está definido por dos curvas de control a ambos lados de la curva ajustada, que
determinan los límites de confianza dentro de los cuales las observaciones se consideran aceptables.
Si x1 es la magnitud del evento que corresponde al tiempo de retorno R (x1), los límites de confianza
inferior y superior serán:
LSx1 x1 x1
Lix1 x1 x1
Donde Δx1 se obtiene en función de: el nivel de confianza adoptado, el factor de frecuencia, el desvío
standard de la serie y el número de años de registro:
Página 56 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
S
x1 t ( ) . . 1 1,14 K 1,1K 2
N
Donde:
K = factor de frecuencia
α = 60 % t (α) = 1,000
α = 80 % t (α) = 1,282
α = 90 % t (α) = 1,645
α = 95 % t (α) = 1,960
Tabla 7-6. Resultados del test de bondad analítico por el método de Chi-Cuadrado.
Teniendo en cuenta lo anterior, la función LogNormal resultó seleccionada como mejor distribución de
ajuste de valores extremos para la precipitación máxima diaria.
Figura 7-8. Valores de PMD (mm) en función del periodo de retorno, según ajuste tipo LogNormal.
Página 57 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Página 58 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Debido a la falta de registros pluviográficos con extensión suficiente para su tratamiento estadístico,
las curvas Intensidad-Duración-Frecuencia (IDF) empleadas para la determinación de caudales
máximos de crecidas fueron estimadas a partir de las siguientes metodologías semiempíricas:
Según la publicación realizada por Dirección Nacional de Vialidad, la intensidad correspondiente a una
recurrencia de 25 años (I25) puede calcularse mediante la siguiente expresión:
a
I 25
t c
b
Donde:
Para determinar los valores de los coeficientes “a”, “b” y “c” Rühle sugiere las siguientes expresiones:
a = 31 * IH + 0,023 * IH 2,295
IH = 60 mm/hora (válida para la zona en estudio según el mapa de isohietas provisto en la publicación
del Ing. Federico Rühle)
a = 2.137
b = 0,82
c = 6,9
Para duraciones de lluvias superiores a 2 horas, el autor recomienda calcular directamente I25
mediante la siguiente expresión simplificada:
I25 = a / tb
50 T
I I 25 * (1 0,44 * * log )
(50 I 25 )
0.5
25
Página 59 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 7-12. Intensidad media de precipitación que cae durante una hora, para un intervalo de recurrencia de
25 años, Dirección Nacional de Vialidad (1966).
7.3.2.2 Curvas I-D-F para el pedemonte del Gran Mendoza (INA-CRA, 2008)
En base al trabajo “Tormenta de Proyecto en Base a Mediciones en el Pedemonte del Gran Mendoza”
efectuado por el Centro Regional Andino (C.R.A.) del Instituto Nacional del Agua (INA) en el año 1997,
se proponen curvas Intensidad-Duración-Frecuencia (I-D-F) para períodos de retorno de 5, 10, 25, 50
100 y 200 años, determinadas sobre la base de las lluvias que ocasionaban escurrimiento superficial.
El trabajo fue modificado posteriormente en el año 2008, considerando todas las precipitaciones
(generen o no derrame). De esta manera, se propone la siguiente expresión analítica:
(Tr )
Id , Tr
(d 0,268) 0,883
Donde:
Página 60 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Tr 2 5 10 25 50 100 200
ω (Tr) 36,049 53,303 64,727 79,161 89,869 100,498 111,088
Tabla 7-7. Coeficiente ω (Tr) en función del periodo de retorno (INA CRA, 2008).
Resultados obtenidos en experiencias realizadas en diversas partes del mundo, señalan que los
valores medios de las relaciones entre la precipitación máxima probable diaria y las correspondientes
a lluvias de duraciones inferiores a 24 horas se asemejan mucho entre sí, en particular en lluvias de
corta duración (d<120 minutos). Ello se debe a que la gran mayoría de las lluvias intensas están
originadas por células convectivas de similares características en distintos sitios del planeta (Bell,
1969).
Existe una amplia gama de coeficientes sugeridos por diversos autores, con valores muy similares
para iguales duraciones. Caamaño y Dasso (2003) proponen a partir de diferentes técnicas de
estimación, los siguientes valores promedios de las relaciones entre lluvias de duraciones inferiores a
24 horas y la diaria.
Tabla 7-8. Coeficientes medios en función de la duración del evento (Caamaño y Dasso, 2003).
Siendo:
Vale señalar que los valores anteriores indican una variación logarítmica de los coeficientes en función
del tiempo de duración de las lluvias. En este sentido, en los primeros 15 minutos precipitaría el 31,5%
de la PMP diaria, en 60 minutos el 54,3%, en 120 minutos el 63,4% y así sucesivamente.
El modelo de predicción de lluvias máximas DIT estima la relación i-d-T entre intensidad (i) de lluvia,
su duración (d) y el período de retorno (T) a partir de series pluviográficas con distribución Log-normal.
Se basa en una estimación algebraica del factor de frecuencia normal e incorpora la duración de la
lluvia en forma analítica, dándole sentido conceptual a sus parámetros y permitiendo su transposición
a estaciones pluviométricas.
Página 61 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
ln id , Tr A. y B. y C
Donde:
id,T: intensidad media máxima de lluvia en mm/h, para una duración d en minutos, esperable en el
período de retorno T en años.
y 2,584458.(ln T )3 / 8 2,252573
y (ln d) q
La ecuación propuesta presenta tres términos independientes entre sí, que provienen de la hipótesis
de independencia entre la duración explicitada en δy y la recurrencia incluida en Φy. Los coeficientes
A, B, C y el exponente q, son los cuatro parámetros resultantes del modelo; A y C incorporan las
características de la lluvia local, mientras que B y q son parámetros zonales.
Estudios antecedentes (Rico et al. 2010, 2011) han evaluado el comportamiento regional de los
parámetros del modelo DIT, a partir de la calibración e interpolación espacial sobre matrices i-d-T
generadas en distintas localidades de la región Centro-Norte de nuestro país. El presente estudio
propone los resultados de dicha investigación como referencia, con el objeto de evaluar las curvas i-d-
T en las inmediaciones de la cuenca del Morro.
El trabajo “Regionalización de lluvias intensas en Argentina” (Devoto, 2002), presenta una metodología
para la estimación de las curvas I-D-F en localidades sin información o con escasez de datos. El
método se basa en la función de distribución de valores extremos de Gumbel. De esta manera, una
variable aleatoria X caracterizada por la función de Gumbel queda expresada como una función doble
logarítmica del período de retorno (Tr), de la media (Xm) y del desvío estándar (σx). La variable X puede
ser expresada por la siguiente formulación:
Con el objeto de caracterizar la variable X para una determinada duración del evento, se definen
valores zonales de la media (X m) y del coeficiente de variación (Cv). A partir de datos pluviométricos
recopilados y analizados por diversos autores correspondientes a estaciones situadas en distintas
zonas del territorio nacional, se proponen planos del país con las isolíneas de las medias y de
coeficientes de variación para duraciones e tormenta de 1, 3, 6 y 12 horas.
Página 62 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 7-13. Isolíneas de Coeficiente de Variación para una duración de 12 horas (Devoto, 2002).
Para la expresión matemática de las curvas I-D-F, Devoto propone dos modelos de tipo hiperbólico,
uno con dos y otro con tres parámetros (aunque en este último uno de ellos es prácticamente
constante).
i = α / tb
α = 60 * x1 / t1 (1-b)
i = a / (tb + c)
a = i1 * (t1b + c)
b = 0,80
Página 63 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
7.3.2.6 Curvas Intensidad – Duración – Recurrencia para la cuenca del Morro (IMASL, 2018)
Las curvas IDF formuladas para la cuenca del Morro se lograron a partir de la media aritmética entre
los resultados obtenidos a partir de cinco metodologías semiempíricas de referencia implementadas
en nuestro país.
Vale señalar que la duración mínima de evento analizada a partir de las curvas propuestas corresponde
a 20 minutos. Esto se explica puesto que las curvas IDF propuestas por Devoto (2002) presentan
intensidades mayores que la media obtenida a partir de las demás metodologías, con lo cual se
recomienda un análisis de mayor detalle para eventos intensos de corta duración.
Figura 7-14. Curvas IDF obtenidas por cinco metodologías semiempíricas (duración<180min).
Figura 7-15. Curvas IDF obtenidas por cinco metodologías semiempíricas (duración>180min).
Estudio Hidráulico Cuenca El Morro - Informe Final.docx
Página 64 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 7-16. Curvas IDF propuestas para la cuenca del Morro (duración<180 minutos).
Figura 7-17. Curvas IDF propuestas para la cuenca del Morro (duración>180 minutos).
Página 65 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Tabla 7-9. Intensidades máximas probables propuestas para la cuenca del Morro.
R 0,162
I 854,20 0,773
d
Donde:
Vale señalar que esta formulación se ajusta para recurrencias de 10 a 200 años y duraciones de
eventos entre 60 minutos y 12 horas.
Por otra parte, para el diseño de obras menores de drenaje y desagüe pluvial se propone la siguiente
expresión, válida para recurrencias de 2 a 5 años y duraciones de eventos entre 20 minutos y 3 horas.
R 0,335
I 394,30 0,677
d
Página 66 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
La predicción del caudal de proyecto en base al análisis de lluvias de diseño requiere definir el valor
medio de ésta sobre la superficie de aporte. La estimación de la precipitación media sobre una cuenca
implica una extrapolación no trivial, dado que los registros puntuales responden a muestras limitadas
de la variable. Asimismo, en cuencas susceptibles a ser afectadas en su totalidad por un evento pluvial,
tener un puesto pluviográfico resulta una circunstancia poco habitual. Esto explica en parte el menor
desarrollo del análisis de frecuencia areal con respecto al de lluvia puntual y la necesidad de
implementar técnicas para transponer la lluvia de diseño.
Se admite que las láminas locales, extraídas de una función intensidad-duración-recurrencia (IDR)
representan la precipitación sobre áreas reducidas en torno al pluviógrafo, cuyas variables se
encuentran sujetas a las características climáticas y topográficas de la región. Para superficies
mayores, dado que el máximo no se produce simultáneamente en todo el sistema, la precipitación
media difiere del dato puntual y debe ser inferida a partir de este valor.
Es habitual evaluar la precipitación media en torno a una estación núcleo, asumida como foco del
evento, aplicando un algoritmo de atenuación el cual es función del área de aporte y de la duración del
evento. Sin embargo, esta técnica resulta objetable puesto que al calibrar las expresiones, la estación
adoptada puede no ser el centro de todos los eventos, los máximos areales no coinciden
necesariamente con los locales y el proceso meteorológico que produce éstos puede diferir del
causante de las grandes crecidas, que en definitiva interesan.
La porción de lámina a utilizar depende de la duración de las lluvias, ya que las que más perduran
localmente (debido al mecanismo propio de precipitación, su superior cobertura areal o desplazamiento
más lento) tienen mayor probabilidad de superponerse en el tiempo con otras en su entorno. Por otra
parte, investigaciones sobre lluvias extremas en diferentes partes del mundo (Occhipinti, 1989) señalan
que la altura pluviométrica máxima varía con la inversa del cuadrado del área afectada. En este sentido,
existen diversos algoritmos y ábacos de aplicación para evaluar el factor de decaimiento areal, los
cuales indican el porcentaje de lluvia local a tomar como promedio sobre la cuenca.
Para el desarrollo de este estudio se propone implementar una expresión analítica de decaimiento-
área-duración obtenida a partir del modelo CoDA, desarrollado por Catalini (2001) y cuyos parámetros
corresponden a la recalibración efectuada por Caamaño Nelli, Rodríguez y López (2004), con una
batería de 10 pluviómetros de la Red Telemétrica del Centro de la Región Semiárida del Instituto
Nacional del Agua (INA-CIRSA).
CoDA asume que la lámina local caracteriza un área de 25 km² en torno al pluviómetro, donde el
coeficiente de decaimiento areal es función de la superficie de la cuenca según la siguiente expresión:
D 25 k .S k
Donde:
k 2,1438.d 0.4774
Página 67 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Estudios antecedentes (Caamaño Nelli et al. 2012) señalan que este enfoque conduce a la
sobrevaloración del factor de decaimiento areal cuando se considera el área total del sistema. Sin
embargo, eso no implica que el enfoque de las expresiones propuestas sea intrínsecamente malo. Su
falla principal, en especial para lluvias breves, es que la atenuación se evalúa a partir de máximos
puntuales, ya que éstos se deben frecuentemente a fenómenos convectivos convencionales que, al
abarcar extensiones reducidas, conducen a decaimientos desmesurados al considerar el área total del
sistema.
De esta manera, la valoración del decaimiento a partir de un valor puntual responde a un camino
forzado debido a la ausencia de información pluviográfica distribuida. Los intervalos de máxima anual
que interesan son los de la precipitación media areal, los cuales dan lugar a crecidas de mayor
magnitud. Éstos pueden provenir de lluvias frontales o de procesos convectivos de mesoescala, que
cubren sectores grandes o la totalidad de la cuenca con intensidades semejantes, lo cual argumenta
que el abatimiento real sea considerablemente menor. Sin embargo, esta metodología exige una buena
cobertura pluviométrica y series simultáneas largas, situación poco corriente. Asimismo, la
extrapolación es raramente aplicable, dado que sus resultados dependen del área de la cuenca donde
se derivaron, pues la atenuación es función de ésta.
Teniendo en cuenta estas consideraciones, se propone evaluar la variación del parámetro k en función
de la recurrencia o del factor de frecuencia normal incluido en el modelo D.I.T.
k 0,03115. y 0,04824
Donde:
y 2,584458.(ln T ) 3 / 8 2,252573
Finalmente, el replanteo de la atenuación con otro enfoque conduce a una metodología extrapolable
regionalmente, la cual aporta valores de reducción más creíbles. El coeficiente de decaimiento resulta
de esta manera función de la superficie de la cuenca y de la recurrencia de la lluvia.
Página 68 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
El escurrimiento superficial sobre cuencas hidrográficas constituye una amenaza crítica para el medio
físico, dado que favorece procesos de anegamiento, erosión hídrica, transporte de sedimentos y
contaminantes. Los caudales de escorrentía pueden cuantificarse mediante la medición directa del
flujo en estaciones de aforo. Sin embargo, este tipo de prácticas implican un costo económico elevado,
dado que demanda inversiones en infraestructura y recursos humanos. En este sentido, una estrategia
alternativa viable para disminuir costos y tiempo resulta en la predicción de escorrentías a partir de
modelos de simulación.
Una herramienta útil consiste en la aplicación de modelos semi-empíricos como el Hidrograma Unitario
desarrollado por el Servicio de Conservación de Suelos (SCS) del Departamento de Agricultura de los
Estados Unidos (Soil Conservation Service, 1972), actualmente bajo el nombre de Natural Resources
Conservation Service (NRCS). Este método, desarrollado inicialmente para pequeñas cuencas
agrícolas, es ampliamente utilizado por el NRCS y muchas otras instituciones nacionales e
internacionales, permitiendo obtener resultados razonablemente aproximados. Asimismo, el
procedimiento requiere un conocimiento moderado del paisaje, con lo cual constituye un recurso
atractivo en términos operativos.
P Ia F Q
Q F
P Ia S
I a .S
Donde:
(P I a ) 2
Q
P Ia S
Página 69 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
25400
S 254
CN
A pesar de las limitaciones señaladas por diversos investigadores, vale señalar que el modelo
propuesto muestra resultados aceptables para cuencas de carácter agrícola (Ponce y Hawkins 1996)
y eventos de lluvia de gran magnitud (>50 mm) (Mishra et al. 2005). Asimismo, resulta importante
destacar que los valores estimados reflejan tendencias promedio respecto a las respuestas de cada
cuenca, dado el espectro posible de variaciones espacio-temporales en las tormentas y de la condición
de humedad antecedente del medio.
Figura 8-1. Hidrograma Unitario estándar propuesto por el Soil Conservation Service (NEH630, 2010)
El hidrograma unitario estándar propuesto por el S.C.S. se define a partir de dos parámetros
principales: tiempo al pico y descarga pico, lo cuales son estimados a partir de relaciones empíricas
en función del área de aporte y el tiempo de concentración de la cuenca. De esta manera, el caudal
pico puede expresarse a partir de la siguiente formulación:
Página 70 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
A.Q
Q p PRF .
Tp
Donde:
El hidrograma unitario original sugerido por el S.C.S. propone un factor de pico (PRF) igual a 484,
asumiendo una distribución triangular del caudal con un Tr igual a 5/3 Tp y un tiempo de base igual a
8/3 Tp. De esta forma, el 37.50% del volumen total de escorrentía se manifiesta en el tramo ascendente
del hidrograma, mientras que el resto se eroga posterior al pico de crecida.
El PRF se encuentra físicamente relacionado con las capacidades de almacenamiento de las cuencas
y puede variar considerablemente en función de las características y extensión del área de estudio
(NRCS, 1972). Estudios antecedentes indican que este parámetro puede variar entre 100 para cuencas
en regiones planas (Capece et al., 1986) hasta 600 en cuencas urbanas y de altas pendientes (NRCS,
1972). Se presenta a continuación una tabla con los valores recomendados por Wanielista et al., 1997
para distintas características del medio físico.
Tabla 8-1. Factores de pico y relación Tr/Tp recomendados para distintos ambientes.
Dado los costos en recursos y tiempo implicados en el desarrollo de modelos matemáticos más
avanzados, los resultados informados se consideran válidos para cuantificar en forma conservadora el
orden de magnitud de los caudales pico asociados a eventos extremos de escorrentía, valorar su
sensibilidad ante la evolución de la cobertura vegetal y servirán de base para la toma de decisiones
respecto a la planificación de infraestructura en el territorio.
Estudio Hidráulico Cuenca El Morro - Informe Final.docx
Página 71 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 8-2. Variación esquemática del hidrograma triangular para distintos PRF (Wanielista y Yousef, 1993)
Tabla 8-2. Caudales máximos probables (m 3/s) estimados por subcuenca y periodo, asociados a un potencial
normal de escorrentía (condición II de humedad antecedente según S.C.S).
Página 72 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Tabla 8-3. Caudales máximos probables (m3/s) estimados por subcuenca y periodo, asociados a distintos potenciales de escorrentía
(condición I, II y III de humedad antecedente según S.C.S)
Página 73 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 8-3. Caudales máximos probables por periodo para la cuenca del Quebrachal. Potencial normal de
escorrentía (condición II de humedad antecedente según S.C.S).
Figura 8-4. Caudales máximos probables por periodo para la cuenca del Río Nuevo. Potencial normal de
escorrentía (condición II de humedad antecedente según S.C.S).
Página 74 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 8-5. Caudales máximos probables por periodo para la cuenca La Guardia. Potencial normal de
escorrentía (condición II de humedad antecedente según S.C.S).
Figura 8-6. Caudales máximos probables por periodo para la cuenca Zanjón del Cerro Negro. Potencial normal
de escorrentía (condición II de humedad antecedente según S.C.S).
Página 75 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 8-7. Caudales máximos probables en función de la condición de humedad antecedente para la cuenca
del Quebrachal (Recurrencia: 25 años).
Figura 8-8. Caudales máximos probables en función de la condición de humedad antecedente para la cuenca
del Río Nuevo (Recurrencia: 25 años).
Página 76 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 8-9. Caudales máximos probables en función de la condición de humedad antecedente para la cuenca
La Guardia (Recurrencia: 25 años).
Figura 8-10. Caudales máximos probables en función de la condición de humedad antecedente para la cuenca
Zanjón del Cerro Negro (Recurrencia: 25 años).
Página 77 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Al momento de definir un escenario de diseño deberán evaluarse los riesgos que son posibles admitir
en función de las características de la obra, con el objeto de obtener una ecuación económica
aceptable. A modo de ejemplo se indican a continuación algunos tiempos de retorno utilizados
comúnmente para proyectos hidráulicos:
Teniendo en consideración los sucesos acontecidos sobre la cuenca del Morro en las últimas dos
décadas en relación a los daños sobre la infraestructura vial, se detallan a continuación los criterios
adoptados por la Dirección Nacional de Vialidad (DNV) para el diseño de obras de drenaje transversal
(puentes y alcantarillas).
En el caso de Rutas Bidireccionales se deberá adoptar como recurrencia de diseño de las alcantarillas
transversales 25 años con una revancha del 20% del tirante de operación de la misma (funcionando al
80 % de su capacidad) y verificación para 50 años bajo el criterio que se produzca escurriendo
libremente por la estructura, o que el desnivel entre el hombro de talud aguas arriba y la cota del pelo
de agua para dicha recurrencia sea igual o mayor a 0.60 metros.
Con respecto a los Puentes el tiempo de recurrencia de diseño será de 100 años, para el cual la
revancha generada entre el fondo de viga y la cota de máxima de creciente será de 1 m, mientras que
el tiempo de recurrencia de verificación será de 200 años, de tal manera que el puente trabaje a
superficie libre.
Sin perjuicio de lo anterior la Dirección Nacional de Vialidad podrá disponer la verificación y aplicación
de otros tiempos de recurrencia según su juicio para casos particulares. En términos del estado de
humedad antecedente, si bien no existe un criterio explícito en las normativas vigentes, se acepta la
definición de caudales máximos de crecida a partir de un potencial normal de escorrentía (condición II
según el S.C.S.).
Página 78 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Vale señalar que la influencia del aumento en los niveles freáticos sobre la capacidad de abstracción
inicial del sistema y su incidencia sobre la condición de humedad antecedente es aún materia de
estudio, requiriendo el desarrollo de modelos matemáticos más avanzados para valorar correctamente
la transformación lluvia-escorrentía en zonas anegadas o con nivel freático a poca profundidad.
Como ejemplo, se muestra a continuación la evolución del nivel freático registrado en parcelas de
pastura y humedales con distintas profundidades del nivel estático. Como puede observarse, se
identifica una respuesta más sensible en relación a la lluvia antecedente en sectores con la napa
cercana a la superficie. Este fenómeno favorecería la transformación lluvia-escorrentía, disminuyendo
la capacidad de abstracción del medio e incrementando así los excedentes potenciales a ser captados
por la red hidrográfica.
Figura 8-11. Dinámica del nivel freático registrada bajo parcelas de pastura y humedales.
Página 79 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Con el objeto de verificar los valores propuestos por el método del S.C.S. y analizar la sensibilidad de
los caudales de diseño ante el incremento de los niveles freáticos y la conectividad del sistema
hidrográfico, se realizaron simulaciones matemáticas a partir del software HEC-HMS, desarrollado por
el Centro de Ingeniería Hidrológica (Hydrologic Engineering Center) del cuerpo de ingenieros de la
armada de los EE.UU. (US Army Corps of Engineers)
La estimación de la escorrentía efectiva fue realizada a partir del modelo de parámetros concentrados
de pérdidas por infiltración del SCS (método CN), con una pérdida inicial del 20% de la abstracción
potencial máxima (S). El hidrograma de escorrentía directa es estimado por convolución y propagación
mediante el concepto de cascada de embalses lineales (Analogía de Nash), con Tlag=0.6 Tc.
Con el objeto de evaluar los valores estimados en la sección 8, se verificaron todos los escenarios
propuestos considerando distintas distribuciones temporales de precipitación. Los caudales pico
obtenidos presentaron diferencias del 5-10% respecto a los calculados por el método del S.C.S,
resultando en una revisión satisfactoria.
Figura 9-1. Hidrograma de crecida obtenido para un evento de lluvia de intensidad constante, R=100 años y
condición II de humedad antecedente. Cuenca del Rio Nuevo. (Qp=70 m3/s)
Estudio Hidráulico Cuenca El Morro - Informe Final.docx
Página 80 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 9-2. Hidrograma de crecida obtenido para un evento de lluvia de intensidad constante, R=25 años y
condición II de humedad antecedente. Cuenca del Rio Nuevo. (Qp=40 m 3/s)
Figura 9-3. Histograma de diseño para la cuenca del Río Nuevo (R=100 años y d=855 minutos)
Página 81 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 9-4. Distribución de intensidades para la cuenca del Río Nuevo (R=100 años y d=855 minutos)
Figura 9-5. Hidrograma de crecida obtenido para un evento de lluvia con distribución temporal por método de
bloques alternos, R=100 años y condición II de humedad antecedente. Cuenca del Rio Nuevo. (Qp=70 m3/s)
Página 82 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 9-6. Hidrograma de crecida obtenido para un evento de lluvia con distribución temporal por método de
bloques alternos, R=25 años y condición II de humedad antecedente. Cuenca del Rio Nuevo. (Qp=40 m3/s)
Figura 9-7. Comparación de hidrogramas de crecida obtenidos por el método de bloques alternos e intensidad
constante para R=100 años y condición II de humedad antecedente. Cuenca del Río Nuevo.
Página 83 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Desde el punto de vista hidrológico se ha observado un ascenso sostenido de los niveles freáticos
generalizado en la región, el cual avanzó progresivamente en las últimas décadas La causante
principal de este proceso se encontraría dada a partir del cambio en el balance hidrológico que sufre
el territorio, motivado fundamentalmente por el cambio en el uso de suelo y el incremento de las
precipitaciones (ver sección 6.2) Por otra parte, vale resaltar que la configuración hidrogeológica de
las áreas afectadas se caracteriza por la presencia de un basamento rocoso cercano a la superficie,
conformando un acuífero libre de baja potencia y en general variable tanto en sentido norte-sur como
en dirección este-oeste. Dicha variabilidad condicionaría de manera significativa la respuesta de los
niveles freáticos frente a la dinámica hídrica existente.
Con el objeto de valorar la incidencia de este fenómeno sobre los caudales pico de diseño, se propone
considerar a estos sectores como áreas impermeables o con capacidad despreciable de infiltración.
Figura 9-8. Hidrograma de crecida obtenido para un área impermeable del 5%, R=25 años y condición I de
humedad antecedente. Cuenca del Rio Nuevo. (Qp=25 m3/s)
Página 84 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Figura 9-9. Comparación de hidrogramas de crecida para distintos porcentajes de área impermeable.
R=25 años y condición II de humedad antecedente. Cuenca del Rio Nuevo.
Figura 9-10. Comparación de hidrogramas de crecida para distintos porcentajes de área impermeable.
R=25 años y condición I de humedad antecedente. Cuenca del Rio Nuevo.
Los resultados indican que el % de áreas impermeables presenta una incidencia significativa sobre los
caudales pico de crecida. Vale destacar los escenarios que presentan una condición de humedad
antecedente seca, en los cuales la única fuente de aporte de escorrentía corresponde a las áreas de
humedales, manifestando caudales del mismo orden de magnitud que los obtenidos para un escenario
en condición húmeda y sin áreas impermeables.
Página 85 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Con el objeto de evaluar la incidencia de este fenómeno sobre los caudales pico de crecida, se propone
estudiar la respuesta hidrológica de la cuenca a partir del tiempo de concentración. Para valorar el Tc
se implementó la expresión propuesta por el S.C.S (1975), la cual surge como una suma individual de
tiempos de viaje calculados a partir de velocidades medias de escurrimiento.
L
Tc
60.V
Donde:
Se adopta 0.50 m/s para el escurrimiento en flujo concentrado sobre el terreno y un valor de 1.20 m/s
para el flujo en cursos de agua. Ambos valores fueron calculados según recomendaciones y
procedimientos del National Engineering Handbook 630.15 (2010). De esta manera, es posible definir
el tiempo de concentración para distintos estados de desarrollo de la red hidrográfica, expresando el
grado de evolución como la relación entre la longitud de flujo en secciones de encauzadas respecto a
la longitud del recorrido hidrológicamente más distante.
El hietograma evaluado corresponde a una lluvia de R=25 años y duración=860 minutos, con
distribución temporal según el método de bloques alternos. Vale señalar que esta tormenta es
compatible con el escenario de diseño aplicado a la cuenca del Río Nuevo.
Figura 9-11. Comparación de hidrogramas de crecida para distintos porcentajes de desarrollo de la red
hidrográfica (R= 25 años, d=860 min y condición de humedad antecedente tipo II). Cuenca del Rio Nuevo.
Estudio Hidráulico Cuenca El Morro - Informe Final.docx
Página 86 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
10 CONCLUSIONES FINALES
La cuenca hidrográfica El Morro ha sufrido un importante proceso de erosión hídrica en las últimas
décadas, manifestando un crecimiento exponencial de su red de drenaje superficial a partir de la
generación abrupta de nuevos cauces sobre el territorio y su posterior desarrollo. Los impactos
documentados sobre la región son variados y de diversa magnitud. Los mismos incluyen daños
puntuales tales como la obstrucción de rutas Nacionales y Provinciales, fragmentación de áreas
productivas, destrucción de caminos vecinales y daños sobre diversa infraestructura pública y privada.
Por otro lado, se han logrado identificar impactos difusos a mayor escala como el anegamiento de
superficies y la salinización de suelos. Resulta enfático destacar que el proceso descripto es de
carácter reciente en la región y no muestra señales de estabilización, observándose que la densidad
de drenaje ha mantenido un crecimiento exponencial en las últimas tres décadas, y que hasta el
presente no ha cambiado su tendencia.
En función de los severos daños registrados en la región a lo largo de las últimas dos décadas, se
efectuó un análisis histórico de eventos con el objeto de caracterizar las condiciones climáticas
asociadas al desarrollo de estos fenómenos. En función del presente análisis se logró identificar una
concordancia entre la apertura de cauces/explayamientos aluvionales y la manifestación de eventos
extremos de escorrentía. Vale señalar que dichos procesos no solo corresponden a lluvias máximas
históricas, sino que cubren un espectro amplio definido por la condición de humedad antecedente, el
estado de conectividad/consolidación de la red hidrográfica y las características del evento de
precipitación.
Por su parte, se evaluó el impacto del cambio de uso de suelo, el aumento de áreas anegadas y la
evolución de la red hídrica sobre los eventos de crecida a partir de técnicas de modelación matemática
y el estudio de las características hidrológicas y climáticas locales. En este sentido, se pueden
establecer las siguientes conclusiones:
- Procesos de erosión fluvial catalizados por eventos extremos de escorrentía, en relación a las
características de la lluvia (duración e intensidad) y las condiciones de humedad antecedente.
Página 87 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
3. A partir de las modelaciones realizadas, se verifica la incidencia del cambio de uso de suelo sobre
dos aspectos principales:
- En relación al balance hídrico, se observa que el reemplazo de las coberturas nativas (bosque
y pasturas/pastizales) por cultivos anuales de secano (soja y maíz) ha provocado un incremento
de las recargas. Esto se explica dada la mayor capacidad de las coberturas nativas para captar
y evapotranspirar el agua de lluvia gracias a una mayor profundidad de raíces. En este sentido,
vale señalar que el aumento del módulo pluviométrico en el último siglo también ha desarrollado
un papel importante sobre los volúmenes anuales de recarga y escorrentía.
- En términos de los eventos de corta duración, el cambio de cobertura vegetal ha demostrado
un impacto relativo sobre los caudales pico de crecida, en función de la distribución de
cobertura vegetal de cada subcuenca.
4. Se estima que las áreas con napa cercana a la superficie presentan una influencia significativa
sobre los caudales pico de crecida, aún para pequeñas extensiones dentro de la cuenca. Sin
embargo, la influencia del aumento en los niveles freáticos sobre la capacidad de abstracción inicial
del sistema y su incidencia sobre la condición de humedad antecedente es aún materia de estudio,
requiriendo el desarrollo de modelos hidrológicos distribuidos más avanzados.
Con el objeto de definir valores de diseño, se aplicó el método propuesto por el Soil Conservation
Service del U.S.D.A. para la estimación de los caudales máximos probables, el cual ha demostrado
valores aceptables para un amplio espectro de cuencas de carácter agrícola y eventos de lluvia
extremos. Por su parte, todos los escenarios fueron verificados con apoyo del software HEC-HMS
encontrando un buen ajuste entre ambos métodos. Teniendo en consideración que las magnitudes
estimadas reflejan tendencias promedio respecto al potencial de escurrimiento, los resultados
informados se consideran válidos para cuantificar en forma conservadora el orden de magnitud de los
caudales pico asociados a eventos extremos de escorrentía y servirán de base para la toma de
decisiones respecto a la planificación del territorio y tendencias de crecimiento en la infraestructura.
Página 88 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
11 REFERENCIAS BIBLIOGRÁFICAS
Assadi-Langroudi A., Ng´ambi S., Smalley I., 2017. Loess as a collapsible soil: Some basic particle
packing aspects. Quaternary International.
Barbeito O., Beltramone C., Contreras P., 2009. Hidrográfico El Morro, provincia de San Luis.
Activación de procesos fluviales en el sistema. Resúmenes del Taller sobre “Manejo sustentable
del agua en Argentina: agua superficial y subterránea”, provincia de Mendoza.
Barbeito, O., Beltramone, C., Ambrosino S., Contreras, P., Quintana Salvat, F., Moya, G., Barrera,
F., 2008. Estudio geomorfológico de la Cuenca del Morro, Departamento Pedernera, San Luis.
Instituto Nacional del Agua, 56 pp.
Becker, A., 2006. Evaluación del proceso de degradación de suelos por erosión hídrica en una
subcuenca representativa de la región pedemontana del suroeste de la provincia de Córdoba,
Argentina. Tesis Doctoral Facultad de Ciencias Exactas Físico-Químicas y Naturales.
Universidad Nacional de Río Cuarto.
Bell, F.C., 1969. Generalized rainfall duration frequency relationships. Journal of Hydraulic Div., ASCE,
95(1), 311-327.
Blarasin, Cabrera, Matteoda, Frontera, et al., 2014. Recursos hídricos subterráneos en la provincia
de Córdoba. Relatorio del XIX Congreso Geológico Argentino.
Bradshaw, Sodhi, Peh, Brook, 2007. Global evidence that deforestation amplifies flood risk and
severity in the developing world. Global Change Biology, 13:2379-2395.
Brown, Zhang, Mcmahon, Western & Vertessy, 2005. A review of paired catchment studies for
determining changes in water yield resulting from alterations in vegetation. Journal of Hydrology,
310:28-61.
Caamaño Nelli, G.; Rico, A. y Dasso, C. 2012. El Modelo DIT 3p para Predicción de Lluvias Máximas.
Libro de Trabajos del III Taller sobre Regionalización de Precipitaciones Máximas. (ISBN 978-
950-673-953-9) pp 71-82. UNR Editora. Rosario, Argentina.
Catalini, C.G., 2001. Altura Areal de Lluvia para Diseño en la Cuenca del Río San Antonio. Trabajo
Final de Ingeniería Civil, Facultad de Ingeniería, UCC. Córdoba, RA.
Contreras S., Santoni C.S., Jobbágy E.G., 2012. Abrupt watercourse formation in a semiarid
sedimentary landscape of central Argentina: the roles of forest clearing, rainfall variability and
seismic activity. Ecohydrology, Volume 6, Issue 5, Pág. 794–805.
Centro Regional Andino (C.R.A.) del Instituto Nacional del Agua (INA), 1997. Tormenta de
Proyecto en Base a Mediciones en el Pedemonte del Gran Mendoza. Mendoza, Argentina.
Página 89 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Chow, V. T., 1951. A General Formula for Hydrologic Frequency Analysis. Transactions American
Geophysical Union. Vol. 32, Nº 2, pp. 231-237. USA.
Chow, V.T., Maidment, D.R. and Mays, L.W., 1988. Applied hydrology. New York: McGraw-Hill.
Degiovanni, S., N. Doffo y M. Villegas, 2008. Rejuvenecimiento de la red de drenaje del Sur de
Córdoba como consecuencia del cambio climático durante la última centuria. Ed. Especial sobre
Geomorfología y Cambio Climático. Sayago y Collantes (Eds), INGEMA.
Devoto, G.A., 2002. Regionalización de Lluvias Intensas en Argentina. Universidad Católica Argentina.
Facultad de Ciencias Fisicomatemáticas y de Ingeniería.
Dooge, J.C.I., 1973. Linear theory of hydrologic systems. US Department of Agriculture, Technical
Bulletin 1468.
Farley, Jobbágy y Jackson, 2005. Effects of afforestation on water yield: a global synthesis with
implications for policy. Global Change Biology, 11:1565-1576.
Galván M.J., Collado A.D., 2009. Escurrimientos hídricos superficiales en la cuenca hidrográfica de
“El Morro”, provincia de San Luis. Centro Regional La Pampa - San Luis – INTA: provincia de
San Luis.
Galván M.J., 2012. Manejo integrado de los recursos hídricos: cuenca hidrográfica El Morro. Proyecto
de restauración, Centro Regional La Pampa - San Luis – INTA: provincia de San Luis.
Grimaldi, Petroselli, Tauro, Porfiri, 2012. Time of concentration: a paradox in modern hydrology.
Hydrological Sciences Journal, 57 (2).
HEC-HMS User´s Manual – Version 2016. U.S. Army Corps of Engineers. Institute for Water
Resources. Hydrologic Engineering Center (CEIWR-HEC).
Huxman, Wilcox Breshears, Scott; Snyder, et al. 2005. Ecohydrological implications of woody plant
encroachment. Ecology, 86:308-319.
Jackson, Carpenter, Dahm, Mcknight, Naiman, et al. 2001. Water in a changing world. Ecological
applications, 11:1027-1045.
Jobbágy, Nosetto, Santoni, Baldi, 2008. El desafío ecohidrológico de las transiciones entre sistemas
leñosos y herbáceos en la llanura chaco-pampeana. Ecología austral, 18:305-322.
Jobbágy, Nosetto, Villagra, Jackson, 2011. Water subsidies from mountains to deserts: their role in
sustaining groundwater-fed oases in a sandy landscape.Ecological applications, 21:678-694.
Jobbágy E.G., Nosseto M.D., Bernasconi H.O., Colazo J.C., Galván M.J., Mercau J.L., Saenz S.A.,
Colazo E.R., Larrusse C.E., Marchi A.A., Barbosa O.A., Giaccardi A., Hellmers M.M.,
Martínez Álvarez D., Alfonsina Tripaldi, 2015. Los Nuevos Cursos de Agua en la Cuenca El
Morro: Descripción del Proceso y Pautas para su Gestión. Consejo Provincial de Ciencia,
Técnica, Desarrollo e Innovación. Comisión de Trabajo sobre la Cuenca El Morro, provincia de
San Luis.
Página 90 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Kirpich Z. P., 1940. Time of concentration in small agricultural watersheds. Civil Engineering, 10 (6):
362-.
Mishra, S. K., Jain, M. K., Pandey, R. P., and Singh, V. P., 2005. Catchment area-based evaluation
of the AMC-dependent SCS-CN-inspired rainfall-runoff models. J. Hydrological Processes
19(14), 2701–2718.
Mockus, V. 1949. Estimation of total (and peak rates of) surface runoff for individual storms. Exhibit A
in Appendix B, Interim Survey Report, Grand (Neosho) River Watershed, USDA, Washington
DC.
Nosetto, Jobbágy, Brizuela, Jackson, 2012. The hydrologic consequences of land cover change in
central Argentina. Agriculture, Ecosystems and Environment, 154:2-11.
Ponce, V. M. and Hawkins, R. H., 1996. Runoff curve number: Has it reached maturity? J. Hydrol.
Engrg. ASCE 1(1), 11–19.
Rocca R.J., Redolfi E.R., Terzariol R.E., 2006. Características geotécnicas de los loess de Argentina.
Revista Internacional de Desastres Naturales, Accidentes e Infraestructura Civil. pp 149 – 166.
Rühle, Federico G. O., 1966. Determinación del derrame máximo superficial de las cuencas
imbríferas. Dirección Nacional de Vialidad.
Qi X., Xu Q., Liu F., 2017. Analysis of retrogressive loess flowslides in Heifangtai, China. Engineering
Geology.
Sáenz C.A., Rusoci N., Colazo J.C., 2016. Balance hídrico de diferentes escenarios en la cuenca El
Morro. Información Técnica 192. EEA San Luis. ISSN 0327-425X.
Sánchez V. H., Salvioli G., Guimaraes R., Di Chiacchio J. C., 2013. Investigación hidrogeológica en
el área de Colonia Los Manantiales. Convenio Fundación Argentina - Ministerio del Campo:
provincia de San Juan.
Santoni C.S., 2004. Circulación vertical del agua y su relación con la vegetación en zonas áridas y
semiáridas. Tesis doctoral. Grupo de Estudios Ambientales, Instituto de Matemática Aplicada
San Luis (CONICET – Universidad Nacional de San Luis).
Scanlon, Keese, Flint, Gaye, et al, 2006. Global synthesis of groundwater recharge in semiarid and
arid regions. Hydrological Processes, 20:3335-3370.
Shao X., Zhang H., Tan Y., 2017. Collapse behavior and microstructural alteration of remolded loess
under granded wetting tests. Engineering Geology.
Soil Conservation Service (SCS), 1972. National Engineering Handbook, Section 4: Hydrology.
Department of Agriculture, Washington DC, 762 p.
Soil and Water Assessment Tool Theoretical Documentation - Version 2009. Texas A&M
University.
Página 91 de 92
ESTUDIO HIDRÁULICO CUENCA EL MORRO
MODELACIÓN HIDROLÓGICA
INFORME FINAL
Soil and Water Assessment Tool Input/Output Documentation - Version 2012. Texas Water
Resources Institute.
Viviroli y Weingartner, 2004. The hydrological signifance of mountains: From regional to global scale.
Hydrology and Earth System Sciences, 8:1016-1029.
Wanielista, Kersten, Eaglin, 1997. Hydrology: Water Quantity and Quality Control. 2nd Edition, Wiley
and Sons, Inc., New York.
Williams J.R., 1995. The EPIC model. Computer Models of Watershed Hydrology. V.P. Singh (Ed.),
Water Resources Publications, CO.
Xie W.-L, Li P., Vanapalli S.K., Wang J.-D, 2017. Prediction of the wetting-induced collapse behaviour
using the soil-water characteristic curve. Journal of Asian Earth Sciences.
Xu L., Richard Coop M., Zhang M., Wang G., 2017. The mechanics of a saturated silty loess and
implications for landslides. Engineering Geology.
Zhang, Dawes, Walker, 2001. Response of mean annual evapotranspiration to vegetation changes at
catchment scale. Water Resources Research, 37:701-708.
Página 92 de 92