1. Antecedentes

Los asentamientos humanos se han establecido históricamente en torno a los recursos hídricos superficiales: costas, ríos, lagos y manantiales (^{Chapelle, 1997}; ^{Pastoreet al., 2010}). A su vez, los asentamientos humanos han afectado el ciclo hidrológico a través del cambio de uso de la tierra, el cambio climático, el uso excesivo del agua, la redirección del agua y los cambios en la calidad del agua (^{Pastoreet al., 2010}; Vörösmartyet al., 2013). Por lo tanto, una buena comprensión de la interacción agua subterránea - agua superficial ayuda a comprender, analizar y trazar no solo la historia pasada, sino también la futura.

Los usuarios de los recursos hídricos aún visualizan los sistemas de aguas superficiales y subterráneas como elementos independientes. Para efectos legales y de política pública, el agua superficial está desvinculada del agua subterránea (^{Carmona-Laraet al., 2017}).

Es común que las autoridades que administran los recursos hídricos se refiere a la transferencia de agua de un río a un acuífero como “pérdidas” y, en el caso contrario como “recuperaciones”. Esto pone de manifiesto una evidente desvinculación entre ambos sistemas que causa problemas y conflictos en la administración y manejo de los recursos hídricos. Si bien toda Ley Nacional de Aguas debe perseguir la gestión integral de las aguas superficiales y subterráneas, en la práctica no siempre es posible el utilizar de manera coordinada las distintas fuentes de agua para satisfacer las demandas hídricas en algunas regiones del mundo (^{Levy y Xu, 2012}).

El tema de la interacción del agua superficial con el agua subterránea ha recibido mucha atención como pieza clave de investigación ya que incrementa el entendimiento del ciclo hidrológico en una región (^{Linet al., 2018}; ^{Barthel y Banzhaf, 2016}; ^{Levy y Xu, 2012}). En lugar de estudiar cada una de estas fuentes de agua de manera separada, entender su interacción da una mejor idea de las respuestas tanto a la dinámica natural como al cambio global (^{Sahaet al., 2017}; ^{Fleckensteinet al., 2010}).

Asimismo, entender el ciclo hidrológico implica estudiar tanto los receptores del recurso como el agua superficial, el agua subterránea y la atmósfera, y los intercambios de un receptor a otro. Se ha definido al agua subterránea como aquella que satura los poros en el subsuelo. El flujo del agua entre el subsuelo y la superficie está gobernado por la ecuación de Darcy para flujo en medios porosos saturados (ecuación 1).

Donde:qes la descarga específica [L/T],Kes la conductividad hidráulica en el lecho del cuerpo de agua superficial [L/T], eies el gradiente hidráulico en la interfase [-].

Aunque esta ecuación tiene potencial para cuantificar las interacciones, los flujos entre el agua subterránea y el agua superficial no están bien caracterizados (^{Winteret al., 1998}). Por definición, la frontera se encuentra bajo el agua, dificultando las mediciones que al ser puntuales, son insuficientes para caracterizar los flujos de agua al recargar al acuífero o descargar el agua subterránea.

La interfase entre el agua subterránea y el agua superficial ocurre cuando el nivel freático (el límite superior de la zona saturada de un acuífero libre) intersecta la superficie del terreno. En otras palabras, la interfase agua subterránea - agua superficial ocurre en las corrientes de los ríos, lagos, humedales y océanos. Las descargas ocurren debajo del nivel superficial del agua y prácticamente sólo se ven en los manantiales que son descargas puntuales. Independientemente del caso, se sabe que la descarga ocurre porque el agua superficial es una manifestación de la descarga del agua subterránea. Muchas corrientes se secarían entre períodos de lluvia si no recibieran descargas de aguas subterráneas, sin considerar aquellas que reciben descargas inducidas por el hombre o deshielo de glaciares. Al descender el nivel freático, los lagos y humedales pueden secarse. Aún en los grandes cuerpos de agua superficiales, la descarga de agua subterránea se ha reconocido como un componente importante del ciclo hidrológico (^{Moore, 2010}). Así, el agua subterránea es clave para la sustentabilidad de todos los recursos de agua superficiales (^{Alleyet al., 1999}).

La interfase puede involucrar flujos de agua en ambas direcciones, es decir que los cuerpos de agua superficiales pueden ganar así como perder agua proveniente del acuífero. Por lo general, las descargas de agua subterránea ocurren en la orilla de estos receptores, según el concepto de sistema de flujo (^{Winteret al., 1998}; ^{Moore, 2010}; ^{Alleyet al., 1999}; ^{McBride y Pfannkuch, 1975}), que pueden ganar agua (recibiendo flujo del agua subterránea) en una orilla y perder agua (recargar flujo al acuífero) en la otra. Las corrientes también pueden ganar o perder agua y ambos procesos pueden ocurrir a lo largo de diferentes tramos de un río.

Las estimaciones de la descarga de agua subterránea son muy variables para los ríos debido a la amplia variedad de escenarios naturales en los que se encuentran. ^{Zektser y Loaiciga (1993)} proporcionaron estimaciones del porcentaje de la escorrentía de los ríos como flujo base de agua subterránea en todo el mundo que oscilan entre el 35 y el 55 %.

En 24 regiones de los Estados Unidos, delineadas por ^{Winteret al. (1998)}, la contribución del agua subterránea al caudal fluvial osciló entre el 14 y el 90 % con una media del 55 %. Por otro lado, ^{Rudolphet al. (2007)} determinaron que la disminución progresiva de la extensión superficial de la zona lacustre en el sur de la Cuenca del Valle de Toluca, México, podría estar asociada a la explotación intensiva del agua subterránea en las proximidades. Así, la disminución de la recarga de agua superficial y el flujo lateral de las montañas circundantes serían los factores principales que incrementan el drenaje superficial y/o aumentan las pérdidas por infiltración, como resultado de la despresurización de los sedimentos subyacentes a la zona lacustre.

Otro caso de explotación intensiva del acuífero y sus consecuencias ocurre en el área de Los Sotos, del Parque Nacional Doñana (Huelva), España, que ha sufrido cambios notables. y que fueron simulados en función de los efectos a mediano y largo plazo (^{Juárezet al., 2012}). ^{Arumiet al., (2012)} presentan casos de estudio en Chile en los que se enfatiza que el uso de modelos simples como la transposición de caudales, crecidas o supuestos de igualdad de rendimiento específico pueden generar resultados erróneos cuando existen cuencas donde las condiciones geológicas producen importantes trasvases de agua subterránea. Durante el período de flujo base el agua del río generalmente se deriva del agua subterránea, a menos que existan descargas de aguas por actividades humanas.

Debido a la evidente predominancia de los flujos superficiales en el Estado de Tabasco, es común encontrar una amplia variedad de investigaciones enfocadas al estudio de ríos y lagos. En contraparte, la hidrología subterránea ha sido relegada a un segundo plano, incluso para las instituciones públicas relacionadas con la administración del agua. Se conoce poco acerca de la dinámica que presenta el agua subterránea de la región.

Gracias a la abundancia y proximidad de los ríos y lagos, la población tabasqueña no necesitó explotar el recurso del subsuelo sino hasta las décadas de 1970 y 1980 debido al incremento de la población. Este fenómeno obligó a los asentamientos a aumentar el número y el volumen de los pozos de extracción para abastecerse de agua (^{INEGI, 2001}).

La fama de Tabasco como un territorio hídricamente abundante se debe a sus condiciones geológicas - climáticas que resultan favorables para la renovación del agua (^{CONABIO, 2019}; ^{Ramos-Herreraet al., 2012}). En las porciones sur - sureste (en el municipio de Tenosique), y suroeste (en el municipio de Huimanguillo), es donde se presentan las principales elevaciones de terreno que conforman una zona de recarga importante debido a las altas precipitaciones que registran (^{INEGI, 2001}).

Esta investigación se llevó a cabo para comprender mejor el movimiento del agua subterránea y su interacción con el agua superficial. La información y conocimiento generados apoyarán a las autoridades a nivel municipal, estatal y federal en la gestión del agua de la región “Los Ríos”, al poniente del municipio Balancán en Tabasco. Se pretende que este estudio sea útil en la toma de decisiones durante contingencias hídricas y ambientales como son inundaciones por desborde de ríos y el transporte de contaminantes a los cuerpos de agua superficiales y al acuífero.

1.1. ZONA DE ESTUDIO

Tabasco es uno de los territorios más importantes de la República Mexicana en materia hídrica (^{CONABIO, 2019}). Sus numerosos cuerpos de agua superficial, vegetación y consecuentes paisajes naturales le han valido el sobrenombre de “El Edén”. La región de Los Ríos, se ubica en la porción oriental del estado de Tabasco, hacia el poniente del municipio de Balancán (Figura 1). La zona de estudio se encuentra en el área que corresponde en su totalidad a la provincia fisiográfica Llanura Costera del Golfo Sur; estructura llana formada por grandes cantidades de aluvión acarreado desde la Sierra de Chiapas por el Río Usumacinta y sus afluentes (^{CONAGUA, 2015}). Geomorfológicamente la zona está constituida por una llanura costera de relieve escaso con altitudes menores a los 100 m.s.n.m. (metros sobre el nivel del mar) como resultado de la acumulación de grandes depósitos fluviales que la han moldeado. Geológicamente está representada por formaciones sedimentarias del Terciario, como calizas y areniscas, y del Cuaternario, como conglomerados y depósitos palustres, lacustres y aluviales (^{CONAGUA, 2015}).

Figura 1 Localización de la zona de estudio en el municipio de Balancán, Tabasco, México. 

El clima de la región de Los Ríos se caracteriza por ser cálido con régimen normal, debido a sus parámetros de temperatura y humedad (^{CONAGUA, 2015}). La temperatura en Balancán oscila anualmente entre los 21.1 y los 32.6 °C, registrándose mínimas promedio de hasta 18.2 °C y máximas promedio de 36.3 °C. Además, se identifican regímenes anuales promedio para la precipitación (1,670.8 mm) y evaporación (1,158 mm). El proceso de evapotranspiración está principalmente regulado por la temperatura, su valor promedio anual ronda los 1,628 mm; se sugiere que la evapotranspiración real sobrepasa este valor debido a la deficiencia de humedad en los meses de abril y mayo (^{CONAGUA, 2015}).

La Figura 2 muestra un climograma representativo de Balancán, obtenido de un análisis estadístico realizado a los parámetros registrados en estaciones climáticas de la zona de estudio, el cual corrobora que este municipio cumple con el patrón estatal de lluvias y temperatura.

Figura 2 Climograma característico Municipal de Balancán, Tabasco. (Elaboración propia con datos de ^{Díaz-Padillaet al., 2006}). 

Se observa la presencia de precipitación todo el año, siendo más abundante en el período junio - octubre, mismo que el ^{INEGI (2019)} indica como temporada de lluvias; además, se ratifica que la evaporación más pronunciada se registra en los meses con mayores temperaturas (abril y mayo), los cuales coinciden con la temporada de estiaje (noviembre - mayo) en la zona de interés (^{CONAGUA, 2015}).

Según el Instituto Nacional de Estadística y Geografía (INEGI), el territorio de la entidad representa apenas el 1.3% de la superficie nacional y, sin embargo, cuenta con 129,331.2 hectáreas de cuerpos de agua divididas entre 33 ríos y 17 lagunas (^{INEGI, 2014}), que ocupan aproximadamente 5.25 % de la superficie de Tabasco. Entre los ríos principales se encuentran el Grijalva, el Usumacinta y el Mezcalapa (Figura 3).

Figura 3 Distribución espacial de los cuerpos y corrientes de agua en la zona de estudio. 

El sistema geohidrológico subyacente abarca una superficie de 1,329 km²en su totalidad, y cubre parcialmente los municipios de Balancán, Centla, Jonuta y Emiliano Zapata (^{CONAGUA, 2015}).

En su mayoría, los aprovechamientos que han sido monitoreados por la Comisión Nacional del Agua (CONAGUA) desde la década de 1990 se encuentran en el municipio de Balancán (Figura 1). En términos generales, Tabasco puede considerarse un solo acuífero debido a la uniformidad de sus condiciones, con excepción de las zonas con altas elevaciones al sur de la entidad. Debido a la gran capacidad de infiltración del subsuelo, las aguas precipitadas son transportadas a través del acuífero después de la recarga hasta desembocar en el Golfo de México (^{INEGI, 2001}).

2. Métodos

Para el análisis de las variables involucradas en esta investigación se requirieron herramientas computacionales entre las cuales se distinguen un software de interpolación y representación gráfica de los parámetros piezométricos y un sistema de información geográfica (SIG) para la ubicación espacial de los sitios de interés.

La Subgerencia Técnica de la Dirección Local de CONAGUA, en Villahermosa, Tabasco, proporcionó datos del nivel freático históricos del acuífero que representan escenarios temporales en intervalos de tres años para el período que va de 2002 - 2017. La dependencia también suministró la ubicación georreferenciada de cada uno de los pozos censados en la zona.

En cuanto al marco cartográfico, se recuperó información ráster y vectorial del portal electrónico del INEGI. Entre las principales contribuciones figuran: el modelo de elevación digital (DEM, por sus siglas en inglés) de tipo LiDAR (Light Detection and Ranging, por sus siglas en inglés) en alta resolución (cinco metros) a escala 1: 10,000; límites geopolíticos del municipio de Balancán y sus similares circundantes; y mapas de características físicas y demográficas.

Los registros históricos de 49 puntos de monitoreo del agua subterránea fueron complementados con otros 284 puntos propuestos que se localizaron en las márgenes de cuerpos de agua superficial. Al tratarse de un acuífero libre, los puntos complementarios corresponden a puntos de nivel freático cuyos valores fueron estimados mediante el DEM.

Con los datos de profundidad de agua en pozos de la red piezométrica existente, los potenciales hidráulicos del acuífero para cada uno de los escenarios temporales fueron calculados. Los contornos equipotenciales para cada año fueron definidos mediante el análisis de los potenciales hidráulicos. Los datos fueron dispuestos en el área de estudio mediante un mallado que, al ser interpolado, genera las curvas de nivel freático. Adicionalmente, las cuadrículas se promediaron para obtener las curvas características de la zona para el período 2002 - 2017. Los vectores de dirección de flujo fueron determinados mediante la interpretación de cada uno de los escenarios modelados. Los gradientes hidráulicos resultan del cálculo de las pendientes de cada uno de los modelos piezométricos generados. Se han representado gráficamente en mapas para los escenarios correspondientes y para los valores característicos del lapso analizado.

La relación del acuífero con el sistema hidrográfico suprayacente se determinó analizando las diferencias de potenciales hidráulicos entre escenarios anuales consecutivos y comparando el resultado con la localización de los ríos y lagos circundantes. La estimación de los caudales de ingreso y egreso se realizó empleando una simplificación de la Ley de Darcy. Esta primera aproximación de los gastos se calculó mediante datos conocidos de los pozos “Provincia” (caudal de entrada) y “Chablé” (caudal de salida). Éstos corresponden a los puntos extremos de la red de flujo de agua subterránea, tanto aguas arriba como aguas abajo, respectivamente. La transmisividad fue tomada de las pruebas de bombeo ejecutadas por ^{CCT (1991)}, el gradiente hidráulico fue determinado al aproximar sus valores con los mapas correspondientes.

3. Resultados y discusión

Los datos obtenidos de la Subgerencia Técnica de la Dirección Local de CONAGUA, en Villahermosa, Tabasco, fueron depurados para seleccionar una red conformada por pozos de monitoreo ubicados en la zona de estudio. Los registros de la red piezométrica se realizaron de manera intermitente y discontinua, sin un patrón establecido en tiempo, pero contenían suficiente información para el análisis realizado. Es así que se establecieron escenarios temporales clave, según la frecuencia de lecturas, considerando los meses correspondientes a la época de lluvias. Se trató de tener una buena cobertura de la zona de estudio. Los pozos seleccionados fueron ubicados mediante elSIG Global Mapper(^{Blue Marble Geographics, 2016}). Adicionalmente se ubicaron puntos auxilia- res en los márgenes de los ríos y cuerpos de agua perennes, según se muestra en la Figura 4. Estos puntos auxiliares complementaron la información para mejorar la precisión de la modelación ya que se ubicaron considerando una interacción directa entre los sistemas hidrológicos superficial y subterráneo. Los potenciales hidráulicos del acuífero se calcularon restando a la altitud del brocal del pozo de monitoreo la profundidad a nivel estático registrada. Con los valores resultantes, se generaron contornos equipotenciales anuales con el apoyo del programaSurfer, utilizando el interpoladorKriging.

Figura 4 Ubicación de pozos de monitoreo en la zona de estudio. 

Los niveles piezométricos máximos de cada escenario se presentaron en el extremo noreste de la zona de estudio; éstos tienen valores de hasta poco más de 45 m.s.n.m. Los valores mínimos se presentaron en la porción más occidental del municipio de Balancán coincidiendo con las llanuras de menor altitud y con las márgenes del Río Usumacinta, hacia el cual migra el agua subterránea. La Figura 5 muestra el mapa promedio para el período 2008 - 2017.

Figura 5 Potencial hidráulico promedio del Acuífero “Los Ríos”, Balancán, Tabasco (Período 2008 - 2017). 

Se contrastaron las diferencias entre los contornos de curvas equipotenciales para diferentes periodos de tiempo que, al analizar los mapas, ponen de manifiesto una reducción paulatina del nivel freático. En particular, esto puede observarse en el extremo noreste de la zona de estudio hasta aproximarse a las márgenes del Río Usumacinta, principal cauce de la zona, lo cual sugiere que el acuífero aporta agua al sistema superficial suprayacente (Figura 6).

Figura 6 Diferencia del nivel freático entre: (A) 2008 y 2011 y (B) 2014 y 2017. 

De manera similar a los contornos equipotenciales, se generaron mapas con los valores de gradientes hidráulicos calculados para los registros en los años considerados para el estudio. Una síntesis de lo anterior se presenta en la Figura 7, como un mapa de contornos de gradientes hidráulicos promedio del acuífero para el período 2008 - 2017. Se resalta que los gradientes hidráulicos más pronunciados se presentan en la zona noreste donde se encuentran los relieves topográficos con mayor altitud.

Figura 7 Gradientes hidráulicos promedio del Acuífero “Los Ríos”, Balancán, Tabasco (Período 2008 - 2017). 

Se distinguieron dos zonas principales, una al Este, con valores mayores a 50 m/km que presenta variaciones con tendencia a incrementar sin exceder 250 m/km; y otra al oeste, donde se observan las pendientes menos pronunciadas y que también permanecen presentes en algunas zonas que corresponden al Río Usumacinta y demás cuerpos de agua perennes.

Los mapas de curvas equipotenciales permiten identificar las zonas de entrada y salida natural de agua subterránea en la zona de estudio. ^{Peñuela Arévalo y Carrillo Rivera (2013)} aplicaron la teoría de sistemas de flujo como una herramienta para delimitar las zonas de recarga y descarga en la porción centro-sur de la Mesa-Central, México. Mostraron la utilidad de la herramienta analizando indicadores para obtener la localización de áreas prioritarias, aportando una aproximación del funcionamiento del acuífero.

El ingreso de volúmenes importantes de agua subterránea en la zona se da principalmente de tres maneras: (1) El agua subterránea que proviene de porciones circundantes al acuífero que aporta por flujo regional; (2) El agua superficial que escurre en la cuenca hidrológica suprayacente que aporta mediante el cauce al acuífero; y (3) La precipitación pluvial de la zona que se infiltra a través de las unidades hidrogeológicas más permeables. Los egresos del agua subterránea se dan principalmente mediante el flujo natural del acuífero, tal como se ha caracterizado anteriormente (Figuras 5 y 6).

El análisis histórico de los quince años correspondientes al período estudiado (2002 - 2017) indica que el agua subterránea fluye del noreste hacia el oeste.

^{INEGI (2012)} sugiere la existencia de una unidad geohidrológica predominante en la zona, de material no consolidado con posibilidades altas. Ésta se compone de tres unidades geohidrológicas que difieren entre sí según el rendimiento del material: rendimiento bajo (8B, rendimiento menor a los 10 litros por segundo), rendimiento medio (7M, cuyo gasto oscila entre los 10 y los 40 litros por segundo), y rendimiento alto (6A, cuyo rendimiento es mayor a los 40 litros por segundo). A grandes rasgos, se trata de unidades con permeabilidad variable, pero con capacidad para almacenar agua. De acuerdo a esta descripción, se puede interpretar que prácticamente toda la zona de estudio tiene condiciones favorables para la recarga y almacenamiento del agua.

Los suelos que componen esta unidad, en la zona de estudio, son principalmente aluviales en las áreas cercanas al Río Usumacinta, conglomerados en las zonas de mayor altitud, y la presencia esporádica de areniscas. Existen también suelos lacustres cercanos a cuerpos de agua menores. No se encontraron fallas geológicas lo suficientemente cercanas a la zona de estudio que puedan afectar el comportamiento hidrodinámico del acuífero. Se determinaron los caudales de entrada y salida de agua subterránea utilizando valores de transmisividad determinados en la zona de estudio mediante pruebas de bombeo, reportadas por CCT (1991). Aplicando la ecuación de Darcy (ecuación 1) se obtienen los resultados que se presentan en la tabla 1. Como primera aproximación, la entrada del agua subterránea al acuífero “Los Ríos” es de 5.57 x 10^-5m³/s por cada metro de longitud transversal al flujo, y la salida es de 3.03 x 10^-4m³/s por metro. Las redes de flujo constituyen una gran herramienta en el análisis de la interacción agua subterránea - agua superficial. La complejidad de este tipo de estudios radica en que:(i)Monitorear bajo el agua es difícil; la evaluación visual para ayudar a identificar la recarga o descarga es limitada. (ii) La interfase de recarga y descarga de agua subterránea es espacialmente heterogénea; la heterogeneidad geológica es muy difícil de caracterizar, pero puede afectar los flujos de agua a través de ésta. Pueden ocurrir variaciones de un orden de magnitud (o más) en el flujo en la escala de metros. Además, los sistemas de aguas superficiales pueden tener capas finas que ejercen control sobre los flujos verticales y son particularmente difíciles de identificar. (iii) La naturaleza transitoria de la interfaz. Las zonas de recarga y descarga pueden variar con el tiempo, por lo que se necesitan datos continuos para capturar esta variación. Las zonas de recarga pueden convertirse en zonas de descarga y viceversa. Estas variaciones en el tiempo son causadas por diferentes factores, desde la variación natural de la precipitación hasta las alteraciones inducidas por el hombre en el uso de la tierra y el agua.

Un monitoreo de los flujos en la interfase de las aguas subterráneas y superficiales podría proporcionar un indicador de los impactos de la climatología. Es probable que el ciclo hidrológico se intensifique en escenarios de cambio climático al manifestarse tormentas más fuertes y sequías más prolongadas. Estos impactos combinados con posibles cambios en el uso de la tierra podrían tener efectos superpuestos y potencialmente amplificados (^{Juckemet al., 2008}). En la década pasada ^{Carrillo-Sosaet al., (2012)} estudiaron las inundaciones de la región con el objetivo de dar una solución con base en análisis hidráulicos realizados por medio de modelación matemática. Se comprobó que los esfuerzos dedicados a controlar las inundaciones no serían suficientes ya que ninguna de las alternativas, planteadas como medidas de mitigación de los niveles de inundación, las evitarían, por lo que entonces se concluyó rescatar la filosofía de convivencia con el agua (^{Alatriste-Domínguez, 2019}). El cambio climático no representa la única presión sobre los recursos naturales y la abundancia del agua; el sistema hidrológico de la zona de estudio ha estado sujeto a varios procesos de degradación ambiental, una comprensión más profunda de estos temas a escala regional se puede consultar en ^{Andrade-Velázquez y Medrano-Pérez (2020)}.

^{Nooren (2017)} presenta un estudio sobre la influencia del clima y las actividades humanas enfocado a las cuencas de los ríos Usumacinta y Grijalva. Estudios previos han enfatizado el impacto de las sequías prolongadas y su posible relación con el declive social, así como el colapso de la cultura Maya en el preclásico y clásico dando como resultado el abandono de las ciudades (^{Hoggarthet al., 2016}; ^{Medina-Elizaldeet al., 2016}; ^{Lentzet al., 2014}), pero se ha prestado muy poca atención a episodios de eventos de precipitación extraordinarios e inundaciones que también pu- dieron haber impactado la sociedad Maya. Esto puede probarse por el hecho de que las inundaciones, así como las sequías, son un tema importante representado en los antiguos códices mayas (^{Thompson, 1972}); sin embargo, la incertidumbre en el balance hidrológico y las diferentes escalas de los modelos climáticos dificultarían una evalua- ción de los impactos, y por ende, una predicción (^{Loaicigaet al., 1996}; ^{Greenet al., 2011}). ^{Noorenet al., (2017)} plantearon la hipótesis de que el clima durante el Período Preclásico Medio pudo haber sido desfavorable para la intensificación de la agricultura basada en el maíz, que formó la economía de subsistencia responsable del desarrollo de muchas sociedades mesoamericanas vecinas, lo que explica el desarrollo relativamente retrasado de la civilización Maya en esta área.

4. Conclusiones

Históricamente el acuífero de la región de “Los Ríos”, en su porción analizada, no presenta variaciones bruscas en cuanto a su dirección de flujo, la cual va de este - noreste a oeste. A pesar de que existen variaciones en los niveles freáticos, éstos suelen oscilar dentro de los mismos rangos de elevación.

Debido al comportamiento piezométrico en las márgenes del Río Usumacinta, se concluye que el agua subterránea recarga los cuerpos de agua superficiales en la zona.

Los gradientes hidráulicos han sido caracterizados, notándose un paralelismo entre sus zonas con mayor valor y las elevaciones más pronunciadas tanto de los niveles freáticos como del terreno natural. Existen zonas bajas que corresponden a los puntos de egreso del agua tanto superficial como subterránea, hacia el oeste - noroeste.

Las recargas del acuífero en la porción analizada se dan mediante precipitación y sus flujos de ingreso. Las descargas se dan mediante su aporte a los cuerpos de agua, los volúmenes subterráneos egresantes y en menor proporción por actividades antropogénicas.

El análisis de la evolución de los niveles piezométricos permite contar con las herramientas técnicas necesarias para la toma de decisiones en casos de contingencias ambientales, también provee información útil para la reglamentación del uso racional de las aguas subterráneas.

Se sugiere establecer programas de monitoreo que garanticen la cobertura total del área para mejorar la cantidad y calidad de los registros históricos. Es necesario ubicar estaciones piezométricas en las márgenes de los principales cuerpos de agua perennes para poder describir con mayor exactitud la relación que éstos guardan con el acuífero. De igual manera, el establecimiento de estaciones hidrométricas es importante para medir de manera indirecta el caudal base de las corrientes principales en la zona de estudio.