Introducción

En la zona sur del estado de Quintana Roo, México, al igual que en toda la península de Yucatán, el acuífero es la principal fuente de abastecimiento de agua para todas las actividades humanas; no obstante, es altamente vulnerable a la contaminación debido a sus características cársticas, que permiten infiltración y arrastre de contaminantes. Aunada a esta situación, las perspectivas de crecimiento demográfico y económico impulsado por la industria turística han traído consigo una creciente demanda de agua y, por consiguiente, una creciente generación de aguas residuales y desechos sólidos (^{CNA, 1997}).

Aparentemente, las características propias de la península deberían proporcionar abundantes suministros de agua; sin embargo, la combinación de factores, como el clima y la hidrogeolo-gía, han generado un sistema hidrológico en el que el agua dulce es escasa, teniéndose inclusive un suministro restringido en algunas regiones (^{Ward, Weidie, & Back, 1985}).

Dadas estas condiciones, se han empezado a desarrollar diversos programas de ordenamiento, con la firme intención de garantizar el abastecimiento de agua potable para zonas urbanas, así como establecer las bases para que el agua que se extraiga para uso público en zonas rurales se mantenga en condiciones adecuadas de calidad y en cantidad suficiente como para satisfacer la demanda de manera oportuna y permanente.

No obstante, para impulsar la gestión adecuada del agua es fundamental contar con información necesaria sobre el estado de este recurso, pues así será posible guiar la toma de decisiones y definir acciones. En este sentido, el conocimiento de la composición química y su distribución espacial a través de la calidad del agua es fundamental desde el punto de vista de su aptitud para el consumo humano, al igual que por la información que pueden aportar acerca de la dirección y extensión de los sistemas de flujo subterráneo (^{CCPY, 2012}); más aún, puede indicar cambios estacionales y temporales en la medida que sus constituyentes químicos estén presentes en el acuífero (^{Sajil & James, 2013}). El índice de calidad del agua (ICA) es una de las herramientas más efectivas para resumir los datos de calidad de las aguas subterráneas en un formato fácilmente expresable y comprensible, capaz de ser utilizado por los responsables políticos y ciudadanos interesados (^{Lateef, 2011}). Se define como una técnica de clasificación que engloba la influencia compuesta de cierto grupo de parámetros de calidad del agua en un solo valor general (^{Aly, Al-Omran, & Alharby, 2014}).

En el mundo se han desarrollado diferentes ICA para aguas superficiales y subterráneas, entre los que destaca el de la National Sanitation Foundation (ICA-NSF), desarrollado por ^{Brown, McClelland, Deininger y Tozer (1970)} para ríos de Estados Unidos, y ampliamente empleado en diferentes estudios internacionales. ^{Dinius (1987)} desarrolló un ICA similar, teniendo como valor agregado el planteamiento de rangos de clasificación basados en usos específicos, entre los que destaca el del consumo humano. ^{Backman, Bodis, Lahermo, Rapant y Tarvainen (1998)} crearon un índice para evaluar y asignar la cantidad de contaminación del agua subterránea, y lo aplicaron en el suroeste de Finlandia y Eslovaquia central. De forma más reciente se han planteado ICA para la evaluación del recurso hídrico a ser destinado al consumo humano, previo tratamiento, que incluye parámetros fisicoquímicos y microbiológicos relacionados con el nivel de riesgo sanitario presente en el agua. ^{Avvannavar y Shrihari (2008)} mostraron los resultados de su intento por desarrollar un índice mediante seis parámetros fisicoquímicos para evaluar la calidad del agua superficial del río Netravathi en India. ^{Saeedi, Abessi, Sharifi y Meraji (2010)} aplicaron ocho parámetros diferentes, incluyendo K+, Na+, Ca²+, Mg²+, SO₄^2-, Cl^-, pH y STD, como los componentes más importantes del agua para el desarrollo de un índice de calidad de aguas subterráneas en Irán. El trabajo de ^{Ketata-Rokbani, Gueddari y Bou-hlila (2011)} informa de la creación de ICA para las aguas subterráneas y los resultados de su aplicación, a fin de evaluar la calidad del agua subterránea en El-Khairat, acuífero profundo de Túnez.

En México, a principios de la década de 1970, se desarrolló un índice de calidad del agua compuesto por 18 parámetros como sistema cualitativo para hacer comparaciones de los niveles de contaminación en diferentes áreas. Este ICA define el grado de contaminación existente en el agua a la fecha de un muestreo, expresado como un porcentaje de agua pura, donde 0% del ICA representa agua altamente contaminada y 100%, agua en excelentes condiciones. Este ICA sigue vigente en la actualidad. No obstante, para fines exclusivos de evaluación, se ha decidido utilizar tres parámetros indicadores: demanda bioquímica de oxígeno (DBO5), demanda química de oxígeno (DQO) y sólidos suspendidos totales (SST), que muestran la influencia antropogénica desde el punto de vista de afectación por la presencia de centros urbanos e industriales.

Si bien este ICA incluye gran cantidad de parámetros indicadores elegidos a partir de las estaciones de medición de la Red Nacional de Monitoreo de la Comisión Nacional del Agua, es de mencionar que los distintos pesos relativos (W_i) asignados a cada parámetro para su sistematización no corresponden con las características geohidrológicas de cada región; no se toma en cuenta la injerencia de cada parámetro en su entorno, como puede ser la geología del suelo e incluso el uso de suelo. Se sustenta más bien en una estandarización lineal homogénea, tomando como base que el recurso hídrico para el abastecimiento de poblaciones se obtiene de fuentes superficiales y no subterráneas, como es el caso de la península de Yucatán.

En este estudio se propone un sistema estimativo de la calidad del agua subterránea a través del ICA basado en la medición física de siete parámetros fisicoquímicos. Cada uno de ellos con diferentes pesos relativos (W_i), enmarcados según el grado de contaminación del agua del acuífero y los diversos impactos sobre los usos de este recurso hídrico.

Se ilustra y explica además el comportamiento de siete parámetros fisicoquímicos, donde se pretende realizar una valoración cuantitativa y cualitativa de la calidad de las aguas subterráneas de manera espacial y temporal en el acuífero del sureste de Quintana Roo.

Materiales y métodos

Area de estudio

El área de estudio forma parte de dos subregio-nes hidrológicas (RH): la 33-A, conocida como Bahía de Chetumal y Otras, y la 33-B, conocida como Cuencas Cerradas (Figura 1). Se localiza geográficamente al sureste de la península de Yucatán, en las coordenadas 17° 49' a 20° 12' de latitud norte, y entre los 87° 30' y 89° 27' de longitud oeste, cubriendo una extensión total de 37 896 420 km². Colinda hacia el nororiente con el mar Caribe y hacia el suroriente por una línea que está definida por la bahía de Chetumal, por el límite nacional con Belice en un tramo de límite natural constituido por el río Hondo y enseguida por el límite político con Belice (^{UQROO, 2002}). Esta área tiene un notable contraste topográfico. Hacia el norte y la franja costera se observa un relieve topográfico plano con elevaciones que van de 0 a 20 msnm, mientras que al sur-suroccidente se observa un relieve ondulado con elevaciones que van de 50 a 340 msnm (^{Bauer-Gottwein, Gondwe, Char-vet, & Marín, 2011}). La temperatura promedio anual es de 26.15 °C, con temperaturas máximas entre 36 y 38 °C, y mínimas entre 12 y 14 °C (^{CCPY, 2012}). De acuerdo con la clasificación climática de Kóppen-Geiger, la región cuenta con un clima aw2 (i); cálido subhúmedo con abundantes lluvias en verano y parte de invierno (^{CNA, 2006}). La precipitación media anual es de 1 365 mm, aproximadamente 600 mm más alta que el promedio nacional, el cual es 772 mm (^{JICA, 2004}; ^{Alcocer, Lugo, Sánchez, & Esco-bar,1999}). La temporada de lluvias se presenta en los meses de mayo a octubre, con una lámina media mensual de 150 a 220 mm, mientras que la temporada de estiaje o secas comprende los meses de noviembre a abril, durante los cuales las láminas llovidas son menores que 60 mm (^{UQROO, 2002}).

Figura 1 Ubicación geográfica del área de estudio. La figura (a) muestra los sitios de monitoreo del año 2012; la figura (b) muestra los sitios de monitoreo del año 2002. El fondo en sombreado gris indica el Modelo Digital de Elevación (DEM), modificado de ^{INEGI (2014)}. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 

Marco hidrogeológico

El área de estudio se define como una unidad geológica cubierta casi en su totalidad de sedimentos calcáreos de origen marino del Terciario reciente (Butterlin & Bonet, 1960; ^{López-Ramos, 1972}). Se sabe por ^{Isphording (1975)} que la geología estructural de la región de estudio está conformada por una secuencia de rocas de edad Terciaria, constituida de evaporitas intercaladas con rocas calizas y minerales, como dolomita, aragonita, yeso y halita (Figura 2). Estas rocas carbonatadas son altamente permeables y solubles, por lo que al ser disueltas dejan escaso residuo, razón por la cual se presentan espesores de suelo muy reducidos, que a su vez resultan en la presencia de un sistema acuífero cárstico maduro con amplias fisuras, fracturas y cuevas (^{Beddows, Smart, Whitaker, & Smith, 2002}; ^{Hausman, 2009}; ^{JICA, 2004}).

Figura 2 Mapa geológico de la península de Yucatán, modificada de ^{SGM (2007)}. El contorno en negro representa el área de estudio. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q) 

Fisiográficamente, el área de estudio se ubica en la unidad "Planicie Central", caracterizada por tener una superficie plana, ligeramente rugosa y ondulada (^{López-Ramos, 1972}).

Por sus características físicas y litológicas, las rocas presentes hacia la carretera de Ucum-La Unión y González Ortega (Figura 2) forman parte de un acuífero de calizas de características variadas y depósitos aluviales, cuya permeabilidad primaria es mayor en las porciones que contienen fracturamiento y carsticidad, y menor en los estratos de caliza masiva (^{CAPA, 2012a}; ^{Perry, Paytan, Pedersen, & Velázquez-Oliman, 2009}).

La mayor parte de la región carece de un sistema de drenaje superficial y sólo en algunas áreas se ha formado una red hidrológica desintegrada no muy bien definida, que desaparece en resumideros o cuerpos de agua superficial. El acuífero se explota por medio de múltiples alumbramientos; los tipos de captación son norias, pozos someros y pozos profundos, que se utilizan para diferentes usos, como público urbano, agrícola y de servicios, entre otros (^{JICA, 2004}).

La recarga del acuífero tiene lugar durante los meses de mayo a octubre (^{UQROO, 2002}). De esta manera, gran parte de la precipitación pluvial se infiltra al subsuelo a través de fracturas, oquedades y conductos cársticos en las calizas y evaporitas (^{CCPY, 2012}). El agua fluye por el subsuelo, alcanzando las costas por donde es descargada al mar; el resto se pierde por evapotranspiración. El espesor saturado de agua dulce decrece desde el interior hacia la costa (^{CNA, 2006}, ²⁰⁰⁹).

En el área de estudio se identifican dos regiones de acuíferos: la primera corresponde a un acuífero compuesto por rocas carbonatadas del periodo del Eoceno-Paleoceno, ubicado hacia las porciones con mayor altitud del sur-suroccidente; la segunda está compuesta por rocas carbonatadas del periodo del Mioceno-Plioceno, distribuida a largo de la costa (^{JICA, 2004}). De acuerdo con ^{Marín, Perry, Essaid y Steinich (2004)}, ambas regiones se pueden describir como acuíferos de porosidad triple, donde el flujo del agua subterránea circula a través de la matriz, fracturas y un sistema de cavernas interconectadas con orientación NNO-SSE, casi perpendicular a la costa.

En el subsuelo, el agua sigue diferentes trayectorias de flujo controladas por el desarrollo o evolución del carst profundo (^{Velázquez, 1986}). La dirección del flujo de agua subterránea se presenta de SO-NE en las zonas de mayor altitud y de O-E hacia la costa (^{Marín, Perry, Essaid, & Steinich, 2001}; ^{Gondwe, Hong, Wdowinski, & Bauer-Gottwein, 2010}; ^{Charvet, 2009}; ^{JICA, 2004}; ^{CNA, 2001}, ²⁰⁰⁷, ²⁰⁰⁹; ^{Perry, Velazquez-Oliman, & Marin, 2002}). La tendencia es en relación con la mayor densidad de fracturamiento de las formaciones geológicas que lo conforman (Figura 3), moviéndose de las zonas de mayor precipitación hacia la costa, constituyendo hacia esta zona su área de descarga, alimentando a su vez al sistema Bacalar-Río Hondo, la bahía de Chetumal y mar Caribe (^{Velázquez, 1986}; ^{CNA, 2007}).

Figura 3 Mapa de flujo de agua subterránea a escala regional modificado de ^{Sánchez-Sánchez, Álvarez-Legorreta, Pacheco-Ávila, González-Herrera y Carrillo-Bibriezca (2015)}. Las flechas indican las direcciones de flujo de agua subterránea según la referencia respectiva. El contorno en negro representa el área de estudio. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 

Metodología

Los datos fisicoquímicos de la calidad de las aguas subterráneas corresponden a dos periodos de muestreo (años 2002 y 2012) realizados durante la temporada de lluvias. La información del periodo 2002 forma parte del estudio de Aguas Subterráneas Serie II, realizado por el Instituto Nacional de Estadística y Geografía (^{INEGI, 2002}), del cual se obtuvieron 103 datos fisicoquímicos cubriendo la región de estudio. La información del periodo 2012 fue obtenida a partir un muestreo realizado por ^{CAPA (2012b)}, cubriendo un total de 120 datos fisicoquímicos.

Según los informes de ^{INEGI (2002)} y ^{CAPA (2012b)}, todas las muestras fueron colectadas de pozos de abastecimiento uniformemente distribuidos en el área de estudio. Todos los sitios contaban con su ubicación geográfica (coordenadas X, Y), obtenidas a partir del Sistema de Posicionamiento Global (GPS, por sus siglas en inglés) y representadas en el Sistema de Coordenadas Métricas UTM (del inglés, Universal Transverse Mercator). Además, estos informes mencionan que todas las mediciones de temperatura (°C), pH y conductividad eléctrica (CE) se verificaron en campo, usando un Hydrolab Quanta G, y los análisis químicos se determinaron en laboratorios certificados empleando las técnicas analíticas que establecen las Normas Oficiales Mexicanas (Cuadro 1). Las muestras analizadas en el laboratorio para sus constituyentes químicos se determinaron de la siguiente manera: alcalinidad total y alcalinidad a la fenolftaleína por titulación con ácido sulfúrico 0.02N (NMX-AA-36-2001), dureza total y dureza de calcio por titulación con EDTA (NMX-AA-072-SCFI-2001), cloruros por titulación con nitrato de plata (NMX-AA-073-SCFI-2001), sólidos disueltos totales por el método gravimétrico (NMX-AA-034-SCFI-2001), nitratos por el método de reducción con cadmio cuperizado (NOM-AA-079-SCFI-2001), sulfatos por el método turbidimétrico (NMX-AA-074-SCFI-1981), sodio por el método potenciométrico (ion selectivo).

Cuadro 1 Técnicas analíticas según las Normas Oficiales Mexicanas. 

Índices de calidad del agua (ICA)

El índice de calidad del agua fue desarrollado por ^{Brown et al. (1970)} y mejorado por Deininger (^{Scottish Development Department, 1975}). El ICA se define como una técnica de clasificación que engloba la influencia compuesta de cierto grupo de parámetros de calidad del agua en un solo valor general (^{Aly et al., 2014}). El método ICA ha sido utilizado debido a su capacidad de expresar de forma integral la información de la calidad del agua subterránea de alguna región determinada (^{Sadat-Noori, Ebrahimi, & Liaghat, 2014}).

El ICA fue calculado con el fin evaluar la calidad del agua subterránea de la región de estudio de acuerdo con los límites máximos permisibles establecidos en la Norma Oficial Mexicana (NOM-127-SSA1-1994). Se utilizaron siete parámetros químicos de calidad de agua (pH, SDT, dureza total, Na+, Cl^-, SO₄^2- y NO₃^-) y el cálculo consistió de tres etapas. En la primera, a cada uno de los siete parámetros se les asignó una ponderación (w_i) acorde con su importancia relativa en la calidad del agua para consumo humano, así como de las condiciones geohidrológicas de la región de estudio (zona cárstica). La ponderación máxima fue de cinco y se asignó a sólidos totales disueltos y nitratos, en función de los efectos adversos en la salud humana (^{Camargo & Alonso, 2007}; ^{Pacheco, Cabrera, & Pérez, 2004}). A los iones sodio y cloruro se les asignó una ponderación de cuatro, ya que son indicativos de mezcla con agua marina (intrusión salina) en zonas costeras. La ponderación de tres fue asignado a la dureza total y al ion sulfato; al pH se le asignó la ponderación mínima de uno debido a su variación mínima.

En la segunda etapa se determinó el peso relativo (W_i ) de cada parámetro de acuerdo con la ecuación del método del peso aritmético (^{Brown et al., 1970}; ^{Horton, 1965}):

1

Donde W_i es la ponderación relativa, w. es el peso de cada parámetro y n es el número de parámetros.

En la tercera etapa se estableció una escala de valoración de la calidad para cada parámetro (qi). Para esto, se dividió la concentración resultante de cada muestra de agua analizada por su respectivo límite máximo permisible de acuerdo con NOM-127-SSA1-1994 y el resultado fue multiplicado por 100:

2

Donde q_i. es la calificación de la calidad, C_i es la concentración de cada parámetro químico en cada muestra de agua en mg/l y S _j es el valor estándar de agua potable en mg/l (Cuadro 2).

Cuadro 2 Pesos relativos de los parámetros considerados para el cálculo del índice de calidad de agua subterránea según su importancia para la salud humana. 

Finalmente, para calcular el ICA se determinó el Sli usando la ecuación (3):

3

Donde Sli es el subíndice de cada parámetro y q_i es la calificación de la calidad:

4

Donde ICA es el índice de calidad del agua y Sli es el subíndice de cada parámetro.

El resultado final son valores clasificados en cinco categorías de calidad del agua: excelente (< 50); buena (de 50 a < 100); pobre (de 100 a < 200); muy pobre (200-300), y no apta para consumo humano (> 300).

Se realizó un análisis estadístico descriptivo para comparar los valores de las medias, desviaciones estándar y varianza de cada parámetro (Cuadro 3). Para determinar si existe diferencia en la calidad del agua entre las muestras de agua correspondientes a los años 2002 y 2012, se hizo la prueba de Kolmogorov-Smirnov debido a que los valores no se ajustaron a una distribución normal. Dichos análisis estadísticos se efectuaron a través del programa estadístico Statgraphics Centurión, v. XV.

Cuadro 3 Estadística descriptiva del agua subterránea durante la temporada de lluvias de 2002. 

Uso de los sistemas de información geográfica (SIG)

En este estudio, la modelación de los patrones de distribución espacial de cada parámetro de calidad del agua se desarrolló con base en técnicas geoestadísticas, usando la interpolación Kriging. A través de estas técnicas se pueden estimar los valores espaciales en aquellos sitios no muestreados o en el que no se ha hecho algún tipo de medición. El método Kriging se sustenta sobre modelos estadísticos que incluyen la autocorrelación y generación de superficies de predicción con elevado grado de certeza (^{Goovaerts, 1997}). Kriging comprende tres etapas: análisis estadístico-exploratorio de los datos, modelo del variograma (ajuste teórico) y creación de la superficie de predicción (^{García, 2004}). Su ecuación general es la siguiente:

5

Donde Z^(0) es la predicción de la variable de interés, Z, en el punto X₀; Z(X_i) es el valor observado en el punto (base) X. con i = 1, 2,.. N, y λ_i es un coeficiente de ponderación con i = 1, 2,...N. Estas ponderaciones se estiman considerando la correlación espacial y configuración de observaciones a través de los modelos del variograma equiparados a los variogramas experimentales (^{Castiglioni, Castellarin, & Montanari, 2009}).

Se elaboraron los mapas de distribución espacial de cinco variables hidrogeoquímicas (dureza total, Na+, Cl^-, SO₄^2-, NO₃^-). Los mapas se crearon utilizando el programa ArcGIS 10.3. El ajuste teórico del variograma se realizó a través del programa GS plus, versión 9.

Resultados y discusiones

Análisis estadístico descriptivo

El resumen estadístico de la composición química de las muestras del agua subterránea para los años 2002 y 2012 se tiene en los Cuadros 3 y 4. Las muestras de agua subterránea presentaron valores de pH de 6.71 a 8.23, indicando que el agua fue moderadamente alcalina para los periodos estudiados. No obstante, ninguna de las muestras analizadas estuvo fuera de los límites permisibles de la NOM-127-SSA1-1994.

Cuadro 4 Estadística descriptiva del agua subterránea durante la temporada de lluvias de 2012 

Los STD para el periodo 2002 variaron entre 228.2 y 2 877.0 mg/l, con una media y desviación estándar de 847.15 ± 448.51 mg/l. Para el periodo 2012, el contenido de STD varió entre 288.6 y 7 466.85 mg/l, con una media de 1 182.04 ± 880.26 mg/l (Figura 4a). Además de exhibir una gran dispersión, se observó para el periodo 2002, que más de 22% de las muestras rebasó el límite máximo permisible establecido en la NOM-127-SSA1-1994 (1 000 mg/l), mientras que para 2012, más de 42% de las muestras excedió dicho límite (Figura 4b). En general, las menores concentraciones de SDT se localizaron en la planicie interna del área de estudio, con valores entre 500 y 1 000 mg/l, incrementándose gradualmente hacia las costas hasta alcanzar concentraciones mayores a 2 800 mg/l para 2002 y mayores a 7 400 mg/l para 2012 (Figura 4b). La magnitud y diferencia de los valores encontrados demuestran el efecto de salinización que da origen a aguas de mala calidad en algunas regiones del área de estudio, en especial en la zona norte. Esto coincide con lo reportado por ^{CCPY (2012)}, en el que reportan que la gran dispersión del contenido de STD entre la porción norte y sur del área de estudio sugiere que el emplazamiento de las aguas subterráneas se da en materiales litológicamente diferentes, unos más solubles que otros; sin embargo, hacia la porción noreste, la menor variación indica la presencia de materiales muy solubles.

Figura 4 Concentraciones de sólidos totales disueltos (a) 2002 y (b) 2012) en el agua subterránea de la zona de estudio. Los puntos representan los sitios monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q) . 

La temperatura (T °C) no presentó mayores cambios en ambos años, teniendo valores medios alrededor de 27 °C, con una desviación estándar de ±1 aproximadamente. El orden de abundancia de los aniones (SO₄^2- > Cl^- > NO₃^-) fue similar para ambos años. La amplitud encontrada en los valores máximos y mínimos de los cationes y aniones revela heterogeneidad espacial del medio acuífero debido al cambio dinámico entre los diferentes procesos hidro-geoquímicos presentes en el área de estudio, como la disolución de minerales de rocas carbonatadas, dolomías y evaporitas, y la mezcla de agua dulce-marina.

Análisis geoestadístico

La prueba Kolmogorov-Smirnof (K-S) indicó que los valores de dureza total, Na⁺, Cl^- y SO₄^2-no presentaron diferencias estadísticamente significativas para ambos periodos, por lo que se agruparon para formar un solo conjunto de datos. Por el contrario, los valores de NO₃^- exhibieron diferencias estadísticamente significativas (p > 0.05). Se obtuvo la media aritmética para aquellos valores que coincidieron en un mismo sitio en ambos periodos.

A partir de los resultados de la prueba de K-S se construyeron los variogramas omnidireccio-nales para cada variable química con ajuste a un modelo teórico isotrópico (Cuadro 5). En ningún caso se consideró la existencia de anisotropía debido al número de datos disponibles.

Cuadro 5 Modelo teórico de variograma con mejor ajuste para los parámetros hidroquímicos. 

Los modelos teóricos a los que se ajustaron los semivariogramas experimentales fueron el modelo esférico, modelo exponencial y modelo gaussiano. De éstos, el parámetro con mejor coeficiente de determinación (R²⁾ fue el ion cloruro con 89%, mientras que el ion NO₃^- del periodo 2002 fue el que menor coeficiente de determinación mostró con 67%. Se utilizó el método de interpolación Kriging ordinario, para la elaboración de los mapas de distribución espacial para cada variable química.

Mapas de distribución nitratos

La presencia de nitratos en el agua subterránea puede tener diversos orígenes relacionados con actividades humanas, como la agricultura; disposición de los desechos de humanos y de animales; infiltración de aguas residuales sin tratamiento, y en pequeña proporción la precipitación (^{Pacheco et al., 2004}; ^{Price, 2003}). En contraste con el ion amonio, el nitrato puede moverse con libertad a través del acuífero, por lo que difícilmente precipita (^{Hounslow, 1995}). Diversos estudios acerca de la contaminación por nitratos en aguas subterráneas han mostrado que por lo general existe una relación inversa entre las concentraciones de nitratos y profundidad por debajo del nivel freático (^{Freeze & Cherry, 1979}; ^{Ritter & Chirnside, 1984}). ^{Pacheco (1985)} reportó para una región rural al norte del estado de Yucatán que la tendencia de las concentraciones de este contaminante en el agua subterránea es de naturaleza estacional, con las mayores concentraciones asociadas con los centros poblacionales.

Las Figuras 5a y 5b muestran que existió una diferencia en la distribución espacial de NO₃^-para los periodos 2002 y 2012, como resultado posible de las actividades agrícolas, por lo que sus cantidades pueden variar de una temporada a otra o de un año a otro (^{Velasco, 2012}). Todos los valores de NO₃^- del periodo 2002 y 2012 se encontraron por debajo del límite máximo permisible (10 mg/l) establecido en la NOM-127-SSA1-1994.

Figura 5 Concentraciones de nitratos (a) 2002 y b) 2012) en el agua subterránea de la zona de estudio. Los puntos representan los sitios monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 

Es importante señalar que cuando el nitrato está en concentraciones de 10 mg/l o superiores, el agua es indeseable para consumo humano debido al posible efecto tóxico que puede tener sobre los niños (cianosis) (^{Hernández, 2008}; ^{Mora & Alfaro, 1999}; ^{Pacheco & Cabrera, 2003}).

Dureza total

La dureza total representa la suma de las durezas individuales debidas a los iones de calcio, magnesio, y es expresado como mg/l de CaCO₃ (^{Appelo & Postma, 1996}; ^{Manahan, 1991}). Las aguas duras son por lo general incrustantes, mientras que las blandas suelen ser agresivas (^{Lillo, 2007}). Dependiendo de varios aspectos, como pH del agua y concentración de carbo-natos de calcio y magnesio, pueden producir procesos de corrosión o incrustación en las tuberías utilizadas para el abastecimiento de aguas domésticas (^{Mora & Alfaro, 1999}).

Los valores de dureza total de las aguas subterráneas en el área de estudio variaron entre 122.5 y 2 304.84 mg/l como CaCO₃. Los valores obtenidos se agruparon en tres categorías: las dos primeras cumplen con lo establecido en la NOM-127-SSA1-1994 (< 500 mg/l) y la tercera presentó valores que excedieron dicho límite (Figura 6).

Figura 6 Concentración de dureza total en el agua subterránea de la zona de estudio. Los puntos representan los sitios monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 

En el área de estudio, 60% de las muestras rebasó la norma, siendo el valor máximo encontrado de 2 304.84 mg/l como CaCO₃. Con base en la clasificación de ^{Custodio y Llamas (1976)}, el agua subterránea del área de estudio se clasificó como de dura a muy dura, lo que pudiera explicarse por la presencia de minerales disueltos, como carbonatos (CaCO₃), dolomita (CaMg (CO₃)₂) y anhidrita (CaSO₄), característicos del sistema cárstico de la península de Yucatán (^{CCPY, 2012}; ^{Pacheco et al., 2004}).

Es importante señalar que el problema que afronta la población al utilizar aguas muy duras no sólo radica en una mala cocción de los alimentos, o en la generación de sarro e incrustaciones en los utensilios de cocina, y en la poca disolución del jabón, sino además predispone a los pobladores a padecer ciertas enfermedades. Por ejemplo, en el estado de Yucatán, ^{Medina, Zaidi, Real y Orozco (2002)} reportaron que el antecedente familiar y la dureza extrema (> 400 mg/l de CaCO₃) del agua empleada para consumo humano son factores que predisponen el desarrollo de la enfermedad de urolitiasis, la cual se incrementa de forma notable con la edad.

Sodio

El sodio no se encuentra libre en la naturaleza, la fuente común de los niveles de sodio en el agua subterránea se debe, en algunos casos, a la presencia de rocas evaporitas, al efecto de la intrusión marina, o bien a la inadecuada disposición de desechos de humanos y animales (^{Basavaraddi, Kouser, & Puttaiah, 2012}; ^{Lillo, 2007}; ^{Pacheco et al., 2004}). Los valores de las concentraciones del sodio de las aguas subterráneas en el área de estudio oscilaron entre 1.0 y 1 060 mg/l. Con base en el límite máximo permisible de 200 mg/l, establecido en la NOM-127-SSA1-1994, se estimó que 9.8% de los aprovechamientos rebasó la norma (Figura 7).

Figura 7 Concentración de sodio en el agua subterránea de la zona de estudio. Los puntos representan los sitios monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 

Con el fin de identificar los procesos geoquímicos que originan la presencia del sodio en el agua subterránea del área de estudio, se elaboró un diagrama con la relación Na+/Cl^- en mEq/l (Figura 8). De manera general, se observó que las aguas subterráneas del área de estudio se encontraban ligeramente enriquecidas en sodio. Además, la mayoría de las muestras estaba por debajo de 10 y 15 mEq/l para Cl^- y Na+, respectivamente.

Figura 8 Relación del Na+ y Cl-. Los puntos representan las muestras. 

Algunos estudios de toxicidad del sodio en las aguas subterráneas y su efecto adverso en la salud humana revelan que la ingestión de sodio en altas concentraciones puede causar hipertensión, falla del corazón, enfermedades renales, cirrosis hepática, toxemias del embarazo, enfermedad de Menier (afección del oído interno que causa vértigo), e hipernatremia (elevado nivel de sodio en la sangre) en los niños (^{Catalán, 1981}). Asimismo, como cloruro de sodio puede causar náuseas, vómitos, inflamación del tracto gastrointestinal, conmutación muscular, convulsiones e incluso la muerte (^{Basavaraddi et al., 2012}).

Cloruros

La mayor parte del cloruro que se encuentra en las aguas subterráneas suele proceder de cuatro fuentes naturales: 1) agua marina de sedimentación atrapada en el interior de sedimentos; 2) disolución de la sal común y minerales afines que existen en formaciones de evaporitas; 3) concentración por evaporación de los cloruros existentes en agua de lluvia, y 4) disolución de las partículas de material sólido existentes en la atmósfera (^{Davis & De-Wiest, 1971}). De forma antropogénica, el aumento de cloruros puede ser motivado por ciertos desechos industriales, aguas residuales domésticas, retornos de aguas de riego o bien residuos de actividades en zonas pecuarias (^{Catalán, 1981}; ^{CCPY, 2012}). El ion cloruro no forma sales de baja solubilidad, no se oxida ni se reduce en aguas naturales, tampoco es adsorbido de modo significativo ni entra a formar parte de procesos bioquímicos, por lo que se puede considerar como una propiedad conservativa del agua, resultando en un buen trazador (^{Manahan, 1991}).

La distribución espacial de los cloruros en el agua subterránea del área de estudio (Figura 9) mostró un rango de concentraciones de 14.0 a 1 578 mg/l. El 9.9% de las muestras excedió los 250 mg/l, límite establecido en la NOM-127-SSA1-1994. Estos valores se ubican geográficamente hacia las porciones noreste y noroeste de la región de estudio. Esta distribución puede estar asociada con la influencia de aguas marinas, aunque no se descartan la relación con aguas derivadas de retornos de riego o aguas residuales derivadas de asentamientos humanos.

Figura 9 Concentración de cloruros en el agua subterránea de la zona de estudio. Los puntos representan los sitios monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 

Sulfatos

El origen de los sulfatos está asociado principalmente con la disolución de sales sulfatadas de yeso y anhidrita (^{Hounslow, 1995}; ^{Lillo, 2007}). El sulfato está sujeto a procesos de reducción, en especial en presencia de bacterias y materia orgánica; además, posee la capacidad de asociarse con facilidad con otros iones (^{Hem, 1986}). En ambientes reductores (pH < 7), la forma reducida estable es el ácido sulfhídrico (H₂S), mientras que en soluciones alcalinas predomina el ácido sulfúrico (HS) (^{Lillo, 2007}).

La concentración del ion sulfato en las aguas de bebida ha sido motivo de atención debido a su acción catártica, habiéndose encontrado un efecto laxante de las aguas que contienen más de 750 mg/l (^{Cabrera 1986}; ^{Catalán, 1981}). En la península de Yucatán, los sulfatos son uno de los principales iones que se encuentran en las aguas subterráneas y además difícilmente alcanzan la saturación debido a que su producto de solubilidad es muy alto (^{CCPY, 2012}).

En el área de estudio, los valores del ion sulfato se encontraron en el rango de 0.01 a 2 000 mg/l. Dentro de este rango, 27.8% de los aprovechamientos rebasó el límite de 400 mg/l que establece la NOM-127-SSA1-1994. La Figura 10 muestra que la zona suroeste y la franja central del área de estudio son las que presentaron menores concentraciones de sulfatos; sin embargo, existen áreas dentro de estas zonas con concentraciones superiores a la establecida en la NOM. No obstante, estas últimas presentan de forma sorprendente concentraciones bajas en cloruros; ello sugiere que sus características provienen en específico de la disolución de evaporitas (yesos/ anhidritas) y no de la mezcla con agua de mar. Por otra parte, se observó un incremento gradual hacia las zonas noreste y noroeste, donde se llegaron a alcanzar concentraciones de hasta 2 000 mg/l de sulfatos. Este incremento puede estar asociado con tres factores: 1) la presencia de rocas evaporitas (yesos y anhidritas), lo que genera un incremento sustancial en la salinidad del agua subterránea (^{Perry, Velazquez-Oliman, & Wagner, 2011}); 2) la descarga de las aguas de las zonas acuíferas de las porciones altas, y 3) la presencia de aguas marinas.

Figura 10 Concentración de sulfatos en el agua subterránea de la zona de estudio. Los puntos representan los sitios monitoreados. 

Índice de calidad del agua subterránea

Para los dos periodos analizados, se encontró una gran variación en los índices de calidad del agua. Respecto a 2002, el índice más bajo fue de 14 y el máximo de 174 (Figura 11a). El valor medio del ICA es de 53 y la mediana de 47, lo que indica que para este periodo aproximadamente la mitad de los pozos se encontró en la categoría de "excelente". De los aprovechamientos muestreados, 58% se ubicó en la categoría de "excelente"; 36%, en la categoría "buena", y el 6% en la categoría "media". En cuanto a 2012, el índice más bajo encontrado fue de 18 y el máximo de 372 (Figura 11b). El valor medio del ICA es de 74 y la mediana de 69. Estos valores permitieron constatar una degradación en la calidad del agua, pues la mitad de los pozos se clasificó en las categorías "excelente" y "buena". De los aprovechamientos muestreados, 34% se ubicó en la categoría de "excelente"; 48%, en la categoría "buena"; 15%, en la categoría "media"; 2% en la categoría "pobre", y 1% en la categoría "extremadamente pobre". Estas dos últimas categorías se ubicaron geográficamente en la porción central este del área de estudio. Tal distribución puede estar asociada con el crecimiento poblacional de las dos comunidades que se localizan en esa región, ya que pasó de 870 habitantes en el año 2000 a 1 500 en 2010 (^{INEGI, 2014}). Este crecimiento dio origen a la generación de más aguas residuales domésticas. Otra causa de dicha distribución es el uso de agroquímicos en actividades agrícolas (producción de piña) que se desarrollan en estas entidades. Respecto a los procesos hidrogeoquímicos que pudieran complementar el elevado ICA encontrado, los valores concuerdan con lo reportado por ^{CCPY (2012)}, en el que se identificaron elevadas concentraciones de sólidos disueltos totales, sulfatos y cloruros en el pozo de abastecimiento de una localidad de esa región, con la posible influencia de aguas marinas, y la asociación con aguas derivadas de retornos de riego y aguas residuales.

Figura 11 Índice de calidad del agua (a) 2002 y (b) 2012 de la zona de estudio. Los puntos representan los sitios monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 

De manera general, el mapa de distribución del ICA mostró un patrón espacial cualitativamente similar para ambos periodos. La calidad del agua subterránea disminuye desde el suroeste hacia el norte, exceptuando la parte sur suroeste, así como algunas zonas de la parte noreste del área de estudio, en donde se encontró la zona más segura en términos de calidad de agua para consumo humano. En estas regiones, casi todos los valores del ICA se encontraron en la categoría de "excelente". La disminución de la calidad se asocia de manera directa con el aumento gradual en las concentraciones de sulfatos, dureza total, cloruros y sodio.

Conclusiones

Los constituyentes químicos que mostraron concentraciones que excedieron el límite permisible de la NOM-127-SSA1-1994 fueron sólidos totales disueltos (22% en 2002 y 42% en 2012), sodio (9.8%), cloruros (9.9%) y dureza total (60%).

El uso de los sistemas de información geográfica en la elaboración de los mapas de distribución espacial de las concentraciones permitió observar que las menores concentraciones para los sólidos totales disueltos, cloruros, dureza total y sodio se localizaron en la planicie interna sur del área de estudio, incrementándose sus concentraciones gradualmente hacia la zona noroeste. El mapa de las concentraciones de nitratos mostró una diferencia en la distribución espacial para los periodos 2002 y 2012.

El ICA permitió demostrar que la calidad química del agua subterránea con fines de abastecimiento es aceptable para la mayoría de los sitios estudiados. No obstante, el análisis comparativo de los periodos 2002 y 2012 indica de manera general una degradación en la calidad del agua subterránea con los máximos valores del ICA ubicados hacia la parte noreste del área de estudio, que corresponde con zonas que presentaron un crecimiento de la población y de zonas agrícolas.

Este estudio demuestra con claridad que el uso de métodos como el valor ICA podría ser de gran utilidad para evaluar la calidad de las aguas subterráneas en sistemas cársticos, proporcionando valiosas herramientas a los tomadores de decisiones, a fin de comprender de forma integral el estado de la calidad del agua subterránea y tener la oportunidad de tomar decisiones adecuadas para un mejor uso en el futuro.

Por último, este estudio demuestra la necesidad de identificar zonas alternas para abastecer agua en aquellas poblaciones donde se reportan los valores más altos del ICA, pues sin duda se están suministrando aguas que sobrepasan los límites permisibles de la NOM-127-SSA1-1994.