Tesis: Coordinación de Desarrollo Profesional E Institucional
Tesis: Coordinación de Desarrollo Profesional E Institucional
Tesis: Coordinación de Desarrollo Profesional E Institucional
PROFESIONAL E INSTITUCIONAL
SUBCOORDINACIÓN DE POSGRADO
TESIS
presenta
Jesús Arturo Rabadán Ramírez
F.C5.02.44.1
Con fundamento en los artículos 21 y 27 de la Ley Federal del Derecho de Autor y como titular
de los derechos moral y patrimoniales de la obra titulada “EVALUACIÓN HIDROGEOLÓGICA
PARA DETERMINAR EL COMPORTAMIENTO DE LOS SISTEMAS DE FLUJO
SUBTERRÁNEO: OBRAJUELO-QUERÉTARO”, otorgo de manera gratuita y permanente al
Instituto Mexicano de Tecnología del Agua, autorización para que fijen la obra en cualquier
medio, incluido el electrónico, y la divulguen entre su personal, estudiantes o terceras personas,
sin que pueda percibir por tal divulgación una contraprestación.
__________________________________________________
Jiutepec, Morelos a 17 de diciembre del 2019
2
F.C5.02.43.1
DEDICATORIA
Dedico esta tesis de investigación a mis padres Rosalía Ramírez Lagunas y Crescencio
Rabadán Delgado por todo el apoyo incondicional en todos los aspectos de mí vida y
haber aportado parte de sus enseñanzas para hacerme de una persona de la que soy
actualmente.
A mis hermanas y demás familia en general por el apoyo motivacional que me brindaron
A todas las personas, tutor, comité, profesores y amigos que durante este proceso de
3
AGRADECIMIENTOS
Doy infinitamente gracias a la divinidad por haber dado la oportunidad de tener esta
experiencia de realizar una maestría y por toda la ayuda y el inmenso apoyo que me ha
del posgrado, por haberme aceptado y ser parte de esta gran Institución y darme la
sociedad.
Le agradezco al Dr. Manuel Martínez Morales, en primera por haberme aceptado como
toda la atención y darme una enseñanza tanto para la vida profesional como personal.
Agradezco a mi comité tutoral a la Dra. Sofía E. Garrido Hoyos, Dra. Edith Rosalba
Salcedo Sánchez, Dra. María del Pilar Saldaña Fabela y al M.I. Carlos Gutiérrez Ojeda,
Gracias a todos.
4
Lo que sabemos es una gota de agua; lo que ignoramos es el océano.
Isaac Newton
5
CONTENIDO
6
3. MÉTODOS Y MATERIALES ............................................................................................................. 58
DELIMITACIÓN DEL ÁREA DE ESTUDIO................................................................................................. 59
MEDICIÓN DEL NIVEL PIEZOMÉTRICO .................................................................................................. 60
MUESTREO ...................................................................................................................................... 63
INTEGRACIÓN Y PROCESAMIENTO DE LA INFORMACIÓN. ...................................................................... 67
ANÁLISIS DE DATOS .......................................................................................................................... 67
3.5.1 Balance iónico ....................................................................................................................... 67
3.5.2 Análisis de datos hidrogeoquímicos ...................................................................................... 68
3.5.3 Análisis estadístico descriptivo y multivariado ...................................................................... 69
3.5.4 Análisis de la calidad del agua .............................................................................................. 70
EVALUACIÓN HIDROGEOLÓGICA ........................................................................................................ 73
CAPITULO IV. ............................................................................................................................................ 74
4. DESCRIPCIÓN DEL ÁREA DE ESTUDIO ........................................................................................ 74
MARCO GEOGRÁFICO........................................................................................................................ 74
4.1.1 Localización y Extensión del Área......................................................................................... 74
4.1.2 Clima ..................................................................................................................................... 75
4.1.2.1 Temperatura ................................................................................................................. 76
4.1.2.2 Precipitación ................................................................................................................. 77
4.1.3 Uso de suelo y vegetación .................................................................................................... 78
4.1.4 Zonas Agrícolas..................................................................................................................... 79
4.1.5 Zonas industriales ................................................................................................................. 80
MARCO FISIOGRÁFICO ...................................................................................................................... 81
4.2.1 Geomorfología ....................................................................................................................... 83
MARCO HIDROLÓGICO ...................................................................................................................... 85
4.3.1 Hidrología superficial. ............................................................................................................ 85
4.3.1.1 Región hidrológica- administrativa: VIII. Lerma- Santiago-Pacífico. ............................ 85
4.3.1.2 Región Hidrológica: No. 12, Lerma-Santiago. .............................................................. 85
4.3.1.3 Cuenca Rio la Laja ....................................................................................................... 86
4.3.1.4 Subcuenca Río Apaseo ................................................................................................ 86
4.3.2 Hidrología subterránea. ......................................................................................................... 88
MARCO GEOLÓGICO ......................................................................................................................... 91
4.4.1 Estratigrafía ........................................................................................................................... 91
4.4.2 Geología Estructural .............................................................................................................. 99
4.4.2.1 Sistema de fallas de Querétaro .................................................................................. 100
4.4.2.2 Sistema de Fallas NE-SW .......................................................................................... 102
CAPITULO V. ........................................................................................................................................... 104
5. RESULTADOS ................................................................................................................................. 104
ANÁLISIS PIEZOMÉTRICO ................................................................................................................. 104
5.1.1 Profundidad del nivel estático y dirección del flujo.............................................................. 104
5.1.1.1 Profundidad del nivel del agua subterránea ............................................................... 105
5.1.1.2 Dirección de flujos subterráneos ................................................................................ 106
RESULTADOS Y UBICACIÓN DE LAS MUESTRAS DE AGUA.................................................................... 107
HIDROGEOQUÍMICA ......................................................................................................................... 111
5.3.1 Balance iónico (cálculo del error) ........................................................................................ 111
5.3.2 Familias de agua ................................................................................................................. 112
5.3.2.1 Familia Bicarbonatada sódica .................................................................................... 113
5.3.2.2 Familia Bicarbonatada Cálcica ................................................................................... 113
5.3.3 Índice de Saturación ............................................................................................................ 117
ANÁLISIS ESTADÍSTICO Y CALIDAD DEL AGUA .................................................................................... 121
5.4.1 Parámetros fisicoquímicos .................................................................................................. 121
5.4.2 Elementos mayores ............................................................................................................. 123
5.4.3 Elementos minoritarios ........................................................................................................ 124
7
5.4.4 Elementos traza................................................................................................................... 124
ÍNDICE DE LA CALIDAD DEL AGUA ..................................................................................................... 125
ANÁLISIS ESTADÍSTICO BIVARIADO ................................................................................................... 129
5.6.1 Coeficiente de correlación ................................................................................................... 129
5.6.2 Índices hidrogeoquímicos.................................................................................................... 132
ANÁLISIS DE LOS PARÁMETROS FISICOQUÍMICOS E HIDROGEOQUÍMICOS. ........................................... 139
5.7.1 Temperatura ........................................................................................................................ 139
5.7.2 Potencial de Hidrogeno (pH) ............................................................................................... 140
5.7.3 Potencial de óxido reducción (Eh)....................................................................................... 142
5.7.4 Conductividad Eléctrica (CE) .............................................................................................. 143
5.7.5 Oxígeno Disuelto (OD) ........................................................................................................ 145
5.7.6 Sólidos totales disueltos (STD) ........................................................................................... 146
5.7.7 Iones mayoritarios ............................................................................................................... 148
5.7.7.1 Sodio (Na+) ................................................................................................................. 148
5.7.7.2 Potasio (K+) ................................................................................................................. 149
5.7.7.3 Calcio (Ca2+) ............................................................................................................... 150
5.7.7.4 Magnesio (Mg2+) ......................................................................................................... 151
5.7.7.5 Cloruros (Cl-) ............................................................................................................... 152
5.7.7.6 Sulfatos (SO42-) ........................................................................................................... 153
5.7.7.7 Bicarbonatos (HCO3-) ................................................................................................. 154
5.7.7.8 Nitratos (NO3-)............................................................................................................. 155
5.7.8 Iones minoritarios ................................................................................................................ 156
5.7.8.1 Boro (B3+) .................................................................................................................... 156
5.7.8.2 Fluoruros (F-) .............................................................................................................. 157
5.7.8.3 Hierro (Fe) .................................................................................................................. 158
5.7.8.4 Litio (Li+) ...................................................................................................................... 159
5.7.9 Elementos traza................................................................................................................... 160
5.7.9.1 Aluminio (Al3+) ............................................................................................................. 160
5.7.9.2 Arsénico (As) .............................................................................................................. 161
5.7.9.3 Cadmio (Cd2+) ............................................................................................................. 162
5.7.9.4 Plomo (Pb2+) ............................................................................................................... 163
ANÁLISIS ESTADÍSTICO MULTIVARIADO ............................................................................................. 164
5.8.1 Análisis de conglomerados, clúster o jerárquico ................................................................. 164
ANÁLISIS DE LA VARIABILIDAD TEMPORAL DE LA COMPOSICIÓN HIDROGEOQUÍMICA DEL AGUA
SUBTERRÁNEA ........................................................................................................................................ 167
5.9.1 Elementos mayoritarios ....................................................................................................... 167
5.9.2 Elementos minoritarios y traza ............................................................................................ 168
ANÁLISIS DE LA RELACIÓN DE LA PROFUNDIDAD Y COMPOSICIÓN HIDROGEOQUÍMICA DEL AGUA
SUBTERRÁNEA ........................................................................................................................................ 169
ORIGEN DE LA COMPOSICIÓN QUÍMICA DEL AGUA SUBTERRÁNEA ................................................... 171
ANÁLISIS HIDROGEOLÓGICO ....................................................................................................... 172
MODELO HIDROGEOLÓGICO CONCEPTUAL ................................................................................... 175
DISCUSIONES............................................................................................................................. 177
CAPITULO VI. .......................................................................................................................................... 184
6. CONCLUSIONES............................................................................................................................. 184
7. RECOMENDACIONES .................................................................................................................... 186
8. BIBLIOGRAFÍA................................................................................................................................ 187
9. ANEXOS .......................................................................................................................................... 195
8
ÍNDICE DE FIGURAS
9
Figura 24. Mapa de la ubicación de los 19 pozos de agua y un manantial termal del área de
estudio (Fuente: elaboración propia con base a datos de INEGI, 2010 a y b) ---------------------- 65
Figura 25. A) Multiparamétrico (Multi-parameter Multiline P4) para la medición del pH,
conductividad eléctrica, temperatura, oxígeno disuelto, potencial redox y sólidos disueltos totales.
B) Pozo muestreado, encargado de la Comisión Estatal del Agua de Querétaro. C) Muestreo del
manantial termal Salitre. ----------------------------------------------------------------------------------------- 66
Figura 26. Representación del diagrama triangular de Piper (A) y diagrama de Stiff (B) (Fuente:
Freeze y Cherry, 1978). ------------------------------------------------------------------------------------------ 68
Figura 27. Mapa geográfico de la ubicación del área de estudio. ------------------------------------- 74
Figura 28. Mapa del clima del área de estudio (Fuente: elaboración propia de acuerdo a datos
vectoriales de CONABIO, 1998; García et al., 1999; INEGI, 2010 a y b y CONAGUA, 2015b). 75
Figura 29. Mapa de la temperatura media anual del área de estudio (Fuente: elaboración propia
con base a datos vectoriales de Vidal-Zepeda, 1990b). ------------------------------------------------- 76
Figura 30. Mapa de precipitación media anual (1990-1992) del área de estudio. ----------------- 77
Figura 31. Tipo de suelo del área de estudio. -------------------------------------------------------------- 78
Figura 32. Mapa de las zonas agrícolas. -------------------------------------------------------------------- 79
Figura 33. Mapa de las zonas industriales del área de estudio. --------------------------------------- 80
Figura 34. Mapa de las provincias fisiográficas de México (Fuente: Raisz, 1964). --------------- 82
Figura 35. Mapa geomorfológico del área de estudio. --------------------------------------------------- 84
Figura 36.Modelo digital de elevación de la zona del valle de Obrajuelo-Querétaro (Fuente propia
con base a datos de INEGI, 2010 a y b; CONAGUA, 2015b). ------------------------------------------ 84
Figura 37. Subcuenca Río Apaseo ubicando el área de estudio. ------------------------------------- 87
Figura 38. Ubicación de los acuíferos administrativos. -------------------------------------------------- 88
Figura 39. Usos del agua subterránea. ---------------------------------------------------------------------- 90
Figura 40. Mapa geológico de la zona de Obrajuelo Guanajuato que muestran las principales
unidades litológicas. (Fuente propia: con base a datos de Alanís et al. 2001 y Consejo de
Recursos Minerales, 1999). ------------------------------------------------------------------------------------- 96
Figura 41. Columna estratigráfica de la zona de Obrajuelo-Querétaro que muestran las
principales unidades litológicas (Fuente: elaboración propia con base a datos de Alanís et al.
2001 y Consejo de Recursos Minerales, 1999).------------------------------------------------------------ 97
Figura 42. Mapa geológico-estructural de la zona de estudio que muestran las principales fallas
y fracturas de la zona de estudio (Fuente propia: con base a datos de Alanís et al. 2001 y Consejo
de Recursos Minerales, 1999). -------------------------------------------------------------------------------- 103
Figura 43. Niveles de profundidad de los 30 pozos monitoreados en el área de estudio.------ 105
Figura 44. Niveles de profundidad de los pozos muestreados para el análisis hidrogeoquímico.
------------------------------------------------------------------------------------------------------------------------ 105
Figura 45. Mapa piezométrico de la zona comprendida de Obrajuelo-Querétaro. --------------- 106
Figura 46. Mapa de la ubicación espacial de los pozos. ----------------------------------------------- 110
Figura 47. Diagrama de Piper. -------------------------------------------------------------------------------- 114
Figura 48. Diagramas de Stiff de las muestras dl agua subterránea de la zona de estudio. -- 115
Figura 49. Mapa de la variación espacial de los diagramas de stiff obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 116
Figura 50. Índice de saturación de la Albita en función de la conductividad eléctrica. ---------- 118
Figura 51. Índice de saturación del feldespato k en función de la conductividad eléctrica. ---- 118
10
Figura 52. Índice de saturación de la mica k en función de la conductividad eléctrica (Fuente:
Elaboración propia). ---------------------------------------------------------------------------------------------- 119
Figura 53. Índice de saturación de la calcita en función de la conductividad eléctrica (Fuente:
Elaboración propia). ---------------------------------------------------------------------------------------------- 120
Figura 54. Índice de saturación de la fluorita en función de la conductividad eléctrica (Fuente:
Elaboración propia). ---------------------------------------------------------------------------------------------- 120
Figura 55. Diagramas de cajas y bigotes para los parámetros fisicoquímicos de las 21 muestras
de agua (Fuente: elaboración propia). ----------------------------------------------------------------------- 123
Figura 56. Diagramas de cajas y bigotes de los elementos mayoritarios de las 21 muestras de
agua subterránea (Fuente: elaboración propia). ---------------------------------------------------------- 123
Figura 57. Diagramas de cajas y bigotes para los elementos minoritarios de las 21 muestras de
agua subterránea (Fuente: elaboración propia). ---------------------------------------------------------- 124
Figura 58. Diagramas de cajas y bigotes de los elementos traza de las 21 muestras de agua
subterránea (Fuente: elaboración propia).------------------------------------------------------------------ 124
Figura 59. Mapa de la distribución de la calidad del agua subterránea de los pozos muestreados
(Fuente propia: de acuerdo a datos vectoriales INEGI, 2010 a y b;). -------------------------------- 127
Figura 60. Diagrama de la relación de temperatura/conductividad eléctrica (T/CE). ------------ 130
Figura 61. Diagrama de dispersión de la relación CE/SD. --------------------------------------------- 131
Figura 62. Relación iónica calcio/magnesio (rCa+2/rMg+2) de las muestras del agua. ----------- 132
Figura 63. Relación iónica sodio/calcio (rNa+/rCa+2) de las muestras del agua. ------------------ 133
Figura 64. Relación iónica sodio/calcio (rNa+/rCa+2+rMg+2) de las muestras del agua. --------- 134
Figura 65. Relación iónica sulfato/cloro (rSO42-/rCl+2) de las muestras del agua. ---------------- 134
Figura 66. Relación iónica cloro y bicarbonato (rCl+2/rHCO42-) de las muestras del agua. ----- 135
Figura 67. Relación iónica cloro y sodio (rCl-/rNa+) de las muestras del agua. ------------------- 136
Figura 68. Relación iónica de la conductividad eléctrica y nitratos (rCE/rNO3-) de las muestras
del agua. ------------------------------------------------------------------------------------------------------------ 136
Figura 69. Relación iónica litio y boro (rLi+/rB+3) de las muestras del agua. ----------------------- 137
Figura 70. Relación iónica litio y flúor (rLi+/rF-) de las muestras del agua. ------------------------- 138
Figura 71. Gráfica de los datos de temperatura. --------------------------------------------------------- 139
Figura 72. Mapa de la variación espacial de los datos de temperatura obtenidos de los pozos
muestreados (Fuente: de acuerdo a datos vectoriales de INEGI, 2010 a y b). -------------------- 140
Figura 73. Gráfica de los valores de pH de los pozos muestreados. -------------------------------- 141
Figura 74. Mapa de la distribución espacial de los datos de pH. ------------------------------------- 141
Figura 75. Gráfica de los valores de Eh de los pozos muestreados. -------------------------------- 142
Figura 76. Mapa de la variación espacial de los datos de Eh obtenidos de los pozos muestreados
(Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------------------------- 143
Figura 77. Gráfica de los valores de CE de los pozos muestreados.-------------------------------- 144
Figura 78. Mapa de la variación espacial de los datos de CE obtenidos de los pozos muestreados
(Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------------------------- 144
Figura 79. Gráfica de los valores de oxígeno disuelto. ------------------------------------------------- 145
Figura 80. Mapa de la variación espacial de los datos de OD obtenidos de los pozos muestreados
(Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------------------------- 146
Figura 81. Gráfica de los valores de STD. ----------------------------------------------------------------- 147
11
Figura 82. Mapa de la variación espacial de los datos de SDT obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 147
Figura 83. Mapa de la variación espacial de los datos de sodio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 148
Figura 84. Mapa de la variación espacial de los datos de potasio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 149
Figura 85. Mapa de la variación espacial de los datos de calcio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 150
Figura 86. Mapa de la variación espacial de los datos de potasio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 151
Figura 87. Mapa de la variación espacial de los datos de cloruros obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 152
Figura 88. Mapa de la variación espacial de los datos de sulfatos obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 153
Figura 89. Mapa de la variación espacial de los datos de bicarbonato obtenidos de los pozos
muestreados (Fuente: de acuerdo a datos vectoriales de INEGI, 2010 a y b). -------------------- 154
Figura 90. Mapa de la variación espacial de los datos de nitratos obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 155
Figura 91. Mapa de la variación espacial de los datos de boro obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 156
Figura 92. Mapa de la variación espacial de los datos de fluoruros obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 157
Figura 93. Mapa de la variación espacial de los datos de hierro obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 158
Figura 94. Mapa de la variación espacial de los datos de litio obtenidos de los pozos muestreados
(Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------------------------- 159
Figura 95. Mapa de la variación espacial de los datos de aluminio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 160
Figura 96. Mapa de la variación espacial de los datos de arsénico obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 161
Figura 97. Mapa de la variación espacial de los datos de cadmio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 162
Figura 98. Mapa de la variación espacial de los datos de plomo obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b). ----------- 163
Figura 99. Dendograma de las 21 muestras del agua subterránea del área de estudio. ------- 164
Figura 100. Mapa de la distribución espacial de la clasificación del agua subterránea del área de
estudio (Fuente propia: de acuerdo a datos vectoriales INEGI, 2010b). ---------------------------- 166
Figura 101. Diagramas de cajas y bigotes de los elementos mayoritarios de los estudios de 2007,
2008 y 2018 (Fuente: elaboración propia con base a datos de Serrano, 2007 y Pérez, 2008).
------------------------------------------------------------------------------------------------------------------------ 168
Figura 102. Diagramas de blox-plot de los elementos minoritarios y traza para los años 2007,
2008 y 2018 (Fuente: elaboración propia con base a datos de Serrano, 2007 y Pérez, 2008).
------------------------------------------------------------------------------------------------------------------------ 169
Figura 103. Gráfica de la relación de la profundidad y el análisis de conglomerados. ---------- 170
12
Figura 104. Mapa hidrogeológico (Fuente: elaboración con base a datos de Alaniz-Álvarez et al.,
2001 y Consejo de Recursos Minerales, 1999).----------------------------------------------------------- 172
Figura 105. Perfiles litológicos de los pozos de la zona de Obrajuelo (Fuente: elaboración con
base a datos de Alaniz-Álvarez et al., 2001; Carreón et al. 2016; Suez-Lesser, 2017; IMTA, 2018).
------------------------------------------------------------------------------------------------------------------------ 173
Figura 106. Perfiles litológicos de los pozos de la zona de Querétaro (Fuente: elaboración con
base a datos de Alaniz-Álvarez et al., 2001; Carreón et al. 2016; Suez-Lesser, 2017; IMTA, 2018).
------------------------------------------------------------------------------------------------------------------------ 173
Figura 107. Sección hidrogeológica de norte a sur, integrando datos piezométricos y dirección
del flujo (Fuente: elaboración con base a datos de Alaniz-Álvarez et al., 2001; Carreón et al. 2016;
Suez-Lesser, 2017; IMTA, 2018). ----------------------------------------------------------------------------- 174
Figura 108. Sección geológica de oeste a este, integrando datos piezométricos y dirección del
flujo (Fuente: elaboración con base a datos de Alaniz-Álvarez et al., 2001; Carreón et al. 2016;
Suez-Lesser, 2017; IMTA, 2018). ----------------------------------------------------------------------------- 174
Figura 109. Modelo hidrogeológico del acuífero del Valle de Obrajuelo-Querétaro. Instituto
Mexicano de Tecnología del Agua (IMTA). ----------------------------------------------------------------- 176
13
ÍNDICE DE TABLAS
Tabla 1. Abatimientos por año del agua subterránea de la zona de Querétaro. ------------------- 22
Tabla 2. Clasificación de aguas subterráneas de acuerdo a la temperatura. ----------------------- 41
Tabla 3. Clasificación de aguas subterráneas de acuerdo a la temperatura. ----------------------- 42
Tabla 4. Componentes más importantes encontrados en la presente investigación. ------------- 45
Tabla 5. Índices hidrogeoquímicos. --------------------------------------------------------------------------- 49
Tabla 6. Ubicación de los 30 pozos monitoreados. ------------------------------------------------------- 60
Tabla 7. Parámetros físico químicos e hidrogeoquímicos. ---------------------------------------------- 63
Tabla 8. Localización de los pozos muestreados. --------------------------------------------------------- 64
Tabla 9. Parámetros y objetivos considerados para la evaluación del índice CCME-WQI para el
uso de abastecimiento y consumo humano (NOM-127-SSA1-1994). -------------------------------- 70
Tabla 10. Criterios para la clasificación de la calidad del agua subterránea. ----------------------- 73
Tabla 11. Datos geográficos y socioeconómicos por región hidrológico-administrativa. --------- 85
Tabla 12. Característica de la región hidrológica Lerma-Santiago. ----------------------------------- 85
Tabla 13. Características del acuífero administrativo 2201. -------------------------------------------- 89
Tabla 14. Usos del agua subterránea de pozos activos del área de estudio.----------------------- 90
Tabla 15. Composición mineralogía para los diferentes tipos de rocas. ----------------------------- 98
Tabla 16. Descripción de fallas del área de estudio. ------------------------------------------------------ 99
Tabla 17. Resultados de las mediciones piezométricas tomadas en campo.---------------------- 104
Tabla 18. Resultados de los parámetros fisicoquímicos tomados in situ. -------------------------- 107
Tabla 19. Resultados de los análisis de laboratorio de los elementos mayoritarios. ------------- 108
Tabla 20. Resultados de los análisis de laboratorio de los elementos minoritarios y traza. ---- 109
Tabla 21. Resultados del porcentaje de error del balance iónico de las 21 muestras del área de
estudio. -------------------------------------------------------------------------------------------------------------- 111
Tabla 22. Familia de aguas y tipos de flujo. ---------------------------------------------------------------- 112
Tabla 23. Resumen estadístico de la composición química de las muestras del agua subterránea.
------------------------------------------------------------------------------------------------------------------------ 121
Tabla 24. Clasificación de la Calidad del Agua. ----------------------------------------------------------- 125
Tabla 25. Resultados de los índices de calidad de los pozos de agua muestreados. ----------- 126
Tabla 26. Evaluación de la calidad del agua subterránea con respecto a los criterios de la NOM-
127-SSA1-1994 y la OMS. ------------------------------------------------------------------------------------- 128
Tabla 27. Matriz de correlación entre las variables ------------------------------------------------------ 129
Tabla 28. Elementos de origen natural de acuerdo a las características del área de estudio. 171
Tabla 29. Elementos que tienen un origen antropogénico de acuerdo a las características del
área de estudio. --------------------------------------------------------------------------------------------------- 171
14
RESUMEN
15
CAPITULO I.
1. INTRODUCCIÓN
16
embargo, se han encontrado el deterioro de la calidad del agua provocado por factores
antropogénicos (CONAGUA, 2017).
17
Planteamiento del problema
18
Figura 1. Localización hidrológica del área de estudio.
19
Antecedentes
20
Figura 3. A) Comportamiento piezométrico de 45 pozos ubicados al poniente de la Falla
Central para el periodo de 1970 al 2003. B) Comportamiento piezométrico de 24 pozos
ubicados al oriente de la Falla Central del periodo de 1970 al 2003 (Fuente: Carreón et
al., 2005).
A B
21
Pérez (2008), menciona que las características estratigráficas del subsuelo como
fallas, fracturas y discontinuidades geológicas actúan como barreras o canales del flujo
subterráneo. Por consiguiente, analizó que los niveles piezométricos de los últimos 40
años medidos en el acuífero del Valle Querétaro, registraron un descenso promedio de
120 m. Definió el acuífero Valle de Querétaro como multicapa, profundo, con
interacciones de secuencias granulares y fracturadas.
22
1.2.2 Calidad del agua subterránea
23
Nitratos (NO3-) Flúor (F-)
120 30
100
A 25
B
80 20
mg/L
mg/L
60 15
40 10
NOM 127
20 5 NOM 127
0 0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 0 3 6 9 12 15 18 21 24 27 30
No. de Pozos No. de Pozos
Pérez (2008), caracterizó la composición química del agua subterránea del Valle
de Querétaro. Los elementos más significativos para este estudio que sobrepasan los
límites permisibles establecidos por la norma (NOM-127-SSA1-1994), se muestran en la
Figura 7, A, B, C y D.
60 Nitratos (NO3-) 10 Flúor (F-)
50
A B
8
40
mg/L
mg/L
6
30
4
20
NOM 127 2 NOM 127
10
0 0
0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40
No. de Pozos No. de Pozos
0.040
Arsénico (As) 0.02
Plomo (Pb2+)
0.035 C D
0.030 NOM 127 0.02
mg/L
0.025
mg/L
NOM 127
0.020 0.01
0.015
0.010 0.01
0.005
0.000 0.00
0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40
No. de Pozos No. de Pozos
24
1.2.3 Caracterización Hidrogeoquímica
Figura 8. Diagrama de los diferentes tipos de agua determinados por cada autor.
A B C D
Figura 9. Diagramas de Piper muestran las facies químicas del agua subterránea: A)
Simbología de las facies de agua. B) Serrano, 2007. C) Pérez, 2008 y D) Suez-Lesser;
CEAQ, 2017.
25
1.2.4 Hidrogeología
López (1996) realizó un modelo geológico conceptual de un sistema hidrotermal no
volcánico, de oeste a este. Menciona que es una zona hidrotermal formada por un
sistema de cuencas escalonadas e interconectadas hidrológicamente, limitadas por
fallas regionales donde los fluidos son forzados a ascender por zonas de debilidad,
generadas principalmente por la falla San Bartolomé al oeste y la falla Querétaro al este,
que funcionan como fronteras de la cuenca y produciendo una zona de convección
forzada. Estas manifestaciones termales no están asociadas a la actividad volcánica
reciente, por lo tanto, se debe a la circulación profunda de fluidos, movidos por el
gradiente hidráulico regional (Figura 10).
Suez-Lesser; CEAQ (2017) realizó un estudio hidrogeológico del acuífero del valle
de Querétaro con el propósito de establecer el comportamiento del acuífero y determinar
posibles fuentes de abastecimiento para la zona metropolitana de Querétaro. Su sección
geológica va de oeste-este y fue elaborada con base a cortes litológicos de pozos, donde
se destacan las fallas geológicas de tipo normal dando origen a fosas y pilares tectónicos.
El acuífero está constituido por materiales granulares (arenas y gravas) con horizontes
de arcillas que se encuentra limitados principalmente por basaltos y andesitas. Presenta
un nivel estático entre 130 a 140 metros de profundidad (Figura 13).
26
Sección A
Figura 10. Modelo geológico conceptual regional de Celaya a Querétaro (Fuente: López, 1996).
W Sección B E
Figura 11. Modelo geológico de la sección del valle de Querétaro (Fuente: Pacheco, 2007).
Sección C
Figura 12. Sección geológica del acuífero de Querétaro (Fuente: Carreón et al., 2016).
Sección D
Figura 13. Sección geológica del valle de Querétaro (Suez-Lesser; CEAQ, 2017).
27
En la Figura 14, se muestran las ubicaciones de las secciones geológicas dentro
del polígono del área de estudio. Para la sección A, corresponde al modelo geológico
conceptual regional de Celaya a Querétaro de López (1996). La sección B, modelo
geológico de la sección del valle de Querétaro de Pacheco (2007). En la sección C,
pertenece al modelo geológico del acuífero de Querétaro de Carreón et al. (2016). Y
finalmente en la sección D corresponde al modelo hidrogeológico del valle de Querétaro
de Suez-Lesser; CEAQ, (2017).
Sección A:
Sección B:
Sección C:
Sección D:
Figura 14. Mapa del área de estudio mostrando la ubicación de las secciones geológicas
de los antecedentes históricos.
28
1.2.5 Isotopía
Carreón-Freyre et al. (2005) realizó un estudio isotópico del oxígeno 18 contra
deuterio (δ18O vs δD). Demostraron que el agua subterránea no proviene directamente
de la infiltración y de la recarga vertical (Figura 15). En este trabajo menciona que el
posible origen del agua subterránea se debe a la evaporación bajo las condiciones
climáticas previas a la infiltración en aguas antiguas o por la mezcla de agua no
meteóricas.
Figura 15. Gráfica de isotopos de δ18O contra δD del agua subterránea (Fuente:
Carreón-Freyre et al., 2005).
Serrano (2007) en su investigación sobre la caracterización hidrogeoquímica e
isotópica en el valle de Santiago de Querétaro, demostró que el agua subterránea no
proviene de la precipitación. Y que las posibles causas se deben a la composición
producida por la evaporación durante la recarga, efectos de altitud, paleoclimáticos y por
infiltración preferencial de la precipitación durante la época del año. El comportamiento
isotópico de las muestras de los pozos, demostró evidencias de fuentes hidrotermales
(Figura 16).
Figura 16. Gráfica de isotopos estables del agua subterránea de la ciudad de Santiago
Querétaro comparadas con la línea meteórica mundial y local (Fuente: Serrano, 2007).
29
Suez-Lesser; CEAQ (2017) en su estudio hidrogeológico del acuífero del Valle de
Querétaro determinó que el comportamiento isotópico del acuífero del valle de Querétaro
presenta evidencias de zonas de recarga de agua meteórica para la zona norte y este
del acuífero y zonas de evaporación para el área centro, norte, sur y este, como se
observa en la Figura 17.
Figura 17. Gráfica de deuterio vs oxígeno 18 para el Valle de Querétaro (Fuente: Lesser,
CEAQ, 2017).
Figura 18. Composición isotópica del agua subterránea del Valle de Querétaro (Fuente:
Mejía et al., 2018).
30
Justificación
31
Hipótesis
32
Objetivos
Querétaro.
33
CAPITULO II.
2. MARCO TEÓRICO
Hidrología
Hidrogeología
34
Hidrogeología: Actualmente se aceptan como términos adecuados como
hidrología subterránea, hidrogeología o geohidrología. Por lo tanto, se encarga del
estudio del almacenamiento, circulación, distribución, composición de las aguas en las
zonas saturadas de las formaciones geológicas, teniendo en cuenta sus propiedades
físicas y químicas, sus interacciones con el medio físico, biológico y sus reacciones a la
acción del hombre (Custodio y Llamas, 2001).
Acuífero
La palabra acuífero proviene del latín y está compuesta por dos palabras, aqua
que significa agua y fero que significan llevar, es decir, es una unidad geológica
permeable saturada que puede transmitir cantidades significativas de agua (Custodio y
Llamas, 1983; Freeze y Cherry, 1979).
35
Zona de saturación y circulación profunda: zona asociada a flujos profundos, y
se caracteriza por presentar elevadas concentraciones de minerales y gases
disueltos. Emergen al exterior a través de manantiales o pozos perforados por el
hombre (Fagundo, 2005).
36
Clasificación según al tipo de porosidad o textura como: kárstica, fisurada y porosa.
Kárstica: Son aquellas formaciones geológicas constituidas por rocas sedimentarias
o susceptibles a generar procesos kársticos, generalmente calizas, dolomías, yesos
y, en general rocas evaporíticas (Custodio y Llamas, 1983; Werner, 1996).
Fisurados o fracturados: son aquellos acuíferos que circula a través de fallas,
fisuras o grietas dentro de las formaciones de rocas ígneas (plutónicas o volcánicas),
metamórficas y sedimentarias consolidadas (muy cementadas) (Custodio y Llamas,
1983; Werner, 1996).
Porosa: Es una unidad acuífera constituida por materiales no consolidados, formada
por sedimentos granulares como arenas, areniscas o gravas, en el cual, el agua
ocupa los poros para su circulación o almacenamiento (Custodio y Llamas, 1983;
Werner, 1996).
Hidrodinámica
2.3.1 Nivel del agua subterránea
Existen diferentes herramientas para efectuar los niveles del agua subterránea,
utilizando una cinta de acero graduada, marcada con gis, sondas eléctricas, las cuales
tienen un amperímetro, voltímetro o un foco, donde se activa una vez que la sonda entra
en contacto con el agua. Los transductores de presión electrónicos se utilizan para medir
los niveles del agua con una gran precisión (25 mm de tolerancia). Estos aparatos se
conectan a una computadora para el registro de las mediciones (Martínez y Mejía, 2000).
37
Gradiente Hidráulico: Se define como la pérdida de energía experimentada por
unidad de longitud recorrida por el agua, es decir, representa la pérdida o cambio de
potencial hidráulico por unidad de longitud, medida en el sentido del flujo del agua.
Permite establecer la dirección del flujo del agua subterránea (Freeze y Cherry,
1979).
Red de Flujo: Los sistemas de flujo del agua subterránea se pueden representar en
diferentes dimensiones por un conjunto de líneas ortogonales (perpendiculares entre
sí) que se conocen como equipotenciales y líneas de flujo. Las redes de flujo son la
superposición de líneas de flujo y líneas de carga hidráulica para ciertas condiciones
de flujo, con la finalidad de representar la dirección e intensidad del flujo del agua en
dos dimensiones (Martínez y Mejía, 2000).
Dirección del Flujo Subterráneo: La determinación de la dirección del flujo permite
establecer las zonas de recarga y descarga de los acuíferos, entre otras cosas. Para
determinar la dirección del flujo del agua subterránea se requiere de la nivelación del
brocal del pozo con respecto al nivel del mar (o algún otro nivel de referencia),
distancia entre los pozos (cuando menos tres) y la medición de la carga hidráulica
en los pozos (Martínez y Mejía, 2000).
38
Los sistemas de flujo se distinguen por la profundidad del recorrido y sus
características químicas, físicas y cinéticas. El contenido de elementos químicos en el
agua subterránea dependerá del medio geológico y generalmente aumenta al ampliar la
escala espacial y temporal. En la Figura 19, se observa una sección transversal de la
distribución del flujo subterráneo, en el cual, se ubican las zonas de recarga, tránsito y
descarga (Tóth, 1999).
39
Figura 19. Esquema de los sistemas de flujo (Tóth, 1999; tomada de Peñuela–Arévalo, 2007).
Gradiente Geotérmico
La temperatura aumenta gradualmente con la profundidad a un ritmo conocido como
gradiente geotérmico. El gradiente geotérmico varía considerablemente de un lugar a
otro. En la corteza, las temperaturas aumentan, a una media de 20°C a 30°C por
kilómetro (1°C / 33 metros) (Tarbuck y Lutgens, 2005).
40
Hidrogeoquímica
Werner (1996), la define como la química de las aguas subterráneas que se ocupa
de la calidad química del agua y que proporciona conocimiento sobre el origen, evolución
y entender el funcionamiento de los sistemas de flujo subterráneos.
a) Temperatura
La temperatura del agua subterránea en zonas más superficiales, depende de las
variaciones de la temperatura ambiente o por el régimen de flujo geotérmico (1ºC cada
33 metros). Las fallas y fracturas pueden generar actividad hidrotermal, aumentando la
temperatura del agua subterránea. La temperatura afecta alguna de las propiedades
geoquímicas del agua subterránea (solubilidad de gases y minerales). Sus mediciones
se determinan directamente en campo (Deutsch, 1997).
41
Existe otra clasificación según Schoeller (1962), en el cual propone una
clasificación de aguas de acuerdo a la temperatura como se muestra en la Tabla 3.
42
Los sistemas redox están en relación directa con el pH, el oxígeno disuelto, la
presión y la temperatura. Por ejemplo, el hierro se disuelve en las aguas subterráneas
en los estados redox (valencia) Ferroso [Fe(II)] y Férrico [Fe(III)]. Mediante el Eh, es
posible calcular a partir de una concentración de hierro disuelto total [Fe(ll)+Fe(lll)] la
cantidad de hierro que está presente en las formas ferrosas y férricas, asumiendo el
equilibrio redox. Muchos de los elementos disueltos en el agua subterránea ocurren en
más de un estado redox como, por ejemplo: Mn, N, S, O y Cr (Deutsch, 1997).
43
f) Sólidos disueltos totales (SDT)
Los sólidos disueltos totales en el agua consisten en la concentración de iones
mayores que ocupan un 99%, los menores un 1% y los traza alrededor de 0.1%. La
abundancia relativa de los elementos que componen a los sólidos disueltos en el agua
subterránea depende de la movilidad química, esto debido a que conforme el agua
subterránea se infiltra se van disolviendo los materiales que encuentra en su camino, y
así, aumentando la cantidad de sólidos disueltos (Freeze y Cherry, 1979; Deutsch, 1997).
44
Tabla 4. Componentes más importantes encontrados en la presente investigación.
Elemento Fuente u origen Características
45
Cationes mayoritarios ( > 5 mg/L)
Calcio Se encuentra frecuentemente en
(Ca2+) Se origina a partir de minerales
estado de saturación y su
(calcita, plagioclasas, piroxenos,
estabilidad en solución depende del
anfíboles, apatita y yeso) y
equilibrio CO2, HCO3- y CO32-.
especialmente de rocas básicas,
Se precipita con facilidad y es muy
calizas, dolomías y anhidritas.
afectado por el cambio iónico.
Magnesio Se disuelve más lento, pero es
Se debe a la disolución de
(Mg2+) más soluble que el calcio, y tiende
minerales de rocas cristalinas
a permanecer en solución cuando
(olivino, biotita piroxenos,
este se precipita. Es afectado por el
anfíboles, silicatos magnésicos y
cambio iónico. La concentración de
ferromagnésicos) y
Mg2+ en aguas subterráneas es
especialmente en rocas básicas,
menor que la del calcio; y por lo
evaporitas. También puede
general, presenta una
tener un origen de mezcla con
concentración de 1 a 40 mg/L y en
agua de mar, contaminación
acuíferos de rocas ricas en Mg
industrial y minera.
hasta 100 mg/L.
Sodio El sodio procede de la
(Na+) desintegración de minerales de Es el catión más abundante en el
feldespatos, feldespatoides, agua de mar y es muy soluble. Es
halita y otros silicatos; y de rocas afectado por el cambio de bases.
de origen marino. Usualmente es el catión dominante
También se originan del en aguas geotermales de alta
cambio de bases de arcillas del temperatura.
mismo origen, intrusiones de La NOM-127-SSA1-1994
agua de mar y agua de lluvia. Y establece un límite permisible de
raras veces por contaminación 200 mg/L para agua potable.
urbana e industrial.
Elementos minoritarios (0.01 – 1.0 mg/L)
Boro Su origen es de aguas El ácido bórico es muy débil y se
(B) termales y puede proceder de disuelve fácilmente en agua
algunos minerales como la caliente (aguas termales).
turmalina. También puede tener Concentraciones altas afectan a la
un origen por la utilización de salud, que de acuerdo a la OMS el
determinados fertilizantes. límite máximo permisible para uso
En aguas subterráneas se y consumo humano es de 0.5 mg/L.
encuentra disuelto como H3BO4 Las concentraciones normales
(ácido bórico) van de los 0.01 - 1.0 mg/L.
Es un buen trazador para En aguas termales de ambientes
aguas usadas y aguas lixiviadas volcánicos y en salmueras de
de basureros, ya que el boro se yacimientos de hidrocarburos se
encuentra en muchos encuentran concentraciones
detergentes. elevadas de 10 - 100 mg/L.
46
Fluoruro Presenta una solubilidad baja, por
(F-) lo que las concentraciones en
aguas subterráneas es reducida,
Procede del mineral fluorita
pero si rebasa los límites
(CaF2) y generalmente de aguas
permisibles de la NOM-127-SSA1-
termales o ambientes
1994 que establece un límite de 1.5
volcánicos.
mg/L para agua potable, podría
afectar a la salud (huesos y
dientes).
Litio Es un elemento poco Es soluble y tiende a permanecer
(Li+) abundante en la corteza en disolución.
terrestre, pero se encuentra en Sus concentraciones normales son
aguas termales de origen 0.001 – 0.5 mg/L, llegando hasta 5
volcánico. mg/L en salmueras.
Potasio Proviene de minerales como la
Presenta una alta solubilidad y
(K+) ortoclasa, biotita, leucita y de
fijación en la estructura de las
rocas cristalinas y suelos
arcillas.
arcillosos. Sin embargo, puede
La concentración natural en
tener origen de actividades
aguas subterráneas es menos de
agrícolas (uso de fertilizantes), o
10 mg/L.
por actividades industriales.
Estroncio Pueden ser insoluble o soluble en
(Sr2+) agua, de acuerdo a las reacciones
químicas. Los compuestos solubles
Procede principalmente de
en agua constituyen una amenaza
minerales como la celestita,
para la salud humana.
(SrSO4), y la estroncianita,
La concentración natural está
(SrCO3).
limitada por el efecto del
intercambio iónico con las arcillas
ricas en calcio y contribuye a la
dureza del agua.
Hierro
(Fe)
Su estabilidad depende del
potencial redox; se solubiliza y
precipita con facilidad. En medios
oxidantes solo se tiene cantidades
Procede de la desintegración
mínimas de Fe disuelto.
de minerales de silicatos
Afecta la potabilidad, lo cual, la
ferriferos, sulfuros, óxidos de
NOM-127-SSA1-1994 establece un
hierro y de rocas sedimentarias.
límite de 0.3 mg/L. Produce
manchas, incrustaciones y
precipitados. El agua se vuelve
amarillenta cuando se oxida.
47
Elementos traza (< 0.1 mg/L)
Aluminio El aluminio se encuentra en
(Al3+) rocas ígneas, que contienen Es un elemento muy abundante
aluminio en forma de minerales pero poco móvil. Es difícil de poner
de alúminosilicato. Cuando en solución, y queda retenido por
estos minerales se disuelven, las arcillas.
según las condiciones químicas, El límite máximo permisible para
es posible precipitar el aluminio uso y consumo humano es de 0.02
en forma de arcillas u hidróxidos mg/L.
de aluminio o ambos.
Arsénico Proceden principalmente de
(As) minerales como arsenopirita
(FeAsS); y otros arseniuros
metálicos como la löllingita Afecta la salud y el límite máximo
(FeAs2), nicolita (NiAs), cobalto permisible para uso y consumo
brillante (CoAsS) y esmatita humano es de 0.025 mg/L.
(CoAs2).
Y puede tener un origen de los
agroquímicos.
Cadmio Se encuentra en la naturaleza
por afinidad química asociado a Es un metal tóxico para la salud
metales como el zinc, plomo y humana y el límite máximo
cobre. Puede tener origen permisible para uso y consumo
industrial y por causa de humano es de 0.005 mg/L.
fertilizantes.
Plomo
El plomo por lo general, se
(Pb2+) Es un elemento fuertemente
encuentra en los sulfuros y el
dañino para la salud humana.
mineral más común es la
El límite máximo permisible para
galeana. Sin embargo, puede
uso y consumo humano es de 0.01
tener un origen antropogénico de
mg/L.
zonas industriales.
Manganeso Afecta la potabilidad del agua, por
(Mn2+) Es similar al hierro, se da por la lo cual, produce manchas negras y
disolución de minerales de precipitados.
manganeso en ambientes El límite máximo permisible para
ácidos y reductores. uso y consumo humano es de 0.15
mg/L.
Fuente: elaborado con base a Custodio y Llamas, 1983; Hem, 1985; Werner, 1996;
Pinto, 2009; Fernández, 2009.
48
2.4.3 Relaciones iónicas
49
2.4.4 Procesos hidrogeoquímicos que determinar la composición del agua
subterránea.
Figura 20. Procesos de interacción en un sistema acuoso (Fuente: Merkel y Planer 2005).
a) Disolución- precipitación
En la hidrogeoquímica, la disolución de sustancias sólidas en el agua subterránea
es un proceso muy importante, conocido también como la mineralización del agua
subterránea. La solubilidad del mineral o cantidad de sustancias disueltas depende de
las siguientes variables: temperatura del agua, solubilidad de los minerales, tiempo del
contacto entre agua y roca, superficie del contacto, concentración de sustancias ya
disueltas, pH y potencial redox (Eh).
50
Sin embargo, la disolución o la precipitación de un mineral a partir de una solución
saturada se rige bajo la “ley de acción de masas” que se define con la siguiente Ecuación
1: (Werner, 1996; Fagundo, 2005; Appelo y Postma, 2005).
Ecuación 1
𝑎𝐴 + 𝑏𝐵 ↔ 𝑐𝐶 + 𝑑𝐷
En una disolución, cuando los productos de una reacción sean todas especies
iónicas, se dice que la disolución es congruente (Fagundo, 2005).
Ejemplo:
Mineral Reacción
Halita NaCl = Na++ Cl-
Calcita CaCO3 = Ca2+ (Ca)+ CO32-
Anhidrita CaSO4 = Ca2+ + SO42-
(Fuente: Fagundo, 2005).
51
Ejemplo:
Mineral Reacción
Dolomita CaMg (CO2)3 →Ca2+ + Mg2+ + 2(CO3)
(Fuente: Fagundo, 2005).
b) Procesos de óxido-reducción
52
c) Sorción-Desorción
Adsorción y absorción
Adsorción:
Absorción:
Intercambio
iónico:
Intercambio iónico
53
El cambio de bases afecta principalmente a los cationes más frecuentes en el agua
subterránea y forman la serie de facilidad de intercambio catiónico siguiente:
Li+ > Na+ > K+ > Mg+2 > Ca+2 > Sr+2 > Ba+2
Por su adherencia más fuerte, los cationes más a la derecha de esta serie
desplazan a aquéllos más a la izquierda. El intercambio catiónico más importante en la
hidrogeoquímica es el sodio (Na+) por el calcio (Ca+2) (Custodio y Llamas, 1983; Werner,
1996).
54
Los modelos consisten en un arreglo de ecuaciones de balance de masa para cada
elemento, ecuaciones derivadas de la ley de acción de masas y sus constantes de
equilibrio para la formación de complejos y ecuaciones que definen coeficientes de
actividad de iones individuales (Fagundo, 2005).
Ecuación 3
𝐼𝐴𝑃
𝑆𝐼 = 𝑙𝑜𝑔
𝐾𝑠𝑝
Donde:
IS es el índice de saturación.
IAP es el producto de actividad iónica.
Kps es la constante de solubilidad a una temperatura y presión específica.
55
El índice de saturación SI indica, si una solución está en equilibrio con una fase
sólida o si está sub-saturada y supersaturada en relación con una fase sólida,
respectivamente. Un valor de 1 significa una sobresaturación diez veces mayor, un valor
de -2 de insaturación cien veces en relación con una cierta fase mineral. En la práctica,
el equilibrio se puede asumir para un rango de -0.2 a 0.2. Si el valor de SI determinado
es inferior a -0.2, se entiende que la solución no está saturada en relación con el mineral
correspondiente, si SI excede de +0.2, se supone que el agua está sobresaturada con
respecto a este mineral (Merkel y Planer 2005).
Modelo hidrogeológico
Un modelo conceptual es una representación gráfica del sistema de flujo del agua
subterránea que se realiza incorporando toda la información hidrogeológica disponible
en un simple esquema conceptual (Anderson y Woessner, 1992).
56
2.5.2 Modelo hidrogeológico matemático
Los modelos matemáticos son los más utilizados en la actualidad ya que son los
que utilizan métodos numéricos para resolver la ecuación diferencial del movimiento del
agua subterránea. Son más efectivos para predecir el comportamiento de los diversos
sistemas subterráneos (Leyva, 2010).
Los modelos pueden ser utilizados para predecir las consecuencias de una acción
propuesta o en un sentido interpretativo para conocer las variables que controlan la
dinámica de un sitio específico. Por lo tanto, los modelos son utilizados, para determinar
las características de los límites del acuífero y realizar evaluaciones de la recarga y
descarga del acuífero (Leyva, 2010).
57
CAPITULO III.
3. MÉTODOS Y MATERIALES
58
Para alcanzar los objetivos antes mencionados y comprobar la hipótesis propuesta
en la presente investigación, se llevaron a cabo las siguientes actividades que describen
los métodos, técnicas e instrumentos a utilizar en el proceso de la investigación.
Se recopiló toda información útil y necesaria que definió la situación histórica hasta
la actual del área de estudio. Para el contexto hidrológico, se recopiló toda la información
disponible con la finalidad de describir brevemente la red hidrográfica de la zona, los
principales escurrimientos, entre otras características a nivel cuenca y subcuenca. Se
analizaron los datos históricos de la calidad del agua subterránea como los parámetros
físico-químicos e hidrogeoquímicos y así también los niveles piezométricos. Por otra
parte, se recopiló y describió la información espacial disponible relacionada al clima,
vegetación, suelo, topografía y geología del área de estudio. La información se obtuvo
de la Comisión Nacional del Agua (CONAGUA), Sistema Nacional de Información del
Agua (SINA), Comisión Nacional para el Conocimiento y Uso de la Biodiversidad
(CONABIO), Universidad Autónoma de México (UNAM), Instituto Nacional de
Estadísticas (INEGI), Servicio Geológico Mexicano (SGM) y del Gobierno del Estado e
inclusive del Municipio de Obrajuelo, así como en diversos bancos de información
consultados a través del Internet.
59
Medición del nivel piezométrico
En la etapa de muestreo se realizó la medición del nivel piezométrico de 30 pozos
distribuidos en la ciudad de Querétaro y de la localidad de Obrajuelo, Guanajuato (Tabla
6). Los pozos fueron seleccionados de acuerdo a la accesibilidad y permisos. Esto con
la finalidad de generar mapas piezométricos y compararlo con estudios precedentes. Se
utilizó un Sistema de Posicionamiento Global (GPS) marca GARMIN, empleando el
sistema de coordenadas geográficas WGS 84 son las siglas en inglés de World Geodetic
System 84 (que significa Sistema Geodésico Mundial 1984).
60
En el mapa de la Figura 22, se muestra la ubicación espacial de los pozos. El
monitoreo se realizó a finales del temporal de lluvias del año 2017, que consistió en la
medición del nivel estático de 30 pozos.
Figura 22. Mapa de la ubicación de los 30 pozos de agua del área de estudio (Fuente:
elaboración propia con base a datos de INEGI, 2010 a y b).
61
En la Figura 23 se muestra el trabajo de campo realizando la medición del nivel
piezométrico y posteriormente el muestreo hidrogeoquímico.
A B
C
Figura 23. A) Equipo de medición. B) Pozo el Castillo. C) Medición del nivel piezométrico
mediante una sonda.
62
Muestreo
Durante los recorridos de campo en el periodo de octubre del 2017 y abril del 2018,
se colectaron 21 muestras de 19 pozos de agua y un manantial termal, de los cuales se
midieron 6 parámetros físico-químicos in situ y 27 parámetros hidrogeoquímicos
procesados en el laboratorio de La Unidad Académica Juriquilla del Instituto de Ingeniería
de la UNAM. Utilizando los respectivos métodos analíticos para cada especie química
como se muestra en Tabla 7.
63
En la Tabla 8 se muestra la ubicación geográfica y las características de los 19
pozos y un manantial termal.
13 San Pedro Mártir 1-BIS 346894 2280236 116 Agua Potable abr-18
16 San José de los Olvera 2 352506 2273910 134 Agua Potable abr-18
64
En la Figura 24, se muestra el mapa de la ubicación de los 19 pozos muestreados
y el manantial termal Salitre, colectados en el periodo de octubre del 2017 y abril del
2018, final del temporal de lluvias y principios de estiaje, respectivamente.
Figura 24. Mapa de la ubicación de los 19 pozos de agua y un manantial termal del área
de estudio (Fuente: elaboración propia con base a datos de INEGI, 2010 a y b)
65
El muestreo se realizó de acuerdo a la norma NOM-230-SSA1-2002, salud
ambiental (Figura 25). En ella se plasman los requisitos sanitarios que se deben cumplir
en los sistemas de abastecimiento públicos y privados durante el manejo del agua, así
también, los procedimientos sanitarios para el muestreo.
A B
C
Figura 25. A) Multiparamétrico (Multi-parameter Multiline P4) para la medición del pH,
conductividad eléctrica, temperatura, oxígeno disuelto, potencial redox y sólidos
disueltos totales. B) Pozo muestreado, encargado de la Comisión Estatal del Agua de
Querétaro. C) Muestreo del manantial termal Salitre.
66
Integración y procesamiento de la información.
Análisis de datos
Una vez obtenidos los resultados de los análisis químicos de las muestras de agua,
se determinó la calidad de las muestras colectadas mediante el cálculo del valor del error
analítico del balance iónico de manera que la suma (en miliequivalentes) de cationes
debe ser igual a la suma de aniones, como se describe en la ecuación de Freeze y Cherry
(1979). Y de acuerdo a Custodio y Llamas (2001), un error de equilibrio de carga menor
o igual al 10 % se considera aceptable.
Ecuación 4
∑𝑎𝑛𝑖𝑜𝑛𝑒𝑠 − ∑𝑐𝑎𝑡𝑖𝑜𝑛𝑒𝑠
% de Error Analítico = ∗ 100
∑𝑎𝑛𝑖𝑜𝑛𝑒𝑠 + ∑𝑐𝑎𝑡𝑖𝑜𝑛𝑒𝑠
67
3.5.2 Análisis de datos hidrogeoquímicos
Familias de agua
Figura 26. Representación del diagrama triangular de Piper (A) y diagrama de Stiff (B)
(Fuente: Freeze y Cherry, 1978).
68
3.5.3 Análisis estadístico descriptivo y multivariado
Una vez obtenidos los resultados del laboratorio de las muestras de agua
recolectadas de los pozos y del manantial, se continuó con su procesamiento y análisis
de los elementos químicos, para identificar los procesos que controlan la evolución
hidrogeoquímica y comportamiento de las aguas subterráneas, procesos de interacción
agua-roca, mezclas de aguas, procesos hidrotermales, contaminación antropogénica,
transporte de contaminantes, entre otros procesos.
Análisis de correlación
consistió en el análisis los datos muéstrales para saber el grado de asociación o
correlación entre dos o más variables de una población. El grado de correlación se
expresa como un número comprendido entre -1 y +1 y se le conoce como coeficiente de
correlación.
69
3.5.4 Análisis de la calidad del agua
Para determinar la calidad del agua subterránea se utilizó el índice propuesto por
Canadian Council of Ministers of the Environment, conocido como CCME-WQI (por sus
siglas en ingles). Se consideró que el principal uso del agua extraída de los pozos en la
zona de estudio es el abastecimiento público, por lo tanto, se consideraron los criterios
para uso y consumo humano establecidos por la normatividad mexicana NOM-127-
SSA1-1994. Salud ambiental. Agua para uso y consumo humano. Considerando la
flexibilidad del índice se consideraron 17 parámetros o variables a evaluar de dicha
norma (Tabla 9).
Tabla 9. Parámetros y objetivos considerados para la evaluación del índice CCME-WQI para el
uso de abastecimiento y consumo humano (NOM-127-SSA1-1994).
No. Parámetro Unidades Limite
1 pH Log10 6.5 - 8.5
2 Sólidos Disueltos Totales mg/L 1000
3 Dureza (CaCO3) mg/L 500
4 Aluminio mg/L 0.2
5 Arsénico mg/L 0.025
6 Bario mg/L 0.7
7 Cadmio mg/L 0.005
8 Cloruros mg/L 250
9 Cobre mg/L 2
10 Fluoruros mg/L 1.5
11 Fierro mg/L 0.3
12 Manganeso mg/L 0.15
13 Sodio mg/L 200
14 Nitratos mg/L 10
15 Plomo mg/L 0.01
16 Sulfatos mg/L 400
17 Zinc mg/L 5
70
Descripción del método de cálculo
El índice de la calidad del agua se calcula mediante tres factores, que son los
componentes del índice que representa el alcance, frecuencia y amplitud de cada
variable; el cálculo de F1 y F2 es relativamente sencillo; pero en el caso de F3 se
requieren algunos pasos adicionales. A continuación, se describen los pasos para el
cálculo de los componentes y finalmente del índice.
Factor 1: Alcance
F1 (alcance o ámbito de aplicación) representa el porcentaje de variables (parámetros)
que no cumplen con los objetivos por lo menos una vez durante el periodo de tiempo
considerando (variables que no cumplen), relativo al número total de variables analizadas
(CCME, 2001) (ecuación 5):
Ecuación 5
Número de variables que no cumplen el objetivo
F1 = ( ) 𝑥 100
Número total de variables
Donde las variables indican la calidad del agua para los objetivos definidos durante el
periodo de tiempo que se calcula el índice.
Factor 2: Frecuencia
F2 (frecuencia) representa el porcentaje de cada una de las pruebas (aprovechamiento)
que no cumplen los objetivos (pruebas que no cumplen) (CCME, 2001) (ecuación 6):
Ecuación 6
Número de pruebas fallidas
F2 = ( ) 𝑥 100
Número total de pruebas
Factor 3: Amplitud
F3 (amplitud) representa la suma total de cada una de las pruebas (aprovechamiento)
que no cumplen los objetivos. F3 se calcula en tres pasos (CCME, 2001).
71
(i) El número de veces que una concentración individual es mayor (ecuación 7).
Ecuación 7
Número de pruebas que no cumplen
Desviación anómala𝑖 = ( ) -1
𝑂𝑏𝑗𝑒𝑡𝑖𝑣𝑜𝑠𝑗
Ecuación 8
𝑂𝑏𝑗𝑒𝑡𝑖𝑣𝑜
𝑗
Desviación anómala𝑖 = (Número de pruebas que ) − 1 -1
no cumplen
(ii) La suma total de los datos de las pruebas individuales que están fuera del rango, y
que son calculados al sumar las desviaciones anómalas de las pruebas individuales para
los objetivos de la evaluación, se divide entre el número total de pruebas (tanto las que
cumplen los objetivos como los que no los cumplen). Esta variable se denomina suma
normalizada o suma de la desviación anómala, o NSE, y se calcula como (CCME, 2001):
Ecuación 9
∑𝑛𝑖=1 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 𝑎𝑛ó𝑚𝑎𝑙𝑎
𝑁𝑆𝐸 =
𝑛𝑜. 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑝𝑟𝑢𝑒𝑏𝑎𝑠
(iii) F3 se calcula por medio de una función asintótica de escalas normalizadas y la suma
de la desviación anómala de objetivos (NSE), para obtener un rango entre 0 y 100
(CCME, 2001).
Ecuación 10
𝑁𝑆𝐸
𝐹3 =
0.01 𝑁𝑆𝐸 + 0.01
Ecuación 11
√𝐹12 + 𝐹22 + 𝐹32
ICA = 100 − ( )
1.732
Ecuación 11. Índice de calidad del agua 1.0, tomado del Consejo Canadiense de
Ministros de Medio Ambiente (CCME, 2011)
72
Una vez obtenido el valor del índice, se emplearon los criterios de la Tabla 10 para
determinar las condiciones y características del agua subterránea (CCME, 2001).
Evaluación Hidrogeológica
73
CAPITULO IV.
4. DESCRIPCIÓN DEL ÁREA DE ESTUDIO
Marco geográfico
4.1.1 Localización y Extensión del Área
El área de estudio se ubica entre los estados de Guanajuato y Querétaro, abarca
desde la comunidad de Obrajuelo municipio de Apaseo el Grande del estado
de Guanajuato a la ciudad de Santiago Querétaro. La zona de estudio se encuentra entre
los paralelos 20°45' y 20°28' de latitud norte y 100°37' y 100°15' de longitud oeste. Cubre
un área de 763 km2. Las principales vías de acceso son: carretera Querétaro
Celaya/México 45D para después continuar por la carretera Corregidora Obrajuelo-Ixtla
hasta llegar a la comunidad de Obrajuelo; saliendo de la ciudad de Querétaro por la
carretera estatal No. 11 e incorporándose al libramiento Norponiente y desviándose a la
derecha por la vía que conduce a Obrajuelo (Figura 27).
74
4.1.2 Clima
Figura 28. Mapa del clima del área de estudio (Fuente: elaboración propia de acuerdo a
datos vectoriales de CONABIO, 1998; García et al., 1999; INEGI, 2010 a y b y
CONAGUA, 2015b).
75
4.1.2.1 Temperatura
En la Figura 29, se muestran las condiciones de la zona centro del área de estudio,
presenta un clima semicálido, con una temperatura media anual mayor de 18ºC, la
temperatura del mes más frio es menor a 18ºC y la temperatura del mes más caliente es
mayor a 22ºC. En las zonas altas presenta climas templados-subhúmedos, con una
temperatura media anual entre 12ºC y 18ºC, la temperatura del mes más frio esta entre
-3ºC y 18ºC y la temperatura del mes más caliente es de 22 ºC. En la zona noroeste del
área de estudio se presentan climas semiáridos-templados, con una temperatura media
anual entre 12ºC y 18ºC, la temperatura del mes más frio es de -3ºC y 18ºC, la
temperatura del mes más caliente es menor a 22ºC (CONABIO, 1998; García et al.,
1999). La temperatura media anual oscila entre 12 y 20°C (CONAGUA, 2010 a y b).
Figura 29. Mapa de la temperatura media anual del área de estudio (Fuente: elaboración
propia con base a datos vectoriales de Vidal-Zepeda, 1990b).
76
4.1.2.2 Precipitación
Figura 30. Mapa de precipitación media anual (1990-1992) del área de estudio.
(Fuente: elaboración propia de acuerdo a los datos de Vidal-Zepeda, 1990).
77
4.1.3 Uso de suelo y vegetación
En el área de estudio existen dos provincias fisiográficas que influyen sobre el
comportamiento de los suelos: la Mesa Central y el Eje Neovolcánico (Figura 31).
78
4.1.4 Zonas Agrícolas
79
4.1.5 Zonas industriales
80
Marco Fisiográfico
Mesa Central
El área que ocupa la Mesa Central comprende a la subprovincia Llanuras y sierras
de Querétaro e Hidalgo. Es una zona de gran variedad morfológica, ya que presenta
llanuras y sistemas montañosos con pequeñas elevaciones aisladas, presentando
alturas que oscilan entre 2 000 y 2 960 msnm. Está compuesta por una gran variedad de
rocas de origen distinto, pudiéndose identificar grandes extensiones de relleno aluvial y
secuencias de rocas ígneas intrusivas y extrusivas, sedimentarias y metamórficas
(INEGI, 1998).
Eje Neovolcánico
El Eje Neovolcánico ha sido definido como una provincia del cenozoico continental
de origen volcánico, es una región activa de volcanes, con una orientación este-oeste
(E-W), sus dimensiones varían entre 20 y 200 km de ancho por 100 km de longitud, cuyas
altitudes oscilan entre 2000 y 3280m (INEGI, 1998). Está compuesta por rocas
volcánicas de composición variada y presenta edades que van del Mioceno medio hasta
el Holoceno. Se ha interpretado como un arco volcánico originado por la subducción de
la placa Rivera y Cocos bajo la Placa Norteamericana (Alanís-Álvarez et al., 2001). Es
una zona con gran actividad tectónica, provocando fracturas y fallas que dan origen a
pilares y fosas tectónicas que al mismo tiempo son rellenados por material clástico y/o
derrames lávicos (CONAGUA, 2015 a y b).
81
Figura 34. Mapa de las provincias fisiográficas de México (Fuente: Raisz, 1964).
82
4.2.1 Geomorfología
La zona de estudio está limitada por sierras y conos volcánicos, en los cuales se
distinguen 4 unidades, y se describen a continuación: valles, conos volcánicos, sierras
volcánicas y mesetas. En la zona se destaca el valle de Obrajuelo y el valle de Querétaro
presentando una forma muy plana con una ligera pendiente hacia el oeste. En la porción
norte, sur y oriental del valle, están limitadas por elevaciones topográficas que alcanzan
hasta 800 metros sobre el nivel de valle (Figura 35 y Figura 36).
El valle de Obrajuelo se ubica al oeste del área de estudio, limitado por elevaciones
topográficas de suave relieve, con topografía plana de 1760 msnm y está drenado por el
río Querétaro que presenta una dirección de flujo este a oeste. Los conos volcánicos se
orientan de norte y sur, constituidos por rocas extrusivas de composición basáltica, con
un drenaje de tipo radial. En las sierras volcánicas de los extremos norte y sur, alcanzan
una altura máxima de 2550 msnm y una mínima de 800 sobre el valle, mismas que
exponen laderas con fuertes pendientes y arroyos que forman parte de un drenaje de
tipo radial–dendrítico (Figura 35 y Figura 36).
83
Figura 35. Mapa geomorfológico del área de estudio.
(Fuente propia con base a datos de INEGI, 2010 a y b; CONAGUA, 2015b).
84
Marco Hidrológico
85
4.3.1.3 Cuenca Rio la Laja
La región hidrológica del río Lerma-Santiago tiene una extensión de 2496 km2 en el
estado de Querétaro. Representa el 21% de la superficie estatal. Los municipios que la
conforman son: Querétaro, Corregidora, El Marqués, Huimilpan y una pequeña porción
al sur de Amealco. En esta región se encuentra ubicada la cuenca del río Laja, con una
superficie de 2274 km2. La cuenca Laja tiene como afluentes a los ríos Pueblito y
Querétaro, que escurren finalmente hacia la zona de Lerma (CQRN, 2001).
En la sierra del sur se generan los escurrimientos del río El Pueblito, con dirección
de sur a norte y llegan y cruzan al municipio de Corregidora, donde cambia de rumbo de
sureste-noroeste y atraviesan la llanura hasta llegar a las Adjuntas, lugar donde se une
con el arroyo Madalena. En la porción norte del área en estudio se encuentran dos
arroyos de gran importancia que aportan al río la Laja y son: el arroyo San Isidro que
cruza la población de Santa Rosa Jáuregui, y el arroyo Jurica, ambos recorren una
dirección general oeste-este (CQRN, 2001; CONAGUA, 2015b).
Los principales cuerpos de agua presentes en la subcuenca Río Apaseo son: Presa
el Batán, Santa Catarina, el Carmen, Juriquilla o lago el Cajón, San Pedro Huimilpan y el
lago Salitre. El la Figura 37, se muestra el mapa de la hidrología de la zona de estudio.
86
Figura 37. Subcuenca Río Apaseo ubicando el área de estudio.
(Fuente propia: de acuerdo a datos de SIALT; INEGI, 2010 a y b)
87
4.3.2 Hidrología subterránea.
88
En la Tabla 13, se describen las características principales del acuífero administrativo
valle de Querétaro.
89
Uso del agua subterránea
Tabla 14. Usos del agua subterránea de pozos activos del área de estudio.
Censo de pozos
Uso Pozos Porcentaje
Potable 127 60
Agrícola 46 22
industrial 37 17
abrevadero 3 1
Total 213 100
(Fuente: elaboración con base a datos de Suez-Lesser 2017).
1%
17%
22%
60%
90
Marco Geológico
4.4.1 Estratigrafía
91
Edad: Su edad ha sido asignada mediante al contenido fósil encontrado en Guanajuato
del eoceno y por el fechamiento isotópico (49.3 ± 1 Ma, K-Ar roca total) (Álaniz-
Álvarez et al., 2001; Nieto-Samaniego et al., 1999).
Ignimbrita (Tig)
Litología y espesor: En esta unidad (Tig en la Figura 40) se agrupan ignimbritas de
composición ácida que presentan distintos grados de piroconsolidación, aunque en
su mayoría se observa poco soldada. Las ignimbritas están intercaladas con
depósitos de toba y vitrófidos. Pertenece a un ambiente continental. El espesor va
de los 100 a 300 m aproximadamente.
Distribución: Los afloramientos de esta unidad aparecen al norte de la localidad de
Juriquilla. Las ignimbritas yacen discordantemente sobre el basamento mesozoico,
o bien sobre la Andesita El Cedro.
Edad: Las rocas de esta unidad pertenecen al Oligoceno tardío de acuerdo a un
fechamiento isotópico de la secuencia ignimbritica realizado por el método K-Ar en
sanidinos, dando una edad de 28.6 ± 0.7 Ma. (Nieto-Samaniego et al., 1996; Álaniz-
Álvarez et al., 2001).
92
Dacita Obrajuelo (TDo)
Litología y espesor: La unidad Dacita Obrajuelo (TDo en la Figura 40) se compone de
rocas volcánicas. Son de color rosa y ricas en fenocristales de andesina, contienen
además feldespato potásico, biotita, anfíbol y bajo contenido de cuarzo. Los domos
alcanzan alturas entre 150 y 200 m desde su base hasta la cima.
Distribución: Afloran al noroeste de Querétaro a manera de domos y forman un amplio
tren estructural de dirección NE-SW, desde el volcán San Pedro hasta Santa Rosa
de Jáuregui. A esta unidad se incluye un cuerpo ígneo de textura porfídica y
composición intermedia que aflora en Juriquilla, debido a su posición estratigráfica
y mineralógica.
Edad: Se le asigna una edad entre el Oligoceno y el Mioceno medio (Álaniz-Álvarez et
al., 2001).
93
Querétaro y el Basalto Querétaro. También se puede observar derrames de basalto
perteneciente a esta unidad intercalados con la unidad Arenisca y Conglomerado
cenozoicos en la depresión.
Edad: Su edad se ubica en el Mioceno medio ya que las edades conocidas varían de
12.4 a 9.9 ± 0.4 Ma. (Pérez-Venzor et al., 1996; Álaniz-Álvarez et al., 2001).
94
Distribución: La secuencia se encuentra intercalada con la unidad vulcaniclastica, estos
sedimentos se depositaron en depresiones topográficas, generando lentes
sedimentarios.
Edad: Presentan una edad del oligoceno-mioceno (Suez-Lesser y CEAQ, 2017).
95
También se encuentran intercaladas con tobas y material piroclástico. El espesor
de esta unidad es variable según la cuenca de que se trate, en los afloramientos de
los alrededores de Querétaro, los espesores reportados son de hasta 350 m.
Distribución: Estos sedimentos afloran en las partes bajas del valle. Estos depósitos
sedimentarios constituyen el relleno de las depresiones o cuencas continentales del
área de estudio.
Edad: La edad de la unidad abarca del Oligoceno al Pleistoceno (Nieto-Samaniego et al.,
1996; Álaniz-Álvarez et al., 2001).
Figura 40. Mapa geológico de la zona de Obrajuelo Guanajuato que muestran las
principales unidades litológicas. (Fuente propia: con base a datos de Alanís et al. 2001 y
Consejo de Recursos Minerales, 1999).
96
Se realizó una recopilación de la información estratigráfica de la región de estudio
con la finalidad de conocer las capas litológicas del subsuelo. En la Figura 41, se
observan dos tipos de rocas, la primera de origen sedimentario constituida por areniscas,
conglomerados, sedimentos lacustres y secuencias vulcanoclásticas; y la segunda por
rocas de origen volcánico compuesta por basaltos, andesitas, ignimbritas y dacitas.
97
En la Tabla 15, se describe la edad, litología, simbología y mineralogía para cada
tipo de roca presente en el área de estudio.
98
4.4.2 Geología Estructural
Falla San
x 80 N05°W/E 27
Bartolomé
Sistema de fallas NE-SW
99
4.4.2.1 Sistema de fallas de Querétaro
El sistema de fallas de Querétaro está compuesto por cuatro fallas paralelas con
orientación aproximada norte-sur (N-S). Las fallas principales de este sistema son: la
falla Querétaro, la falla 5 de Febrero, la falla Tlacote y la falla San Bartolomé.
La falla Querétaro-norte, se presenta por el límite poniente del alto topográfico entre
San José Iturbide y Santa Rosa de Jáuregui. El escarpe tiene una orientación N-S y se
extiende por más de 10 km, se encuentra en contacto con ignimbritas riolíticas de 29.3
Ma y basaltos del Mioceno tardío (10 y 7 Ma), con depósitos sedimentarios
pertenecientes a la unidad de Arenisca y Conglomerados del cenozoico. La actividad de
la falla fue posterior al Mioceno medio, que es la edad de los basaltos desplazados
(Alaniz-Álvarez et al., 2001; López, 1996).
La falla Querétaro-sur, con una traza de 25.5 km de largo, muestra una orientación
N22°W. Esta falla corta a los Basaltos Querétaro y expone al volcanoclásticos Querétaro
a lo largo de su escarpe. La roca más joven se desplazada por esta falla que corresponde
a las lavas del volcán Cimatario. El desplazamiento vertical promedio es de 80 m. La falla
Querétaro-sur está cortada y limitada por la falla Huimilpan de rumbo N80°E. Esta
relación se evidencia por el truncamiento que presenta el escarpe N-S en el costado este
del volcán Cimatario. La actividad de la falla fue en el Mioceno tardío (Alaniz-Álvarez et
al., 2001).
100
4.4.2.1.2 Falla 5 de Febrero
Es una falla de tipo normal, su escarpe se observa por más de 20 km. Está
compuesta por dos segmentos separados por la falla Central, con orientación ENE-
WSW, cuya expresión topográfica y relaciones de corte indican que es anterior a la falla
Tlacote. El desplazamiento vertical fue de 80 m. Tanto en el bloque levantado como en
el bloque hundido se encuentra depositado el Basalto Querétaro. La edad de actividad
de esta falla es más joven que 7.5 Ma y no afecta a los derrames del dique El Patol de
6.2 Ma. (Alaniz-Álvarez et al., 2001).
101
4.4.2.2 Sistema de Fallas NE-SW
La falla Ixtla tiene una estructura con dirección N50°E y está compuesta por zonas
escarpadas que cortan al Basalto y Dacita Querétaro. Sobre este conjunto de fallas se
observa un alineamiento de emisiones volcánicas de composición dacítica con dirección
NE-SW (Figura 42). El bloque hundido generado por esta falla se depositan sedimentos
volcanoclásticos Querétaro, sedimentos aluviales y lacustres (Alaniz-Álvarez et al.,
2001).
La falla Sur tiene una orientación N70°E y longitudes desde 5 km, estas estructuras
cortan al volcán Cimatario y a la falla Tlacote (Figura 42). En el bloque hundido se
depositaron sedimentos lacustres y aluviales. La falla Sur desplaza a la falla Querétaro-
sur (Alaniz-Álvarez et al., 2001).
102
Figura 42. Mapa geológico-estructural de la zona de estudio que muestran las
principales fallas y fracturas de la zona de estudio (Fuente propia: con base a datos de
Alanís et al. 2001 y Consejo de Recursos Minerales, 1999).
103
CAPITULO V.
5. RESULTADOS
Análisis piezométrico
5.1.1 Profundidad del nivel estático y dirección del flujo
104
5.1.1.1 Profundidad del nivel del agua subterránea
Del monitoreo de los 30 pozos, se realizó un mapa piezométrico a nivel regional
con la finalidad de determinar la profundidad y la dirección del flujo subterráneo (Figura
45). De la línea de color café de la Figura 45 y de los datos de la Tabla 17 se realizó el
gráfico (Figura 43) donde se aprecia una profundidad mínima de 43 m, una máxima de
167 m y una profundidad media de 117 m.
No. de Pozos
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
0
20
Porfundidad (m)
40
60
80
100
120 Línea de
140 tendencia del
160 nivel freático
180
0
20
40
Nivel estático (m)
60
80
100
120
140
160
180
200
Pozos Obrajuelo-Querétaro
Figura 44. Niveles de profundidad de los pozos muestreados para el análisis hidrogeoquímico.
105
5.1.1.2 Dirección de flujos subterráneos
Con base a los datos de la Tabla 17, se realizó el mapa piezométrico del monitoreo
del año 2017 de la zona de estudio, con la finalidad de demostrar la dirección de los flujos
subterráneos. En la Figura 45, se observa un cono de abatimiento hacia las localidades
de Santa María Magdalena y el Pueblito en el Estado de Querétaro y el Castillo en el
Estado de Guanajuato.
Sección
106
Resultados y ubicación de las muestras de agua
14 Ejido Romeral 344153 2272218 29.8 7.8 256 730 4.8 283.7
Pueblito 2
15 349448 2272279 29.5 7.7 296 665 5.7 210
Garza
San José de los
16 352506 2273910 31.6 7.7 175 510 4.6 178.3
Olvera 2
San Francisco
18 353751 2267773 27.9 8.0 193 501 5.2 133.5
Corregidora
19 Poniente 1 347775 2275764 30.1 7.7 172 575 3.2 170.1
20 Zapata 349406 2274513 29.4 7.8 175.8 617 4.5 211.3
107
Tabla 19. Resultados de los análisis de laboratorio de los elementos mayoritarios.
Ca2+ K+ Mg2+ Na+ HCO3- SO42- Cl- NO3-
No POZO
mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L
1 San Miguel Ixtla 60.27 7.45 31.61 49.51 336.7 18.6 11.68 56.4
Castillo 150m 27.96 8.31 3.4 87.53 265.8 27.8 8.12 4.18
2
Castillo 200m 17.35 6.18 1.93 80.56 198.3 26.2 10.67 0.28
3 Granja Quanni 13.1 5.75 1.52 89.23 179 42.1 11.33 11.5
Manantial termal
10 0.8 12.29 0.36 200.4 267.5 40.6 89.68 0.29
Salitre
14 Ejido Romeral 44.14 15.44 19.9 78.45 235.2 90.5 25.47 32.9
15 Pueblito 2 Garza 58.01 14.72 18.18 50.68 269.7 32.6 22.42 52.9
17 Los Olvera 2 17.33 5.23 3.8 112.2 58.92 59.9 21.83 44.9
San Francisco
18 31.44 14.18 9.73 54.76 248.6 15.7 5.63 3.51
Corregidora
108
Tabla 20. Resultados de los análisis de laboratorio de los elementos minoritarios y traza.
No POZO Si4+ S Al3+ As B3+ Ba2+ Be2+ Cd2+ Cu2+ F- Fe Li+ Mn2+ Ni2+ Pb2+ Sr2+ V Zn2+
mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L
San Miguel
1 Ixtla
36.25 6.38 0.007 <l.c. 0.027 0.011 <l.c. <l.c. 0.002 0.161 <l.c. 0.026 <l.c. 0.018 0.011 0.617 0.030 0.513
Castillo 150m 44.21 9.90 0.047 <l.c. 1.024 0.019 <l.c. 0.387 0.090 0.437 0.479 0.264 0.013 0.023 0.018 0.213 0.020 0.154
2
Castillo 200m 28.32 9.37 0.018 <l.c. 1.142 0.016 <l.c. 0.553 0.089 0.693 1.102 0.211 0.032 0.019 0.016 0.127 <l.c. 0.206
Granja
3 Quanni
37.67 15.28 0.021 <l.c. 1.116 0.004 <l.c. <l.c. 0.003 0.835 0.013 0.159 <l.c. 0.018 0.013 0.090 0.034 0.064
4 La Estancia 35.83 11.60 0.005 <l.c. 0.356 0.005 <l.c. <l.c. <l.c. 0.516 <l.c. 0.094 <l.c. 0.017 0.010 0.240 0.025 0.277
Fuentes
5 Balvanera
37.81 27.39 0.007 <l.c. 11.370 0.010 <l.c. <l.c. <l.c. 1.848 0.006 0.282 <l.c. 0.017 0.011 0.723 0.042 0.236
Calera
6 Obrajuelo
35.02 28.68 0.005 <l.c. 2.579 0.006 <l.c. <l.c. <l.c. 0.748 <l.c. 0.196 <l.c. 0.017 0.009 0.404 0.026 0.547
Punta de
7 Obrajuelo
34.24 13.49 0.008 <l.c. 0.553 0.010 <l.c. <l.c. 0.020 0.424 0.002 0.095 <l.c. 0.016 0.011 0.373 0.037 0.187
8 La Norita 22.26 19.86 0.019 0.053 17.840 0.129 0.0008 <l.c. <l.c. 11.36 0.637 0.915 0.035 0.010 <l.c. 0.225 <l.c. 0.123
9 Don Manuel 33.22 22.05 0.015 <l.c. 0.471 0.020 <l.c. <l.c. 0.002 0.312 0.010 0.114 <l.c. 0.015 0.010 0.491 0.031 0.101
Manantial
10 termal Salitre 100.33 275.27 0.342 <l.c. 19.480 0.002 0.001 <l.c. <l.c. 23.09 0.012 1.399 <l.c. 0.045 0.023 0.057 <l.c. 0.295
11 Santa Fe 36.20 14.49 0.010 <l.c. 0.360 0.019 <l.c. <l.c. 0.002 0.774 0.004 0.218 <l.c. 0.017 0.012 0.341 0.022 0.621
12 Poniente 4 38.39 13.96 0.015 <l.c. 0.315 0.005 <l.c. <l.c. 0.002 0.587 0.009 0.158 <l.c. 0.018 0.010 0.287 0.023 0.506
San Pedro
13 Mártir 1-BIS 32.34 12.94 0.017 <l.c. 0.303 0.018 <l.c. <l.c. 0.003 0.430 0.006 0.098 <l.c. 0.015 0.008 0.412 0.036 0.388
Ejido
14 Romeral 34.43 34.79 0.009 <l.c. 0.776 0.010 <l.c. <l.c. <l.c. 0.648 0.003 0.122 <l.c. 0.015 0.010 0.392 0.036 0.357
Pueblito 2
15 Garza 38.31 11.74 0.007 <l.c. 0.184 0.010 <l.c. <l.c. <l.c. 0.221 <l.c. 0.078 <l.c. 0.018 0.009 0.476 0.020 0.411
San José de
16 los Olvera 2 36.98 17.52 0.005 <l.c. 0.240 0.003 <l.c. <l.c. <l.c. 0.573 <l.c. 0.091 <l.c. 0.018 0.010 0.222 0.029 0.329
17 Los Olvera 2 33.89 23.07 0.007 <l.c. 0.760 0.001 <l.c. <l.c. <l.c. 0.481 0.002 0.064 <l.c. 0.014 0.009 0.133 0.035 0.405
San
18 Francisco 33.95 5.23 0.009 <l.c. 0.093 0.022 <l.c. <l.c. 0.003 0.256 <l.c. 0.036 <l.c. 0.017 0.010 0.333 0.030 0.432
Corregidora
19 Poniente 1 37.91 10.86 0.004 <l.c. 0.285 0.010 <l.c. <l.c. 0.002 0.536 <l.c. 0.120 <l.c. 0.017 0.011 0.259 0.027 0.584
20 Zapata 37.64 19.52 0.005 <l.c. 0.281 0.008 <l.c. <l.c. 0.002 0.548 <l.c. 0.117 <l.c. 0.018 0.010 0.285 0.026 0.437
* <l.c: límite de cuantificación
109
El muestreo se realizó en el periodo de octubre del 2017 y abril del 2018. Donde 19
fueron de pozos de agua y un manantial termal. Los puntos de muestro del agua
subterránea se ubican en los municipios de Apaseo el Alto y Apaseo el grande en el
Estado de Guanajuato; y en el municipio de Santiago de Querétaro y Corregidora para
el Estado de Querétaro. En la Figura 46 se muestra el mapa de la ubicación geográfica
de los puntos de muestreo.
110
Hidrogeoquímica
5.3.1 Balance iónico (cálculo del error)
∑𝒂𝒏𝒊𝒐𝒏𝒆𝒔 − ∑𝒄𝒂𝒕𝒊𝒐𝒏𝒆𝒔
% 𝐝𝐞 𝐄𝐫𝐫𝐨𝐫 = ∑𝒂𝒏𝒊𝒐𝒏𝒆𝒔 + ∑𝒄𝒂𝒕𝒊𝒐𝒏𝒆𝒔 ∗ 𝟏𝟎𝟎
Tabla 21. Resultados del porcentaje de error del balance iónico de las 21 muestras del
área de estudio.
Suma de Suma de % Error
No. Pozo
cationes aniones (*100)
1 San Miguel Ixtla 7.95 7.15 5.3
2 Castillo 150m 5.69 5.23 4.2
3 Castillo 200m 4.69 4.10 6.6
4 Granja Quanni 4.81 4.32 5.4
5 La Estancia 5.17 4.55 6.4
6 Fuentes Balvanera 10.22 9.68 2.7
7 Calera Obrajuelo 7.40 6.64 5.4
8 Punta de Obrajuelo 5.28 5.01 2.6
9 La Norita 10.17 8.97 6.3
10 Don Manuel 6.40 5.64 6.3
11 Manantial termal Salitre 9.10 7.76 7.9
12 Santa Fe 7.21 6.50 5.2
13 Poniente 4 6.48 5.90 4.7
14 San Pedro Mártir 1-BIS 5.31 5.12 1.8
15 Ejido Romeral 7.65 6.99 4.5
16 Pueblito 2 Garza 6.97 6.58 2.9
17 San José de los Olvera 2 4.95 4.45 5.3
18 Los Olvera 2 6.19 3.55 27.1
19 San Francisco Corregidora 5.11 4.62 5.1
20 Poniente 1 5.80 5.56 2.1
21 Zapata 6.13 5.51 5.4
111
5.3.2 Familias de agua
112
De acuerdo al diagrama de Piper, se determinaron dos familiar de agua:
Bicarbonatada sódica (Na-HCO3-) y Bicarbonatada cálcica (Ca-HCO3-) como se muestra
en la Figura 47.
El orden de abundancia de los cationes mayoritarios fue de Na+ > Ca2+ > Mg2+ > K+,
mientras que para los aniones mayoritarios fue de HCO3- > SO42- > Cl- > NO3-., estos
valores permiten observar el tiempo de residencia y la evolución del agua subterránea,
que, de acuerdo a la literatura, indica una secuencia normal de un agua de reciente
infiltración, de: HCO3- - > SO42- > Cl-.
113
Figura 47. Diagrama de Piper.
(Fuente: Elaboración propia).
Diagramas de stiff
De acuerdo a los resultados del laboratorio (Tabla 19) se realizaron los diagramas
de Stiff como se muestran en la Figura 48. Para la realización de estos diagramas se
utilizaron los elementos mayoritarios en concentraciones de meq/L. La interpretación se
basa en la similitud o forma de los diagramas y se pueden observar cuatro formas
distintas. En el grupo A, se presentan dos muestras de agua de reciente infiltración: el
pozo San Miguel Ixtla y Pueblito 2 Garza. Para el grupo B y C, las muestras representan
flujos intermedios o zonas de mezcla de aguas y para el grupo D, son muestras de flujos
regionales con altos contenidos de sodio asociado a fuentes hidrotermales (Figura 48).
114
1 12
2 13
3 14
4 15
5 16
6 17
7 18
8 19
9 20
10 21
11
Figura 48. Diagramas de Stiff de las muestras dl agua subterránea de la zona de estudio.
(Fuente: Elaboración propia).
115
En la Figura 49, se muestra el mapa de la distribución espacial de los diagramas
stiff, se puede observar que para el grupo A, el pozo de San Miguel Ixtla y Pueblito 2
Garza son lugares de recarga o agua de reciente infiltración. En el grupo B y C,
representan flujos intermedios que se caracterizan por aguas tipo bicarbonatadas
sódicas. Y para el grupo D, se caracteriza por ser flujos regionales debido a que la zona
suroeste del área de estudio presenta evidencias de fuentes hidrotermales. Estos flujos
profundos ascienden a través de la falla Obrajuelo y migran hacia la dirección del flujo
subterráneo (suroeste-noreste) modificando la composición química del agua del pozo
Fuentes Balvanera y La Norita.
Figura 49. Mapa de la variación espacial de los diagramas de stiff obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
116
5.3.3 Índice de Saturación
Plagioclasas (Albita)
En la Figura 50, se observan las muestras en un rango de los 0.31 a -2.04, lo que
indica que el acuífero está subsaturada (en disolución) con respecto al mineral albita.
Para la muestra de San Miguel Ixtla está en equilibrio con respecto a la Albita, sin
embargo, para el agua del pozo el Castillo de 150 m y 200 m de profundidad se encuentra
sobresaturada y subsaturada respectivamente; debido que a mayor profundidad el
mineral la albita tiende a la disolución provocado principalmente por el aumento de la
temperatura, entre otras propiedades fisicoquímicas del ambiente.
117
2.5
0.5 CA SMI
Equilibrio
P4
-0.5 CA SF ER
DM FB
SFC
GQ SPM
LEPO ZAP2G
SJO2 P1 LO2 CO
-1.5 MTS LN
Disolución
-2.5
0 100 200 300 400 500 600 700 800 900 1000 1100 1200
CE (μS/cm)
Feldespatos K
En la Figura 51, el 80% de las muestras de agua de pozo están en equilibrio con
respecto al feldespato potásico. La muestra la Norita y el Manantial termal Salitre el agua
se encuentra subsaturada en feldespato k, y la muestra del pozo el Castillo a 150 m de
profundidad el agua se encuentra en precipitación el feldespato potásico.
2.5
Precipitación
Índice de saturación de
1.5 CA
Feldespato K
P4
SFC SF ER
DM
0.5 CA SPM ZA P2G
LE FB
Equilibrio POP1
SJO2 SMICO
GQ
-0.5
LO2
LN
-1.5 MTS
Disolución
-2.5
0 100 200 300 400 500 600 700 800 900 1000 1100 1200
CE (μS/cm)
118
Mica K
12.00
Índice de saturación de Mica K
ZA
Precipitación
9.00
SFC P1SJO2
LO2
MTSER
DM SPM P2G
P4 LN SF
6.00 FB
GQ CO PO
CA LE
CA
3.00
0.00 Equilibrio SMI
-3.00
-6.00
-9.00
-12.00
Disolución
0 100 200 300 400 500 600 700 800 900 1000 1100 1200
CE (μS/cm)
Calcita
119
4.00
2.00
1.00
SFC P1 ZAP2G SMI
Equilibrio GQ LE
SPM
SJO2
PO SF ERCO
P4DM MTS FB
0.00
CA LO2 LN
CA
-1.00
-2.00
-3.00
Disolución
-4.00
0 100 200 300 400 500 600 700 800 900 1000 1100 1200
CE (μS/cm)
Fluorita
4.00
Índice de saturación de Fluorita
Precipitación
3.00
2.00
1.00
Equilibrio LN
0.00
MTS FB
-1.00 SF ERCO
CA SJO2
LE P4
P1 ZA
GQ SPM
PO
-2.00 CA SFC DM
LO2
P2G
SMI
-3.00
Disolución
-4.00
0 100 200 300 400 500 600 700 800 900 1000 1100 1200
CE (μS/cm)
120
Análisis estadístico y calidad del agua
Tabla 23. Resumen estadístico de la composición química de las muestras del agua subterránea.
No. de Desviación NORMA
Parámetro Unidades Mínimo Máximo Media
muestras estándar 127
Temperatura ºC 19 27.3 43.7 35.48 16.21
pH 19 6.97 9.06 7.7 0.40 6.5 - 8.5
C.E. μS/cm 19 472 1090 682.95 187.84
OD mg/L 19 0.8 6.4 4.021 1.47
SDT mg/L 20 133.5 451.8 223.5 82.81 1000
Eh mV 15 -370 519.1 171.2 213.28
Dureza mg/L 20 38.97 280.7 140.09 64.06 500
2+
Ca mg/L 21 0.8 60.27 34.26 15.16
K+ mg/L 21 5.23 18.04 11.649 3.44
2+
Mg mg/L 21 0.36 31.61 11.668 8.06
Na+ mg/L 21 49.51 201.4 83.617 43.80 200
-
HCO3 mg/L 21 58.92 336.7 228.64 59.71
2-
SO4 mg/L 21 15.7 93.9 45.20 21.58 400
Cl- mg/L 21 5.63 121.2 29.95 31.57 250
3-
NO mg/l 21 0.28 56.4 22.52 16.55 10
Si4+ mg/L 21 22.26 100.331 38.34 14.83
Sr2+ mg/L 21 0.057 0.723 0.319 0.17
Fe mg/L 13 0.002 1.102 0.176 0.35 0.3
B mg/L 21 0.005 3.406 0.496 1.01
F- mg/L 21 0.161 23.090 2.166 5.35 1.5
3+
Al mg/L 21 0.004 0.342 0.028 0.07 0.2
As mg/L 1 0.054 0.054 0.054 0.00 0.025
Cu2+ mg/L 12 0.002 0.090 0.018 0.03 2
Pb2+ mg/L 20 0.008 0.023 0.011 0.004 0.01
Mn2+ mg/L 3 0.013 0.036 0.027 0.01 0.15
Zn2+ mg/L 21 0.064 0.621 0.342 0.16 5
2+
Ba mg/L 21 0.001 0.129 0.016 0.03 0.7
Li+ mg/L 21 0.026 1.399 0.231 0.32
Be2+ mg/L 2 0.001 0.001 0.001 0.0001
2+
Cd mg/L 2 0.387 0.553 0.470 0.12 0.005
2+
Ni mg/L 21 0.010 0.045 0.018 0.01
V mg/L 18 0.020 0.042 0.029 0.01
S mg/L 21 5.234 275.271 28.73 56.99
5.4.1 Parámetros fisicoquímicos
121
Los resultados de los parámetros fisicoquímicos se analizaron mediante diagramas
estadísticos de cajas y bigotes para identificar los datos atípicos, mediana y valores
mínimos y máximos. En la Figura 55, se visualizan los diagramas para cada parámetro
indicando la dispersión de los datos y características del agua subterránea.
122
Leyenda
Media x
Figura 55. Diagramas de cajas y bigotes para los parámetros fisicoquímicos de las 21
muestras de agua (Fuente: elaboración propia).
123
5.4.3 Elementos minoritarios
Figura 57. Diagramas de cajas y bigotes para los elementos minoritarios de las 21
muestras de agua subterránea (Fuente: elaboración propia).
Figura 58. Diagramas de cajas y bigotes de los elementos traza de las 21 muestras de
agua subterránea (Fuente: elaboración propia).
124
Índice de la calidad del agua
El 62% de las muestras de agua tiene un rango del 80-94 lo que significa que están
en una buena condición a pesar de que presentan algunas amenazas o daños de poca
magnitud; el 19% de los pozos muestreados presentan una condicion regular lo que
indica que el agua es ocasionalmente amenazada, contaminada o dañada; y un 19% de
los pozos muestreados (entre ellos el manantial termal) se encuentran frecuentemente
contaminados, amenazados o dañados. Existen diversos factores que pudieran afectar
la calidad del agua, entre los factores que pudieran influir son las descargas de aguas
residuales y la fuente del manantial termal Salitre.
125
Tabla 25. Resultados de los índices de calidad de los pozos de agua muestreados .
Número
Número Número real Número Número de CCM
Alcance Frecuencia Amplitu de
Pozo total de de variables de pruebas E
(F1) (F2) d (F3) pruebas WQI
variables probadas pruebas aprobadas
fallidas
S. Miguel
Ixtla
23.1 23.1 27.3 17 13 13 3 10 75.4
Granja
Quanni
21.4 21.4 4.2 17 14 14 3 11 82.3
Fuentes
Balvanera
23.1 23.1 3.5 17 13 13 3 10 81.0
Calera
Obrajuelo
16.7 16.7 12.6 17 12 12 2 10 84.6
Punta de
Obrajuelo
21.4 21.4 12.5 17 14 14 3 11 81.1
Mtal termal
Salitre
45.5 45.5 60.5 17 11 11 5 6 49.0
San Pedro
Mártir 1
14.3 14.3 10.0 17 14 14 2 12 87.0
Ejido
Romeral
15.4 15.4 16.1 17 13 13 2 11 84.4
Pueblito 2
Garza
16.7 16.7 27.2 17 12 12 2 10 79.2
S. José de
los Olvera 2
16.7 16.7 14.0 17 12 12 2 10 84.2
S. Fco
Corregidora
7.7 7.7 1.7 17 13 13 1 12 93.6
126
La Figura 59, muestra el mapa de la distribución espacial de la calidad del agua de
los pozos, donde evidentemente la muestra del agua del manantial termal Salitre se
encuentra en una condición contaminada debido al origen geotermal, sin embargo, el
pozo la Norita, Balvanera y el Castillo se encuentran influenciados por flujos regionales
provenientes del manantial termal presentando correlaciones en aluminio, cobre, plomo,
manganeso y cadmio.
Figura 59. Mapa de la distribución de la calidad del agua subterránea de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales INEGI, 2010 a y b;).
127
En la Tabla 26, se muestran los resultados de los 17 parámetros analizados y comparados con la norma NOM-127-
SSA1-1994 y la OMS. De acuerdo al análisis hidrogeoquímico del agua subterránea se tiene evidencias de contaminación
de origen antropogénico y la modificación química de fuentes naturales de origen geotermal debido a la zona de estudio.
Tabla 26. Evaluación de la calidad del agua subterránea con respecto a los criterios de la NOM-127-SSA1-1994 y la OMS.
SDT Dureza Na+ SO42- Cl- NO3- Fe B3+ F- Al3+ As Cu2+ Pb2+ Mn2+ Zn2+ Ba2+ Cd2+
No POZO Uso pH
mg/l mg/l mg/L mg/L mg/L mg/l mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L
1 San Miguel Ixtla Agua Potable 7.4 236.9 280.7 49.51 18.6 11.68 56.4 <l.c. 0.027 0.161 0.007 <l.c. 0.002 0.011 <l.c. 0.513 0.011 <l.c.
Castillo 150m Riego 170 83.82 87.53 27.8 8.12 4.18 0.479 1.024 0.437 0.047 <l.c. 0.090 0.018 0.0131 0.154 0.019 0.387
-
2
Castillo 200m Riego 146.8 51.27 80.56 26.2 10.67 0.28 1.102 1.142 0.693 0.018 <l.c. 0.089 0.016 0.0325 0.206 0.016 0.553
-
3 Granja Quanni Agua Potable 7.9 177 38.97 89.23 42.1 11.33 11.5 0.013 1.116 0.835 0.021 <l.c. 0.003 0.013 <l.c. 0.064 0.004 <l.c.
4 La Estancia Agua Potable 7.7 181.8 127.3 53.3 32.7 15.58 22.4 <l.c. 0.356 0.516 0.005 <l.c. <l.c. 0.010 <l.c. 0.277 0.005 <l.c.
5 Fuentes Balvanara Agua Potable 7.5 451.8 223.3 121.7 93.9 121.2 9.74 0.006 11.370 1.848 0.007 <l.c. <l.c. 0.011 <l.c. 0.236 0.010 <l.c.
6 Calera Obrajuelo Agua Potable 7.4 295 158.9 89.96 82 45.05 25.8 <l.c. 2.579 0.748 0.005 <l.c. <l.c. 0.009 <l.c. 0.547 0.006 <l.c.
7 Punta de Obrajuelo Agua Potable 7.5 166.3 123.8 57.61 35.4 11.13 27.9 0.002 0.553 0.424 0.008 <l.c. 0.020 0.011 <l.c. 0.187 0.010 <l.c.
8 La Norita Agua Potable 7.0 413.5 49.01 201.4 53.5 92.98 0.5 0.637 17.840 11.360 0.019 0.0538 <l.c. <l.c. 0.0358 0.123 0.129 <l.c.
9 Don Manuel Agua Potable 7.7 220.2 163.7 64.67 55.3 22.63 30 0.010 0.471 0.312 0.015 <l.c. 0.002 0.010 <l.c. 0.101 0.020 <l.c.
10 Manantial termal Salitre Manantial 9 - - 200.4 40.6 89.68 0.29 0.012 19.480 23.090 0.342 <l.c. <l.c. 0.023 <l.c. 0.295 0.002 <l.c.
11 Santa Fe Agua Potable 7.5 223.6 185.3 73.47 40.2 27.32 23.4 0.004 0.360 0.774 0.010 <l.c. 0.002 0.012 <l.c. 0.621 0.019 <l.c.
12 Poniente 4 Agua Potable 7.6 203.2 168.4 64.52 37.5 24.31 22.6 0.009 0.315 0.587 0.015 <l.c. 0.002 0.010 <l.c. 0.506 0.005 <l.c.
13 San Pedro Mártir 1-BIS Agua Potable 7.5 165.1 126.4 57.36 35.7 9.96 24 0.006 0.303 0.430 0.017 <l.c. 0.003 0.008 <l.c. 0.388 0.018 <l.c.
14 Ejido Romeral Agua Potable 7.8 283.7 192.2 78.45 90.5 25.47 32.9 0.003 0.776 0.648 0.009 <l.c. <l.c. 0.010 <l.c. 0.357 0.010 <l.c.
15 Pueblito 2 Garza Agua Potable 7.7 210 219.7 50.68 32.6 22.42 52.9 <l.c. 0.184 0.221 0.007 <l.c. <l.c. 0.009 <l.c. 0.411 0.010 <l.c.
16 San José de los Olvera 2 Agua Potable 7.7 178.3 118.5 52.76 49.1 16.03 27.6 <l.c. 0.240 0.573 0.005 <l.c. <l.c. 0.010 <l.c. 0.329 0.003 <l.c.
17 Los Olvera 2 Agua Potable 7.8 232.4 58.92 112.2 59.9 21.83 44.9 0.002 0.760 0.481 0.007 <l.c. <l.c. 0.009 <l.c. 0.405 0.001 <l.c.
18 San Francisco Corregidora Agua Potable 8.0 133.5 118.6 54.76 15.7 5.63 3.51 <l.c. 0.093 0.256 0.009 <l.c. 0.003 0.010 <l.c. 0.432 0.022 <l.c.
19 Poniente 1 Agua Potable 7.7 170.1 146.1 59.05 29.4 13.13 16.2 <l.c. 0.285 0.536 0.004 <l.c. 0.002 0.011 <l.c. 0.584 0.010 <l.c.
20 Zapata Agua Potable 7.8 211.3 167 56.84 50.6 22.88 35.9 <l.c. 0.281 0.548 0.005 <l.c. 0.002 0.010 <l.c. 0.437 0.008 <l.c.
NOM-127-SSA1-1994 Limite permisible 6.5 - 8.5 1000 500 200 400 250 10 0.3 1.5 0.2 0.025 2 0.01 0.15 5 0.7 0.005
OMS-2017 Limite permisible 6.5 - 8.5 1000 100 200 250 250 50 0.3 2.4 1.5 0.2 0.01 2 0.01 0.4 3 0.7 0.003
128
Análisis estadístico bivariado
5.6.1 Coeficiente de correlación
Mg2+ -0.17 0.31 0.36 0.25 0.49 0.45 0.97 0.90 0.37 1.00
Na+ 0.62 0.04 -0.67 0.75 -0.56 0.19 -0.56 -0.65 0.20 -0.42 1.00
HCO3
- 0.08 -0.05 -0.24 0.41 -0.30 0.15 0.35 0.30 0.45 0.41 0.17 1.00
SO42- 0.16 0.26 0.18 0.56 0.19 0.65 0.16 0.11 0.42 0.22 0.31 -0.25 1.00
Cl- 0.54 0.23 -0.45 0.93 -0.27 0.51 -0.08 -0.19 0.58 0.07 0.79 0.28 0.57 1.00
NO3- -0.09 0.37 0.48 -0.11 0.67 0.11 0.65 0.66 -0.13 0.59 -0.51 -0.16 0.07 -0.35 1.00
Si4+ 0.75 0.19 -0.49 0.30 -0.30 -0.57 -0.33 -0.38 0.00 -0.24 0.46 0.13 -0.07 0.34 -0.23 1.00
Sr2+ -0.21 0.22 0.34 0.32 0.50 0.56 0.87 0.81 0.47 0.89 -0.32 0.39 0.30 0.20 0.43 -0.31 1.00
Fe -0.40 -0.82 -0.37 -0.07 -0.72 0.06 -0.43 -0.41 -0.22 -0.44 0.30 0.12 -0.19 0.05 -0.51 -0.21 -0.34 1.00
B3+ 0.68 0.14 -0.70 0.83 -0.49 0.20 -0.40 -0.51 0.41 -0.25 0.94 0.34 0.28 0.89 -0.51 0.51 -0.13 0.21 1.00
F- 0.83 0.19 -0.71 0.62 -0.49 -0.23 -0.52 -0.60 0.21 -0.39 0.85 0.29 0.03 0.64 -0.44 0.77 -0.37 0.12 0.89 1.00
Al3+ 0.78 0.11 -0.61 0.37 -0.43 -0.53 -0.49 -0.55 0.00 -0.38 0.63 0.17 -0.08 0.42 -0.37 0.95 -0.39 -0.01 0.66 0.90 1.00
As 0.14 0.00 -0.39 0.51 -0.34 0.49 -0.28 -0.31 0.35 -0.23 0.62 0.35 0.09 0.46 -0.30 -0.25 -0.13 0.43 0.59 0.39 -0.03 1.00
Cu2+ -0.58 -0.98 -0.19 -0.39 -0.61 -0.22 -0.32 -0.25 -0.44 -0.37 -0.04 0.02 -0.32 -0.26 -0.40 -0.06 -0.28 0.78 -0.13 -0.13 0.00 -0.09 1.00
Pb2+ 0.22 -0.33 -0.36 -0.15 -0.43 -0.67 -0.30 -0.31 -0.38 -0.25 0.11 -0.01 -0.24 -0.04 -0.28 0.76 -0.33 0.09 0.10 0.35 0.69 -0.58 0.46 1.00
Mn2+ -0.27 -0.65 -0.42 0.13 -0.68 0.23 -0.44 -0.43 -0.05 -0.42 0.45 0.22 -0.12 0.20 -0.51 -0.27 -0.31 0.95 0.38 0.23 -0.02 0.70 0.58 -0.14 1.00
Zn2+ 0.01 0.31 0.58 -0.03 0.27 -0.06 0.50 0.51 0.08 0.46 -0.34 0.13 -0.08 -0.18 0.41 0.01 0.21 -0.41 -0.32 -0.20 -0.13 -0.30 -0.36 -0.08 -0.43 1.00
Ba2+ 0.02 -0.10 -0.33 0.44 -0.33 0.48 -0.22 -0.24 0.37 -0.19 0.53 0.43 0.01 0.38 -0.35 -0.32 -0.04 0.48 0.52 0.31 -0.09 0.97 0.01 -0.57 0.73 -0.31 1.00
Li+ 0.73 0.07 -0.71 0.65 -0.59 -0.12 -0.54 -0.61 0.24 -0.42 0.90 0.33 0.08 0.70 -0.53 0.70 -0.37 0.24 0.92 0.98 0.84 0.48 -0.02 0.32 0.36 -0.26 0.41 1.00
Be2+ 0.75 0.16 -0.71 0.64 -0.50 -0.10 -0.53 -0.60 0.26 -0.41 0.88 0.34 0.02 0.64 -0.44 0.60 -0.37 0.21 0.91 0.97 0.77 0.61 -0.13 0.15 0.38 -0.25 0.53 0.97 1.00
Cd2+ -0.56 -0.97 -0.20 -0.36 -0.65 -0.20 -0.32 -0.27 -0.44 -0.37 0.00 -0.01 -0.28 -0.21 -0.41 -0.08 -0.30 0.85 -0.10 -0.10 0.01 -0.07 0.97 0.45 0.64 -0.31 0.01 0.00 -0.10 1.00
Ni2+ 0.63 0.00 -0.50 0.18 -0.41 -0.65 -0.37 -0.41 -0.11 -0.29 0.40 0.14 -0.17 0.25 -0.30 0.98 -0.36 -0.06 0.45 0.73 0.94 -0.29 0.12 0.85 -0.16 -0.03 -0.35 0.66 0.55 0.11 1.00
V -0.21 0.38 0.51 -0.21 0.72 0.24 0.50 0.47 0.01 0.49 -0.56 -0.33 0.30 -0.26 0.46 -0.34 0.51 -0.68 -0.53 -0.64 -0.51 -0.48 -0.37 -0.25 -0.71 0.13 -0.46 -0.70 -0.69 -0.46 -0.42 1.00
1.0
S 0.85 0.24 -0.55 0.46 -0.32 -0.44 -0.42 -0.49 0.09 -0.30 0.65 0.11 0.08 0.49 -0.29 0.94 -0.33 -0.10 0.69 0.90 0.98 -0.04 -0.13 0.62 -0.10 -0.07 -0.12 0.83 0.77 -0.11 0.90 -0.44
0
129
Se realizó una matriz de coeficientes de correlación para los 33 parámetros
analizados en el agua subterránea del área de estudio. Se consideró una correlación
positiva a los coeficientes mayores a 7 (r>0.7) resaltados en color rojo y para las
correlaciones negativas se establecieron los coeficientes superiores a -0.7 y son
resaltados en color azul, estos valores representan una buena correlación entre las
variables (Tabla 27).
1000 B
800 C
LO2
(μS/cm)
600 GQ
400
A
200
0
0 20 40 60 80 100 120
Temperatura (°C)
130
Las correlaciones mayores a r>0.7 que ayudaron en la interpretación de la
evolución y procesos hidrogeoquímicos del agua subterránea son: la relación de la
conductividad eléctrica y los sólidos disueltos totales (CE/SDT), esto es evidente debido
al aumento de las concentraciones iónicas aumenta la conductividad eléctrica (Figura
61).
500
Solidos disueltos totales (mg/L)
FB
400 LN
300 ELCO
SMI
200
CA SFC
100
r=0.97
0
0 200 400 600 800 1000 1200
Conductividad Eléctrica (μS/cm)
131
5.6.2 Índices hidrogeoquímicos
En la Figura 62, se observa que las muestras del agua subterránea presentan un
coeficiente de correlación positiva (r=0.86) de Ca/Mg y una relación mayor a 1. El
incremento del calcio puede relacionarse a procesos de disolución de rocas
carbonatadas. Para el caso del magnesio se debe a la influencia de terrenos ricos en
silicatos magnésicos como los gabros y basaltos.
3.0
1:1
SMI
2.5
FB
Mg2+ (meq/L)
2.0 r= 0.86
ER
1.5 SF P2G
P4
DM
1.0
SFC
PO P1 COZA
SPM
SJO2 LE
0.5
LN LO2 CA
MTS GQ CA
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
Ca2+ (meq/L)
Figura 62. Relación iónica calcio/magnesio (rCa+2/rMg+2) de las muestras del agua.
132
En la Figura 63, se muestra la relación de sodio/calcio, en el cual el 90% de las
muestras están por arriba de la relación a 1 y para la muestra de San Miguel Ixtla (SMI)
y Pueblito 2 Garza (P2G) tienen una relación menor a 1, lo que significa que son muestras
de flujos locales.
3.5
1:1
SMI
3.0
P2G
2.5
FB
Ca2+(meq/L)
ZADM SFER
P4 CO
2.0
SJO2 P1 Enriquecimiento
LE
PO
1.5 SPM SFC de sodio
CA
1.0 LO2
r= -0.651 LN
0.5 GQ
MTS
0.0
0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0
Na+ (meq/L)
Figura 63. Relación iónica sodio/calcio (rNa+/rCa+2) de las muestras del agua.
133
6.0
SMI 1:1 r= -0.55
5.0
P2G FB
Ca+2 + Mg+2 (meq/L)
4.0 SFER
P4
ZADM CO
3.0 P1
SPM
LE
PO
SJO2
SFC
2.0
CA LO2
LN
1.0 GQ
MTS
0.0
0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0
Na+ (meq/L)
Figura 64. Relación iónica sodio/calcio (rNa+/rCa+2+rMg+2) de las muestras del agua.
4.00
FB r=0.56
3.50
3.00 MTS LN
rCl+2(meq/L)
2.50 1:1
2.00
CO
1.50
SF ER
1.00 P2GP4 ZA DMLO2
SMI CA LE
P1 SPM SJO2
0.50 SFC CA PO GQ
0.00
0.00 0.50 1.00 1.50 2.00 2.50
rSO42-(meq/L)
Figura 65. Relación iónica sulfato/cloro (rSO42-/rCl+2) de las muestras del agua.
134
En la Figura 66, se presenta la relación del cloro y bicarbonato (Cl/HCO), que de
acuerdo a los resultados se puede determinar que mediante este índice hidrogeoquímico
la relación iónica va de los 0.04 a 0.82 y se ubica en un rango de 0<rCl/rHCO<5, lo que
significa que no existe la evidencia de intrusión salina en el acuífero, lo que con lleva a
que las muestras del agua subterránea tienen un origen continental. A pesar de ello se
pueden diferenciar 2 tipos de aguas. En el grupo A se aprecia aguas ricas en bicarbonato
y para el grupo B se presentan aguas ricas en cloruros y se encuentra el pozo La Norita
(LN), Fuentes Balvanera (FB) y el manantial termal Salitre (MTS). Estas muestras están
influenciadas por un flujo regional que asciende desde las profundidades a través de
fallas y fracturas presentes en el área de estudio.
6
SMI
LN
5
CA P1 P2G SF MTS
SFC P4 FB
HCO3-(meq/L)
4 SPM ER
PO DM
CA ZA CO 1:1
3 GQ LE
SJO2
2
1 LO2
0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0
rCl+2(meq/L)
Figura 66. Relación iónica cloro y bicarbonato (rCl+2/rHCO42-) de las muestras del agua.
135
4.0
1:1
3.5 FB
3.0
LN
Cl- (meq/L)
2.5
MTS
2.0
1.5 CO
1.0 P4 SF ER
P2GZADM LO2
SJO2
LE
SMI P1
0.5 PO
SPM
SFC GQ
CACA
0.0
0.0 2.0 4.0 6.0 8.0 10.0
Na+(meq/L)
Figura 67. Relación iónica cloro y sodio (rCl-/rNa+) de las muestras del agua.
60
P2G SMI
50 LO2
40 ZA
NO3- (mg/L)
EL
DM
PO
SJO2
30 SPM CO
LE P4 SF
20 P1
GQ FB
NOM 127
10
CA
MTS LN
CA SFC
0
0 200 400 600 800 1000 1200
Conductividad Electrica
Figura 68. Relación iónica de la conductividad eléctrica y nitratos (rCE/rNO 3-) de las
muestras del agua.
136
El índice de litio/boro permite determinar posibles fuentes hidrotermales. En la
Figura 69, se observa como punto de referencia el manantial termal Salitre (MTS), que
se encuentra modificando a la composición química del agua subterránea de los pozos
La Norita (LN) y Fuentes Balvanera (FB). Las concentraciones de boro van desde 0.01
meq/L a 5.41 meq/L, asociados a fuentes geotermales presentes en la zona suroeste del
área de estudio, estos flujos regionales están modificando las características químicas
de los pozos LN y FB.
6.0
5.0 LN MTS
Flujo Regional
4.0
B+3(meq/L)
FB
3.0
Flujos intermedios
2.0
ZA
SFCER CO
1.0 GQCACA
PODM
LE
SPM
SJO2
P2G
SMI P4SF
Flujos locales
0.0 LO2 P1
0.00 0.05 0.10 0.15 0.20 0.25
Li+(meq/L)
Figura 69. Relación iónica litio y boro (rLi+/rB+3) de las muestras del agua.
137
25
MTS
20 Flujo Regional
de origen termal
F-(mg/L)
15
LN
10
Flujos locales o
intermedios
5
FB NOM 127
LEGQ
SMIPO CA
CA
0
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6
Li+(mg/L)
Figura 70. Relación iónica litio y flúor (rLi+/rF-) de las muestras del agua.
138
Análisis de los parámetros fisicoquímicos e hidrogeoquímicos.
En este apartado, se describen los resultados de los parámetros fisicoquímicos de
temperatura, pH, conductividad eléctrica, Eh, solidos disueltos totales y oxígeno disuelto
de las 21 muestras de agua tomadas in situ, de las cuales 19 corresponden a pozos y un
manantial termal. En la Tabla 18, se muestran los resultados de los parámetros
fisicoquímicos y en la Figura 46, se ubican geográficamente los datos.
5.7.1 Temperatura
Debido a la presencia de fuentes hidrotermales y por el gradiente geotérmico, la
temperatura del agua subterránea es muy variable, tanto espacial como temporal. En el
área de estudio se reportaron temperaturas que van de los 27.3°C a los 43.7°C y se
clasifican como aguas moderadamente tibias a calientes a excepción de la muestra del
manantial termal Salitre (MTS) que se reportó la mayor temperatura por arriba de los 47
ºC y que se clasifican como aguas muy calientes (Figura 71).
100
MTS
90
80
Temperatura (ºC)
70
60
LN
GQ
LO2
SJO2
50
SPM
P2G
CO
SMI
FB
SFC
PO
DM
P1
EL
SF
P4
40
LE
30
ZA
10
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
No. de Pozos
139
Figura 72. Mapa de la variación espacial de los datos de temperatura obtenidos de los
pozos muestreados (Fuente: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
140
14.0
MTS
SJO2
SFC
LO2
P2G
SPM
GQ
DM
EL
LE
PO
CO
P1
ZA
SMI
FB
SF
P4
pH (Log10)
LN
7.0
0.0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
No. de Pozos
141
5.7.3 Potencial de óxido reducción (Eh)
Los 13 pozos restantes, presentan valores positivos que van desde los 105 mV
hasta los 519 mV. Esto indica que las condiciones son oxidantes o existe perdidas de
electrones. Estos pozos se localizan en zonas alejadas del manantial termal Salitre o
alejados de la zona de falla Obrajuelo (Figura 76).
600
CO
SPM
P2G
SF
SJO2
400
EL
SFC
P4
PO
LO2
ZA
P1
DM
200
Eh (mV)
0
LN
MTS
-200
-400
-600
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
No. de Pozos
142
Figura 76. Mapa de la variación espacial de los datos de Eh obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
Sin embargo, se puede notar una gran diferencia en el pozo Fuentes Balvanera
(FB), La Norita (LN) y el manantial termal Salitre (MTS) que presentaron conductividades
de 1090 μS/cm 413.5 μS/cm y 989 μS/cm respectivamente. El aumento de la
conductividad eléctrica depende de la interacción de agua y roca (disolución-
precipitación), pero también puede ser provocada por la mezcla de aguas o flujos
regionales. En la Figura 78 se muestra el mapa de la distribución espacial de las
muestras de la conductividad eléctrica del agua subterránea en la zona de estudio.
143
1200
FB
LN
MTS
1000
CO
SMI
P2G
EL
LO2
DM
SF
800
CE (μS/cm)
P4
SJO2
SPM
ZA
P1
SFC
PO
CA
LE
GQ
600
CA
400
200
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
No. de Pozos
Figura 78. Mapa de la variación espacial de los datos de CE obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
144
5.7.5 Oxígeno Disuelto (OD)
El agua subterránea presentó valores de oxígeno disuelto que oscilan entre 0.8
mg/L a 6.4 mg/L y una media de 4 mg/L (Tabla 23 y Figura 79). Estas concentraciones
de oxígeno disuelto en aguas subterráneas representan valores normales que están
entre 0 y 5 mg/L (Custodio y Llamas, 1983; Werner, 1996).
Se puede observar que en el pozo La Norita (LN) y el manantial termal Salitre (MTS)
presentaron concentraciones de oxígeno disuelto de 0.9 mg/L y 0.8 mg/L
respectivamente, estos valores se correlacionan con los datos del potencial redox. Esto
nos indica que el agua tiene mayor tiempo de residencia en el subsuelo y existe un
ambiente reductor. En la Figura 80 se presenta la distribución espacial del oxígeno
disuelto y se observa que las concentraciones menores a 1.0 mg/L se ubican al suroeste
del valle y en los límites de la falla Obrajuelo.
P2G
PO
SFC
LO2
6
Oxígeno Disuelto (mg/L)
SJO2
DM
SPM
EL
ZA
P4
5
CO
FB
SMI
LE
SF
P1
GQ
3
MTS
2
LN
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
No. de Pozos
145
Figura 80. Mapa de la variación espacial de los datos de OD obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
Los valores de los sólidos totales disueltos medidos en el agua subterránea van
desde los 133.5 mg/L a 451.8 mg/L y una media de 223.52 mg/L (Tabla 23 y Figura 81).
Estas concentraciones de SDT en el agua subterránea representan valores que están
por debajo de la NOM-127-SSA1-1994 que estable un límite permisible de 1000 mg/L.
El 86 % de las muestras oscilan entre el 133.5 mg/L a 295 mg/L, sin embargo, en
la zona suroeste del área de estudio el pozo Fuentes Balvanera (FB) y La Norita (LN)
presentaron valores de 451.8 mg/L y 413.5 mg/L respectivamente. La concentración de
SDT de estos pozos están por arriba de la media, posiblemente están influenciados por
el manantial termal Salitre y por la mezcla de flujos de tipo regional (Figura 82).
146
500
Solidos disueltos totales (mg/L)
FB
LN
400
CO
EL
LO2
SMI
300
P2G
DM
SF
SJO2
P4
ZA
SPM
GQ
LE
CA
PO
P1
SFC
CA
200
100
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
No. de Pozos
Figura 82. Mapa de la variación espacial de los datos de SDT obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
147
5.7.7 Iones mayoritarios
5.7.7.1 Sodio (Na+)
Las concentraciones de sodio van de los 49.5 mg/L a los 201.4 mg/L, con una media
de 83.6 mg/L. En la Figura 83, se muestra el mapa de la distribución espacial de sodio,
las mayores concentraciones se ubican en los pozos la Norita y Fuentes Balvanera
asociados a la presencia del manantial termal Salitre por el ascenso de flujos regionales.
Estas concentraciones de sodio en el agua subterránea representan valores que están
por el límite de la NOM-127-SSA1-1994 que establece un límite permisible de 200 mg/L.
Figura 83. Mapa de la variación espacial de los datos de sodio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
148
5.7.7.2 Potasio (K+)
Los valores de potasio van de los 5.23 mg/L a los 18.04 mg/L, con una media de
11.65 mg/L. Los valores más altos de potasio se localizan en la zona al suroeste del área
de estudio y se encuentra el pozo la Norita, Fuentes Balvanera y Ejido Romeral. Esto es
debido a la presencia geotermal o disolución de feldespatos alcalinos como la ortoclasa
y sanidino o material arcilloso de composición potásica. La distribución de potasio se
muestra en la Figura 84.
Figura 84. Mapa de la variación espacial de los datos de potasio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
149
5.7.7.3 Calcio (Ca2+)
Las concentraciones de calcio van de los 0.8 mg/L a los 60.27 mg/L, con una media
de 34.2 mg/L. Los valores más altos y representativos de calcio se localizan al noroeste
y sureste del área de estudio y representan aguas de reciente infiltración debido a la
poca actividad en el intercambio iónico y aún se conserva el calcio proveniente de la
infiltración, estas muestras son: San Miguel Ixtla y Pueblito 2 como se muestra en el
mapa de la Figura 85.
Figura 85. Mapa de la variación espacial de los datos de calcio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
150
5.7.7.4 Magnesio (Mg2+)
Los valores de magnesio van de los 0.36 mg/L a los 31.6 mg/L, con una media de
11.6 mg/L. Las concentraciones magnesio son menores a las del calcio, sin embargo,
tienen una íntima correlación. Esto por la posible disolución de rocas carbonatadas o
volcánicas (magnésicas o ferromagnésicos). En el mapa de la Figura 86, se muestra la
distribución espacial del magnesio.
Figura 86. Mapa de la variación espacial de los datos de potasio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
151
5.7.7.5 Cloruros (Cl-)
Las concentraciones de cloruros van de los 5.6 mg/L a los 121.2 mg/L, con una
media de 29.9 mg/L. Los valores más altos de cloruros no rebasan el límite permisible
de la norma NOM-127-SSA1-1994 que establece un límite permisible de 250 mg/L. Sin
embargo, las altas concentraciones se localizaron al suroeste del área de estudio. El
origen del ion cloruro es causado por la presencia de fuentes hidrotermales, en especial
la muestra del manantial termal Salitre, Fuentes Balvanera y la Norita (Figura 87). La
relación iónica de cloruro y bicarbonatos va de los 0.04 a 0.82 y se ubica en un rango de
0<rCl/rHCO<5, lo que significa que las muestras del agua subterránea tienen un origen
continental.
Figura 87. Mapa de la variación espacial de los datos de cloruros obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
152
5.7.7.6 Sulfatos (SO42-)
Las muestras del agua subterránea presentan valores de sulfatos que van de los
15.7 mg/L a los 93.9 mg/L, con una media de 45.2 mg/L. Estas concentraciones de
sulfatos no rebasan el límite permisible de la norma NOM-127-SSA1-1994 que establece
un límite de 400 mg/L. Las altas concentraciones de sulfatos registradas en el área de
estudio van de los 40 mg/L a los 93 mg/L, corresponden al pozo Fuentes Balvanera, Ejido
Romeral, Calera Obrajuelo, Don Manuel, la Norita, Zapata, los Olvera 2 y San José de
los Olvera 2. El posible origen de los sulfatos de debe a la gran actividad agrícola que
existe en la zona (Figura 88).
Figura 88. Mapa de la variación espacial de los datos de sulfatos obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
153
5.7.7.7 Bicarbonatos (HCO3-)
Figura 89. Mapa de la variación espacial de los datos de bicarbonato obtenidos de los
pozos muestreados (Fuente: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
154
5.7.7.8 Nitratos (NO3-)
Figura 90. Mapa de la variación espacial de los datos de nitratos obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
155
5.7.8 Iones minoritarios
5.7.8.1 Boro (B3+)
Las concentraciones de boro en el agua subterránea van de los 0.005 mg/L a los
3.4 mg/L, con una media de 0.5 mg/L. Se encontraron 3 muestras de agua donde las
concentraciones fueron las más altas de boro y que corresponden al manantial termal
Salitre, Fuentes Balvanera y la Norita que van de los 11.3 mg/L a los 19.4 mg/L (Figura
91). El boro tiene diversos orígenes, sin embargo, los valores por arriba de los 10 mg/L
son representativos para aguas termales de ambientes volcánicos que corresponde a la
zona suroeste del área de estudio.
Figura 91. Mapa de la variación espacial de los datos de boro obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
156
5.7.8.2 Fluoruros (F-)
Las concentraciones de fluoruro en el agua subterránea van de los 0.16 mg/L a los
23 mg/L, con una media de 2.1 mg/L. Se encontraron 3 muestras de agua que
sobrepasan el límite permisible de flúor de la norma NOM-127-SSA1-1994 que establece
un límite de 1.5 mg/L. Las altas concentraciones de flúor corresponden al manantial
termal Salitre con 23 mg/L, la Norita con 11.36 mg/L y el pozo Balvanera con 1.8 mg/L.
El posible origen del fluoruro está asociado a la disolución de minerales como la fluorita
(CaF2), fluoroapatita (Ca5(PO4)3F), criolita (Na3AlF6), villiaumita, (NaF), topacio (Al2(SiO4)
F2) y gases volcánicos. La Figura 92, se muestra la distribución espacial del fluoruro,
donde se puede observar en la zona suroeste los valores más altos de fluoruro que se
asocian al ascenso de flujos regionales a través de la falla Obrajuelo.
Figura 92. Mapa de la variación espacial de los datos de fluoruros obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
157
5.7.8.3 Hierro (Fe)
Figura 93. Mapa de la variación espacial de los datos de hierro obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
158
5.7.8.4 Litio (Li+)
El mapa de la Figura 94, se muestra la distribución espacial del litio, cuyos valores
van de los 0.2 mg/L a los 1.4 mg/L y una media de 0.2 mg/L.
Figura 94. Mapa de la variación espacial de los datos de litio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
159
5.7.9 Elementos traza
Figura 95. Mapa de la variación espacial de los datos de aluminio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
160
5.7.9.2 Arsénico (As)
Figura 96. Mapa de la variación espacial de los datos de arsénico obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
161
5.7.9.3 Cadmio (Cd2+)
Figura 97. Mapa de la variación espacial de los datos de cadmio obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
162
5.7.9.4 Plomo (Pb2+)
Figura 98. Mapa de la variación espacial de los datos de plomo obtenidos de los pozos
muestreados (Fuente propia: de acuerdo a datos vectoriales de INEGI, 2010 a y b).
163
Análisis estadístico multivariado
5.8.1 Análisis de conglomerados, clúster o jerárquico
Figura 99. Dendograma de las 21 muestras del agua subterránea del área de estudio.
164
A continuación, se describen los 5 grupos a partir de la posición de la línea
fenotípica que se dividió en el dendograma.
165
Son aguas bicarbonatadas cálcicas (Ca-HCO3-) y bicarbonatada sódica (Na-HCO3-)
representando un flujo local e intermedio respectivamente. Las profundidades de los
pozos están en 160, 160, 139, 161, 146, 175, y 51 metros respectivamente.
Grupo E: Para este grupo solamente se ubicó la muestra del pozo San Miguel Ixtla,
ubicada en la zona noroeste del área de estudio que representa la zona de recarga. Se
caracteriza por presentar la máxima concentración en dureza, calcio, magnesio,
bicarbonatos y nitratos. Presenta concentraciones mínimas de flúor, litio y valores
significativos de estroncio. Es un agua de tipo bicarbonatada cálcica (Ca-HCO3-)
evidenciando un agua de reciente infiltración o zona de recarga. La profundidad estimada
es de 30 m.
El agrupamiento de las muestras con base a su similitud química permitió identificar
la distribución espacial de los procesos hidrogeoquímicos. En el mapa de la Figura 100,
se muestra la distribución geográfica de los grupos determinados mediante el análisis de
conglomerados.
0
0
0
0
166
Análisis de la variabilidad temporal de la composición hidrogeoquímica del
agua subterránea
167
Leyenda
Media x
mg/L
Límite de
la NOM
127
Figura 101. Diagramas de cajas y bigotes de los elementos mayoritarios de los estudios
de 2007, 2008 y 2018 (Fuente: elaboración propia con base a datos de Serrano, 2007 y
Pérez, 2008).
168
mg/L
Leyenda
Media x
mg/L
Límite de
la NOM
127
Figura 102. Diagramas de blox-plot de los elementos minoritarios y traza para los años
2007, 2008 y 2018 (Fuente: elaboración propia con base a datos de Serrano, 2007 y
Pérez, 2008).
169
El grupo B, corresponden a la muestra del pozo el Castillo (150m) y los Olvera 2.
Este grupo se ha caracterizados por tener similitudes en los valores de la conductividad
eléctrica, sodio, potasio, hierro, cobre, plomo, manganeso, litio, cadmio, boro y zinc. Los
niveles del agua subterránea están a 150 m y 146 m respectivamente.
El grupo C, se encuentra el pozo Punta Obrajuelo, San Pedro Mártir y San José de
lo Olvera 2. Están caracterizados por presentar similitudes en los siguientes elementos:
dureza, calcio, potasio, magnesio, estroncio, hierro, cobre, plomo y litio. Presentan un
agua de tipo bicarbonatada sódica (Na-HCO3-), y se encuentran a unas profundidades
de 110, 116 y 134 metros respectivamente.
Aun nado a lo anterior, existen una serie de factores que determinan la composición
química del agua subterránea, por ejemplo: la distancia entre cada uno de los pozos, la
profundidad, la dirección del flujo subterráneo, la distribución espacial de los pozos, tipo
de roca, la edad del agua subterránea, entre otras variables ambientales.
No. de Pozos
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
0
MTS SFC
20 Leyenda
40 SMI Grupos
CO
A
Profunidad (m)
60
B
80 GQ C
100 DM LN D
FB
120 PO SF SPM E
SJO2
140 P4
LO2 LE CA
160 P2G ZA P1 ER
180
CA
200
Figura 103. Gráfica de la relación de la profundidad y el análisis de conglomerados.
170
Origen de la composición química del agua subterránea
Tabla 28. Elementos de origen natural de acuerdo a las características del área de
estudio.
Elementos Origen
Concentraciones de iones
Son elementos que se encuentran de manera natural
mayoritarios: calcio, potasio,
debido a la disolución de rocas y minerales presentes
magnesio, sodio, bicarbonatos,
en el área de estudio.
sulfatos y cloruros.
Las concentraciones significativas de estos
elementos presentes en el área de estudio, son
Concentraciones de elementos
evidencias de fuentes hidrotermales provenientes de
minoritarios: Boro, Flúor, Hierro
flujos regionales y profundos y lo confirma la
y Litio
presencia del manantial termal Salitre al suroeste de
la zona de estudio.
El posible origen de estos elementos trazas están
Concentraciones de elementos
relacionados a flujos regionales de origen
traza: Aluminio, arsénico,
hidrotermal provenientes del manantial termal
cadmio y Plomo.
Salitre.
Tabla 29. Elementos que tienen un origen antropogénico de acuerdo a las características
del área de estudio.
Elementos Origen
El nitrato tiene un origen antropogénico, provocado
Concentraciones de ion por la intensa actividad agrícola provocado por la
mayoritario: Nitratos utilización de fertilizantes y por la posible infiltración
de aguas residuales provenientes de la superficie.
171
Análisis Hidrogeológico
B B'
A'
Figura 104. Mapa hidrogeológico (Fuente: elaboración con base a datos de Alaniz-
Álvarez et al., 2001 y Consejo de Recursos Minerales, 1999).
172
En la Figura 105 y Figura 106, se muestran los perfiles litológicos de 17 pozos y un
manantial termal, en los cuales, se describe el nivel de profundidad del agua subterránea
expresado en metros sobre el nivel del mar (msnm). El acuífero fluye a través de rocas
volcánicas (basaltos) y rocas granulares (arenas). En cuanto a la composición química
de los pozos tenemos la presencia de 7 parámetros (As, Cd2+, F-, Fe, Na+, NO3- y Pb2+)
que sobrepasan los límites permisibles de la NOM-127-SSA1-1994.
Figura 105. Perfiles litológicos de los pozos de la zona de Obrajuelo (Fuente: elaboración con
base a datos de Alaniz-Álvarez et al., 2001; Carreón et al. 2016; Suez-Lesser, 2017; IMTA, 2018).
Figura 106. Perfiles litológicos de los pozos de la zona de Querétaro (Fuente: elaboración con
base a datos de Alaniz-Álvarez et al., 2001; Carreón et al. 2016; Suez-Lesser, 2017; IMTA, 2018).
173
La Figura 107 y Figura 108 se muestran la sección hidrogeológica de norte-sur y
oeste-este respectivamente. Se puede observar la dirección del agua subterránea y la
zona de abatimiento. La estratigrafía está compuesta principalmente por rocas
volcánicas, secuencias vulcanoclásticas, derrames andesíticos y basálticos.
174
Modelo Hidrogeológico conceptual
175
Figura 109. Modelo hidrogeológico del acuífero del Valle de Obrajuelo-Querétaro.
Instituto Mexicano de Tecnología del Agua (IMTA).
176
Discusiones
177
Mediante el análisis de los elementos mayoritarios de 21 muestras de agua, se
determinaron dos tipos de familias de agua: bicarbonatada cálcica (Ca-HCO3-) y
bicarbonatada sódica (Na-HCO3-). Lo que se interpreta que la evolución del agua
subterránea demuestra una secuencia normal de un flujo de reciente infiltración a un flujo
intermedio a regional. Por otro lado, las investigaciones de López (1996), Serrano (2007),
Pérez (2008) y Suez-Lesser (2017), determinaron las dos mismas familias de aguas;
bicarbonatada cálcica y bicarbonatada sódica. Por lo tanto, se deduce que el acuífero
permanece con la misma composición hidrogeoquímica, que significa no se han
presentado variaciones importantes en la composición mayoritaria del agua subterránea.
Para el análisis del índice de saturación, con base al tipo de rocas, se consideraron
cinco minerales presentes en el agua subterránea que corresponde a rocas volcánicas
de composición andesítica y basáltica, constituidas por los siguientes minerales:
plagioclasas (albita), feldespatos, micas, calcita y fluorita. Se pudo constatar que las
plagioclasas (albita) y la fluorita se encuentra en disolución, sin embargo, las muestras
de agua de La Norita y el Manantial termal el salitre presentan los feldespatos en
disolución. La mica de composición potásica se encuentra en precipitación y para la
calcita se encuentra en equilibrio en el acuífero.
178
En el análisis estadístico bivariado de coeficientes de correlación, se obtuvieron
resultados tanto positivos como negativos, los cuales permitieron observar procesos
hidrogeoquímicos que ocurren en la interacción agua-roca, procesos de mezclas de
aguas, disolución de minerales y observar el aporte de componentes de origen
antropogénico, así como la evolución del agua subterránea.
Los parametros mas significativos que determinaron las calidad del agua subterranea se
describen a continuación:
179
pozos de agua, estas concentraciones van de los 0.28 mg/L a los 53 mg/L rebasando en
algunas muestras de agua, los límites máximos permisibles de la NOM-127-SSA1-1994
que establece un límite permisible de 10 mg/L y que modifican la calidad del agua
subterránea, sin embargo, se han reportado concentraciones de hasta 117.1 mg/L para
el pozo Alameda de cuerdo a Serrano (2007). Para el estudio de Pérez (2008), su registro
máximo fue de 148 mg/L, pozo ubicado en la zona metropolitana de Querétaro. Estas
concentraciones se atribuyen por la intensa actividad agrícola e infiltración de aguas
residuales que ocurren en la zona.
Así mismo, las investigaciones realizadas por Serrano (2007) y Pérez (2008),
reportaron concentraciones de flúor por arriba de la norma NOM-127-SSA1-1994, sin
embargo, en la presente investigacion se encontraron valores que van de los 0.16 mg/L
a los 23 mg/L rebasando el limite permisible de 1.5 mg/L provocado por la mezcla de
flujos regionales provenientes del manantial termal Salitre.
De los datos recopilados de los años 2007, 2008 y comparándolos con los datos
del 2018, se realizó el análisis de la variabilidad temporal de la composición de los
elementos mayoritarios, en el cual, dio como resultado poca variabilidad en las medias
de las concentraciones de calcio, potasio, magnesio, sodio, bicarbonatos, sulfatos,
cloruros y nitratos. Es decir, la composición mayoritaria del agua subterránea se ha
mantenido.
180
Para los elementos minoritarios y traza más significativos presentes en el agua
subterránea, se analizaron para los años 2007, 2008 y 2018, que se describen a
continuación:
181
El grupo C, se encuentra el pozo Punta Obrajuelo, San Pedro Mártir y San José de
lo Olvera 2. Están caracterizados por presentar similitudes en los siguientes elementos:
dureza, calcio, potasio, magnesio, estroncio, hierro, cobre, plomo y litio. Presentan un
agua de tipo bicarbonatada sódica (Na-HCO3-), y se encuentran a unas profundidades
de 110, 116 y 134 metros respectivamente.
Aun nado a lo anterior, existen una serie de factores que determinan la composición
química del agua subterránea, por ejemplo: la distancia entre cada uno de los pozos, la
profundidad, la dirección del flujo subterráneo, la distribución espacial de los pozos, tipo
de roca, la edad del agua subterránea, entre otras variables ambientales.
Para los iones mayoritarios: calcio, potasio, magnesio, sodio, bicarbonatos, sulfatos
y cloruros; son elementos que se encuentran de manera natural debido a la disolución
de rocas y minerales presentes en el área de estudio.
Las concentraciones de los elementos minoritarios: boro, flúor, hierro y litio; son
evidencias de fuentes hidrotermales provenientes de flujos regionales y profundos y se
confirma con la presencia del manantial termal Salitre al suroeste de la zona de estudio.
182
Las concentraciones de nitratos, tienen un origen antropogénico, provocado por la
intensa actividad agrícola provocado por la utilización de fertilizantes y por la posible
infiltración de aguas residuales provenientes de la superficie.
183
CAPITULO VI.
6. CONCLUSIONES
La profundidad media del agua subterránea se encuentra a los 120 metros, las
máximas profundidades se registraron en los pozos: el Castillo, Ejido Romeral, Poniente
1 y Pueblito 2 Garza. Estos pozos se ubican en las localidades del Castillo en el Estado
de Guanajuato; Sta. María Magdalena y el Pueblito en el Estado de Querétaro. Lo que
se infiere que la dirección preferencial del flujo subterráneo es de tipo radial, con dirección
al centro del valle de Querétaro. Dicha explotación es inducida por la intensa actividad
agrícola, industrial y urbana de la región. De acuerdo al análisis hidrogeológico, los
sistemas de flujo del agua subterránea se encuentran circulando por dos unidades
litológicas: granular (arenas y conglomerados) y volcánicas fracturadas (basaltos y
andesitas).
184
Mediante el análisis estadístico permitió simplificar, agrupar y correlacionar las 33
variables fisicoquímicas e hidrogeoquímicas. Dichos resultados ayudaron a definir
procesos hidrogeoquímicos que ocurren en la interacción agua-roca, procesos de
mezclas de aguas, disolución de minerales y observar el aporte de componentes de
origen antropogénico.
Los resultados del cálculo del índice de saturación demuestran que el mineral la
albita y la fluorita se encuentran en disolución en el acuífero, mientras que el feldespato
potásico y la mica potásica se encuentran precipitadas en el agua subterránea. Para esta
investigación, se tomó en cuenta la geología del área de estudio, en este caso se
consideraron rocas volcánicas de composición andesítica y basáltica.
185
7. RECOMENDACIONES
Una de las recomendaciones planteadas para este estudio, es que las autoridades
correspondientes realicen monitoreos periódicos para vigilar la calidad del agua de los
pozos circundantes al manantial termal Salitre.
Para extraer agua de los mantos acuíferos de mejor calidad en el área de estudio,
es necesario tomar en cuenta las siguientes recomendaciones: no perforar en zonas de
abatimiento (no mayor a 150 m de profundidad); estar a una distancia considerable del
manantial termal Salitre y de fallas que tengan una orientación norte-sur por la posible
mezcla de flujos regionales.
186
8. BIBLIOGRAFÍA
Amézaga, B.S (2016). Actualización hidrogeológica del acuífero del Valle de Celaya.
Análisis y determinación por concentraciones de Arsénico; propuesta de sitio para
la perforación de un pozo de adsorción. Tesis de Licenciatura en Ingeniero Geólogo.
Universidad Nacional Autónoma de México. México.
Bear, J. y Cheng, A.H.D. (2010). Modeling Groudwater Flow and Contaminant Transport.
Springer.
Campa, M.F. (1983). The Tectonostratigraphic Terranes and the Thrust Belt in Mexican
Teritory. Proceedings of the Circum-pacific Terrane Conference, Stanford University
Publications Geological Sciences, 18, 44-46.
Campos, D. (1998). Procesos del ciclo Hidrológico. (3°ed.) México, D.F: Universidad
Autónoma de San Luis Potosí.
Carreón-Freyre D., Cerca, M., Luna, l., Ochoa-González G., Tatini P., Zuñiga F. R.
(2016). Shearing along faults and stratigraphic joints controlled by land subsidence
in the Valley of Queretaro, Mexico. Springer, 24, 657-674. DOI 10.1007/s10040-
016-1384-0.
187
CCME (2001). Canadian Water Quality Guidelines for the Protection of Aquatic Life:
CCME Water Quality Index 1.0, Technical Report, Canadian Council of Ministers of
the environment winnipeg, MB, Canada. [en línea].
http://www.ccme.ca/sourcetotap/wqi.html.
Centro Queretano de Recursos Naturales (CQRN), (2001). Uso del suelo y vegetación
de la zona Sur de Estado de Querétaro. Consejo de Ciencia y Tecnología del
Estado de Querétaro. Reporte Técnico.
Cervera, L.A. (2007). Indicadores de uso sustentable del agua en Ciudad Juárez,
Chihuahua. Scientific Electronic Library Online, 8.
Comisión Estatal del Agua de Querétaro (CEAQ), 2011. Zona Metropolitana de la ciudad
de Querétaro (ZMQ), México. Kit de capacitación SWITCH, estudio caso Querétaro
(ZMQ), México.
CONAGUA (Comisión Nacional del Agua) (2016). Estadísticas del Agua en México.
Edición 2016. Secretaría de Medio Ambiente y Recursos Naturales. México.
188
CONAGUA (Comisión Nacional del Agua) (2017). Estadísticas del Agua en México.
Edición 2016. Secretaría de Medio Ambiente y Recursos Naturales. México.
Consejo de Recursos Minerales (CRM), (1999). Monografía Geológico Minera del Estado
de Querétaro: Secretaría de Comercio y Fomento Industrial, Coordinación General
de Minería.
Córdoba, D.F., Meffe, R., Martínez, V., Segura, M. (2012). Metodología para el desarrollo
de un modelo hidrogeológico mediante herramientas de Información Geográficas.
Centro de investigación y Desarrollo en Recursos Hídricos (CIDERH), Chile.
Instituto Madrileño de Estudios Avanzados (IMDEA) Agua, España.
Custodio, E. y Llamas, M.R. (1983). Hidrología subterránea. (2a ed.). Barcelona: Omega.
Ferrari, L., Garduño, V.H., Innocenti, F., Manetti, P., Pasquarè, G., Vaggelli, G., (1994).
A widespread mafic volcanic unit at the base of the Mexican Volcanic Belt between
Guadalajara and Queretaro: Geofísica Internacional, 33, 107-124.
Ferrari, L. (2004). “Slab detachment control on mafic volcanic pulse and mantle
heterogeneity in central México”. En: Geology, 32 (1), 77-80.
189
Fetter, C. (1988). Applied Hydrogeology. New York: Company Mac Millam Publiling.
Foster S., Kemper, K., Garduño, H., Nanni, N., y Tuinhof. (2006). Gestión
sustentable del agua subterránea: Amenazas naturales a la calidad del agua
subterránea. Washington, D.C., United State of America: Banco mundial.
Freeze, A., & Cherry T. A. (1979). Groundwater. New Jersey, USA: Prentice - Hall.
Gleik, P., H., Sandra P. L. & Jason M. (1996), The Sustainable Use of Water in the Lower
Colorado River Basin, The Pacific Institute for Studies in Development, Environment
and security, Oakland, California.
Güler, C., Thyne, G., McCray, J., y Turner, A. (2002). Evaluation of graphical and
multivariate statistical methods for classification of water chemistry data.
Hydrogeology Journal, 10, 455-474.
Gutiérrez Carrillo, N., Palacios Vélez, E., & Peña Díaz, S., & Palacios Vélez, O. (2002).
Escenarios para el aprovechamiento sustentable del acuífero del Valle de
Querétaro. Agrociencia, 36 (1), 1-10.
Hem, J.D. (1985). Study and Interpretation of the Chemical Characteristics of Natural
Water, Alexandria, Virginia, U.S. Geological Survey.
190
Instituto Nacional de Estadística, Geografía e Informática (INEGI) (2010b). Compendio
de información geográfica municipal de los Estados Unidos Mexicanos Querétaro,
Querétaro clave geoestadística 22014. México.
Loucks, P. Daniel & John Gladwell S. (1999). Sustainability Criteria for Water Resource
Systems, International Hydrology Programme, Cambridge University Press.
Lugo, J. H. (2011). Diccionario geomorfológico. Geografía para el siglo XXI, serie: Textos
universitarios, 7 Instituto de Geografía UNAM, México, 480.
191
Mejía M., Martínez M., González L., Aparicio J.C. y Mañón D.S., 2018. Evidencias
isotópicas e hidrogeoquímicas de flujos profundos en el acuífero de Obrajuelo,
Guanajuato. XXV Congreso Nacional de Hidráulica. 1-4.
Melloul, A., Collin, M., (1992). The principal components statistical method as a
complementary approach to geochemical methods in water quality factor
identification; application to the Coastal Plain aquifer of Israel. Journal of Hidrology,
140, 49-73
Nieto-Samaniego A. F., Ferrari L., Alaniz-Álvarez S.A., Labarthe Hernández, G., Rosas-
Elguera, J., (1999). Informe de la Carta Geológico-Minera y Geoquímica Hoja
Celaya F14C64 escala 1:50000 Estado de Guanajuato. Pachuca, Hidalgo, México:
Consejo de Recursos Minerales (Servicio Geológico Mexicano).
Norma Oficial Mexicana (NOM), 2000, Modificación a la Norma Oficial Mexicana NOM-
127-SSA1-1994, Salud ambiental. Agua para uso y consumo humano. Límites
permisibles de calidad y tratamientos a que debe someterse el agua para su
potabilización: México, Secretaría de Salud, publicado en el Diario Oficial de la
Federación el 22 de noviembre de 2000.
Ortega, F., Mitre, L.M., Roldán, J., Aranda, J.J., Morán, D., Alaniz, S.A. y Nieto, A.F.
(1992). “Carta Geológica de la República Mexicana Escala 1: 2.000.000”.
Universidad Nacional Autónoma de México UNAM, Instituto de Geología y
Secretaría de Energía, Minas e Industria. 1 Mapa con texto explicativo, 74.
192
Paces, J.B., Ludwig, K.R., Peterman, Z.E., Neymark, L.A., (2002). 234U/238U evidence
for local recharge and patterns of groundwater flow in the vicinity of Yucca Mountain,
Nevada, USA. Applied Geochemistry, 17, 751-779.
Pérez-Venzor, J.A., Aranda-Gómez, J.J., McDowell, F., Solorio Munguía, J.G., (1996).
Geología del volcán Palo Huérfano, Guanajuato, México: Revista Mexicana de
Ciencias Geológicas, 13 (2), 174-183.
Pérez V., J. (2008). Caracterización Hidrogeoquímica del Agua Subterránea del Valle de
Querétaro para el Estudio de Sistemas de Flujo. Tesis de Maestría en Ciencia
Ambientales. Centro Universitario Santiago de Querétaro, Querétaro. México.
Shiklomanov. (1993). World fresh water resources, In: Peter HG (ed), Water in crisis: a
guide to the world’s fresh water resources, Oxford University Press, New York.
Tóth J, 1962. A theory of groundwater motion in small drainage basins in Central Alberta,
Canada. J. of Geophysical Res. 67(11), 4375-4387.
Tóth, J. (2000), "Las aguas subterráneas como agente geológico: causas, procesos y
manifestaciones", Boletín Geológico y Minero, 111, (4), 9–26
193
Tóth, J. (2009). Gravitational Systems of Groundwater Flow: Theory, Evaluation,
Utilization. Cambridge University Press.
194
9. ANEXOS
Índice de saturación
No Pozo ID
Albita Feldespato K Mica K Calcita Fluorita
1 San Miguel Ixtla SMI 0 0.14 6.34 0.26 -2.78
2 Castillo 150m CA 0.31 1.63 3.70 -0.89 -2.05
3 Castillo 200m CA -0.73 0.50 10.17 -0.58 -1.87
4 Granja Quanni GQ -1.16 -0.19 5.13 0.02 -2.01
5 La Estancia LE -1.27 0.40 5.96 0.07 -1.90
6 Fuentes Balvanara FB -1.08 0.34 6.07 0.15 -0.83
7 Calera Obrajuelo CO -1.38 0.00 5.73 -0.07 -1.58
8 Punta de Obrajuelo PO -1.34 0.23 6.26 -0.05 -2.15
9 La Norita LN -2.03 -0.98 5.17 -0.62 0.17
10 Don Manuel DM -1.01 0.55 6.79 0.23 -2.30
11 Manantial termal Salitre MTS -2.04 -1.63 0.24 -0.04 -0.89
12 Santa Fe SF -0.93 0.58 6.98 0.15 -1.50
13 Poniente 4 P4 -0.75 0.82 7.36 0.20 -1.75
14 San Pedro Mártir 1-BIS SPM -1.15 0.40 7.00 0.00 -2.14
15 Ejido Romeral ER -1.07 0.50 6.15 0.36 -1.69
16 Pueblito 2 Garza P2G -1.21 0.54 6.19 0.45 -2.49
17 San José de los Olvera 2 SJO2 -1.50 0.09 5.28 0.02 -1.89
18 Los Olvera 2 LO2 -1.53 -0.68 4.08 -0.49 -2.36
19 San Francisco Corregidora SFC -1.12 0.60 6.03 0.47 -2.54
20 Poniente 1 P1 -1.43 0.16 5.29 0.31 -1.86
21 Zapata ZA -1.32 0.31 5.49 0.32 -1.78
195