Modflow PDF

Modelación de Flujo Subterráneo.
Aplicación en la cuenca del Milagro-Ecuador y análisis de pozo en

subcuenca del Río Suquía
Facultad de Ciencias Exactas Físicas y Naturales
Universidad Nacional de Córdoba
Maestría en Ciencias de la Ingeniería – Mención en Recursos Hídricos
Modelacion de Flujo Subterraneo.

Aplicacion en la cuenca del Milagro-
Ecuador y analisis de pozo en
Ing. Omar Yagual Muñoz
Director:
Dra. Ing. Teresa María Reyna
Córdoba, 2016
1
Modelación de Flujo Subterráneo. Aplicación en la cuenca del Milagro-Ecuador y análisis de pozo en

Ecuador y analisis de pozo en
Contenido
1 Introducción al Problema ............................................................................................................ 11

1.1 Introducción y Planteamiento del Problema ...................................................................... 11
1.2 Cuenca del Río Milagro ....................................................................................................... 12
1.3 Subcuenca del Río Suquía ................................................................................................... 13
1.4 Objetivos ............................................................................................................................. 14
1.4.1 Objetivo general ........................................................................................... 14
1.4.2 Objetivos Específicos ................................................................................... 15
1.5 Metodología ........................................................................................................................ 15
2 Base Conceptual .......................................................................................................................... 17
2.1 Ciclo del Agua. ..................................................................................................................... 17
2.2 Parámetros hidrogeológicos ............................................................................................... 18
2.3 Hidráulica de Captaciones de Agua ..................................................................................... 20
2.4 Explotación y sobre explotación de acuíferos. .................................................................... 21
3 Modelos Matemáticos de Simulación de Agua Subterránea ...................................................... 23
3.1 Tipos de Modelos ................................................................................................................ 23
3.2 MODFLOW 3-D .................................................................................................................... 24
3.2.1 Características Importantes .......................................................................... 24
3.2.2 Solución de la ecuación de flujo por diferencias finitas. ................................ 25
3.2.3 Datos de Salida MODFLOW ......................................................................... 27
3.2.4 Fortalezas y Debilidades del Programa ........................................................ 27
3.3 Estructura de Modelos Hidrogeológicos ............................................................................. 29
3.4 Etapas en el desarrollo de un modelo de Simulación Hidrogeológico ................................ 30
3.5 Modelos de Simulación ....................................................................................................... 31
3.5.1 Modelos de medio poroso............................................................................. 32
2
3.6 Modelos matemáticos para la simulación del flujo subterráneo en un medio saturado ... 32
3.7 Flujo en sistemas de acuíferos ............................................................................................ 33
3.7.1 FLOWPATH.................................................................................................. 33
3.7.2 MODFLOW ................................................................................................... 34
3.7.3 VISUAL MODFLOW ..................................................................................... 34
3.7.4 FRAC3DVS .................................................................................................. 35
3.7.5 FEFLOW ...................................................................................................... 35
3.7.6 FEMWATER ................................................................................................. 36
3.7.7 ASMWIN....................................................................................................... 36
3.8 Transporte de contaminantes ............................................................................................. 37
3.8.1 MODPATH ................................................................................................... 37
3.8.2 PATH3D ....................................................................................................... 37
3.8.3 MT3D ........................................................................................................... 38
3.9 Validación y Análisis de Sensibilidad. .................................................................................. 38
3.10 Hipótesis de simulación del modelo ................................................................................... 38
4 Caracterización Regional de Cuenca del Milagro y Suquía ......................................................... 40
4.1 Cuenca del Río Suquía. Provincia de Córdoba-Argentina ................................................... 40
4.2 Límites de la Zona de estudio .............................................................................................. 43
4.2.1 Cuenca del Río Suquía-Argentina ................................................................ 43
4.2.2 Evaluación de Antecedentes – Cuenca del Río Suquía ................................ 45
4.2.3 Datos y Antecedentes. .................................................................................. 45
4.2.4 Historia y Problemática Sistema Laguna Mar Chiquita ................................. 46
4.2.5 Clima ............................................................................................................ 48
4.2.6 Temperatura Media Mensual ........................................................................ 51
4.2.7 Precipitación en la cuenca ............................................................................ 52
4.2.8 Evapotranspiración ....................................................................................... 53
4.3 Suelos .................................................................................................................................. 54
4.4 Cobertura del Suelo............................................................................................................. 58
4.5 Consideraciones Hidrológicas ............................................................................................. 60
4.5.1 Balance Hidrológico ...................................................................................... 60
4.6 Aspectos Geológicos e Hidrogeológicos.............................................................................. 61
4.6.1 Formaciones Geológicas .............................................................................. 63
4.6.2 Plano Geológico de Superficie ...................................................................... 65
4.6.3 El agua subterránea en las Sierras ............................................................... 68
4.6.4 El agua subterránea en la zona de piedemonte ............................................ 69
3
4.6.5 Calidad del agua subterránea ....................................................................... 70
4.7 Cuenca del Río Milagro. Milagro-Ecuador........................................................................... 73
4.7.1 Precipitación media anual ............................................................................. 76
4.7.2 Temperatura Media Anual ............................................................................ 76
4.7.3 Evaporación media anual ............................................................................. 77
4.7.4 Evapotranspiración media anual ................................................................... 77
4.8 Geología de la Cuenca del Río Milagro................................................................................ 80
4.8.1 Rasgos estructurales e Historia Tectónica ................................................... 81
5 Trabajos de Campo Sub Cuenca del Milagro............................................................................... 85
5.1 Metodología y Materiales. .................................................................................................. 85
5.1.1 Equipos utilizados. ........................................................................................ 85
5.1.2 Obtención de datos de los instrumentales de Campo ................................... 86
5.1.3 Metodología .................................................................................................. 86
5.2 Modelación Numérica ......................................................................................................... 87
5.3 Calibración del modelo y de parámetros ............................................................................ 99
6 Trabajos de campo Sub Cuenca del Suquía............................................................................... 101
6.1 Metodologìa y Materiales ................................................................................................. 101
6.1.1 Equipos utilizados ....................................................................................... 102
6.2 Afectación del recurso en la zona ..................................................................................... 103
6.3 Factores a monitorear ....................................................................................................... 104
6.3.1 Humedad .................................................................................................... 105
6.3.2 Nivel freático ............................................................................................... 105
6.4 Modelación de transporte de Contaminantes en área circundante Urbanización Cañuelas.
126
7 Análisis de Resultados ............................................................................................................... 128
7.1 Cuenca del Milagro............................................................................................................ 128
7.1.1 Registros obtenidos en los pozos Cuenca del Milagro. ............................... 133
7.1.2 Condiciones iniciales y de frontera ......................................................... 135
7.1.3 Mediciones de Calidad del Agua................................................................. 136
7.2 SubCuenca del Suquìa ....................................................................................................... 137
8 Conclusiones y recomendaciones ............................................................................................. 139
8.1 Introducción ...................................................................................................................... 139
8.1.1 Conclusiones Generales ............................................................................. 139
8.1.2 Conclusiones sobre la Cuenca del Milagro ................................................. 140
8.1.3 Conclusiones sobre el sistema Cuenca del Suquía. ................................... 141
4
8.2 Recomendaciones ............................................................................................................. 141
8.2.1 Cuenca del Milagro ..................................................................................... 141
8.2.2 Cuenca del Suquía ..................................................................................... 142
9 Bibliografía ................................................................................................................................ 144
10 Anexos ................................................................................................................................... 150
Índice de Ilustraciones:
Ilustración 1-1 Cuenca del Milagro y área de estudio Unemi .................................................... 13

Ilustración 2-1 Tipos de Agua disponible en el Mundo. (de Marsily, 1986)............................. 17
Ilustración 3-1 : Sistema acuífero hipotético discretizado (Fuente: (Harbaugh, A.W., 2005)).
............................................................................................................................................................. 25
Ilustración 3-2 Formulación de Bloque centrado (Fuente: (Harbaugh, A.W., 2005)) ........... 26
Ilustración 3-3 Formulación de Punto central (Fuente: (Harbaugh, A.W., 2005)) ................. 26
Ilustración 3-4 .- Visual Modflow Menú Principal ......................................................................... 29
Ilustración 3-6 Etapas en el desarrollo de un modelo Hidrogeológico (Fuente: (Facultad de
Ciencias Fisicas y Matematicas de la Universidad de Chile, 2010) ......................................... 31
Ilustración 3-7 Registro de monitoreo pozo .................................................................................. 38
Ilustración 4-1.-Cuencas Hidrogeológicas – Pcia. de Cba. Fuente (Vasquez, Miatello,
Roque, & otros, 1979) ...................................................................................................................... 41
Ilustración 4-2 Cuenca Hidrográfica Río (Suquía). .................................................................... 43
Ilustración 4-3: Clasificación climática en la llanura chaco pampeana (Fuente: INTA, 1995
citado en Revista “Visión Argentina”. Publicación N°5). ............................................................ 49
Ilustración 4-4 Isohietas medias anuales República Argentina (Fuente: http://inta.gob.ar/).
(Fulginiti F., 2013)............................................................................................................................. 53
Ilustración 4-5 Distribución de Loess y Loessoides (Fuente: (Rocca R. J., 2005), adaptado
de Zárate, 2003). ( (Fulginiti F., 2013)) ......................................................................................... 55
Ilustración 4-6 Suelos en la provincia de Córdoba (Fuente: http://inta.gob.ar/ - Elaborado
en base al Atlas de suelos de la Rep. Argentina). ( (Fulginiti F., 2013)) ................................. 56
Ilustración 4-7 Cobertura del Suelo en la macrocuenca de Mar Chiquita (Fuente:
Elaborado en base a información del SIG del IGN. ( Web www.ign.gob.ar-Instituto
Geografico Nacional de la Republica Argentina). ....................................................................... 59
Ilustración 4-8 Imagen satelital Río Primero – Suquía (Fuente: Elaboración propia - Imagen
tomada de Google Earth). (Fulginiti F., 2013).............................................................................. 60
Ilustración 4-9 Perfil esquemático del subsuelo de la Pampa, desde las montañas
cordobesas hasta el Uruguay (Fuente: según Stappenbeck, 1965 en Vazquez, 1979).
(Fulginiti, 2013) ................................................................................................................................. 63
Ilustración 4-10 Carta Geológica del Segemar. (Fulginiti F., 2013) ........................................ 66
Ilustración 4-11 Plano con la geología superficial de la cuenca (Fuente: elaborado en base
a Cartas Geológicas del Segemar). (Fulginiti F., 2013) ............................................................. 67
Ilustración 4-12 Perfil Hidrogeológico en la zona serrana de la provincia de Córdoba.
(Fuente: Mapa Hidrogeológico de la pcia. de Córdoba – 1978 – Dirección Provincial de
Hidráulica). (Fulginiti F., 2013) ....................................................................................................... 69
Ilustración 4-13.- Cuenca del Río Guayas-Ecuador (CELEC, 2012) ....................................... 74
5
Ilustración 4-14.-Cuencas Hidrográficas del Ecuador-Cuenca 130 Rio Guayas (INAMHI,
2008)................................................................................................................................................... 75
Ilustración 4-15 Precipitación Media Anual Estaciones circundantes. (INAMHI, 2013) ....... 76
Ilustración 4-16 Temperatura Media Anual. (INAMHI, 2013) .................................................... 77
Ilustración 4-17 Evaporación Media Anual. (INAMHI, 2013) ..................................................... 77
Ilustración 4-18 Evapotranspiración Anual ................................................................................... 78
Ilustración 4-19 Cuenca del Rio Milagro (SNI, 2009) ................................................................ 78
Ilustración 4-20 Mapa geomorfológico de la Cuenca del Río Milagro (RemotosCLIRSEN,
2009)................................................................................................................................................... 79
Ilustración 4-21 Subcuencas del Rio (SNI, 2009) ....................................................................... 80
Ilustración 4-22 Mapa Geológico del Ecuador (Ministerio de Recursos Naturales y
Energeticos, 1982) ........................................................................................................................... 83
Ilustración 4-4-23 Perfil Geologico Ecuador (Ministerio de Recursos Naturales y
Energeticos, 1982) ........................................................................................................................... 84
Ilustración 4-4-24 Zona Geologica donde se encuentra la zona de estudio (INAMHI, 2008)
............................................................................................................................................................. 84
Ilustración 5-1 Mapa de Izopiezas Cuenca del Guayas ( Instituto Nacional de Meteorología
e Hidrologìa, 2013) ........................................................................................................................... 88
Ilustración 5-2 Mapa Geomorfologico de la zona de estudio. ................................................... 90
Ilustración 5-3 Datos de Precipitación Diaria (INAMHI, 2013) .................................................. 91
Ilustración 5-4 Líneas equipotenciales escenario actual. ......................................................... 92
Ilustración 5-5 Balance de masas, modelo de masas en escenario actual. ........................... 93
Ilustración 5-6 líneas equipotenciales escenario futuro uno. .................................................. 94
Ilustración 5-7 Balance de Masas escenario futuro uno ............................................................ 95
Ilustración 5-8 Líneas equipotenciales escenario futuro dos..................................................... 95
Ilustración 5-9 Balance de Masa Escenario Futuro dos ............................................................. 96
Ilustración 5-10 Líneas equipotenciales y trayectoria de partículas de escenario futuro tres
............................................................................................................................................................. 96
Ilustración 5-11 Frecuencia Temporal Niveles Calculados....................................................... 97
Ilustración 5-12 Niveles Observados............................................................................................ 97
Ilustración 5-13 Validación Recarga .............................................................................................. 98
Ilustración 5-14 Validación de Conductividad. ............................................................................. 99
Ilustración 6-1 Ubicación del conducto de salida de Aguas Cordobesas en el predio ........ 101
Ilustración 6-2 Laguna al sur del Club House que se encuentra ubicada por encima del
conducto de Aguas Cordobesas en el predio ............................................................................ 102
Ilustración 6-3 Registro de caudales de ingreso a la planta de tratamiento de la
urbanización .................................................................................................................................... 103
Ilustración 6-4 Ubicación de los sistemas de monitoreo .......................................................... 104
Ilustración 6-5 Grafico de monitoreo Pozo Cañuelas Córdoba ............................................... 122
Ilustración 6-6 Histograma Pozo Cañuelas ................................................................................ 124
Ilustración 6-7 Gráfico de datos probabilidad Normal ............................................................. 124
Ilustración 6-8 Análisis de Residuales ........................................................................................ 125
Ilustración 6-9 Curva de regresión ajustada ............................................................................. 125
Ilustración 7-1 Mapa satelital de la ubicación de pozos de estudio (Foto Google Earth) ... 129
Ilustración 7-2 Diagramas Ingreso-Egresos T= 500 días ......................................................... 132
6
Ilustración 7-3 Diagramas Ingreso-Egresos T= 1000 días ...................................................... 132
Ilustración 7-4 Nivel freático Pozo 1 Garita ................................................................................ 133
Ilustración 7-5 Nivel Freático Pozo 2 UNEMI ............................................................................. 134
Ilustración 7-6 Niveles simulados escenario actual .................................................................. 134
Ilustración 7-7 Ubicación Pozo de Observación Loteo Cañuelas Córdoba-Argentina ........ 137
Ilustración 7-8 Zona de Interés Loteo Cañuelas Superficie 60 Ha. ........................................ 138
Ilustración 10-1 Serie de tiempo pluviometría 1981-2014 (INAMHI, 2013)........................... 151
Ilustración 10-2 Datos de Estaciones Meteorológicas cercanas al lugar de estudio. (Instituto
Nacional de Meteorología e Hidrología(INAMHI), 2009) .......................................................... 152
Ilustración 10-3 (INAMHI, 2013)................................................................................................... 152
Ilustración 10-4 Descripción de Mapa Izopiezas (INAMHI, 2013) .......................................... 153
Ilustración 10-5 Ubicación de la Cuenca Hidrográfica del Guayas (INAMHI, 2008)............ 153
Ilustración 10-6 Precipitación plurianual SERIE 1965-1999 (INAMHI, 2008) ...................... 154
Ilustración 10-7 Detalle según colores de los rangos de precipitaciones (INAMHI, 2008) . 154
Ilustración 10-8 Perforación estratigráfica próxima a los pozos de estudio. (B.S.L.C
Consultores Cía. LTda., 2013) ..................................................................................................... 155
Índice de Tablas:
Tabla 4-1 Parámetros climáticos promedio de Córdoba. (SMN, 2013)

(www.freemeteo.com) ...................................................................................................................... 42
Tabla 4-2 Precipitación –Temperatura Anual Años 1961-1990 (SMN, 2013) ........................ 42
Tabla 4-3 Meteorología habitual en CORDOBA OBS. Lat.-31.24 Lon:-64.11. Período 1961-
1990. (Fuente: (www.freemeteo.com)) ......................................................................................... 52
Tabla 4-4 Meteorología habitual en CORDOBA AERO. Lat.-31.24 Lon:-64.11. Período
1961-1990. (Fuente: (www.freemeteo.com)) ............................................................................... 52
Tabla 4-5 Estratigrafía sector occidental de la cuenca en la provincia de Córdoba .............. 62
Tabla 4-6 Estratigrafía sector central y oriental de la cuenca en la provincia de Córdoba.
(Fulginiti F., 2013)............................................................................................................................. 62
Tabla 4-7 Expone las principales características hidrogeológicas a nivel general que
distinguen la zona de sierra de la de piedemonte según (Cioccale M. A., 1998). (Fulginiti F.,
2013)................................................................................................................................................... 70
Tabla 4-8 Valores de Estaciones Meteorológicas cercanas al sitio de Estudio. (INAMHI,
2013)................................................................................................................................................... 76
Tabla 5-1 Pozos Registrados Universidad Estatal de Milagro .................................................. 87
Tabla 5-2 Perfil Estratigráfico de una perforación cercana a los pozos referidos. (Moreno L.,
2013)................................................................................................................................................... 89
Tabla 6-1 Parámetros de calidad del agua de la perforación .................................................. 104
Tabla 7-1 Pozos del Área de Estudio .......................................................................................... 130
Tabla 7-2 Parámetros ajustados y métodos utilizados ............................................................. 131
Tabla 7-3 Geometría y Límites del Área de Estudio Zona 17-UNEMI Milagro ..................... 135
Tabla 7-4 Parámetros de Calidad de Agua de los pozo Garita. ............................................. 136
Tabla 7-5 Parámetros de Calidad de Agua de los pozo Bloque K ......................................... 136
7

Ecuador y estudio de pozo en
Resumen
El conocimiento del flujo subterráneo es parte fundamental de los estudios de
cuencas, en particular si se desea realizar una planificación integral del recurso. En
la actualidad los planificadores y profesionales del área, buscan información
integral de las cuencas subterráneas que les permitan tomar decisiones referidas
no sólo a la disponibilidad del recurso, sino también a la calidad del mismo. Dentro
de este contexto los estudios generales de las cuencas, las mediciones específicas
continuas y los modelos matemáticos permiten mejorar el conocimiento del recurso
hídrico subterráneo. Dentro de este panorama el uso de modelos matemáticos
permite el análisis de diferentes escenarios que favorecen la toma de decisiones
que sean ambientalmente sustentables.
En el presente trabajo, se seleccionó una subcuenca donde se modelaron

escenarios potenciales, como el incremento del consumo, cambio en el uso del
agua subterránea, contaminación del acuífero. Además se planteó un escenario
para remediación. El área de estudio se encuentra en la Subcuenca del Río
Yaguachi, que a su vez tiene las subcuencas del Milagro y Estero Los Monos en la
cual están los predios de la Universidad Estatal de Milagro. Se consideró en la
calibración del modelo, el criterio cuantitativo en el cual, se debe tener en cuenta
los datos de niveles observados y contrastar con los niveles que se simulan, el
parámetro calibrado fue la conductividad hidráulica. Es preciso indicar que la
recarga se estimó con los datos técnicos recolectados de la estación meteorológica
cercana al sitio de estudio.
Y por otra parte se analizaron las fluctuaciones y ventajas del monitoreo de

acuíferos no confinados para esto se consideró la subcuenca, ubicada en un área
situada en una urbanización del Sur de la ciudad de Córdoba, donde se realizó el
monitoreo de los niveles freático a través de instrumental colocado en un pozo de
monitoreo del acuífero no confinado. Además se consideró, una modelación de
transporte de contaminantes, que describe en la zona circundante a la
Urbanización Cañuelas un derrame de combustible desde un camión repartidor de
gasolina, que a través de las cunetas de la vía, llevaron al hidrocarburo hasta una
capa freática.
8
Los pozos monitoreados en Ecuador, demostraron una gran capacidad de

recuperación, luego de estar sometidos a bombeos de largos períodos, en donde
en la perforación estratigráfica hasta los 23 m efectuada en la zona, se muestra
que existe un predominio de limos de baja plasticidad, por lo cual la información
resultante es de gran interés para posteriores trabajos de investigación que
permitan cuidar de mejor manera el recurso hídrico subterráneo.
En Córdoba, el pozo que se analizó, demuestra datos de una menor capacidad de

permeabilidad que los analizados en Milagro, hay que tomar en cuenta que ambos
están en entornos diferentes. Sin embargo ante el escenario del vertido de un
contaminante a su capa freática, presenta poca capacidad de dispersión.
Palabras Claves: agua subterránea, monitoreos, modelos matemáticos de

acuíferos.
Abstract
The Knowledge of the groundwater flow is an essential part of studies of

watersheds, particularly if you want to perform a comprehensive resource planning.
Nowadays planners and practitioners, searching comprehensive information on the
groundwater basins to make informed referring not only to the availability of the
resource, but also the quality of the decisions. Within these context general studies
of watersheds, specific continuous measurements and mathematical models help
improve knowledge of groundwater resources. Within this scenario the use of
mathematical models allows the analysis of different scenarios that favor decisions
that are environmentally sustainable.
A sub-basin where potential scenarios as the increase in consumption, changes in
groundwater use, aquifer pollution were modeled were selected. Also he raised a
remediation scenario. The study area is located in the basin of the Rio Yaguachi,
which in turn has the watersheds of the Miracle and Estero Los Monos in which are
the premises of the State University of Milagro.
It was considered in the model calibration, the quantitative criterion in which should
take into account the observed data and contrast levels to levels that are simulated,
the parameter was calibrated hydraulic conductivity. It should be noted that the
recharge was estimated with the technical data collected from the study site near
the weather station.
On the other hand it analyses fluctuations and advantages of monitoring confined
aquifers us for this sub-basin, located in an area located in a residential area of the
south of the city of Cordoba, where monitoring of groundwater levels was made
through considered analyzed Instrumental placed in a monitoring well unconfined
aquifer.
9
In addition it was considered a modeling contaminant transport, which describes in

the area surrounding the urbanization Cañuelas, a fuel spill from a delivery truck
gasoline, which through the gutters of the road, led the hydrocarbon to a water
table.
Wells monitored in Ecuador, showed great resilience after being subjected to
pumping long periods, where the stratigraphic drilling up to 23 m made in the area
shows that there is a predominance of silt of low plasticity, whereby the resulting
information is of great interest for further research to enable better care for
underground water resources.
In Cordoba, well that was analyzed, shows data from a smaller capacity
permeability than those discussed in Milagro, take into account that both are in
different environments. However at the stage of discharge of a pollutant to the water
table, it has little capacity for dispersion.
10
1 Introducción al Problema
1.1 Introducción y Planteamiento del Problema
Según datos publicados por el Programa Mundial de Evaluación de los Recursos

Hídricos de la UNESCO (WWAP, 2003), en el siglo pasado la demanda mundial de
agua se multiplicó por más de seis mientras que la población del planeta se triplicó.
De no mejorar la gestión de los recursos hídricos y los ecosistemas conexos, en el
2025, dos tercios de la humanidad padecerán problemas de penuria de agua grave
o moderada. Como se desprende de las variadas Conferencias, Cumbres y
Encuentros Internacionales sobre el tema: “El acceso al agua, una fuente de vida
no sustituible, es un derecho humano individual, colectivo y fundamental, y su
conservación es una responsabilidad fundamental colectiva de toda la humanidad.
Todos los seres vivos tienen derecho al agua. El agua es parte de la herencia
común de la humanidad.”
Esta creciente preocupación por el agua a nivel mundial, tanto como elemento de
consumo así como factor integrante de nuestro hábitat, justifica sobremanera el
estudio de la hidrogeología dentro de cualquier zona en que se encuentre presente.
(Fulginiti F., 2013)
Es importante mencionar que los usos del agua son básicamente los siguientes:
 Uso Doméstico
 Uso Agrícola (riego )
 Uso Industrial
El uso agrícola es el que demanda la mayor cantidad de agua. Cuantificar que

porcentaje corresponde a cada uso, a fin de planificar y gestionar el uso adecuado,
será siempre una prioridad, es por esto que este estudio se basará en el análisis
del acuífero no confinado. Como es conocido un modelo, es una representación
simplificada y seleccionada de un sistema real, la que reproduce en forma
aproximada las relaciones causa-efecto de interés (Espinoza C., 2012). Un modelo
de agua subterránea es una representación física, analógica, o matemática que
simplifica la compleja hidrología y química de un sitio seleccionado. Estos modelos
utilizan sistemas de tamaño más reducido para estudiar temas de escurrimiento en
aguas subterráneas, por ejemplo columnas de suelo, analogías con procesos o
fenómenos diferentes para estudiar sistemas hidrogeológicos, las ecuaciones que
gobiernan el flujo y conservación de la masa para simular flujo de agua y transporte
de sustancias contaminantes. Dicho modelo se basa en observaciones del “mundo
real” y percepciones del sitio seleccionado.
En esta propuesta las dos Cuencas que se estudiaron están ubicadas en dos
regiones geográficamente distintas. Una de ellas ubicada en la Provincia del
Guayas, Ecuador y la otra en la Provincia de Córdoba, Argentina.
11
La cuenca del Suquía, está ubicada en la zona Mediterránea Argentina mientras

que la cuenca del Milagro se ubica en la zona litoral, región Costa ecuatoriana.
Sobre estas cuencas existen pluviometrías distintas.
Uno de los mayores problemas que tiene el uso de modelos matemáticos para el
estudio de agua subterránea, es la escasa información que existe sobre la misma,
es por esto que cuando se quiere iniciar un estudio, primero se debe instrumentar
la zona para poder recabar información.
Por otra parte la mayoría de los estudios de perforaciones se realizan a través de

modelos de flujo en pozos en régimen impermanente a través de los resultados de
los ensayos de bombeo de pozos los cuales no permiten considerar toda la
variabilidad que existe en los medios porosos saturados o no saturados.
1.2 Cuenca del Río Milagro

En relación a la cuenca del Río Milagro, la corriente hidrológica más notable en el
área de estudio, es el río del mismo nombre, cuyo origen se sitúa en las cercanías
del poblado de Barraganetal a 110 msnm, el cual recorre 80 kilómetros hasta el
Cantón Milagro. Aguas arriba de este sitio, recibe recarga de los principales ríos
tributarios por la margen derecha, el río Venecia y en la ciudad las aguas del estero
Las Damas, para terminar su recorrido en el Río Yaguachi, a unos 14 kilómetros
aguas debajo de Milagro aproximadamente. Tiene una área 226,74 Km², la zona de
estudio del acuífero descrito queda dentro de los predios de la Universidad Estatal
de Milagro, en donde el mayor uso en esta zona es para riego de áreas verdes y
abastecimiento de agua de las Edificaciones Administrativas y Académicas, que
3
actualmente bombean 814.574,88 m /año, esta cantidad de extracción ha
provocado descensos leves en los niveles de agua subterránea, por lo tanto es
importante predecir el comportamiento del flujo. Las estimaciones de la extracción
de agua por bombeo suelen ser datos muy importantes para el modelado de aguas
subterráneas, especialmente para las predicciones de reducción de niveles
piezométricos. El valor de extracción se tomó como un gasto constante los 365
días del año.
12
Ilustración 1-1 Cuenca del Milagro y área de estudio Unemi
Se planeó la proyección del bombeo durante los próximos 10 años, tomando como
base las condiciones iniciales del año 2013,estas condiciones iniciales se refieren a
los dos pozos con bombeo, que se denominó escenario de modelación uno.
Se planteó otro escenario de simulación, en el cual partiendo de las condiciones

iniciales al año 2013, descrita en el párrafo anterior, se llegue a tener cinco pozos
en total con funcionamiento y gasto constante durante los 365 días del año, que se
denominó modelación dos.
Otro escenario de simulación que se consideró, es el de una pluma de

contaminación debido a un vertido cloacal de un pozo séptico defectuoso, que se
denominó escenario tres.
1.3 Subcuenca del Río Suquía

La Cuenca del río Suquía tiene su inicio en las Sierras Pampeanas de Córdoba.
Una vez que cruza las Sierras Chicas fluye por la planicie y cuando esta próximo al
final de su curso, el cauce principal se divide en varios ramales que según las
oscilaciones temporales del balance de aguas unas se infiltran, otras causan
humedales o desembocan en la laguna.
Este río tiene un sistema discontinuo, en donde existen picos de caudales máximos
en la época de precipitaciones, en cambio en la época de estiaje sus caudales se
reducen de manera importante.
La zona que tiene como plataforma de progreso a las aguas del río que se
menciona en líneas más arriba, también llamado Primero, es comparativamente
pequeña, aunque en su interior se encuentra la Ciudad de Córdoba, segunda en
13
importancia de la República Argentina después de su capital Buenos Aires, con un

importante movimiento industrial consolidado en los últimos años, además de
poseer un polo de desarrollo macro turístico nacional y sin dejar de lado el gran
poder agrícola y ganadero.El área hidrográfica del curso hídrico mencionado es de
aproximadamente 7.500 Km2 (Fulginiti F., 2013), donde una pequeña parte está en
la zona montañosa, es decir, 1900 Km², el resto del área se sitúa en la llanura
hasta su llegada a la Laguna Mar Chiquita.
Existen dos principales ríos que se forman en las montañas Sierras Grandes y
Sierras Chicas que dan inicio al Suquía, estos son el río Cosquín y el San Roque.
El río Cosquín lo forman el río Yuspe y el Grande de Punilla y un poco antes de la
ciudad de Cosquín estos se unen. En el trayecto hasta el embalse San Roque el
Cosquín recibe a los arroyos las Mojarras y Chorrillos. En la ciudad de Villa Carlos
Paz, donde inicia prácticamente el río Suquía su carrera hacia la laguna Mar
Chiquita, recibe al río San Antonio. En la ciudad de Córdoba, alimenta sus aguas
del arroyo La Cañada, famosa vertiente que a su paso por la urbe fue encajonada
en piedra para evitar desastres en sus crecidas. Es importante mencionar que
entre el dique San Roque y la Mar Chiquita, tiene el río Suquía un recorrido de
aproximadamente de 200 km. (Gobierno de la Provincia de Córdoba, 2012).
El estudio se realiza dentro de la cuenca del río Suquía en una subcuenca ubicada
al sur este de la Ciudad de Córdoba fuera de la Avenida de Circunvalación. El
estudio planteó el seguimiento de los niveles piezométricos en un pozo ubicado al
sur este de la ciudad de Córdoba donde el acuífero freático se haya utilizado
especialmente como destino final de efluentes.
1.4 Objetivos
A continuación se plantean los objetivos generales y particulares de la presente

tesis.
1.4.1 Objetivo general
El objetivo general es:

- Evaluar la importancia del conocimiento hidrogeológico para la definición de
manejos ambientalmente sustentables en zonas con posibles riesgos
ambientales.
- Desarrollar estudios sobre la modelación de flujo en acuíferos no confinados
sometidos a explotación a través de pozos.
- Evaluar la modelación del flujo en la cuenca subterránea del Río Milagro,
Ecuador.
- Estudiar la evolución de los niveles freáticos en los predios de la
Urbanización Cañuelas ubicada en la Ciudad de Córdoba, en una
subcuenca del río Suquía parte de la cuenca endorreica de Mar Chiquita,
Argentina.
14
- Conocer características del acuífero en la zona puntual de estudio, así como

los movimientos del agua subterránea y el medio que la contiene, para
predecir el comportamiento para años futuros, así mismo prevenir zonas de
abatimientos y reducir sus efectos a causa del bombeo en exceso.
1.4.2 Objetivos Específicos
Como objetivos específicos se contemplan:
- El trabajo en la Cuenca del río Milagro permitió estudiar el comportamiento

de un área bajo bombeo ubicada en los predios de la Universidad Estatal de
Milagro
- Y por otra parte se analizaron las fluctuaciones y ventajas del monitoreo de
acuíferos no confinados para esto se consideró un área situada en una
urbanización del Sur de la ciudad de Córdoba, donde se realizó el monitoreo
de los niveles freático a través de instrumental colocado en un pozo de
monitoreo de tipo antes mencionado.
Además se proyectan los siguientes objetivos para poder realizar los estudios:
- Establecer la geometría del área de estudio de ambas cuencas.

- Determinar la calibración a fin de que se reproduzcan las variables de
estado medidas en terreno.
- Evaluar la calidad del agua en cuanto a Conductividad, Salinidad y de
Solidos Disueltos Totales de los pozos.
- Analizar los niveles piezométricos de los pozos estudiados.
- Determinar la dirección del flujo.
- Predecir los niveles para diferentes tiempos.
- Determinar las características del acuífero.
- Usar modelos computacionales exclusivos para el análisis de flujo
subterráneo y transporte de contaminantes.
1.5 Metodología
La metodología para obtener los objetivos propuestos, será la que se presenta a

continuación:
1. Análisis y comprensión de los procesos matemáticos generales en el flujo de

agua subterránea.
2. Monitoreo de pozos, en donde se registre en una serie de tiempo definida, la
variación de niveles y temperatura a lo largo del tiempo.
3. Obtener información pluviométrica de la Estación INAMHI Milagro, de
acuerdo al periodo de monitoreo.
4. Procesamiento de la información obtenida en el monitoreo de pozos y
contrastarla con la conseguida desde los datos de precipitación de lluvias.
15
5. Tomar las mediciones de calidad de agua con el equipo respectivo, tales

como sólidos totales disueltos, conductividad y salinidad.
6. Obtener el caudal de bombeo del equipo instalado, verificar el nivel dinámico
de los pozos.
7. Comparación de los datos necesarios para la entrada y calibración del
modelo referido.
Los puntos anteriormente citados, se resumen en los capítulos que tiene este
trabajo de tesis y son:
1.- Introducción al Problema.
2.- Base Conceptual.
3.- Modelos Matemáticos.
4.- Caracterización Regional de la Cuenca del Río Suquía y Milagro.
5.- Materiales y Métodos para el estudios para la zona sur este de la ciudad de
Córdoba.
6.- Materiales y Métodos Sub Cuenca del Río Milagro.
7.- Análisis de los Resultados.
8.- Conclusiones y Recomendaciones.
Finalmente en el capítulo 9 se encuentra la bibliografía y en el 10 los anexos.
16
2 Base Conceptual
En este capítulo se presentan los aspectos conceptuales vinculados al desarrollo

de la tesis.
2.1 Ciclo del Agua.
El agua es un elemento vital para la vida de los seres vivos, como tal se debe
entender el mecanismo natural por el cual ésta circula, este mecanismo se conoce
como el ciclo del agua.
1400.0 Agua Disponible a nivel mundial

1200.0
1000.0
Millones de Km³
800.0
1320.0
600.0
400.0
200.0
38.2
0.0
agua salada agua dulce
Ilustración 2-1 Tipos de Agua disponible en el Mundo. (de Marsily, 1986)
En la Ilustración 2-1, se aprecia la cantidad de agua salada a nivel mundial, en

océanos, mares y lagunas, y la cantidad de agua dulce en lagos, ríos, arroyos,
atmósfera y bajo la superficie terrestre, es decir, el agua subterránea.
De lo antes mencionado, se nota que el agua salada es una gran mayoría,

mientras que el agua dulce es un porcentaje menor, sin embargo esta es la que se
puede utilizar para la mayoría de usos, sea doméstico, agrícola e industrial. Una
definición del ciclo del agua es aquella en la cual se considera que la cantidad o
masa total de agua es en esencia constante en el tiempo, al hacer esta
consideración nos lleva a tomar en cuenta aunque no sea necesario expresarlo, la
del movimiento de masa de agua de lugar, de estado o ambos a la vez. Los
mencionados movimientos o cambios de estado vienen impulsados por dos fuerzas
motrices principales:
 La energía térmica, que es generada principalmente por la radiación solar,

hace que se desencadene los cambios de estado del líquido (líquido a
gaseoso, sólido a gaseoso o líquido) y de la circulación atmosférica.
17
 La gravedad, que supone el desplazamiento de las masas de agua desde

las cotas más altas a las cotas más bajas. (Fraile, J.; Escuder, R.; Jordana,
S.; Ribera, F.; Sanchez-Vila, X.; Vásquez-Suñé, E., 2009)
2.2 Parámetros hidrogeológicos
Es importante entender cada uno de los parámetros de un acuífero. Para esto es

necesario considerar los acuíferos como sistemas que se caracterizan por un
funcionamiento, que se ajusta o regulariza mediante la recarga y las extracciones,
es decir lo que entra y sale del sistema. Las características o parámetros de un
acuífero permiten definir el funcionamiento de un acuífero frente a unas
determinadas acciones exteriores.
Porosidad (total y efectiva).- La mayor parte de rocas y suelos naturalmente

contienen cierto porcentaje de espacios vacíos, que pueden ser ocupados por agua
o fluidos. En el estudio de la porosidad, se puede distinguir entre(1) la existencia de
vacíos y su geometría adecuada, la definición de la porosidad total y (2) la manera
en que el fluido se distribuye en estos vacíos y las consiguientes relaciones fluido-
sólido, que nos permiten definir la porosidad cinemática o porosidad efectiva. (de
Marsily, 1986)La porosidad total se define por el cociente entre volumen de vacíos
y Volumen del sólido
Volumen de vacíos
Porosidad total
Volumen de sólido
Por su parte la porosidad efectiva o cinemática es el cociente entre volumen de

agua capaz de circular y el volumen total de roca.
Porosidad Volumen de agua capaz de circular

cinemática Volumen total de roca
Conductividad hidráulica.- La conductividad hidráulica K es una constante que

depende del medio poroso o algunas veces denominada coeficiente de
permeabilidad. Es usualmente expresada en metros por segundo. (de Marsily,
1986).
La conductividad es el valor encontrado por Darcy en las experiencias realizadas a

través de filtros de arena para purificar el agua. De estos estudios encontró que la
relación del caudal era linealmente proporcional a la sección y al gradiente
hidráulico. La ley de Darcy expresada en medios saturados o no saturados
unidimensionalmente, tiene en cuenta que la conductividad es función del estado
de humedad como:
18
H
q = -K(  ).
L
donde q es el caudal por unidad de área de flujo o velocidad aparente, K() es la

conductividad hidráulica a humedad del suelo  y :
H = h(  )  z
donde h() es la altura de presión en el suelo con humedad  y z es la distancia

medida positivamente hacia abajo desde la superficie del suelo.
Actualmente al término permeabilidad se deja para identificar la permeabilidad

intrínseca (k).
K.
k=
 .g
donde:
K = conductividad hidráulica
 = viscosidad absoluta del fluído
 = densidad del fluído
g = aceleración de la gravedad
Coeficiente de Transmisividad.- Representa la capacidad del medio acuífero para

transmitir agua en toda su altura, en otras palabras, el caudal que atraviesa una
faja de base unitaria y altura igual del acuífero.
Es el parámetro que nos muestra la facilidad del agua para circular horizontalmente
por una formación geológica
Coeficiente de almacenamiento.- Es el volumen unitario de agua descargado por

un prisma vertical de base unitaria y altura del acuífero, cuando desciende una
unidad de longitud de altura piezométrica media. (Saenz, 1999)
Coeficiente Volumen de agua liberado

Almacenamiento Volumen total que ha bajado la superficie piezométrica
Tanto para determinar el coeficiente de almacenamiento, como el coeficiente de

transmisividad, existen los métodos de Theis y de Jacob, entre los más conocidos.
El coeficiente de almacenamiento es, como la porosidad eficaz, adimensional

(volumen/volumen).
19
Tipos de Acuíferos Valores Típicos

Acuíferos Libres 0,3 a 0.001
Acuíferos Semiconfinados 10-3 a 10-4
Acuíferos Confinados 10-4 a 10-5
Tabla 2-1 Coeficientes de almacenamiento según tipo de acuífero
2.3 Hidráulica de Captaciones de Agua
Una captación consiste en extraer a la superficie un caudal determinado de agua

presente en el subsuelo, de tal forma que, se obtenga un rendimiento satisfactorio
desde el punto de vista de la relación beneficio-costo. (Fraile, J.; Escuder, R.;
Jordana, S.; Ribera, F.; Sanchez-Vila, X.; Vásquez-Suñé, E., 2009)
Las captaciones de aguas subterráneas tienen como fin definido, satisfacer una
demanda de agua en un sitio concreto. Cuando la zona saturada se pone en
contacto con la superficie del suelo, presentándose como un manantial o una
fuente, se dice que hay una captación natural. Las fuentes o manantiales brotan
generalmente en la zona de contacto entre una formación geológica permeable y
otra de baja permeabilidad, tienen además un caudal variable en función de las
condiciones hidrogeológicas existentes en cada período en el acuífero. No siempre
las captaciones naturales, están donde la demanda lo amerita, es aquí que ante
esta necesidad de satisfacer la mencionada demanda, corresponde realizar una
captación artificial.
Estas captaciones artificiales constan de 2 partes, y son las siguientes:
 La obra de acceso al acuífero

 La obra de captación del agua subterránea, propiamente dicha.
Las captaciones artificiales más conocidas son:
 Zanjas
 Galerías
 Pozos
Las zanjas, son excavaciones horizontales de poca profundidad, se establecen en

acuíferos de pequeño espesor, donde el nivel freático está a poca profundidad.
(Fraile, J.; Escuder, R.; Jordana, S.; Ribera, F.; Sanchez-Vila, X.; Vásquez-Suñé,
E., 2009)
Las galerías son extracciones horizontales, se construyen en acuíferos en medios

fracturados, en los que se aprovecha el corte de la mayor cantidad de fracturas, y
en acuíferos en medios no consolidados para poder llegar fácilmente en sentido
lateral al nivel freático. (Benitez, A., 1972).
20
Los pozos son hoyos verticales, que mediante bombas se extrae fácilmente el
agua.
Los objetivos que tiene la hidráulica de captaciones se pueden reunir en tres

grupos:
El problema de simulación. Una vez conocido el medio geológico y las

propiedades hidráulicas del acuífero, el problema de la simulación directa consiste
en obtener los niveles piezométricos en cualquier dominio del acuífero, y en
cualquier tiempo, cuando dicho acuífero se somete a un bombeo, constante o
variable en el tiempo.
El problema de estimación. El problema principal de un estudio de hidráulica de

captaciones es la dificultad de conocer los parámetros hidráulicos. Tratar de
obtener los mencionados parámetros es lo que se conoce con el nombre de
problema inverso. Este problema consiste en que dada la geología, a través de
mapas geológicos, muestras de perforación, etc., si se puede conocer la respuesta
del acuífero, es decir, la variación de los niveles a un registro determinado de
caudales, que puede ser un registro histórico, en una serie de puntos y para una
serie de tiempos, entonces será posible obtener los parámetros que rigen el
sistema acuífero.
Comportamiento de las captaciones a escala local. Constantemente se debe

evaluar la eficiencia y correcta construcción de las captaciones, de la misma forma
su comportamiento a corto y largo plazo. Esto nos permite evaluar y saber cómo
evolucionan los caudales extraíbles con el tiempo, pues, en caso que estos
disminuyan por debajo del valor de interés, se puede tomar las medidas de
corrección oportunas. (Fraile, J.; Escuder, R.; Jordana, S.; Ribera, F.; Sanchez-
Vila, X.; Vásquez-Suñé, E., 2009).
Para obtener los parámetros hidráulicos de un acuífero, se deben realizar ensayos

hidráulicos.
El ensayo hidráulico consiste, en bombear agua desde un pozo y observar la

evolución en el tiempo de los descensos de nivel en dicho pozo de bombeo y en
algunos puntos de observación, llamados pozos de observación. Luego se realiza
el ensayo de recuperación y consiste en medir la recuperación de los niveles del
pozo una vez que se ha parado el bombeo. Es importante que en todo ensayo
hidráulico se midan las evoluciones de caudales bombeados, que generalmente se
los realizan a caudal constante.
2.4 Explotación y sobre explotación de acuíferos.
La explotación de acuíferos, data desde la época en que el hombre preciso extraer

agua del subsuelo, en sitios en donde las aguas superficiales eran escasas o
nulas. Poco a poco se fueron desarrollando técnicas de extracción de agua
21
subterránea, inicialmente para el consumo doméstico, luego cuando se desarrolla

la agricultura, se usa para el riego de inmensas zonas con cultivos. Es importante
que dicha explotación de acuíferos, tenga en cuenta extraer la cantidad necesaria y
tratar de que en este proceso se observe asegurar la buena calidad de la misma.
Al ser el agua subterránea un recurso vulnerable, que al afectarla se crean

consecuencias no deseables en el acuífero que pueden subsistir largos periodos
de tiempo. El concepto de calidad está ligado profundamente a dos aspectos
fundamentales, como son: los condicionantes naturales que alteran o controlan la
composición química de las aguas subterráneas ; así como el uso que se pretende
dar a dicha agua con las políticas de gestión y planificación asociadas a dicho
recurso. Conjuntamente, desde el punto de vista de la gestión del recurso hídrico,
la valoración de la calidad del agua está muy ligada al uso que se le quiera dar.
Porque, no es igual el enfoque de un estudio que trate de valorar el agua para
suministro urbano como el que se realizaría para un uso industrial (Fraile, J.;
Escuder, R.; Jordana, S.; Ribera, F.; Sanchez-Vila, X.; Vásquez-Suñé, E., 2009)
Una de las principales ventajas del agua subterránea respecto a la superficial

circundante reside en que, generalmente, su calidad química y bacteriológica suele
ser superior. (Benitez, A., 1972). Esto se debe principalmente a la capacidad de
retardo de la velocidad de propagación, de la disminución de la concentración o de
la disminución de gran parte de los compuestos disueltos o transportados por el
agua. En concreto, la cuantificación del grado de autodepuración del agua al
incorporarse al medio subterráneo depende de la interrelación de diversos
procesos, entre los que se destacan los derivados de la propia composición
química del agua y de sus características físicas, de las características geológicas
de la zona donde se produce la recarga y de la historia de usos del suelo en ese
sector. El aspecto más interesante del mecanismo de autodepuración del terreno
es que éste se produce de forma natural, durante el flujo del efluente a través de la
zona no saturada (ZNS) y, en menor grado, en la zona saturada (ZS) del acuífero.
Una cuantificación más precisa de la capacidad de autodepuración requiere
estudios específicos que tendrían que ser considerados en las etapas más
avanzadas de cada proyecto, en concreto. (Fraile, J.; Escuder, R.; Jordana, S.;
Ribera, F.; Sanchez-Vila, X.; Vásquez-Suñé, E., 2009). Entonces se debe tener en
cuenta al hacer la explotación de los acuíferos, el caudal del que se puede
disponer y sus reservas, para no sobre explotarlos.
22
3 Modelos Matemáticos de Simulación de Agua Subterránea
Los modelos matemáticos son fundamentales para la resolución de problemas de

agua subterránea, dado que no todos los problemas de aguas subterráneas y flujo
en medios porosos, pueden ser resueltos a través de soluciones cerradas.
3.1 Tipos de Modelos
Según la implementación de los modelos, los mismos pueden clasificarse en:

FÍSICO: representan a escala la geometría y propiedades de los materiales de un
sistema: maquetas en arquitectura y modelos hidráulicos de un canal o una central
hidroeléctrica a escala.
ANALÓGICO: representan las propiedades del sistema a través de fenómenos
diferentes pero que tienen un comportamiento similar.
CONCEPTUAL: identifican las características más relevantes de un sistema, sin
llegar a representarlas cuantitativamente. Se utilizan para tener una mejor
comprensión del sistema real.
MATEMÁTICO: plantean relaciones principales entre elementos relevantes de un
sistema, se reemplazan por relaciones matemáticas que las representan en forma
aproximada.
Resumiendo, las potenciales aplicaciones de modelos de agua subterránea son:
PREDICTIVAS: Permiten predecir cambios futuros. Deben ser calibrados para
representar zona de interés.
INTERPRETATIVAS: Se utilizan para interpretar la dinámica de un sistema y para
mejorar conocimiento del mismo.
GENERICAS. Colaboran en el análisis de condiciones hipotéticas para comprender
fenómenos. (Espinoza C., 2012)
Una numeración de ejemplos concretos de aplicaciones se presenta a

continuación:
 Manejo o gestión de recursos hídricos.

Determinar tasas de bombeo máximas en un sector.
Calcular áreas de protección del agua subterránea.
 Predecir impactos sobre el agua subterránea a través de casos estudio.
Influencia de una presa sobre el flujo regional de agua subterránea.
Analizar drenaje de exacavacion de una mina.
 Diseñar y optimizar sistemas de remediación.
 Analizar riesgos ambientales.
Determinar dirección más probable de una pluma de contaminante.
 Mejorar conocimiento de la hidrogeología.
(Reyna, Reyna, & Labaque, 2010)
23
A continuación se describe el modelo de diferencias finitas MODFLOW.
3.2 MODFLOW 3-D
Preliminares.
Es importante describir el programa computacional que se utiliza en la presente

investigación, de esta manera el lector tiene una noción de las propiedades y
limitaciones del programa y como este efectúa las aproximaciones, en base a los
datos de entrada registrados, así como también interpreta los resultados.
MODFLOW muestra una estructura mediante sistemas patrones, compuestos de

un esquema principal y una sucesión de subrutinas en paquetes. Cada paquete de
MODFLOW representa alguna característica específica del sistema hidrológico,
como recarga, evapotranspiración, pozos, ríos y drenes. Además, desde su
desarrollo original han sido incorporados al mismo, un amplio rango de módulos
adicionales, tales como anisotropía variable en el espacio y re-humedecimiento de
celdas.
La investigación pormenorizada de este modelo y sus referidos paquetes se

presenta en (Mc Donald, M.G., y Harbaugh, A.W., 1988), (Harbaugh, A.W., 2005), y
en escritos técnicos de libre acceso desarrollados por el USGS para MODFLOW, o
Schlumberger en su adaptación comercial. En el presente subtítulo se desarrollan
en términos generales, la estructura del programa principal y los paquetes de
cálculo que son empleados en las modelaciones, además se resumen las
fortalezas y debilidades que deben considerarse en su utilización.
Generalizaciones
El programa que se hace referencia, es el modelo hidrogeológico con mayor

utilización en el mundo científico. Ofrece una simulación versátil, ya sea con flujo
estacionario o transitorio, en escenarios saturados confinados, no confinados o
ambos.
3.2.1 Características Importantes
• Representa la salida de agua subterránea en (Plano X, Plano Y,

Plano Z) en un medio heterogéneo y anisotrópico.
• Admite representar flujo en estado permanente y transitorio en
acuíferos del tipo no confinados, confinados y modelos confinantes
• Admite simular los efectos de pozos, recargas, ríos, drenes y
evapotranspiración.
• Muestra una distribución por módulos en el programa y en la
documentación que facilita una mejor comprensión del modelo.
24
• Usa diferencias finitas (calcula de manera aproximada las soluciones

usando ecuaciones diferenciales finitas para aproximar derivadas).
• Permite la utilización de condiciones de frontera de Cauchy,
condiciones de Dirichlet y de Neuman.
• El módulo MODPATH, acepta simular el rastreo de partículas y el
módulo MT3D, permite simular el transporte de contaminantes.
3.2.2 Solución de la ecuación de flujo por diferencias finitas.
Un procedimiento de aplicación de patrones matemáticos en hidrogeología, es el

método de diferencias finitas, que no es otra cosa que la determinación de valores
finales aproximados, partiendo de valores de inicio comprobados, que aplican a
una función por intermedio del cálculo en fases finitas. (Fulginiti F., 2013)
La discretización espacial en MODFLOW se obtiene tomando en cuenta acuíferos
reales, para esto se elabora una red conformada por celdas, cada celda está
caracterizada en términos de filas (I), columnas (J), y capas (K) (Ilustración 3-1). La
mencionada celda se considera uniforme con respecto a sus propiedades. Existe
un punto por celda denominado nodo, en el cual es calculada la altura de presión.
COLUMNAS (J)
1 2 3 4 5 6 7 8 9
1
2 Límite del Acuífero
Celda Inactiva
FILAS (I) 3 Celda Activa
vk
4 rj Dimensión de la celda a lo largo
5 de la dirección de las filas
ci Dimensión de la celda a lo largo
1 de la dirección e las columnas
2 vk Dimensión de la celda a lo largo
CAPAS (K) 3 de la dirección vertical
ci
4
5
rj
Ilustración 3-1 : Sistema acuífero hipotético discretizado (Fuente: (Harbaugh, A.W., 2005)).
Según la ubicación de los nodos, se obtienen dos clasificaciones: Celda Central(los

nodos están centrados equidistantes de los contornos de las celdas) y punto
central(los nodos se encuentran en donde se intersectan las celdas).
El flujo de agua subterránea es simulado por MODFLOW usando solo la
configuración de bloque central de diferencias finitas (Mc Donald, M.G., y
Harbaugh, A.W., 1988).
En la Ilustración 3-2 y la Ilustración 3-3 se pueden apreciar ambas configuraciones
de resolución en un sistema bidimensional.
25
rj
ci
Ilustración 3-2 Formulación de Bloque centrado (Fuente: (Harbaugh, A.W., 2005))
rj
ci
Ilustración 3-3 Formulación de Punto central (Fuente: (Harbaugh, A.W., 2005))
Tanto en la formulación de bloque central como en la de punto central, el espacio

entre nodos tiene que ser seleccionado de tal forma que las propiedades
hidráulicas sean uniformes en el bloque.
La elección de como transcurre el tiempo y la elaboración de la red son secuencias
críticas en el esquema del modelo, en razón de que los valores de la discretización
espacial y temporal, poseen un predominio importante en los resultados numéricos
(Anderson, M. P. y Woessner, W. W.,, 1992).
En este programa computacional, se pueden hacer corridas (run) de modelos en
estado estacionario, como también en estado transitorio. Cuando se eligen
simulaciones en estado estacionario, son para saber cómo transitan las
aproximaciones de las condiciones de equilibrio en el largo plazo. Cuando se eligen
simulaciones en estado transitorio, se efectúan a partir de un tiempo inicial e
inclusive un tiempo final, para esto se tiene como patrón una serie de sucesiones
en el tiempo. Los niveles del medio saturado, cambian en función de las
variaciones que presenten las condiciones de frontera, las afectaciones externas, la
transmisividad y las propiedades de almacenamiento del medio.
26
3.2.3 Datos de Salida MODFLOW
El programa arroja la distribución del nivel piezomètrico, para el intervalo de tiempo

analizado. El programador tiene la posibilidad de controlar la frecuencia, con la que
los datos de salida se presentan. Existen también otras elecciones de datos de
salida que contienen descensos de nivel y el balance de masas entre regiones que
han sido definidas con anterioridad por el usuario.
3.2.4 Fortalezas y Debilidades del Programa
El mencionado programa, como la mayoría de programas computacionales

presenta propiedades que son fortalezas y limitaciones que son debilidades que se
describe brevemente.
Fortalezas del Programa:
Se ha comprobado la eficiencia en la resolución de problemas de flujo en las

dimensiones X, Y,Z.
Entre las principales ventajas del mismo se enfatizan:
 El esquema numérico de resolución por diferencias finitas, entrega

soluciones aceptables a los problemas modelados, a través de una
metodología numéricamente sencilla.
 Permite su aplicación a una gran variedad de condiciones del mundo real.
 La representación gráfica de la modelación, facilita el ingreso de los datos y
la visualización de resultados.
 MODFLOW puede aplicarse como un modelo unidimensional, bidimensional,
o cuasi o totalmente tridimensional
 Al ser un modelo cuasitridimensional, Visual MODFLOW está capacitado
para representar de una manera más realista los problemas que se
presentan en el mundo real.
 El diseño del programa MODFLOW, (en su versión de código abierto)
permite que nuevas características de la simulación sean agregadas con
relativa facilidad.
 MODFLOW ha sido aceptado en muchos casos de la Corte en los Estados
Unidos, como una herramienta legítima para el análisis de sistemas de agua
subterránea.
 Una amplia variedad de programas, están disponibles para leer los datos de
salida desde MODFLOW y presentar gráficamente los resultados del modelo
de manera que se entiendan fácilmente.
Debilidades del Programa:
Entre las debilidades que presenta el programa se pueden enunciar las siguientes:
27
 El esquema numérico de resolución, (como en la mayoría de los programas

similares) presenta limitaciones en el espaciado y la densidad de la malla, lo
que en modelos de gran extensión es un factor de importancia.
 La resolución numérica al utilizar el método de diferencias finitas, trabaja
mejor con acuíferos rectangulares o prismáticos y de composición uniforme,
por lo que los modelos de acuíferos irregulares o curvos, con propiedades
de anisotropía y heterogeneidad son difíciles de resolver mediante esta
metodología.
 Utilizar diferencias finitas implica una menor precisión de las soluciones
entregadas, en relación a la utilización de técnicas como los elementos
finitos.
 Una debilidad de Visual MODFLOW, producto de su formulación numérica,
radica en la imposibilidad de rehumedecer celdas que se vacían (o secan)
durante el proceso de cálculo. Esto lleva a cometer errores de cálculo al
momento de cuantificar caudales de extracción o direcciones de flujos. El
paquete de rehumedecimiento (REWET) puede salvar esta dificultad si es
aplicado de forma adecuada.
 El agua debe tener una densidad constante, viscosidad dinámica y en
consecuencia temperatura igual durante todo el modelo (SEAWAT-2000 es
una versión modificada de MODFLOW-2000) la cual está diseñada para flujo
y transporte de agua subterránea dependiente de la densidad.
 Las opciones para el cálculo de la conductancia horizontal bajo condiciones
de acuífero freático, la corrección del flujo vertical, y la conversión del
término de almacenamiento, se desarrollaron para el supuesto de que cada
capa del modelo corresponde a distintos acuíferos o distintos horizontes
permeables, y que estos horizontes se encuentran separados por unidades
distintas de baja permeabilidad. El uso de estas opciones donde estas
condiciones no se cumplen, puede conducir a una variedad de problemas e
imprecisiones en la simulación.
 Por ejemplo, si la opción para el cálculo de la conductancia horizontal bajo
condiciones de la tabla de agua, se utiliza cuando un acuífero freático está
representado por varias capas del modelo, y la tabla de agua se espera que
atraviese más de una capa durante la simulación, la conversión incorrecta
de las células a un no flujo de condición puede presentarse como iteraciones
se llevan a cabo. Si la opción de humectación no se utiliza, la conversión a
ningún flujo es irreversible. La opción de humectación permite la conversión
a ningún flujo a ser invertida, sin embargo, la conversión hacia atrás y
adelante a veces pueden producirse varias veces de modo que la
convergencia solucionador no se logra. Por lo tanto, se debe tener cuidado
en la decisión de utilizar estas opciones. (Fulginiti F., 2013)
28
Como sucede con todo modelo numérico el objetivo de la modelación de un

sistema hídrico subterráneo con Visual MODFLOW es la obtención de los
siguientes resultados:
– evolución de los niveles de agua subterránea, h(x, y, z, t).
– velocidades de flujo en el entorno definido como zona de interés, v(x, y, z, t).
ENTRADA DE EJECUTAR EL SALIDAS-

DATOS MODELO RESULTADOS
Ilustración 3-4 .- Visual Modflow Menú Principal
3.3 Estructura de Modelos Hidrogeológicos
Los modelos hidrogeológicos pretenden representar de manera física, matemática,

en su forma más simple, lo complicado que son los procesos hidrológicos,
hidrogeológicos e hidro geoquímicos que se despliegan en un sitio determinado.
Usando las ecuaciones que gobiernan el flujo y la conservación de la masa en un

medio poroso, expresadas en diferenciales, representan los variados escenarios de
estimación para el flujo de agua y transporte de sustancias contaminantes en el
subsuelo.
El Centro Internacional de Modelismo de Agua Subterránea (IGWMC, 1993)define

un modelo de simulación hidrogeológico como “una descripción matemática, no
única, simplificada de un sistema de aguas subterráneas existente, codificado en
un lenguaje de programación, junto con una cuantificación del sistema acuífero que
incluya las condiciones de borde relevantes, los parámetros del sistema, así como
las presiones o acciones sobre éste”. La solución de un modelo de simulación
hidrogeológico, corresponde a la distribución de los niveles de agua (acuífero libre)
29
o niveles de energía o piezométricos (acuífero confinado) a través del espacio y

tiempo. (Fulginiti F., 2013)
3.4 Etapas en el desarrollo de un modelo de Simulación

Hidrogeológico
Estas pueden a la vez agruparse en tres macroactividades o etapas. La primera

etapa se orienta a elaborar un modelo conceptual de la zona de estudio (A). La
segunda etapa es la construcción del modelo propiamente dicho (B), y la tercera
etapa incluye su uso posterior (C). El modelo se construye en base a
observaciones reales y percepciones sobre el sistema hidrogeológico y el problema
a analizar (conceptualización).¡Error! No se encuentra el origen de la referencia.
30
IDENTIFICAR EL PROBLEMA PASO 1
DEFINIR ÁREAS DE INTERÉS PASO 2 A
ELABORAR MODELO CONCEPTUAL PASO 3
DATOS
SELECCIONAR MODELO MATEMÁTICO PASO 4
IDEALIZAR SISTEMA FÍSICO PASO 5
PROPIEDADES HIDROGEOLÓGICAS Y CONDICIONES

DE BORDE PASO 6
B
ACCIONES EXTERNAS PASO 7
CALIBRACIÓN / VALIDACIÓN / ANÁLISIS DE

SENSIBILIDAD PASO 8
VALIDACIÓN
OPERACIÓN DEL MODELO PASO 9
C
POSTAUDITORIA PASO 10
Ilustración 3-5 Etapas en el desarrollo de un modelo Hidrogeológico (Fuente: (Facultad de

Ciencias Fisicas y Matematicas de la Universidad de Chile, 2010)
3.5 Modelos de Simulación
En las simulaciones del agua subterránea, se utilizan diferentes tipos de modelos

según se pretende conocer el comportamiento en la zona saturada o no saturada,
según la extensión de los límites de simulación, los procesos intervinientes o los
parámetros involucrados.
En términos generales, los que más se han usado son los modelos físicos,
eléctricos, analógicos y matemáticos, siendo estos últimos los que mayor
relevancia presentan en la actualidad.
31
Existen diversos modelos que permiten predecir el movimiento del agua y los
químicos a través de la zona saturada o no-saturada del suelo. Estos modelos
numéricos necesitan de la estimación de numerosos parámetros.
3.5.1 Modelos de medio poroso.
El modelo de medio poroso o modelo a escala pequeña, es uno de los modelos

ampliamente utilizados en ingeniería hidráulica. El sistema de agua subterránea se
representa en una escala apropiada en el laboratorio con sus fronteras incluidas.
En los modelos de medio poroso, las propiedades hidrogeológicas se distribuyen
en espacio y magnitud de acuerdo al sistema natural. Cuando se estudia el flujo
libre del agua subterránea, son necesarias correcciones para el ascenso por
capilaridad, ya que este es desproporcionalmente grande en el modelo, comparado
con las condiciones de campo. Las líneas de corriente se pueden producir por
medio de inyección de tinta. Los modelos de medio poroso, como las cajas de
arena, aún se usan para investigación. En el estudio de problemas de campo, los
modelos de medio poroso han sido reemplazados casi totalmente por los modelos
numéricos.
3.6 Modelos matemáticos para la simulación del flujo subterráneo en

un medio saturado
Los modelos matemáticos se han convertido en una herramienta esencial para el
estudio de los problemas complejos del flujo subterráneo y el transporte de solutos
en el subsuelo, ya que proporcionan un marco cuantitativo que sintetizan un gran
número de parámetros que describen la variabilidad espacial dentro del subsuelo.
Los modelos matemáticos son empleados en distintas aplicaciones o códigos para

computadoras, con el fin de resolver los sistemas de ecuaciones que definen los
procesos físicos y químicos dentro del medio poroso.
Los códigos más usados son aquellos que simulan el flujo, con o sin transporte de
solutos, en un medio poroso continuo (en contraste con el medio fracturado).
Existen códigos analíticos que permiten evaluar el comportamiento de un acuífero

de características uniformes ante situaciones locales de bombeo, como
AQTESOLV ( (AQTESOLV), o para la evaluación de ensayos de bombeo, como
TECTYPE o PUMPTEST. Sin embargo la mayoría de los análisis de flujo
subterráneo son estudiados mediante códigos numéricos debido a la alta
complejidad que presentan.
Entre los principales códigos numéricos existentes se destacan:
• MODFLOW, (United State Geological Survey, USGS)
• FLOWPATH (Waterloo Hydrogeologic, Inc.)
• FRAC3DVS (Waterloo Hydrogeologic, Inc.)
32
• FEFLOW (DHI-WASY)
• 3DFEMWATER (US.EPA)
• ASMWIN
• PLASM (Organismo de Control de Aguas del Estado de Illinois, EEUU).
• AQUIFEM-1 (Instituto Tecnológico de Massachussets).
Estos códigos requieren generalmente como entrada la distribución espacial de las

propiedades del acuífero, en general, conductividad hidráulica, transmisividad y
parámetros de almacenamiento, además del potencial inicial y las diferentes
aportaciones o extracciones que se producen en el sistema.
Por lo general los mismos paquetes de software incluyen códigos para simular el
transporte de contaminantes. Estos simulan los procesos de advección, dispersión
y las reacciones químicas entre los diferentes solutos disueltos en el subsuelo.
Entre los principales códigos se encuentran RNDWALK (Centro de Investigación y
Estudios en Medio Ambiente, 2003), como módulo de PLASM y MT3D como
módulo de MODFLOW. Generalmente estos modelos requieren como entradas la
velocidad, el coeficiente de dispersión longitudinal, la concentración inicial, la
porosidad efectiva, el factor de retardo y la vida media de los contaminantes, entre
otros.
Teniendo en cuenta tanto los códigos que simulan el flujo subterráneo, como los
que contemplan el transporte de solutos, se pueden citar entre los principales los
siguientes:
3.7 Flujo en sistemas de acuíferos
3.7.1 FLOWPATH
FLOWPATH (Waterloo Hydrogeologic, Inc., 1996) es un modelo de líneas de

corriente y flujo de agua subterránea bidimensional que trabaja bajo estado
estacionario. El código puede simular acuíferos resquebrajados, confinados y no
confinados en medios heterogéneos y anisotrópicos. Permite simular condiciones
de contorno complejas. La información de salida del programa incluye niveles
piezométricos, líneas de corriente, tiempos de viaje, velocidades y balances
másicos. Las paredes de las pantallas, pueden ser simuladas construyendo una
rejilla modelo con tamaños de celda variables. No puede ser utilizado si el flujo de
agua subterránea del lugar a simular es muy complejo debido a flujos verticales o si
han de ser simulados campos de flujo transitorio.
33
3.7.2 MODFLOW
Es uno de los códigos de modelización de agua subterránea más versátil, de

amplia utilización y difusión. Desarrollado por el United State Geological Survey, se
trata de un modelo en 3-D que trabaja bajo el método de diferencias finitas
(MacDonald & Harbaugh., 1988,1996). MODFLOW simula el flujo transitorio de
agua subterránea en dos dimensiones y casi o completamente en 3 dimensiones
en sistemas de acuíferos anisotrópicos, heterogéneos y estructurados por capas.
Está escrito en Fortran 77 (Mc Donald, M.G., y Harbaugh, A.W., 1988). Permite
calcular distribuciones de niveles piezométricos, de velocidades de flujo y balances
de agua. Incluye módulos para flujos hacía pozos, a través de cauces y en
desagües. Otros módulos se ocupan de la evapotranspiración y recarga.
Disponibles en el mercado se encuentran pre- y post-procesadores que facilitan el
uso del código así como el análisis de los resultados. Entre estos se encuentra
Visual MODFLOW desarrollado por Waterloo Hydrogeologic, Inc. (1996).
3.7.3 VISUAL MODFLOW
Existen otras versiones que, utilizando el motor de cálculo original, le añaden una
interfaz cómoda para el usuario, tanto para introducir las características del modelo
como para observar al final los resultados. La más extendida es Visual Modflow,
incluye además a otros programas que fueron creados como complementos de
MODFLOW; utilizan la salida de MODFLOW y proporcionan los siguientes
resultados: MODPATH describe la trayectoria descrita por las partículas de fluido
que se haya indicado; ZONE BUDGET calcula el balance hídrico de cada zona que
se establece dentro del modelo, MT3D calcula el transporte de solutos. (J. Sánchez
San Román, 2008)
Visual Modflow es una interfaz gráfica de Modflow, dicho software fue desarrollado
por el U. S. Geological Survey (MacDonald & Harbaugh., 1988,1996). El programa
fue implementado en 1980 y ha sido continuamente modificado con el desarrollo de
muchos paquetes y programas relacionados con el estudio de las aguas
subterráneas. (Herrera- Carbajal, Villalba-Atondo, & Molinar-Tabares, 2006)
Este programa es frecuentemente utilizado para la simulación del comportamiento
de las aguas subterráneas sobre la base del método de las diferencias finitas en 1,
2 ó 3D.
En este método, se divide el acuífero por medio de una malla en bloques
rectangulares, los cuales se organizan en filas, columnas, y capas, denominándose
a cada bloque “celda”. (Romero & Mercado, 2007)
El modelo, Visual MODFLOW, es una herramienta computacional muy útil y fácil de
operar, por sus características visuales y de administración de datos. Como
descripción general de este modelo se puede mencionar que es capaz de abordar
simulaciones del flujo de aguas subterráneas en tres dimensiones, las cuales en
conjunto con una serie de paquetes que posee el programa (M3TD y RT3D)
34
permiten abordar problemas de transporte de contaminantes. Presenta ventajas en

la representación gráfica de los sistemas a modelar, facilitando la entrada de datos
para definir el escenario de simulación. Este modelo tridimensional de diferencias
finitas, permite hacer una estimación de la geometría del área y de sus límites de
estudio.
Para definir dentro de Visual MODFLOW la geometría del sistema, se debe indicar
las dimensiones en planta de la zona y luego el número de columnas, filas y las
capas. Esta generalmente es utilizada principalmente para desarrollar modelos
conceptuales y teóricos, modelos de carácter local como el diseño de pruebas de
bombeo y sistemas de drenaje.
3.7.4 FRAC3DVS
FRAC3DVS es un modelo en 3D, basado en el método de los elementos finitos

para simular flujos de agua subterránea en estado estacionario o transitorio,
saturado o variablemente saturado, así como el transporte de solutos, advectivo-
dispersivo, en medios porosos o discretamente fracturados. El código fue
desarrollado en la Universidad de Waterloo (Therrien R., 1992) y está siendo por
Waterloo Hydrogeologic, Inc. El código incluye pre-procesadores para mallas de
rejillas y generación de archivos de entrada y post-procesadores para la
visualización de los resultados de la simulación. Este programa presenta un gran
número de características avanzadas en relación al resto de los programas, a la
vez que se incrementan la cantidad de parámetros necesarios.
FRAC3DVS pueden utilizar un bloque de 8 nodos o 6-nodos de elementos finitos o

prisma de 7 puntos de diferencia finita para la discretización del medio poroso y
elementos planos rectangulares o triangulares que representen fracturas. Los
elementos de bloque se pueden subdividir en tetraedros de igual modo que las
redes compuestas de elementos de bloque no ortogonales.
Incluye los Procesos de transporte, advección, dispersión hidrodinámica, sorción

según una isoterma de lineal de Freundlich, y el transporte de múltiples especies
de cualquiera de las cadenas de desintegración, rectas o ramificadas.
3.7.5 FEFLOW
FEFLOW fue desarrollado por DHI-WASY filial alemana de DHI Group.
Los componentes primarios incluyen:
- Una completa selección de herramientas gráficas que permiten construir la

malla de elementos finitos, asignar propiedades y definir condiciones de
borde.
35
- Se integra totalmente con ArcView y ArcInfo (ESRI) GIS, permitiendo

compartir los datos de modelación de aguas subterráneas, con una base de
datos GIS. Es compatible con WGEO, una excelente herramienta para la
georreferenciación fácil y rápida, procesamiento de datos de trama
especialmente relacionados.
- Integra técnicas de solución y algoritmos muy robustos.
- Permite abordar problemas de variada indole como: Flujos de densidad

variable (salt water intrusion).Flujo permanente e impermanente. Flujo
saturado y no saturado. Múltiples superficies libres (acuíferos colgados).
Transporte de masas y calor
- Permite interpretar datos en tiempo real y visualizarlos en 3D.
(https://www.mikepoweredbydhi.com/products/feflow)
3.7.6 FEMWATER
FEMWATER (desarrollado por El Cuerpo de Ingenieros de los Estados Unidos) es

un modelo de elementos finitos 3D, totalmente utilizado para simular el flujo
acoplado densidad impulsada y transporte de contaminantes en las zonas
saturadas e insaturadas. FEMWATER permite el modelado de la intrusión de la
salinidad y otros contaminantes dependientes de la densidad. Estratigrafía
compleja puede desarrollarse en GMS y directamente representada en el modelo.
Las soluciones se pueden visualizar mediante gráficos realistas 3D y secuencias
de animación.
GMS proporciona una interfaz personalizada para el modelo FEMWATER, que

ofrece una forma sencilla de configurar los parámetros del modelo y una interfaz
gráfica de usuario para ejecutar el modelo y visualizar los resultados. Recopilar
datos sobre los antecedentes de una variedad de fuentes de SIG a CAD y acceder
a los datos en línea a partir de numerosas bases de datos de mapas, imágenes y
datos de elevación. GMS le permite interactuar con los modelos en 3D real que se
aprovechan de gráficos OpenGL y optimizados para crear representaciones
fotorrealistas y animaciones para presentaciones de PowerPoint, impresión y web
3.7.7 ASMWIN
ASMWIN (modelo de simulación de acuíferos para Windows) es una versión

mejorada de ASM 5.0 que corría bajo MS-DOS (Kinzelbach W. Y Rausch R.,
1995), implementa un modelo completo de flujo subterráneo bidimensional y
modelo de transporte. ASMWIN posee una interfaz gráfica para el usuario, un
modelo de flujo por diferencias finitas, una herramienta para la calibración
automática del modelo de flujo, un modelo de rastreo de partículas, un modelo de
36
transporte aleatorio, un modelo de transporte por diferencias finitas y muchas otras

herramientas de modelaje.
Las ecuaciones de flujo discretizado se resuelven por el método de los gradientes

conjugados precondicionados (PCG) (módulo ASMSIM1) con elección de la
precondición diagonal y de Cholesky. Para campos de flujos estacionarios, un
procedimiento automático de calibración de modelos utilizando el algoritmo de
Marquardt-Levenberg está disponible en el módulo ASMOPTI.
El módulo de rastreo de partículas ASMPATH ofrece muchos métodos de

interpolación de velocidades y utiliza la integración de Euler para calcular las
trayectorias y tiempos de recorrido del flujo estacionario o transitorio. Se disponen
dos módulos de simulación de transporte. El primero usa un esquema de
diferencias finitas (Módulo ASMT2SIM) y el segundo un método aleatorio basado
en la teoría de Ito-Fokker-Planck (Módulo ASMWALK).
3.8 Transporte de contaminantes
3.8.1 MODPATH
Los resultados de MODFLOW pueden ser aplicados en códigos de rastreo de

partículas, como MODPATH (Pollock D. W., 1989), para calcular sendas de agua
subterránea y tiempos de viaje. MODPATH es un código de postprocesado,
utilizado para computar líneas de sendas de agua subterránea tridimensionales,
basadas en la información de salida de simulaciones en estado estacionario
obtenidas por el código de modelización MODFLOW. MODPATH usa un esquema
semi-analítico de rastreo de partículas, basándose en la hipótesis de que cada
componente de velocidad direccional varía linealmente dentro de una célula de
rejilla en su propia dirección de coordenada.
3.8.2 PATH3D
PATH3D (Zheng C., 1989) emplea también los resultados de MODFLOW para
calcular sendas de agua subterránea y tiempos de viaje. PATH3D es un programa
de rastreo de partículas general, para calcular sendas de agua subterránea y
tiempos de viaje en campos de flujo transitorios y tridimensionales. El programa
incluye 2 principales segmentos, un interpolador de velocidad, el cual convierte los
niveles piezométricos generados por el MODFLOW en un campo de velocidad, y
un programa de resolución numérica de cuarto orden, basado en el método
numérico de Runge-Kutta con ajuste del tamaño de intervalo de tiempo automático
para el rastreo del movimiento de las partículas de fluido. (Van der Heijde,
Elnawawy, 1993)
37
3.8.3 MT3D
Es un módulo que forma parte de Modflow. Este modelo 3D se emplea en la

simulación del transporte de contaminantes. Permite simular procesos de
advección, dispersión, sumideros, fuentes, procesos de mezcla y reacciones
químicas de disolución de constituyentes en sistemas acuíferos.
3.9 Validación y Análisis de Sensibilidad.

Al realizar estos elementos, se debe tener en cuenta algunas consideraciones para
lo cual se debe analizar los siguientes aspectos:
¿Son los supuestos del modelo válidos?

¿Son razonables los datos básicos utilizados?
¿Cuál es el efecto de pequeñas variaciones en los datos de entrada?
¿Se puede defender el modelo? (Espinoza, 2012)
3.10 Hipótesis de simulación del modelo
Ilustración 3-6 Registro de monitoreo pozo
De acuerdo a Ilustración 3-6, que muestra un registro de monitoreo de pozo, en

donde están las variaciones de profundidad de dicho pozo a lo largo de un tiempo t
determinado, la curva de datos medidos y los calculados, lo ideal es que los datos
calculados se acercan con mayor precisión a los datos medidos, es decir, la
calibración del modelo, está simulando la realidad, que es lo que se quiere, para a
partir de eso, simular los escenarios futuros y que estos sean confiables.
38
Las hipótesis básicas del modelo de flujo en medios porosos son:

La matriz sólida está saturada de fluido.
Puede haber o no difusión del fluido.
La matriz sólida está en reposo.
La matriz sólida es elástica.
El fluido es comprensible.
La velocidad del fluido está dada por la ley de Darcy. (Suarez E. L., 2010)
39
4 Caracterización Regional de Cuenca del Milagro y Suquía

A continuación se presenta la caracterización de la Cuenca del Milagro y del Río
Suquía, dentro de estas cuencas se encuentran los estudios de flujo en acuíferos
que se midieron y analizaron luego con el uso de modelos matemáticos
4.1 Cuenca del Río Suquía. Provincia de Córdoba-Argentina
En una clasificación general y esquematizada de la totalidad de la provincia de

Córdoba, la llanura del Este pertenece a lo que en el país se considera la gran
Cuenca Chaco-Paranense, y las del Noroeste y Oeste, forman parte de los Llanos
occidentales. La cuenca Chaco-Paranense se divide en la Provincia en dos
unidades bien definidas: la de Mar Chiquita y la de los ríos Tercero, Cuarto y
Quinto. De acuerdo a estas consideraciones, se definen las siguientes cuencas
hidrogeológicas (Ilustración 4-1) (Vasquez, Miatello, Roque, & otros, 1979) :
I) Cuenca de Mar Chiquita (subregiones Río Dulce, Sistema Río Suquía y Río
Xanaes)
II) Cuenca de los ríos Tercero (Ctalamochita), Cuarto (Chocancharava) y
Carcarañá
III) Cuenca de las Salinas Grandes (subregiones Río Cruz del Eje, Río Soto, Río
Pichanas, Sistema Noroeste y Guasapampa)
IV) Cuenca del Conlara
V) Cuencas Intermontanas
VI) Cuenca de la Llanura Medanosa, subregión Río Quinto (Popopis)
40
Figura I.1: Cuencas subterráneas: I, Cuenca de Mar Chiquita; II, Cuenca de los Ríos Tercero
Ilustración 4-1.-Cuencas
III, Cuencayde las Salinas Grandes; IV, Cuenca –
Hidrogeológicas
Cuarto, del Pcia.
Conlara; V,de Cba. Fuente (Vasquez, Miatello,
Cuenca
Intermontanas; VI Cuenca del Ríos Quinto, (o)
Roque, & otros, 1979)
Perforaciones Perforaciones. (Vázquez y otros, 1979).
En la República Argentina, la importancia de las aguas subterráneas es

fundamental, porque las fuentes de agua superficiales se encuentran mal
distribuidas y en algunas regiones son escasas.
La cuenca hidrológica de Mar Chiquita limita al Oeste con los cordones de las
Sierras Pampeanas; al sur con la dorsal, que con rumbo Este-Oeste separa la
llanura del Río Segundo o Xanaes con la del Río Tercero o Ctalamochita siguiendo
una línea que pasaría por las localidades de San Agustín, Oliva, El Fortín, es decir,
aproximadamente siguiendo el paralelo de los 32°. Al Este, limita por el bloque
elevado, en la línea San Francisco – Selva - Tostado, en el límite con la provincia
de Santa Fe. El límite Norte, impreciso aún, está dentro de la provincia de Santiago
del Estero (C.A.A.A.S., 1975)
La cuenca de interés es la del Río Primero o Suquía, en esta cuenca el promedio
de lluvias anuales varía de 600 a 800 mm por año. Los meses de mayo a
septiembre, conforman el período seco y los de octubre a marzo el de lluvias,
siendo mayores las precipitaciones en los meses de diciembre, enero y febrero.
41
Tabla 4-1 Parámetros climáticos promedio de Córdoba. (SMN, 2013)

(www.freemeteo.com)
Datos Estación Meteorológica Córdoba

900 35
800 30
700
Precipitación(mm)
Temperatura(°C)
25
600
500 20
Precipitación
400 15
300
10 Temperatura Diaria
200 Minima
100 5
Temperatura Diaria
0 0
Máxima
Julio
Mayo
Diciembre
Agosto
Anual
Enero
Abril
Febrero
Octubre
Marzo
Junio
Septiembre
Noviembre
Meses
Tabla 4-2 Precipitación –Temperatura Anual Años 1961-1990 (SMN, 2013)
42
4.2 Límites de la Zona de estudio
4.2.1 Cuenca del Río Suquía-Argentina
Ilustración 4-2 Cuenca Hidrográfica Río (Suquía).
El Río Suquía o Primero nace de las confluencias de los ríos San Antonio y
Cosquín a los que se les une el de Los Chorrillos y el Arroyo de Las Mojarras. En la
actualidad sus aguas se encuentran represadas por el dique San Roque. Su
cuenca activa está limitada al Norte por el dorso de La Cumbre; al Sur, por la
cresta de la Sierrita o Cordón de Santiago, que la separa de la cuenca del Río
Anizacate; al Oeste, por las divisorias de las aguas del Río Pintos, de la Pampa de
San Luis y de la Sierra Grande. La región que tiene como base de desarrollo a las
aguas del Río Suquía, es relativamente pequeña, pero dentro de ella se encuentra
ubicada la segunda concentración demográfica de la República (Ciudad de
Córdoba), con una actividad industrial de gran solidez e importancia, un importante
centro turístico, y con una no despreciable potencialidad agrícola-ganadera y
minera. Antes de su paso por la ciudad de Córdoba, se encuentra regulado aguas
arriba con el embalse San Roque, además sufre varias obstrucciones a lo largo de
su recorrido con el azud El Diquecito y diversos puentes fuera y dentro de la
metrópolis.
Sus afluentes generan con sus crecidas, problemas importantes aguas arriba del
Dique San Roque y el mismo río los genera aguas abajo en su tramo por la ciudad.
43
La presa San Roque es del tipo presa de gravedad y de planta curva, la altura de
su cierre es de 51,30 m, la longitud de coronamiento es de 145 m, su volumen
embalse a cota labio vertedero es de: 201 Hm3, y la superficie del lago a cota de
labio del vertedero es de 1.501 Ha. (Lábaque, M., 2010).
Aguas abajo del lago no existe ningún mecanismo de regulación eficiente. El

Primero nace de la unión de los arroyuelos que descienden de la dorsal de La
Cumbre entre Santa Rosa y el Pungo, por el flanco occidental de la Sierra Chica y
de los que bajan por la cuenca oriental de La Cumbre del Perchel y del Mogote de
Flor, como los arroyos de las Chacras y el de Piedras Blancas. A pocos kilómetros
se reúnen al Oeste de Huerta Grande, en el fondo del valle, y constituyen el río
Grande o de La Punilla, luego cambia por San Francisco, hasta su confluencia con
el Yuspe.
El río Yuspe nace en las pendientes de Los Gigantes, cuyas aguas se encauzan en
los arroyos de Las Vacas Muertas y de la Esquina. A partir de dicha confluencia,
corre hacia el Noreste, atraviesa el valle de Las Cuchillas donde se incorpora el río
Negro. Cruza por el Lavadero, continúa hacia el Este y luego cambia de rumbo
hacia el Nor-noreste. En El Real se suma el río Carnerío (o Carnerillo), que se
inicia en la confluencia del Arroyo Horcosún y de los Piquillines.
El primero de los mencionados, es decir, el Arroyo Horcosùn junta las aguas del
Cerro Ensenada y de la pampa del mismo nombre; recoge al arroyo Quebrada del
Tigre, colector de la Mesa de la Mula Muerta y de las vertientes de la Pampa de las
Ensenadas; corre hacia el Noreste y finalmente toma hacia el Sureste. Luego forma
la Cascada de Olaén, tras lo cual cambia su nombre por el de arroyo de la
Cascada y adelante se une al arroyo de los Piquillines. A unos tres kilómetros de la
afluencia del Carnerío el río Yuspe gira hacia el Este y converge con el San
Francisco. Aquí se inicia el río Cosquín, que transita hacia el Sur, e incorpora a los
arroyos El Rosario y Suncho Huayco o Santa María, y ya al entrar en el lago, al río
Las Mojarras. Éste, escurre una extensa área por medio de sus afluentes, los ríos
Tanti y Mallín.
El sector meridional del Valle de Punilla es avenado por la cuenca del río San
Antonio o San Roque, que nace del encuentro de los ríos Icho Cruz y Malambo. El
Icho Cruz nace en las vertientes de la Pampa de Achala. A la altura de Cerro Bayo,
estas aguas, se reúnen y forman el colector Icho Cruz que se dirige hacia el
Noreste. Tras unos quince kilómetros y donde aumenta su caudal con el aporte de
arroyitos como el de Paso de los Bayos, Paso de la Esquina y de Copina, recoge
las aguas del río de los Sauces, varía su trazo hacia el oriente, rodea el Cerro
antes mencionado entre los dos ríos por el Sudeste y se une al Malambo.
Las aguas de la pendiente oriental y sudoriental de los Gigantes, son reunidas por
el río Malambo, en el que se vierten torrentes y arroyitos como El Mogotes, El
Sauce, Las Barrancas, etc. Corre hacia el Sudeste y concurre con el Icho Cruz
para formar el río San Antonio. El nuevo colector se dirige con rumbo al Este
44
recogiendo las aguas de pequeños arroyos, tales como el de las Achiras, Cuesta
Blanca, etc., destacándose entre los caudalosos el río del Cajón, llamado así por
atravesar una profunda garganta. Desde esta afluencia toma hacia el Noreste. En
Santa Rosa incorpora las aguas del flanco Norte de la Sierrita o Cordón de
Santiago por intermedio del arroyo llamado también San Antonio. Luego transita al
Norte unos ocho kilómetros y alcanza el Lago San Roque.
Atraviesa Carlos Paz y recibe al río de Los Chorrillos, que escurre las pendientes
orientales de la Pampa del Matadero entre los cerros Blanco, Chorrillos y Potrero.
Reunidos en el valle, los ríos Cosquín y San Antonio dan nacimiento al Suquía, que
luego de su travesía por el Dique San Roque, circula hacia la localidad de La
Calera. Aguas abajo del Dique Mal Paso, parten los dos Canales Maestros de
distribución de agua de riego, y recoge al arroyo Saldán. Luego entra en la llanura
y su valle se ensancha hasta alcanzar casi cuatro kilómetros en la ciudad de
Córdoba, con un cauce de aproximadamente cien metros de ancho, pero hoy
reducido a un estrecho canal de hormigón. En pleno centro urbano incorpora, por el
Sur (margen derecha), al arroyo de La Cañada, que desagua el área de La
Lagunilla. (Vázquez, Miatello, & Roqué, 1979)
4.2.2 Evaluación de Antecedentes – Cuenca del Río Suquía
Para un mayor detalle de la Cuenca Suquía y para explicar algunos factores que
influyen en las variables de análisis, reviste particular valor una adecuada y
extensa apreciación de las referencias históricas. Por ser esta cuenca una de los
subsistemas de la Cuenca de Laguna Mar Chiquita tiene un extenso interés
hidrológico y ambiental, tanto a nivel nacional como internacional, por lo que los
antecedentes existentes, aunque dispersos, son de considerable extensión. Este
capítulo se dedica a la recopilación y análisis de información relativa a la
mencionada cuenca, profundizando el aspecto hidrogeológico y los relacionados al
mismo.
A continuación se enumeran los datos y antecedentes en que se basan los

siguientes ítems.
4.2.3 Datos y Antecedentes.

Es importante destacar que salvo algunos estudios puntuales, en la actualidad no
se realizan estudios hidrogeológicos a nivel regional, ni controles por parte de las
reparticiones encargadas del área. (Fulginiti F., 2013)
Los estudios que se destacan, principalmente en el área de hidrología, fueron

realizados por organismos nacionales y provinciales como:
• Consejo Federal de Inversiones (C.F.I.) en el año 1983,
• Universidad Nacional de Córdoba (U.N.C.) en el año 1998 y durante el periodo

2000-2003,
45
• Universidad Nacional de Santiago del Estero (U.N.S.E.) en el año 1997,
• La Agencia Córdoba Ambiente S.E. y la Agencia Córdoba Ciencia que

cofinancian proyectos de investigación desde el año 1998 al igual que la
Secretaría de Ciencia y Tecnología (Se.C.yT.) y la Dirección Provincial de Agua
y Saneamiento (Di.P.A.S. hoy Ministerio de Agua, Ambiente y Servicios
Públicos), a través de Convenios con la U.N.C. desde el año 2000.
• El INTA realiza controles de calidad de agua de las perforaciones existentes en

la llanura pampeana utilizadas para riego complementario.
Los estudios y publicaciones que se consultaron en la presente recopilación son

entre otros:
• “Estudio Hidrológico – Ambiental del Sistema Mar Chiquita - Bañados y Tramo

Inferior del Río Dulce”. Informe Final UNC-UNSE a SSRH (Rodríguez, A. y otros,
2007.)
• Recursos de aguas Subterráneas y su Aprovechamiento en la Llanura

Pampeana y en el Valle del Conlara (C.A.A.A.S., 1975).
• Datos obtenidos del freatímetro ubicado en las proximidades de la laguna.
• Datos de perforaciones registradas en DIPAS.
• Evaluación Hidrogeológica - Diagnóstico Ambiental de Las Aguas Subterráneas

de la Ciudad de Córdoba (Ferral A.;Zappino R.;Alaniz E.y Cheniuk S., 2007).
• Tesis de maestría en Ciencias de la Ingeniería (Plencovich, 2010). “Simulación

Hidrológica de los humedales de Mar Chiquita con apoyo satelital”.
• Informe para el ENOHSA “El agua potable en la República Argentina -

Diagnóstico de fuentes de agua en algunas localidades” (Pflüger J. C., 2007).
 Modelación Hidrogeológica en Cuencas Endorreicas - Aplicación a la Cuenca

de la Laguna Mar Chiquita y la Laguna Llancanelo. Tesis de Maestría. (Fulginiti
F., 2013)
4.2.4 Historia y Problemática Sistema Laguna Mar Chiquita
Antes se explica brevemente la historia y los temas ligados a la Cuenca de la

Laguna Mar Chiquita, del cual la cuenca del río Suquía forma parte.
En los últimos veinte años, se han ejecutado simulaciones hidrológicas de este

sistema constituyendo, algunas de ellas, el inicio para esbozar escenarios más
complicados y edificar las bases de futuras políticas de manejo de los recursos
hídricos.
46
La cuenca de la Laguna Mar Chiquita y Bañados del Río Dulce alcanza cinco
provincias de Argentina y son las siguientes: Salta, Catamarca, Tucumán, Córdoba
y Santiago del Estero. Las tres primeras provincias conforman la cuenca alta que
recolectan sus aportes hídricos en la presa de Río Hondo ubicado en el límite entre
la Provincia de Tucumán y la Provincia de Santiago del Estero. A partir de este
dique nace el Río Dulce que escurre (formando humedales o bañados en épocas
de excedentes hídricos) hacia la Laguna Mar Chiquita, que recibe finalmente el
agua de toda esta cuenca endorreica. (Plencovich, 2010)
Los cambios atmosféricos que desde hace unos 25 años afectan al país y que
quebrantaron el ecosistema de la laguna Mar Chiquita, no sólo acarrearon
aspectos negativos como las inundaciones; además esto mejoró los escenarios de
fertilidad de los terrenos. Anteriormente predominaban las salineras y hoy los
sembrados sustituyeron a los bosques de chañares. El gran espejo de agua
favoreció la presencia de más lluvias y optimizó la calidad de la tierra. Se
aumentaron las aves playeras.
Las constantes oscilaciones, impactaron directamente tanto en el ecosistema

regional, como en el adelanto social y económico de las sociedades asentadas en
las proximidades de la Laguna y en la superficie de los humedales.
La Laguna sufrió impresionantes cambios en su nivel de agua; a partir de los años

’70, se registró una evolución significativa pasando de un área de 2000 km2 a 6500
km2 en el año 2003, cuando se registró el mayor nivel en la costa Sur frente a
Miramar. Esta variación representó un aumento mayor al 300% para dicho período.
La zona es un lugar de paso de aves que emigran desde Canadá hasta el sur del
país. Hoy por hoy la zona transita por una etapa de abundante desarrollo turístico,
aunque las crecidas registradas provocan temor y paran las inversiones por
desconocerse la cota máxima que puede alcanzar la laguna.
Históricamente el mayor desastre se inició a mediados de 1978 cuando la laguna

empezó a elevar su nivel. En 1985, 37 manzanas edificadas quedaron bajo las
aguas, se perdieron 12 Hectáreas de construcciones que incluían el 90 por ciento
de las instalaciones turísticas, de edificios públicos y viviendas. Unos 60 negocios y
200 casas de familia fueron devastados. Para revertir esa situación y bajo el plan
llamado Recuperación de Miramar, el Comando del Tercer Cuerpo del Ejército se
encargó en 1992 de derrumbar con explosivos todas las ruinas que afloraban sobre
la superficie del agua, quitar los escombros de la ciudad sumergida y así devolver
el horizonte limpio que atrajo a turistas de todo el mundo en décadas anteriores.
Los humedales que conforman la zona de bañados, paleocauces y lagunas

dispersas comparten una propiedad primordial: el agua. Esta variable juega un rol
fundamental en el ecosistema, en la determinación de la estructura y las funciones
ecológicas del humedal haciendo que los mismos tengan características diferentes
de los ecosistemas terrestres, como lo es la gran variabilidad espacio - temporal.
47
Esto tiene efectos muy importantes sobre la diversidad biológica que habita en los
humedales, la que debe desarrollar adaptaciones para sobrevivir a estos cambios
que pueden llegar a ser muy extremos, por ejemplo, ciclos hidrológicos de gran
amplitud con períodos de gran sequía y períodos de gran inundación.
En este contexto surge la necesidad de un sistema de planificación multisectorial,

que considere el incremento sustentable de la región y la conservación de la
biodiversidad.
Una de las herramientas de mayor difusión en la que se apoya el proceso de

planificación es la simulación. Los actuales estudios y modelos se limitan al cuerpo
lagunar y a los bañados, independizándose de los procesos en el resto de la
cuenca. Esto permite simplificar los estudios limitando la extensión de los modelos
de funcionamiento propuestos. Sin embargo, el avance en el conocimiento de las
distintas variables y sus interrelaciones obliga a extender las fronteras de
simulación incorporando nuevos elementos. Uno de los elementos de difícil
valoración tanto a escala local como regional es el componente hidrogeológico. La
presente tesis trata de aportar una herramienta adicional a los actuales modelos de
balance, que permita justificar y orientar las decisiones de los organismos
intervinientes. (Fulginiti F., 2013)
4.2.5 Clima
Entre la variedad de climas que presenta la República Argentina (Ilustración 4-3), la

gran extensión de la cuenca superficial de la laguna Mar Chiquita cubre cuatro de
los diferentes tipos en que el territorio fue dividido por el INTA (http://inta.gob.ar/).
Estos son:
1°) Climas cálidos: Dentro de los denominados climas cálidos, la cuenca se

encuentra ubicada en una región que abarca las siguientes variedades:
Clima cálido tropical con estación seca: Se extiende por la llanura chaqueña central
y occidental, abarcando el Oeste de Formosa y Chaco, el Este de Salta y
Tucumán, casi todo Santiago del Estero y el Noroeste de Santa Fe.
Se identifica por las elevadas temperaturas, con pequeñas oscilaciones mensuales

y anuales, y precipitaciones abundantes superiores a 1000 mm anuales. La
temperatura media anual es superior a 20° y no existe ningún mes con promedio
inferior a 15°. Los ciclos fríos no son prolongados,pueden registrarse heladas entre
mayo y agosto.
Hay tres meses con temperaturas medias superiores a 25°, y en enero se registran
temperaturas máximas superiores a 40°. Hay altas temperaturas de verano que
coinciden con la estación de las lluvias, en cambio los inviernos son secos, y en
ellos se registran temperaturas bajas, hasta por debajo de 0°.
48
La insolación es fuerte durante las horas del día, mientras que las noches son
frescas. Esta sequedad se nota en toda la región chaqueña, pero especialmente en
Santiago del Estero. Hay pues una fuerte insolación diaria, relacionada con la gran
diafanidad del aire, sobre todo en el verano, pues generalmente es baja la
humedad relativa; de ahí las fuertes amplitudes térmicas. Entre mayo y septiembre
son frecuentes las heladas.
Los vientos más comunes son los del Norte, Nordeste y Este, y luego los del Sur;
estos últimos provocan los cambios de tiempo acompañados por fuertes descargas
pluviales. Las calmas suelen ser prolongadas durante el verano.
La característica principal es la existencia de un período seco, que abarca casi la

mitad del año, durante el cual los ríos disminuyen su caudal, desciende el nivel de
las aguas subterráneas y la vegetación denota la aridez. Sobre los 539 mm de
lluvia anual que recibe al año Santiago del Estero, entre los meses de mayo y
septiembre cae solamente el 10%.
Ilustración 4-3. Clasificación climática en la llanura chaco pampeana (Fuente: INTA, 1995
citado en Revista “Visión Argentina”. Publicación N°5).
Clima cálido tropical serrano. Comprende las sierras Subandinas del noroeste, los
valles bajos que conducen a las quebradas de Humahuaca y del Toro, y la ladera
oriental del Aconquija.
Las condiciones climáticas varían con la altitud y la dirección de las sierras,

quebradas y valles. Las cuestas que miran hacia el Este y Sur son las más
lluviosas, y la época de las precipitaciones corresponde al verano, coincidiendo con
la culminación del Sol en el trópico de Capricornio, y con el alcance del aire cálido y
húmedo del anticiclón del Atlántico Sur, cuando por su aridez el Noroeste de la
Argentina se convierte en una zona de baja presión.
49
Las laderas serranas que miran al Oeste y Norte, contrastan de las anteriores por
su aridez, en particular las más altas crestas occidentales. Los valles bajos, como
el del río San Francisco, tienen un clima tropical, con un período de heladas
sumamente breve, reducido en general al mes de julio. En los valles algo más altos
se hallan las condiciones típicas del clima tropical serrano.
El régimen pluviométrico sigue siendo de tipo tropical; las zonas de pie de monte
se confunden con las precipitaciones del chaco occidental, con promedios anuales
dentro 600 y 1200 mm anuales, de los cuales el 80% corresponde al verano, entre
noviembre y marzo. El período de mayo a septiembre es netamente seco, pero en
él tienen lugar las heladas, favorecidas por la transparencia del aire, a lo cual
contribuye también el pronunciado enfriamiento nocturno.
2°) Climas templados: Son los climas que reúnen las mejores condiciones para el
desarrollo de los cultivos, que incluyen en el país el Sur de la Mesopotamia y la
llanura pampeana oriental. Dentro de esta extensa región, se notan algunas
diferencias que permiten establecer tres variedades: templado pampeano o de
transición, templado con influencia oceánica y templado serrano. A la cuenca en
estudio le corresponden estas variedades, salvo el templado con influencia
oceánica.
Clima templado pampeano o de transición: Comprende la región pampeana

oriental. Es el más extenso y característico de los climas templados de la
Argentina. Se destaca por los veranos calurosos e inviernos templados, con olas
de aire frío, por su posición abierta a las invasiones de aire polar antártico. Las
estaciones intermedias no están bien definidas, o bien duran poco tiempo, siendo
la más estable el otoño.
El período caluroso comprende los meses de noviembre a marzo. El invierno se

inicia a fines de mayo y dura hasta agosto, con promedios mensuales superiores a
10°, pero se registran con frecuencia temperaturas inferiores a 0°. Los días de mal
tiempo son frecuentes al comienzo del período invernal, pero no faltan días de sol,
y noches de intenso frío.
La elevada humedad de ciertos días de verano, hace asfixiante el ambiente. El

coeficiente de humedad relativa experimenta una resta constante hacia el interior.
Las épocas de menor quietud atmosférica son la primavera y el verano. A fines del
verano y en el otoño, el tiempo es más estable. En noviembre y diciembre, con los
nacientes calores fuertes, tienen lugar tormentas violentas, con brusca caída
termométrica, vientos huracanados y diluvios. Los vientos más frecuentes son los
del norte, nordeste y sudeste, siendo escasos los del noroeste y oeste. Las
quietudes son raras, sobre todo en verano.
Los vientos locales característicos son el Norte, el pampero y la sudestada. El

viento “norte”, cálido y húmedo, trae tiempo caluroso e inestable, que origina los
50
días sofocantes del verano y las bajas presiones. El “pampero” sopla del Sudoeste,
a veces con violencia, y hace posible, al provocar los cambios de temperatura, las
condiciones normales para la vida humana. La atmósfera posee un bajo coeficiente
de diafanidad, y el promedio de heliofonía es más alto en verano que en invierno.
El total anual de lluvias, se halla próximo a los 1000 mm a lo largo de la costa del
río de la Plata. Las lluvias disminuyen paulatinamente hacia el oeste.
En Córdoba hay un período seco entre junio y septiembre, en el cual cae

solamente el 8% del total anual disminuyendo el mismo hacia la faja ribereña
donde prácticamente desaparece.
Circundando por el Oeste al clima templado de transición, de la pampa oriental, se

nota en una franja un aumento creciente de la aridez y una mayor oscilación en las
temperaturas, que señalan el pasaje paulatino a la zona árida del Oeste, con sus
travesías arenosas. Este clima es el existente en la Cuenca del Suquía en su
mayor parte.
Clima árido de las sierras y llanos inter montañosos: Al este de los Andes Áridos se
extiende una gran zona de clima árido, que coincide aproximadamente con la
región de las sierras Pampeanas, con excepción de las sierras de Aconquija,
Córdoba y San Luis. El rasgo saliente de este tipo de clima es la aridez, pero la
presencia de las sierras determina variantes climáticas acentuadas, y se puede
distinguir un clima de los llanos, un clima de las sierras y un clima de los valles
altos. Fuente: Revista “Visión Argentina”. Publicación N°5. (Fulginiti F., 2013)
4.2.6 Temperatura Media Mensual
A continuación se presentan los valores registrados en la zona de interés.
Las tablas siguientes Tabla 4-3 y Tabla 4-4 presentan los valores de temperatura,
precipitaciones y velocidad del viento registrados en las estaciones Córdoba
Observatorio Lat.-31.24 Lon:-64.11 y Córdoba Aeropuerto Lat.-31.24 Lon:-64.11
Mes Temp Temp mín Temp máx Media de Mínima de Máxima de Media
Media Media Media precipitaciones precipitaciones precipitaciones mensual
mensuales mensuales mensuales de la
velocidad
del viento
°C °C °C mm mm mm m/sec
Enero 24.1 18.1 31.1 121.7 24.2 311.5 5.2

Febrero 23.1 17.4 30.1 99.8 20.8 169.9 4.9
Marzo 20.9 15.6 27.6 110.3 5.9 324.5 4.7
Abril 17.9 12.3 24.9 52.2 1.1 263.5 4.3
Mayo 14.9 9.3 22.0 18.9 0.4 79.5 4.4
Junio 11.3 5.7 18.5 11.4 0.0 58.1 3.9
51
Julio 11.3 5.5 18.6 12.8 0.0 111.4 4.8

Agosto 13.2 6.7 21.0 9.7 0.0 40.0 5.8
Septiembre 15.6 9.1 23.3 33.8 0.0 161.9 6.7
Octubre 18.9 12.6 26.1 66.4 11.3 192.7 7.4
Noviembre 21.3 15.2 28.4 96.6 29.2 224.9 7.3
Diciembre 23.2 17.3 30.3 136.9 33.3 286.8 6.4
Tabla 4-3 Meteorología habitual en CORDOBA OBS. Lat.-31.24 Lon:-64.11. Período 1961-
1990. (Fuente: (www.freemeteo.com))
Mes Temp Temp mín Temp máx Media de Mínima de Máxima de Media
Media Media Media precipitaciones precipitaciones precipitaciones mensual
mensuales mensuales mensuales de la
velocidad
del viento
°C °C °C mm mm mm m/sec
Enero 23.5 17.1 29.8 133.3 45.0 302.4 10.6

Febrero 22.5 16.5 28.9 109.7 20.2 209.7 10.1
Marzo 20.3 14.8 26.6 116.1 19.1 257.3 9.8
Abril 17.2 11.5 24.2 51.1 4.2 199.8 9.6
Mayo 14.1 8.5 21.2 21.5 0.3 133.4 10.1
Junio 10.5 4.7 17.9 11.8 0.0 48.9 9.7
Julio 10.4 4.4 17.9 14.3 0.0 124.4 10.7
Agosto 12.3 5.6 20.3 11.4 0.0 42.4 12.5
Septiembre 14.9 7.9 22.4 38.2 0.0 153.3 13.8
Octubre 18.2 11.4 25.2 73.1 6.1 187.9 14.4
Noviembre 20.8 14.2 27.2 107.3 33.7 260.5 13.9
Diciembre 22.7 16.3 29.0 148.7 31.5 296.1 12.8
Tabla 4-4 Meteorología habitual en CORDOBA AERO. Lat.-31.24 Lon:-64.11. Período
1961-1990. (Fuente: (www.freemeteo.com))
4.2.7 Precipitación en la cuenca
Dado que la precipitación constituye una de las principales fuentes de recarga de

los acuíferos, el análisis de las estadísticas de precipitación toma vital importancia
en la estimación de los aportes al sistema subterráneo.
En la Ilustración 4-4 se observa el mapa de isohietas medias anuales para la

República Argentina, realizado por el INTA en el marco del proyecto de
ecorregiones.
El mapa de isohietas evidencia la relación entre la precipitación y la topografía,

aumentando la intensidad de las mismas en la proximidad de la zona costera y
ante ciertos grupos montañosos que influyen en la generación de precipitaciones
orográficas.
52
Ilustración 4-4 Isohietas medias anuales República Argentina (Fuente: http://inta.gob.ar/).

4.2.8 Evapotranspiración
La evapotranspiración, es definida como la cantidad de agua que es transpirada

por la vegetación y evaporada desde la tierra. Esta depende de múltiples
componentes como lo son el tipo de vegetación, las condiciones climáticas y la
clase de suelo entre otros. Su medición directa implica sustanciales dificultades,
por lo que en la generalidad de los casos se recurre a fórmulas experimentales que
la correlacionan con variables de más fácil determinación. La cantidad de
53
elementos que se consideran en cada método, se relacionan por lo general con la

precisión obtenida en los mismos.
Todos los valores disponibles en las fuentes consultadas, hacen uso de este tipo
de metodologías para la determinación de la evapotranspiración, por este motivo
dichos valores se asocian al período de tiempo al que corresponde la información
de base empleada para el cálculo.
La evapotranspiración, es un elemento primordial en la determinación de la recarga

por precipitación al acuífero.
4.3 Suelos
Se deduce por suelos, a la primera capa de aproximadamente 0.5 m de material no

consolidado, ubicado inmediatamente por debajo de la superficie y que se
aprovecha como medio natural para el desarrollo de la vegetación. La magnitud de
la infiltración y escorrentía producto de las precipitaciones está directamente
relacionada con las características del suelo. ( (Ci66j, Apunte de Cátedra, 2006))
Someramente gran parte del suelo que cubre los depósitos aluviales de los ríos
que desembocan en la laguna Mar Chiquita, corresponden a depósitos de tipo
loésico.
La zona ocupada por loess en la República Argentina se extiende en las planicies

situadas desde 23° a 38° S, cubriendo más de 600,000 km2 (Ilustración 4-5). Hacia
el Oeste, está contenida por varias cadenas montañosas que ocupan una faja que
se extiende en dirección N-S, a lo largo de Argentina. Recientemente se ha
propuesto la subdivisión entre loess pampeano (al sur de 30° S) y loess chaqueño,
subtropical (al norte de 30° S). Hay otros depósitos loésicos en valles montañosos
y planicies de altura del Oeste del país.
En el caso del Tipo de loess Pampeano, predominan partículas de origen volcánico

originadas en los Andes Centrales y en la Patagonia. Dichas partículas fueron
transportadas por agentes erosivos tales como glaciaciones y condiciones fluviales
hacia el norte de la Patagonia. Desde esta región las partículas fueron
transportadas por la acción de los vientos del S - SW hacia el centro del país.
(Rocca R. J., 2005)
54
Ilustración 4-5 Distribución de Loess y Loessoides (Fuente: (Rocca R. J., 2005), adaptado
de Zárate, 2003). ( (Fulginiti F., 2013))
Los suelos loésicos exhiben características en cuanto a la anisotropía de la

conductividad hidráulica, la conductividad hidráulica vertical (Kv) es mayor que la
horizontal (Kh). Los niveles superiores de limos arcillosos, presentan
estructuralmente una red de fisuras y microfisuras visibles. Esta red incrementa el
valor de la conductividad hidráulica.
A menor escala y focalizando en las provincias que conforman la mayor parte de la

cuenca, se presenta la clasificación de sus suelos según el sistema "Soil
Taxonomy" (SSS- USDA, 1975) representado en el "Atlas de Suelos de la
República Argentina" elaborado por el INTA. Esta clasificación hace posible la
subdivisión sistemática en diferentes tipos de suelos, caracterizados por sus
propiedades específicas e identificadas por nombres de clases.
La clasificación taxonómica según el sistema norteamericano denominado "Soil

Taxonomy" (SSS- USDA, 1975), muestra un patrón de relaciones genéticas entre
los suelos y su estructura bajo un esquema de seis categorías o niveles sucesivos.
55
Los niveles superiores poseen pocas divisiones y se encuentran definidos por

características generales, a medida que se progresa hacia los niveles inferiores las
divisiones aumentan y las propiedades son definidas en forma específica por un
gran número de propiedades del suelo.
En el orden de mayor a menor nivel, las categorías que conforman el sistema de

clasificación "Soil Taxonomy" son: 1) Orden; 2) Suborden; 3) Gran Grupo; 4)
Subgrupo; 5) Familia y 6) Serie
El plano siguiente presenta los órdenes y subórdenes de suelos, según este criterio
para la provincia de Córdoba (Ilustración 4-6). El mencionado gráfico junto a los
planos geológicos nos podrán dar un mejor concepto para tener una mejor
definición de las propiedades hidráulicas superficiales al instante del desarrollo del
modelo.
Ilustración 4-6 Suelos en la provincia de Córdoba (Fuente: http://inta.gob.ar/ - Elaborado

en base al Atlas de suelos de la Rep. Argentina). ( (Fulginiti F., 2013))
Los Molisoles ocupan áreas importantes en la planicie Chaco-Pampeana y

constituyen los suelos dominantes, entre los que tienen la mejor aptitud para la
agricultura.
La Región Pampeana, tanto húmeda como semiárida, se caracteriza

respectivamente por Udoles y Ustoles, con ocurrencia menor de Acuoles en áreas
planas utilizadas para la producción de ganado.
56
Los Molisoles son suelos pardos o negros, que se han desarrollado a partir de
sedimentos minerales en climas templado- húmedo a semiárido, no obstante
también pueden presentarse en climas fríos y cálidos cubiertos con una vegetación
de gramíneas.
Las propiedades que caracterizan a los Molisoles son: la presencia de un epipedón

mólico, por la incorporación de residuos vegetales y su mezcla con la parte mineral,
que con el paso del tiempo produce el oscurecimiento del suelo por la
incorporación de materia orgánica; la estructura granular o migajosa moderada y
fuerte que facilita el movimiento del agua y del aire; el dominio del catión calcio en
el complejo de intercambio catiónico, que favorece la fluctuación de los coloides; la
dominancia de las arcillas, moderada a alta capacidad de intercambio y elevada
saturación con bases.
Los suelos desarrollados en las planicies de inundación, derrames, deltas de ríos y

sobre sedimentos depositados recientemente por las aguas son afectados con
frecuencia por inundaciones pero no llegan a saturarse con agua.
En el perfil de suelos es común la presencia de materiales provenientes de la

erosión de áreas altas conteniendo una cantidad considerable de materia orgánica.
La característica común a todos los órdenes fluventes, es la ausencia de rasgos
asociados a la humedad, los que están presentes en suelos más mojados como los
Acueptes.
El suelo de la región de bañados y costas de la Laguna Mar Chiquita es

francamente salino, particularmente en áreas bajas, exhibiendo eflorescencias de
sales, a veces de gran extensión, sobre todo en periodos secos. El carácter salino
de esta cuenca no parece ser el mismo que el de las Salinas Grandes, si se mira
las diferencias de vegetación.
Tratándose de una región de clima seco abundan los suelos salinos en general
intrazonales; las Salinas Grandes y la cuenca cerrada de Mar Chiquita representan
los tipos más importantes, en ambas el suelo es por demás salino, sin vegetación y
con eflorescencias de sales, hasta poco salino con monte; como se trata de
depresiones, son suelos pobremente drenados y de constitución en general areno -
arcilloso con elevada cantidad de sodio intercambiable distribuido en el perfil.
El suelo en general es limo – arenoso, con mayor permeabilidad cuando hay mayor
pendiente, mientras que en las depresiones predomina la arcilla aunque no son
raros aquí los lugares completamente arenosos.
La parte oriental es algo salina y más húmeda que la occidental, con mayor
salinidad en las cercanías del escalón (Borde de Los Altos).
El suelo es pardo - rojizo, desarrollado sobre el loess o material limo - loésico. Su

relieve es plano a plano - cóncavo con drenaje moderado ( (Rodríguez, A. y otros,
2007.))
57
4.4 Cobertura del Suelo
La cobertura del suelo influye de forma directa, en la evapotranspiración como en

la capacidad de infiltración del terreno. Los principales tipos de cobertura
(Ilustración 4-7) en la macro cuenca que incluye la cuenca del Río Salado en
Santiago del Estero según información obtenida del SIG (sistema de información
geográfica) del IGN. ( Web www.ign.gob.ar-Instituto Geografico Nacional de la
Republica Argentina).
TIPO USO DE TIERRA MIXTO: Área con cultivos diversos como cereales,
pasturas y silvicultura alternados.
TIPO PAJONAL, JUNCAL, MALEZAL, MATORRAL Y JARILLAL: Terreno cubierto

con vegetación de distinta clase como el pajón, totoras, juncos, etc. en general se
trata de crecimientos en terrenos húmedos.
TIPO BAÑADOS: Terreno anegado que contiene agua durante todo el año. Esta
zona es de piso generalmente fangoso, pantanoso o cenagoso, cubierto de agua
de poca profundidad, proveniente del desborde de un lago, río, arroyo o por
invasión de las aguas del mar.
TIPO BOSQUE, SELVA, FORESTA, PARQUE NATURAL TRANSITABLE: terreno

cubierto de árboles de diferentes especies perennes y de no más de tres o cuatro
metros de altura.
TIPO BOSQUE, SELVA, FORESTA, PARQUE NATURAL INTRANSITABLE:

Terreno cubierto de árboles de diferentes especies perennes y de más de cuatro
metros de altura.
TIPO PLANTACIONES PERENNES: Área cubierta por siembra sistemática de

vegetación persistente, como ser: árboles, yerba mate, café, té y frutales, etc.
TIPO AFLORAMIENTO ROCOSO: Manifestación superficial de estratos, filón o

masa de cualquier tipo de roca que asoma a la superficie del terreno. Pueden
afectar áreas más o menos grandes.
58
Ilustración 4-7 Cobertura del Suelo en la macrocuenca de Mar Chiquita (Fuente:

Elaborado en base a información del SIG del IGN. ( Web www.ign.gob.ar-Instituto
Geografico Nacional de la Republica Argentina).
59
4.5 Consideraciones Hidrológicas
El área de estudio se encuentra en la región central de Argentina, en la cuenca

hídrica del Río Suquía, provincia de Córdoba.
El río Suquía a la altura de la ciudad de Córdoba, presenta un módulo de 9,7 m³/s

con un valor máximo de 24 m³/s y un mínimo de 2 m³/s. La pendiente media del
cauce varía desde 6 m/km (6 ‰) a 3 m/km (3 ‰), en la zona de desembocadura en
la Laguna Mar Chiquita.
Ilustración 4-8 Imagen satelital Río Primero – Suquía (Fuente: Elaboración propia - Imagen
tomada de Google Earth). (Fulginiti F., 2013)
En la Ilustración 4-8 se observa el Río Primero o Suquía y su contexto geográfico.
4.5.1 Balance Hidrológico
La cuantificación del ciclo hidrológico, puede ser representado por una ecuación de
balance de masa, donde las entradas menos las salidas son iguales a los cambios
en el almacenamiento. A pesar de la simplicidad de este concepto, la elaboración
de modelos matemáticos que lo reflejen no resulta fácil, en especial para simular
sistemas de gran escala como es el caso de cuencas de gran extensión y escasa
instrumentación.
Estos modelos necesitan ser nutridos con series de datos hidrometeorológicos, que
necesariamente deben ser medidos in-situ y a su vez es necesario poder contrastar
los resultados obtenidos del modelo con datos reales.
60
La inclusión de la variable hidrogeológica a escala regional, al ampliar las fronteras

del modelo, dificulta los procesos de calibración y validación al incrementarse el
número de datos hidrometeorológicos necesarios para este fin.
4.6 Aspectos Geológicos e Hidrogeológicos
En el territorio de la provincia de Córdoba el perfil geológico generalizado es el

siguiente: Basamento Cristalino, Paleozoico (Pérmico), Cretácico, Terciario y
Cuartario (Schlagintweit O., 1946). El basamento aflorante de las sierras se hunde
hacia el este con una estructura de fosas y pilares. La zona central corresponde a
una fosa que superficialmente se manifiesta de forma homóloga (posibilitando una
cuenca hidrográfica de tipo endorreico).
La estructura en bloques, fosas y pilares del Basamento Cristalino y el afloramiento

del mismo en las sierras, permite distinguir varias cuencas hidrogeológicas, cuyos
límites son un tanto imprecisos por la falta de información del ambiente profundo.
Los cordones montañosos, con un rumbo aproximado Nor noreste, dividen a la
Provincia en tres unidades características: la montañosa, la llanura del este y las
del Noroeste y Oeste. Desde un punto de vista geológico, se encuentran en toda la
región tres formaciones bien diferenciadas, formación Paraná, formación Puelches
y formación Pampeana.
En el sector de las cuencas intermontanas, aparte de las pequeñas acumulaciones

de sedimentos modernos que rellenan depresiones y valles intermontanos, hacia el
Este y bordeando las sierras, hay un conjunto de conglomerados y areniscas rojas
(Cretácico), que afloran en varias localidades, como Saldán, Ongamira,
Despeñaderos, etc. Este complejo, con buzamiento suave al Este, sirve de base
impermeable para las aguas que circulan dentro de los sedimentos más recientes,
que los cubren, al hundirse en la llanura. (Vasquez, Miatello, Roque, & otros, 1979)
Si bien los bloques del basamento, regulan la acumulación y circulación de las

aguas subterráneas, la heterogeneidad de los sedimentos que lo cubren, es la
causante de la distribución de los distintos acuíferos. Se puede decir que hacia el
Este, los acuíferos profundos se manifiestan con mayor caudal, incluso, muchos
son surgentes, pero altamente salinizados.
En profundidad, esta área presenta una estructura fallada de bloques del

Basamento Cristalino, que originaron fosas y pilares. Sobre el mismo se hallan
sedimentitas rojizas continentales con carbonato de calcio y yeso. Por encima se
disponen los materiales sedimentarios marinos de la Formación Paraná (Mioceno)
constituida por arcilitas verdes y arenas presentando un espesor de 60 a 80 m. La
Formación Santiago Temple, de origen continental sólo existe en el área occidental
de la llanura (área pedemontana), siendo reemplazada lateralmente hacia el este
por las Formaciones Paraná y Puelches. Está constituida por sedimentitas arcillo-
limosas calcáreas de coloración rojiza y castaño con capas aisladas de arena
61
gruesa y gravas, presentando espesores variables entre 100 y 400 metros.

Estratigráficamente se pueden sintetizar las principales características de la zona
occidental y oriental de la cuenca en la Tabla 4-5 y Tabla 4-6:
Estratigrafía Hidroestratigrafía Acuíferos Principales

Sedimentos actuales o
Epiparaneano Acuíferos Aluviales
aluvio
Sedimentos Pampeanos Epiparaneano Acuífero Pampeano
Formación Santiago
Epiparaneano Acuífero / Acuitardo
Temple
Formación Saldán Paraneano Acuífero Paraná
Basamento Cristalino Basamento Hidrogeológico Acuífugo
Tabla 4-5 Estratigrafía sector occidental de la cuenca en la provincia de Córdoba
Estratigrafía Hidroestratigrafía Acuíferos Principales

Sedimentos Pampeanos Epiparaneano Acuífero Pampeano
Formación Puelches Epiparaneano Acuífero / Acuitardo
Formación Paraná Paraneano Acuífero Paraná
Formación Arcillas Pardas Hipoaraneano Acuitardo

Formación Serra Geral Basamento Hidrogeológico Acuífero / Acuífugo
Tabla 4-6 Estratigrafía sector central y oriental de la cuenca en la provincia de Córdoba.
En general, el agua subterránea de la Llanura Chaco-Pampeana aumenta su

contenido salino a medida que crece la aridez ambiental. También, incrementa su
salinidad en profundidad, debido a la presencia de sedimentos arcillosos marinos
(depósitos paranianos o Formación Paraná). Localmente, hacia las zonas de
descarga, en ámbitos mal drenados, con predominio de sedimentos pelíticos y
niveles freáticos cercanos a la superficie (zonas de lagunas, esteros o bañados),
también se produce un notorio incremento salino (Laguna Mar Chiquita - Córdoba).
Otro sedimento que comúnmente integra esta vasta llanura aluvial es el loess (limo
– arenoso de origen eólico) que posee permeabilidades más bajas (1 m/día) que
las que presentan las Arenas Puelches, mientras que la salinidad de sus aguas
normalmente está controlada por la profundidad del nivel freático y por los
caracteres climáticos dominantes (en regiones húmedas aguas de buena calidad y
en regiones áridas aguas con elevados tenores salinos).
El Loess Pampeano contiene a la capa freática y a veces, debido a variaciones en

su permeabilidad tanto verticales como horizontales, puede presentar uno o más
acuíferos bajo presión, constituyendo sistemas semiconfinados o un acuífero
multiunitario.
62
4.6.1 Formaciones Geológicas
En la Ilustración 4-9 se puede apreciar un perfil esquemático (Stappenbeck, R,

1965) del subsuelo desde las montañas cordobesas hasta Uruguay. Las
formaciones predominantes y claramente diferenciadas como ya se expresó en el
ítem anterior son la formación Paraná, formación Puelches y formación Pampeana.
Ilustración 4-9 Perfil esquemático del subsuelo de la Pampa, desde las montañas
cordobesas hasta el Uruguay (Fuente: según Stappenbeck, 1965 en Vazquez, 1979).
(Fulginiti, 2013)
Formación Paraná
La formación Paraná es la más profunda de las tres (Paraná, Puelches y

Pampeana). Constituida por sedimentos marinos de facies costeras y deltaicas de
posición cronoestratigráfica no aclarada, variando entre Oligoceno y Pleistoceno.
Consta de arenas de clastos de cuarzo, feldespatos y minerales livianos de
coloración rojiza a verdosa, intercalando con arcillas de color verde oscuro no
calcáreas (Moyano Miguel A., Leynaud José M., 1994).
Por lo general sobre ésta, se encuentra la Formación Puelches y por encima la

unidad más reciente que generalmente forma la superficie terrestre, la formación
Pampeana. La formación Paraná contiene capas de arena cuyo espesor varía
entre los 2 y 5 m aumentando su participación en la formación hacia el este.
Muchos perfiles de suelo de esta región tienen en común la siguiente distribución:
en la parte inferior, encima de las arcillas de la Formación Paraná, la Formación
Puelches (C.A.A.A.S., 1975).
63
Formación Puelches
Sobre la Formación Paraná yacen sedimentos litológicamente muy parecidos, pero

los cuales ofrecen pocos o ningún indicio sobre su origen marino. Predominan aquí
las arenas sobre las capas arcillo-limosas.
La parte superior de esta formación está formada por una capa continua y extensa
de arena, la cual reviste gran importancia como acuífero en vastas regiones de la
pampa argentina. El espesor total de la formación oscila entre 20 y 80 m.
Posiblemente las diferencias de espesor no revistan mayor importancia en el
aspecto hidrogeológico, ya que en lugares de mayor espesor se depositaron más
sedimentos de grano fino (capas de arcilla, limo y arena fina) (C.A.A.A.S., 1975)
Contrariamente a las capas de arena de la Formación Paraná, estas arenas

contienen esencialmente más cuarzo, cuyos granos son más redondeados dentro
del mismo tamaño de granos. Las capas superiores de arena presentan
generalmente un color amarillo claro hasta pardo, debido a que los granos
cuarzosos tienen un recubrimiento limolítico; contienen menos feldespato, pero
más ágata y carniola. Entre las capas arenosas, estrechamente interdigitadas con
ellas, yacen sedimentos limoso-arcillosos de color verde claro y, a veces verde
oscuro.
Estas capas arcillosas no se encuentran en determinados horizontes

correlacionables, sino que forman junto con las capas arenosas un complejo que
difiere de lugar a lugar. Pero, lo que muchos perfiles tienen en común, es que en la
parte inferior encima de las arcillas de la Formación Paraná, la Formación Puelches
comienza con una arena de grano relativamente grueso, luego se desarrolla como
secuencia alternada de arcilla y arena, terminando hacia arriba con una capa de
arena relativamente pura. El cambio hacia la Formación Pampeana a menudo está
claramente marcado por el límite inferior de una capa arcillosa característica de
color amarillo claro, pardo o verde oliva. Pero también existen lugares con una
zona limítrofe en la cual las arenas cuarzosas de la Formación Puelches forman
intercalaciones con las arcillas multicolores (C.A.A.A.S., 1975)
La extensión de la Formación Puelches coincide casi completamente con la de la

Formación Paraná y sobrepasa tan solo un poco el límite occidental de esta última.
Como característica diferencial con la Formación Pampeana suprayacente al
Oeste, se registra en todos los casos la existencia de una capa de arena de alto
contenido cuarzoso y marcado redondeamiento de los granos (Facies Puelches).
Formación Santiago Temple
Coetanea a la Formación Paraná (Oligo-Pleistoceno) pero formada en ambiente

continental. Está constituida por rocas arcillosas limosas de color rojo a pardo y
capas aisladas de arenas gruesas y gravas. Esta formación solo existe en el área
occidental de la llanura (piedemonte) siendo reemplazada hacia el Este por las
formaciones Puelches y Paraná. Su espesor es de 100 a 400 m tanto en el techo
64
como en el piso no son claramente distinguibles. (Moyano Miguel A., Leynaud José
M., 1994).
Formación Pampeana
Siendo la unidad más reciente, la Formación Pampeana forma generalmente la

superficie terrestre. (González Bonorino F., 1965), le adjudica a esta formación
todos los sedimentos que yacen encima de la capa superior de la facies Puelches.
La Formación Puelches está cubierta por una capa de arcilla poco permeable, que
forma en esta zona la base de la formación más reciente, la Formación Pampeana.
La Formación Pampeana está compuesta exclusivamente por depósitos eólicos, el
loess. El espesor de la formación Pampeana es de 80 a 120 m.
La capa de arcilla, muy poco permeable en la base de la Formación Pampeana,

provoca una separación en dos cuerpos del agua subterránea: uno inferior en las
arenas de la Formación de Puelches, y otro superior en la Formación Pampeana.
El Acuífero Pampeano se caracteriza como de baja a mediana productividad. La

recarga se realiza por infiltración directa de la lluvia y además de sus propias
características hidrogeológicas, se destaca por constituir la fuente de recarga del
Acuífero Puelche, mediante el proceso de filtración vertical descendente (Auge M.
P. y Hernández M. A., 1984).
4.6.2 Plano Geológico de Superficie
El disponer de información hidrogeológica de cada sector, que constituye la cuenca

en estudio, suele ser una limitante en la generación de modelos numéricos. Por
este motivo los mapas geológicos de superficie suelen resultar de suma utilidad
para establecer las bases de las condiciones hidrogeológicas, debido a que en los
mismos se representan los afloramientos de las formaciones geológicas y los
depósitos superficiales. Como ya se estableció una formación agrupa a un conjunto
de unidades litológicas semejantes. Por ello, con la base de un mapa geológico, un
reconocimiento de campo y datos de perforaciones, se pueden establecer las
propiedades hidráulicas de dichas unidades y así identificar rápidamente los
ambientes de mayor interés hidrogeológico.
Desde el punto de vista hidrogeológico interesa conocer el comportamiento de las

unidades aflorantes en función de su capacidad para almacenar y transmitir agua.
Por ello es necesario destacar que un mapa hidrogeológico puede diferir
sustancialmente del geológico de la misma zona. El mapa geológico puede agrupar
en una misma formación dos coladas de lava superpuestas, tengan o no alvéolos y
estén o no fracturadas. El mapa hidrogeológico deberá diferenciar aquellos
basaltos integrados por coladas masivas, de los que presenten oquedades o
alvéolos comunicados y fracturas o diaclasas, pues pueden darle elevada
permeabilidad secundaria.
65
A escala regional, las propiedades diferenciales de las mismas formaciones no

suelen revestir gran trascendencia, por lo que el mapa geológico superficial es una
buena aproximación en la definición inicial de las propiedades intervinientes.
Con este fin Fulginiti generó un plano que sintetiza la geología superficial a partir
de las cartas geológicas elaboradas por el Instituto de Geología y Recursos
Minerales de la República Argentina (Segemar).
Para su desarrollo se empleó el siguiente procedimiento:
• Digitalización de las cartas a escala 1:50000.
• Digitalización y ubicación de los perfiles transversales.
• Georreferenciación de las distintas unidades mediante la identificación de puntos

singulares.
• Unificación de las referencias geológicas.
• Identificación de zonas de interés hidrogeológico.
Ilustración 4-10 Carta Geológica del Segemar. (Fulginiti F., 2013)
66
Ilustración 4-11 Plano con la geología superficial de la cuenca (Fuente: elaborado en base
a Cartas Geológicas del Segemar). (Fulginiti F., 2013)
67
4.6.3 El agua subterránea en las Sierras
Las condiciones de drenaje subterráneo en la zona de las sierras están vinculadas

al relieve, al espesor de la capa meteorizada, a las características de las diaclasas
y fisuras de las rocas, a la cobertura vegetal y a las condiciones climáticas. El
sector serrano constituye la zona de infiltración y alimentación de los cauces que
llegan hasta el piedemonte. En el se encuentra un solo acuífero, generalmente
ubicado entre el sedimento o capa meteorizada y la roca del basamento. (Cioccale
M. A., 1998)
Las oscilaciones del nivel de agua subterránea están reguladas por las variaciones
estacionales de las precipitaciones. En la estación lluviosa se produce la saturación
de los suelos, el aporte permite el ascenso del nivel freático y la ocurrencia de
surgencias, al aumento de los caudales de los cauces permanentes y la activación
de los cauces temporarios. Durante la época invernal se seca, con un marcado
déficit hídrico, los ríos y arroyos serranos mantienen su caudal casi exclusivamente
a expensas de las aguas almacenadas durante el verano.
La Ilustración 4-12 muestra en forma esquemática la posición del nivel freático

respecto al relieve. Las quebradas más profundas siempre tienen en sus fondos
cauces permanentes, mientras que los valles situados a mayor altitud y con
pendientes menos pronunciadas presentan cauces temporarios o efímeros.
En las sierras es muy frecuente durante la época estival la ocurrencia de vertientes

fuertemente vinculadas con el relieve, el espesor de la cubierta y las diferentes
estructuras del macizo rocoso. Se ha observado que estas vertientes se producen
en la media ladera (Quintana Salvat F. y Barbeito O. L., 1998), donde el nivel
freático intersecta a la superficie. Esta situación es bastante característica en
aquellas laderas cóncavo – convexas que presentan en su parte media un punto de
inflexión. En este sector, de máxima pendiente se intensifican los procesos
erosivos que dejan al descubierto los afloramientos rocosos impermeables.
68
Ilustración 4-12 Perfil Hidrogeológico en la zona serrana de la provincia de Córdoba.

(Fuente: Mapa Hidrogeológico de la pcia. de Córdoba – 1978 – Dirección Provincial de
Hidráulica). (Fulginiti F., 2013)
4.6.4 El agua subterránea en la zona de piedemonte
En el piedemonte, el nivel freático sufre un marcado descenso, y se ubica muy por

debajo de las corrientes superficiales, aunque en épocas de lluvias incrementa su
69
nivel, rápidamente desciende y la mayor parte del año el movimiento de agua se

efectúa desde el canal hacia el acuífero libre.
El piedemonte presenta varios niveles de aguas subterráneas que constituyen un

sistema de acuíferos de diferentes edades, los cuales aumentan conforme se
incrementa la profundidad y disminuye la velocidad de circulación. La salinización
es el proceso indeseable más frecuente en las aguas subterráneas del piedemonte
y de las sierras.
Sierra Piedemonte
Solo se presenta el acuífero libre Presencia de más de un acuífero
El nivel freático está en contacto y
El nivel freático esta por debajo de los
vinculado con los cauces. El agua
cauces. El agua durante todo el año se
almacenada durante el verano mantiene los
infiltra, los ríos alimentan el acuífero libre.
cauces durante el invierno.
Es común la presencia de vertientes debido
a la geometría de la ladera y las diferencias No existen fuentes emergentes naturales.
de permeabilidad de las rocas
En los movimientos del agua subsuperficial, En los movimientos del agua subsuperficial,
la circulación paralela a la superficie la circulación de ascenso y descenso
domina sobre los movimientos de ascenso (verticales) domina sobre los movimientos
y descenso (verticales) paralelos a la superficie
Tabla 4-7 Expone las principales características hidrogeológicas a nivel general que
distinguen la zona de sierra de la de piedemonte según (Cioccale M. A., 1998). (Fulginiti F.,
2013)
4.6.5 Calidad del agua subterránea
Con respecto al uso del agua subterránea, la misma representa en la Provincia un

recurso natural inapreciable. La explotación para consumo se realiza mayormente
de los acuíferos profundos, que por lo general están aislados del proceso
bacteriológico y protegidos de la contaminación resultante de la actividad antrópica,
que se presenta muy concentrada en los núcleos urbanos. No obstante una gran
cantidad de poblados y particulares se nutren del acuífero freático para riego,
ganado e incluso consumo humano con los riesgos y costos asociados a la
potabilización del mismo.
La explotación sustentable de los acuíferos profundos, requiere un posterior

análisis contemplando las obras de captación, con el fin de evitar que un uso
impropio produzca la infiltración del agua salobre circundante.
El uso excesivo de fertilizantes y pesticidas constituye la causa más importante de

contaminación para los cursos de ríos y acuíferos.
70
A fines del año 2003, se produjo una gran mortalidad de peces en acumulaciones
de agua al Norte de la Laguna (a la altura de la Laguna Los Porongos), a ciencia
cierta no se conoce cuáles fueron las causas pero las dos hipótesis más fuertes
fueron: (a) el derrame masivo de insecticidas de alta toxicidad, provenientes de las
áreas de agricultura que fueron lavadas por una lluvia, en el instante posterior a su
colocación o (b) un nivel muy bajo del río que llega a interrumpir totalmente el flujo
(Plencovich, 2010).
Los elementos inorgánicos diluidos en el agua, proceden de tres fuentes distintas:

la atmósfera, la biosfera y la litosfera.
La precipitación que cae de la atmósfera hacia la tierra, ya contiene sustancias

diluidas. Provienen de la evaporación del mar (especialmente NaCl), son
absorbidas por el aire como elemento de forma gaseosa y transmitida en forma de
iones (HCO3, NO3), o bien son componentes disueltos de partículas de polvo, las
cuales - amén de las cenizas volcánicas - vienen de la superficie terrestre (Ca, K,
SO4, HC O3, SiO2). Cerca de la costa puede llegar a ser más de 30 mg/l, y se
reduce rápidamente hasta 1 mg/l y menos a medida que la distancia hacia la costa
se hace mayor.
Las porciones de iones HCO3, NO3 y SO4 provenientes de los gases aéreos, no
desempeñan un rol importante en el balance hídrico del agua subterránea, ya que,
o bien en la biosfera son muy influenciados en cuanto a su cantidad (HCO 3, NO3), o
bien como partículas de tierra tienen más importancia (SO4). La proporción de
iones en el agua de lluvia, que se han diluido de partículas sólidas de polvo (sobre
todo Ca y SO4), aumenta naturalmente en dirección del mar hacia tierra adentro, y
depende de la vegetación y de la composición de las rocas de la superficie
terrestre.
En la biosfera se gastan (NO3) o se acumulan (HCO3) iones que juegan un papel

en el metabolismo de las plantas y de los animales. En un promedio a través de
largos años, la interacumulación tiene tan solo un papel reducido. El ion nitrato,
originariamente contenido en el agua de lluvia en igual cantidad que el ion cloruro,
en el suelo es absorbido por las plantas y debido a la descomposición
microbiológica de partes vegetales muertas, vuelve al aire como N 2. Por lo tanto los
altos contenidos de nitrato con alto contenido total de sales, no tienen que deberse
forzosamente a contaminaciones por materias fecales humanas o animales, sino
que pueden ser la expresión de una descomposición biológica deficiente del
contenido de nitrato natural bajo condiciones climáticas áridas.
Durante la infiltración de la precipitación al suelo, el contenido de HCO 3 es

enriquecido por el aire del suelo, que con una vegetación activa, contiene mucho
CO2. (Los ácidos humíferos no tienen papel importante alguno en las aguas
subterráneas de la zona de trabajo). El anhídrido carbónico absorbido ocasiona
una reducción del valor PH y con ello la disolución del calcio existente en el suelo,
con lo cual llegan a diluirse aniones adicionales HCO3.
71
Las influencias de la litosfera se generan por causas alternantes del agua

subterránea y del agua de infiltración y de sus elementos contenidos en rocas, que,
como dilución o precipitación, parcialmente pueden describirse como procesos
Redox.
De tales cifras se desprende, a grandes rasgos, que con una total meteorización y
disolución de la roca existirían en sus aguas capilares todos los iones mencionados
menos Na y Cl, en estado completamente sobresaturado.
En las aguas de La Pampa, generalmente ya se alcanzó una saturación en cuanto

a Ca, HCO3 y Fe. Por lo tanto, una meteorización paulatina de los silicatos del
suelo, ya no puede modificar los contenidos en cuanto a esos iones, pero si los
contenidos de iones en SO4, Mg, Na, CL
De estos últimos, el SO4, en muchos casos, se acerca ya a la saturación, y podría

encontrarse en equilibrio con el yeso del suelo. Sin embargo, en la mayoría de los
casos aún, existe una falta de SO4, es decir, el contenido de este ion puede
aumentar todavía mediante hidrólisis y oxidación de silicatos con contenido S.
Empero, puede demostrarse que la concentración de SO4 alcanzada de esta
manera, es muy reducida en las aguas subterráneas profundas; ya que los
oxidantes son aportados por el mismo flujo de agua que se lleva al ion SO 4
formado. Exclusivamente en el área de fluctuación del nivel de agua subterránea,
puede alcanzarse una alta concentración de SO 4, bajo influencia del oxígeno
atmosférico.
Paralelamente a las dependencias de solubilidad de los iones individuales con sus

sales en las aguas subterráneas, desempeñan un rol importante los procesos de
absorción de minerales arcillosos, recristalización y la nueva formación de
minerales arcillosos. Estos procesos se refieren ante todo a los cationes Ca, Mg, K,
Na y ocasionan modificaciones de sus contenidos relativos, siendo así que Na es el
que menos, y Ca el que más se ligan.
En cambio apenas puede apreciarse la porción de iones procedentes de rocas en

el agua subterránea, en comparación a aquellos provenientes de otras fuentes.
Empero, puede establecerse que una constante meteorización, influiría de la

siguiente manera sobre las existencias de iones del agua subterránea en la zona
comprendida desde la Laguna de Mar Chiquita al Sur:
• Ca, HCO3 y Fe están mayormente saturados y se modifican en dependencia

recíproca.
• SO4 sólo excepcionalmente está en equilibrio con el yeso de los suelos;

mediante incremento de la meteorización en superficie, o por procesos de
concentración, puede aún subir el contenido en el agua subterránea freática,
pero no así en acuíferos confinados.
72
• Pero como las relaciones de los iones entre sí (Mg con Ca, Cl con SO 4) en el
agua del cuerpo inferior de agua subterránea, apenas difieren de aquellas
del agua subterránea freática, un tiempo de estadía prolongado dentro de un
acuífero confinado, tendría - en su totalidad - tan sólo un efecto reducido
sobre su reserva de elementos.
• En el agua subterránea freática estancada, pueden enriquecerse los cationes

Mg, Na, K, a raíz de una larga meteorización; pero predomina aparentemente
la cantidad aportada por la precipitación, ya que cloruro también aumenta en
la misma relación. (Fulginiti F., 2013)
4.7 Cuenca del Río Milagro. Milagro-Ecuador
En el Ecuador se estudiará la zona acuífera que rodea a la ciudad de Milagro, en la

Provincia del Guayas.
Se ubica en la cuenca del río Guayas, que es la más grande de Ecuador. La
Cuenca del Río Guayas, tiene una extensión de 34.500 km², el 13 por ciento del
territorio nacional del Ecuador, y en esta habita el 40 por ciento de la población
total del país.
La cuenca del Guayas tiene las siguientes Subcuencas:
1. Subcuenca Río Daule
2. Subcuenca Río Macul
3. Subcuenca Río Vinces
4. Subcuenca Rio Babahoyo
5. Subcuenca Río Jujan
6. Subcuenca Río Yaguachi
La subcuenca del río Milagro, es parte de la Subcuenca Río Yaguachi que posee
una extensión de 316.4 Km².
73
Ilustración 4-13.- Cuenca del Río Guayas-Ecuador (CELEC, 2012)
El promedio de lluvias que cae anualmente en la región va desde el extremo de

3000 mm. , en el límite Noreste (Santo Domingo de los Colorados), disminuyendo a
2.200 mm., en Quevedo, 1.800 mm., en Daule-Peripa, y 1.200 mm., en Guayaquil,
750 mm. ,en Chongón, hasta llegar a los 50 mm., en el extremo occidental de la
Península de Santa Elena, en Salinas.
En lo que respecta al temporal lluvioso en la región, vale destacar que el 85 por

ciento de las precipitaciones se registran en los cuatro primeros meses de año,
durante el invierno ecuatorial, mientras que el restante 15 por ciento, cae en los
otros ocho meses del año. (Benites Cardenas, A., 2002)
Con respecto a la ubicación geológica del área de estudio la misma se ubica en el

litoral ecuatoriano, al Oeste de la cordillera occidental de Los Andes. Pertenece a la
llanura aluvial que forma una zona distal de una serie de abanicos aluviales de
material cuaternario, a los que se les ha designado como el cono de deyección del
rio Chimbo, un tributario del rio Babahoyo, uno de los dos principales afluentes del
Rio Guayas.
74
Ilustración 4-14.-Cuencas Hidrográficas del Ecuador-Cuenca 130 Rio Guayas (INAMHI,

2008)
Las temperaturas máximas y mínimas media anual de Milagro, basadas en la

estación que posee el Instituto Nacional de Meteorología e Hidrología en los
predios de la Compañía Azucarera Valdez, indican que son 29.9 y 21.7 °C.
(Benites Cardenas, A., 2002).
La humedad relativa es del orden de 75% y 87 %, donde el periodo con mayor
humedad es el que va desde el mes de enero a julio. El promedio anual de
precipitación para Milagro, tomado de la estación meteorológica es de 1.361 mm.
Entre Enero y Mayo llueve aprox. el 95% de la media anual. La evaporación anual
es de 1.346 mm. (CELEC, 2012).
75
Estacion Precipitacion Temperatura( Evaporacion( Evapotranspiracion(
meteorologica Latitud Longitud (mm) °C) mm) mm)
Milagro 02°06´56´´ 79°35´57´´ 1526.10 25.8 1331.2 735
San Carlos 02°12´00´´ 79°26´35´´ 1843.50 27.3 1000.6 764
Guayaquil 02°09´30´´ 79°53´02´´ 1242.40 26.7 1703.5 739
Babahoyo 01°49´41´´ 79°33´49´´ 2108.60 26.0 1300.9 784
Tabla 4-8 Valores de Estaciones Meteorológicas cercanas al sitio de Estudio. (INAMHI,
2013)
La Tabla 4-8 presenta las estaciones climatológicas con valores de temperatura,

precipitación, evaporación y evapotranspiración media anual de las cuatro
estaciones circundantes al área de estudio, donde la estación de Milagro se
encuentra sobre la zona de interés. El análisis de los datos mostrados son entre
los años 1981 al 2013, es decir, la serie tiene 33 años.
4.7.1 Precipitación media anual
El tipo de lluvias en el área de estudio es característico de invierno que es la

estación lluviosa, similar a la mayor parte de la región Litoral del país. En invierno
la precipitación es producida por el calentamiento de las aguas del Océano Pacifico
frente a las costas Ecuatorianas, debido a la corriente cálida marina procedente del
Norte. En verano las lluvias cesan, esta estación es seca, por la influencia de la
corriente fría del sur Humboldt. A continuación se colocan los datos de las
estaciones cercanas al área de estudio. Ilustración 4-15.
Precipitación
2500.00
2108.60
2000.00 1843.50
Precipitacion(mm)
1526.10
1500.00 1242.40
1000.00 Precipitación
500.00
0.00
Milagro San Carlos Guayaquil Babahoyo
Estaciones Meteorológicas
Ilustración 4-15 Precipitación Media Anual Estaciones circundantes. (INAMHI, 2013)
4.7.2 Temperatura Media Anual
Las temperaturas en las estaciones cercanas, fluctúan entre 25ºC a 27ºC de valor
promedio de temperatura media anual, según se muestra la Ilustración 4-16.
76
Temperatura
27.5 27.3
27.0
Temperatura (°C) 26.7
26.5
26.0
26.0 25.8
Temperatura
25.5
25.0
Ilustración 4-16 Temperatura Media Anual. (INAMHI, 2013)
4.7.3 Evaporación media anual
A continuación, se presentan los datos de Evaporación de las estaciones

próximas al sitio de estudio. (Ilustración 4-17).
Evaporación
1800 1703.5
1600
1331.2 1300.9
1400
Evaporacion(mm)
1200 1000.6
1000
800
600 Evaporación
400
200
0
Ilustración 4-17 Evaporación Media Anual. (INAMHI, 2013)
4.7.4 Evapotranspiración media anual
Se muestran las estaciones cercanas a la zona de estudio con los valores de

Evapotranspiración. (Ilustración 4-18).
77
Evapotranspiración
770 764 764
760
Evapotranspiración
750
739
740 735 Evapotranspiración
730
720
Ilustración 4-18 Evapotranspiración Anual
El río Milagro posee una longitud de 80 kilómetros aproximadamente, inicia su

recorrido en las cercanías del poblado conocido con el nombre de Barraganetal,
perteneciente al Cantón Naranjito, a una altitud de 110 msnm aproximadamente, al
pasar por la ciudad de Milagro con una altitud aproximada de 12.55 msnm, la
atraviesa de Este a Oeste, y desemboca en el Río Yaguachi
Ilustración 4-19 Cuenca del Rio Milagro (SNI, 2009)
En la Ilustración 4-19 de color marrón al centro de la gráfica, se observa la cuenca

del Río Milagro.
78
Ilustración 4-20 Mapa geomorfológico de la Cuenca del Río Milagro (RemotosCLIRSEN,

2009)
En la Ilustración 4-20, se puede distinguir que el área de estudio en los predios de

la Universidad Estatal de Milagro, corresponde a una unidad ambiental
denominada como Llanura Aluvial, constituida por Depósitos Aluviales de textura
variable (arcillas, limos y arenas no consolidadas), con una morfología de origen
Deposicional es un Dique o Banco Aluvial, cuya morfometría tiene una pendiente
dominante de 0-2 %, que corresponde a un desnivel relativo entre 0-5 metros.
79
Ilustración 4-21 Subcuencas del Rio Guayas (SNI, 2009)
4.8 Geología de la Cuenca del Río Milagro

La cuenca baja del Guayas, cercana a la desembocadura del Océano Pacífico,
tiene una extensión aproximada de 4.200 Km² y está formada por las subcuencas
hidrográficas de los ríos Milagro, Chimbo-Chanchán, Barranco Alto, Culebras,
Estero Verde, Bulubulu, Ruidoso y Churute.
Estratigrafía de la zona Cuenca del Río Milagro
Dos tipos de rocas sedimentarias reposan sobre el basamento oceánico de la

región: Rocas tipo flysch, con edades que van desde el Cretáceo superior hasta el
Eoceno inferior y del tipo molásico, desde el Oligoceno superior hasta el
Cuaternario. (Dugas F., 1987)
Depósitos Aluviales (Holoceno)
Lo integran conglomerados, arenas, limos y arcillas depositadas a lo largo de los

actuales ríos, dispuestos en terrazas emplazadas sobre las formaciones anteriores,
excepto los depósitos estuarinos que son contemporáneos. (Dugas F., 1987)
80
4.8.1 Rasgos estructurales e Historia Tectónica
Regionalmente en el Sur de la costa ecuatoriana, se observan rasgos

estructurales, que nos indican una tectónica de bloques fallados.
La Cuenca del Guayas perteneciente al Mio-Plioceno, rellenada por espesores

importantes del Plio-Cuaternario. Finalmente el Golfo de Guayaquil, que representa
la depresión actual.
Dos rasgos tectónicos importantes se identifican en la costa ecuatoriana: uno,

relacionado con la tectónica de fallamientos sucesivos de bloques inclinados
diferentemente, que a lo largo de la historia estructural de la región, ha creado un
sistema de horsts y grabens; y otro, relacionado con el desplazamiento de grandes
masas, originando cuerpos alóctonos, de controvertida interpretación.
A continuación se describirá la secuencia de eventos que tuvieron lugar durante la

evolución geológica de la costa del Ecuador:
Una fase de compresión, iniciada por un desplazamiento del piso oceánico, en el

Cretáceo medio. A partir del Turoniano superior, en las partes deprimidas de la
dorsal y lejanas a zonas continentales, se acumularon sedimentos pelágicos de
consistencia silícea y carácter turbidítico profundo que constituyeron los Miembros
Cayo y Guayaquil de edad Senoniano y Maestrichtiano, respectivamente; a esta
fase de compresión algunos investigadores atribuyen el metamorfismo de la
corteza oceánica, que produjo los metaesquistos verdes de la Formación Churute,
así como también algunas rocas sedimentarias (Formación Punta de Piedra). Las
intrusiones granodioríticas existentes en el área, han sido interpretadas como
consecuencia de la continentalización del piso oceánico. (Barragan Chang E.,
2002)
Un mecanismo de acreción de corteza oceánica, tuvo lugar al fin del Cretáceo

(Formación Macuchi), así como la formación de una fosa y un talud profundo,
emplazados al Sur de la Cordillera de Chongón-Colonche, la cual se encontraba en
proceso de levantamiento. La depresión fue rellenada en principio por flysch
turbidítico del Grupo Azúcar inferior (Formación Estancia), de edad Paleoceno y
posteriormente por descargas de areniscas y conglomerados en depósitos de
canales (Formación Chanduy), gracias a la pendiente SE-NW que existía desde los
cerros de Amotape en el Norte de Perú hacia la Península de Santa Elena, a través
del Golfo de Guayaquil.
En el Eoceno medio, un flysch turbidítico calcáreo (Formación San Eduardo), se

depositó en los bordes de la Cordillera Chongón-Colonche y del Horst Azúcar-
Playas; siendo sus aportes calcáreos, provenientes de la Cordillera Occidental de
los Andes, donde se reportan arrecifes. (Barragan Chang E., 2002)
Al fin del Eoceno medio un evento tectónico importante, levantó la pendiente que
se mantenía en el sentido SE-NW, movilizando el talud hasta adquirir una
81
inclinación hacia el SW y formar la Cuenca Ancón. Como consecuencia de este

levantamiento, en el Eoceno superior se originaron deslizamientos del Grupo
Azúcar y en algunos sectores de la Cordillera de Chongón-Colonche, del Miembro
Guayaquil Chert y Piñón, hacia la Depresión Ancón. De esta manera, bloques de
diversos tamaños, relacionados con las formaciones pre-existentes, se
desplazaron, usando como matriz un flysch silíceo, con turbiditas distales, que
rellenaban la depresión, constituyendo así la Formación Socorro, denominada Wild
Flysch, de la Península de Santa Elena y de la región litoral, hasta Puerto López al
Norte.
Capas masivas de grauvacas y mantos delgados de areniscas y lutitas, se

acumularon sobre Socorro, constituyendo la formación superior del Grupo Ancón
(Formación Seca). En el norte de la Cordillera de Chongón-Colonche, la formación
San Mateo es el equivalente directo de Socorro. (Barragan Chang E., 2002)
Una laguna estratigráfica regional, separa el Eoceno superior y el Oligoceno

superior; haciéndose evidente por la ausencia de depósitos de estas edades.
Durante el Oligoceno superior se intensificó el levantamiento de la cordillera de

Chongón-Colonche, dando origen al hundimiento de algunas áreas, como la
cuenca de Manabí, Cuenca del Guayas y el Golfo de Guayaquil. Este evento se
relaciona directamente con un movimiento dextral de la falla transcurrente
Guayaquil-Dolores, que con rumbo paralelo a los andes se emplazó
perpendicularmente al eje longitudinal de la Cuenca Progreso.
Al fin del Oligoceno superior y durante el Mioceno inferior, se produjo una

transgresión marina, dando origen a la acumulación de un relleno tipo molásico, de
más de 5.000 m. de espesor (Formación Tosagua, Subibaja y Progreso en la
Cuenca Progreso y sus equivalentes Angostura Onzole y Borbón en la Cuenca de
Manabí).
Al fin del Mioceno superior y principios del Plioceno, un evento tectónico hizo
emerger la cuenca de Manabí y Progreso, mientras que una sedimentación más
litoral siguió en la Cuenca del Guayas y Golfo de Guayaquil, al pie mismo de la
Cordillera Occidental. Por efecto de este nuevo movimiento de bloques, se levantó
definitivamente la Cordillera de Chongón-Colonche; y los escalones del flanco Sur,
se curvaron y deslizaron parcialmente las capas de la Formación Subibaja y del
Miembro Bellavista.
En el estuario de la alta cuenca del Guayas y durante el Plioceno se sedimentó la

Formación Balzar, mientras que al mismo tiempo y en el Golfo de Guayaquil se
depositó la Formación Puná.
82
Ilustración 4-22 Mapa Geológico del Ecuador (Ministerio de Recursos Naturales y

Energeticos, 1982)
Durante el Cuaternario y después de un débil hundimiento de toda la región, que

permitió la acumulación de los Tablazos en el Litoral y la Formación Pichilingue, en
la alta cuenca del Guayas, se produjo el levantamiento general de la región, que
hizo emerger estas formaciones, mientras que la sedimentación continuó en el
Golfo de Guayaquil. Al final el resultado de todo este tectonismo y acumulación de
sedimentos dio a la costa su forma actual, que es como se conoce.(Barragan
Chang E., 2002)
83
Ilustración 4-23 Perfil Geologico Ecuador (Ministerio de Recursos Naturales y Energeticos,

1982)
Ilustración 4-24 Zona Geologica donde se encuentra la zona de estudio (INAMHI, 2008)
84
5 Trabajos de Campo Sub Cuenca del Milagro

Para poder realizar el modelo hidrogeológico fue necesario ejecutar trabajos de
campo.
5.1 Metodología y Materiales.
5.1.1 Equipos utilizados.
Se utilizaron una sonda y dos tipos de sensores con un data loggers de marca
SOLINST, entre otros equipos que a continuación se describen:
1. Una sonda freatimetrica Modelo 102 MINI. Longitud del cable es 25 metros
y dimensiones de sonda 6.4mm x 38mm.
2. Un sensor de Nivel y Temperatura Levelogger Edge-M30. Exactitud del

sensor de nivel +- 0.05%FS.Exactitud de sensor temperatura +- 0.05°C.
3. Un sensor de Presión Barométrica Barologger M 1.5 compensador de

mediciones del Levelogger.
4. Cable de Acero para Barologger de 15 metros de longitud.
5. Cable de Acero para Levelogger de 30 metros de Longitud.
6. Set de Comunicación USB para Levelogger.
7. CD con Software de Instalación.
8. Trazador Testeador de bolsillo Marca La Motte, EC/TDS/SALT, mide la

Conductividad, el TDS (Sólidos Totales disueltos) y la salinidad. Rango de
medición TDS= 0 a 9,999 ppm, Salinidad 0 a 9,999 ppm, Temperatura 0°C a
65°C.
9. Computador portátil que efectuara las funciones de: configuración, captura,

análisis y compensación preliminar in situ de las mediciones de los sensores
Marca Solinst.
10. GPS Marca Garmin.
11. Cinta Métrica 30 metros
85
5.1.2 Obtención de datos de los instrumentales de Campo
Se procede a instalar el software en la pc, posteriormente se conecta y sincroniza

el sensor, configurando las unidades apropiadas, el tiempo en que el sensor realiza
las lecturas, que en este caso es un intervalo de 15 minutos.
El sensor Levelogger se lo coloca a una profundidad de 22 metros y el Barologger

se ubica a una profundidad de 6 metros de profundidad, es decir, algo distante del
nivel del agua.
Es importante sincronizar ambos sensores con las mismas unidades y los mismos
intervalos de tiempo, para que posteriormente haya una precisión en la efectividad
de compensación entre los valores que registran tanto el Levelogger y Barologger,
pues como se ha descrito anteriormente, el Barologger registra los datos de
presión atmosférica, mientras que el Levelogger registra la presión total, es decir,
presión barométrica + presión de agua. En otras palabras, la lectura del Levelogger
menos la lectura del Barologger se obtiene la altura de columna de agua, que se
ejecuta con la mencionada compensación que el mismo software de los equipos
efectúa.
Una vez que se ha determinado el tiempo de registro en cada pozo, se retira el

sensor y se descarga la información a la PC para su posterior análisis. Vale indicar
que se debe considerar la información de precipitación del intervalo de tiempo que
se está registrando.
5.1.3 Metodología
Para esto se definió un terreno de flujo tridimensional con una longitud de 1500 m.,
ancho de 1200 m. y una profundidad de 44 m.
En la zona que se definió para el presente estudio existen 2 capas, en las cuales
se representa la zona no saturada y la saturada. La capa que corresponde a la
zona no saturada tiene un espesor aproximado de 14 m. La capa que corresponde
a la zona saturada, que es de interés de este estudio, tiene un espesor de 30 m.,
que en total nos da una profundidad total de 44 m.
Una vez que se precisan las dimensiones y algunas características del nivel
freático existentes en la Ciudadela Universitaria UNEMI, se procede a detallar los
procedimientos que se ejecutaron en los pozos existentes.
Se realizó el monitoreo de los pozos mencionados y con fecha 31 de enero del
2013, se instalaron los sensores respectivos en el pozo ubicado en la plazoleta de
la referida institución de educación superior, situado en las siguientes coordenadas:
S 02°09’01.0”; W 79°36’07.2”, posteriormente con fecha 28 de Febrero del 2013 se
retiraron del mencionado pozo los instrumentos de medición .Este pozo
actualmente está en producción y se extrae agua mediante una Bomba sumergible
de 2200 Watios de Potencia.
86
Luego en el pozo situado atrás del Bloque K, ubicado en las siguientes

coordenadas: S 02°09’06.3”; W 79°36’11.0”, se instaló los equipos el día 15 de
Abril del 2013 y fueron retirados del pozo el 18 de mayo del 2013. Previamente
estos equipos se calibran, de acuerdo al manual del fabricante.
Como paso previo a la instalación de los sensores, se realiza con el freatimetro, la

medición del nivel estático del pozo. Posteriormente, se midieron los indicadores de
calidad de agua con el Conductimetro, Marca La Motte, EC/TDS/SALT TRACER,
que son la Conductividad, el TDS (Sólidos Totales disueltos) y la salinidad.
Fecha Registro
Nombre del Tipo de Fecha
Pozo Latitud Longitud Sensor Inicio Fecha Fin
Pozo 1(Plazoleta) 655465S 9762252E Freatígrafo 31-ene-13 28-feb-13
Pozo 2(Bloque K) 655346S 9762106E Freatígrafo 15-abr-13 18-may-13

Tabla 5-1 Pozos Registrados Universidad Estatal de Milagro
5.2 Modelación Numérica
El esquema de resolución numérico empleado por el modelo, corresponde a un

esquema de diferencias finitas, el cual entrega soluciones aceptables a los
problemas modelados, a través de una metodología numéricamente sencilla. La
representación gráfica de la modelación facilita el ingreso de los datos, la
visualización de resultados, etc. Las dos acciones externas que se debe considerar
en un modelo de simulación de aguas subterráneas son “la recarga y la
extracción”.
La recarga se obtiene mediante el balance hídrico. Se debe evaluar, no sólo la

recarga natural, sino además el volumen recargado por riego o por mecanismos
artificiales. No se incluye el efecto de recargas en cursos de agua, ya que estas se
indican en las condiciones de borde. (Espinoza C., 2012)
La extracción incluye aquella realizada mediante pozos de bombeo, así como

también aquella que ocurre por otros mecanismos como la evapotranspiración.
87
Ilustración 5-1 Mapa de Izopiezas Cuenca del Guayas ( Instituto Nacional de Meteorología
e Hidrologìa, 2013)
MUESTRA
DESCRIPCION ESTRATIGRAFIA
# PROF.
1 1.50
LIMO DE ALTA PLASTICIDAD COLOR GRIS
2 2.50
OSCURO DE CONSISTENCIA MUY COMPACTA
3 3.00
LIMO DE BAJA PLASTICIDAD GRIS VERDOSO 4 4.50
LIMO GRIS OSCURO DE CONSISTENCIA MUY
5 5.00
BLANDA
ARCILLA GRIS OSCURO DE CONSISTENCIA
6 6.00
MUY BLANDA
88
TURBA 7 7.00
LIMO ARCILLOSO CON INTERCALACIONES DE
8 8.00
TURBA
ARCILLA DE ALTA PLASTICIDAD GRIS
9 9.50
VERDOSA DE CONSISTENCIA BLANDA
LIMO DE BAJA PLASTICIDAD COLOR GRIS
10 10.00
VERDOSO DE CONSISTENCIA MUY BLANDA
LIMO DE ALTA PLASTICIDAD COLOR GRIS
11 11.00
VERDOSA CON ALGO DE YESO
LIMO DE BAJA PLASTICIDAD DE COLOR CAFÉ
CLARO CON ARCILLA GRIS VERDOSA Y ALGO 12 12.00
DE YESO
LIMO DE ALTA PLASTICIDAD CAFÉ CLARO 13 12.50
CON INTERCALACIONES GRIS VERDOSO DE
14 14.00
CONSISTENCIA MUY COMPACTA
LIMO DE BAJA PLASTICIDAD CAFÉ CLARO CON
LIMO GRIS CLARO DE CONSISTENCIA MUY 15 15.00
COMPACTA
LIMO DE ALTA PLASTICIDAD COLOR CAFÉ
16 16.00
CLARO DE CONSISTENCIA MUY COMPACTA
LIMO DE BAJA PLASTICIDAD COLOR CAFÉ
17 17.00
OSCURO DE CONSISTENCIA DURA
LIMO DE BAJA PLASTICIDAD COLOR GRIS
18 18.50
CLARO DE CONSISTENCIA DURA
19 19.00
20 20.00
21 21.00
22 22.00
23 23.00
FIN DE LA PERFORACION
Tabla 5-2 Perfil Estratigráfico de una perforación cercana a los pozos referidos. (Moreno L.,
2013)
En la decada de 1980, el INAMHI elabora el Mapa Hidrogeológico Nacional del

Ecuador, que abarca una compilación de cerca 5.000 puntos de agua
inventariados, que sirvió como base fundamental para formular y sobretodo definir
las Unidades Hidrogeologicas del Ecuador, cuya conceptualización es actuar como
medios acuìferos continuos dotados de cierta homogeneidad, que se constituyen
en unidades naturales de planeamiento y gestion del recurso. ( Instituto
Ecuatoriano de Recursos Hidraulicos-Inerhi, 1986)
El objetivo principal, fue definir para cada unidad hidrogeològica, sus

caracterìsticas bàsicas, que se cuantifique en la medida de lo posible sus
elementos principales, como son geometría, heterogeneidad geológica,
productividad acuìfera, piezometría y profundidad del agua, calidad quimica,
recarga media interanual.
89
Con referencia a la unidad hidrogeologica correspondiente a la zona de estudio,

que es la Unidad Hidrogeologica del Guayas, tiene una superficie aproximada de
28.000 KM2, se encuentra cubierta por sedimentos del Cuaternario; en la parte
oriental se destacan rocas cretacicas, hacia el Oeste afloran rocas de origen
marino del Terciario.Los caudales reportados de los pozos perforados varían de
0,0002 m³/s a 0,096 m³/s, la profundidad del nivel freatico oscila entre 3m a 30m.La
calidad del agua en general es muy buena.Como caracterizaciòn hidrogeològica de
la cuenca del rìo Guayas, donde se encuentra la zona de estudio, se puede indicar
que el acuìfero principal esta relacionado con los sedimentos
cuaternarios(aluviales). A manera de referencia en los cantones Pedro Carbo,
Santa Lucia y subcuencas de los rios Macul, Vinces y Baba, los pozos pueden
alcanzar caudales de 0,015 m³/s, estos sitios se encuentran a unos 68 km. del sitio
de estudio.Ver Ilustración 5-1
La subcuenca del rio Milagro, en el sector de Chobo; se estima que la producciòn

diaria es de 60.000 m³, además en el sector de Naranjito,cercano al Cantòn
Milagro,especìficamente en el Ingenio Azucarero San Carlos, existen alrededor de
64 pozos con fines agroindustriales que proveen un caudal de 11.4 m³/s.
Ilustración 5-2 Mapa Geomorfologico de la zona de estudio.
(SISTEMA NACIONAL DE INFORMACION-SNI, 2009)
90
Estación Milagro Enero-Junio (Inamhi)

140.0
120.0
Precipitación(mm)
100.0 Enero
80.0 Febrero
60.0 Marzo
40.0 Abril
Mayo
20.0
Junio
0.0
0 10 20 30 40
Días
Ilustración 5-3 Datos de Precipitación Diaria (INAMHI, 2013)
En la Ilustración 5-3, se observa los meses donde se presentan los picos de

precipitación, que son Febrero y Marzo.
Los datos de la conductividad de la zona de interés, se obtiene de un estudio

realizado a una finca bananera en la parroquia Mariscal Sucre, 15 Km a Nor-Este
de Milagro y cuyo valor es: 10,1m/d.
Utilizando los datos en el subcapítulo 5.1.3 y usando las 3D, para representar el
cubo a modelar, se obtiene un volumen de 79’200.000 de m³ de terreno de estudio.
Precisando un poco más los escenarios de simulación, descritos en el Subcapítulo

1.2, estos son:
• Escenario Actual.- Que es el estado actual del sistema a analizar con dos
pozos con bombas sumergibles y un pozo intermedio de observación.
• Escenario Futuro uno.-Que es el estado proyectado con dos pozos de

bombeo más al escenario actual.
• Escenario Futuro dos.- Escenario proyectado con tres pozos de bombeo

más al escenario actual.
• Escenario Futuro tres.- Se proyecta con tres pozos de bombeo adicionales al

estado actual más la contaminación de un tanque séptico.
91
Ilustración 5-4 Líneas equipotenciales escenario actual.
En Ilustración 5-4 , se observa las líneas equipotenciales del área de estudio con
una conductividad de 10 m/d y recarga de 743 mm/a con el escenario actual.
Con ajuste que da una raíz cuadrada media de 0,178 y error estándar de
estimación de 0,038; se observa la mejor aproximación.
92
Ilustración 5-5 Balance de masas, modelo de masas en escenario actual.
Ingresos por 1’145.109 m³ y 1’145.238m³ en salidas, interpretando esto significa

que lo que ingresa como recarga es menor a lo que por los sistemas de bombeo de
los pozos instalados tiene salida. Se obtiene así -129 m3, debido a lo antes
expuesto.
93
Ilustración 5-6 Líneas equipotenciales del escenario futuro uno.
En Ilustración 5-6 , se observa las líneas equipotenciales del área de estudio con
una conductividad de 10 m/d y recarga de 743 mm/a.
En este escenario futuro uno, se incrementa dos pozos de bombeo, uno llamado
polideportivo y el otro Bloque J, esta calibración arrojó una raíz cuadrada media de
0,818 y error estándar de estimación de 0,117.
94
Ilustración 5-7 Balance de Masas escenario futuro uno
Ingresos por 1’321.491 m³ y 1’321.559m³ en salidas, interpretando esto significa

que lo que está siendo bombeado por los pozos es mayor que el valor de recarga
del área de estudio. De la suma algebraica nos da -68 m³.
Ilustración 5-8 Líneas equipotenciales escenario futuro dos
Se observa en la Ilustración 5-8 , se observa las líneas equipotenciales del área de

estudio. En este escenario futuro dos, se incrementa tres pozos de bombeo, el
situado en el polideportivo, el otro en el Bloque J y uno llamado estándar futuro,
esta calibración arrojó una raíz cuadrada media de 1,76 y error estándar de
estimación de 0,22.
95
Ilustración 5-9 Balance de Masa Escenario Futuro dos
Ingresos por 1’563.673 m³ y 1’563.737m³ en salidas. Se obtiene -64 m³ que arroja

el balance de masas.
Ilustración 5-10 Líneas equipotenciales y trayectoria de partículas de escenario futuro tres
El escenario tres mostrado en la Ilustración 5-10, tiene cinco pozos de bombeo y

además se agregó unas partículas de residuos cloacales, producto de un pozo
séptico que derrama los líquidos al acuífero libre en las cercanías del pozo del
polideportivo.
Tomando en cuenta que el suelo es limo-arenoso, es un volumen representativo

respecto del tamaño del grano de suelo.
96
Ilustración 5-11 Frecuencia Temporal Niveles Calculados
Ilustración 5-12 Niveles Observados
97
En Ilustración 5-11 e Ilustración 5-12 se puede observar los valores observados y

los calculados con el mejor ajuste.
Validación
Primero se debe establecer el rango a validarse, que partirá del valor que mejor
ajuste mostró, entonces posteriormente se ingresara los valores restantes del tal
manera que, exista una gráfica que nos demuestre cual es la mejor aproximación.
Para esto se elabora una tabla con los valores de Recarga y conductividad
ajustados con los respectivos valores obtenidos por raíz cuadrada media.
Validación recarga
1.8
1.6
1.4
Raiz Cuadrada Media
1.2
1
0.8 Validación Parametro
0.6 Recarga
0.4
0.2
0
0 200 400 600 800 1000 1200
Valores de recarga
Ilustración 5-13 Validación Recarga
Recarga que muestra el valor calibrado y que de acuerdo a la gráfica se valida.
A partir de los resultados del modelo calibrado en estado estacionario donde

resultó una conductividad hidráulica de 10 m/d y las recargas mostradas en las
figuras respectivas se efectuó el análisis de sensibilidad de los parámetros
manejados para observar la sensibilidad del modelo frente a la variación de los
mismos. En una primera instancia se dejó constante la conductividad hidráulica y
se procedió a variar la recarga en un rango de 400 a 1100 mm/año. Posteriormente
sobre las recargas calibradas, se varió la conductividad en el rango de 4 a 20 m/d.
Así se obtuvieron varias nuevas simulaciones para los procedimientos
mencionados, donde se evaluó la influencia de estas transiciones en los medidas
sobre los niveles afectados que se ajustó en el pozo de observación. De esta
manera, el estudio de sensibilidad es índice del parámetro que altera al nivel
supuesto.
98
En las figuras referidas a las validaciones se observan los resultados de la

magnitud de la raíz cuadrada media con cada cambio en los valores de la recarga y
conductividad hidráulica.
Validación conductividad
0.5
0.45
0.4
Raiz cuadrada media
0.35
0.3
0.25
0.2 Validación Parametro
Conductividad
0.15
0.1
0.05
0
0 5 10 15 20
Valores de conductividad
Ilustración 5-14 Validación de Conductividad.
5.3 Calibración del modelo y de parámetros
La calibración radica en obtener un conjunto de valores de los parámetros del

modelo donde se incluye las propiedades de acuífero, los términos fuente-
sumideros, las condiciones iniciales y condiciones de frontera que permiten que el
modelo represente los valores de observaciones (datos cualitativos) dentro de un
error aceptable. Existen dos procedimientos generales para la calibración: el ajuste
manual por ensayo y error, y la estimación automatizada de parámetros (Win
PEST).
El ajuste manual suele ser las más utilizada hoy en día para los modelos, sin
embargo suele ser muy subjetivo y conducente a resultados difícil de evaluar
(Carrera J.y Neuman, 1986). Por otra parte, el uso de la calibración automatizada
con Win PEST consiste en asignar un límite de valores de calibración a partir de los
datos alimentados al modelo. Así mismo el procedimiento automatizado calcula
directamente la sensibilidad del modelo hacia sus parámetros y la correlación de
estos mismos.
Esta información permite determinar si los parámetros del modelo y las

predicciones fueron calculados de manera confiable con los datos disponibles, y
que datos adicionales son requeridos para mejorar el modelo (Poeter, 1997)
Resumiendo en sí, la calibración es el ajuste de los parámetros o un grupo de
99
propiedades hidráulicas o condiciones de frontera ajustadas, que concuerden con

los valores medidos en campo dentro de un rango de error preestablecido. WinPest
se emplea como un módulo de ModFlow que ayuda al proceso de calibración por
medio de algoritmos matemáticos complejos, PEST puede ayudar a determinar si
el modelo conceptual es razonable para encontrar el peor y el mejor caso de
predicción bajo condiciones calibradas. Sin embargo PEST no garantiza que los
resultados sean 100% razonables, si esto ocurre, podría ser un fuerte indicador de
que nuestro modelo conceptual es incorrecto. (Carrera J.y Neuman, 1986)
El modelo implementado en la zona de estudio de Ecuador, se utilizó el criterio

cuantitativo, basado en comparar los valores medidos y los valores simulados.
Es por esto que, se tomó en cuenta el cálculo de residuales hasta encontrar un

valor razonable que permita el ajuste apropiado.
La conductividad hidráulica fue el valor calibrado, con los valores obtenidos de los
pozos ya mencionados que existen en los predios de la Universidad Estatal de
Milagro. El valor de la recarga, parámetro importante en una modelación
hidrogeológica, se obtuvo recopilando los datos de la Estación Meteorológica
Inamhi de la ciudad de Milagro, tales como precipitación y evapotranspiración(ET).
Gráfico Calibración
20
18
16
Valor de Conductividad(m/d)
14
12
10
Gráfico Calibración
8
0
0 0.1 0.2 0.3 0.4 0.5
Raíz cuadrada del Residual
Ilustración 5-15 Grafico Calibración Conductividad Hidráulica
100
6 Trabajos de campo Sub Cuenca del Suquía

6.1 Metodologìa y Materiales
Los estudios en la cuenca del sur de la ciudad de Córdoba que corresponden al

Río Suquía se detalla a continuación.
Para ello se utilizó de base, el pozo de monitoreo que se ejecutó dentro del
convenio que mantiene Aguas Cordobesas S.A., con la empresa desarrollista de la
zona, dicho convenio de monitoreo se realizó porque el conducto que sale desde la
planta potabilizadora Los Molinos atraviesa el terreno y era de interés poder
monitorear las afectaciones que el conducto podría sufrir. Sobre la traza del
conducto no hay edificaciones, solo se emplazan un espacio recreativo (cancha de
golf) y dos lagunas de pequeñas dimensiones. El objetivo de este convenio es
realizar un monitoreo que permita identificar desperfectos de la cañería ante un
eventual hecho.
Ilustración 6-1 Ubicación del conducto de salida de Aguas Cordobesas en el predio
El predio del loteo, donde se encuentran ubicados los instrumentales, presentan un

total de 355 lotes con superficies individuales del orden de los 1.500 m², con
amplios espacios verdes y una cancha de golf de 9 hoyos en su interior.
101
Ilustración 6-2 Laguna al sur del Club House que se encuentra ubicada por encima del
conducto de Aguas Cordobesas en el predio
6.1.1 Equipos utilizados

Dentro de este convenio se realizan mediciones de humedad en el suelo,
precipitación y monitoreos de freática, para ello se encuentran instalados en el
predio diverso instrumental. El equipo que se instaló en el año 2011 consta de:
1) Tres Sensores ECH2O de Humedad, Temperatura y Conductividad

Eléctrica (5-TE) del Suelo. Rango 0-100% humedad volumétrica. Longitud
parte sensor 5 cm.
2) Cuatro Sensores ECH2O de Humedad de Suelo (EC-5). Rango 0-100%
humedad volumétrica. Longitud parte sensor 5 cm.
3) Tres Sensores ECH2O de Humedad de Suelo (10 HS). Rango 0-100%
humedad volumétrica. Longitud parte sensor 10 cm.
4) Pluviómetro marca Decagon modelo ECRN-50, resolución 1 mm, tipo de
sensor de cuchara simple, salida plug stereo, material de carcasa ABS,
temperatura óptima de trabajo 0° a 60°C, tamaño 5 cm x 10 cm, compatible
con datalogger Em 50, y otros.
5) Sensor de temperatura del aire marca Decagon, modelo ECT, precisión
entre 5° y 40°C de temperatura de trabajo: ±0,5°C, fuera de este rango la
precisión es de ±1,0°C
Equipo instalado en Diciembre de 2014:
6) Un freatigrafo modelo LT F300/M100 y barógrafo modelo LT F5/M15 ambos

de la marca Solinst, para monitoreo de niveles de aguas subterránea. El
rango de precisión de la lectura para una profundidad entre 20- 50 metros de
la presión es de +/- 2,0 cm de columna de agua (cmH2O) con una resolución
de 0.4 cmH2O.
102
6.2 Afectación del recurso en la zona

La urbanización dispone una ocupación de alrededor del 50 % (200 viviendas) y
tiene instalada una planta de tratamiento de líquidos cloacales, de acuerdo a las
lecturas recogidas por el caudalímetro electromagnético que se encuentra a la
salida de la cámara de desbaste de rejas, se consumen diariamente un promedio
de 1,20 m3/d por lote, lo que arroja en la actualidad un consumo medio diario de la
urbanización de 240 m3/d. El caudalímetro mencionado, dispone además un
sistema para la medición y almacenamiento en una memoria interna de la
conductividad del afluente a la planta de tratamiento, variable esta de fundamental
importancia al momento de cuantificar la aptitud del efluente. De acuerdo a los
registros el valor referenciado se encuentra en el orden de entre los 400 a 600
s / cm , rango este correspondiente con los valores admisibles para agua potable
(400 a 800 s / cm) , presentando en algún momento del mes picos eventuales,
pero de muy corta duración, como resultado del vertido de ciertos productos al
interior del sistema cloacal por parte de los residentes de la urbanización. En la
Ilustración 6-3, se puede advertir el registro correspondiente a dos meses de
registro utilizado durante la etapa de proyecto del sistema de desinfección UV,
donde queda evidenciado los valores de conductividad señalados anteriormente.
Ilustración 6-3 Registro Salinidad vs tiempo de caudales de ingreso a la planta de

tratamiento de la urbanización
103
Parámetro Unidad Valor
pH 8,90
Conductividad
eléctrica
s / cm 2220
Salinidad ppm 1020
TDS ppm 1400
Tabla 6-1 Parámetros de calidad del agua de la perforación
Con el fin de disponer de un panorama claro, se presenta a continuación una

planimetría general de la urbanización con la indicación de las unidades de
monitoreo existentes.
Ilustración 6-4 Ubicación de los sistemas de monitoreo
6.3 Factores a monitorear

Los factores que se monitorean son humedad del suelo, nivel freático, precipitación
y presión atmosférica. La litología que predomina en la zona son depósitos de
limos, limos arcillosos, limo arenosos intercalados con arenas limosas de la
Formación General Paz, y en forma subyacente se encuentra la Formación Río
Primero compuesto por arenas, gravas y pelitas.
104
El perfil geológico en el terreno, obtenido de la perforación de estudio es el

siguiente:
Profundidad (m.) Litología

0,00 -4,00 Limo arenoso color castaño, macroporoso
4,00- 8,00 Limo arcilloso, ligeramente plástico
8,00 -12, 50 Limo arenoso color castaño, en parte cementado
12,50 -19,00 Limo arenoso, fracción arena abundante, húmedo.
19,00 -22.00 Limo arenoso intercalado con arena limosa, muy húmedo
Tabla 6-2 Perfil Geológico Zona Loteo Cañuelas
6.3.1 Humedad
Dentro de las tareas de monitoreo que se realizan una de las variables es la

humedad del suelo, lo que permitiría realizar una detección temprana de las
eventuales pérdidas de agua tanto de la laguna como del caño.
La dirección de escurrimiento superficial en la zona de estudio es O – E. El terreno

se presenta plano sin grandes ondulaciones.
La posición del cuerpo de agua freático en el mismo se ubica a la profundidad de

30m., según datos conseguidos en la zona y verificado en censo de pozos (P1).
Por encima de esta profundidad el suelo se encuentra con mucha humedad, lo que
puede formar capas locales denominadas colgadas o postizas.
6.3.2 Nivel freático
Los movimientos del nivel freático (tanto ascensos como descensos) son más
funciones regionales que locales. Para identificar si el cambio en la humedad del
suelo se ha producido por algún cambio estacional y regional o es producto de un
fenómeno local, es necesario monitorear los niveles freáticos de la zona.
Los instrumentos para realizar este tipo de medición son los piezómetros o los
sensores. Para la obtención en continuo del nivel piezométrico, existen diferentes
propuestas tecnológicas en función del sensor utilizado y del principio físico en el
que se basan.
A los fines de determinar la dirección del flujo subterráneo en la zona de estudio se

efectuó un censo de pozos en el área de influencia, donde se realizaron las
mediciones conducentes a determinar la profundidad del cuerpo de agua freático.
P1 perforación Cañuelas 30 m. ; P2 perforación Nº1 Country Jockey Club: 24m. ;
P3 perforación Inaudi: 29.00 m.; y P4 pozo molino Cno. San Carlos Km 6.5: 37 m.
105
De acuerdo a estos datos el movimiento del cuerpo de agua freático en la zona

tiene circulación NO – SE, coincidente con el escurrimiento subterráneo regional,
diferente de la dirección de escurrimiento superficial.
Las cuatro unidades de monitoreo han sido evaluadas, a continuación se detalla la

información que se obtuvo del instrumental instalado.
El período monitoreado del que se tienen registros son Enero, Abril 2013, Enero
2014, Enero 2015.
Enero 2013
Unidad 1
Unidad 1: Humedad y Lámina de Lluvia

0.320 30
0.300
0.280
0.260 25
0.240
Lámina de lluvia (mm)

0.220 20
0.200
Humedad
0.180
0.160 15
0.140
0.120
0.100 10
0.080
0.060 5
0.040
0.020
0.000 0
27-dic 1-ene 6-ene 11-ene 16-ene 21-ene 26-ene 31-ene 5-feb
Tiempo (días)
Humedad 10-HS Humedad EC-5b adicional Humedad 5-TE Lámina Lluvia
Ilustración 6-5 Medición de los sensores EC-5b y 5-TE a 0,30 m, 10-HS a 1m. 01-01-13
hasta 31-01-13
El perfil de humedad que se muestra en la Ilustración 6-5 pertenece a las

mediciones tomadas en la unidad 1, se observa que el sensor 10-HS (Ilustración
6-5) registró un promedio de humedad del 27,03% con desviación de 2,20%.
Presentando valores bastantes estables lo que descarta una posible pérdida de la
conducción de Aguas Cordobesas en este sector de la conducción.
Los sensores EC-5b y el 5-TE (Ilustración 6-5) ubicados a 0.3 m de la superficie del
suelo muestran ligeras fluctuaciones en la variación del contenido de humedad del
suelo. En lo referente a los valores propiamente dichos, el contenido de humedad,
presentó promedios de 17,40% y 14,30% y desviaciones de 3,19% y 0,56% de
humedad durante el periodo monitoreado. Alcanzado el sensor EC-5b un pico
máximo de humedad de 19,30% mientras que el sensor 5-TE alcanza un valor de
máximo de 15,00%.
106
Unidad 2
Esta unidad durante el mes de Enero de 2013 pese a los ajustes técnicos
efectuados, siguió presentando inconvenientes técnicos que no permitieron la
recolección normal de lecturas de los sensores ubicados en este sector del
emprendimiento.
Unidad 3
Los sensores 10-HS y EC-5 ubicados a 1 m y 1,80 m, muestran muy ligeras

variaciones de humedad, presentando el sensor 10-HS, una humedad promedio de
13,90% con desviación de 0,30%, como se puede observar este sensor
experimentó un ligero pico en su humedad luego de la precipitación del día 5 de
enero retornando luego a valores normales.
Para el sensor EC-5 presenta una humedad promedio de 12,50% y una desviación
de 0,10%. Estos valores bajos se deben a que la zona de la unidad 3 presenta un
suelo antrópico, muy compactado y con poca vegetación, manteniendo humedades
casi residuales donde las variaciones de humedad son muy escasas.

0.320 30
0.300
0.280
0.260 25
0.240

0.220 20
0.200
Humedad
0.180
0.160 15
0.140
0.120
0.100 10
0.080
0.060 5
0.040
0.020
0.000 0
Tiempo (días)
Humedad 10-HS Humedad EC-5 Lámina Lluvia
Ilustración 6-6 . Medición de los sensores, 10HS a 1m y EC-5 1.80 m. 01-01-13 hasta 31-
01-13
107
Temperatura Máxima y Mínima diaria

45
40
35
Temperatura en º Celsius
Min
30
25
Max
20
15
10
5
0
1-ene
6-ene
21-ene
11-ene
16-ene
26-ene
31-ene
Tiempo (días)
Ilustración 6-7 Medición del sensor de Temperatura ECT 01-01-13 hasta 31-01-13
En la Ilustración 6-7 se muestran los valores de temperatura diarios máximos y

mínimos durante el 01 de Enero y 31 de enero del 2013, con una media máxima
mensual de 32,83 ºC y media mínima mensual de 18,36 ºC; el día 24 de enero
ocurre la temperatura más elevada del mes de 40,66 °C y el 02 de enero se
presenta la mínima del mes de 11,16 ºC.
Unidad 4
Nivel freático
387.92 30
387.90 25
Altura del nivel freático (msnm)
387.88
20
Lámina lluvia (mm)
387.86
15
387.84
10
387.82
5
387.80
387.78 0
387.76 -5
Nivel Freático Lámina Lluvia
Ilustración 6-8 Variación del nivel freático en la unidad 4. 01-12-13 hasta 31-01-13
108
En la Ilustración 6-8 se registra las variaciones del nivel freático, al respecto se

debe mencionar que la unidad registró valores desde el 01 de Enero hasta el 31 de
Enero de 2013, alcanzando un valor máximo de 387,90 m y un valor mínimo de
387,77 m, la media calculada es de 387,82 m.
Abril 2013
Unidad 1

0.320 25
0.300
0.280
0.260 20
0.240

0.220
0.200 15
Humedad
0.180
0.160
0.140
0.120 10
0.100
0.080
0.060 5
0.040
0.020
0.000 0
27-mar 1-abr 6-abr 11-abr 16-abr 21-abr 26-abr 1-may 6-may
Tiempo (días)
Humedad 10-HS Humedad EC-5 Humedad 5-TE Lámina Lluvia
Ilustración 6-9 Medición de los sensores EC-5b y 5-TE a 0,30 m, 10-HS a 1m. 01-04-13
hasta 30-04-13

conducción de Aguas Cordobesas en este sector.
Los sensores EC-5 y el 5-TE (Ilustración 6-9) ubicados a 2,70 y 0.3 m

respectivamente de la superficie del suelo muestran ligeras fluctuaciones en la
variación del contenido de humedad del suelo. En lo referente a los valores
propiamente dichos, el contenido de humedad, presentó valores medios de 22,70%
y 16,80% de humedad durante el periodo monitoreado. Alcanzado el sensor EC-5
un pico máximo de humedad de 23,90% mientras que el sensor 5-TE alcanza un
valor de máximo de 17,70%.
109
Unidad 2
Debido a los problemas que presentaron los sensores de esta estación, se ha

procedido a extraerlos y comprar nuevas unidades similares a las instaladas,
encontrándose en la actualidad en la espera del arribo de dichos sensores para su
instalación.
Unidad 3
Los sensores 5 TE y EC-5 ubicados a 030 m y 2,00 m, muestran muy ligeras

variaciones de humedad, presentando el sensor 5 TE, una humedad promedio de
14,60% con desviación de 0,48%.
Para el sensor EC 5 presenta una humedad promedio de 12,40% y una desviación

de 0,09%. Estos valores bajos se deben a que la zona de la unidad 3 presenta un
suelo antrópico, muy compactado y con poca vegetación, manteniendo humedades
casi residuales donde las variaciones de humedad son muy escasas.

0.180 25
0.160
0.140 20

0.120
15
Humedad
0.100
0.080
10
0.060
0.040 5
0.020
0.000 0
Tiempo (días)
Humedad EC-5 Humedad 5-TE Lámina Lluvia
Ilustración 6-10 Medición de los sensores, 10HS a 1m y EC-5 1.80 m. 01-04-13 hasta 30-
04-13
110
Temperatura Máxima y Mínima diaria del

aire y Temperatura del suelo
40 40
Temperatura del suelo en º Celsius

Temperatura del aire en º Celsius
35 35
30 30
25 25
20 20
15 15
10 10
5 5
0 0
11-abr
16-abr
21-abr
26-abr
1-abr
6-abr
Tiempo (días)

mínimos del aire durante el 01 de Abril y 30 de Abril del 2013, con una media
máxima mensual de 28,13 ºC y media mínima mensual de 12,80 ºC; el día 09 de
Abril ocurre la temperatura más elevada del mes de 34,00 °C y el 13 de Abril se
presenta la mínima del mes de 4,18 ºC.
En la Ilustración 6-11 se presenta conjuntamente la temperatura del suelo, medida

a una profundidad de 30 cm, observándose fluctuaciones normales con la variación
de la temperatura del aire debido a lo superficial de la lectura, con un promedio de
22,36Cº.
111
Unidad 4
Nivel freático
388.17 40
388.16 35
388.15 30
388.14 25
Lámina lluvia (mm)

388.13 20
388.12 15
388.11 10
388.10 5
388.09 0
388.08 -5
Ilustración 6-12 . Variación del nivel freático en la unidad 01. 04-13 hasta 30-04-13

debe mencionar que la unidad registró valores desde el 01 de Abril hasta el 30 de
Abril de 2013, alcanzando un valor máximo de 388,16 m y un valor mínimo de
112
Enero 2014
Unidad 1

0.320 0.09
0.300
0.280 0.08
0.260
0.240 0.07

0.220 0.06
0.200
Humedad
0.180 0.05
0.160
0.140 0.04
0.120
0.100 0.03
0.080 0.02
0.060
0.040 0.01
0.020
0.000 0
Tiempo (días)
Ilustración 6-13 Medición de los sensores EC-5 y 5-TE a 0,30 m, 10-HS a 1m. 01-01-14
hasta 31-01-14

Presentando algunas ligeras fluctuaciones pero estabilizándose en valores
normales a los presentados en los meses anteriores lo que descarta una posible
pérdida de la conducción de Aguas Cordobesas en este sector.
Los sensores EC-5 y el 5-TE (Ilustración 6-13) ubicados a 2,70 y 0,30 m

respectivamente de la superficie del suelo muestran ligeras fluctuaciones en la
variación del contenido de humedad del suelo. En lo referente a los valores
propiamente dichos, el contenido de humedad, presentó valores medios de 20,10%
y 13,40% de humedad durante el periodo monitoreado. Alcanzado el sensor EC-5
un pico máximo de humedad de 26,40% mientras que el sensor 5-TE alcanza un
valor de máximo de 14,40%.
113
Unidad 2

0.320 0.09
0.300
0.280 0.08

0.260
0.240 0.07
0.220 0.06
Humedad
0.200
0.180 0.05
0.160
0.140 0.04
0.120
0.100 0.03
0.080 0.02
0.060
0.040 0.01
0.020
0.000 0
Tiempo (días)
Humedad 10-HS Humedad 5TE Lámina Lluvia
Ilustración 6-14 Medición de los sensores, 10HS a 1m y 5TE 0.30m. 01-01-14 hasta 31-
01-14
Debido a los problemas que presentaron los sensores de esta estación, se ha

procedido a extraerlos y se ha realizado la compra de nuevas unidades y se la
instaló y calibró, habiendo comenzado a tomar lecturas el día 15 de Enero de 2014.
conducción de Aguas Cordobesas en este sector. El sensor 5-TE (Ilustración 6-15)
muestra ligeras fluctuaciones en la variación del contenido de humedad del suelo.
En lo referente a los valores propiamente dichos, el contenido de humedad,
presentó un valor medio de 4,80% durante el periodo monitoreado. Alcanzado el
sensor 5 TE un pico máximo de humedad de 5,00%.
114
Unidad 3

0.260 0.09
0.240 0.08
0.220 0.07
0.200

0.180 0.06
0.160
Humedad
0.140 0.05
0.120 0.04
0.100 0.03
0.080
0.060 0.02
0.040 0.01
0.020
0.000 0
Tiempo (días)
Ilustración 6-15 Medición de los sensores, 10HS a 1m y EC-5 1.80 m. 01-01-14 hasta 31-
01-14
Los sensores 5 TE y EC-5 ubicados a 0,30 m y 2,00 m, muestran muy ligeras

variaciones de humedad, presentando el sensor 5 TE, una humedad promedio de
15,50% con desviación de 0,24%.
Para el sensor EC 5 presenta una humedad promedio de 11,90% y una

desviación de 0,39%. Estos valores bajos se deben a que la zona de la unidad 3
presenta un suelo antrópico, muy compactado y con poca vegetación, manteniendo
humedades casi residuales donde las variaciones de humedad son muy escasas.
115
Temperatura Máxima y Mínima diaria del

aire y Temperatura del suelo
45 40

40 35
35 30
30
25
25
20
20
15
15
10 10
5 5
0 0
1-ene
6-ene
11-ene
16-ene
21-ene
26-ene
31-ene
Tiempo (días)

mínimos del aire durante el 01 de Enero y 31 de Enero del 2014, con una media
máxima mensual de 33,68 ºC y media mínima mensual de 19,10 ºC; el día 23 de
Enero ocurre la temperatura más elevada del mes de 41,17 °C y el 26 de Enero se
presenta la mínima del mes de 10,26 ºC.En la figura 4 se presenta conjuntamente
la temperatura del suelo, medida a una profundidad de 30 cm, observándose
fluctuaciones normales con la variación de la temperatura del aire debido a lo
superficial de la lectura, con un promedio de 31,47Cº.
Unidad 4
Nivel freático
0.090
0.080
387.66
0.070
387.64 0.060
Lámina lluvia (mm)
387.62 0.050
0.040
387.60
0.030
387.58 0.020
387.56 0.010
0.000
387.54
-0.010
387.52 -0.020
Ilustración 6-17 Variación del nivel freático en la unidad 01. 01-14 hasta 31-01-14
116

debe mencionar que la unidad registró valores desde el 01 de Enero hasta el 31 de
Enero de 2014, alcanzando un valor máximo de 387,66m y un valor mínimo de
Enero 2015
Unidad 1
0.320 50
0.300
0.280 45
0.260 40
0.240

0.220 35
0.200 30
Humedad
0.180
0.160 25
0.140
0.120 20
0.100 15
0.080
0.060 10
0.040 5
0.020
0.000 0
Tiempo (días)
Ilustración 6-18 Medición de los sensores EC-5 2,70m 5-TE a 0,30 m, 10-HS a 1m. 01-01-
15 hasta 31-01-15

Presentando algunas ligeras fluctuaciones propias de las precipitaciones
acontecidas, pero estabilizándose en valores normales a los presentados en los
meses anteriores lo que descarta una posible pérdida de la conducción de Aguas
Cordobesas en este sector.
Los sensores EC-5 y el 5-TE (Ilustración 6-18) ubicados a 2,70 y 0,30 m

respectivamente de la superficie del suelo muestran fluctuaciones temporales en la
variación del contenido de humedad del suelo producto de las precipitaciones. En
lo referente a los valores propiamente dichos, el contenido de humedad, presentó
valores medios de 21,70% y 20,00% de humedad durante el periodo monitoreado.
Alcanzado el sensor EC-5 un pico máximo de humedad de 23,00%, mientras que el
sensor 5-TE alcanza un valor de máximo de 21,00%.
117
Unidad 2

0.320 50
0.300
0.280 45

0.260 40
0.240
0.220 35
Humedad
0.200 30
0.180
0.160 25
0.140
0.120 20
0.100 15
0.080
0.060 10
0.040 5
0.020
0.000 0
29-dic 6-ene 13-ene 21-ene 28-ene 4-feb
Tiempo (días)
Humedad 10-HS Humedad 5TE Lámina Lluvia
Ilustración 6-19 Medición de los sensores, 10HS a 1m y 5TE 0.30m. 01-01-15 hasta 31-
01-15

Presentando valores con ciertas fluctuaciones propios de la época estival de lluvias
pero sin valores extremos, lo que descarta una posible pérdida de la conducción de
Aguas Cordobesas en este sector.
El sensor 5-TE (Ilustración 6-19) presenta fluctuaciones similares al sensor 10-HS

propio de esta época del año. En lo referente a los valores propiamente dichos, el
contenido de humedad, presentó un valor medio de 15,20% durante el periodo
monitoreado. Alcanzado el sensor 5 TE un pico máximo de humedad de 17,40%.
118
Unidad 3

50
0.26
0.24 45
0.22 40
0.2

0.18 35
0.16 30
Humedad
0.14 25
0.12
0.1 20
0.08 15
0.06 10
0.04
0.02 5
0 0
Tiempo (días)
Ilustración 6-20 . Medición de los sensores, 10HS a 1m y EC-5 1.80 m. 01-01-15 hasta 31-
01-15
Los sensores 5 TE y EC-5 ubicados a 0,30 m y 2,00 m, muestran ligeras

variaciones de humedad, producto de las precipitaciones acontecidas durante el
periodo de mediciones, presentando el sensor 5 TE, una humedad promedio de
22,70% con desviación de 1.16%.
Para el sensor EC 5 presenta una humedad promedio de 16.24% y una

desviación de 0,94%. Estos valores bajos se deben a que la zona de la unidad 3
presenta un suelo antrópico, muy compactado y con poca vegetación, manteniendo
humedades casi residuales donde las variaciones de humedad son muy escasas.
119
Temperatura Máxima y Mínima diaria del aire y Temperatura del

suelo
40 40
35 35

30 30
25 25
20 20
15 15
10 10
5 5
0 0
6-ene
1-ene
26-ene
11-ene
16-ene
21-ene
Tiempo (días)

mínimos del aire durante el 01 de Enero y 31 de Enero del 2015, con una media
máxima mensual de 30,13ºC y media mínima mensual de 16,84 ºC; el día 06 de
Enero de 2015 ocurre la temperatura más elevada del mes de 36,00 °C y el 02 de
Enero se presenta la mínima del mes de 9 ºC.
En la Ilustración 6-21 se presenta conjuntamente la temperatura del suelo, medida

a una profundidad de 30 cm, observándose fluctuaciones normales con la variación
de la temperatura del aire debido a lo superficial de la lectura, con un promedio de
25,85Cº.
Conclusiones de los sistemas de monitoreos en las 4 Unidades del

Loteo Residencial Cañuelas.
Entre las mediciones del mes de Enero y Abril 2013, se puede observar en la
Unidad 1, la humedad de acuerdo a los sensores respectivos, que miden valores
similares, con ciertos picos causados por la lámina de lluvia que en ambos meses
registró un pico cercano a los 25 mm.
La unidad 2, no se obtuvo registros por cuanto el instrumental presento problemas

técnicos.
La unidad 3, los valores de registro de humedad se mantuvieron estables con una

tendencia neutral, que obviamente registró el mismo valor de lámina de lluvia
referido en la unidad 1.
120
En cuanto a las mediciones de temperatura diaria máxima y mínima, en Enero

2013, registró una media máxima superior a los 30ºC y la media mínima inferior a
los 20ºC, exactamente de 18º,36 C.
En la Unidad 4, el nivel freático que muestra la gráfica en el mes de Enero, no fue

inferior a los 386 msnm. En cambio en el mes de Abril, el valor más bajo fue de 388
msnm. En el mes de Enero, la lámina de lluvia registró picos de 10 y 26 mm,
mientras que en Abril, la mencionada lámina de lluvia registró picos de 30 y 35 mm.
Con lo cual, se deduce que el aumento de la cota freática en Abril se debió a la

lámina de lluvia que presentó dos picos importantes en el mes.
En el mes de enero 2014, en la Unidad 1, la humedad registra una leve tendencia

negativa, que fue distinta a ese mismo mes del año 2013, donde existió una
tendencia neutral, es decir, ni positiva ni negativa.
La unidad 2, debido a la compra del instrumental nuevo, recién empezó a registrar

lecturas a partir del 15 de enero, por lo que no amerita comparaciones.
En la unidad 3, conservó una tendencia estable con valores promedios cercanos al

12%.
La temperatura, registró una media máxima mensual de 33ºC, mientras que la

media mínima mensual, apenas estuvo por debajo de los 20ºC. Valores de
temperatura mayores al año 2013.
La unidad 4, el nivel freático mostro minimos cambios, que no llegaron a completa

la unidad. Esto debido a que la lámina de lluvia, fue muy escasa no llegando ni a
valores de 1mm.
En Enero del 2015, en la unidad 1, la humedad registra una tendencia positiva

poco pronunciada, cuyo origen se debe a que la lámina de lluvia registró 3 picos
importantes cuyos valores son superiores a 35 mm en los tres casos.
En la unidad 2, de igual forma, se registra una leve tendencia positiva en los

valores de humedad con un 14%.
En la unidad 3, los registros muestran valores de 22%, que hacen de los graficos
de humedad tengan de igual forma, un leve crecimiento.
En lo que respecta a la temperatura, los valores de la media mensual, indican una

temperatura de 30ºC, como media máxima y 16ºC como media mínima.
Comparando con los anteriores años, 2013 y 2014, en el 2015 existieron valores
mas bajos.
121
El conjunto pozo de observación y sensor denominado Unidad N°4 fue instalado en

la zona cercana al Club House y su ubicación georreferenciada es latitud
31°28'23,71" S, longitud 64°9'05,30"O.
Cota Nivel Freàtico Unidad Nº4 Cañuelas

Nùmero de Registros
0 100 200 300 400 500 600 700 800 900
30 418.6
Nivel Freàtrico(msnm)
25 418.4
Precipitaciòn(mm)
20 418.2
15 418
10 417.8
5 417.6
0 417.4
Días
Lluvia Diciembre 2012 nivel
Ilustración 6-22 Gráfico de monitoreo Pozo Cañuelas Córdoba
En la Ilustración 6-22 se observa los valores del nivel freático en la zona. Cerca de
la zona de ubicación del pozo de monitoreo de freática se encuentra ubicada la
planta de tratamiento de líquidos cloacales dicha planta durante el año 2012
volcaba los líquidos tratados en dos pozos que se ubicaban en la cercanía de la
misma.
Para evaluar la vulnerabilidad del acuífero se utiliza para la determinación del

acuífero libre la metodología de Foster e Hirata mediante la aplicación del “Sistema
DiOS para la Evaluación del Índice de Vulnerabilidad del Acuífero”.
Para los autores la vulnerabilidad a la contaminación del acuífero es la sensibilidad

de un acuífero a ser adversamente afectado por una carga contaminante. Esta
metodología permite evaluar la vulnerabilidad de un acuífero mediante una
selección de parámetros de campo disponibles y que le otorgan al cálculo un alto
grado de certeza.
122
La vulnerabilidad del acuífero es una función de:
1- La inaccesibilidad de la zona saturada en un sentido hidráulico, a la penetración

de contaminantes.
2- La capacidad de atenuación de los estratos encima de la zona saturada del

acuífero como resultado de su retención física y reacción química con
contaminantes.
Estos dos componentes de la vulnerabilidad del acuífero interactúan con los

siguientes elementos correspondientes de la carga contaminante al subsuelo:
1.- El modo de disposición del contaminante en el subsuelo.
2.- La clase de contaminante en términos de su movilidad y persistencia.
Esta interacción determinará el tiempo de residencia en la zona no saturada y la

demora de la llegada del contaminante al acuífero. Se realizó una evaluación
tomando en cuenta las características del sitio y la posición del acuífero libre en el
predio.
Para ello se consideran las siguientes características:
a) Distancia al agua: 30 m.
b) Ocurrencia del agua subterránea: acuífero no confinado libre.
c) Sustrato litológico: Sedimentos no consolidados: Limos arenosos.
Mediante el sistema DiOS, se obtiene un índice de 0,18; lo que indica que la

vulnerabilidad de contaminación del acuífero libre es baja.
Por otro lado al analizar los datos obtenidos del monitoreo se encontra los
siguientes:
Cota del Pozo: 440 msnm, se obtienen 809 registros separados por intervalos de
15 minutos, el nivel estático mínimo está situado a 417,51msnm y el nivel estático
máximo se ubica a 418,49 msnm.
A los resultados obtenidos con los sensores del pozo referido, se realizaron los
análisis de datos correspondientes con la finalidad de interpretar de mejor forma el
monitoreo de sus niveles estáticos.
123
Nivel Freático %Acumulado

120.00%
100.00%
Acumulado(%)
80.00%
60.00%
40.00%
20.00%
0.00%
418.2 418.3 418.1 418 417.7 417.6 417.8 417.9 418.5 417.3 417.4 417.5 418.4
Niveles(msnm)
Ilustración 6-23 Distribución Normal Datos Pozo Cañuelas
En esta ilustración se observa el nivel freático vs el porcentaje acumulado de los

datos obtenidos en el pozo de Observación en el Loteo Cañuelas, en donde se
agregan las cotas de niveles freáticos.
De acuerdo a esto, indica una muestra de datos están cercanos a la normal.
A continuación se muestra los gráficos que se obtienen del análisis de datos por
regresión lineal.
Gráfico de probabilidad normal

418.6
418.4
418.2
Y(msnm)
418
417.8
417.6
417.4
0 20 40 60 80 100 120
Muestra percentil
Ilustración 6-24 Gráfico de datos probabilidad Normal
124
Variable X 1 Gráfico de los residuales

0.4
0.3
0.2
0.1
0
Residuos
-0.1 0 100 200 300 400 500 600 700 800 900
-0.2
-0.3
-0.4
-0.5
-0.6
-0.7
Variable X 1
Ilustración 6-25 Análisis de Residuales
Variable X 1 Curva de regresión ajustada

418.6
418.4
418.2
Y(msnm)
418
Y
417.8 Pronóstico para Y
417.6
417.4
0 200 400 600 800 1000
Variable X 1
Ilustración 6-26 Curva de regresión ajustada
En la Ilustración 6-24 está el grafico de la probabilidad normal para analizar valores

no ordenados de un grupo de datos que se ha obtenido en el pozo Nº 4 de
Cañuelas. Este gráfico se usa para auxiliar a clasificar si una muestra de datos
numéricos proviene o no de una distribución normal. De no ser el caso,
comúnmente se puede determinar el tipo de alejamiento de la normalidad
examinando la forma en la que los datos se desvían de la línea de referencia
normal. ( StatPoint, Inc., 2007)
125
En la Ilustración 6-25 el gráfico de los residuales vs cada variable, en este caso X1

se utilizan para detectar valores atípicos y la no semejanza de varianzas. Este tipo
de gráficos permite detectar si una variable en modelo debe incluirse de manera
lineal o no. Cuando los datos se desvían de la linealidad, la variable no debe entrar
en el modelo de manera lineal y en este caso el gráfico indicará como transformar
los datos para alcanzar la linealidad.
La Ilustración 6-26 la curva “de regresión ajustada” se observa que la correlación

es negativa ya que al aumentar X, Y disminuye, cabe mencionar que el coeficiente
de correlación calculado por el sistema siempre es positivo, por lo cual se debe
basar la gráfica de regresión para determinar el signo.
6.4 Modelación de transporte de Contaminantes en área circundante

Urbanización Cañuelas.
Ilustración 6-27 Área de modelación de Transporte
En este trabajo se desea resolver un problema de transporte de contaminantes,

simulando el efecto de una fuente de descarga puntual de combustible, como
podría ser el caso de un accidente con un camión repartidor de gasolina en malas
condiciones, que derrame el contenido en la vía y por las cunetas, llega hacia el
acuífero adyacente.
En este trabajo se estudia un acuífero de dimensiones 1500m x 1500m, que dan un

producto de 2,25 Km² de superficie. El acuífero tiene una conductividad (K) de 0,22
m/día, un espesor saturado inicial (b) de 30 m, porosidad total (n) de 0,52 y un
gradiente hidráulico (iH) de 0,5%. Este gradiente hidráulico puede ser simulado con
126
dos condiciones de nivel constante (30 m en el lado izquierdo y 23 m en el lado

derecho).
El vertido descarga al acuífero una concentración de 500 mg/L. Otro dato

importante es que se supone un acuífero inicialmente sin contaminación y la pluma
se desplazará por efectos de advección y difusión. La dispersividad longitudinal se
asume conocida (aL = 10 m), mientras que la relación entre las dispersividades
horizontal y longitudinal es uno (1). Se ubica al norte un pozo de observación de
calidad de agua subterránea a una distancia de 300 m.
Ilustración 6-28 Concentración-tiempo
Ilustración 6-29 Pluma contaminante-Descarga Puntual
127
7 Análisis de Resultados
7.1 Cuenca del Milagro

Proceso de Levantamiento de información de pozos existentes
En la Universidad Estatal de Milagro, existen 4 pozos de los cuales, el pozo

ubicado en las coordenadas UTM 655465S y 9762252E, posee bomba sumergible
con un caudal 0,002 m³/s. que sirve para regar las áreas verdes cercanas a dicho
pozo, que son aproximadamente 500 m². Se denomina para efectos de
identificación del proyecto, Pozo 1 Garita.
Este pozo está en la cota topográfica 10,88 msnm, el nivel estático del pozo se
encuentra al momento del estudio, en 14,95 m. de profundidad, lo que nos da:
Z=10,88; menos h=14,95; total H= - 4,07m. Este nivel piezométrico obviamente
varía y en su respectivo cuadro de registro se ha visto, realizando la prueba de
bombeo, desciende 0,55 metros respecto del nivel estático del agua. En el
monitoreo se puede apreciar los valores de descenso al momento del bombeo y el
tiempo de recuperación de dicho nivel estático, tomando en cuenta intervalos de 15
minutos que es el tiempo que se calibró para el registro de los sensores.
Existe otro pozo a una distancia en línea recta de 197 m. aproximadamente, que no
tiene uso específico, se ubica en las coordenadas UTM 655346S y 9762106E, que
será para efectos del estudio un pozo de observación. Este se denomina dentro del
proyecto como pozo 2 Bloque K y está en la cota topográfica 9,60 msnm, el nivel
estático del agua está a la fecha del estudio a 13,17 m, calculando se obtiene: Z=
9,60m, h=13,17; que nos da un H= -3,57m.
128
Ilustración 7-1 Mapa satelital de la ubicación de pozos de estudio (Foto Google Earth)
Para efectos de incluir otra referencia, se mide el nivel estático del pozo ubicado
junto al Bar del Polideportivo de la Universidad, de coordenadas UTM 655611S
9762552E con un nivel estático de 11,80 m, la cota es de 9,83msnm, entonces se
calcula el H:
Z= 9.83m, h=11,80; que nos da un H=-1.97m.
En el pozo número 4 existe una estructura de un tanque elevado de 15 m. de

altura que abastece de agua a la Universidad, se extrae el mencionado líquido con
una bomba sumergible de 11.000 Vatios, con un caudal de 0,02 m³/s, está ubicado
en la cota 9,72 msnm, cuyo nivel estático 13,39, entonces tiene un H=-3,67 msnm.
Como el pozo está en producción, su nivel dinámico es de 14,93; es decir, que
1,54m. es el descenso(s).En la zona de estudio se bombean en 2 pozos, el pozo 1
Garita y el pozo 4 Tanque Elevado, con un total de 525.600 m³ al año
aproximadamente.
129
Cota Prof. Nivel Cota

Nombre del Pozo Topográfica estático Piezométrica Diferencia
Pozo 1(Pileta
Ornamental) 10,88 14,95 -4,07
Pozo 2(Bloque K) 9,60 13,17 -3,57 -0,50
Pozo 3
(Polideportivo) 9,83 11,8 -1,97 -2,10
Pozo 4(Tanque
elevado) 9,72 13,39 -3.67
Tabla 7-1 Pozos del Área de Estudio
Inmediatamente se muestran los resultados obtenidos del modelo tras la

calibración, cuyo ajuste puede ser considerado como aceptable y responde
satisfactoriamente a los objetivos propuestos.
El transcurso de la calibración, fue repeticiones sucesivas y radicó en alcanzar el

balance entre los ingresos, la variación del almacenamiento en el suelo que se
manifiesta en los cambios que soporta el nivel freático y los egresos para cada
hipótesis. Los parámetros de calibración utilizados para alcanzar el ajuste entre los
valores de carga hidráulica observados (o considerados de acuerdo a las
condiciones superficiales del flujo) y simulados fueron la conductividad hidráulica y
la recarga.
Los parámetros calculados son coherentes con el modelo conceptual propuesto y

con las estimaciones previas, los valores de conductividad hidráulica estimados
inicialmente no difieren de forma sustancial de los que surgen de la calibración;
están dentro del mismo orden de magnitud y mantienen la misma relación de
anisotropía.
En relación a las condiciones de frontera consideradas se verificaron y contrastaron

las diferentes hipótesis de funcionamiento:
A - no existe flujo a través del perímetro del área de estudio.
B- frontera del tipo Recarga y Evapotranspiración en la zona de estudio, basado en

los datos de la estación meteorológica de Milagro, administrada por el Instituto
Nacional de Meteorología e Hidrología (INAMHI).
Haciendo un detalle de los rasgos geológicos de la cuenca referida y su zona de

influencia se considera que al acumularse los sedimentos aluviales, son suelos
con gran capacidad de almacenar volúmenes de agua importantes, por ser de gran
permeabilidad. En relación de este hipotético escenario, el nivel del acuífero
freático dentro de la subcuenca y por consiguiente área de estudio, presenta gran
130
capacidad de recuperación el acuífero. Por este motivo se adopta el modelo con

las cargas hidráulicas iniciales, que han sido promediadas de los niveles estáticos
medidos en sitio. Para poder encontrar el valor de la conductividad, entre el ajuste
inicial y el ajuste final existió una variación de alrededor 10 puntos en los
parámetros de calibración del modelo (medianamente significativa) que permite
representar las situaciones con caudales altos que atraviesan la frontera de
simulación.
Conductividad Kx Ky Kz
(m/s) (m/s) (m/s)
arena-limosa 0.00011 0.00011 0.00011
Recarga RCH en el area de estudio

Condiciones de contorno Evapotranspiracion EVT promediada
Rehumedecimiento de celdas secas activo

Umbral de humecimiento 0.2
Intervalo de Humedimiento 5
Metodo de humedimiento Desde abajo WETDRY }0
Carga de humedecimiento A partir del umbral IHDWET }0
Factor de humedimiento 1
Carga en celdas secas -1.00E+30
Layer N1 type No confinado 1
Layer N2 type Confinado/No confinado var S y T 3

Tipo de Corrida Transitoria
Modflow 2000
Solver WHS
Max outer iterations(MXITER) 50
Max inner iterations(ITER1) 25
Head change criterion(HCLOSE]) 0.01
Residual criterion(RCLOSE) 0.01
Damping factor(DAMP) 1
Relative residual criterion(RSCRIT) 0
Factorization level 0
Tabla 7-2 Parámetros ajustados y métodos utilizados
131
Ilustración 7-2 Diagramas Ingreso-Egresos T= 500 días
Con referencia al balance de entradas y salidas del modelo (Ilustración 7-2 e

Ilustración 7-3), se presta atención bajo esta hipótesis, que el mayor intercambio de
caudales ocurre principalmente a través de la recarga que alimenta el acuífero y
pasan a engrosar el volumen de almacenamiento de este y el bombeo de los pozos
de la Unemi.
Ilustración 7-3 Diagramas Ingreso-Egresos T= 1000 días
132
7.1.1 Registros obtenidos en los pozos Cuenca del Milagro.
Cotas Nivel FreáticoPozo 1 Garita
Número de registros
0 500 1000 1500 2000 2500
60.0 -3.0000
50.0 -3.2000
Nivel Freatico(msnm)
-3.4000
Precipitación(mm)
40.0
-3.6000
30.0
-3.8000
20.0
-4.0000
10.0 -4.2000
0.0 -4.4000
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31
Días
Lluvia Febrero Niveles piezometricos
Ilustración 7-4 Nivel freático Pozo 1 Garita
El pozo 1 Garita fue monitoreado con los sensores ya previamente citados con
intervalos de 15 minutos, en total son 2.695 registros. El gráfico demuestra una
tendencia positiva, es decir, se sostuvo un crecimiento de los niveles monitoreados.
También indica cuantos metros desciende, mientras está el pozo en bombeo y
luego de que dicho bombeo cesa, y así se puede obtener el tiempo que tarda en
recuperarse el nivel estático. Se recopiló información pluviométrica, para lo cual se
tomó los datos correspondientes a la Estación Meteorológica INAMHI Milagro.
133
Cotas nivel Freàtico Pozo 2 Bloque K

Numero de registros
0 500 1000 1500 2000 2500
12 -2.3
10 -2.35
Nivel Freatico (msnm)

Precipitación(mm)
8 -2.4
6 -2.45
4 -2.5
2 -2.55
0 -2.6
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31
Días
Mayo COTAS(msnm)
Ilustración 7-5 Nivel Freático Pozo 2 UNEMI
En el pozo 2 ubicado en las coordenadas UTM 655346 y 9762106, parte posterior

del Edificio de Aulas K. Tiene una tendencia más pronunciada que el pozo uno,
pero en un menor rango de tiempo.
Simulación situación actual
Datos: Conductividad: 10 m/d; Recarga: 760mm; Nivel Inicial:-2.61
Ilustración 7-6 Niveles simulados escenario actual
134
Elaboración del modelo computacional Modflow.

Se muestra una síntesis de los pasos seguidos para la reproducción del modelo
computacional por medio del software Visual Modflow. Como primera acción, se
muestran las unidades que se utilizan en el modelo de simulación, para luego
concretar el dominio del sistema de acuíferos que integra el área de estudio y la
extracción de un determinado número de puntos en el espacio tanto en planta
como a lo largo del eje vertical. El modelo de flujo sólo simula el comportamiento
del Acuífero libre, la malla de elementos finitos es tridimensional. En la malla de
diferencias centradas en los bloques, se debe especificar que es el tipo de malla
que utiliza Modflow, la frontera del modelo se sitúa en los bordes de las celdas, el
régimen temporal es transitorio cubriendo un periodo de 1.000 días.
En la Tabla 7-3 se sintetizan los límites de simulación y las características
geométricas de la malla de simulación.
Malla Coordenadas UTM

(x1,y1) límite inferior izquierdo x1 655079.05 m E
y1 9762206.94 m S
(x2,y2) límite superior derecho x2 655649.45 m E
y2 9762646.66 m S
Delta x(Ancho,m) 1000
Delta y(Largo,m) 1000
Delta z 50
ancho de celdas(m) 50
largo de celdas(m) 50
# de col 20
# de filas 20
Tabla 7-3 Geometría y Límites del Área de Estudio Zona 17-UNEMI Milagro
7.1.2 Condiciones iniciales y de frontera
En este ítem, se detallan las condiciones iniciales (de contorno) en las

demarcaciones de la franja de estudio, que fueron uno de los componentes cuya
incidencia es importante a tener en cuenta en las sucesivas pruebas en el
programa, es decir ,cuando se introduce todos los datos y luego se oprime RUN.
También, se debe detallar que las condiciones de frontera se fueron haciendo más
complejas de manera gradual, de acuerdo a como se obtenían escenarios de
equilibrio en el sistema después de la calibración de los distintos parámetros.
La primera hipótesis consistió en considerar que no existe flujo del acuífero a
través del perímetro del Área de estudio.
La segunda hipótesis consistió en imponer una frontera del tipo Recarga y
Evapotranspiración en toda la zona.
Las mediciones realizadas del nivel estático fueron tomadas en cuenta como nivel
inicial.
135
7.1.3 Mediciones de Calidad del Agua
De lo observado en las gráficas de los resultados, la dirección del flujo es en

sentido Noroeste a la zona de flujo, con dirección al cauce principal del Río Milagro
que dista aproximadamente unos 3 Kilómetros de la zona.
A continuación los resultados de las mediciones con el conductímetro es la

siguiente:

Conductividad μS 1174
Salinidad ppm 470
Solidos totales Disueltos ppm 820
Tabla 7-4 Parámetros de Calidad de Agua de los pozo Garita.

Conductividad μS 255
Salinidad ppm 100
Solidos totales Disueltos ppm 170
Tabla 7-5 Parámetros de Calidad de Agua de los pozo Bloque K
Estos resultados indican que el pozo uno, que está siendo utilizado, tiene valores
aceptables para uso doméstico, estos valores son tomados en la parte superior del
nivel del acuífero libre.
136
7.2 SubCuenca del Suquìa
Ilustración 7-7 Ubicación Pozo de Observación Loteo Cañuelas Córdoba-Argentina
De lo mostrado en el gráfico de nivel freático mostrado en el capítulo 6 con los

sensores conocidos, instalados en la unidad 4 , se puede indicar que estos se
mantienen sin variación alrededor de la cota 418msnm, excepto un pico negativo
de alrededor de 1 metro, es decir, cota 417msnm que puede ser causa de la falta
de lluvias en el periodo anterior a este evento, pero luego vuelve a recuperar el
nivel inicialmente mencionado. Ante esto se puede indicar que en el periodo de
monitoreo no existe ningún problema o fallo del conducto Los Molinos que
atraviesa este lote.
137
Ilustración 7-8 Zona de Interés Loteo Cañuelas Superficie 53 Ha aproximadamente.
138
8 Conclusiones y recomendaciones
8.1 Introducción
El actual capítulo abrevia las principales conclusiones generales y particulares,

conseguidas a partir de los modelos desplegados.
Consecutivamente se exhiben recomendaciones generales, particulares y futuras
líneas de estudio que lograrán continuarse en base al presente trabajo.
8.1.1 Conclusiones Generales
Tras los análisis y modelaciones surgen conclusiones generales relacionadas tanto

con el proceso de recolección de información, de conceptualización y simulación
propiamente dicha.
Entre estas se destacan las siguientes:
 Es de suma importancia la unificación de los conocimientos hidrológicos e

hidrogeológicos existentes de cada ambiente, sistematizando la información
de base para hacerla accesible y comparable con futuras investigaciones. El
desarrollo de modelos hidrogeológicos de mayor precisión requiere datos de
campos actualizados y representativos de los ambientes de estudio.
 La información de los pozos existentes debe estar georreferenciada y en un
sistema unificado de manera que permita que se puedan incorporar a los
modelos matemáticos.
 La construcción de los modelos y la simulación de distintos escenarios, permite

estudiar el comportamiento hidráulico de los acuíferos y establecer criterios
para la planificación del aprovechamiento sostenible del agua subterránea.
 En los sistemas acuíferos las condiciones pueden variar notablemente en

forma espacial (variación de conductividad, porosidad, grado de saturación,
etc.) o en forma temporal, grado de saturación y condiciones de borde, en
consecuencia para descifrar o predecir el comportamiento de un área o
situación particular, se deben asumir esquemas o modelos con niveles de
detalle para los que, en general, no se posee tal grado de conocimiento.
 Los modelos de simulación constituyen herramientas muy poderosas para

el manejo y gestión de las aguas subterráneas. Un modelo intenta reproducir
la realidad de la forma más acertada posible en función de la información
disponible, lo que permite realizar estimaciones o análisis comparativos de
distintos comportamientos posibles. La aplicación de MODFLOW para los
casos que se presentan en esta tesis muestra que esta herramienta puede
aplicarse para estudios locales como para estudios más generales y que los
resultados que se obtienen permiten adoptar acciones tendientes a un
correcto manejo de las cuencas subterráneas.
 Los modelos de flujo propuestos resultan en términos generales coherentes

con la información existente y reproducen cualitativamente las curvas
isofreáticas de los sectores donde existe información.
139
 Las limitaciones observadas en los modelos surgen de la escasez de

información y un bajo grado de precisión en muchos de los datos existentes.
Existe necesidad de aumentar la información hidrogeológica disponible.
8.1.2 Conclusiones sobre la Cuenca del Milagro
De acuerdo a los resultados del modelo en el capítulo respectivo se llegaron

a las conclusiones que inmediatamente se reproducen:
 Los resultados del modelo numérico, representan adecuadamente el

sistema hidrogeológico.
 Los escenarios planteados muestran que el acuífero tiene buena respuesta

al bombeo en cuanto a intervalos largos, como sucede con el pozo de
Plazoleta. El pozo Bloque K actúa como pozo de observación. Un tercer
pozo llamado Tanque Elevado se tomó la medición de los niveles estáticos
para comprobar y sacar conclusiones, respecto de hacia dónde drena o
circula el agua subterránea. Después de haber monitoreado los pozos, es
que tienen una buena capacidad de recuperación luego de los ciclos
bombeos realizados. Además el pozo que abastece al tanque elevado es
donde se presentan el nivel estático más bajo por ser una bomba
sumergible de mayor capacidad que sirve para que haya agua en todas las
edificaciones de los predios universitarios, por lo tanto, las medidas del nivel
estático indican que el pozo de observación está ubicado a una cota más
superficial que el nivel descrito del pozo de bombeo, sin embargo el pozo
cercano al polideportivo tiene una mayor cota sobre el nivel del pozo de
tanque elevado mencionado.
Los escenarios que se plantean adicionalmente, son el aumento de pozos

en la zona, agregando 2 bombas de 0,002 m³/s de caudal, dado que la
ciudadela universitaria continuará creciendo y se proyecta al 2050, tener una
infraestructura que doble a la actual. Hay que aclarar que, la mayor
demanda que va a tener el agua extraída de los pozos será para riego de
sus áreas verdes y de sus canchas de césped natural, que es precisamente
las áreas que se acrecentaran.
Adicionalmente, otro escenario se daría si las aguas residuales de algún

pozo o tanque séptico contamina el manto freático, que al realizar
excavaciones y calicatas se encuentra en la mayoría de los casos a 2.80 y
3.00 metros de la superficie. Existe la simulación de este escenario, que nos
indica que la pluma contaminante predomina en la dirección del acuífero
libre.
 Frente a extracciones de gran proporción mencionadas, que es el escenario

futuro, ocurre que los niveles descienden aproximadamente 1,5 m. en el
140
caso del pozo del Tanque elevado y 0,50 metros en el pozo ubicado en la
garita de ingreso a los predios universitarios.
8.1.3 Conclusiones sobre el sistema Cuenca del Suquía.
Al analizar este pozo, que como se ha descrito está ubicado en una

urbanización, se observa que tiene una tendencia estable desde la fecha
que se lo registra.
 El análisis de los registros revelan que son homogéneos, inclusive en un par

de días que registro un pico negativo con una variación de 1,20 m. pero
aceptables, lo que no es suficiente para hacer cambiar radicalmente la
tendencia.
 Estos mismos datos en cambio no son independientes, el histograma indica
que los valores con mayores repeticiones son los situados entre los valores
418,20 y 418,30
 El gráfico de residuales muestra uniformidad excepto unos valores que
causan un pico negativo, ya descrito en párrafos anteriores.
 Los resultados muestran que el monitoreo de niveles freáticos en esta zona
y su modelación permiten considerar los futuros escenarios vinculados a los
efectos de las afectaciones antrópicas.
8.2 Recomendaciones
8.2.1 Cuenca del Milagro
De acuerdo a lo que se ha analizado, la zona de estudio posee gran espacio de

áreas verdes, áreas que son permeables y en donde la precipitación que cae se
infiltra fácilmente al subsuelo, lo ideal sería tener más pozos de observación en
donde se pueda medir y registrar cómo evolucionan los niveles estáticos. Esta
Ciudadela Universitaria posee pocas áreas impermeables, se puede describir que
de las 19 Ha. en total que tiene la Universidad, un 20 % es área impermeable y el
restante 80% es área permeable, por lo que existen áreas de potencial infiltración.
Se observa en la simulación en la cual existe un incremento de los pozos de

bombeo, el nivel estático baja notoriamente. Para esto se recomendaría
incrementar un solo pozo adicional con una bomba que de acuerdo a la magnitud
del riego, su caudal sea el más óptimo. Con esto se mejora el tiempo de
recuperación de los niveles estáticos del agua.
Como propuestas generales de programas y acciones a implementar para alcanzar

el objetivo del desarrollo sustentable de las aguas subterráneas, se indican las
siguientes:
• Construcción y operación de una red de monitoreo hidrogeológico en cada

subcuenca.
141
Con este propósito deberán distribuirse los elementos de medición de forma

racional en la superficie de cada subcuenca, de forma que sus mediciones ayuden
a determinar de manera inequívoca y representativa de vastos sectores, las
direcciones del flujo, la variación del gradiente y sirvan al discernimiento de las
condiciones de contorno admitidas.
• En lo referente a información hidrogeológica se considera conveniente la creación

de un banco con datos de: perforaciones (coordenadas geográficas, perfiles,
caudales, uso, niveles, composición del agua, rendimientos, etc); clima
(precipitación, temperatura, humedad, viento, insolación, radiación, etc); geología
(informes, mapas, perfiles); geomorfología (geoformas); hidrología superficial
(caudales y derrames de los ríos y arroyos y alturas de las lagunas); edafología
(tipos y distribución de los suelos); riego (volúmenes aplicados, origen del agua,
pérdidas); piezometría (niveles de agua subterránea referidos al sistemas
georreferenciados y su evolución en el espacio y en el tiempo); uso del agua
(consumo humano, riego, ganadero, industrial).
• Se considera apropiado la realización de estudios que permitan ajustar los

parámetros de recarga desde los cauces de los ríos. Esto reduciría la fluctuación
asociada al valor de la conductividad de la capa que vincula el acuífero con el lecho
de los cauces.
• Anticipar sobre políticas y estrategias para el manejo del agua subterránea fijando
guías sobre la protección natural de los acuíferos y los criterios más apropiados
para el manejo sustentable del recurso.
• Evaluar en ambos sistemas la vulnerabilidad de los acuíferos estudiados

mediante las metodologías existentes (DRASTIC, GODS, AVI, etc…). Dado el bajo
ritmo de renovación de las aguas subterráneas, su depuración se torna muy
embarazosa. Es necesario en consecuencia evitar el manejo incorrecto de las
aguas subterráneas y el conocimiento de los sectores con mayor vulnerabilidad
ayuda a fijar normas y medidas de conservación.
• Se propone la generación de un mapa hidrogeológico para cada sistema que

junto con un banco de datos dinámico de libre acceso en Internet, que sea una
herramienta valiosa, para facilitar la toma de decisiones.
8.2.2 Cuenca del Suquía
En la urbanización Loteo Cañuelas, es importante ampliar mucho más tiempo el

monitoreo del pozo y de pozos que estén cercanos, para contrastar los datos,
sobre todo en meses donde haya precipitaciones pluviométricas que permitan
sacar conclusiones del comportamiento del acuífero libre mientras existe la
temporada lluviosa.
De igual manera es interesante monitorear en la temporada seca y registrar los

niveles para analizar todos estos datos.
142
Estos monitoreos se continúan realizando pero debería incorporarse la extracción

de muestras y el análisis de las mismas para verificar si las acciones antrópicas
modifican los parámetros dado que este sector se encuentra muy afectado por la
expansión de las urbanizaciones.
143
9 Bibliografía
Instituto Ecuatoriano de Recursos Hidraulicos-Inerhi. (1986). Plan Nacional de Recursos Hidraulicos
del Ecuador-Recursos Subterráneos.
Instituto Nacional de Meteorología e Hidrologìa. (2013). Cuenca del Guayas Hidrogeología_Mapa

de Izopiezas. Quito.
StatPoint, Inc. (04 de Abril de 2007). http://www.statgraphics.net/. Recuperado el 16 de julio de

2015
Web www.ign.gob.ar-Instituto Geografico Nacional de la Republica Argentina. (s.f.). Instituto

Geografico Nacional de la Republica Argentina. Obtenido de www.ign.gob.ar
A., Burgos J. J. y Vidal. (1951). “Los climas de la República Argentina según la nueva clasificación de
Thornthwaite”Meteoros I (1). Buenos Aires.
Anderson, M. P. y Woessner, W. W.,. (1992). “Applied Groundwater Modeling”. Academic Press,

San Diego.
AQTESOLV. (s.f.). AQTESOLV. Obtenido de www.aqtesolv.com
Argentina, M. d. (s.f.).
(http://www.mininterior.gov.ar/provincias/cuencas_riosali.php?idName=provincias&idNa
meSubMenu=provinciasCuencas&idNameSubMenuDer=provinciasCuencasCooComites.
Recuperado el 07 de 09 de 2013, de
(http://www.mininterior.gov.ar/provincias/cuencas_riosali.php?idName=provincias&idNa
meSubMenu=provinciasCuencas&idNameSubMenuDer=provinciasCuencasCooComites
Auge M. P. y Hernández M. A. (1984). “Características geohidrológicas de un acuífero

semiconfinado (Puelche) en la Llanura Bonaerense. Su implicancia en el ciclo hidrológico de
las llanuras dilatadas”. Coloquio Internacional sobre Hidrología de Grandes L.
Auge, M. (2004). Regiones Hidrogeológicas de la República Argentina y Provincias de Buenos Aires,

Mendoza y Santa Fé. La Plata.
Auge, M. P. ( 2009). “Hidrogeología de Llanuras” Apunte Cátedra de Hidrogeología. Universidad de

Buenos Aires.
B.S.L.C Consultores Cía. LTda. (2013). Consultoria para estudios definitivos para la construccion del
Edificio de Idiomas Unemi. Guayaquil.
Balek, J. (1988). Groundwater recharge concepts. Reidel.
Banco Mundial, T. A. (2000). Water Resources Management Policy and Issues Thematic Annexes.
Barragan Chang E. (2002). analisis y Evaluacion de la Salinidad en la Zona de Taura. Guayaquil:

Ecuela Superior Politecnica del Litoral.
144
Bear, J. (1988). Dynamics of fluids porous media.
Benites Cardenas, A. (2002). investigacion del acuifero de chobo. En Mapa de superficies

piezometricas en el Area Milagro-Chobo y sus aplicaciones en la determinacion de las zonas
de captacion de agua subterranea. Guayaquil.
Benitez, A. (1972). Captación de Aguas Subterráneas. Madrid: Dossat S.A.
C.A.A.A.S., C. A. (1975). Recursos de Agua Subterránea y su Aprovechamiento en la Llanura

Pampeana y en el Valle del Conlara. Tomos I y II. Buenos Aires-Argentina: Instituto Nacional
de Ciencia y Técnica Hídricas -Instituto Federal de Geociencias y Recursos Naturales.
Carrera J.y Neuman. (1986). Estimation of aquifer parameters under transient and steady state
conditions:.
CELEC. (2012). www.hidronacion.org. Recuperado el 4 de Abril de 2013, de Sitio Web de CELEC.
Center, International Groundwater Modelling. (1993). IGWMC. Obtenido de

www.igwmnc.mines.edu
Centro de Investigaciòn y Estudios en Medio Ambiente. (2003). "Curso de Modelaciòn de acuiferos".

Universidad Nacional de Ingenieria de Nicaragua.
Ci66j, Apunte de Cátedra. (2006). “Modelación numérica de aguas subterráneas - Modelación

numérica acuífero zona norte de Santiago”. Santiago de Chile: Universidad de Chile -
Facultad de Ciencias Físicas y Matemáticas - Departamento de Ingeniería Civil.
Cioccale M. A. (1998). “Geomorfología de las Sierras Chicas”. Tesis Doctoral.
Comision Ambientalista Xanaes. (s.f.). Obtenido de www.coamxa.fullblog.com.ar
de Marsily, G. (1986). Quantitative Hydrogeology Groundwater Hidrology for Engineers. London:

Academic Press Inc.
Dinesh Kumar, M., & Singh. (2008). How Serious are Groundwater Over-Exploitation Problems in
India ? A Fresh Investigation Into An Old Issue. Nueva Delhi.
Donado Garzón, L. (2002). Aplicación de Visual Modflow para el modelamiento numérico del flujo
de agua subterránea en el acuífero Morroa(Sucre-Córdoba). VI Congreso Colombiano de
elementos finitos y modelamiento numérico., Bogota.
Dugas F. (1987). Guía Geologica del Suroeste de la Costa Ecuatoriana.
Espinoza C. (2012). Curso electivo de Magíster en Recursos y Medio Ambiente Hídrico. Seminario:
Modelación de Aguas Subterráneas. Santiago De Chile: División de Recursos Hídricos y
Medio Ambiente. Ingeniería Civil de la Facultad de Ciencias Físicas y Matemáticas.
Universidad de Chile.
Espinoza, C. (2012). Curso electivo de Magíster en Recursos y Medio Ambiente Hídrico. Seminario:
Modelación de Aguas Subterráneas. Santiago De Chile: División de Recursos Hídricos y
145
Medio Ambiente. Ingeniería Civil de la Facultad de Ciencias Físicas y Matemáticas.
Universidad de Chile.
Facultad de Ciencias Fisicas y Matematicas de la Universidad de Chile. (2010). Modelaciòn de aguas

subterraneas-Introducciòn a los modelos de simulaciòn. Universidad de Chile.
Ferral A.;Zappino R.;Alaniz E.y Cheniuk S. (2007). Evaluaciòn Hidrogeologica- Diagnostico ambiental
de las aguas subterràneas de la Ciudad de Còrdoba. Parana: V Congreso Argentino de
hidrogeologia.
Fraile, J.; Escuder, R.; Jordana, S.; Ribera, F.; Sanchez-Vila, X.; Vásquez-Suñé, E. (2009).
Hidrogeología Conceptos Básicos de Hidrología Subterránea. Barcelona: Fundación Centro
Internacional de Hidrología Subterránea.
Fulginiti F. (2013). Modelación Hidrogeológica en Cuencas Endorreicas - Aplicación a la Cuenca de la

Laguna Mar Chiquita y la Laguna Llancanelo. En F. Fulginiti, Modelación Hidrogeológica en
Cuencas Endorreicas - Aplicación a la Cuenca de la Laguna Mar Chiquita y la Laguna
Llancanelo. Córdoba: Tesis Maestría en Ciencias de la Ingeniería - Mención en Recursos
Hídricos. U.N.C. Argentina.
García-Arostegui, J., Heredia, J., Murillo, J., Rubio-Campos, J., González-Ramon, A., & Lopez.Geta, J.
(2000). Contribución desde la modelización del flujo subterráneo al conocimiento del
acuífero del Río Verde(Granada). Madrid.
Gobierno de la Provincia de Córdoba. (2012). Obtenido de

http://web2.cba.gov.ar/actual_web/estadisticas/informes_departnuevos/capital/ambiente
/clima.htm
González Bonorino F. (1965). “Mineralogía de las fracciones arcilla y limo del Pampeano en el área
de la ciudad de Buenos Aires y su significado estratigráfico y sedimentológico”. Revista de la
Asociación Geológica Argentina 20(1): 67-148.
Harbaugh A. W.;Banta E. R.;Hill M.C.;Mc Donald M. G. (1992). Modflow 2000,The U.S Geological
Survey Modular Groundwater-User guide modularization Concepts and the Groundwater
Flow Process.
Harbaugh, A.W. (2005). “MODFLOW-2005, The U.S. Geological Survey modular ground-water model
the Ground-Water Flow Process”. . U.S. Geological Survey Techniques and Methods 6-A16.
Hernández, J., Fornero, L., & Martinis, N. (s.f.). Aplicación del Modelo del Flujo de Aguas
Subterráneas Modflow al Acuífero de la Subcuenca El Carrizal y su Extensión al Oasis Norte.
Herrera- Carbajal, S., Villalba-Atondo, A., & Molinar-Tabares, M. (2006). Modelación Hidrogeológica
y Evaluación de la Sustentabilidad del Acuífero Transfronterizo del Río San Pedro en Sonora.
México.
Hill, M. C. (1990). “Preconditioned Conjugate-Gradient 2 (PCG2). A computer program for solving

ground-water flow equations”. U.S. Geological Survey Water-Resources Investigations.
Report 90-4048.
146
IGWMC. (1993). International GroundWater Modelling Center. Obtenido de
http://igwmc.mines.edu/
INAMHI. (2008). CUENCAS HIDROGRAFICAS.
INAMHI. (2013). Informes climatologicos.
Instituto Nacional de Meteorología e Hidrología(INAMHI). (2009). Caracterizacion Hidrogeologica

de la zona sur de la Cuenca Baja del Río Guayas. Quito.
Issar, A. (1990). Regional hydrogeological. Hannover.
J. Sánchez San Román. (Abril de 2008). Departamento de Geologia Universidad de Salamanca.

Obtenido de http://web.usal.es/javisan/hidro
Kinzelbach W. Y Rausch R. (1995). Grundwassermodellierung-Eine Einführung mit Übungen

Gebrüder Borntraeger. Berlin.
Lábaque, M. (2010). PROPUESTA PARA LA CONFECCIÓN DE PLANES DE EMERGENCIA POR ROTURA

DE PRESAS EN LA PROVINCIA DE CÓRDOBA. APLICACIÓN AL DIQUE SAN ROQUE. Cordoba-
Argentina.
Linsley, R. (1882). Hydrology for Engineers. Mc Graw Hill.
MacDonald, M. G., & Harbaugh., A. (1988,1996). En Modflow,A modular three dimensional finite –
difference groundwater flow model (pág. 568). U.S. Geological Survey Techniques of Water-
Resources Investigations,book 6,chapter A1.
Madrid, F. d. (2009-2010). Hidraulica de Captaciones. En Curso de Hidrogeología. Madrid.
Martín, A. P.; Storniolo, A.; Castellano, J. A.; Bejarano, R.; Schejtman, W.,. ( 1997). “Carta
Hidrológica de la Provincia de Santiago del Estero”. I Congreso Nacional de Hidrogeología,
III Seminario Hispano Argentino sobre temas actuales de hidrología subter.
Mc Donald, M.G., y Harbaugh, A.W. (1988). “A modular three-dimensional finite-difference ground-

water flow model”. Techniques of Water-Resources Investigations of the United States
Geological Survey, Book 6.
Ministerio de Recursos Naturales y Energeticos. (1982). Mapa Geologico del Ecuador. Quito.
Moreno L. (2013). Estudio de suelo Edificio Idiomas Universidad Estatal de Milagro. Milagro.
Moyano Miguel A., Leynaud José M. (1994). “Aguas subterráneas para Riego en la provincia de
Córdoba”. Dirección Provincial de Hidráulica.
Pflüger J. C. (2007). El agua potable en la Republica Argentina-Diagnostico de fuentes de agua en

algunas localidades. Informe para el ENOHSA.
PICTOR. (2004). Programa: "Bases para la gestión de los recursos hídricos del NE de la Pcia de Cba.".
Cordoba: Grupo compuesto por UNC, UCC, CEPROCOR e IUA.
147
Plencovich. (2010). “Simulación Hidrológica de los humedales de Mar Chiquita con apoyo satelital”.
Tesis de Maestría en ciencias de la Ingeniería – Mención en Recursos Hídricos. UNC.
Córdoba.
Poeter, E. P. (1997). Inverse models: A necessary next step in groundwater modeling.
Pollock D. W. (1989). Documentation of computer programs to complete and display pathlines using
results from the USGS modular three dimensional finite difference groundwater model .
Open File Report .
Prudic, D. E. (1989). “Documentation of a computer program to simulate stream-aquifer relations

using a modular, finite difference, ground-water flow model”. U.S. Geological Survey Open-
File Report 88–729.
Prudic, D. E., Konikow, L.F., y Banta, E. R. (2004). “A new streamflow-routing (SFR1) package to
simulate streamaquifer interaction with MODFLOW–2000”. U.S. Geological Survey Open-
File Report 2004–1042.
Pulido Bosch, A. (2000). Sobre Explotación de Acuíferos y Desarrollo Sostenible. Publicación en 3er
Taller Internacional sobre Gestión y Tecnologías de Suministro de Agua Potable y
Saneamiento Ambiental, La Habana.
Quintana Salvat F. y Barbeito O. L. (1998). “Carta Geológica Geomorfológica de la Ciudad de

Córdoba y su Entorno” Fotointerpretación. Córdoba, Argentina.: ISSN/ISBN: 03277410.
RemotosCLIRSEN, C. d. (2009). Sistema Nacional de Información. Recuperado el 2 de julio de 2013,

de Sitio Web Sistema Nacional de Información: www.sni.gob.ec
Reyna, S., Reyna, T., & Labaque, M. (2010). Recurso Hídrico Subterráneo. Apuntes de Cátedra,
Córdoba.
Rocca R. J., R. E. (2005). “Características Geotécnicas de los Loess de Argentina”. Rev. Internacional
de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 6(2).
Rodríguez, A. y otros. (2007.). “Estudio Hidrológico – Ambiental Del Sistema Mar Chiquita - Bañados
y Tramo Inferior Del Río Dulce”.Informe Final UNC-UNSE a SSRH. Cordoba.
Rodríguez, F. A. (2007). Modelación del Comportamiento Hidrodinámico del agua subterránea de la

zona comprendida entre Carbó, Pesqueira y Zamora, Sonora. Mexico.
Romero, M., & Mercado, G. (2007). Evaluación del Impacto de la explotación de acuíferos en el
flujo subterráneo en la zona de "El Paso" .
Rushton, K. (1988). Numerial and conceptual models for recharge estimation in arid and semi-arid
zones. Rótterdam.
Saenz, G. M. (1999). Hidrología en la Ingeniería. Mexico: Alfaomega.
Schlagintweit O. (1946). “El Subsuelo de la Llanura Cordobesa en Base a la Perforación S.T.1, con
Referencia al Gondwana”. Dirección Prov. de Minería de Córdoba. Pub. Nº 10. Argentina.
148
SISTEMA NACIONAL DE INFORMACION-SNI. (2009). MAPA DE LA DIVISION HIDROGRAFICA
MILAGRO. El Centro de Levantamientos Integrados de Recursos Naturales por Sensores
Remotos (CLIRSEN).
SMN. (01 de mayo de 2013). www.smn.gov.ar. Recuperado el 2013, de www.smn.gov.ar
SSS- USDA. (1975). Soil Taxonomy.
Stappenbeck, R. (1965). “Geología y Aguas subterráneas de la Pampa. Argentina”.
Suarez E. L. (2010). Acuiferos Semiconfinados y su modelación:Aplicaciones al acuifero de la zona

metropolitana de la ciudad de Mexico. Mexico: Tesis Programa de Postgrado en Ciencias de
la Tierra UNAM.
Subsecretaría de Recursos Hídricos de Argentina. (2004). Estadística Hidrológica de la República

Argentina. Buenos Aires: Evarsa.
Suquía, R. (s.f.). Rìo Suquìa, Provincia de Còrdoba. Obtenido de web2.cba.gov.ar
Therrien R. (1992). Three Dimensional analisys of variably satured flow and solute transport in
discretely fractured porous media. PHD thesis University Waterloo, Waterloo Ontario.
Van der Heijde, Elnawawy. (1993). Compilation of Ground Water Models.
Vasquez, J. B., Miatello, R., Roque, M., & otros. (1979). Geografía Física de la Provincia de Córdoba.
Córdoba,Argentina: Boldt.
Vázquez, J., Miatello, R., & Roqué, M. (1979). Geografía Física de la Provincia de Córdoba.
Villarroya, F. (1994). Regulatory issues mainly about aquifer over-exploitation within the scope of
sustainable development. Publicación, Helsinki.
WWAP, T. W. (2003). Obtenido de sitio web de WWAP: http://webworld.unesco.org/water/wwap/
www.freemeteo.com. (s.f.). Freemeteo. Obtenido de www.freemeteo.com
Zheng C. (1989). PATH3D. Rockville, Maryland.
149
10 Anexos
150
Modelación de Flujo Subterráneo. Aplicación en la cuenca del Milagro-Ecuador y análisis de pozo en subcuenca del Río Suquía
Serie 33 años
Registro de estación Meteorológica Milagro 1981-2014
Año Enero Febrero Marzo Abril Mayo Junio Julio Agosto Septiembre Octubre Noviembre Diciembre Suma promedio
1981 130.00 370.00 316.00 150.00 6.00 2.00 4.00 0.00 0.00 0.00 0.00 0.00 978.00 81.50
1982 140.00 130.00 40.00 30.00 15.00 12.00 7.00 0.00 0.00 51.00 0.00 0.00 425.00 35.42
1983 420.00 410.00 660.00 790.00 425.00 435.00 0.00 0.00 0.00 6.00 0.00 0.00 3146.00 262.17
1984 620.00 305.00 270.00 10.00 12.00 0.00 0.00 0.00 0.00 0.00 0.00 1217.00 110.64
1985 130.00 60.00 210.00 40.00 6.00 4.00 0.00 0.00 0.00 0.00 0.00 0.00 450.00 37.50
1986 550.00 220.00 70.00 160.00 6.00 12.00 0.00 0.00 0.00 0.00 0.00 0.00 1018.00 84.83
1987 490.00 700.00 525.00 390.00 100.00 105.00 0.00 0.00 0.00 4.00 0.00 0.00 2314.00 192.83
1988 420.00 270.00 30.00 260.00 120.00 125.00 0.00 0.00 0.00 0.00 0.00 0.00 1225.00 102.08
1989 510.00 470.00 305.00 230.00 35.00 30.00 0.00 0.00 0.00 0.00 0.00 0.00 1580.00 131.67
1990 56.30 302.70 177.60 110.80 28.40 43.80 0.20 0.10 0.40 4.50 0.50 67.90 793.20 66.10
1991 60.50 448.40 222.40 10.80 25.70 1.40 1.40 0.40 0.90 0.60 1.70 58.80 833.00 69.42
1992 127.30 628.70 687.50 506.40 382.40 11.80 5.90 0.30 0.10 0.70 3.70 4.20 2359.00 196.58
1993 237.60 668.50 389.50 275.20 59.00 0.50 1.10 0.00 0.00 2.00 0.40 147.00 1780.80 148.40
1994 255.00 318.10 264.20 198.90 43.90 0.20 0.00 0.00 0.20 0.60 4.40 139.00 1224.50 102.04
1995 308.10 203.20 159.70 108.00 15.60 0.00 3.90 0.10 0.60 0.20 2.10 17.60 819.10 68.26
1996 79.50 355.70 338.60 42.00 0.80 0.10 0.50 0.00 0.10 0.00 1.00 8.60 826.90 68.91
1997 208.60 487.20 680.00 381.50 99.50 200.50 156.30 31.70 127.10 30.50 551.20 595.50 3549.60 295.80
1998 712.70 846.00 966.40 957.30 517.80 91.40 33.80 0.00 0.50 0.40 0.90 3.20 4130.40 344.20
1999 113.50 519.80 0.00 207.50 60.30 2.10 0.00 0.00 2.90 3.50 3.80 82.00 995.40 82.95
2000 110.10 221.40 284.90 164.10 134.10 5.90 0.00 0.00 3.50 2.00 0.10 4.00 930.10 77.51
2001 376.00 386.60 748.40 179.50 16.70 0.00 0.50 0.00 0.00 0.10 0.00 0.70 1708.50 142.38
2002 57.10 590.20 594.00 278.80 6.50 5.80 0.00 0.00 0.00 5.30 0.00 55.90 1593.60 132.80
2003 115.20 337.00 173.10 90.70 51.90 0.80 0.40 0.00 0.00 1.30 0.00 28.20 798.60 66.55
2004 101.70 275.40 332.70 145.60 66.00 2.30 1.00 0.00 5.90 1.50 0.40 20.60 953.10 79.43
2005 131.00 201.50 263.80 120.00 25.00 0.00 0.00 0.10 0.00 0.00 1.10 44.00 786.50 65.54
2006 316.90 573.80 339.00 20.00 56.00 2.10 0.00 1.40 0.00 0.10 2.70 23.20 1335.20 111.27
2007 211.10 97.50 481.60 162.00 15.90 3.00 1.20 0.00 0.00 2.90 5.10 3.00 983.30 81.94
2008 468.00 669.90 619.20 272.90 48.40 1.00 0.20 0.60 2.90 2.50 0.40 0.10 2086.10 173.84
2009 456.90 380.80 187.60 49.10 66.10 2.40 0.30 0.20 0.00 0.10 0.00 18.00 1161.50 96.79
2010 163.50 464.90 238.40 330.50 34.80 1.30 3.30 0.70 0.00 0.20 7.10 124.80 1369.50 114.13
2011 158.10 263.20 39.20 393.80 0.60 18.30 31.50 0.10 0.00 0.00 0.50 6.80 912.10 76.01
2012 457.50 806.20 494.10 223.10 85.00 0.40 0.00 0.00 0.00 4.90 0.40 6.10 2077.70 173.14
Ilustración 10-1 Serie de tiempo pluviometría 1981-2014 (INAMHI, 2013)
151
Estacion Precipitacion(m Temperatura(° Evaporacion(m Evapotranspiracion(m

meteorologica Latitud Longitud m) C) m) m)
Milagro 02°06´56´´ 79°35´57´´ 1526.10 25.8 1331.2 735
San Carlos 02°12´00´´ 79°26´35´´ 1843.50 27.3 1000.6 764
Guayaquil 02°16´01´´ 79°54´03´´ 1242.40 26.7 1703.5 739
Babahoyo 01°49´41´´ 79°33´49´´ 2108.60 26.0 1300.9 764
Ilustración 10-2 Datos de Estaciones Meteorológicas cercanas al lugar de estudio.

(Instituto Nacional de Meteorología e Hidrología(INAMHI), 2009)
Ilustración 10-3 (INAMHI, 2013)
152
Ilustración 10-4 Descripción de Mapa Izopiezas (INAMHI, 2013)
Ilustración 10-5 Ubicación de la Cuenca Hidrográfica del Guayas (INAMHI, 2008)
153
Ilustración 10-6 Precipitación plurianual SERIE 1965-1999 (INAMHI, 2008)
Ilustración 10-7 Detalle según colores de los rangos de precipitaciones (INAMHI, 2008)
154
NOMENCLATURA
OBRA: CONSTRUCCION, BLOQUE DE AULAS, UBICACIÓN: CANTON MILAGRO, PROVINCIA DEL GUAYAS.
W HUMEDAD NATURAL qu COMPRESION SIMPLE PERFORACIÓN N°: 1
UBICADA EN LA UNIVERSIDAD DE MILAGRO. Wl LIMITE LIQUIDO ε DEFORMACION UNIT.
Wp LIMITE PLASTICO N NUMERO DE GOLPES VERIFICADO POR: ING. LUCRECIA MORENO
IP INDICE PLASTICO Su.tor Su.torvane FECHA INICIO: 24 DE AGOSTO DE 2013 COORDENADAS:
g DENSIDAD HÚMEDA S-3 TUBO SHELBY NIVEL FREATICO: SI NORTE: 9762102
CONSULTOR: BSLC CONSULTORES CIA. RELLENO NO ESTE: 655377
LTDA.
0 50 100 150 200 250 300
0
MUESTRA W Wl Wp IP g qu Su.tor.
ESTRATI e # 4 # 40 # 200
DESCRIPCION GRAFIA
SUCS SPT
kg/cm3 Ton/m2 % T/m2 % % % 1
# PROF. % % % %
2
1 1.50 MH 33.05 71,30 38,70 32,60 100.0 98.38 95.87 (4-8-9)
LIMO DE ALTA PLASTICIDAD COLOR GRIS OSCURO 3

DE CONSISTENCIA MUY COMPACTA 2 2.50 MH 39.56 53,30 29,56 23,74 100.0 100.00 98.20 (9-9-10)
4
3 3.00 MH 25.98 58,70 33,36 25,34 100.0 97.39 88.89 (6-7-10)
5
LIMO DE BAJA PLASTICIDAD GRIS VERDOSO 4 4.50 ML 34.77 39,64 27,75 11,89 11.97 6.5 99.6 97.62 82.80 SHELBY
LIMO GRIS OSCURO DE CONSISTENCIA MUY

BLANDA 5 5.00 ML 38.93 36,76 26,43 10,33 100.0 98.92 66.20 (1-1-1) 6
ARCILLA GRIS OSCURO DE CONSISTENCIA MUY

BLANDA 6 6.00 CL 30.77 39,52 24,12 15,40 2.14 100.0 98.52 61.95 SHELBY 7
TURBA 7 7.00 PT 293.15 162,15 111,39 50,76 100.0 56.80 15.49 SHELBY 8
LIMO ARCILLOSO CON INTERCALACIONES DE

TURBA
8 8.00 MH 73.72 179,28 120,88 58,40 3.88 100.0 99.19 91.47 SHELBY 9
ARCILLA DE ALTA PLASTICIDAD GRIS VERDOSA DE

CONSISTENCIA BLANDA 9 9.50 CH 66.32 65,72 29,19 36,53 8.05 4.5 100.0 96.63 86.97 SHELBY 10
LIMO DE BAJA PLASTICIDAD COLOR GRIS VERDOSO

DE CONSISTENCIA MUY BLANDA
10 10.00 MH 40.33 71,59 41,04 30,55 22.18 100.0 96.68 90.22 SHELBY 11
LIMO DE ALTA PLASTICIDAD COLOR GRIS VERDOSA

CON ALGO DE YESO
11 11.00 MH 33.00 67,29 37,72 29,57 98.1 93.77 90.05 (7-12-20) 12
LIMO DE BAJA PLASTICIDAD DE COLOR CAFÉ CLARO

CON ARCILLA GRIS VERDOSA Y ALGO DE YESO
12 12.00 ML 37.27 0.00 0.00 0.00 100.0 95.84 84.58 (11-11-13) 13
LIMO DE ALTA PLASTICIDAD CAFÉ CLARO CON 13 12.50 MH #¿NOMBRE? 0.00 0.00 0.00 100.0 0.00 0.00 (11-8-10) 14
INTERCALACIONES GRIS VERDOSO DE
CONSISTENCIA MUY COMPACTA 14 14.00 MH 64.49 69,73 55,13 14,60 100.0 99.16 84.86 (13-13-16) 15
LIMO DE BAJA PLASTICIDAD CAFÉ CLARO CON LIMO

GRIS CLARO DE CONSISTENCIA MUY COMPACTA
15 15.00 ML 36.54 43,72 32,06 11,66 100.0 99.63 77.49 (11-11-17) 16
LIMO DE ALTA PLASTICIDAD COLOR CAFÉ CLARO

DE CONSISTENCIA MUY COMPACTA
16 16.00 MH 40.50 50,43 37,89 12,54 100.0 99.91 98.58 (10-12-18) 17
LIMO DE BAJA PLASTICIDAD COLOR CAFÉ OSCURO

DE CONSISTENCIA DURA 17 17.00 ML 38.90 48,58 35,36 13,22 100.0 99.11 93.24 (7-12-20) 18
LIMO DE BAJA PLASTICIDAD COLOR GRIS CLARO

DE CONSISTENCIA DURA
18 18.50 ML 32.61 49,31 31,69 6,23 100.0 99.60 95.28 (8-14-21) 19
19 19.00 ML 31.10 44,38 31,69 12,69 100.0 99.10 92.65 (14-20-25) 20
20 20.00 ML 30.4 33,72 27,49 6,23 100.0 93.7 76.3 (20-22-27) 21

21 21.00 ML 29.7 32,44 26,32 6,12 100.0 94.7 78.5 (25-28-32) 22
22 22.00 ML 28.3 29,82 24,25 5,57 100.0 95.8 77.1 (27-30-33) 23
SPT H WL IP
23 23.00 ML 34.3 37,48 32,24 5,24 100.0 92.3 73.0 (26-31-35)
FIN DE LA PERFORACION
Ilustración 10-8 Perforación estratigráfica próxima a los pozos de estudio. (B.S.L.C Consultores Cía. LTda., 2013)
155
156

Modflow PDF

Cargado por

Copyright:

Formatos disponibles

Modflow PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Modflow PDF

Cargado por

Copyright:

Formatos disponibles

Modelación de Flujo Subterráneo.

Aplicación en la cuenca del Milagro-Ecuador y análisis de pozo en

Facultad de Ciencias Exactas Físicas y Naturales

Universidad Nacional de Córdoba

Maestría en Ciencias de la Ingeniería – Mención en Recursos Hídricos

Modelacion de Flujo Subterraneo.

Ing. Omar Yagual Muñoz

Dra. Ing. Teresa María Reyna

Modelacion de Flujo Subterraneo.

1 Introducción al Problema ............................................................................................................ 11

Ilustración 1-1 Cuenca del Milagro y área de estudio Unemi .................................................... 13

Tabla 4-1 Parámetros climáticos promedio de Córdoba. (SMN, 2013)

Modelacion de Flujo Subterraneo.

En el presente trabajo, se seleccionó una subcuenca donde se modelaron

Y por otra parte se analizaron las fluctuaciones y ventajas del monitoreo de

Los pozos monitoreados en Ecuador, demostraron una gran capacidad de

En Córdoba, el pozo que se analizó, demuestra datos de una menor capacidad de

Palabras Claves: agua subterránea, monitoreos, modelos matemáticos de

The Knowledge of the groundwater flow is an essential part of studies of

In addition it was considered a modeling contaminant transport, which describes in

1.1 Introducción y Planteamiento del Problema

Según datos publicados por el Programa Mundial de Evaluación de los Recursos

El uso agrícola es el que demanda la mayor cantidad de agua. Cuantificar que

La cuenca del Suquía, está ubicada en la zona Mediterránea Argentina mientras

Por otra parte la mayoría de los estudios de perforaciones se realizan a través de

1.2 Cuenca del Río Milagro

Ilustración 1-1 Cuenca del Milagro y área de estudio Unemi

Se planteó otro escenario de simulación, en el cual partiendo de las condiciones

Otro escenario de simulación que se consideró, es el de una pluma de

1.3 Subcuenca del Río Suquía

importancia de la República Argentina después de su capital Buenos Aires, con un

A continuación se plantean los objetivos generales y particulares de la presente

1.4.1 Objetivo general

El objetivo general es:

- Conocer características del acuífero en la zona puntual de estudio, así como

1.4.2 Objetivos Específicos

Como objetivos específicos se contemplan:

- El trabajo en la Cuenca del río Milagro permitió estudiar el comportamiento

- Establecer la geometría del área de estudio de ambas cuencas.

La metodología para obtener los objetivos propuestos, será la que se presenta a

1. Análisis y comprensión de los procesos matemáticos generales en el flujo de

5. Tomar las mediciones de calidad de agua con el equipo respectivo, tales

1.- Introducción al Problema.

2.- Base Conceptual.

3.- Modelos Matemáticos.

4.- Caracterización Regional de la Cuenca del Río Suquía y Milagro.

6.- Materiales y Métodos Sub Cuenca del Río Milagro.

7.- Análisis de los Resultados.

8.- Conclusiones y Recomendaciones.

Finalmente en el capítulo 9 se encuentra la bibliografía y en el 10 los anexos.

En este capítulo se presentan los aspectos conceptuales vinculados al desarrollo

2.1 Ciclo del Agua.

1400.0 Agua Disponible a nivel mundial

Ilustración 2-1 Tipos de Agua disponible en el Mundo. (de Marsily, 1986)

En la Ilustración 2-1, se aprecia la cantidad de agua salada a nivel mundial, en

De lo antes mencionado, se nota que el agua salada es una gran mayoría,

 La energía térmica, que es generada principalmente por la radiación solar,

 La gravedad, que supone el desplazamiento de las masas de agua desde

2.2 Parámetros hidrogeológicos