facultad de ciencias exactas y naturales universidad ... · factores como la topografía,...
TRANSCRIPT
FACULTAD DE CIENCIAS EXACTAS y NATURALES
UNIVERSIDAD NACIONAL DE LA PAMPA
TESINA PRESENTADA PARA OBTENER
EL GRADO ACADEMICO DE
INGENIERO EN RECURSOS NATURALES Y MEDIO AMBIENTE
“Estudio de la recarga puntual en las depresiones del centro-este de la
provincia de La Pampa”
Autor: Roberto Alejandro POCHETTI
Director: Dr. Pablo Fernando DORNES
SANTA ROSA (LA PAMPA)
ARGENTINA
2016
2 3 de junio de 2016
PREFACIO
Esta tesis es presentada como parte de los requisitos para alcanzar el grado académico de
Ingeniero en Recursos Naturales y Medio Ambiente, de la Universidad Nacional de La
Pampa y no ha sido presentada previamente para la obtención de otro título en esta
Universidad ni en otra Institución Académica. Se llevó a cabo en la cátedra de
“Climatología e Hidrología General” dependiente del Departamento de Recursos
Naturales, durante el periodo comprendido entre el 01 de Septiembre de 2014 al 15 de
Abril de 2016. Bajo la dirección del Dr. Pablo Fernando Dornes.
3 3 de junio de 2016
AGRADECIMIENTOS:
Agradezco a la Universidad Nacional de La Pampa y la nación Argentina por
brindarme educación pública y de calidad durante todos estos años.
A mis padres, mi tía Mónica, a mi hermana Mayra y a mi novia Yesica por toda su
ayuda y contención en todo momento.
A mi director Pablo Dornes por sus explicaciones y tiempo dedicado durante la
realización de esta tesis.
A Roció Comas, Pablo Diaz y Ramiro Páez por su colaboración activa durante las
salidas a campo sin importar las condiciones climáticas. Además de sus consejos y
recomendaciones.
A Claudia Montoya por su colaboración en el procesamiento de muestras de suelo.
A la Facultad de Ciencias de Exactas y Naturales por el espacio brindado en para
trabajar en la tesis.
A la Facultad de agronomía y la subsecretaria de recursos hídricos por el apoyo
recibido.
Y a todas las personas que me rodean, las cuales siempre me ayudaron a resolver
inconvenientes que surgieron durante la tesina.
4 3 de junio de 2016
INDICE GENERAL
Lista de Figuras 5
Lista de Tablas 6
RESUMEN
1. INTRODUCCIÓN 8
2. HIPÓTESIS Y OBJETIVOS 9
2.1 Hipótesis 9
2.2 Objetivo General 10
2.3 Objetivos Específicos 10
3. MATERIALES Y MÉTODOS 10
3.1 Caracterización del área de estudio 10
3.1.1 Ubicación 10
3.1.2 Clima 11
3.1.3 Geomorfología 11
3.1.4 Suelos 12
3.1.5 Hidrología 13
3.2 Sitio de Estudio 14
3.3 Metodología 15
3.3.1 Caracterización topográfica 15
3.3.2 Caracterización hidroquímica 17
3.3.3 Caracterización edáfica 18
3.3.4 Determinación de superficie piezométrica 19
3.3.5 Determinación de uso de suelo y cobertura 21
3.3.6 Información Meteorológica 22
4. RESULTADOS 24
4.1 Caracterización topográfica 24
4.2 Caracterización hidroquímica 32
4.3 Caracterización edáfica 33
4.4 Determinación de superficie piezométrica 36
4.5 Determinación de uso de suelo y cobertura 40
4.6 Información Meteorológica 43
5. DISCUSIÓN 45
6. CONCLUSIONES 56
7. BIBLIOGRAFÍA 57
8. ANEXOS 59
5 3 de junio de 2016
LISTA DE FIGURAS:
Figura 3.1: Ubicación del área de estudio. Fuente imagen satelital: Google Earth (2013).
Figura 3.2: Precipitación media mensual para Winifreda, La Pampa. (1921-2011) en barras
azules. Fuente: Administración Provincial del Agua (APA).
Figura 3.3: Ubicación de las subregiones geomorfologicas en la provincia de La Pampa.
Fuente: Inventario de Recursos Naturales.
Figura 3.4: Posición de los suelos en la subregión de las colinas y lomas. Fuente: Inventario
de Recursos Naturales.
Figura 3.5: Ubicación de la región hídrica de los valles transversales en la provincia de La
Pampa.
Figura 3.6: Ubicación del sitio de estudio. Fuente: SAS Planet.
Figura 3.7: Ubicación del área sujeta a nivelación y puntos acotados.
Figura 3.8: Ubicación de los puntos muestreados en la zona de estudio y cuenca
hidrográfica Pichi Luan. Fuente Google Earth.
Figura 3.9: Muestreo de suelo en la loma, en inmediaciones de la Estación Meteorológica.
Figura 3.10: Puntos muestreados para la determinación piezométrica.
Figura 3.11: Índice de vegetación diferencial normalizado. Fuente: Li et., al 2013.
Figura 3.12: Ecuación de Penman-Monteith (FAO-56). Fuente Organización de las Naciones
Unidas para la Alimentación y la Agricultura.
Figura 3.13: Ubicación de la Estación Meteorológica.
Figura 4.1: Cuenca hidrográfica de la Laguna Pichi Luan y su Modelo Digital de Elevación.
Figura 4.2: Cuenca hidrográfica de la Laguna Pichi Luan y sus Curvas de Nivel equidistancia
5m.
Figura 4.3: Cuenca hidrográfica de la Laguna Pichi Luan y drenajes preferenciales.
Figura 4.4: Ubicación de las transectas niveladas en la parte baja de la cuenca.
Figura 4.5: Perfil topográfico de la Transecta A (SSO-Freatímetro).
Figura 4.6: Perfil topográfico de la Transecta B (SO-Freatímetro).
6 3 de junio de 2016
Figura 4.7: Perfil topográfico Transecta C.
Figura 4.8: Perfil topográfico de la Transecta D.
Figura 4.9: Mapa de iso-salinida correspondiente a la cuenca hidrográfica de la Laguna
Pichi Luan.
Figura 4.10: Distribución de partículas en la muestra LPEG1-10cm.
Figura 4.11: Distribución de partículas en la muestra LPEG1-60cm.
Figura 4.12: Distribución de partículas en la muestra LPII Lag-P.Luan-10cm.
Figura 4.13: Distribución de partículas en la muestra LPII Lag-P.Luan-30cm.
Figura 4.14: Mapa Piezométrico de la cuenca hidrográfica de la laguna Pichi Luan.
Figura 4.15: Mapa Piezométrico y flujo de agua subterránea de la cuenca hidrográfica de la
laguna Pichi Luan.
Figura 4.16: Nivel piezométrico en el área de la laguna Pichi Luan, fecha: 7 de octubre del
2014.
Figura 4.17: Nivel piezométrico en el área de la laguna Pichi Luan, Fecha 29 de octubre del
2014.
Figura 4.18: Evolución de la cobertura vegetal año 2014.
Figura 4.19: Precipitación estación meteorológica de El Guanaco, 2014.
Figura 4.20: Evapotranspiración calculada para el año 2014.
Figura 5.1: Evolución de la cobertura vegetal la estación de verano año 2014.
Figura 5.2: Evolución de la cobertura vegetal la estación de otoño año 2014.
Figura 5.3: Evolución de la cobertura vegetal la estación de invierno año 2014.
Figura 5.4: Evolución de la cobertura vegetal la estación de primavera año 2014.
Figura 5.5: Comparación topografía y nivel Piezométrico de la Transecta A, posterior a la
precipitación del 2-5 de Octubre, iniciando en el freatímetro extendiéndose al S.S.O. véase
figura 4.4.
7 3 de junio de 2016
Figura 5.6: Comparación topografía y nivel piezométrico de la Transecta B, posterior a la
precipitación del 2-5 de Octubre, iniciando en el freatímetro extendiéndose al S.O. véase
figura 4.4.
Figura 5.7: Comparación topografía y nivel piezométrico de la Transecta A en su sección
más próxima al freatímetro, posterior a la precipitación del 27 de Octubre, iniciando en el
freatímetro extendiéndose al S.O. Véase figura 4.4.
Figura 5.8: Comparación topografía y nivel piezométrico de la Transecta C, posterior a la
precipitación del 27 de Octubre, iniciando en el freatímetro extendiéndose al Norte. Véase
figura 4.4.
Figura 5.9: Comparación topografía y nivel piezométrico de la Transecta D, posterior a la
precipitación del 27 de Octubre, iniciando en la intersección con la transecta C y
terminando al Este de la mima. Véase figura 4.4.
LISTA DE TABLAS
Tabla 3.1: Comparación de los sensores del Landsat 8 OLI/TIRS y el sensor del Landsat 7
ETM+. Fuente: Li et., al 2013.
Tabla 4.1: Comparación de las propiedades físicas de muestras extraídas en la Loma.
Tabla 4.2: Comparación de las propiedades físicas de muestras extraídas en el bajo.
Tabla 4.3: Extensión de las diferentes coberturas en el año 2014.
Tabla 4.4: Balance hídrico del año 2014.
Tabla 5.1: Detalle de eventos de precipitación de Verano.
Tabla 5.2: Detalle de eventos de precipitación en otoño año 2014.
Tabla 5.3: Detalle de eventos de precipitación en invierno año 2014.
Tabla 5.4: Detalle de eventos de precipitación en primavera año 2014.
8 3 de junio de 2016
1. INTRODUCCIÓN
La región oriental de la provincia de La Pampa se caracteriza por poseer un régimen
hídrico subhúmedo seco. Geomorfológicamente la zona de estudio se sitúa en la
subregión de colinas y lomas, caracterizada por zonas deprimidas alargadas donde se
acumula el excedente pluvial. Desde un punto de vista hidrológico se encuentra en una
zona de transición entre la región de la planicie pampeana central y los valles
transversales, en el área no existen patrones de drenajes definidos como arroyos y ríos
(Salazar Lea Plaza, 1980). Los cuerpos de agua pueden ser agrupados hidrográficamente
en el grupo de lagunas de El Guanaco o Luan Lauquén, caracterizado por depresiones con
orientación SO-NE que forman pequeños valles separados por lomas de origen eólico con
la misma orientación las cuales generan una desconexión superficial propiciando la
generación de cuencas endorreicas que presentan acumulación de excedentes de agua en
épocas ricas en precipitaciones (Cazenave Et. al., 1987).Los cuerpos de agua son en
general de carácter temporal. El balance hídrico del área de estudio es negativo durante la
mayor parte del año. Sin embargo existen períodos, principalmente en las épocas de
transición, donde las precipitaciones exceden a la evaporación, dando lugar a un balance
positivo.
El comportamiento de la dinámica de agua en laderas es complejo y difícil de caracterizar,
debido a que los procesos de generación de escorrentía y almacenamiento de agua en el
suelo están fuertemente influenciados por la heterogeneidad espacial que presentan las
áreas de aporte. La distribución de la humedad del suelo depende principalmente de
factores como la topografía, propiedades físicas del suelo, y la vegetación o cobertura. La
topografía a su vez ejerce un control primario sobre la partición de volúmenes de
infiltración y escorrentía. La textura y la estructura del suelo, afectan a la capacidad de
almacenar y trasmitir agua, en tanto que la vegetación, cobertura y/o manejo afecta a la
conductividad hidráulica, intercepción y generación de encostramientos. (Famiglietti et.
al., 1998). En estos ambientes semiáridos las respuestas suelen ser no lineales, siendo
determinante en los procesos hidrológicos la condición previa del sistema, por ejemplo la
humedad antecedente. Otros controles secundarios en la generación de escorrentía, son
la duración de la tormenta y la intensidad de la precipitación (Beven, 2002; Castillo et. al.,
2003).
Antecedentes en el área de estudio mostraron que existen fenómenos de escorrentía
superficial y que los mismos eran respuesta a factores climáticos, de relieve, de manejo y
a las condiciones antecedentes. La combinación o presencia de dichos factores resultó en
la activación de distintos procesos de generación de escorrentía (Secco et. al., 2013).
9 3 de junio de 2016
La mayoría de descripciones hidrológicas en áreas de llanura con nivel freático somero,
caracterizan a los bajos por su comportamiento efluente es decir que constituyen zonas
de descarga del agua freática, e identifica como zonas de recarga, a las áreas
topográficamente elevadas o intermedias (Auge, 2008). Las lagunas o humedales que
forman parte del paisaje de llanura han sido identificados como sitios ecológicos
importantes y dependen en gran manera del aporte de agua subterránea (Custodio,
2010). Aunque estudios en ambientes semiáridos de praderas de Norte América
demostraron que los humedales desempeñan un rol importante en la recarga del agua
subterránea, siendo más evidente su efectos a nivel local, pero podían llegar a influir en
los sistemas hidrológicos regionales (Van der Kamp, Hayashi, 1998).Así mismo, se
comprobó que debajo de humedales y lagunas temporales existe una lixiviación de sales
debido a fenómenos de recarga puntual en depresiones, indistintamente de su tamaño
(Berthold et. al., 2004).
En los ecosistemas semiáridos de la llanura pampeana si bien las precipitaciones son la
principal fuente de ingreso de agua a los sistemas hidrológicos (Secco et. al., 2012), la
presencia de niveles freáticos profundos pone en evidencia la existencia de otros procesos
de recarga. Existen diferentes mecanismos que explican el proceso de recarga. La misma
puede darse en un proceso espacial diferido y difuso a través de la zona no saturada o
deberse a un proceso instantáneo o localizado donde los flujos preferenciales permitan al
agua circular rápidamente por el perfil (Lerner., 1990). La recarga por lluvia es la situación
más común y regionalmente más importante, sin embargo la acumulación de agua en
depresiones debido a la concentración de escorrentía superficial, puede constituir una
fuente de recarga adicional y puntual. De esta manera, la existencia de flujos superficial en
laderas o pendientes, son los que producen la redistribución local de la recarga y la
acumulación de agua en las áreas deprimidas del relieve (Custodio., 1997).
En este contexto, se caracterizó el rol hidrológico y su relación con la recarga, de bajos o
lagunas con concentración de flujos de escorrentía.
2. HIPOTESIS Y OBJETIVOS:
2.1 Hipótesis
Típicamente se describe la dinámica de las lagunas o depresiones de la región pampeana
en función de la descarga acuífera, lo que determina que su ocurrencia está
mayoritariamente vinculada a fluctuaciones de los niveles saturados y sólo indirectamente
a eventos pluviométricos. Por otro lado, el proceso de recarga se asocia a geoformas
eólicas que facilitan el proceso de infiltración. Sin embargo, tales conceptualizaciones en
la región pampeana occidental con niveles saturados profundos no explican la ocurrencia
de ciertos cuerpos de agua o lagunas.
10 3 de junio de 2016
La eventual presencia de estos sistemas hidrológicos desconectados del nivel freático
estaría relacionada a procesos tales como la escorrentía directa y la redistribución local de
la precipitación, ambos procesos son favorecidos por la presencia de geoformas que
actuarían como cuencas endorreicas de acumulación de excedente pluvial. Estos procesos
darían lugar a la eventual recarga puntual del acuífero libre ubicado debajo de estos bajos
o lagunas temporales.
2.2 Objetivo General
El objetivo de este estudio es identificar el proceso de recarga a partir de los volúmenes
de agua almacenados en depresiones del centro-este de la provincia de La Pampa.
2.3 Objetivos Específicos
En particular se propone:
a) Evaluar el carácter intermitente de lagunas o depresiones a partir del aporte de la
escorrentía superficial y del análisis de características climáticas, topográficas, edáficas y
de manejo o cobertura.
b) Vincular la dinámica hidrológica superficial con las fluctuaciones del nivel freático y su
relación con el ciclo hidrológico.
3. MATERIALES Y MÉTODOS
3.1 Caracterización del área de estudio
3.1.1 Ubicación
El trabajo se llevó a cabo en la cuenca de la laguna Pichi Luan ubicada a 35 km al norte de
la ciudad de Santa Rosa, en las coordenadas (36º 17' 59.9" S, 64º 17' 1.6" O). Se sitúa en el
departamento Capital, lote 2, fracción A, sección XXII. El área de estudio se encuentra
inmersa en una matriz agrícola-ganadera. (Figura 3.1)
Figura 3.1: Ubicación del área de estudio. Fuente imagen satelital: Google Earth (2013).
11 3 de junio de 2016
3.1.2 Clima
El área de estudio se encuentra en el régimen de humedad subhumedo seco y se
encuentra entre la isohietas de 600 y 700 mm (Casagrande y Conti, 1980). Las
precipitaciones son características de los meses de transición entre el verano y el invierno
debido a influencia de los anticiclones ubicados en el Océano Pacífico y el Océano
Atlántico. La estación estival presenta eventos de precipitaciones generados por
condiciones ciclónicas o por el avance de frentes fríos sobre masas de aire caliente, estos
eventos poseen la particularidad de ser muy intensos y de corta duración. Así también
presenta una estación invernal típicamente seca donde de la presencia de eventos
pluviométricos está caracterizada por la baja intensidad y larga duración. Los patrones
pluviométricos poseen mucha variabilidad espacial, como así también interanual (Gatto
Cáceres y Dornes, 1996). (Figura 3.2).
Figura 3.2: Precipitación media mensual para Winifreda, La Pampa. (1921-2011). Fuente:
Administración Provincial del Agua (APA).
La zona de estudio se encuentra dentro de las isotermas medias anuales de 15.5 °C y 16
°C, con temperaturas máximas en verano cercano a 24 ° C y mínimas en julio entre las
isotermas de 7 y 8 °C. Los vientos predominantes son del S-SO y N-NE, con una velocidad
media de 10-12 km/h (Casagrande y Conti, 1980).
3.1.3 Geomorfología
Geomorfológicamente la zona de estudio se sitúa entre la subregión de las planicies con
tosca y la subregión de colinas y lomas (Figura 3.3). El relieve primitivo de la región fue el
de una pediplanicie, modificada en el tiempo por procesos morfogenéticos de origen
0
20
40
60
80
100
120
Pre
cip
itac
ión
(m
m)
Tiempo (meses)
12 3 de junio de 2016
hídrico y eólicos. En la litología de la región se distinguen esencialmente dos tipos de
materiales. El primero de ellos se dispone superficialmente y posee una textura arenosa a
franco arenosa, correspondiente a sedimentos depositados por acción del viento en
orientación SE-NO. Debajo de dichos materiales se disponen los sedimentos
correspondiente a la Formación Cerro Azul, estos poseen una coloración rojiza a pardo
rojiza. La constitución de los depósitos está dada por una sucesión materiales loéssicos
aunque es necesario destacar la presencia niveles arenosos. (Salazar Lea Plaza, 1980)
Figura 3.3: Ubicación de las subregiones geomorfologicas en la provincia de La Pampa.
Fuente: Inventario de Recursos Naturales
3.1.4 Suelos
El área de estudio posee suelos de texturas que van desde franco limoso a franco arenosa.
(Peña Zubiate et. al., 1980) (Figura 3.4). Los suelos descriptos para la subregión de colinas
y lomas son los siguientes:
-Suelos dispuestos en lomas: Son haplustoles énticos, familia franco gruesa, mixta,
térmica. Tiene poca evolución edafogenética, presenta un perfil sencillo del tipo A-AC-C; el
Horizonte superficial es profundo y tiene buena provisión de materia orgánica.
-Suelos dispuesto en las pendientes: Son Ustipsamente típico, familia silícea (calcárea),
térmica. Son suelos muy arenosos, moderados a fuertemente inclinados, pobres en
materia orgánica, secos, susceptible a la erosión hídrica y eólica.
13 3 de junio de 2016
-Suelos dispuestos en los bajos: Son haplustoles énticos, familia franco fina, mixta,
térmica. Tienen características muy similares a la de las lomas, pero más húmedos, con
mayor cantidad de materia orgánica y textura más fina (franco limoso)
En el inventario de suelos del país realizado por el INTA, se caracteriza los suelos como no
sódicos, además se reconoce en la zona de trabajo dos limitantes para el uso agropecuario
la profundidad efectiva y el riesgo de erosión hídrica (GeoINTA). Esta última adquiere
mayor importancia en zonas de pendientes, presentes en el área estudiada.
Figura 3.4: Posición de los suelos en la subregión de las colinas y lomas. Fuente: Inventario
de Recursos Naturales.
3.1.5 Hidrología
El área de estudio se ubica en la región hídrica de los valles transversales, la misma
presenta un relieve de orientación SO-NE, conformado por valles, depresiones, cordones
medanosos y otras geoformas de origen eólico (Giraut et. al., 2006). El cuerpo de agua en
estudio, es considerado parte del sistema hidrográfico de grupos de lagunas del Guanaco.
(Cazenave., Et. al. 1987). Este sistema posee las misma orientaciones predominantes en
los valles transversales, las lagunas se ubican en pequeños valles separados por lomas de
origen eólico presentando una desconexión superficial. La profundidad de los cuerpos de
agua es somero presentado en promedio de 1 metro. La región presenta niveles saturados
profundos que rondan los 10 metros de profundidad, siendo variable dependiendo de la
14 3 de junio de 2016
ubicación en el relieve. No presenta redes de drenaje superficiales como arroyos o
manantiales, sin embargo la presencia de estos últimos existe solo en referencias
históricas de la zona (Figura 3.5).
Figura 3.5: Ubicación de la región hídrica de los valles transversales en la provincia de La
Pampa.
3.2 Sitio de estudio
Las actividades se realizaron dentro de la cuenca hidrográfica de la laguna temporaria
Pichi Luan perteneciente al grupo de lagunas denominados Luan Lauquen en el
departamento Capital, provincia de La Pampa. (Figura 3.6)
15 3 de junio de 2016
Figura 3.6: Ubicación del sitio de estudio. Fuente: SAS Planet.
3.3 Metodología
3.3.1 Caracterización Topográfica
A los efectos de poder delimitar el sitio de estudio se procedió a la confección de un
modelo digital de elevación (MDE). Se decidió utilizar las imágenes satelitales de la SRTM
(Shuttle Radar Topography Mission) por su cobertura global y acceso libre. Dentro de los
productos de esta misión se optó por los archivos SRTM Void Filled de resolución espacial
de 90m. Este tipo de imágenes resultaron ser mejores respecto a la representación del
terreno, debido a que para afrontar áreas con faltantes de datos se utilizaron algoritmos
de interpolación dando como resultado imágenes sin celdas con datos nulos. Dichas
imágenes fueron adquiridas mediante la utilización de la herramienta virtual Earth
Explorer del USGS (United States Geological Survey). En base a estas imágenes se elaboró
un modelo raster de celdas igualmente espaciadas con un único valor en el eje (z) o MDE.
La determinación y delimitación de la cuenca hidrográfica de la laguna Pichi Luan se
obtuvo mediante la aplicación de algoritmos de estimación de redes de drenaje, basados
en un modelo determinístico de direccionamiento de flujo D8, el cual parte de un área de
9 celdas en el que se evalúa el comportamiento de la celda central respecto a las
16 3 de junio de 2016
circundantes. Otro algoritmo genera áreas de aporte de dimensiones variables
relacionadas a un punto de salida o drenaje identificado por el modelo D8. Dichas áreas
son polígonos delimitados vectorialmente, los cuales al ser combinados permiten la
creación de cuencas de aporte. Es así que se logra delimitar la cuenca hidrográfica y
obtener las redes de drenajes teóricas.
Con el fin de realizar estudios en detalles en la zona más baja de la cuenca, se realizaron
cuatro perfiles altimétricos mediante el uso de un nivel óptico montado en un trípode y
una regla graduada. Dichos perfiles consistieron en la determinación de cota de diferentes
puntos ubicados en una transecta respecto a un plano de referencia, el datum del mismo
fue definido en base a la altura del freatímetro ubicado en la laguna Pichi Luan a
145msnm. (Figura 3.7).
Figura 3.7: Ubicación del área sujeta a nivelación y puntos acotados.
Las características físicas son determinantes en los procesos hidrológicos y en la
concepción de la magnitud de los mismos, es por eso que se analizaran los siguientes
parámetros morfológicos:
Área: el área de una cuenca es la superficie, delimitada por la línea de divisorias de aguas.
17 3 de junio de 2016
Desnivel altitudinal: es el valor de diferencia entre la cota más alta de la cuenca y la más
baja de la misma. (DA= HM – Hm).
Índice de compacidad: relaciona el perímetro de la cuenca con el perímetro de un círculo
de igual área. Este índice es mayor a medida que la irregularidad de la cuenca aumenta,
así mismo si la cuenca fuera un círculo perfecto el índice sería igual a 1.
𝐾𝑐 = 0.282 𝑃
√𝐴
P= Perímetro A= Área KC= Índice de compacidad
Índice de alargamiento: Este índice relaciona la longitud máxima encontrada en la cuenca,
medida en la dirección principal del drenaje, y el ancho de la misma medida
perpendicularmente. Cuando esta variable toma valores mayores a la unidad se trata de
cuencas alargadas, mientras que valores tendientes a 1 son característicos de cuencas con
drenajes centrales cortos y formas de abanico.
𝐼𝑎 =𝐿𝑚
𝑙
Ia= Índice de alargamiento Lm= Longitud máx. de la cuenca l=ancho máx. de la cuenca.
3.3.2 Caracterización Hidroquímica
Para poder caracterizar la hidroquímica de la zona de estudio, se procedió a realizar un
muestreo de las perforaciones y cuerpos de aguas ubicados dentro y fuera de la cuenca
hidrogeográfica de la laguna Pichi Luan. Los datos fueron obtenidos a campo mediante la
utilización de una sonda multi-paramétrica obteniéndose así valores de PH, Conductividad
y Sales Disueltas Totales. Las campañas fueron realizadas a fines del año 2014 y comienzos
del 2015 relevándose 17 sitios de muestreo. En base a los datos obtenidos se confeccionó
un mapa de iso-conductividad, mediante el uso del método de interpolación Inverse
Distance Weighting. El mismo se basa en la ponderación de los pesos y su relación con la
distancia. Se utilizó el software libre System for Automatet Geoscientific Analyses (SAGA)
(Figura 3.8).
18 3 de junio de 2016
Figura 3.8: Ubicación de los puntos muestreados en la zona de estudio y cuenca
hidrográfica Pichi Luan. Fuente Google Earth.
3.3.3 Caracterización Edáfica
Con el fin de describir la variación de las características del suelo que presenta el área de
estudio, se decidió tomar muestras en la zona alta y baja del relieve de la cuenca de
aporte a Pichi Luan. Las muestras de la loma fueron extraídas a 10cm (LP-EG1-10cm) y a
60cm (LP-EG1-60cm) de profundidad. En el bajo las muestras fueron tomadas a 10cm
(LPII-LAG-Pichi Luan-10cm) y a 30 cm (LPII-LAG-Pichi Luan-30cm). Se determinó para cada
una el porcentaje de carbono orgánico y carbono inorgánico, mediante la utilización de un
horno-mufla SIMCIC HM-3. Para poder conocer la texturas de las muestras se utilizó un
equipo de difracción laser Mastersizer 2000 con una unidad de dispersión Hydro 2000MU
(A) obteniéndose la distribución y el tamaño de las partículas. La determinación de
conductividad se realizó con un conductímetro digital. (Figura 3.9).
19 3 de junio de 2016
Figura 3.9 A: Muestreo de suelo en la loma, en inmediaciones de la Estación
Meteorológica
Figura 3.9 B: Muestreo de suelo en la loma, en inmediaciones de la Estación
Meteorológica.
20 3 de junio de 2016
3.3.4 Determinación de superficie piezométrica
La determinación del flujo de agua subterráneo en la zona se obtuvo mediante la
construcción de un mapa piezométrico. El mismo fue confeccionado en base a la
profundidad de los niveles freáticos obtenidos en catorce perforaciones, ubicadas dentro
y fuera de la cuenca hidrográfica de la laguna Pichi Luan (Figura 3.8). Las mediciones
fueron realizadas mediante la utilización de una sonda piezométrica que permitió la
determinación puntual de la profundidad del nivel de agua. Las alturas piezométricas de
cada perforación fueron obtenidas a partir de los valores del MDE, previamente creado
para caracterización topográfica. Una vez obtenidas dichas alturas se procedió a la
generación de una imagen raster de valores de altura mediante el método de
interpolación Inverse Distance Weighting con el software libre Qgis, en base a estos
valores se generaron las curvas de equipotenciales para la zona de estudio utilizando
geoalgoritmos.
Con el fin de obtener mayor detalle en la zona de la laguna, se midieron los niveles
freáticos de 22 pozos de observación generados mediante la utilización de barrenos
(figura 3.10), en campañas de medición realizadas días posteriores a eventos de
precipitación. La creación de los perfiles piezométricos correspondientes se llevó a cabo
mediante la aplicación de la misma metodología utilizada para productos de nivel
regional; la diferencia radicó en la utilización de alturas obtenidas mediante la nivelación
realizada en la caracterización topográfica.
21 3 de junio de 2016
Figura 3.10: Puntos muestreados para la determinación piezométrica.
3.3.5 Determinación de uso de suelo y cobertura
Para poder determinar el uso del suelo y su cobertura se utilizaron imágenes satelitales
del satélite orbital Landsat 8 (Tabla 3.1). Las imágenes fueron adquiridas mediante la
plataforma virtual Earth Explorer del USGS (United States Geological Survey). Se utilizó un
índice de vegetación diferencial normalizado que permite estimar la calidad, cantidad y
desarrollo de la vegetación NDVI. Se implementó el mismo mediante la fórmula propuesta
por Li et., al (2013) (Figura 3.11).
Figura 3.11: Índice de vegetación diferencial normalizado. Fuente: Li et., al 2013.
22 3 de junio de 2016
Tabla 3.1: Comparación de los sensores del Landsat 8 OLI/TIRS y el sensor del Landsat 7
ETM+. Fuente: Li et., al 2013.
3.3.6 Información Meteorológica
Los datos fueron obtenidos a partir de la estación instalada en la parte superior de la
cuenca, la misma posee registros de humedad de suelo, radiación solar, precipitación,
velocidad del viento, temperatura y presión atmosférica cada 10 minutos (figura 3.13). A
partir de estos datos se realizó un balance de entradas y salidas para poder explicar las
fluctuaciones del nivel freático. Para la evaporación se utilizó el modelo estandarizado de
Penman-Monteith (FAO-56) (figura 3.12).
Figura 3.12: Ecuación de Penman-Monteith (FAO-56). Fuente Organización de las Naciones
Unidas para la Alimentación y la Agricultura.
Donde Rn es la radiación neta en superficie en (Mj/m2 * día), G es el flujo del calor en el
suelo (Mj/m2 * día), es es la presión de vapor de saturación en kPa, ea es la presión real de
vapor en kPa, (es – ea) es el déficit de presión de vapor en kPa, u2 es la velocidad del viento
a 2 metros de altura en (m/s), Δ representa la pendiente de la curva de presión de vapor
(kPa/ºC) y ϒ es la constante psicométrica (kPa/ºC).
Esto permitió definir las perdidas por evaporación en sistema. Las entradas fueron
obtenidas de la precipitación medida en dicha estación automática.
23 3 de junio de 2016
Figura 3.13: Ubicación de la Estación Meteorológica.
24 3 de junio de 2016
4. RESULTADOS
4.1 Topografía
El análisis del MDE permitió la delimitación de la cuenca hidrográfica de Pichi Luan (figura
4.1), la determinación de la curvas de nivel del terreno (Figura 4.2) y redes de drenaje
preferenciales (Figura 4.3). La cuenca posee un perímetro de 31,63 km y un área de
aporte de 1649,5 hectáreas. Mediante la aplicación de los cálculos correspondientes se
pudo determinar que la cuenca presenta un índice de compacidad de 2,21 característico
de cuencas ovaladas alargadas. Así mismo los resultados obtenidos mediante el cálculo
del índice de alargamiento adquiere el valor de 2,77 siendo clasificada como una cuenca
moderadamente alargada. En términos de altura la cuenca presenta un desnivel de 41 m
siendo la cota más alta 184 m. y la cota más baja 143m. La morfología que presenta la
cuenca es cóncava, las redes de drenaje son dendríticas y convergen en la parte más baja
generando una concentración del escurrimiento superficial en las zonas más deprimidas
del terreno.
Figura 4.1: Cuenca hidrográfica de la Laguna Pichi Luan y su Modelo Digital de Elevación.
25 3 de junio de 2016
Figura 4.2: Cuenca hidrográfica de la Laguna Pichi Luan y sus Curvas de Nivel equidistancia
5m.
26 3 de junio de 2016
Figura 4.3: Cuenca hidrográfica de la Laguna Pichi Luan y drenajes preferenciales.
Los puntos nivelados en la parte baja de la cuenca permitieron la realización de cuatro
transectas o perfiles altimétricos (Figura 4.4) (Anexo I), las cuales representan el micro
relieve de la zona nivelada. Mediante el análisis de los perfiles topográficos se demostró
que el nivel de base de la cuenca está en proximidades del freatímetro. Las pendientes
son muy suaves siendo menores al 0,7 % poniendo al descubierto una superficie llana, la
cual se corresponde con el área afectada al anegamiento, producto de la dinámica
presentada por la laguna temporal de Pichi Luan.
27 3 de junio de 2016
Figura 4.4: Ubicación de las transectas niveladas en la parte baja de la cuenca.
28 3 de junio de 2016
Transecta A (Figura 4.5): Se extiende desde el Molino C al SSO de la laguna temporal hasta
el freatímetro ubicado en la parte baja de la misma.
Longitud: 516 metros.
Desnivel altitudinal: 3,59 metros.
Altura máxima: 148,59 MSNM (extremo Sur).
Altura mínima: 145 MSNM (Freatímetro).
Pendiente: 0,69 %.
Figura 4.5: Perfil topográfico en metros de la Transecta A (SSO-Freatímetro).
29 3 de junio de 2016
Transecta B (Figura 4.6): Se extiende desde una zona anegable al SO hasta el freatímetro
ubicado en la zona central de la laguna temporal.
Longitud: 409 metros.
Desnivel altitudinal: 0,82 metros.
Altura máxima: 145,82 MSNM (extremo SO).
Altura mínima: 145 MSNM (Freatímetro)
Pendiente: 0,2 %.
Figura 4.6: Perfil topográfico en metros de la Transecta B (SO-Freatímetro).
30 3 de junio de 2016
Transecta C (Figura 4.7): Se extiende desde el freatímetro hacia el Norte.
Longitud: 413 metros.
Desnivel altitudinal: 0,17 metros.
Altura máxima: 145,17 MSNM (extremo N).
Altura mínima: 145 MSNM (Freatímetro)
Pendiente: 0,04 %.
Figura 4.7: Perfil topográfico en metros de la Transecta C.
31 3 de junio de 2016
Transecta D (Figura 4.8): Se extiende desde una zona anegable a NE del freatímetro hasta
su intercepción con la transecta B 194 metros al Norte del freatímetro.
Longitud: 276 metros.
Desnivel altitudinal: 0,32 metros.
Altura máxima: 145,23 MSNM (72 metros al O del extremo E).
Altura mínima: 144,91 MSNM (Extremo E)
Pendiente: 0,04 %.
Figura 4.8: Perfil topográfico en metros de la Transecta D.
32 3 de junio de 2016
4.2 Hidroquímica
El análisis de los valores de conductividad eléctrica medida en las perforaciones
muestreadas permitió la confección de un mapa de iso-salinidad del acuífero libre (Figura
4.9) (Anexo II). En el mismo pueden observarse la variabilidad de la calidad de agua,
presente en el área de estudio los valores de conductividad eléctrica fluctúan entre los
1,98 mS/m y los 9,98 mS/m. Las zonas que presentan valores inferiores a 5 mS/m se
ubican en el sector central de la cuenca hidrográfica coincidiendo con las zonas
deprimidas del relieve en la mayoría de los casos. A su vez por lo general, las zonas con
mayores valores de conductividad se ubican en las partes elevadas del terreno.
Figura 4.9: Mapa de iso-salinidad correspondiente a la cuenca hidrográfica de la Laguna
Pichi Luan.
33 3 de junio de 2016
4.3 Edafología
Las muestras de suelos analizadas en la loma presentan texturas franco-arenosas, siendo
característico una distribución equilibrada de las fracciones granulométricas con un leve
predominio de las fracciones más gruesas. Poseen un drenaje moderado y buena
aireación, a su vez la presencia de limo y arcillas lo hace un suelo que puede conservar la
humedad necesaria para la radicación de especies vegetales. En términos de salinidad
estos suelos presentan valores bajos siendo clasificados como suelos no salinos. Poseen
bajos contenidos de materia orgánica por ser suelos de textura gruesa.
Tabla 4.1: Comparación de las propiedades físicas de muestras extraídas en la Loma.
Muestra LPEG1-10cm (Fig. 4.10) LPEG1-60cm (Fig. 4.11)
Textura Franco-Arenoso Franco-Arenoso
Arcilla 9,18 % 11,12 %
Limos 26,87 % 35,08 %
Arenas 63,95 % 53,8 %
Conductividad Eléctrica 81 µS/cm 95,9 µS/cm
Materia Orgánica (LOI 105C°) 0,91 % 1,25 %
Carbono Orgánico (LOI 550C°) 3,46 % 3,12 %
Carbono Inorgánico (LOI 950C°) 0,66 % 0,82 %
Figura 4.10: Distribución de partículas en la muestra LPEG1-10cm.
34 3 de junio de 2016
Figura 4.11: Distribución de partículas en la muestra LPEG1-60cm.
Las muestras de suelos analizadas en el bajo presentan suelos de textura Franco-Limosa,
estos poseen mayor contenido de limos, pero con presencia de las demás fracciones
granulométricas en menor proporción. Presentan menor drenaje, menor aireación y
mayor capacidad de contener humedad que los suelos descriptos en la zona de la Loma.
Estos suelos son propensos a la compactación por el laboreo agrícola dado su elevado
contenido de limos, estas prácticas provocan un aumento en la densidad. Posee buenos
contenidos de materia orgánica los cuales favorecen la retención de humedad, y la
radicación de especies vegetales. En términos de salinidad no son considerados suelos
salinos, pero cabe destacar que presentan mayor contenidos de sales que los suelos
muestreados en la loma. Esto podría deberse a procesos que permiten la concentración
de sales tales como la acumulación de excesos de precipitación y su posterior
evaporación. Por lo general los suelos que se ven afectados por el ascenso del nivel
freático presentan valores más elevados de salinidad, generando costras salinas en
superficie, características que los suelos de la laguna temporal de Pichi Luan no presentan.
A su vez las muestras posen un tono grisáceo característicos de condiciones reductoras y
ambientes anaeróbicos, saturados de agua.
35 3 de junio de 2016
Tabla 4.2: Comparación de las propiedades físicas de muestras extraídas en el bajo.
Muestra LPII Lag-P.Luan-10cm (Fig. 4.12)
LPII Lag-P.Luan-30cm (Fig. 4.13)
Textura Franco-Limoso Franco-Limoso
Arcilla 18,1 % 19,95 %
Limos 71,22 % 66,83 %
Arenas 10,67 % 13,22 %
Conductividad Eléctrica 367 µS/cm 148,1 µS/cm
Materia Orgánica (LOI 105C°) 4,82 % 3,4 %
Carbono Orgánico (LOI 550C°) 14,52 % 4,68 %
Carbono Inorgánico (LOI 950C°) 1,95 % 1,73 %
Figura 4.12: Distribución de partículas en la muestra LPII Lag-P.Luan-10cm.
36 3 de junio de 2016
Figura 4.13: Distribución de partículas en la muestra LPII Lag-P.Luan-30cm.
4.4 Superficie Piezométrica
El análisis del mapa de la superficie piezométrica (Figura 4.14) creado mediante la
interpolación de perforaciones de la zona permitió inferir el movimiento del agua
subterránea de la cuenca. Se identificaron los flujos de agua subterránea
correspondientes (Figura 4.15), los mismos son de carácter concéntrico y fluyen hacia la
zona de la laguna temporaria. También se identificó otra zona deprimida al oeste de la
cuenca hacia donde se dirigiría parte del flujo del agua subterránea, existiendo en el
centro de la cuenca una divisoria de aguas subterráneas.
37 3 de junio de 2016
Figura 4.14: Mapa Piezométrico de la cuenca hidrográfica de la laguna Pichi Luan.
38 3 de junio de 2016
Figura 4.15: Mapa Piezométrico y flujo de agua subterránea de la cuenca hidrográfica de la
laguna Pichi Luan.
En lo que respecta a la piezometría de la zona de la laguna temporaria, la misma presenta
mayor dinámica. Asociada a la ocurrencia de eventos pluviométricos y a su intensidad.
Para evaluar esta variación puntual de los niveles piezómetros se analizaron en dos fechas
posteriores a eventos de precipitación el 7 de octubre del 2014 (Figura 4.16) (Anexo III A5)
y el 29 de octubre del 2014 (figura 4.17) (Anexo III A6).
39 3 de junio de 2016
Figura 4.16: Nivel piezométrico en el área de la laguna Pichi Luan, fecha: 7 de octubre del
2014.
Figura 4.17: Nivel piezométrico en el área de la laguna Pichi Luan, Fecha 29 de octubre del
2014.
40 3 de junio de 2016
4.5 Uso del suelo y cobertura vegetal
La interpretación de los índices de vegetación aplicados a las imágenes satelitales Landsat
8, permitieron la creación categorías que representan la cobertura en vegetal a lo largo
del año. Es así que para cada imagen se crearon tres clases de cobertura vegetal:
Clase A: Valores del índice NDVI comprendidos entre 0 y 0,2 estos valores se
corresponden con suelo desnudo o baja cobertura de vegetación. Dicha clase en la
representación gráfica adopta el color rojo.
Clase B: Valores del índice NDVI comprendidos entre los valores 0,2 y 0,4 estos valores se
corresponden con vegetación en crecimiento o cobertura de vegetación media. Esta clase
en la representación gráfica adopta color amarillo.
Clase C: Valores del índice NDVI comprendidos entre los valores de 0,4 y 0,6 estos valores
se corresponden con vegetación desarrollada o coberturas vegetales altas. Esta clase en la
representación gráfica adopta color verde.
Estas clases permitieron la creación de una secuencia de imágenes (Figura 4.18) en las que
es posible identificar la variación de la cobertura vegetal en el área de la cuenca,
identificando periodos de máxima cobertura en la primavera y principios del otoño.
Mediante la vectorización de estas imágenes de formato rasters se pudo cuantificar el
área correspondiente a cada clase y su porcentaje del total, a lo largo del año en la cuenca
hidrográfica de la laguna Pichi Luan (Tabla 4.3).
41 3 de junio de 2016
Tabla 4.3: Extensión de las diferentes coberturas en el año 2014.
Fecha 06/01/2014 07/02/2014 23/02/2014 11/03/2014 12/04/2014 30/05/2014 01/07/2014
% Clase A
74,49 50,08 27,04 25,46 56,55 88,73 91,95
% Clase B
25,14 49,31 56,58 64,70 41,99 10,48 7,93
% Clase C
0,34 0,61 16,38 9,84 1,46 0,79 0,12
Ha. Clase A
1228,7 826,06 446,02 419,96 932,79 1463,7 1516,8
Ha. Clase B
415,2 813,36 933,29 1067,22 692,62 172,8 130,8
Ha. Clase C
5,6 10,08 270,19 162,3 24,09 13 1,9
Fecha 17/07/2014 02/08/2014 18/08/2014 21/10/2014 06/11/2014 08/12/2014 24/12/2014
% Clase A
92,08 86,28 78,48 22,27 19,68 63,98 65,69
% Clase B
7,83 12,97 20,08 47,73 55,22 34,91 32,68
% Clase C
0,12 0,75 1,45 29,50 25,11 1,11 1,63
Ha. Clase A
1518,9 1423,2 1294,5 375,6 324,6 1055,3 1083,06
Ha. Clase B
129,1 214 331,1 787,3 910,8 575,8 539,1
Ha. Clase C
1,4 12,3 23,9 486,6 414,1 18,4 26,9
42 3 de junio de 2016
Figura 4
.18
: Evolu
ción
de la co
bertu
ra vegetal año
20
14.
43 3 de junio de 2016
La utilización de estos índices permitieron identificar dos periodos en el año en los que la
cobertura seria de regular (Clase B) a buena (Clase C). El primero se ubica en los meses de
febrero y marzo el cual se corresponde con el crecimiento máximo vegetativo de cultivos y
pasturas denominadas como cultivos de verano o plantas C4, lo que vulgarmente es
conocido en el ámbito agrario como la gruesa. El segundo periodo va desde octubre a
noviembre, este se corresponde con la máxima acumulación vegetal de los cultivos de
invierno o plantas C3, conocidas por los productores como la fina. El resto del año se
caracteriza por presentar baja (Clase A) a regular (Clase B) cobertura.
4.6 Información Meteorológica
La precipitación sigue los patrones de distribución normales de la zona de estudio.
Durante el periodo de estudio la misma se concentra en los meses de transición siendo el
mes de octubre y el de abril los de mayor acumulación (Figura 4.19). El total de la
precipitación en el año 2014 fue de 733,4 mm.
Figura 4.19: Precipitación estación meteorológica de El Guanaco, 2014.
La evapotranspiración obtenida mediante la aplicación el método de Penman-Monteith
(FAO-56) permitió determinar una evaporación para la zona de estudio de 1446,3
milímetros. La distribución de la misma fue normal, adquiriendo valores máximos en los
meses de verano y valores mínimos en los meses de invierno (Figura 4.20).
0,0
20,0
40,0
60,0
80,0
100,0
120,0
140,0
160,0
Pre
cip
itac
ión
[mm
]
Año 2014
E.M. El Guanaco
44 3 de junio de 2016
Figura 4.20: Evapotranspiración calculada para el año 2014.
Al contrastar los valores de precipitación y evapotranspiración se pudo identificar
periodos en los que existe un balance positivo, particularmente en los meses de abril y
octubre; contrario a esto del resto de los meses presentaron balances negativos (Figura
4.34). A su vez el balance anual resulto ser negativo, adquiriendo valores que se acercan a
la precipitación total (Tabla 4.4).
Tabla 4.4: Balance hídrico del año 2014.
0,0
50,0
100,0
150,0
200,0
250,0
Evap
otr
ansp
irac
ión
[m
m]
Año 2014
E.M. El Guanaco
Mes Enero Febrero Marzo Abril Mayo Junio Julio Agosto Septiembre Octubre Noviembre Diciembre Total
Precipitacion [mm] 74,6 104,4 109,8 122,6 31,6 2 39 13,8 21 156,6 33 25 733,4
Evaporacion [mm] 250,0 143,7 124,7 69,1 49,7 49,1 46,6 79,1 92,5 135,1 178,3 228,4 1446,3
Balance [mm] -175,4 -39,3 -14,9 53,5 -18,1 -47,1 -7,6 -65,3 -71,5 21,5 -145,3 -203,4 -712,9
Estacion Meteorologica Guanaco
45 3 de junio de 2016
Figura 4.34: Balance hídrico del año 2014.
5. DISCUSIÓN
La caracterización topográfica, hidroquímica y edáfica sumado a la determinación de las
superficie piezométrica y las clases de uso del suelo permitieron comprender los
parámetros físicos y químicos de la cuenca hidrogeográfica de la laguna Pichi Luan. Al
relacionar estos datos con la información meteorológica proporcionada obtenida de la
estación ubicada en la parte alta de la cuenca, se pudo analizar la respuesta de la misma a
diferentes eventos de precipitación. Esto permitió evaluar si la laguna Pichi Luan presenta
comportamiento de recarga puntual del acuífero ubicado en la zona.
Respuesta hidrológica a los eventos de verano
La precipitación acumulada desde el 1 de enero al 21 de marzo del 2014 fue de 236,2
milímetros. Los mismos se distribuyeron en 12 eventos de precipitación (Tabla 5.1), de los
cuales solo dos eventos superaron los 50 mm estos presentaron baja intensidad debido a
que ambos se dieron en un lapso de 48 hs. En esta estación el freatímetro ubicado en la
laguna no presento variaciones, encontrándose el nivel freático en todo momento a
profundidades mayores a los 3 metros bajo el nivel del suelo. Esta respuesta se explica por
el hecho de que durante esta estación se registraron las mayores tasas de evaporación
sumada a la baja intensidad de las precipitaciones; en cuanto a la cobertura vegetal al
inicio de la estación era baja, aumentando hacia fines de la misma (figura 5.1). Bajo estas
condiciones no se pudo comprobar la presencia de procesos de recarga puntual.
-220,0-200,0-180,0-160,0-140,0-120,0-100,0
-80,0-60,0-40,0-20,0
0,020,040,060,0
[mm
]
Año 2014
Balance P vs E
46 3 de junio de 2016
Tabla 5.1: Detalle de eventos de precipitación de Verano.
Fecha Precipitación
17 de Enero 5,2 mm
20 de Enero 10,4 mm
22 y 23 de Enero 59 mm
5 de Febrero 0,6 mm
8 de Febrero 25,2 mm
11 de Febrero 23,8 mm
16 y 17 de Febrero 40 mm
20 y 21 de Febrero 13 mm
24 de Febrero 1,8 mm
8 y 9 de Marzo 4,2 mm
13 y 14 de Marzo 52,8 mm
19 de Marzo 0,2 mm
Figura 5.1: Evolución de la cobertura vegetal la estación de verano año 2014.
0
10
20
30
40
50
60
70
80
90
100
% d
e C
ob
ertu
ra
Verano año 2014
Evolución de la cobertura vegetal
Clase A
Clase B
Clase C
47 3 de junio de 2016
Respuesta hidrológica a los eventos de otoño
La precipitación acumulada entre el 21 de marzo y el 21 de junio del 2014 fue de 208,8
milímetros. Estos se distribuyeron en 12 eventos de precipitación (Tabla 5.2), la
particularidad de esta estación fue que el 78,73 % de las lluvias se concentraron en un
lapso de 15 días acumulando 164,4 mm; puntualmente solo el evento comprendido entre
el 2 y el 8 de abril acumulo 111,8 mm. Posterior a estas precipitaciones el freatímetro de
Pichi Luan mostro un ascenso del nivel freático, pasando de profundidades mayores a 3
metros a 1,2 metros bajo el nivel del terreno. Esta respuesta se explica debido a que en
esta época del año, más precisamente en el mes de abril el balance hídrico fue positivo
superando la precipitación a la evaporación (Figura 4.33). A su vez la cobertura vegetal
comienza a disminuir producto de la maduración de los cultivos de verano disminuyendo
su foliación y así la intercepción de la precipitación (Figura 5.2). Teniendo en cuenta estas
variables la rápida repuesta del nivel freático frente a estos eventos de precipitación
permitió inferir la presencia de procesos de concentración de la escorrentía superficial en
la parte baja de la cuenca y su posterior infiltración.
Tabla 5.2: Detalle de eventos de precipitación en otoño año 2014.
Fecha Precipitación
25 de marzo 23,4 mm
30 y 31 de marzo 29,2 mm
2-8 de Abril 111,8 mm
10 de Abril 1,4 mm
12 de Abril 3 mm
27-29 de Abril 6,4 mm
5 de Mayo 5 mm
10-12 de mayo 16,6 mm
20 y 21 de Mayo 10 mm
7 de junio 1 mm
12 de junio 0,6 mm
15 y 16 de junio 0,4 mm
48 3 de junio de 2016
Figura 5.2: Evolución de la cobertura vegetal la estación de otoño año 2014.
Respuesta hidrológica a los eventos de invierno
La precipitación acumulada desde el 21 de junio al 21 de septiembre del 2014 fue de 73,6
milímetros. Los mismos se distribuyeron en 7 eventos de precipitación (Tabla 5.3). Estas
lluvias de invierno se caracterizaron por su baja intensidad, siendo el mayor evento el del
14 y 15 de julio con 24,4 mm. Durante esta estación el freatímetro ubicado en la laguna no
presento variaciones de importancia iniciando y finalizando el ciclo con 1,5 metros bajo
nivel del suelo. Esta estabilidad del nivel freático existió debido a que si bien el balance
hídrico para esta estación resulto negativo, la mayor parte del agua precipitada fue
incorporada al complejo del suelo. Esto fue posible debido a la baja cobertura vegetal que
caracteriza este ciclo (Figura 5.3), la cual es producto de las siembras de cultivos de
invierno trigo y cebada; como así también materia seca producto de los residuos de
cultivos de verano. Bajo estas condiciones no fue posible identificar procesos de recarga
puntual en la laguna de Pichi Luan.
0
10
20
30
40
50
60
70
80
90
100
11-mar 12-abr 30-may
% d
e C
ob
ertu
ra
Otoño año 2014
Evolución de la cobertura vegetal
Clase A
Clase B
Clase C
49 3 de junio de 2016
Tabla 5.3: Detalle de eventos de precipitación en invierno año 2014.
Fecha Precipitación
4-6 de Julio 6,4 mm
12 de Julio 8,2 mm
14 y 15 de Julio 24,4 mm
22 y 23 de Agosto 11,6 mm
25 de Agosto 2,2 mm
2-5 de Septiembre 20,6 mm
9 de Septiembre 0,2 mm
Figura 5.3: Evolución de la cobertura vegetal la estación de invierno año 2014.
Respuesta hidrológica a los eventos de primavera
La precipitación acumulada desde el 21 de septiembre al 21 de diciembre del 2014 fue de
205,2 milímetros. Los mismos se distribuyeron en 14 eventos de precipitación (Tabla 5.4).
La particularidad de estas lluvias es que solo en el mes de octubre se concentró el 76,31 %
de las mismas, acumulando un total de 156,6 mm. Puntualmente las precipitaciones del 2
al 5, y del 27 al 29 de Octubre generaron un aumento importante del nivel freático,
observado en le freatímetro ubicado en la laguna. El mismo alcanzo niveles que
alcanzaron la superficie aumentando desde 1,5 metros bajo el nivel del suelo a estar a
nivel con el mismo, siendo estos los inicios de la formación de la laguna temporal de Pichi
Luan. En cuanto a la cobertura vegetal al inicio de la estación presentaba valores bajos, la
0
10
20
30
40
50
60
70
80
90
100
01-jul 17-jul 02-ago 18-ago
% d
e C
ob
ertu
ra
Invierno año 2014
Evolución de la cobertura vegetal
Clase A
Clase B
Clase C
50 3 de junio de 2016
misma aumento en los meses de octubre y noviembre llegando a su máxima expresión (en
correspondencia con la época de mayor precipitaciones), para luego disminuir hacia
finales de diciembre (Figura 5.4). Bajo estas condiciones se evaluó la presencia de un
proceso de recarga puntual en las zonas bajas de la cuenca. El mismo fue comprobado
mediante la determinación piezométrica en las transectas ubicadas en la laguna temporal
posterior a eventos de precipitación. Donde se pudo comprobar la formación de un domo
temporal (Figura 5.5, 5.6, 5.7, 5.9), ubicado debajo de la laguna el cual tiene su origen en
la concentración de la precipitación captada por la cuenca y su posterior infiltración.
Tabla 5.4: Detalle de eventos de precipitación en primavera año 2014.
Fecha Precipitación
29 de Septiembre 0,2 mm
2-5 de Octubre 89,4 mm
15 de Octubre 14,4 mm
17 de Octubre 0,6 mm
23-25 de Octubre 13 mm
27 de Octubre 32,2 mm
30 de Octubre 1 mm
2 y 3 de Noviembre 7,8 mm
21 de Noviembre 9 mm
23 de Noviembre 1 mm
27-29 de Noviembre 15,2 mm
12 de Diciembre 0,2 mm
15 de Diciembre 15 mm
19 de Diciembre 0,2 mm
51 3 de junio de 2016
Figura 5.4: Evolución de la cobertura vegetal la estación de primavera año 2014.
Figura 5.5: Comparación topografía (Marrón) y nivel Piezométrico (Azul) de la Transecta A,
posterior a la precipitación del 2-5 de Octubre, iniciando en el freatímetro extendiéndose
al S.S.O. Escala en metros [m], véase figura 4.4.
0
10
20
30
40
50
60
70
80
90
100
18-ago 21-oct 06-nov 08-dic
% d
e C
ob
ertu
ra
Primavera año 2014
Evolución de la cobertura vegetal
Clase A
Clase B
Clase C
52 3 de junio de 2016
Figura 5.6: Comparación topografía (Marrón) y nivel piezométrico (Azul) de la Transecta B,
posterior a la precipitación del 2-5 de Octubre, iniciando en el freatímetro extendiéndose
al S.O. Escala en metros [m], véase figura 4.4.
53 3 de junio de 2016
Figura 5.7: Comparación topografía (Marrón) y nivel piezométrico (Azul) de la Transecta A
en su sección más próxima al freatímetro, posterior a la precipitación del 27 de Octubre,
iniciando en el freatímetro extendiéndose al S.O. Escala en metros [m], Véase figura 4.4.
54 3 de junio de 2016
Figura 5.8: Comparación topografía (Marrón) y nivel piezométrico (Azul) de la Transecta C,
posterior a la precipitación del 27 de Octubre, iniciando en el freatímetro extendiéndose
al Norte. Escala en metros [m], Véase figura 4.4.
55 3 de junio de 2016
Figura 5.9: Comparación topografía (Marrón) y nivel piezométrico (Azul) de la Transecta D,
posterior a la precipitación del 27 de Octubre, iniciando en la intersección con la transecta
C y terminando al Este de la mima. Escala en metros [m], Véase figura 4.4.
56 3 de junio de 2016
6. CONCLUSIONES
Este estudio evidenció la presencia de procesos hidrogeológicos que favorecerían la
recarga puntual del acuífero freático en las lagunas temporales del centro-este de la
provincia de La Pampa. Esto es posible debido a la redistribución de los excedentes
pluviales a zonas deprimidas por efectos del relieve de la cuenca y su posterior infiltración
en el área central de las lagunas temporales.
Si bien el balance hídrico es negativo durante la mayor parte del año, existe un balance
positivo en los meses de transición donde las precipitaciones son intensas. Es en estos
meses donde la precipitación excede la evaporación dando lugar a los procesos de
recargar puntual. Bajo estas condiciones los procesos de redistribución de la precipitación,
se hacen evidentes y generan una mayor dinámica en el nivel freático en inmediaciones
de la laguna temporal de Pichi Luan.
La cuenca hidrográfica de la laguna temporal de Pichi Luan, se caracteriza por el drenaje
centrípeto de carácter endorreico de un área 1649,5 ha. La misma se encuentra aislada
superficialmente de las lagunas circundantes, aunque subterráneamente puede tener
interacciones.
La hidroquímica de la cuenca se caracteriza por la presencia de zonas de baja salinidad
ubicadas en la parte cóncava de la misma. Especialmente las muestras en el área de la
laguna temporal presentan menor concentración de sales, esto es producto de la
lixiviación de las mismas generada por procesos de recarga puntual.
Los suelos de la cuenca se clasifican en francos-arenosos en las partes altas y franco
limosos en las depresiones. Esta clase de suelos se caracteriza por sus buenos parámetros
de infiltración del agua precipitada, es por eso que la misma debe ser intensa para que
existan procesos de escorrentía que generen recarga preferencial.
La piezometría de la zona permite inferir un flujo de agua subterránea hacia la parte
inferior de la cuenca. Este flujo se revierte en épocas donde el balance hídrico es positivo,
formándose en la zona central de la laguna un domo subterráneo, producto de la recarga
puntual.
La cobertura del suelo presentó sus máximos en los meses de transición, conjuntamente
con los meses de mayor precipitación atenuando la escorrentía que pudieran generar las
mismas. A su vez la cobertura en la parte baja de la cuenca fue superior a la de las zonas
periféricas. Esto se debe a que producto de la topografía existe un balance hídrico positivo
que permite inferir la presencia de procesos de redistribución de la precipitación pieza
clave que favorece la recarga puntual.
57 3 de junio de 2016
7. BIBLIOGRAFIA
Auge M., (2008) “Apuntes de Hidrología – Curso 2008”. Departamento de Geología. FCEN-
UBA 1-372. Ined. Bs. As.
Berthold, S., Bentley, L.R., Hayashi, M., (2004). Integrated hydrogeological and geophysical
study of depression-focused groundwater recharge in the Canadian prairies. Water
Resources Research Vol. 40.
Bronstert, A. & Plate, E.J., (1997). Modelling of runoff generation and soil water dynamics
for hillslope and micro-catchments. Journal of Hydrology, 198, 177-195.
Beven, K. (2002). Runoff production in semiarid areas. In: Dryland rivers: hydrology and
geomorphology of semi-arid channels. Bull, L. and Kirkby, M. (Eds.). 57- 116. John Wiley,
Chichester, U.K.
Casagrande, G y H. Conti (1980) “Clima”. En: Inventario Integrado de los Recursos
Naturales de la Provincia de La Pampa. INTA, Gob. DeLa Pampa y UNLPam (Eds.). p. 19-32.
Buenos Aires, Argentina
Castillo, V.; A. Gómez Plaza y M. Martinez Mena. (2003). The role of antecedent soil water
content in the runoff response of semiarid catchaments: a simulation approach. Journal of
Hydrology 284:114-130
Cazenave, H., Dillon, B., Sevilla, J. (1987) El grupo de lagunas del guanaco considerado
como posible sistema hidrográfico. D.A.Y.C.H., A.P.A.
Custodio, E., (1997). Recarga a los acuiferos: aspectos generales sobre el proceso, la
evaluación y la incertidumbre, en La evaluación de la recarga a los acuíferos en la
planificación hidrológica. E. Custodio, M.R. Llamas, J. Samper. Pag.19-39.
Custodio, E., (2010). Aguas Subterráneas en zonas llanas. I Congreso Internacional de
Hidrología de Llanuras. Tomo I: 21-32, Azul, Argentina.
Famiglietti, J.S., Rudinicki, J.W. & Rodell, M., (1998). Variability in surface moisture content
along a hillslope transect: Rattlesnake Hill, Texas. Journal of Hydrology, 210, 259-281.
Peña Zubiate C., Maldonado Piñedo D., Martínez H. y Hevia R. (1980). “Suelos”. En:
Inventario Integrado de los Recursos Naturales de la Provincia de La Pampa. INTA, Gob.
deLa Pampa y UNLPam (Eds.). p. 300-301. Buenos Aires, Argentina.
58 3 de junio de 2016
Salazar Lea Plaza J. (1980). “Geomorfología” En: Inventario Integrado de Recursos
Naturales de la Provincia de La Pampa. INTA, Gob. De La Pampa y UNLPam (Eds.). p. 295.
Buenos Aires, Argentina.
Secco N., Buss G., Dornes P., (2012) Análisis de la respuesta hidrológica en distintos
ambientes pampeanos a la dinámica espacio temporal de la precipitación. Primer
encuentro de investigadores en formación en recurso hídricos I.N.A. Ezeiza, Buenos Aires,
Argentina. 16 pag.
Secco N., Dornes P., Aimar S., (2013) Analisis de la respuesta hidrológica bajo distintos
manejos en una zona agrícola de la provincia de La Pampa. XXIV Congreso Nacional del
Agua. Anales Congreso, 38:12pp. San Juan, Argentina
Taro Uchida ,Ilja Tromp-van Meerveld , Jeffrey J. McDonnell The role of lateral pipe flow
in hillslope runoff response: an intercomparison of non-linear hillslope response Journal of
Hydrology, Volume 311, Issues 1–4, 15 September (2005), Pages 117–133
Van der Kamp G., Hayashi M., (1998), The groundwater recharge function of small
wetlands in the semi-arid northern prairies, Great Plains Res., 8, 39–56.
ZUO, H.J., ZHANG, Q., MA, L.Y., HARTMANN, H., ZHAI, J.Q. & XU, L.G., 2010. Evaluation of
Soil Water Percolation in Response to Different Rainfall Conditions using Hydrus-1D Model
in the Low Hill Red Soil Region of Jiangxi, China. Research Journal of Soil and Water
Management 1(3-4), 76-84.
59 3 de junio de 2016
8. ANEXOS
Anexo I. Tablas de alturas de las transectas topográficas.
Tabla A1. Alturas transecta topográfica A.
Tabla A2. Alturas transecta topográfica B.
Transecta B: Sur Oeste - Freatímetro
Puntos Altura [m] X Y
68 145.82 4384311 5982034
67 145.84 4384380 5982083
66 145.63 4384445 5982119
65 145.31 4384499 5982157
64 145.17 4384543 5982183
63 145.14 4384580 5982205
62 145.09 4384616 5982225
52 145 4384655 5982255
Transecta A: Sur - Freatímetro
Puntos Altura [m] X Y
M-C 148.59 4384436 5981788
97 147.06 4384456 5981829
98 146.21 4384476 5981871
99 145.16 4384494 5981911
59 145.56 4384504 5981926
100 145.65 4384522 5981970
58 145.45 4384549 5982023
57 145.23 4384567 5982063
61 144.99 4384596 5982128
56 144.8 4384613 5982170
55 144.95 4384633 5982210
54 145.02 4384645 5982235
52 145 4384655 5982255
60 3 de junio de 2016
Tabla A3. Alturas transecta topográfica C.
Transecta C: Norte - Freatímetro
Puntos Altura [m] X Y
82 145.17 4384654 5982669
81 145.16 4384654 5982522
80 145.14 4384655 5982449
79 145.15 4384655 5982379
78 145.11 4384656 5982303
52 145 4384655 5982255
Tabla A3. Alturas transecta topográfica D.
Transecta D: Este - Transecta C
Puntos Altura [m] X Y
85 144.91 4384929 5982484
84 145.23 4384867 5982457
83 145.14 4384718 5982454
80 145.14 4384655 5982449
Anexo II. Tabla de conductividad eléctrica.
Tabla A4. Conductividad eléctrica de las perforaciones.
Valores de conductividad eléctrica
Puntos X Y CD [mS/m]
1 4386427 5984801 3.28
2 4381252 5983183 8.42
3 4380581 5982872 3.41
4 4383772 5981645 7.1
5 4384436 5981788 9.8
6 4385047 5981322 10
7 4384561 5982839 4.35
8 4384658 5982259 2.49
9 4386536 5984566 1.92
10 4385562 5983252 5.5
11 4384478 5983696 9.85
12 4383393 5982749 8.91
13 4385906 5981450 9.16
14 4386752 5980952 5.87
15 4381175 5981783 9.34
16 4382502 5981531 1.68
17 4381306 5980794 9.6
61 3 de junio de 2016
Anexo III. Tablas de nivel piezométrico.
Tabla A5. Nivel piezométrico de la laguna temporal 07/10/14.
Nivel Piezométrico: 07/10/2014
Puntos Nivel Estático Ajustado
[m] X Y
52 0 4384655 5982255
54 0.42 4384645 5982235
55 0.58 4384633 5982210
56 0.72 4384613 5982170
57 1.3 4384567 5982063
61 0.8 4384596 5982128
62 0.51 4384616 5982225
63 0.67 4384580 5982205
64 0.91 4384543 5982183
65 1.11 4384499 5982157
66 1.21 4384445 5982119
67 1.6 4384380 5982083
68 1.55 4384311 5982034
M-C 5.84 4384436 5981788
Tabla A6. Nivel piezométrico de la laguna temporal 29/10/14.
Nivel Piezométrico: 29/10/2014
Puntos Nivel Estático Ajustado
[m] X Y
78 0.56 4384656 5982303
79 0.42 4384655 5982379
80 0.26 4384655 5982449
81 0.4 4384654 5982522
82 0.74 4384654 5982669
83 0.18 4384718 5982454
84 0.68 4384867 5982457
85 0 4384929 5982484
52 0.92 4384655 5982255
54 0.92 4384645 5982235
55 0.98 4384633 5982210
56 1.06 4384613 5982170
57 1.59 4384567 5982063
61 1.06 4384596 5982128
62 3 de junio de 2016