Comportamiento hidrogeoquímico de flujos subterráneos en acuíferos cársticos fracturados, aplicando modelación inversa: Caso Huasteca Potosina

Hydrogeochemical behavior of underground flow in fractured karstic aquifers, by means of inverse modeling: the Huasteca Potosina case

Janete Morán-Ramírez², José Alfredo Ramos-Leal^1,*, Briseida López-Álvarez^2,3, Simón Carranco Lozada², Germán Santacruz-De León³

¹ Instituto Potosino de Investigación Científica y Tecnológica, A.C. (IPICYT), C. Presa San José # 2055, Lomas 4ª sección, C.P. 78216, San Luis Potosí, San Luis Potosí, México. *jalfredo@ipicyt.edu.mx

³ Programa Agua y Sociedad. Colegio de San Luis, A. C. (COLSAN). Parque de Macul No. 155, Fracc. Colinas del Parque, San Luis Potosí, S.L.P. México, C.P. 78299.

Manuscrito recibido: Junio 3, 2011.
Manuscrito corregido recibido: Mayo 5, 2012.
Manuscrito aceptado: Junio 27, 2012.

Resumen

El área de estudio se encuentra dentro de la Plataforma Valles-San Luis Potosí la cual pertenece a la Sierra Madre Oriental (SMO), cuya porción montañosa formada por anticlinales y sinclinales tienen una orientación preferencial NW-S y un sistema muy importante de fallas y fracturas con orientación E-W que controlan el flujo subterráneo hacia el Golfo de México. Para demostrar la conectividad hidráulica del sistema de fallas y fracturas E-W, se aplicó modelación inversa en cuatro secciones, otras pruebas realizadas a posibles conexiones entre fracturas no fueron capaces de reproducir estos procesos, lo que se interpreta como puntos sin conexión hidráulica. Se identificaron dos principales tipos de agua: uno bicarbonatado cálcico debido a interacción con rocas carbonatadas como las calizas y otro sulfatado cálcico que se relaciona con un proceso de interacción agua-roca, principalmente de origen evaporítico como los yesos. En el trayecto de Rioverde a Santa Anita, el agua subterránea disminuye su composición química debido a un proceso de dilución por mezcla con otras aguas de reciente infiltración o menos evolucionadas por lo que sus parámetros fisicoquímicos disminuyen. El flujo subterráneo de Apesco a Huichihuayan muestra un incremento en sus parámetros fisicoquímicos debido a la disolución de yeso, calcedonia y dolomita; así como, precipitación de anhidrita, aragonita, calcita, halita y celestita. La evolución hidrogeoquímica de Unión de Guadalupe a Tambaque muestra un incremento de los parámetros fisicoquímicos por la disolución de calcita, fluorita, halita, sílice, azufre; así como, precipitación de dolomita. En la sección El Salto- Minas Viejas, se relaciona con la disolución de calcita, pirita y cuarzo, junto con la precipitación de celestita, dolomita, FeS, yeso, halita y jarosita.

Palabras clave: modelación hidrogeoquímica, modelación inversa, balance de masas, conectividad hidráulica, índices de saturación, karst, precipitación y disolución.

Abstract

The study area is located inside the Valles-San Luis Potosí platform, which belongs to the Sierra Madre Oriental (SMO), whose mountain portion, formed by anticlines and synclines, has a North-South preferential orientation and a very important system of faults and fractures with East-West orientation that control the groundwater flow into the Gulf of México. In order to demonstrate the hydraulic connectivity of the East-West system of faults and fractures, inverse modeling was applied to four sections. Other performed tests in possible connections between fractures were not able to reproduce these processes, which are interpretated as non-hydraulic connection points. Two main types of water were identified: one with calcium bicarbonate because of the interaction with carbonate rocks such as limestone, and other with calcium sulfate that is related with a process of water-rock interaction, mainly from evaporite origin such as plaster. On the flow path to Rioverde-Santa Anita, groundwater decreases its chemical composition due to a process of dilution by mixing with recent water infiltration or less evolved so their physicochemical parameters decrease. The flow path from Apesco to Huichihuayán shows an increase in their physicochemical parameters due to the dissolution of gypsum, chalcedony and dolomite, as well as precipitation of anhydrite, aragonite, calcite, halite and celestite. The hydrogeochemical evolution of Union de Guadalupe to Tambaque shows an increase in the physicochemical parameters for the dissolution of calcite, fluorite, halite, silica, sulfur, as well as precipitation of dolomite. The Salto-Minas Viejas flow section is related to the dissolution of calcite, pyrite and quartz, accompanied by the precipitation of celestite, dolomite, FeS, gypsum, halite, jarosite.

Keywords: hydrogeochemical modeling, inverse modeling, mass balance, hydraulic connectivity, saturation index, karst, hydrogeochemical facies, precipitation and dissolution.

1. Introducción

La composición química del agua subterránea que no se ha mezclado es controlada principalmente por factores como tiempo de residencia, distancia recorrida, mineralogía y geología del acuífero. La distribución espacial de los parámetros fisicoquímicos en el agua subterránea da una idea de la heterogeneidad y la conectividad del subsuelo por donde circula el agua. Ésta puede fluir por acuíferos en medio granular o en medio fracturado karstificados. En el primer caso el flujo es laminar, en tanto que el flujo subterráneo en medios fracturados se caracteriza por tener velocidades altas y heterogéneas, tiempos de residencia relativamente cortos y gran interacción agua-roca (Roback et al., 2001).

Durante el trayecto del flujo subterráneo ocurren procesos hidrogeoquímicos como el intercambio iónico, reacciones ácido-base, procesos óxido-reducción, disolución, precipitación, posibles mezclas de aguas, reducción bacteriana y degradación.

Una forma de conocer los procesos hidrogeoquímicos es mediante la aplicación de modelos numéricos los cuales proporcionan información detallada acerca de su origen, interacción agua-roca y patrones de flujo (Hidalgo y Cruz-Sanjulián, 2001).

La modelación directa y la modelación inversa son dos tipos de modelaje hidrogeoquímicos en donde se aplican los principios físico-químicos para describir las reacciones que ocurren en el sistema hidrológico (Parkhurst et al., 1980). La primera, predice la composición del agua y la transferencia de masa que puede resultar de reacciones químicas hipotéticas; mientras que la segunda, se define como la transferencia de masas a partir de datos químicos, isotópicos y mineralógicos observados, intentando explicar el significado de un conjunto de datos hidrogeoquímicos e identificar de forma cuantitativa la reacciones hidrogeoquímicas que dan origen a dicha composición (Parkhust y Apello, 1999).

Ambos modelos hidrogeoquímicos utilizan reacciones químicas específicas como intercambio iónico, reacciones ácido-base, procesos óxido-reducción, disolución, precipitación, posibles mezclas de agua; así como reducción bacteriana los cuales pueden ocurrir en el trayecto del flujo subterráneo.

Los cálculos inician con la composición inicial del agua, a partir de la cual se realiza una especiación para obtener los índices de saturación presentes en la solución, su composición determina la distribución de las especies presentes en la solución, actividad termodinámica, así como, el estado de saturación de la solución con respecto a las distintas fases minerales. Finalmente se calcula con códigos numéricos la transferencia de masa entre una fase sólida y la fase acuosa en equilibrio, este proceso de interacción agua-roca es reversible (Plummer y Back, 1980; Apello y Postma, 1996). Como resultado se obtiene la composición del agua, identificación de la roca con la que reacciona en el sistema en equilibrio (evolución composicional) y las solubilidades de los minerales, así como, la transferencia de masa entre distintas fases (Helgeson et al., 1970; Hereford et al., 2007). Una manera simple de ejemplificar este proceso es con la siguiente ecuación.

Agua inicial + reactivos = agua final + productos

Agua inicial es la composición del agua en el punto origen, el reactivo es el medio rocoso con el que interactúa, el agua final y los productos son el resultado de los procesos de interacción con el medio.

La modelación inversa se aplica entre dos puntos, los cuales se infiere que están a lo largo de una línea de flujo. Se supone una composición inicial y una final (puntos extremos). Para la aplicación del modelo se parte de la premisa de que existen diferencias en los componentes químicos de ambos puntos (delta iónico). El modelo cuantifica los procesos que son posibles; para explicar los cambios observados sobre una línea de flujo; además, calcula todos los posibles conjuntos de reacciones entre puntos extremos y la fase mineral definida (André et al., 2005). Cabe mencionar que la conectividad entre dos puntos no siempre es posible.

En el presente trabajo se utiliza la modelación hidrogeoquímica inversa para demostrar la conexión hidráulica existente entre pares de manantiales; como Apesco-Huchihuayan, Unión de Guadalupe-Tambaque, Rioverde-Santa Anita, El Salto-Minas Viejas, localizados en la parte cárstica de la Huasteca Potosina.

La Huasteca Potosina se encuentra en la parte central de la Sierra Madre Oriental; está integrada principalmente por calizas plegadas y fracturadas formando valles intermontanos (Figura 1). El área de estudio se considera parte de la zona de descarga de la Sierra Madre Oriental.

En general se conoce poco sobre el funcionamiento hidrodinámico de los sistemas kársticos en la SMO, debido a que son difíciles de estudiar por las condiciones abruptas del terreno y la excesiva vegetación; así como, a la ausencia de pozos de extracción para su estudio.

El área de estudio también conocida como Plataforma Valles-San Luis Potosí (PV-SLP), pertenece a la Sierra Madre Oriental, cuya porción montañosa formada por anticlinales y sinclinales tienen una orientación preferencial NW-SE.

Las edades de las rocas que afloran en el área de la PV-SLP varían del Mesozoico al Reciente (Figura 2). Las rocas más antiguas dentro del Mesozoico corresponden al Cretácico Inferior, constituido por depósitos evaporíticos, (principalmente yesos de la Formación Guaxcama) cabe mencionar que estas rocas afloran al W del área de estudio.

El Cretácico Inferior se depositó fuera de la PV-SLP, consiste de calizas compactas de textura fina, con nódulos de pedernal perteneciente a la Formación Tamaulipas Inferior de facies batial. Al interior de la PV-SLP se desarrolló en un ambiente evaporítico la Formación Guaxcamá, cuyo espesor es mayor a los 2000 m, como lo indica el pozo Lagunillas, perforado por PEMEX (CRM, 1992). En el Cretácico medio fue depositada la Formación El Abra en un ambiente lagunar con su correspondiente frente arrecifal y en los bordes exteriores se acumularon restos del mismo que dieron origen a la Formación Tamabra, el espesor de la Formación El Abra varía de 1500 - 2000 m, (Carrillo y Suter, 1982). El Cretácico Superior involucra un número más variado de formaciones las cuales cambian dependiendo de la posición que guardan dentro de la PV-SLP, la primera es la Formación Agua Nueva la cual está formada por una secuencia de caliza-lutita, por espesores aproximados de 80m (Padilla y Sánchez, 1985); ésta es cubierta por calizas arcillosas y margas de la Formación San Felipe, su espesor varia de 120 a 300m (Muir, 1936; CRM, 1992). Sobre esta última se encuentra la Formación Méndez compuesta por lutitas y capas delgadas de margas grises y azules, con un espesor no mayor de 70m (CRM, 1992).

Durante la formación de la Sierra Madre Oriental las rocas sedimentarias fueron plegadas originando formas alargadas y en posición paralela (anticlinales y sinclinales). El esfuerzo compresivo en las rocas sedimentarias produjo anticlinales con esfuerzos distensivos en sus crestas, generando un intenso fracturamiento. En algunos casos el plegamiento fue tal que produjo grandes cabalgaduras.

En general, existen dos sistemas estructurales (pliegues, fallas y fracturas) que afectan a las unidades cretácicas (Zapata-Zapata y Pérez-Benzor, 1979); uno es longitudinal que forma pliegues con orientación NW-SE y el otro es perpendicular a los ejes del plegamiento (Pérez-Corona, 2010). Las estructuras kársticas (dolinas, úvalas y poljes) se desarrollan en el sistema de pliegues NW-SE (Pérez-Corona, 2010). Posterior a la generación de este tipo de estructuras se produjeron fracturas y fallas normales con orientación E-W, que sirven de conducto entre los valles intermontanos (Ramos-Leal y Santacruz-De León, 2010).

3. Hidrogeología

A nivel regional el flujo subterráneo, va de las regiones de mayor elevación topográfica, localizadas en la Sierra de Álvarez de la zona media (2423 msnmm), a las de menor elevación hacia el Golfo de México en la zona de la Huasteca (Figuras 1 y 3).

En el área de estudio existen unidades con condiciones hidrogeológicas que funcionan como acuitardos (formaciones San Felipe y Méndez) y otros como acuíferos (formaciones Abra y Agua Nueva), esta última tiene poco espesor y no siempre está presente en la secuencia estratigráfica debido a que fue erosionada.

En general el agua fluye a través de las calizas en los valles y es captada principalmente por fallas regionales que descargan hacia el Golfo de México en la zona de la Huasteca (Figuras 1 y 3). Solo en la parte interna de la PV-SLP el flujo subterráneo interacciona con los yesos de la Formación Guaxcamá y es controlada por estructuras como anticlinales y sinclinales con orientación NW-SE y fallas regionales con una disposición E-W. Por lo que es posible suponer que el sistema de fallas E-W podría conectar la Zona Media con la Zona de la Huasteca (Ramos-Leal y Santacruz-De León, 2010).

4. Metodología

El muestreo fisicoquímico de este estudio se realizó en Agosto de 2008. Durante este periodo se tomaron 28 muestras de manantiales en época de estiaje en el sistema kárstico para determinar cationes y aniones mayoritarios (Figura 1). En época de lluvia se colectaron 5 muestras en sitios ya muestreados, para conocer el efecto de la mezcla meteórica en el sistema. Cabe mencionar que estas representan aproximadamente el 98 % de los manantiales más importantes del área de estudio.

Todas las muestras de agua fueron colectadas en botellas de polietileno, (HDLP), lavadas y enjuagadas 7 veces con agua desionizada. Para colectar las muestras para determinar cationes y elementos traza, las botellas se lavaron además con HCl 10 %. Para cada muestra colectada, se midió in situ el pH, conductividad eléctrica, temperatura, potencial redox (PR), oxígeno disuelto (OD) y alcalinidad.

Inmediatamente después de ser colectadas, las muestras para cationes y elementos traza se acidificaron con ácido nítrico puro hasta pH < 2. Todas las muestras se almacenaron a una temperatura menor a 4 °C. Los iones mayores y elementos traza se analizaron en el laboratorio del Centro de Geociencias de la Universidad Nacional Autónoma de México. Las concentraciones de cationes mayores Ca²⁺, Mg²⁺, Na⁺, y K⁺ fueron determinados por ICPOES y los aniones mayores (SO₄^2- y Cl^-) se analizaron usando cromatografía de líquido. La alcalinidad y HCO₃^- se determinaron en el sitio por titulación. El error de equilibrio iónico (electro neutralidad) fue menor a 5 % (Tabla 1).

Los índices de saturación (SI = logIAP / K, donde IAP es el producto de actividad de iónica y K es el producto de solubilidad) con respecto a la calcita (SICa), aragonita (SIAr), dolomita (SIDo), yeso (SIGy), fluorita (SIF), y anhidrita (SIAn), fueron calculados con PHREEQC (Parkhurst y Apello, 1999).

5. Resultados y Discusión

Las facies hidrogeoquímicas, tambien conocidas como familias de agua, pueden ser identificadas a partir de diferentes métodos. Los iones mayores de las muestras tomadas en la zona de estudio fueron gráficados un diagrama de Piper (Figura 4), método gráfico cualitativo de identificación y además, permite establecer la evolución hidrogeoquímica del agua subterránea.

En la figura 4, se identifican tres principales tipos de agua: 1) bicarbonatada cálcico (Ca-HCO₃) característica de la interacción con rocas carbonatadas, 2) sulfatada cálcica (Ca-SO₄) típico de la disolución de rocas evaporíticas con predominancia en yesos y 3) la transición entre las dos anteriores Ca-HCO₃-SO₄ (Figura 4). La facie Ca-HCO₃, se distribuye principalmente en los valles de la zona de la Huasteca, en tanto que la facie Ca-SO₄ se encuentra en la zona montañosa de la Sierra Madre Oriental.

En el área de estudio ocurren principalmente dos procesos, la mezcla con agua meteórica y la interacción agua-roca (Figura 4). El pH del agua subterránea disminuye en época de lluvia y el agua disuelve más carbonatos, esto se observa con un incremento del índice de saturación (IS) para la dolomita y calcita; mientras que para la época de estiaje se tiene una disminución en solución debido a la precipitación de estos carbonatos. Si la mezcla con agua meteórica fuera predominante se tendrían valores bajos en los componentes fisicoquímicos debido a la dilución y los índices de saturación también serian bajos, pero si además se lleva acabo la interacción agua-roca, se tendrá un incremento en los IS, como en el caso del área de estudio.

En el diagrama de Piper se muestran las 5 muestras colectadas durante la época de lluvia. En la figura 4, se observa la variación química del agua en época de lluvia y estiaje, en el periodo de lluvia, las aguas subterráneas son ligeramente más bicarbonatadas (Ecuación 1). Sin embargo, el proceso de mezcla e interacción agua-roca no es homogéneo en toda la región, debido a que es una zona kárstica y estas características dependen del grado de fracturamiento de cada sitio.

5.2. Modelación Hidrogeoquímica

Los principales flujos subterráneos en el área de estudio van de W – E. Debido a que la zona de estudio está sumamente fracturada, es posible que un número importante de estas fracturas estén comunicadas hidráulicamente, como lo evidencian los manantiales de la región. Para comprobar la conectividad hidráulica del sistema de fallas y fracturas en la Sierra Oriental, se eligieron 4 secciones E-W 1) Apesco-Huchihuayan, 2) Unión de Guadalupe-Tambaque, 3) Rioverde-Santa Anita, 4) El Salto-Minas Viejas, en los que se desarrollaron modelos hidrogeoquímicos para entender esta conectividad en la zona kárstica de la Huasteca Potosina. Las secciones modeladas fueron seleccionadas y ubicadas sobre líneas de flujo que son perpendiculares a las líneas equipotenciales (Figura 1).

En la modelación hidrogeoquímica inversa las principales ecuaciones que gobiernan los procesos de interacción agua-roca en la zona kárstica del área de estudio son:

Disolución de carbonatos:

Disolución de Dolomita y formación de Calcita:

Disolución de Yeso:

Reducción de sulfatos con Materia Orgánica:

Oxidación de pirita:

La modelación hidrogeoquímica se realizó con PHREEQC (Parkhurst et al., 1980), que genera cientos de modelos, de los cuales se eligieron solo tres que cumplen con los criterios de selección, la cual se basa principalmente en la mínima incertidumbre (< 0.05) y la cual se condicionó en el ingreso de datos. También considera las reacciones termodinámicas afines a la geología del las secciones modeladas (marcadas en negritas) de las Tablas 2, 3, 4 y 5; así como, la congruencia de estas reacciones con los deltas iónicos (Figuras 5, 6, 7, 8) de los parámetros fisicoquímicos introducidos al modelo.

En el área de estudio se identificaron principalmente tres flujos subterráneos, uno local, otro intermedio y el tercero regional. El primero esta asociado a la recarga local y el tercero es un flujo más evolucionado con mayor tiempo de residencia y mayor distancia recorrida, entre ambos se encuentra el flujo intermedio (Figura 3 y 4).

La primera sección modelada está influenciada por el flujo local (Apesco–Huichihuayan), la sección Unión de Guadalupe–Tambaque alcanza flujos intermedios, y en ambas secciones el proceso dominante es la mezcla y en las secciones Rioverde-Santa Anita y El Salto-Minas Viejas, prevalece el flujo regional, con interacción agua-roca como proceso dominante. Esto se puede observar en la cantidad de especies involucradas en los modelos seleccionados (Tablas 2, 3, 4 y 5).

5.2.1. Secciones modeladas

1) Apesco-Huichihuayan

El flujo local se produce en calizas kársticas de la Formación El Abra, esta sección se localiza al sur de la zona de la Huasteca, en donde el Manantial de Apesco está en la zona montañosa con una elevación de 866 msnmm y el manantial de Huchihuayan en la zona baja a 149 msnmm. La distancia en línea recta entre ambos sitios es de aproximadamente 7 km. En este trayecto se incrementaron todos los parámetros fisicoquímicos lo cual involucra disolución de minerales (Figura 5). Los manantiales de Apesco y Huichihuayan representan flujos locales bajo condiciones oxidantes.

En la Tabla 2, se muestran los tres modelos principales de transferencia de masa, de acuerdo con los resultados obtenidos, el modelo 2 es el que mejor reproduce los procesos hidrogeoquímicos de la sección.

El modelo 2 predice la formación de calcita (-9.320E-04), y celestita (-5.34E-05); la disolución de anhidrita (3.71E-04), calcedonia (8.51E-06), halita (1.66E-07) y dolomita (5.24E-04). La reacción de sulfato (SO₄^2-) con materia orgánica en medio ácido (H⁺) ó bien en presencia de H⁺, produce CO₂, agua y sulfuro.

En el trayecto el sulfato se reduce con materia orgánica (Ec. 5a) esta reacción puede llevarse a cabo posiblemente debido a la acción bacteriana.

2) Unión de Guadalupe – Tambaque

El modelo 2 (Tabla 3) predice la disolución de calcita, fluorita, halita, sílice, azufre y formación de dolomita. Debido a la disolución de carbonatos, en el agua subterránea se tiene un incremento de calcita (1.40E-03), fluorita (7.38E-06), halita (6.78E-05), sílice (6.50E-05), azufre (8.07E-05). La precipitación de lluvia produce un aumento de dolomita en solución, mientras que en época de estiaje existe disminución de dolomita (-5.76E-05) por la formación del mineral, durante su formación se libera CO₂ (2.84E-03), aunque parte de este CO₂ también puede ser aportado por el carbono orgánico procedente del suelo vegetal.

En la Figura 6 se observan las principales variaciones en pH, conductividad eléctrica (CE), sólidos totales disueltos (STD), Alcalinidad, sulfatos (SO₄) y Ca, entre la muestra que va de Unión de Guadalupe a Tambaque, en todos los casos hubo incremento de estos parámetros fisicoquímicos.

En el trayecto entre Unión de Guadalupe a Tambaque, hay incremento de los iones mayores (SO₄^-2, Cl^-, HCO₃^-, NO₃^-, Ca⁺², Mg⁺², Na⁺, K⁺) y parámetros fisicoquímicos (CE, STD, pH, T, OD y ORP), lo cual indica que la disolución de minerales fue el principal proceso.

3) Rioverde-Santa Anita

La sección 3 se localiza en la Zona Media del área de estudio. Esta sección conecta al manantial de Rioverde a una elevación de 1007 msnmm con Santa Anita de menor elevación (178 msnmm). En la Zona Media el agua subterránea del flujo regional interacciona con formaciones evaporíticas de yesos, en tanto que hacia la zona Huasteca la circulación se desarrolla en calizas cársticas de la Formación El Abra.

En esta sección hay ganancia de CO₂ por la formación de HCO₃^- por lo que aumenta la alcalinidad como se muestra en la Figura 7, donde el flujo va de un tipo de agua sulfatada cálcica (Ca-SO₄) a una bicarbonatada cálcica (Ca-HCO₃). En el trayecto de Rioverde a Santa Anita, hay un incremento de los iones mayores (SO₄^-2, Cl^-, HCO₃^-, NO₃^-, Ca⁺², Mg⁺², Na⁺, K⁺) y los parámetros fisicoquímicos (CE, STD, pH, T, OD y PR), lo cual indica que la disolución de minerales fue el principal proceso que ocurre a lo largo del trayecto del flujo subterráneo.

Para esta sección 3, el modelo 2 (Tabla 3) predice la precipitación de yeso, cuarzo, dolomita, estroncianita, y la disolución de halita. La anhidrita y calcita se encuentran en una zona conocida como metaestable, muy cercana a su precipitación (Ec. 1); por lo cual se encuentra anhidrita en disolución y precipitación de calcita.

Debido a que el yeso se encuentra en una zona metaestable, éste puede disolverse o precipitarse, dependiendo de las condiciones fisicoquímicas (Ec. 3). La presencia de calcio y magnesio en forma de carbonatos favorece la precipitación de dolomita, debido a que ésta tiene un índice de formación mayor que otras fases, como la calcita que se forma en menor proporción.

En la trasferencia se perdió yeso (-7.23E-02), CO₂ (-2.56E-03), cuarzo (-4.83E-05), dolomita (-2.87E-03), estroncianita (-4.68E-05) y halita (-1.72E-04). Por otro lado hubo ganancia de calcita (4.39E-03), anhidrita (6.83E-02).

4) El Salto- Minas Viejas

La sección 4, se localiza al norte del área de estudio, conecta al manantial de Minas Viejas en la zona montañosa con el manantial El Salto en el Valle del Naranjo. En esta sección el agua circula a través de rocas carbonatadas de la Formación El Abra. Aunque en la figura de flujo subterráneo entre el Salto y Minas Viejas muestra poco gradiente, hay evidencia en el flujo superficial de tal conexión hidráulica, debido a que el río que nace en el Salto desaparece cerca de Minas Viejas a través de cavidades de disolución, por esta razón se consideró realizar la modelación hidrogeoquímica.

El modelo 2 (Tabla 4) predice la disolución de calcita, pirita y cuarzo; así como, precipitación de celestita, dolomita, FeS, yeso, halita, jarosita.

En la Figura 8 se observan las principales variaciones en pH, CE, STD, Alcalinidad, SO₄^-2 y Ca⁺², entre la muestra que va de El Salto a Minas Viejas, en todos los casos hubo incremento de estos parámetros fisicoquímicos.

En este sitio es donde se puede observar una mayor precipitación de dolomita, yeso, halita, que al incorporar estos elementos ocurren las reacciones señaladas en las ecuaciones 4 b y c, dando lugar a la formación de pirita.

La disolución de la caliza incorpora Ca⁺² y CO₃^-2 al agua, debido a esto es más alto el contenido de cationes y aniones en el Salto que en Minas Viejas. Debido a la disolución e incorporación al agua subterránea se tiene un incremento en calcita (1.11E+02), pirita 1.11E+02 y cuarzo (4.50E-05). Debido a precipitación, la solución pierde celestita (-1.68E-04), halita (-8.36E-04), dolomita (-8.82E-03), FeS (-1.11E+02), yeso (-1.11E+02), halita (-8.36E-04) y jarosita (-3.51E-04). Durante su formación se libera CO₂ (-1.11E+02).

6. Conclusiones

Con la modelación inversa de 4 secciones con posible conexión hidráulica por medio de fracturas, se logran identificar los procesos hidrogeoquímicos que ocurren en el agua subterránea.

En general la modelación hidrogeoquímica sugiere que el principal proceso que domina en el flujo subterráneo es la interacción agua roca y la mezcla con agua de reciente infiltración.

En el trayecto de Rioverde a Santa Anita, la modelación hidrogeoquímica inversa muestra que en el trayecto ocurrió precipitación de yeso, sílice, dolomita, estroncianita y halita en tanto que se disolvió anhidrita y calcita. El modelo predice la disolución de calcita, pirita y cuarzo; así como, precipitación de celestita, dolomita, FeS, yeso, halita y jarosita.

En el flujo subterráneo de Apesco a Huichihuayan el modelo considera precipitación de anhidrita, aragonita, calcita, halita y celestita; así como, disolución de yeso, calcedonia y dolomita, el proceso gana CO₂ (g) debido a la reacción de SO₄^-2 con materia orgánica en presencia de H.

En la evolución hidrogeoquímica de Unión de Guadalupe a Tambaque la modelación inversa indica la disolución de calcita, fluorita, halita, sílice y azufre; así como precipitación de dolomita.

En la sección El Salto- Minas Viejas, el modelo predice la disolución de calcita, pirita y cuarzo; así como, precipitación de celestita, dolomita, FeS, yeso, halita y jarosita.

La modelación hidrogeoquímica resultó una herramienta económica, útil y eficiente para evidenciar la conexión hidráulica entre cuencas y subcuencas hidrogeológicas e identificar los procesos naturales durante su evolución.

Referencias

André, L., Franceschi, M., Pouchan, P., Atteia, O., 2005, Using geochemical data and modelling to enchance the understanding of groundwater flow in a regional deep aquifer, Aquitaine Basin, south-west of France: Journal of Hydrology, 305, 40-62.         [ Links ]

Apello, C.A., Postma, D., 1996, Geochemistry, groundwater and pollution: Rotterdam, A. A. Balkema, 536p.         [ Links ]

Carranco-Lozada, S.E., 2011, Evaluación de la calidad química del agua subterránea en el acuífero cárstico de San Luis Potosí: México, Guanajuato, Universidad de Guanajuato, Tesis de Licenciatura, 79 p.         [ Links ]

Carrillo, M.M., 1981, Contribución al estudio geológico del Macizo Calcáreo el Doctor, Estado de Querétaro: Universidad Nacional Autónoma de México, Revista del Instituto de Geología, 5 (1), 25-29.         [ Links ]

Consejo de Recursos Minerales (CRM), 1992, Monografía geológico-minera del estado de San Luis Potosí: México, SEMIP, Publicación M-7e, 218 p.         [ Links ]

Helgeson, H.C., Brown, T.H., Nigrini, A., Jones, T.A., 1970, Calculations of mass transfer in geochemical processes involving aqueous solutions: Geochimica et Cosmochimica Acta, 34, 569-592.         [ Links ]

Hereford, G., Keating, H., Guthrie, Jr. D., Chen, Z., 2007, Reactions and reaction rates in the regional aquifer beneath the Pajarito Plateau, north-central New Mexico, USA: Environmental Geology, 52, 965-977.         [ Links ]

Hidalgo, C., Cruz-Sanjulián, J., 2001, Groundwater composition, hydrochemical evolution and mass transfer in a regional detrital aquifer (Baza Basin, Southern Spain): Applied Geochemistry, 16, 745-758.         [ Links ]

Muir, J.M., 1936, Geology of the Tampico Region, Mexico: American Association of Petroleum Geologists, Memoir, 280 p.         [ Links ]

Padilla y Sánchez, R.J., 1985, Las estructuras de la curvatura de Monterrey y estados de Coahuila, Nuevo León, Zacatecas y San Luis Potosí: México, D.F., Universidad Nacional Autónoma de México, Revista del Instituto de Geología, 6, 1-20.         [ Links ]

Parkhurst, D.L., Thorstenson, D.C., Plummer, L.N., 1980, PHREEQE- A Comúter program for geochemical Calculations. U.S. Geological Survey, Water Resources Investigations Report, 80-96.         [ Links ]

Parkhurst D.L., Apello, C.A.J., 1999, User guide to PHREEQC (version 2)-a computer program for speciation, batch reaction, one dimensional transport, and inverse geochemical calculations: U.S. Geological Survey Water Resources Investigations Report, 99–4259.         [ Links ]

Pérez-Corona, F.Y., 2011, Caracterización Geofísica de Ambientes Cársticos, Caso: Polje de Joya de Luna, San Luis Potosí, México: México, San Luis Potosí, Instituto Potosino de Investigación Científica y Tecnológica, A.C., Tesis de Maestría, 197 p.         [ Links ]

Plummer, L.N., Back, W., 1980, The mass balance approach: application to interpreting the chemical evolution of hydrologic systems: American Journal of Science, 280, 130-142.         [ Links ]

Roback, R.C., Johnson, T.M, McLing, T.L., Murrell, M.T., Luo, S., Ku, T-L, 2001, Uranium isotopic evidence for groundwater chemical evolution and flow patterns in the eastern Snake River Plain Aquifer, Idaho: Geological Society of America Bulletin, 113 (9), 1133-1141.         [ Links ]

Ramos-Leal , J.A., Santacruz, G., Torres-Gaytan, D., Morán-Ramírez, J., López-Loera, H., 2010, Caracterización terrenos cársticos por métodos hidrogeoquímicos y geofísicos: Caso Zona de La Huasteca, S.L.P. Unión Geofísica Mexicana, Geos, 29 (1), 38.         [ Links ]

Zapata-Zapata, J.L., Pérez-Benzor, J.A., 1979, Cartografía Hoja Peotillos, S.L.P.: México, San Luis Potosí, Universidad Autónoma de San Luis Potosí, Instituto de Geología y Metalurgia, Folleto Técnico 63, 30p.         [ Links ]