Evidencias de interacción entre aguas superficiales y subterráneas a través de las zonas hiporreicas mediante el uso de la hidroquímica y el análisis multivariado en el acuífero de Cuajinicuilapa, Guerrero, México

Edición: Vol.3 Núm.2 mayo-agosto 2012

PDF

 

El acuífero de Cuajinicuilapa en el estado de Guerrero, México, se ubica a la salida de la cuenca hidrográfica del río Quetzala, donde se ha depositado una serie de rellenos aluviales que cubren una fosa tectónica, lo que permite almacenar agua subterránea de donde se extraen importantes gastos para el riego de cultivos; dados los escurrimientos superficiales y el régimen de lluvia, el acuífero no presenta problemas de abatimiento, sin embargo, la calidad fisicoquímica se ha visto ligeramente alterada, en especial por cuerpos de agua superficiales que presentan eutrofización debido a que han quedado aislados de las corrientes superficiales y sólo se renuevan cuando se presentan grandes avenidas, formando en el subsuelo zonas hiporreicas, que influyen en la calidad fisicoquímica del agua subterránea de las regiones aledañas, por lo que fue necesario conocer el modelo geohidrológico del área y la aplicación de herramientas estadísticas para su identificación, apoyados en un sistema de información geográfica (SIG) levantado por el Instituto Nacional de Estadística y Geografía (INEGI)   Palabras clave: zona hiporreica, SIG, fosa tectónica, análisis multivariado		The Cuajinicuilapa aquifer in the state of Guerrero, Mexico, is located off the Quetzala River basin, where they have lodged a series of alluvial fills covering a rift valley, which can store groundwater are extracted substantial volumes for the irrigation of crops, given the runoff and the rainfall, the aquifer is not a problem of depression, but the physicochemical quality has been impacted, particularly for surface water eutrophication presented because have been isolated from the surface currents and are renewed only when there are major avenues, forming hyporheic zones, which influence the chemical quality of groundwater in the surrounding areas, so it was necessary to know the area hydrological modelling and implementation statistical tools for identification, supported by a GIS built by INEGI. Key words: hyporheic zone, GIS, rift valley, multivariate analysis

 

 

Nota: los autores agradecen a la Dirección de Información de Recursos Naturales y Medio Ambiente del INEGI por permitir el uso de la información para la elaboración del presente trabajo, así como a los árbitros anónimos quienes, con sus comentarios, ayudaron a mejorarlo de manera sustancial.

Introducción

Para entender la interacción, de inicio compleja, entre el agua superficial y subterránea, se hace necesario conocer los factores climáticos, morfológicos, geológicos y bióticos a través de un modelo hidrogeoecológico en el que se consideren los mecanismos de recarga del acuífero, sus interacciones, la importancia ecológica y el impacto que pudiera presentarse en la calidad del agua subterránea del área.

El acuífero de Cuajinicuilapa se localiza al sur de la República Mexicana, en la Costa Chica de Guerrero; se encuentra formado por rellenos aluviales que cubren una fosa tectónica; su condición es de tipo libre donde la recarga principal ingresa al sistema de manera directa a través de las corrientes superficiales que discurren por la planicie aluvial por medio de cuerpos de agua (lagos) que se han formado en los meandros abandonados y por flujos horizontales subterráneos que alimentan al acuífero.

Respecto a su calidad fisicoquímica, ésta se ha visto influenciada por la recarga que el acuífero recibe a través de los cuerpos de agua superficiales que han quedado aislados del flujo transitorio (originando la proliferación de lirios acuáticos) y que se manifiesta en las características del agua subterránea de las zonas adyacentes, sin incidir de manera alarmante en el acuífero principal debido a la cantidad de líquido infiltrado a partir de las corrientes superficiales y por flujos subterráneos, que diluyen las concentraciones de aniones y cationes principales, identificado mediante la interpretación de los parámetros fisicoquímicos reportados en un SIG levantado por el INEGI (INEGI, 2010, en proceso de liberación).

Localización y geología

La zona de estudio se localiza entre las coordenadas 16º21’12” y 17º17’33” de latitud Norte y 97º46’14” y 97º52’08” al Oeste del meridiano de Greenwich. Se encuentra inmersa en la provincia fisiográfica de la Sierra Madre del Sur, que se extiende desde el Puerto de San Blas, Nayarit, hasta el Istmo de Tehuantepec en Oaxaca; la parte centro-norte del área o zona serrana está clasificada dentro de la subprovincia Cordillera Costera del Sur y la parte media-baja de la cuenca se ubica dentro de la subprovincia Costas del Sur (INEGI, 2001).

Los sistemas de topoformas presentes en el área son muy variados: destacan las sierras altas complejas, las sierras bajas complejas, el lomerío con llanuras, el lomerío con cañadas, los valles ramificados con lomerío y la llanura costera con lomerío, entre otras (INEGI, 2001).

Las rocas que afloran en la zona de estudio abarcan edades desde el Jurásico hasta el Cuaternario (INEGI, 2010). Dentro del primero, se ha catalogado al basamento de la formación aluvial que aloja al acuífero, representado por un gneis cuarzofeldespático ubicado en la zona centro del área, afectado por intrusiones graníticas- granodioríticas del Terciario; hacia el NO-N y NE, bordeando la cuenca, se localiza un granito- granodiorita del Terciario en contacto hacia el NE con un granito del Jurásico-Cretácico; al SE y centro del área se presentan afloramientos de granitos de edad jurásica-cretácica, infrayaciendo a areniscas del Terciario superior y algunos conglomerados del Cuaternario. En la parte superior de la columna estratigráfica se encuentra el material aluvial clasificado como del Cuaternario, que rellena la parte baja de la cuenca hidrográfica, donde se han depositado los materiales producto de la degradación de las rocas circundantes (limos, arenas y gravas); en menor proporción, hay algunos sedimentos de origen lacustre y litorales, en especial a la salida de la cuenca o en contacto con la zona marina (ver figura 1).

Geología estructural y modelo geohidrológico conceptual

De manera estructural, la zona rocosa muestra un patrón de fracturamiento con una ligera tendencia NO-SE y NE-SO, así como algunos diques que se presentan en estas mismas direcciones; sin embargo, como se sitúan en la parte alta de la sierra, éstas sólo dan origen a pequeños manantiales que alimentan a los ríos Quetzala y Santa Catarina, corrientes principales que se unen al norte de la planicie aluvial para adentrarse en ésta y discurrir por la planicie aluvial (ver figuras 1, 3 y 4), de donde el acuífero recibe el aporte principal de su recarga, sin menospreciar las entradas que se llevan a cabo al SE del área a través de flujos subterráneos y que contribuyen a resarcir la extracción y diluir las concentraciones minerales ocasionadas por la evaporación en los cuerpos de agua superficiales que se infiltran hacia el subsuelo.

En la parte baja de la cuenca, cubierta por el relleno aluvial, no se identifica de forma superficial alguna estructura geológica de trascendencia, pero es probable que la acumulación del relleno aluvial se encuentre asociada a una fosa tectónica, que se deduce a partir de las profundidades y la forma del basamento cristalino identificado en trabajo de la extinta Secretaría de Agricultura y Recursos Hidráulicos (SARH, 1984), actualmente Comisión Nacional del Agua (CONAGUA).

De acuerdo con la metodología para calcular los balances de agua subterránea que lleva a cabo la CONAGUA en los acuíferos nacionales, el de Cuajinicuilapa no se ha visto alterado por las extracciones que hacen los usuarios del agua subterránea debido a que “…existe disponibilidad para nuevas concesiones en la unidad hidrogeológica…”, con poco más de 140 Mm³/año (CONAGUA, 2002); asimismo, los valores de transmisividad calculados para la zona granular oscilan entre 7 y 40 x 10^-3 m²/s, que los clasifica como materiales porosos intergranulares, constituidos por gravas y arenas (Villanueva e Iglesias, 1984) con un coeficiente de almacenamiento estimado en 0.15 (SARH, 1984).

Respecto a la disponibilidad del líquido en la zona de estudio, la CONAGUA (2003) la considera, a nivel nacional, como la tercera con mayor disponibilidad en agua superficial debido a las precipitaciones y volúmenes de escurrimientos, donde sólo se aprovecha 5% del total que circula por la zona, dado que los ríos casi no han sido impactados por desviaciones o extracciones, lo cual origina un aporte constante debido a que sólo se ha construido una presa derivadora sobre el río Cortijo (norte de Cuajinicuilapa) que, por lo suave de la pendiente del terreno, se levantó para desviar una pequeña parte del gasto de escurrimiento que se llena muy rápido y se desborda, por lo que prácticamente es de circulación continua.

Con base en el modelo geohidrológico conceptual del acuífero, por los valores de transmisividad hidráulica reportados y las corrientes superficiales que drenan la zona se deduce que la unidad hidrogeológica es de tipo libre, donde la recarga principal se lleva a cabo a lo largo de los lechos de las corrientes superficiales a través de la superficie de contacto definida entre los materiales granulares permeables y las rocas impermeables, así como por infiltraciones de agua de lluvia precipitada de manera directa en la superficie del valle (ver figura 2); su descarga se produce por evapotranspiración en áreas donde los niveles freáticos se presentan someros y, en menor medida, por bombeo. En reporte de la SARH (1984), las recargas y descargas son numéricamente equivalentes, del orden de 317 Mm³/año.

La creciente demanda social de un medioambiente más limpio ha impuesto en la planificación hidráulica la consideración de que en los cauces circulen, al menos, lo que se ha dado en llamar caudales ecológicos o caudales mínimos medioambientales, en el entendido de que sean capaces de mantener el funcionamiento, composición y estructura del ecosistema fluvial con el fin de sostener la integridad ecológica del sistema, sin disminuir su complejidad estructural y biodiversidad.

En la figura 3 se presentan de manera esquematizada las corrientes principales, la ubicación de estaciones hidrométricas en el área, así como un resumen de los aforos realizados durante 10 años por la SARH (1984) y que confirma la anterior aseveración.

Método de trabajo

Para conocer el comportamiento de los flujos subterráneos en la zona de estudio, se identificaron, geoposicionaron y se les midió el nivel estático a 39 obras de extracción (pozos y norias), así como el registro de la altura de los brocales donde los había para, después — apoyados en los modelos digitales de elevación— identificar la altitud de los sitios en los que se ubicaban éstos, con lo cual se dedujeron las direcciones de los flujos subterráneos (ver figura 4).

Durante la etapa de campo, en las obras de extracción que se encontraban trabajando (previo análisis en gabinete de la ubicación de los sitios a muestrear para que éste resultara lo más sistematizado posible), se colectó una muestra de agua de un litro, con base en el protocolo de muestreo desarrollado en el Laboratorio de Análisis de Materiales del INEGI, que fueron analizadas para determinar el contenido de aniones y cationes principales (ver tabla 1), además de algunos calificadores, como: dureza, relación de adsorción de sodio (RAS), calidad del agua para riego, familia y agresividad, entre otros; asimismo, en las descargas directas se le determinó la temperatura, conductividad eléctrica y pH mediante un conductivímetro portátil, previa calibración del equipo mediante soluciones buffer.

Discusión y resultados

Con base en la información de los niveles estáticos obtenida en campo para la fecha del recorrido, se identificó que la dirección del flujo subterráneo tiende a seguir la pendiente natural del terreno (por lo general NE-SO), como corre el río Santa Catarina, principal afluente en la zona por su desarrollo y el aporte que hace al acuífero y que pudiera ser indicativo de la influencia que tiene sobre éste; aun cuando el nivel estático se ha profundizado de forma ligera, en especial en la zona centro del acuífero aluvial, no se identifican conos de abatimiento que pudieran ocasionar algún problema de intrusión marina debido a la diferencia de cargas hidráulicas (ver figura 4) y que se confirmó mediante el análisis de la posible presencia de algunos elementos, como: el ion Na⁺, Cl^– y la conductividad eléctrica, los cuales pudieran auxiliar en la identificación de alguna anomalía hidroquímica asociada a esta naturaleza (Neri, 2007).

Analizando los resultados reportados por el laboratorio de las 39 muestras ingresadas, se comprobó el balance de electro-neutralidad entre aniones y cationes, cuyos errores resultaron inferiores a 6%, por lo que fueron considerados como aceptables (Hem, 1985), procediendo a la identificación de los sitios donde se presentaban las concentraciones más significativas mediante el uso de la Geoestadística¹ (ver figura 5). El método utilizado para la elaboración de las configuraciones de los parámetros fisicoquímicos fue el Kriggin, que es el más adecuado para la información que deseamos representar.

De las anomalías identificadas en el agua subterránea mediante esta metodología, destacan las asociadas a los cuerpos de agua superficiales formados en los meandros de los ríos, donde se almacena el líquido de los periodos de grandes avenidas que, al quedar aislado de las corrientes superficiales de régimen transitorio, la temperatura y la carga de nutrientes que éste trae origina un proceso de eutrofización² ocasionado por las descargas de aguas residuales domiciliarias sin tratamiento de las localidades ubicadas aguas arriba, el uso de fertilizantes y el lavado de suelos por la deforestación.

Asociado a este fenómeno se encuentra la formación de zonas hiporreicas en el subsuelo, que es la interface que vincula el cuerpo de agua superficial y la zona saturada del acuífero, donde se lleva a cabo la interacción o mezcla entre el agua superficial y subterránea (ver figura 6), cediendo algunas de las características fisicoquímicas de la superficie hacia el acuífero, como se identifica en las configuraciones de algunos de los elementos químicos de la figura 5. Desde el punto de vista biológico, estas zonas generan un gran interés debido a los diferentes ecosistemas que se desarrollan en esta capa que, a su vez, sirven de protección para los organismos que las habitan durante las etapas de grandes avenidas.

Para conocer la tendencia de posibles procesos de mezclas y consecuente evolución química (registrada en las muestras de agua subterránea analizadas por el laboratorio del INEGI), se utilizó el diagrama de Piper, en el cual se puede apreciar que, debido a la posición que guardan en el gráfico, éstas tienden a ubicarse en la zona de intercambio catiónico, así como de procesos de mezclas principalmente, con un ligera tendencia hacia la disolución de yesos-sulfuros presentes en las rocas con las que entra en contacto durante su trayecto (ver figura 7).

Asimismo, mediante estos diagramas triangulares se identificaron las familias de agua (ver tabla 1), donde su agrupamiento está en función del catión y anión predominantes y que son los que le dan nombre y apellido a la muestra de agua. La familia dominante fue la cálcica sódica-bicarbonatada (16 muestras), en segundo lugar se ubicó la mixta-bicarbonatada (15 muestras), sobre todo.

Con el propósito de diferenciar mediante la química del agua la recarga con la que se sustenta el acuífero de Cuajinicuilapa, se utilizó la metodología desarrollada por Mifflin (1968), en la que considera las órdenes de magnitud de los indicadores hidrogeoquímicos Na⁺+K⁺ vs. Cl^–+SO₄^– para clasificar los flujos en: locales, intermedios y regionales.

En el primero se ubicaron 30 muestras, siete de ellas se catalogaron como de flujos intermedios y dos más se posicionaron en el extremo de la gráfica, sea por alimentación de flujos regionales o por fuentes externas que están modificando la calidad del agua subterránea en esos sitios y que no fueron identificadas durante el recorrido de campo (ver figura 8).

Como complemento de la información, se llevó a cabo un análisis multivariado con los resultados de las muestras, que es una herramienta estadística que proporciona resultados importantes en el estudio de variaciones geoquímicas muy sutiles en el agua subterránea, las cuales no pueden ser obtenidas por métodos convencionales y que permite diferenciar acuíferos de origen similar, debido a que se utiliza casi toda la información obtenida en campo y la reportada por el laboratorio.

Para su elaboración, se diseñó una matriz A x B, donde A simboliza los aprovechamientos muestreados en el acuífero de Cuajinicuilapa, representados por 15 norias y 24 pozos y B, las variables fisicoquímicas (aniones y cationes principales: SO₄ ^2-, HCO₃^–, CO₃ ^2-, Cl^–, NO₃ ^–, Ca²⁺, Mg²⁺, Na⁺ y K⁺), temperatura, pH y conductividad eléctrica (ver tabla 1); el software que se utilizó para el análisis de la información fue el PASW Statistics 18 (SPSS).

Mediante la integración de la información, se identificaron dos grupos o facies principales (ver figura 9). En la facie 1 se agruparon 27 muestras y en la 2, sólo 12, con sus respectivos subgrupos (4), que fue hasta donde se llevó la descripción que hacemos y que se indica su ubicación física en la parte superior del dendrograma (ver figuras 9a, 9b, 9c y 9d).

En la figura 9a se aprecia que las muestras se ubican en toda la zona del acuífero aluvial, lo cual pudiera indicar que se encuentran asociadas al desarrollo de los escurrimientos en la planicie aluvial; las agrupadas en la figura 9b se asocian con la llegada del río Quetzala y únicamente una de ellas se ubica en el extremo opuesto, por donde llega el río Santa Catarina, así que pudiera ser una característica del agua que discurre de forma superficial y se infiltra de inmediato a la planicie aluvial; en la 9c, la distribución de las muestras se presenta con la llegada de los ríos Santa Catarina y Cortijo (oriente del área); la 9d representa los pozos que se distribuyen en la zona centro-norte, prácticamente en la cabecera del acuífero aluvial y que pudiera indicar zonas de recarga con una mayor concentración de algunos componentes químicos.

Conclusiones

El acuífero de Cuajinicuilapa ubicado en la Costa Chica del estado de Guerrero recibe su recarga sobre todo a través de las corrientes superficiales que drenan los valles debido a su condición de libre, por lo que el mantenimiento de su buena calidad fisicoquímica redundará en la calidad del agua subterránea.

El modelo geohidrológico conceptual y el graficado de índices hidrogeoquímicos indican la presencia de flujos locales, intermedios y/o regionales que, aunque no son determinantes en la calidad fisicoquímica del agua subterránea debido a procesos de mezclas, se hace necesaria su identificación a mayor detalle para prever su comportamiento dentro del acuífero.

Es de suma importancia identificar en los estudios hidrológicos de aguas subterráneas el papel y extensión de las zonas hiporreicas, ya que éstas pudieran influir en la calidad fisicoquímica del agua subterránea y, de ser posible, que ambas sean manejadas como un mismo sistema hidrogeológico con el fin de proteger los acuíferos y preservar estas zonas debido a su importancia biológica.

La implementación de técnicas estadísticas en la identificación de la evolución fisicoquímica del agua subterránea cumplió con las expectativas de apoyo para la definición de sus orígenes o etapas de circulación, por lo que ésta se pudiera implementar para dar mayor uso a los datos contenidos en el SIG levantado por el INEGI.

Por último, se hace necesario concientizar a las autoridades de las localidades que se ubican dentro de la cuenca hidrográfica para que se comprometan al confinamiento o tratamiento de las aguas residuales domiciliarias con el propósito de que no sean descargadas de manera directa a los escurrimientos superficiales en detrimento de las comunidades que se encuentran en la parte baja o a la salida de la cuenca.

Referencias

Brito-Castillo, L. et al. “Groundwater differentiation of the aquifer in the Vizcaino Biosphere Reserve, Baja California Peninsula, Mexico”, en: Geofísica Internacional. 49 (4), 167-179, 2010.

CONAGUA. Determinación de la disponibilidad de agua en el acuífero Cuajinicuilapa, estado de Guerrero. 2002.

____Estudio de actualización de mediciones piezométricas en los acuíferos de Tehuantepec, Nochixtlán y Cuajinicuilapa del estado de Guerrero. COPEI Ingeniería, 2003, 128 p.

Dahm, C. N. et al. “Nutrient dynamics at the interface between surface waters and groundwaters”, en: Freshwater Biology. 40, 427-451, 1998.

Fox, G. A. & D. S. Durnford. “Unsaturated hyporheic zone flow in stream/aquifer conjunctive systems”, en: Advances in Water Resources. 26: 989-1000,2003.

González, A. Fernando. Guía rápida del programa SPSS ver. 13.0 a 15.0. Curso Intensivo de Posgrado. UACH, México, 2010.

Hem, J. D. “Study and interpretation of the chemical characteristics of natural water”, en: US. Geological Survey Water-Supply. Paper 2254, 1985.

Hayashi, M. and D. O. Rosenberry. “Effects of groundwater exchange on the hydrology and ecology of surface water”, en:Japanese Journal of Groundwater Hydrology. 43, 327-341. Japanese Association of Groundwater Hydrology, 2001.

Hergt, T. et al. Análisis multivariado en la definición de sistemas de flujo de agua subterránea en San Luis Potosí, México. 2007.

Hernández, H. et al. “Determinación de las características del agua subterránea del graben de Aguascalientes, usando análisis multivariado”, ponencia presentada en: Groundwater flow understanding from local to regional scales. XXXIII Congress IAH & 7º Congress ALHSUD, 2004.

INEGI. Sistema de clasificación y levantamiento fisiográfico, escala 1:1 000 000. Aguascalientes, México, INEGI, 2001, 21 pp.

____Conjunto de datos geológicos Chilpancingo, E14-08, escala 1:250 000. En formatos impreso y digital. México, INEGI, 2010.

Mifflin, M. D. “Delineation of groundwater flow systems in Nevada”, en: Report Series H-W. No. 4, p. 111, Desert Research Institute Technical, 1968.

______ “Region 5, Great Basin”, en: Back, W. et al. (eds.). Hydrogeology. Geological Society of America. 69-78, 1988.

Neri Flores, Iris. Identificación de zonas vulnerables a intrusión salina en una zona costera de Veracruz. Universidad Nacional Autónoma de México. Programa de Posgrado en Ciencias de la Tierra. Instituto de Geofísica, 2007, 148 p.

Rajesh Reghunath, T. R. Sreedhara Murthy y B. R. Raghavan. “The utility of multivariate statistical techniques in hydrogeochemical studies: an example from Karnataka, India”, en: Water Research. 36, 2437-2442, 2002.

Salman, S. R., Abu Ruka’h. “Multivariate and principal component statistical analysis of contamination in urban and agricultural soils from north Jordan”, en: Environmental Geology. 38 (3), p. 265-270, 1999.

SARH. Estudio geohidrológico preliminar de la zona de Cuajinicuilapa, estado de Guerrero. Elaborado bajo contrato por Lesser y Asociados, SA de CV, 1984.

Simeonov, V. et al. “Assesment of the surface water quality in Northern Greece”, en: Water Research. 37, 4119-4124, 2003.

Sophocleous, M. “Interactions between groundwater and surface water: the state of the science”, en: Hydrogeology Journal. 10:52-67, 2002.

Tóth, J. “Las aguas subterráneas como agente geológico: causas, procesos y manifestaciones”, en: Boletín Geológico y Minero. Vol. 111-4, 9-26, 2000.

Villanueva Martínez, Manuel y Alfredo Iglesias López. Pozos y acuíferos. Técnicas de evaluación mediante ensayos de bombeo. Instituto Geológico y Minero de España, 1984, 313 p.

Winter, T. C. et al. “Ground water and surface water. A single resource”, en: U.S. Geological Survey. Circular 1139, 1998. Williams, D. D. “Nutrient and flow vector dynamics at the hyporheic/groundwater interface and their effects on the interstitial fauna”, en: Hydrobiologia. 251:185-198, 1993.

Wondzell, S. M. and F. J. Swanson. “Floods, channel change and the hyporheic zone”, en: Water Resources Research. Vol. 35, No. 2, p. 555-567, 1999.

¹ Conjunto de técnicas para resolver problemas de estimación que involucran variables espaciales georrefenciadas, representados por algunos iones, como: los NO₃^–, Cl^–, SO₄=, Na⁺ y los STD, entre otros.

² Proliferación excesiva de plantas acuáticas debido al enriquecimiento del agua superficial con nutrientes primarios, en especial nitrógeno y fósforo presentes en el agua superficial (Winter et al., 1998).