Geoquímica de los sistemas de flujo de agua subterránea en rocas sedimentarias y rocas volcanogénicas de Loreto, BCS, México

Geochemistry of groundwater flux systems in sedimentary and volcanogenic rocks at Loreto, BCS, Mexico

Antalia González-Abraham^1*, Juan Reynerio Fagundo-Castillo², José Joel Carrillo-Rivera³, Ricardo Rodríguez-Estrella¹

¹ Centro de Investigaciones Biológicas del Noroeste S.C., Mar Bermejo No. 195, Col. Playa Palo de Santa Rita Apdo. Postal 128, CP 23090, La Paz, B. C. S., México. *agonzalez@cibnor.mx

³ Instituto de Geografía, Universidad Nacional Autónoma de México Ciudad Universitaria, Coyoacán, 04510, México D. F.

Manuscrito recibido: Junio 15, 2010.
Manuscrito corregido recibido: Septiembre 27, 2010.
Manuscrito aceptado: Diciembre 20, 2010.

Resumen

La composición química del agua subterránea es el resultado de continuos procesos de interacción entre el agua de precipitación (y nieve o granizo), que se infiltra en el terreno, y los minerales de la litología por donde circula. El objetivo del presente estudio es determinar el origen de la composición química del agua subterránea del área de Loreto, Baja California Sur. Para ello se utilizaron modelos hidrogeoquímicos de balance de masa y mezcla de aguas, así como otras herramientas de la hidrogeoquímica. Los procesos relacionados se manifiestan en un aumento en la concentración Cl^-, Br^-, I^-, As_total, B^-, Li y Sr. Se identificaron 7 grupos de agua que corresponden a manantiales someros, pozos, manantiales termales y manantiales submarinos. Los patrones hidrogeoquímicos principales obtenidos se relacionan con flujos específicos de agua subterránea que son típicos de procesos de disolución de plagioclasa del tipo albita, andesina y oligoclasa. Los patrones restantes resultado de la mezcla de agua de diferentes flujos, muestran influencia de procesos de intercambio iónico. De esta manera, la composición química del agua que circula en esta área a través de rocas volcánicas félsicas y sus sedimentos derivados debe su origen a la meteorización de los mismos. La mezcla de agua y los procesos de interacción agua-roca se dan entre flujos de tipo local y otros de jerarquía intermedia o regional. Además, la influencia del clima árido y el intensivo desarrollo agrícola, juegan también un papel importante en el modo en que el agua subterránea adquiere su composición química en el área.

Palabras clave: Patrones hidrogeoquímicos, química de fluidos, zonas áridas, Loreto.

Abstract

The chemical composition of groundwater is the result of continuous water-rock interaction between rain water (and hail or snow) that filters into the ground and the minerals in the lithology it flows through. The objective of this study is to determine the origin of the chemical composition of groundwater in the Loreto area, Baja California Sur. Hydrogeochemical mass balance and water mixture models as well as other hydrogeochemical tools were used. All related processes are marked by an increase in the concentrations of Cl^-, Br^-, I^-, As_total, B^-, Li and Sr. We identified seven groups of water corresponding to shallow springs, wells, thermal springs and submarine springs. The hydrogeochemical patterns of the underground water flows (defined numerically and by graphic expressions of water chemistry) are related to the prevailing lithology and are typical of albite, andesine and oligoclase dissolution. The remaining patterns result from the mixture of different groundwater flows, and show influence of ionic exchange processes. Thus, the chemical composition of groundwater circulating through the felsic volcanic rocks and their derived sediments has its origin in the weathering of those materials. Processes of water-rock interaction and water mixture take place among local flows and others of intermediate or regional hierarchy. Furthermore, the influences of the arid climate and intensive agriculture development also play an important role in how groundwater in the area acquires its chemical composition.

Keywords: Hydrogeochemical patterns, fluid chemistry, arid zones, Loreto.

1. Introducción

La composición química del agua subterránea es el resultado de continuos procesos de interacción entre el agua de precipitación (y nieve o granizo) que se infiltra en el terreno y la litología por donde circula. Parte de los componentes químicos son adquiridos en la zona de recarga, mientras que otros se adquieren en el largo recorrido del flujo del agua subterránea hasta que ésta es captada en pozos o emerge en la zona de descarga a través de un manantial. En ese recorrido se puede originar un flujo de diferente naturaleza, de carácter local, intermedio o regional (Tóth, 2000), lo cual depende de la trayectoria de recorrido del agua subterránea en términos de su longitud y profundidad de recorrido y el tiempo de residencia del agua en el subsuelo. El anión predominante evoluciona, por lo general, desde facies hidrogeoquímicas de tipo HCO₃^-, distintiva en flujos de corto recorrido (local), pasando por facies de tipo SO₄^2-, típica en flujos de mayor recorrido (intermedio), hasta facies Cl^-, característica en la descarga de flujos de muy largo recorrido (los regionales) (Tóth, 2000). El tipo de catión depende de la litología (en el sentido estricto, de la mineralogía) atravesada por el agua: Ca²⁺, en carbonatos donde predominan la calcita o plagioclasas cálcicas, ricas en anortita; Ca²⁺> Mg²⁺, en carbonatos de tipo dolomita y en rocas sedimentarias e ígneas, que contienen minerales con abundante Ca²⁺ y Mg²⁺, Mg²⁺> Ca²⁺, fundamentalmente en alumino-silicatos de tipo básico o ultrabásico. De esta manera, la presente investigación se realizó con base en la colecta de muestras de agua del área de estudio que incorpora dos tipos principales de evidencia de agua subterránea, pozos y manantiales perennes.

El estudio del origen de la composición química del agua, producto de diferentes procesos, ha sido abordado desde distintas perspectivas. Por ejemplo, el origen de la composición del agua de manantiales efímeros y perennes, resultado de procesos de meteorización, ha sido explicado aplicando modelos de balance de masa, por ejemplo en Sierra Nevada, E.U.A., (Garrels y Mackenzie, 1967). En este caso se reconstruyeron los principales minerales originales constitutivos de las rocas ígneas encontradas en campo (plagioclasa, andesina, biotita, K-feldespato, montmorillonita). Estos autores hicieron extensiva esta aproximación para explicar el origen de la composición del agua subterránea que circula en otras rocas ígneas, identificando procesos de intemperismo en que la plagioclasa albita originaba Na⁺ y H₄SiO₄ por disolución incongruente de caolinita y esmectita (Na-beideita). Por otro lado, el origen de la composición del agua de algunos lagos de Sierra Nevada fue analizada mediante un modelo de evaporación, donde fueron identificados procesos de precipitación de yeso, calcita, hidroxisilicatos de Mg²⁺ y SiO₂ amorfo (Garrels y Mackenzie, 1967).

El origen de la composición química del agua puede determinarse mediante modelos hidrogeoquímicos inversos (que es una aproximación de balance de masa), y por medio de modelos de mezcla (i.e., de aguas de diferente calidad). Los primeros se basan en la transferencia de masa a partir de datos químicos, isotópicos y mineralógicos observados. Los mismos fueron desarrollados como una estrategia determinista para interpretar la composición química del agua subterránea, a partir de identificar cuantitativamente las reacciones geoquímicas que dan origen a esa composición. Entre estos modelos se pueden citar BALANCE (Parkhurst et al, 1982) y MODELAGUA (Fagundo-Sierra et al, 2001).

En forma similar, el origen de la composición química del agua que circula por rocas sedimentarias y volcánicas (caliza, caliza dolomitizada, calcarenita, arenisca y serpentinita) se ha determinado utilizando modelos de balance de masa y mezcla de aguas (Fagundo-Sierra et al., 2001) así como modelos hidrogeoquímicos basados en equilibrio químico en la Sierra del Rosario, Cuba (Fagundo et al, 2006; ver recopilación en Fagundo et al., 2008). Los principales procesos de intemperismo, en este caso, se producen por disolución congruente de calcita, dolomita y serpentinita, y por disolución incongruente de la plagioclasa albita y microclina. En el agua que circula en medio cárstico costero, contiguo a un humedal de la Cuenca Sur de La Habana, se identificaron, principalmente, procesos de disolución y precipitación de calcita, reducción anaeróbica de sulfato e intercambio iónico directo e inverso entre Na⁺ y Ca²⁺ así como dolomitización (Fagundo et al., 2004). En estos trabajos se han utilizado, además, las ventajas que ofrece la teoría de los patrones hidrogeoquímicos (Fagundo, 1998) con el fin de interpretar la relación entre la composición química del agua y el medio geológico por el que ésta ha circulado.

El objetivo del presente trabajo es establecer el origen de la composición química del agua, utilizando esencialmente modelos que no han sido aplicados en el área de estudio, como los hidrogeoquímicos de balance de masa, mezcla de aguas y reconocimiento de patrones hidrogeoquímicos, así como otras herramientas hidrogeoquímicas que facilitan el entendimiento de la vinculación entre los sistema de flujo de agua subterránea, la litología por donde circula y el efecto de la actividad humana en la composición química del agua.

2. Área de estudio

El área de estudio está delimitada por límites político-administrativos y corresponde al municipio de Loreto; sin embargo, para poder entender el contexto hidrogeológico del área se utilizaron muestras de agua subterránea de sitios fuera del municipio.

El Municipio de Loreto se localiza en el estado de Baja California Sur, México y abarca una superficie de 4402 km², (entre 26° 33' 11" y 25° 12' 15" latitud norte; 111° 46' 22" y 110° 55' 15" longitud oeste). Al N colinda con el Municipio de Mulegé, al S y E con el Municipio de Comondú; al E con el Golfo de California (Figura 1). Cuenta con una población de 11839 habitantes (INEGI, 2005). Las actividades principales son el turismo y la pesca. El principal uso de suelo en el municipio corresponde a la agricultura, ganadería, minería, pesca y turismo, principalmente.

Según la clasificación de García (1988) presenta clima árido, cálido con régimen de lluvias en verano BW(h') HW(x')(e'). Las lluvias que se producen en verano son esporádicas de baja intensidad y corta duración, debido a que son resultado de ciclones debilitados. El régimen pluvial se caracteriza por presentar un periodo de escasa lluvia entre los meses de febrero a junio y un periodo de lluvias de agosto a septiembre con una media anual de ≈ 115 mm (Steinitz et al., 2005). La temperatura media anual es de 21° C (mínima 9° C y máxima 37.5° C) (CONANP, 2000).

El territorio correspondiente al área de Loreto se encuentra ubicado dentro de la provincia geológica denominada Faja Volcánica de la Sierra La Giganta. Su origen se remite al Mioceno temprano cuando la península se encontraba unida al macizo continental y el terreno era invadido por secuencias estrato-volcánicas y coladas de lava proveniente de volcanes en lo que hoy es la costa del estado de Sinaloa. Las rocas volcánicas que constituyen la Sierra de La Giganta fueron denominadas por Heim (1922) como Formación Comondú (Tco), compuestas por material volcánico extrusivo así como sedimentario con un espesor superior a los 1000 m. Posteriormente, Beal (1948) y Mina-Uhink (1957) limitaron esta definición a una secuencia de rocas volcánicas y vulcano-clásticas en la Sierra de la Giganta, donde el último autor estimó un espesor de 1200 m. El Instituto de Geofísica (1986) sugiere, para esta formación, un espesor de > 3000 m. Las partes topográficamente bajas y la planicie costera están representadas por grava y arena así como por un conglomerado de reciente formación que puede alcanzar en Loreto hasta los 70 m de espesor (Demant, 1975). Afloramientos presentes en el área permiten identificar cuatro unidades básicas: rocas ígneas extrusivas, rocas sedimentarias, depósitos recientes y rocas ígneas intrusivas.

Las rocas ígneas extrusivas están formadas por una secuencia estrato-volcánica y coladas de lava andesítica, toba dacitica, basalto, brecha volcánica, toba ácida y material vulcano-clástico y escasos derrames andesítico y dacítico. Estas rocas muestran fracturas y fallas normales escalonadas formando escarpes y terrazas estructurales; la variabilidad de estructura y composición química en sus afloramientos propicia intemperismo y erosión diferencial. Esta unidad representa el 66% del total de la superficie del área con una geomorfología de montañas y lomeríos con laderas escarpadas y mesetas aisladas con elevaciones de hasta 1680 msnm. Las condiciones físicas de estas rocas las hace de singular importancia en una comparativa rápida filtración y movimiento del agua subterránea que en turno le permiten tener una firma geoquímica particular (Carrillo-Rivera, 2000). Por su parte, las rocas sedimentarias, tienen intercalaciones de sedimentos marino y terrígeno en un ambiente de transgresión durante el Mioceno tardío hasta el Plioceno (Formación Salada, Tps), y se asocia a la apertura del Golfo de California; consiste en un espesor de más de 1200 m compuesto de arenisca, conglomerado y arenisca terrígena, limonita-arenisca y lutita-arenisca de origen marino, con alternancia de limolita y horizontes de coquina y conglomerado polimíctico. Cubriendo a estos materiales se encuentran andesita, basalto y brecha basáltica en un espesor conjunto que varía entre los 100 y 200 m. Las rocas sedimentarias cubren 24.5 % del área, con una geomorfología de pie de monte y escasamente de montaña y lomerío abrupto (Cardona et al., 2004).

Los depósitos recientes son de origen fluvial localizados en lechos de cauce de río (intermitente) resultado de la erosión de unidades adyacentes, se diferencian en terrazas y pie de monte (Qp), y aluviones (Qal) por su posición morfológica. En la parte N del área se encuentra arena media a fina de color gris con pocos clastos redondeados de origen volcánico mientras que al S se compone por arena media a gruesa de color claro. Se encuentra hacia el NO y extremo O del área ocupando el 9% de su superficie.

Las rocas ígneas intrusivas, corresponden a un granito y granodiorita (Kgr) de edad Mesozoica que aflora en la parte media del área con una superficie que incorpora el 0.5 % de la misma. El Kgr se manifiesta como bloques de color rosa y la granodiorita es gris oscuro; éste se presenta como lomeríos de pendiente suave y elevación de < 300 m (Bigioggero et al, 1996).

La geomorfología del área de Loreto corresponde a una angosta planicie costera de entre 1 y 4 km de ancho, la cual se extiende desde el poblado de Loreto hasta el sitio conocido como Ensenada Blanca. En este mismo referente, se encuentra San Juan B. Londó, hacia el límite N del área, el cual corresponde a una planicie aluvial con ligera pendiente al SE, que se eleva 30-40 m sobre la planicie.

Dorsey et al. (1997), estudiaron el proceso estratigráfico, sedimentológico y tectónico de la cuenca Plioceno-Cuaternaria del SE de Loreto encontrando litofacies que incluyen caliza siliclástica marina y caliza bioclásticas, contenidas por la Formación Comondú.

Por otro lado, Bigioggero et al. (1996) expusieron características estratigráficas y petrográficas del cinturón andesítico del Mioceno en el área de estudio que incluye, entre otras, una serie de domos dacítico-andesítico, algunos de ellos ricos en minerales como el piroxeno, anfíbol, biotita, plagioclasas, y raramente cuarzo.

Desde el punto de vista de estructura geológica la disposición de los materiales litológicos en el área de estudio pueden ejemplificarse como una estructura de pilares y fosas como la descrita para San Juan B. Londó (Carrillo-Rivera, 2000), donde la estructura está dominada y delimitada por un sistema de fallas N-NO de unos 100 km de longitud que cruzan varias cuencas superficiales. La roca basamento es el granito-granodiorita cuya topografía está cubierta por rocas de la Formación Comundú y material Cuaternario. El espesor conjunto de material volcánico Terciario y Cuaternario es de 1600 a > 3000 m (Instituto de Geofísica, 1986). Las fosas contienen piroclástos y flujos de lava inter-estratificados con sedimentos no consolidados (arcilla, limo, arenisca, y conglomerado). Es importante enfatizar que a lo largo de la Península los materiales volcánicos de la Formación Comondú están distribuidos más allá de cada cuenca superficial. En el área de estudio se encuentra presente una capa de material de grano fino y compacto (localmente referida como arcilla), de unos 50100 m de espesor a unos 50 m por debajo de la superficie de la planicie; esta capa está cubierta por material aluvial del Cuaternario (Figura 3). La actividad geotérmica más cercana se presenta en Las Tres Vírgenes, unos 250 km al N de San Juan B Londó; es de anotar que la CFE (1998) reporta que la actividad tectónica ha producido un sistema de fallas y fracturas con una orientación N-NE el cual se vincula como una zona de conductividad hidráulica mayor que ha permitido la migración de flujos hidrotermales a la superficie.

2.2. Referente hidrogeológico

Varias de las unidades geológicas arriba descritas forman en el área un sistema acuífero constituido por dos acuíferos (Figura 2): uno menor, libre y colgado que se encuentra sobre otro aquí denominado Acuífero Principal, descritos por Carrillo-Rivera (2000) para San Juan B. Londó. El acuífero colgado es somero en material granular sedimentario y piro-clastos indiferenciados del Cuaternario y Mioceno, su espesor es de 50 m y sobre-yace a un cuerpo de arcilla compacta de 20 a 80 m de espesor que se considera cubre parcialmente la planicie en espesores diversos (> 100 m).

El Acuífero Principal es el que comúnmente se usa para extraer agua en la región; se encuentra en rocas fracturadas vulcano-sedimentarias y material granular (Formación Comondú) con una amplia distribución geográfica y gran espesor (> 3000 m). Sus condiciones hidráulicas, dependiendo del sitio son libres, confinadas o semi-confinadas. Los pozos (con caudal entre 11 y 22 l s^-1) usualmente penetran sólo los primeros 20-80 m de este acuífero y terminan en toba, sedimentos y flujos de lava. Este acuífero está cubierto en lugares por la capa de arcilla compacta reportándose condiciones hidráulicas "semi-confinadas" que bien pueden estar ligadas a la respuesta ascendente al extraer agua termal. La falta de evidencia de subsidencia como de reporte de fracturas en el suelo sugieren que esta arcilla tiene baja compresibilidad. Así, la inferida baja conductividad hidráulica vertical y la baja porosidad de esta capa posiblemente restringen su aporte como semi-confinante y es posible que el caudal de esta fuente (pozo de extracción) en el área sea mínimo.

Debido a que el área de estudio está a unos 250 km de distancia de la zona hidrotermal más cercana, y como no hay geiseres y fuentes hidrotermales con fluido magmático en la vecindad, sugiere que el agua subterránea local no sea afectada por este tipo de fuente, ni en su química ni en su temperatura. Se carece de información sobre la temperatura o química del agua subterránea extraída con anterioridad; sin embargo, los usuarios reportan (en Carrillo-Rivera, 2000) que se ha incrementado la temperatura del agua extraída. En este trabajo se reporta que el agua subterránea tiene alta temperatura (29.0 - 46.5°C), conductividad eléctrica de 1724 - 4300 μS/cm, y concentración de B (1.8 - 9.9 mg/l), F (0.3 - 3.4 mg/l) y Na (283.3 - 436.5 mg/l), características que consideradas en el contexto hidrogeológico, sugieren que un 70 % del agua obtenida en un pozo tiene una firma diferente a aquella de un flujo local (Carrillo-Rivera, 2000), indicando la inducción de flujos de tipo intermedio hacia los pozos de extracción. Comportamiento reportado para partes de México central donde se tienen condiciones hidrogeológicas similares (Carrillo-Rivera et al., 1996; Molina-Maldonado, 1997; Carrillo-Rivera, 2000).

Reportes confiables del caudal de agua subterránea extraída en el tiempo son escasos pese al marcado desarrollo y aumento poblacional. En general, la extracción de agua subterránea en el área se ha acrecentado, y sólo hay algunos datos sobre San Juan B. Londó y el uso agrícola con una extracción de 0.37 m³/s (fines de la década de los 70) y de 0.39 m³/s (para 1985) (Instituto de Geofísica, 1986). La extracción en la década de los 70 incorporó un abatimiento de 6.5 m/mes, y para 1985 de 0.7 m/mes; esto puede sugerir cambio en las condiciones de extracción de confinadas a libres (Instituto de Geofísica, 1986) o inducción de agua termal de menor densidad (Carrillo-Rivera y Cardona, 2008).

El Acuífero Principal surte de agua potable a las ciudades de Loreto y Nopoló (CNA, 1991). El agua circula horizontalmente por estos materiales granulares a una velocidad promedio de 24 m/año (porosidad ≈ 10^-1, conductividad hidráulica horizontal ≈ 10^-4 m/s; gradiente hidráulico horizontal, ≈ 7.6x10^-4) (CONAGUA, 2006).

La profundidad al nivel estático es de 20 m (alrededor y al O del poblado de San Juan B. Londó) y 50 m (al S de este poblado) donde se presenta la profundidad mayor (CONAGUA, 2006). Este estudio considera que el Grupo Salada y el basamento granítico permiten el flujo regional de agua por fracturas y fallas, el cual se manifiesta en manantiales termales submarinos y en los (pocos) reportados en la línea de costa (Prol-Ledesma et al., 2007). La recarga se tiene principalmente en la Sierra de La Giganta y se considera que la descarga se realiza por flujos locales hacia el material granular que forma la planicie costera, y por flujos de mayor jerarquía (intermedio y regional) que descargan hacia el mar en ambas vertientes.

3. Material y Métodos

Se obtuvieron un total de 20 muestras. Las muestras 10-15 (de los pozos CNA3, CNA8, CNA4, San Nicolás, CNA6, CNA15), fueron tomadas en pozos (Figura 2) que usualmente realizan una intensa extracción (uso agrícola) y captan el Acuífero Principal en San Juan B. Londó. La muestra 19 (manantial Juncalito) fue tomada en la cercanía de la costa y corresponde con agua del Acuífero Principal en el área de Loreto; la muestra 20 (manantial hidrotermal Santispac) es del Municipio de Mulegé y fue tomada cerca de la línea de costa. Las muestras 1-9 y 16-18 fueron tomadas en la parte alta de la Sierra de La Giganta. Los datos de referencia (Tabla 3) reportados en Prol-Ledesma et al. (2007) corresponden a muestras de agua termal, manantiales submarinos y agua de mar.

Se midieron in situ en cada sitio, la temperatura, el pH, oxígeno disuelto, potencial redox (Eh) y conductividad eléctrica mediante equipos portátiles precalibrados en laboratorio y calibrados en campo (Tabla 1). Estas mediciones se realizaron en campo con el propósito de adquirir información confiable rápidamente, además de tener un control de referencia adecuado a las mediciones en laboratorio.

La determinación de alcalinidad total (HCO₃^- y CO₃^2-) fue también efectuada en campo mediante el método de Gran (Cardona, 2003). En cada sitio se tomaron 5 botellas de agua en envases de polietileno con capacidad de 60 ml; las cuales fueron lavadas previamente con una solución de ácido clorhídrico al 10 % y enjuagadas perfectamente con agua destilada (Armienta et al., 1997). El contenido de la botella para el análisis de metales pesados se filtró con una membrana de acetato de celulosa con abertura de 0.45 μm para remover las partículas suspendidas de arcilla e hidróxidos de hierro y manganeso coloidales. Posteriormente, se le añadió 10 gotas de ácido nítrico de alta pureza llevando la solución a un pH de ≈ 2 para preservar las especies disueltas. Para su transporte se mantuvieron en hielo (APHA-AWWA-WPCF, 1989) hasta el análisis correspondiente. Los macro y micro componentes fueron determinados mediante ICP-MS (Espectrometría de Emisión Atómica por Plasma acoplado a Espectrometría de Masas), en el laboratorio de la Universidad Autónoma de San Luis Potosí, México y Actlabs, Ontario, Canadá.

La calidad de los datos químicos de componentes mayores fue verificada mediante balance iónico y por la comparación entre la conductividad eléctrica real y la teórica (en ambos casos el error aceptado fue menor de 10 %). Para ello, se utilizó el algoritmo del programa HIDROGEOQUIM (Fagundo et al., 2005). Aplicando este programa se calculó además el contenido de CO₂ en equilibrio del agua. Para la clasificación de datos se utilizó el método de Kurlov (Fagundo, 1996), el cual toma en consideración para la clasificación los iones que pasan de un 20 % de meq/l. La utilidad de este método se debe a que los patrones hidrogeoquímicos (expresados en términos de relación numérica de número entero entre 1 y 8; Tabla 3), se calculan dividiendo por 10 el porcentaje de meq/l de aniones y cationes.

Los patrones hidrogeoquímicos y procesos geoquímicos que explican el origen de la composición química de las aguas se determinó mediante el sistema informático MODELAGUA (Fagundo-Sierra et al, 2001), basado en modelos de reconocimiento de patrones, balance de masa y de mezcla de aguas. El uso de este programa permitió considerar la posibilidad de que, durante la trayectoria del flujo del agua en estudio, haya ocurrido variación de concentración (i.e. evapotranspiración); en el caso de gases no se considera este factor. Cuando se definió que la composición de una muestra fue resultado de mezcla de dos grupos, se calculó el delta iónico, lo cual permite conocer la proporción de mezcla, a partir de un ión conservativo. Como tal, se tomó el ion Cl^-, ya que no se modifica en los procesos de interacción agua-roca.

4. Resultados y discusión

La composición química del agua subterránea en el área de estudio está determinada en parte por procesos de interacción agua-roca. Sin embargo, en la planicie de San Juan B. Londó existe una intensa extracción de agua subterránea para fines agrícolas (CONAGUA, 2006). Los resultados sugieren que intervienen otros procesos, como el ascenso de flujos intermedios o regionales y fuentes de salinidad que están determinando la calidad del agua en esta zona.

Los parámetros medidos en campo y los resultados de los análisis químicos para los componentes mayores son presentados en la Tabla 1 y los componentes menores son presentados en la Tabla 2.

4.1. Grupos de agua y patrones hidrogeoquímicos

Con base en el criterio que considera la determinación de los patrones hidrogeoquímicos (Fagundo, 1998) se determinaron 7 grupos de agua (Tabla 3). Como representación gráfica se utilizaron los diagramas de Stiff (1951) para mostrar los diferentes grupos de agua (Figura 4); también se incluyen el correspondiente patrón hidrogeoquímico. Esto último se determinó a partir del valor medio de la composición química del agua de cada grupo.

Con base en estos resultados (Tabla 3) se analizan las características hidrogeoquímicas de cada grupo de agua:

Grupo 1, está constituido por 2 muestras cuya composición es de tipo bicarbonatada clorurada sódica, patrón hidrogeoquímico 811-271 (Tabla 3). Su rasgo distintivo es Na⁺ > Ca²⁺, siendo despreciable el contenido de Mg²⁺, lo que indica el predominio de contacto con feldespato tipo plagioclasa sódica (albita-oligoclasa) así como disolución equivalente a la presencia de halita. El rasgo distintivo de la segunda muestra que compone este grupo es HCO₃^- > Cl^-, siendo despreciable el contenido de Ca²+. La concentración de sales totales disueltas varía (STD) entre 719 y 755 mg/l.

Grupo 2, está formado por dos muestras de agua de tipo bicarbonatada clorurada cálcica sódica magnésica, patrón hidrogeoquímico 352-271 (Tabla 3). En este grupo se aprecia alto contenido de Na+ y Ca²+ > Mg²+. La segunda muestra del grupo presenta un tipo de agua bicarbonatada clorurada cálcica magnésica sódica, ambas son típicas de la andesina y biotita así como disolución equivalente a la presencia de halita. El contenido de STD varía entre de 449-506 mg/l.

Grupo 4, se integra por 5 muestras de composición clorurada bicarbonatada sódica que tienen un patrón hidrogeoquímico 811-631 (Tabla 3). Esta agua tiene un patrón Na⁺ > Cl^-, lo que pudiera deberse a procesos de intercambio iónico de tipo directo entre Na+ y Ca²+. Los valores de STD son entre 482 y 1540 mg/l.

Grupo 5, está formado por 6 muestras cuya composición química es clorurada bicarbonatada sódica, en algunos casos cálcica (19) o cálcica magnésica, como en las muestras 14 y 15; en otros sódico cálcica como en las muestras 17 y 18, su patrón hidrogeoquímico es 532-721 (Tabla 3). A diferencia del patrón anterior, el Cl^- > Na⁺, posiblemente debido a procesos de intercambio iónico inverso entre Na⁺ y Ca²+. El contenido de STD de estas muestras es entre 371 y 2239 mg/l.

Grupo 6, lo forma una muestra de agua termal, de composición clorurada sódica con un patrón hidrogeoquímico 721-811 (Tabla 3); su contenido de STD es 17209 mg/l.

Grupo 7, está constituido por agua de manantiales termales submarinos y el agua de mar (datos tomados de Prol-Ledesma et al., 2007). Este grupo de agua clorurada sódica, con patrón hidrogeoquímicos 811-811 (Tabla 3). El contenido de STD se encuentra entre 24137 y 35402 mg/l.

En las muestras de los grupos 4 y 5, se aprecia mejor la relación entre el patrón hidrogeoquímico, el tipo de litología y los procesos geoquímicos modificadores de la composición debida a la interacción agua-roca, tales como procesos de intercambio iónico. Para apreciar mejor la relación entre la composición química del agua y la litología por donde circula se utilizó el modelo computacional MODELAGUA y se le extrajo a los datos, el contenido de NaCl, Ca²⁺ y Na⁺ que pudiera estar siendo aportado por procesos de mezcla debido al intercambio iónico. De esta manera se obtuvo un nuevo patrón hidrogeoquímico que nos permitió analizar con mayor claridad los procesos que determinan la composición química del agua subterránea en el área de estudio. Estos son:

El patrón del Grupo 1 es típico de agua en contacto con plagioclasa sódica (albita, oligoclasa) y biotita; el patrón del Grupo 2 el de un agua que ha disuelto plagioclasa sódico cálcica (andesina) y biotita; el patrón del Grupo 3 muestra mayor contenido de Mg²+ que de Ca²+, posiblemente porque el agua circula en un terreno con alto contenido de minerales constitutivos de rocas máficas como basalto (olivino), además de plagioclasa tipo andesina; el patrón del Grupo 4 (excluyendo aporte de NaCl de halita) es propio de agua que circula por minerales con alto contenido de Na⁺ y similar contenido de Ca²⁺ y Mg²⁺, posiblemente diversos minerales como plagioclasa sódico-cálcica (oligoclasa, andesina) y biotita, aunque parte del Na⁺ del patrón original podría deberse al aporte por intercambio iónico directo; el patrón del Grupo 5 (excluyendo el aporte de NaCl de halita), es parecido al del Grupo 2 en cuanto a la relación Ca²⁺/Mg²⁺, pero con déficit de Na+ con respecto al Cl^- (patrón original), lo cual pueda deberse a un proceso de intercambio iónico inverso. El patrón del agua del manantial termal cercano a la costa (Grupo 6) muestra una relación Ca²⁺/Mg²⁺ similar a la del Grupo 2 y 5, pero con apreciable déficit de Na+. Por último, el patrón del Grupo 7, que incluyen los manantiales termales submarinos, si bien muestran más parecido al agua marina, se distingue en el agua de éstos una evolución hidrogeoquímica distinta evidente en elementos traza como Li, As_total, B^-, Br^-, Sr y I^-.

Mediante el uso de modelos de balance de masa y mezcla de aguas, se procederá más adelante a interpretar, con mayor precisión aún los procesos geoquímicos que originan la composición química del agua en esta área de estudio.

4.2. Distribución espacial de los patrones hidrogeoquímicos

La variación en unidades logarítmicas de la concentración de algunos componentes mayores: Cl^-, Na⁺, K⁺ y SO₄^2-(mg/l); y de componentes menores: Br, As, B, Li y Sr (μm/l) de las muestras obtenidas en manantiales y pozos del área de estudio, así como de muestras de manantiales y agua de mar compilados de la literatura (Prol-Ledesma et al., 2007) puede apreciarse mejor en la Figura 4. Los datos están ordenados aproximadamente, de menor a mayor, con respecto al contenido de Cl^-. Esto con la finalidad de representar mejor el aumento en el contenido de algunos componentes. El resultante responde en relación con el incremento del Cl^-. De esta manera es posible mostrar los resultados que indican que los pozos sometidos a extracción agrícola y algunos manantiales incrementan su contenido en componentes que, al parecer proceden de flujos de carácter intermedio o regional (agua termal) como As_total y Li, más que de origen marino. Esto parece estar también respaldado por los valores de temperatura del agua muestreada, la cual es definitivamente cuatro grados mayor a la temperatura media ambiental del área de estudio (muestras 6, 7, 8, 10, 11, 12, 13, 14, 15, 18). En el caso de los pozos, el que presenten alto contenido de elementos traza se debe específicamente al ascenso de agua de flujo de tipo intermedio o regional por efecto de la extracción; este proceso ha sido claramente identificado en rocas similares en el Centro de México (Carrillo-Rivera et al, 1996).

Si se considera la relación de concentración entre el componente estudiado y la concentración del mismo componente correspondiente a la muestra de menor contenido de Cl^- (San José de Comondú), el aumento de concentración (con relación al Cl^-) que experimentan los iones en el agua de los pozos es notable. Esta concentración se obtiene de la división del contenido de Cl^- del agua subterránea entre el agua de lluvia.

Esta magnitud es mucho mayor que la que experimentan los manantiales, con excepción de San Nicolás y Juncalito que también es alta (Tabla 4). Esto puede deberse a la cercanía con el agua de mar. En el agua de los manantiales termales esta relación de aumento de concentración es muy elevada debido al largo recorrido subterráneo así como por el aumento de temperatura que el agua adquiere durante su recorrido y la posible presencia de agua congénita marina aún presente en su recorrido hacia la descarga.

4.3. Factor de concentración (FC) del agua subterránea respecto a la precipitación

El factor de concentración cuantifica los procesos de concentración y puede determinarse a partir de la selección de un ión conservativo. Si este es igual a 1 significa que no se producen fenómenos de concentración. Este ión es un compuesto que no participa en ninguno de los procesos geoquímicos presentes en la región y que por tanto la diferencia de masa entre el valor que presenta el mismo en la muestra y el que presenta en la referencia se debe solamente a procesos de concentración (Schoeller, 1960). Teniendo en cuenta que el Cl^- es un ión conservativo, su incremento en el agua, comparado con el de la precipitación en el área, puede tomarse como un indicador aproximado del proceso de evapotranspiración (Schoeller, 1960). El uso de este método hidrogeoquímico simple de estimación de la recarga al agua subterránea se ha aplicado con efectividad en algunas regiones áridas, donde la humedad del suelo es difícil de medir (Edmunds et al., 1988).

En el presente trabajo se tomó la concentración de agua de lluvia más cercana para la región, que correspondió a Arizona, E.U.A. La concentración de Cl^- del agua subterránea en el área de Loreto con respecto al contenido de este ión en el agua de lluvia varía entre: 9.5 - 12.8 (agua del Grupo 1); 5.9 - 6.3 (Grupo 2) y 3.3 - 10.6 (Grupo 3). Las muestras de agua de los grupos 4 y 5, afectadas por actividad agrícola, y por el recorrido largo en el subsuelo del agua muestreada proponen que el Cl^- tiende a concentrarse, por lo que la magnitud es mayor, con un rango que oscila entre 13.5 y 96.8. En general, estos valores de FC en el agua son coherentes con el carácter árido de Loreto donde la lluvia media anual es de ≈ 115 mm (Steinitz et al., 2005). La tabla 4 muestra el FC por grupo de agua. En regiones húmedas, como en la Sierra del Rosario en Cuba, se han determinado valores de FC del orden de 3.0 (Fagundo et al., 2000), mientras en otras regiones áridas de Cuba (Provincia de Las Tunas), donde el promedio de lluvia anual oscila entre 100 y 180 mm, se obtuvieron valores de FC entre 27.2 y 45.3 (Zamora, 2008).

Esto es similar a lo encontrado en la región árida de Loreto en nuestro trabajo (Tabla 3). Las muestras fueron agrupadas por grupos de agua.

4.4. Origen de la composición química del agua

La selección de minerales constitutivos de las rocas por donde circuló el agua consideró información del referente geológico y aspectos discutidos en el análisis de patrones hidrogeoquímicos originales: plagioclasa (aporte de Na⁺ y Ca²⁺), biotita (aporte de Mg²⁺ y algo de K⁺), olivino (Mg²⁺), microclino (K+), pirita (la oxidación de la pirita, mineral muy distribuido en la naturaleza podría aportar SO₄^2-), cuarzo (SiO₂), así como CO₂ biogénico y el equilibrio CO₂-HCO₃^-. Para las aguas salinas, que experimentan déficit iónico de Ca²⁺, Mg²⁺, Na⁺ y K⁺ en el balance de masa, intervienen equilibrios relacionados con la calcita, la calcita-dolomita, e intercambio iónico entre Ca²⁺ y Na⁺; así como entre Ca²⁺ y K+. Los procesos geoquímicos utilizados en este trabajo para explicar el origen de la composición química del agua se muestran en la Tabla 5 (Appelo y Postma, 1993).

El modelo de balance de masas aplicado para explicar el origen de la composición química del agua en los grupos con bajo contenido salino (1, 2 y 3) utilizó como muestra de referencia el agua de lluvia (Arizona, E.U.A.), mostrando que los principales procesos geoquímicos corresponden al aporte de Na⁺ y Ca²⁺ por la albita-caolinita y la anortita respectivamente, entre los iones principales (Tabla 6). Para los grupos de agua con elevado contenido salino (4 y 5) se tomó como referencia la muestra del manantial Las Parras (para las muestras 9-18 y 20) y de Santa Rosa (para la muestra 19), como referencia de agua dulce. Como miembro extremo de la mezcla el manantial BC4 (Prol-Ledesma et al., 2007) por su alto contenido de NaCl. Los resultados muestran las reacciones de intercambio iónico así como el aporte equivalente a la disolución de halita (Tabla 7).

Para seleccionar la muestra de referencia se consideró que de acuerdo con el tipo de patrón hidrogeoquímico del agua, los grupos 1, 2 y 3 de menor salinidad adquieren su composición mediante procesos de interacción agua-roca entre la lluvia y los minerales incluidos en la litología del terreno (Tabla 6). Los segundos, que corresponden a agua obtenida en pozo, más salina, que adquiere su composición química mediante un proceso combinado de interacción agua-roca entre la lluvia y medio por el que circula el agua así como de un proceso adicional que es el aporte de un flujo de mayor mineralización, posiblemente de tipo intermedio o regional, que es inducido a la zona de extracción de agua del pozo (Tabla 8). En ambos casos se aprecia el aporte de la disolución equivalente a la presencia de halita.

Como se mencionó anteriormente, una forma de interpretar el modo en que el agua adquiere su composición en el caso de las muestras con alto contenido salino es mediante procesos de balance de masa y mezcla con flujos de carácter intermedio o regional. En estos casos se observa la precipitación significativa de calcita, compensándose el contenido de Ca²⁺ adquirido mediante procesos de intercambio iónico entre Na⁺ y Ca²⁺. También intervienen procesos de intercambio iónico entre Na⁺ y Ca²⁺ así como entre K⁺ y Ca²⁺, oxidación de pirita o reducción de SO₄^- y disolución y precipitación de SiO₂. La tabla 8 muestra que el porcentaje de agua marina varía en la mezcla, en los manantiales entre 0.1 y 2.7 %, pero en el manantial termal Santispac alcanza un 64%. La mezcla en este último manantial puede ser debida a que al estar en una oquedad sobre la línea de costa permite la entrada de agua de mar a la misma.

Para determinar el tipo de plagioclasa por la que circula el agua de las muestras pertenecientes a los Grupos 1, 2 y 3, se utilizó el aporte de Na⁺ (albita) y Ca²⁺ (anortita) del agua.

Multiplicando la proporción de moles con que cada mineral (formado, precipitado o transformado) participa en los procesos geoquímicos de interacción agua-roca por el peso molecular de dicho mineral (mg/l), se puede determinar la masa transferida en proceso de intemperismo que dio origen a la composición química del agua (Appelo y Postma, 1993). En estos cálculos se utilizaron las siguientes expresiones: Halita = 58.5 ⋅Cl^- (halita); Plagioclasa = 274 *⋅Na+ (albita) + 286 * ⋅Ca²⁺ (anortita); Microclino = 290 * K⁺ (microclino); Biotita = (172.3)* ⅓ Mg²⁺; Olivino = 120 * Mg²⁺; Pirita = ½ SO₄^2-; Cuarzo = 96 * H₄SiO₄. Los resultados muestran la disolución equivalente a la presencia de halita, disolución de plagioclasa y precipitación de cuarzo como procesos principales. Se aprecia que en el intemperismo originado mediante proceso de interacción agua-roca, se produce una apreciable disolución equivalente a la presencia de halita y transformación de plagioclasa en caolinita en todos los casos. La plagioclasa tipo albita es la fase mineralógica más abundante en la parte alta de la Sierra La Giganta y se considera su disolución como equivale en los manantiales muestreados para este trabajo; puede constituir del 50 % al 90 °% de los fenocristales presentes en las lavas y cristales de la matriz rocosa (Bigioggero et al., 1996). Varios manantiales (muestras 5, 6, 7 y 8), ubicados principalmente en la región más occidental, originan también notables procesos de disolución de olivino presente en el basalto que compone el llamado Cinturón Andesítico Miocenico (CAM) ubicado en el área de estudio (Bigioggero et al., 1996). En todas las muestras se aprecia precipitación abundante de cuarzo, compuesto presente en las rocas de la Formación Comundú (Demant, 1975).

5. Conclusiones

En el área de Loreto, Baja California Sur (México), pueden distinguirse siete grupos de agua subterránea, correspondientes a manantiales de flujo local, pozos, manantiales termales y manantiales submarinos. A partir de su composición se identifican igual número de patrones hidrogeoquímicos, expresados mediante relaciones estequiométricas por medio de la notación Na⁺ + K⁺ : Ca²⁺ : Mg²+ - Cl^- : IICO. : SO².

El origen de la composición química del agua en el área de estudio está controlado, principalmente, por procesos de interacción agua-roca los cuales incluyen disolución de halita, transformación de plagioclasa en caolinita, oxidación de pirita, intemperismo de biotita así como precipitación de cuarzo.

El origen de la composición química que adquiere el agua de los grupos 4 y 5 fue posible confirmarlo y explicarlo mediante procesos de balance de masa y mezcla de aguas. Los principales procesos geoquímicos que controlan dicha composición son: la precipitación de calcita y dolomita (que se justifica por el ascenso de temperatura al mezclarse el agua fría con la caliente), el intercambio iónico entre el Na⁺ y Ca²⁺ así como entre K⁺ y Ca²⁺ (al establecerse el nuevo balance por precipitación de calcita), oxidación de pirita o reducción de sulfato y disolución o precipitación de cuarzo.

El porcentaje de agua marina en la mezcla varía, en los manantiales de flujo local entre 0.1 y 2.7 %, y en el manantial termal alcanza 64 %. Los resultados confirman que la salinidad del agua obtenida en los pozos no está relacionada con la entrada de agua de mar, sino con la inducción de flujos regionales.

El entender el funcionamiento de agua subterránea en el área de estudio tiene varias implicaciones desde la perspectiva del ordenamiento ecológico territorial. Al determinar que sólo se tienen zonas de descarga de flujos locales en el área de estudio es indicativo de que una parte importante del agua de lluvia infiltrada continúa su recorrido subterráneo hacia fuera del territorio de la Península de Baja California como flujo intermedio y regional, es decir, hacia el Golfo de California. Los manantiales submarinos detectados sugieren que esta posibilidad es plausible. El agua subterránea extraída por los pozos es una mezcla de flujos intermedio y regional, siendo posible que la mezcla extraída incremente la componente termal a medida que la extracción aumenta en el tiempo.

Los diferentes tipos de flujos están siendo utilizados por la agricultura; sin embargo, si esta extracción no se modera y al contrario, aumenta como respuesta al incremento poblacional esto implicará que la calidad química del agua obtenida por pozos continuará en su deterioro al ir evolucionando hacia la calidad de los flujos regionales (de mayor salinidad y temperatura). Es decir, hay una disponibilidad baja de agua adecuada para fines de uso agrícola y para asentamientos humanos. Esto da la pauta de la necesidad de entender el funcionamiento de los sistemas de flujo para incluirlos en el manejo del agua subterránea que pueda ser propuesto en los planes de ordenamiento ecológico territorial en regiones áridas como la aquí estudiada.

Agradecimientos

Referencias

Appelo, C.A.J., Postma, D., 1993, Geochemistry, groundwater and pollution: Leiden, The Netherlands, A. A. Balkema Publishers, 536 p.         [ Links ]

American Public Health Association-American Water Works Association-Water Environment Federation (APHA-AWWA-WPCF), 1989, Standard methods for the examination of water and waste water, 17^th Edition: Washington, D.C., USA, American Public Health Association, 1268 p.         [ Links ]

Armienta, M.A., Rodríguez, R., Aguayo, A., Ceniceros, N., Villaseñor, G., Cruz, O., 1997, Arsenic contamination of groundwater at Zimapán, México: Hydrogeology Journal, 5, 39-46.         [ Links ]

Beal, C.H., 1948, Reconnaissance geology and soil possibilities of Baja California Mexico: Geological Society of America Memoir, 31, 138.         [ Links ]

Bigioggero, B., Corona-Chávez, P.C., Tunesi, A., 1996, Características estratigráficas, petrológicas y extinción del cinturón andesítico miocénico de la región de Loreto-Las Parras, Baja California Sur, México: Revista Mexicana de Ciencias Geológicas, 13, 152-173.         [ Links ]

Canet, C., Prol-Ledesma, R.M., Torres-Alvarado, I., Gilg, H.A., Villanueva, R.E., Lozano-Santa Cruz, R., 2005, Silica-carbonates stromatolites related to coastal hydrothermal venting in Bahía Concepción, Baja California Sur, México: Sedimentary Geology, 174, 97-113.         [ Links ]

Cardona, A., 2003, Evaluación de la Contaminación difusa del agua subterránea (resumen), en IV Congreso Nacional de Aguas Subterráneas: San Luís Potosí. S.L.P, México, 26.         [ Links ]

Cardona, A., Carrillo-Rivera, J.J., Huizar-Álvarez, R., Granel-Castro, E., 2004, Salinization in coastal aquifers of arid zones: an example from Santo Domingo, Baja California Sur, Mexico: Environmental Geology, 45, 350-366.         [ Links ]

Carrillo-Rivera, J.J., 2000, Application of the groundwater-balance equation to indicate interbasin and vertical flow in two semi-arid drainage basins, Mexico: Hydrogeology Journal, 8, 503-520.         [ Links ]

Carrillo-Rivera, J.J., Cardona, A., Moss, D., 1996, Importance of the vertical component of groundwater flow: a hydrogeochemical approach in the valley of San Luis Potosi, Mexico: Journal of Hydrology, 185, 23-44.         [ Links ]

Carrillo-Rivera, J.J. y Cardona, A., 2008, Groundwater flow system response in thick aquifer units: theory and practice in Mexico (Trabajos selectos), en XXXIII-IAH Congreso Internacional, Asociación Internacional de Hidrogeólogos: Zacatecas, Zacatecas, México Editorial Balkema, Taylor y Francis Group, 25-46.         [ Links ]

Comisión Nacional del Agua (CNA), 1991, Sinopsis geohidrológica del estado de Baja California Sur: Subdirección General de Administración del Agua, Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación, 81 p.         [ Links ]

Comisión Nacional del Agua (CONAGUA), 2006, Actualización Geohidrológica de los acuíferos de San Juan B. Londó y Loreto, B. C. S, Informe Final: Lesser y Asociados S. A de C. V.         [ Links ]

Comisión Nacional de Areas Naturales Protegidas (CONANP), 2000, Programa de Manejo del Parque Nacional Bahía de Loreto: SEMARNAT, D.F., 182 p.         [ Links ]

Comisión Federal de Electricidad (CFE), 1998. Síntesis Geológica de las Tres Vírgenes, BCS, México, en Marsily. G., 1986., Quantitative hydrogeology: groundwater hydrology for engineers: Academic Press, New York, 440 p.         [ Links ]

Demant, A., 1975, Caracteres químicos principales del vulcanismo terciario y cuaternario de Baja California Sur, Relaciones con la evolución del margen continental Pacifico de México: Revista del Instituto de Geología, UNAM, 75, 19-65.         [ Links ]

Dorsey, R. J., Stone, K. A., Umhoefer, P. J., 1997, Stratigraphy, sedimentology, and tectonic development of the southeastern Pliocene Loreto Basin, Baja California Sur, México, en Johnson, M. E., Ledesma-Vázquez, J. (eds.), Pliocene carbonates and related facies flanking the Gulf of California, Baja California, Mexico: Boulder, Colorado, The Geological Society of America, Special Paper 318, 83-110.         [ Links ]

Drever, J. I. 1982. Geochemistry of natural waters: New Jersey, Prentice-Hall, Inc., 388.         [ Links ]

Edmunds, W. M., Darling, W. G., Kinniburgh, D. G., 1988, Solute profile techniques for recharge estimation in semi-arid and arid terrain, en Simmers, I. (ed.), Estimation of natural groundwater recharge: Dordecht, D. Reidel Publishing Company, 139-157.         [ Links ]

Fagundo J. R., 1996, Química del agua kárstica, Editorial Grupo de Investigación Recursos Hídricos y Geología Ambiental, Universidad de Granada, España: 11-119.         [ Links ]

Fagundo. J. R., 1998, Patrones hidrogeoquímicos y relaciones matemáticas en aguas naturales: Ingeniería Hidráulica, 19, 62-78.         [ Links ]

Fagundo, J. R., González, P., Suárez, M., Fagundo-Sierra, J., Melián, C., Llerena, M., Sánchez, L, 2000, Orígen de la composición química de las aguas naturales y minerales procedentes de acuíferos no carbonatados de la Sierra del Rosario, Contribución a la Educación y la Protección Ambiental: La Habana, Editorial Academia, 1, 204-210.         [ Links ]

Fagundo, J. R., González, P., Suárez, M. Fagundo-Sierra, J., Melián, C, Gimenéz. S., González, A., Romero, E. Orihuela, D. 2004. Procesos geoquímicos en un acuífero costero en contacto con un humedal. En: humedales de Iberoamérica. Red Iberoamericana de humedales. Programa CYTED, J. J. Neiff (editor). 188-203.         [ Links ]

Fagundo, J.R., P. González, M. Suárez Muñoz, J. Fagundo-Sierra, C. Melián, E. Álvarez, HIDROGEOQUIM, 2005, Contribución a la Educación y Protección Ambiental: ISCTN, 6, Q 58-67, La Habana.         [ Links ]

Fagundo, J. R., González, P., Fagundo Sierra, J., Álvarez, E., Suárez, M., Melián, C., 2006. Application of hidrogeochemical modeling to characterization and water quality control of coastal karst aquifer, en Demuth, S., Gustard, A., Planos, E., Scatena, F., Servat, E. (eds.), 2006, Climate variability and change: hydrological impacts: The Netherlands, IAHS Publication 308, 596-600.         [ Links ]

Fagundo, J. R., Carrillo-Rivera, J. J., Antigüedad, I., González, P., Peláez, R., Hernández, R., Cáceres, D., Hernández-Santana, J. R., Suárez, M., Melián, C., Rodríguez-Piña, M., 2008, Chemical and geological control of spring water in Eastern Guaniguanico, Pinar del Río, Cuba: Environmental Geology, Springer-Verlag, 247-267.         [ Links ]

Fagundo-Sierra, J., J.R. Fagundo, P. González, M. Suárez, 2001, Modelación de las aguas naturales Contribución a la Educación y la Protección Ambiental: La Habana, Cuba ISCTN. 2, 8 .         [ Links ]

García E., 1988, Modificaciones al sistema de clasificación de Köppen (para adaptarlo a las condiciones de la Republica Mexicana): Instituto de Geografía, Universidad Nacional Autónoma de México, Distrito Federal, México, 246 p.         [ Links ]

Garrels, R. M., Mackenzie, F.T., 1967, Origin of the chemicals composition of springs and lakes, equilibrium concepts in natural water systems: American Chemical Society, Advances in Chemistry Series, 67, 222-242.         [ Links ]

Heim, A., 1922, Notes of the tertiary of southern lower California: Geological Magazine, 59, 529-547.         [ Links ]

Instituto de Geofísica, 1986, Estudio Geohidrológico-Geofísico en la cuenca de San Juan Bautista Londó, Municipio Comondú, Reporte para SARH, México: Instituto de Geofísica, Universidad Nacional Autónoma de México, Coyoacán, México. 1-122.         [ Links ]

Instituto Nacional de Estadística y Geografía (INEGI), 2005, II Conteo de Población y Vivienda, consultado en línea http://www.inegi.org.mx/sistemas/TabuladosBasicos/Default.aspx?c=10398&s=est.         [ Links ]

Mina-Uhink, F, 1957, Bosquejo geológico del territorio sur de la Baja California, México: Boletín de la Asociación Mexicana de Geólogos Petroleros, 9, 139-270.         [ Links ]

Molina-Maldonado, A., 1997, Diferenciación hidrogeoquímica de los sistemas de flujo de agua subterránea en la zona sur-centro del estado de Aguascalientes, México: B Sc Thesis (Geology), Facultad de Ingeniería, UNAM.         [ Links ]

Parkhurst, D. L., Plummer, L. N., Thordtenson, D. C., 1982, BALANCE - A computer program for calculating mass transfer for geochemical reactions in ground water: U.S. Geological Survey Water-Resources Investigations Report 82-14.         [ Links ]

Prol-Ledesma, R.M., Canet, A., Torres-Vera, M.A., Forrest, M.J., Armienta, M.A., 2007, Vent fluid chemistry in Bahía Concepción coastal submarine hydrothermal system, Baja California Sur, Mexico: Journal of Vulcanology and Geothermal Research, 137, 311-328.         [ Links ]

Schoeller, H. 1960. Salinity of groundwater, evapotranspiration and recharge of aquifers en: Appelo, C.A.J., Postma, D., 1993, Geochemistry, groundwater and pollution: Netherlands, A. A. Balkema Publishers, 536 p.         [ Links ]

Steinitz, C. Faris, R., Vargas-Moreno, J. C., Huang, G., Shiau-Yun, L., Arizpe, O., Angeles, M., Ivanova, A., Gámez, A. E., Baird, K., Maddhock, T., Ajami, H., Huato, L., Haro. M.J., Flaxman, M., Ganster, P., Villegas, A., López, C., 2005, Futuros alternativos para la región de Loreto, Baja California Sur, México: Informe Técnico. Harvard University.         [ Links ]

Stiff, H. A., 1951, The interpretation of chemical water analysis by means of patterns: Journal of Petroleum Technology, 1951, 15-17.         [ Links ]

Tóth, J., 2000, Las aguas subterráneas como agente geológico: causas, procesos y manifestaciones: Boletín Geológico y Minero, 111, 9-25.         [ Links ]

Zamora, V. M., 2008, Geoquímica de las aguas subterráneas de la Provincia de Las Tuna: La Habana, Universidad de La Habana, Cuba, Facultad de Química, tesis doctoral 180 p.         [ Links ]