Introducción

Los recursos hídricos en México se caracterizan por asimetrías territoriales y temporales en la disponibilidad natural del agua, con el agravante de que la concentración espacial de la población y la economía se encuentran en las regiones con menor disponibilidad. A diferencia de las regiones norte y centro, el sureste presenta el mayor escurrimiento, dotando así a esta región de una riqueza hídrica privilegiada. Sin embargo, en contraste con el norte y centro, esta región posee menor población y una alta marginalidad (^{CONAGUA, 2018}).

En México, de acuerdo con ^{CONAGUA (2018)} la disponibilidad natural media anual de agua, o agua renovable, al 2017 alcanza los 451,585 hm³/año (incluye las importaciones y exportaciones con países vecinos) que, al dividirse por la población del 2017, se tiene una disponibilidad per cápita de 3,656 m³/hab./año. Asimismo, del total de agua disponible se extraen 87.84 miles hm³ para uso consuntivo (19.5%), donde el 61% es de fuentes superficiales y el 39% de subterráneas, y los usos agrícolas y abastecimiento público son los principales usuarios. Sin embargo, estas cifras no muestran la dependencia de agua subterránea en diversos territorios que, por sus condiciones climáticas y geomorfológicas o por las presiones demográficas y económicas, ven limitados sus recursos superficiales (^{Peñuela-Arévalo & Carrillo-Rivera, 2013}). Así, ^{Spring & Cohen (2012)} indican que, como consecuencia de la concentración espacial del crecimiento de la población, la producción agrícola e industrial, y la distribución natural del recurso (espacial y temporal), la disponibilidad de agua per cápita es cada vez más limitada. Lo que, en efecto, se manifiesta en las distintas regiones donde estos factores convergen, principalmente, centro y norte del país.

No obstante, de ambos recursos, superficial y subterráneo, el agua subterránea es el de menor comprensión y entendimiento en el país (^{Hatch-Kuri, 2017}; ^{López-Morales, 2017}). De hecho, ^{Peñuela-Arévalo & Carrillo-Rivera (2013)} consideran que el agua subterránea no solo es uno de los elementos que mayor impacto negativo ha sufrido desde la última mitad del siglo XIX, sino también que hasta ahora existe poco interés estatal en el estudio y análisis de su dinámica y funcionamiento. En este contexto, México enfrenta una serie de desafíos hídricos como la contaminación, el impacto del cambio climático en el ciclo hidrológico (^{Ojeda-Olivares et al., 2019}), la escasez de agua, el fortalecimiento de la gestión del agua con la participación de los usuarios, promover una buena planificación ambiental, así como impulsar la investigación y la tecnología (^{Cortés et al., 2012}). En particular, en esto último, el uso de técnicas geofísicas puede contribuir al conocimiento y comprensión de las aguas subterráneas del país (^{CONAGUA, 2007}).

Los métodos geofísicos de exploración son esenciales para inferir y estimar la distribución de las propiedades físicas de los materiales del subsuelo. En particular, para este propósito existen los métodos gravimétricos y magnéticos, conocidos como métodos potenciales. Estos se aplican tanto en áreas urbanas como rurales (^{Henderson, 1992}; ^{Campos-Enríquez et al., 1997}; ^{Huizar-Álvarez et al., 1997}; ^{Campos-Enríquez et al., 1997}; ^{Flores-Márquez et al., 1998}; ^{López-Loera, 2014}), o en sectores como la industria petrolera (^{Ndougsa-Mbarga et al., 2007}), la explotación de recursos minerales, la agricultura, geotermia, geología, arqueología, forense y la hidrología (^{Silva, 1974}; ^{Allred et al., 2008}; ^{Robinson et al. 2008}; ^{Reynolds, 2011}; ^{Alatorre-Zamora et al., 2015}; ^{Tinivella et al., 2013}; ^{Binley et al., 2015}; ^{Romero-Ruiz et al. 2018}; ^{Florsch et al., 2019}). De tal modo, la magnetometría es considerada un método versátil y de fácil operación que puede ser aplicado a distintos problemas de exploración del subsuelo (^{Ernstson, 2006}). En particular, entre otras cosas, puede ayudar a determinar regiones con potencial acuífero por medio de la identificación de estructuras geológicas que podrían estar asociadas con la permeabilidad (fracturamiento), como ha sido reportado por ^{López-Loera (2014)}.

La aplicación de los métodos geofísicos a la hidrología no es reciente, la hidrogeofísica, como disciplina de la geofísica para investigar procesos hidrológicos, se ha fortalecido en los últimos años debido a los avances y la disminución de costos en la tecnología e instrumentación (^{Hubbard & Rubin, 2005}; ^{Robinson et al., 2008}). En los últimos años, desde la hidrogeofísica, el uso de métodos geofísicos ha estado enfocado en investigaciones hidrogeológicas, en el conocimiento de problemas ambientales e infiltración de contaminantes en el suelo, en el mapeo de las características del subsuelo, en la variabilidad espacial y temporal de la humedad del suelo, en el funcionamiento biogeoquímico del sistema agua-suelo, en el estudio de procesos hidrológicos, y en las interacciones de las aguas subterráneas con las superficiales y costeras (^{Hubbard & Rubin, 2005}; ^{Vereecken et al., 2006}; ^{Robinson et al., 2008}; ^{Kirsch, 2009}; ^{Binley et al., 2015}; ^{Brocca et al., 2017}; ^{Romero-Ruiz et al., 2018}; ^{Florsch et al., 2019}).

No obstante, si bien el uso de métodos eléctricos y electromagnéticos en la exploración del agua subterránea (i.e., ^{Belmonte-Jiménez et al., 2012}; ^{Belmonte-Jiménez et al., 2014}) está más extendido en relación a los sísmicos, gravimétricos y magnéticos, el costo, área de estudio y la facilidad tecnológica y operativa, son factores a evaluar para la selección del método a usar (^{Silva, 1974}; ^{Randall-Roberts, 2000}; ^{Ernstson, 2006}; ^{López-Loera & Tristán-González, 2013}; ^{Guevara-Mansilla et al., 2017}).

En el caso de México, destacan los estudios para el levantamiento magnético del lago Chapala, Guadalajara (^{Alatorre-Zamora et al., 2015}); en la caracterización geoeléctrica de un sitio contaminado por hidrocarburos en Tabasco (^{Shevnin et al., 2012}); el método electromagnético para caracterizar el acuífero de la Cuenca de México (^{Arango-Galván et al., 2016}) y el acuífero del altiplano de Tula, Tamaulipas (^{Guevara-Mansilla et al., 2017}); la aplicación de métodos magnetométricos y eléctricos para la localización de agua subterránea en ambientes volcánicos áridos de San Luis Potosí (^{López-Loera, 2014}); la potencialidad acuífera del Graben de Villa de Reyes, San Luis Potosí (^{López-Loera & Tristán-González, 2013}); el estudio hidrogeológico de la cuenca de Villa Hidalgo, San Luis Potosí (^{Pérez-Corona et al., 2017}); y el uso de sondeos eléctricos en la identificación de pozos de infiltración para la mitigación de inundaciones en El Barreal, Chihuahua (^{Dena et al., 2013}). Asimismo, destaca el estudio de ^{CONAGUA (2007)} sobre el uso de los métodos potenciales en la prospección hidrogeológica e ingeniería de construcción de obras hidráulicas en México. De tal manera, existen distintos trabajos que han utilizado los métodos geofísicos para la exploración y entendimiento de la dinámica y funcionamiento de las aguas subterráneas y problemáticas asociadas en el país.

Desde esta perspectiva, el objetivo principal de este manuscrito es la identificación de zonas con potencial acuífero en la cuenca baja del río Grijalva mediante el análisis magnetométrico. En tal sentido, para responder a este objetivo, el presente trabajo está organizado de la siguiente manera. En primer lugar, se muestra una descripción de la zona de estudio. En segundo, se explica la metodología usada para identificar las áreas con potencial acuífero. En tercer lugar, se presentan los resultados y la discusión. Y finalmente, la conclusión. Sobre esta base se concluye que el estudio realizado puede favorecer no sólo la investigación en estos temas, sino también la reflexión y la toma de decisión de los distintos actores en torno a la protección, gestión, aprovechamiento y planificación de las aguas subterráneas en el país.

Área de estudio

El área estudiada se encuentra dentro de la Región Hidrológica Grijalva-Usumacinta Núm. 30 (RH-30), la cual se ubica en la latitud norte 18° 0.106’ y 17° 18.122’ y longitud oeste 92° 52.205’ y 93° 34.427’. El río principal es el Río Grijalva, o Mezcalapa, como es conocido en esta porción de la cuenca Grijalva (^{Álvarez-Soberano & Medrano-Pérez, 2020}). Administrativamente, se localiza en el sector norte del estado de Chiapas cubriendo pequeñas porciones de Tabasco en su parte norte (Figura 1a). Está delimitada por las subcuencas Río Viejo Mezcalapa, Río Mezcalapa, Río Tzimbac, Río Zayula, Río Platanar, Río Paredón, Río Pichucalco (Figura 1b). Posee una superficie de 381.97 km², y alberga una población de 340,160 habitantes, según datos del año 2015 del Instituto Nacional de Estadística y Geografía (INEGI). Asimismo, el área de estudio cubre en su totalidad a la delimitación administrativa del acuífero Reforma, y parcialmente, a los acuíferos Samaria-Cunduacán, La Sierra, Huimanguillo, Ocozocoautla (Figura 1c).

Figura 1 Área de estudio: a) Ubicación; b) Subcuencas; c) Acuíferos; d) Elevación, y e) Panorámica del volcán El Chichón (o Chichonal). 

El área de estudio se caracteriza por valores de elevación de entre 1 y 2,425 m s. n. m., aunque gran parte de la zona puede considerarse de baja elevación (Figura 1d). En cuanto la geomorfología de la zona de estudio comprende una porción de la Planicie Costera del Golfo en su parte norte, encontrado relieves suaves con lomeríos en cadena orientados preferencialmente en dirección noreste-suroeste y elevaciones que alcanzan poco más de los 300 m s. n. m. La porción sur del área de estudio se caracteriza por una serie de sierras y valles. Esta zona está constituida por las denominadas Montañas del Norte, formando una franja con orientación este-oeste y altitudes que rebasan los 1,500 m s. n. m. (^{DOF, 20/04/2015}). Dentro de esta porción, a unos 20 km al suroeste de Pichucalco, se encuentra el volcán Chichonal y sus rocas volcánicas asociadas (Figura 1e) (^{Plascencia et al., 2014}).

En tanto que la geología de la zona, al cubrir el acuífero Reforma prácticamente la totalidad del área estudiada, las características geológicas de este acuífero pueden ser representativas para el área estudiada. En la porción sur del área, correspondiente al límite norte de la Sierra de Chiapas, se encuentran aflorando rocas del Cretácico Inferior y Superior surcadas por una serie de fallas y plegamientos, así como sedimentos Terciarios y rocas ígneas extrusivas. Tienen amplia distribución las rocas calizas y les siguen lutitas y areniscas muy intemperizadas, como consecuencia del clima cálido y las altas precipitaciones de esta zona. El resto del área corresponde a depósitos Cuaternarios y a sedimentos miocénicos. La parte principal del acuífero se localiza en la Planicie Costera constituida por depósitos granulares finos con espesores de alrededor de 30 m, cubriendo los sedimentos del Terciario (lutitas y areniscas) que constituyen la Planicie Costera del Golfo. Esta zona se encuentra configurada por una terraza aluvial del Pleistoceno, actualmente fragmentada y muy desgastada, que se eleva sólo pocos metros arriba de la llanura aluvial del Reciente (el Holoceno) e integrada por una composición de canales de ríos activos e inactivos, diques naturales y pantanos secundarios (^{Plascencia et al., 2014}; ^{DOF, 20/04/2015}). En general, se muestra la Figura 2 una caracterización geolitológica de mayor detalle sobre el área estudiada.

Figura 2 Características geolitológicas del área de estudio: a) tipología de rocas y b) litología. 

Metodología

Para cumplir con el objetivo propuesto para esta investigación, se llevaron a cabo varias etapas fundamentales (Figura 3). En primer lugar, se realizó una recopilación de información relacionada con la geología de la zona de estudio. Posteriormente, se obtuvo la información aeromagnética de las 9 cartas magnéticas a escala 1:50,000 en el Servicio Geológico Mexicano (SGM) en formato PDF (Pichucalco E15-C29, Reforma E15-C19, Raudales Malpaso E15-C48, Villahermosa E15-D11, Huimanguillo E15-C18, Rómulo Calzada E15-C38, Estación Chontalpa E15-C28, Ixhuatán E15-C39, Tecpatán E15-C49), a las cuales se les aplicaron tratamientos para mejorar la calidad de las imágenes para su posterior digitalización y la delimitación del área de estudio en el sistema de información geográfica ArcMap.

Figura 3 Metodología de identificación de áreas con potencial acuífero. Nota: Servicio Geológico Mexicano (SGM); Anomalía Magnética de Campo Total (AMCT) y TILT (Derivada TILT). 

Una vez concluido el proceso de digitalización, se continuó con el preprocesado, el cual consistió en la agrupación de los datos correspondientes a las nueve cartas magnéticas. A continuación, se generó una matriz de datos inicial o grid (archivo GRD) que posteriormente se cargó al software Oasis Montaj, generando así el mapa de Anomalía Magnética de Campo Total (AMCT). Consecuentemente, se aplicaron diversos filtros o algoritmos matemáticos considerados en la metodología para el procesado de los datos, tal como se muestra en la Tabla 1. Derivado de ello, se obtuvieron los distintos mapas correspondientes a los filtros aplicados, los cuales fueron Anomalía Magnética del Campo Total (AMCT), Señal Analítica y Derivada TILT.

Resultados y discusión

En el mapa de la anomalía magnética de campo Total (AMCT) mostrado en la Figura 4a, se observa una anomalía con comportamiento dipolar de gran longitud de onda (dimensiones mayores) que se extiende sobre el área de estudio. La parte con valores de tendencia positiva se ubica al norte, y la parte de valores negativos al sur. La dirección principal de las anomalías presentes en esta configuración es NE-SW. Según la distribución de subcuencas hidrográficas en la zona, los bajos magnéticos están enmarcados en ciertas porciones de las subcuencas R. Zayula, R. Pichucalco, R. Tzimbac y R. Platanar, al sureste de la zona de estudio. Cabe señalar que esta zona está asociada al vulcanismo reciente producido por el volcán Chichonal, por lo que las respuestas magnéticas presentes se ven fuertemente afectadas, debido a que las anomalías termales en profundidad afectan la magnetización de las rocas y, por lo tanto, a las anomalías magnéticas observadas en superficie (^{Carrillo-de la Cruz et al., 2020}).

Figura 4 a) Configuración de anomalía magnética del campo total (AMCT). b) Configuración de Señal Analítica de la AMCT. c) Configuración de derivada TILT de la AMCT. 

En la zona norte, predomina una anomalía positiva de gran longitud de onda. Esta zona coincide con la porción correspondiente a la Planicie Costera del Golfo, donde encontramos los valores altimétricos más bajos. Superpuestas a estas dos anomalías mayores tenemos una serie de anomalías de dimensiones menores. En la porción central, entre los valores positivos mencionados anteriormente, tenemos una zona caracterizada por un fuerte gradiente con contrastes de valores magnéticos de menor extensión que guardan una dirección preferencial noroeste-sureste (NE-SW). Este comportamiento magnético coincide con una zona de transición entre materiales de cobertura sedimentaría y productos volcánicos. Es en esta región en donde la altimetría del relieve comienza a descender.

Asimismo, la configuración de la señal analítica de la AMCT (Figura 4b), al igual que la derivada TILT, es una herramienta que delimita gradientes, ayudando a establecer los límites de los cuerpos que originan las anomalías (i.e., ^{Pérez-Corona et al., 2017}). De misma forma como en la AMCT, la dirección principal de las anomalías presentes en esta configuración es NE-SW. El primer alto magnético está ubicado al norte de la zona de estudio, manteniendo el comportamiento suave como se ha mencionado anteriormente. El bajo magnético presente al sureste en la AMCT está ahora representado por valores positivos intermedios. Estas dos anomalías mayores son separadas por un trend de bajos magnéticos alienados con dirección suroeste-noreste (NE-SW). Respecto al bajo y alto magnético de la gran anomalía dipolar asociada al volcán Chichonal, estos se encuentran claramente definidos por la señal analítica.

La configuración de la derivada TILT de la AMCT, mostrada en la Figura 4c, indica una serie de lineamientos distribuidos en toda la zona de interés, asociados con fuentes superficiales que en un principio estaban superpuestas a las anomalías de carácter regional, por tal razón es posible observar una mayor cantidad de fuentes magnéticas de menor tamaño. Al igual que en las configuraciones anteriores, la dirección principal de las anomalías presentes en esta configuración es NE-SW. Nótese que las repuestas magnéticas muestran tendencias lineales, por lo que pueden representar características estructurales de interés hídrico (zonas de fallas y/o fracturas), puesto que se trata de un algoritmo matemático útil para la detección de contactos entre cuerpos magnéticos, debido a su naturaleza. Respecto a la zona sur, se observan una serie de altos magnéticos con dirección preferencial NE-SW distribuidos por toda la región, los cuales no presentan un interés hídrico debido a que estas respuestas coinciden con la zona con relieves más altos que corresponden a las denominadas Montañas del Norte. Estos rasgos se mantienen presentes hasta la parte central de la zona estudiada, punto a partir del cual la altimetría del relieve comienza a descender. Asimismo, se muestra un alto magnético de dimensiones considerables en la zona norte con dirección NE-SW, sobre la cual se encuentran superpuestas una serie de anomalías de dimensiones menores. Esta zona corresponde al área de mayor interés o potencial acuífero debido a que se encuentra en la porción de menores valores altimétricos, estando influenciada por los escurrimientos de agua provenientes de las partes altas en la zona estudiada. Como se observa en la Figura 4c, la derivada TILT también enfatiza los límites del alto magnético de la anomalía dipolar asociada con el volcán Chichonal.

Según ^{Padilla y Sánchez (2007)}, las rocas de la zona de estudio se encuentran afectada por una serie de rasgos estructurarles (pliegues y fallas) producto de los procesos geológicos que configuraron la geología de la zona de estudio y sus alrededores. Estos rasgos estructurales distribuidos tienen una dirección preferencial NE-SW. Como ya me menciono, las anomalías presentes en las 3 configuraciones mostradas en la Figura 4 presentan principalmente direcciones NE-SW, mostrando una fuerte relación con los rasgos estructurales propuestos por ^{Padilla y Sánchez (2007)}. En concreto, los valores de tendencia positiva distribuidos particularmente al norte de la zona de estudio se encuentran presentes en las tres configuraciones obtenidas (Figura 4). En este sentido, el alto magnético sobresaliente en todas las configuraciones podría estar asociado con un anticlinal.

En relación con la construcción del modelo magnético en la zona estudiada, se aplicó la Deconvolución de Euler (DE) para los índices estructurales 0.5, 1.0, 1.5, 2.0, 2.5 y 3.0 utilizado una ventana de 9x9 datos, y así obtener los rasgos de profundidad asociados a cambios en la geometría de la distribución de los materiales magnéticos que constituyen la base para la construcción del modelo magnético (Figura 5). Todas las soluciones obtenidas con la DE fueron integradas a una misma base de datos. Posteriormente se interpolaron estos valores utilizando el método de vecinos cercanos en el software Voxler 4.1, esto con el objetivo de integrar la respuesta de una fuente magnética equivalente asociada a la capa de roca magnética, que constituye el basamento cristalino impermeable. El uso de un solo índice estructural solo proyecta una parte de la estructura, así que al interpolar las soluciones se genera una imagen que muestra la distribución de estas fuentes magnéticas (Figura 5).

Figura 5 Visualización 3D de las soluciones de la Deconvolución de Euler para los índices estructurales. Nota: 0.5 (Azul marino), 1.0 (Azul celeste), 1.5 (Verde), 2.0 (Amarillo), 2.5 (Rojo) y 3.0 (Morado). Las zonas en blanco en el área estudiada corresponden a las porciones en donde no tienen soluciones de Euler. En la porción inferior izquierda se muestra la vista en planta de los resultados. 

Una vez interpoladas las soluciones, se realizó un proceso de asociación o conexión de fuentes para ayudar visualmente a resaltar la presencia de materiales magnéticos asociados al basamento, contrayendo así un modelo tridimensional (Figura 6, 7a). En las Figuras 7b y 8a, las áreas iluminadas en color azul corresponden a zonas en donde no se encontró una fuente magnética equivalente dentro de la dispersión en la interpolación. Estas zonas son asociadas a espesores de paquetes sedimentarios presentes en la zona de estudio y representan las posibles zonas con mayor potencial para el almacenamiento de fluidos. Los colores oscuros se asocian a la delimitación de lo que suponemos es el basamento magnético. Cabe señalar que estos espesores alcanzan profundidades mayores a los 3,000 metros en esta zona. Asimismo, la Figura 8b corresponde a la superposición de las zonas constituidas por espesores sedimentarios junto con los escurrimientos de las regiones topográficamente más altas, pudiéndose apreciar una buena correlación entre estos dos. En este sentido, se puede observar cómo los escurrimientos provenientes de las partes más altas (sur de la zona estudiada) desembocan en las zonas de espesor sedimentario más profundo (centro y norte de la zona estudiada). Por tanto, se puede inferir, con fundamento en la altimetría y la magnetometría, que la parte central-norte del área estudiada alberga zonas de potencial acuífero.

Figura 6 Visualización 3D del modelo del basamento obtenido a través de la interpolación de las soluciones de Euler. Nota: En la porción inferior izquierda se muestra la vista en planta de los resultados. 

Figura 7 a) Modelo de basamento correspondiente a la zona de estudio (obtenido al interpolar los índices estructurales); b) Zonas de mayor espesor sedimentario (obtenido al interpolar los índices estructurales). Nota: Las barras de colores para a) y b) representan los valores de los índices estructurales interpolados. 

Figura 8 a) Superposición del modelo de basamento y las soluciones correspondientes a las zonas de mayor espesor sedimentario; b) Superposición de las zonas de mayor espesor sedimentario y los escurrimientos potenciales provenientes de las partes más altas [m]. 

Esta zona presenta una compleja evolución sedimentológica y tectónica, pasando de facies carbonatadas del Jurásico-Cretácico a secuencias siliclásticas del Terciario con una gran aportación de sedimentos arenosos con espesores que alcanzan hasta los 6,000 metros hacia la zona costera y que son derivados de la erosión de las zonas continentales altas (^{Velázquez-Aguirre & Ordaz-Ayala, 1993}; ^{Padilla y Sánchez, 2007}). El promedio de la elevación del nivel estático, con base en la red piezométrica de la Comisión Nacional del Agua (CONAGUA), es de 12 m s. n. m. y corresponde a flujos intermedios. Interactuando con la superficie se identifican flujos subterráneos locales en zonas focalizadas de altos topográficos que funcionan como zonas locales de infiltración. Se considera un acuífero libre de buena permeabilidad. En un registro geofísico de un pozo en Cunduacán, Tabasco se interpreta que el acuífero presenta buenas propiedades en una secuencia alternante de arenas, arcillas y limos, donde a partir de los 150-250 metros de profundidad registrados existen capas de arena gruesa, arcillas, y arena media con diferentes espesores (^{Aurora Geofísica, 2016}).

En este contexto, la zona de estudio se encuentra ubicada principalmente en la provincia hidrogeológica Sierras y Valles del Sureste, donde predominan altos niveles de precipitación y riqueza hídrica superficial de abundantes caudales (Grijalva y Usumacinta), siendo las zonas bajas las receptoras de los escurrimientos superficiales y del almacenamiento subterráneo proveniente de las zonas altas (^{Velázquez-Aguirre & Ordaz-Ayala,1993}). A todo ello, respecto a la litología mostrada previamente, el área de estudio se caracteriza por materiales permeables y de alta infiltración, destacando la arenisca-lutita en la parte alta de la zona chiapaneca y una llanura de materiales aluviales en la parte baja del estado de Tabasco.

Con base en los resultados mostrados se evidencia que el método magnetométrico identifica las zonas potenciales acuíferas regionales, por lo que para una mejor caracterización de los acuíferos se recomienda se realicen otros estudios geofísicos como gravimetría y el método transitorio electromagnético para poder corroborar los espesores sedimentarios y el potencial hídrico en la zona, así como una testificación de las zonas con notables características de almacenamiento para un control estratigráfico y para fines de explotación sostenible del agua subterránea de los acuíferos regionales de la zona.

Adicionalmente, cabe señalar que, según datos a febrero ^{2019 del Registro Público de Derechos de Agua (REPDA)} de la CONAGUA, en la zona de estudio se tienen otorgadas 273 concesiones de agua con un volumen total de 62.2 hectómetros cúbicos/año, concentradas fundamentalmente en la agricultura (34%), industria (41%), y tan solo un 15.4% en abastecimiento doméstico. Asimismo, estas concesiones se encuentran ubicadas en el área media-baja de la zona de estudio (Figura 9), por lo que los resultados de la presente investigación confirman la idoneidad de la metodología usada para el análisis de las aguas subterráneas, así como para la identificación preliminar de áreas y estructuras geológicas que podrían albergar un potencial acuífero a nivel cuenca.

Figura 9 Concesiones de agua en el área de estudio (REPDA-CONAGUA) y las zonas de mayor espesor sedimentario. 

Desde esa perspectiva, y teniendo en cuenta los resultados obtenidos, el uso de la magnetometría puede ser clave en la delineación de los rasgos estructurales de las rocas (fallas, intrusiones o la morfología) que configuran a los acuíferos. No obstante, la aplicación geohidrológica de las técnicas magnetométricas no se encuentran muy extendidas en México, fundamentalmente por dos razones: 1) Se trata de un método auxiliar o complementario a otros más determinativos para la detección de rocas que albergan agua y 2) Por el mayor grado de difusión que poseen otros métodos geofísicos, como los eléctricos, en los estudios geohidrológicos (^{CONAGUA, 2007}). Por tanto, más allá de las diferencias entre los métodos geofísicos, la fortaleza de los métodos magnetométricos radica en que permiten analizar grandes superficies, y, por ende, reducir el costo para el uso de sondeos eléctricos verticales (SEV’s) en áreas específicas, teniendo así un conocimiento previo sobre las unidades de rocas que contienen acuífero en la zona estudiada. De tal manera, podrían utilizarse como método complementario en el análisis y conocimiento de las aguas subterráneas.

En tal sentido, los resultados obtenidos en la presente investigación puede representar una herramienta preliminar y de bajo coste para el análisis y conocimiento de las aguas subterráneas, así como para la toma de decisión en torno al ordenamiento, protección y manejo territorial de áreas específicas de las cuencas hidrológicas, con la finalidad, entre otras cosas, de garantizar eventuales o futuras fuentes para el abastecimiento de agua ante la ocurrencia de eventos climáticos extremos y contaminación de las actuales, respectivamente. Por tanto, los resultados de este trabajo se enmarcan como una contribución en esa dirección.

Conclusión

Los resultados muestran las características y rasgos sobresalientes de las principales fuentes magnéticas que fueron interpretados de las diferentes configuraciones obtenidas al aplicar a los datos los algoritmos matemáticos propuestos. Aunado a lo anterior, el análisis del modelo 3D del basamento magnético construido, confirman que el método magnetométrico representa una herramienta factible para la identificación preliminar de zonas constituidas por unidades rocosas con posible potencial acuífero, siendo este el caso para la zona considerada en esta investigación. En la zona de estudio cabe destacar dos puntos importantes:

No obstante, para un mayor detalle y conocimiento de la profundidad de la zona de estudio se requiere la aplicación complementaria de métodos más resolutivos que indiquen la presencia del recurso hídrico, como los sondeos eléctricos verticales (SEV’s). Por tanto, la metodología aplicada podría reducir los costos asociados a la planeación de campañas de exploración por métodos directos e indirectos. En tal sentido, la presente investigación es un aporte al conocimiento y comprensión de las aguas subterráneas, y a su vez, una contribución para la toma de decisiones en los temas geohidrológicos, ordenamiento territorial, protección y conservación del recurso agua, y las políticas públicas asociadas.