Aplicación de los clorofluorocarbonos y el hexafluoruro de azufre como trazadores hidrogeológicos en los acuíferos carbonatados del Parque Natural del Barranco del Río Dulce (Guadalajara)
Using chlorofluorocarbons and sulphur hexafluoride as hydrogeological tracers in the limestone aquifers of Parque Natural del Barranco del Río Dulce (Guadalajara)

M. Martín-Loeches Garrido¹, R. Fisher², S. Díaz Alcaide³, P. Martínez Santos³ 

¹ Departamento de Geografía y Geología de la Universidad de Alcalá. Unidad docente de Geología. Edificio de Ciencias del Campus Universitario. Ctra. N-II, km. 33,600. 28871 Alcalá de Henares, Madrid, España. Email: miguel.martin@uah.es

² NΒΊ 402, Village Trea, Stung MeanChey, MeanChey, Phnom Penh, Camboya

³ Departamento de Geodinámica, Facultad de Ciencias Geológicas, Universidad Complutense de Madrid 28040 Ciudad Universitaria Madrid, España

IntroducciónTOP

Los clorofluorocarbonos (CFCs) son compuestos sintéticos de cloro, carbono y flúor que fueron producidos comercialmente desde la década de los años 30 del pasado siglo para refrigeración y otras aplicaciones industriales. El primero en liberarse a la atmósfera fue el CFC-12, en 1930; desde 1936 se produjo el CFC-11, y finalmente, desde 1944, el CFC-113 (IAEA, 2006). La concentración de CFCs en la atmósfera se incrementó de manera muy notable desde los años 50 hasta el comienzo de los años 90, cuando dejaron de emitirse tras los acuerdos derivados del Protocolo de Montreal. Actualmente, las curvas de concentración de estos gases en la atmósfera se encuentran decreciendo.

Los CFCs se han utilizado con éxito en numerosos estudios como trazadores y para datar aguas subterráneas infiltradas después de 1950 (Darling et al., 2012; Cook & Solomon, 1995 y 1997; Plummer & Busemberg, 1999; Busemberg & Plummer, 1992 y 2000; IAEA, 2006); también para estimar las proporciones de mezcla de agua de diferentes edades en las muestras (Plummer et al., 2003; IAEA, 2006). Indirectamente, estos compuestos se pueden usar para calibrar modelos de flujo y para evaluar la susceptibilidad de una fuente de agua potable a contaminarse (Busenberg & Plummer, 1992).

El hexafluoruro de azufre (SF₆) comenzó a emitirse a la atmósfera a partir de 1953 y es un gas con un claro potencial para los mismos propósitos (Busenberg & Plummer, 2000, 2006 y 2008). A diferencia de lo que ocurre con los CFCs, la concentración en la atmósfera del SF₆ sigue incrementándose, lo que supone un factor importante para el futuro de su uso como trazador.

Ninguno de estos compuestos ha sido utilizado antes en España en el ámbito hidrogeológico.

La datación con CFCs y SF₆ se basa en la comparación de las concentraciones de estos elementos en equilibrio en las muestras de agua subterránea, con el contenido histórico de sus concentraciones en la atmósfera. La Fig. 1 muestra los valores de CFC-11, CFC-12, CFC-113 y SF₆ que se esperarían en las aguas subterráneas del hemisferio norte recargadas entre 1950 y 2010 a 10 °C de temperatura y a nivel del mar. Según los datos de esta figura, en condiciones ideales, el agua recargada en 1985 y correspondiente a un modelo de flujo pistón (PFM) (muestra de agua sin mezcla, recogida de un sondeo equipado con una rejilla estrecha en un acuífero detrítico homogéneo e isótropo, o de un manantial somero), debería contener, 4.5 pmol/l de CFC-11, 2 pmol/l de CFC-12 y 0.3 pmol/l de CFC-113 y 0.65 fmol/L de SF₆ (ver Fig. 1).

El modelo de flujo pistón (PFM) es aquel en el que no existe mezcla de aguas de diferente procedencia o edad; se suele describir como la situación ideal a la que se llegaría si el agua fluyera a través de un tubo de flujo desde la recarga hasta la descarga sin sufrir ninguna interacción con tubos de flujo contiguos. En estos casos, efectivamente, el agua presentaría una edad uniforme.

Es difícil sin embargo recoger muestras que se ajusten a esas condiciones. El agua que se bombea de un pozo o que emerge en un manantial es normalmente una mezcla de aguas procedentes de diferentes líneas de flujo. Para la interpretación de los contenidos en CFCs y SF₆ en estos casos, debemos recurrir a otros modelos; los más utilizados son el modelo de mezcla exponencial (EXM) y el modelo de mezcla binaria (MB) (IAEA, 2006); existen además modelos intermedios, como el flujo pistón exponencial o el de mezcla exponencial (Darling et al., 2012).

En el modelo exponencial, la muestra se interpreta como un número infinito de sub-muestras que presentan una distribución exponencial de edad y se aporta una edad media general. Esta situación puede corresponder con la que sufrirían muestras tomadas de pozos abiertos en una longitud amplia de su tubería en acuíferos detríticos, o en manantiales de descarga regional.

En el modelo de mezcla binaria se considera al agua de la muestra como el resultado de la mezcla entre “agua vieja”, o libre de CFCs (recargada anteriormente a 1950) y “agua joven” o con contenido en CFCs (post 1950). Este modelo es uno de los más simples y quizá el más importante en acuíferos fracturados (Gooddy et al., 2006). Las muestras tomadas en acuíferos carbonatados en los que existe una red de canales de flujo preferente junto con una matriz rocosa saturada donde se produce un flujo mucho más lento, podría ajustarse a este modelo.

Existen además algunos factores que deben tenerse en cuenta ya que pueden distorsionar los valores reales de edad. Entre ellos destacan la temperatura y la altitud en la recarga, el espesor de la zona no saturada, la elevada concentración de estos compuestos en el aire urbano, la degradación microbiana – que afecta de manera diferente a cada compuesto-, la absorción de estos compuestos a la fase sólida del acuífero, si este tiene alto contenido en carbono orgánico, y las mezclas de aguas (Darling et al., 2012).

En el presente trabajo se han analizado los contenidos de CFCs y SF₆ de aguas subterráneas procedentes de los acuíferos carbonatados de edad Jurásico y Cretácico que se encuentran en el entorno del Parque Natural del Barranco del Río Dulce (PNBRD) en la provincia de Guadalajara (Fig. 2) y que pertenecen a la masa de agua subterránea “030.003 Tajuña-Montes Universales” (CHT, 2013) o la Unidad Hidrogeológica homónima 03.02 de IGME (2013), para caracterizar su edad (tiempo transcurrido desde su infiltración hasta su muestreo) como primer objetivo a partir de los diferentes modelos mencionados. Siendo la primera vez que se utilizan en España estos compuestos en el campo de la hidrogeología, otro objetivo del trabajo es valorar su utilidad en acuíferos del centro peninsular.

Material y métodosTOP

Las muestras fueron recogidas siguiendo a Busenberg et al. (2006). Se tomaron dos muestras por cada punto para confirmar los resultados. Se utilizaron botellas de vidrio de 500 ml (suficiente para realizar todas las determinaciones) que fueron llenadas con el agua problema bombeada a través de un tubo de nylon de 10 mm de diámetro. Dicha botella debe transportarse al laboratorio en el interior de una cápsula metálica llena también con el agua muestreada.

Durante el muestreo, tanto la botella como la cápsula deben situarse en el fondo de un cubo; se bombea el agua al interior de la botella a través de la manguera de nylon y se deja que fluya hasta que rebose por encima del cubo. Después, en condiciones sumergidas, se cierra la botella y se introduce ésta en la cápsula metálica que también debe sellarse a través de una cinta de goma que encaja a presión. Esta técnica de muestreo crea una barrera contra la contaminación del agua de la botella interior. Antes de proceder a tomar una muestra, el equipo era lavado con el agua a muestrear.

Se debe evitar el contacto del agua con el aire de la atmósfera, ya que de producirse, la muestra se contaminaría. Esto supuso un factor limitante a la hora de seleccionar los puntos de agua ya que muchos no presentaban las condiciones mínimas que garantizasen un completo aislamiento del agua durante el muestreo. Para evitar esta contaminación, dependiendo del tipo de salida del agua, se introducía un extremo del tubo de nylon hasta asegurarse de que no existiera contacto del agua con el aire atmosférico. El otro extremo terminaba en el fondo de la botella de muestras.

Las muestras fueron tomadas en el mes de octubre, tras el periodo de estiaje, para minimizar la influencia de aguas superficiales.

Se tomaron muestras de 6 puntos de agua; su posición se indica en el mapa de la Fig. 2, y en los cortes de la Fig. 3 con las letras identificativas A, B, C, D, E y F. Los resultados se presentan en la Tabla 1. Corresponden a tres muestras en el Jurásico carbonatado y tres en el Cretácico, también carbonatado.

Las determinaciones analíticas de CFCs y de SF₆ se realizaron por cromatografía de gases obtenidos mediante purga y atrapamiento en el laboratorio del Dr. Oster en Wachenheim, Alemania. Las concentraciones de los gases fueron medidas con un detector de captura de electrones (ECD, en sus siglas en inglés). Una exhaustiva explicación de los métodos de análisis utilizados puede obtenerse en Bullister & Weiss (1988) y en Busenberg et al. (2006).

Los datos fueron interpretados con la ayuda de la hoja de cálculo que para tal propósito ofrece el Servicio Geológico de los Estados Unidos (USGS, en sus siglas en inglés) en su página web (USGS, 2013). Esta herramienta calcula las edades aparentes del agua a partir de los registros históricos de concentración en la atmósfera en CFCs y SF₆ según el modelo pistón, el modelo de mezcla binaria y el modelo exponencial. El programa usa los datos de la estación de Niwot Ridge en Estados Unidos que son válidos para todo el hemisferio Norte a excepción de las áreas industriales y urbanas donde existe una mayor concentración de estos compuestos. Los resultados de estas interpretaciones se muestran en la Tabla 2 y la Fig. 5.

Marco geológico e hidrogeológico y modelo conceptual de flujo de los acuíferos estudiadosTOP

El PNBRD se encuentra en el borde centro-oriental de la Cordillera Ibérica, próximo a la Cuenca terciaria del Tajo. La sucesión estratigráfica en esta región está compuesta por materiales triásicos, jurásicos y cretácicos suavemente plegados en la Orogenia Alpina. La Fig. 2 muestra un esquema de localización y un esquema hidrogeológico de la región en el que se señalan las unidades acuíferas que diferencia Matín-Loeches, (2008a) en base a la estructura tectónica, la presencia de la Formación Arenas de Utrillas, del Cretácico, y de las arcillas, limos y yesos en facies Keuper del Triásico, que son materiales de baja permeabilidad.

La erosión ha aislado algunos macizos calcáreos de edad Jurásico, que constituyen unidades hidrogeológicas independientes como la Unidad “Saúca” (Fig. 2). Las unidades Pelegrina-Torremocha (N y S) (Jurásico), Algora-La Cabrera (Cretácico) y Aragosa (Cretácico) (Fig. 2), forman el sistema hidrogeológico Pelegrina-La Cabrera-Aragosa de Martín-Loeches (2008 a).

El sustrato Jurásico puede alcanzar hasta 375 metros de potencia de materiales margosos intercalados con calcáreos (Goy et al., 1976). Las margas son predominantes en los 100 metros superiores, aunque culminan la serie jurásica 10 metros de calizas micríticas (Ríos et al., 2002). Es frecuente que las calizas y dolomías que dominan en las unidades inferiores se presenten en estratos finos y medios o en tramos brechoides, oquerosos (carniolas) y con un elevado grado de fracturación y diaclasado. No existen grandes bancos de calizas o dolomías continuas de edad Jurásico a través de los cuales se desarrolle un sistema kárstico; la única forma kárstica reconocible en superficie, salvo pequeños torcales, es el valle ruiniforme que se localiza en el cañón del río Dulce entre Pelegrina y Jodra del Pinar sobre la formación Calizas y Dolomías Tableadas de Adell et al. (1981a) (García Hidalgo et al, 2008) cerca del nivel de base actual representado por el cauce del río Dulce.

Se han inventariado 11 pozos perforados en los que se registra una oscilación piezométrica anual cercana a los 20 metros entre 2008 y 2013. Solo en 5 de ellos existe acceso y posibilidad de realizar medidas periódicas y de estos, dos son piezómetros de control de la Confederación Hidrográfica del Tajo (CHT, 2013). Los niveles restantes corresponden a los tomados durante la ejecución de las perforaciones. Los pozos más profundos perforados en estas calizas (entre 200 y 256 metros), raramente resultan negativos aunque aportan caudales siempre inferiores a los 5 l/s. Son frecuentes los manantiales estacionales derivados de la alternancia de calizas y margas situación que define en algunas zonas un nivel piezométrico colgado con interés local (Torremocha del Campo).

No hay evidencias de importantes conductos en las calizas Jurásicas a través de los cuales se defina un flujo veloz del agua subterránea, y divida al conjunto rocoso en grandes bloques. Al menos no se han podido localizar puntos importantes de descarga donde estudiar las variaciones de caudal tras las precipitaciones (Pulido, 2001). Las margas en las unidades geológicas superiores pueden jugar un cierto papel en el almacenamiento del agua que es transferida al medio fisurado inferior donde la transmisión es mayor; en la zona saturada las margas deben retardar la velocidad del flujo. Entre los rangos propuestos para modelos de flujo conceptuales en acuíferos carbonatados por diferentes autores (White, 1969, Atkinson & Samart, 1981; Smart & Hobbs, 1986; Quinlan et al., 1992), deducimos que el material de edad Jurásico se situaría entre aquellos en los que domina el flujo difuso con movimiento lento a través del sustrato, comparable al medio de porosidad intergranular con elevada anisotropía y heterogeneidad, que describe para otros macizos carbonatados españoles Pulido, (1998).

La recarga en estos materiales jurásicos es predominantemente difusa, aunque no podemos descartar que una parte se produzca de forma puntual. También pasa agua de las unidades cretácicas a través de las arenas de Utrillas en una cantidad no cuantificada (Martín-Loeches et al., 2011). Su descarga se debe producir en el paraje de los Heros, a 948 msnm, y al pie del Mirador de Pelegrina a 983 msnm (ver Fig. 2) con unos 2 hm³ anuales en años de precipitación media (Martín-Loeches et al., 2011).

Las líneas isopiezas marcadas en las unidades jurásicas en la Fig. 2 constituyen aproximaciones de la posición del nivel de saturación regional, y son válidas para mostrar la dirección preferente del flujo. Su trazado marca un gradiente hidráulico próximo a 0.01 en las zonas más alejadas de los ríos Dulce y Henares, que provocaría un flujo en la matriz del sustrato muy lento.

Los materiales carbonatados del Cretácico lo forman hasta 400 metros de calizas y dolomías en bancos de espesor medio bien definidos, a veces masivos. Culminan en algunas zonas del área de estudio con unos 60–100 metros de yesos, margas, calizas y dolomías en tránsito hacia el Paleógeno (Ríos et al., 2002). La proporción de materiales margosos es mucho menor que en los materiales jurásicos; también se aprecia en las calizas del cretácico una menor presencia de fracturación. Existen formas exokársticas en superficie como torcales y simas y se desarrollan cavidades visitables (cueva de la Moza o de La Oliva, ver Fig. 2).

En las unidades cretácicas al sur del cauce del río Henares no existen pozos. Nos consta que algunos intentos de perforación en el Cretácico fueron infructuosos, circunstancia frecuente en macizos calcáreos típicamente kársticos (Pulido, 1998). La mayoría de los puntos del inventario corresponde a manantiales de descarga regional de entre los que destaca el manantial de La Cabrera porque registra un caudal mínimo de 30 L/s.

La respuesta de los manantiales tras las lluvias es muy rápida y se superpone durante un breve lapso de tiempo a un flujo de base muy continuo y a otro que refleja las lluvias de primavera y otoño. Las estimaciones realizadas en el manantial de La Cabrera (García-Hidalgo et al, 2008) señalan que el flujo de base puede constituir hasta el 50% de la aportación total (0.85 hm³ de 1.4 hm³ totales en años de precipitación media). Ensayos con trazadores realizados en las simas cercanas a Algora muestran en el Cretácico velocidades máximas próximas a los 4000 m/día.

Estos datos señalan que en las unidades cretácicas debe existir un nivel de saturación general, y que coexisten conductos de alta permeabilidad a través de los cuales el agua circula con rapidez y bajo un régimen turbulento, con una red de discontinuidades menores, fisuras, planos de estratificación e incluso poros en la matriz rocosa en la que se produce flujo lento (Worthington, 2007; Ford & Williams, 2007; Goldscheider & Drew, 2007; Pulido, 2001). Los bloques entre conductos constituyen el elemento capacitivo mientras que los canales serían el elemento transmisivo.

La recarga de las unidades cretácicas es en gran medida puntual, a través de las simas; en los núcleos sinclinales aflora la unidad superior impermeable (Fig. 2) que forma áreas de recarga alogénicas (Goldscheider & Drew, 2007). La descarga se produce al rio Dulce en los manantiales próximos a su cauce con 3.15 hm³ en años de precipitación media y en parte al río Henares.

Hemos marcado en la Fig. 2 líneas isopiezas supuestas en las unidades cretácicas que tratan de indicar de manera aproximada la dirección general del flujo subterráneo y su relación con los ríos; no incluye la red de conductos kársticos de mayor velocidad, ya que no se ha podido investigar en detalle su alcance y naturaleza.

Martín-Loeches et al, (2011) aportan valores de permeabilidad en las calizas jurásicas entre 1.5 y 0.014 m/d y en las cretácicas entre 0.5 y 0.045 m/d; aunque estos rangos de valores son normales en acuíferos carbonatados (Domenico & Schwartz, 1998; Fetter 2001; Pulido, 1998 y 2001) resultarían superiores con mayor número de datos en diferentes puntos de las unidades estudiadas. El único valor que se posee en la zona para las Utrillas es de 0.068 m/d y se encuentra en el rango indicado para arenas arcillosas por Fetter (2001).

El agua subterránea en todas las unidades descritas es HCO₃-Ca-Mg con conductividades muy regulares próximas a los 600 μs/cm salvo en los puntos influenciados por las arcillas del Keuper donde se registran hasta 780 μs/cm (Martín-Loeches, 2008b). La presencia de NO₃⁻ es constante, oscilando en los puntos de descarga regional entre 10 y 19 mg/l. No se ha registrado de forma continua ningún parámetro hidroquímico en los manantiales de descarga de los macizos cretácicos, aunque hemos de esperar variaciones notables con la llegada del agua procedente de los conductos kársticos mayores tras las precipitaciones (Pulido, 2001).

La precipitación alcanza en esta región entre 500 y 600 mm en años medios constituyendo los ríos Dulce y Henares el drenaje principal de los acuíferos carbonatados que estudiamos.

A partir de la descomposición de hidrogramas de la estación de aforos 3 254 (Fig. 2), y considerando 87.5 km² de superficie vertiente a la citada estación de aforos, se estima que la infiltración eficaz oscila entre un 14 y un 20% de las precipitaciones de años medios a húmedos respectivamente.

Durante los estiajes rigurosos algunos tramos del río Dulce sobre los materiales de edad Jurásico se encuentran secos (Fig. 2) con la única aportación continua del manantial situado a 983 msnm en el propio cauce (Fig. 4). El agua de los manantiales de Estriégana y Sauca, parece sufrir infiltración en la unidad Pelegrina-Torremocha Norte cerca de la localidad de Jodra del Pinar (Fig. 2) aunque no se ha localizado un sumidero claro; más bien parece producirse una entrada difusa del agua a lo largo de un tramo corto del cauce. A partir de los macizos cretácicos, en el entorno de la localidad de La Cabrera, se mantiene en el río Dulce un caudal mínimo.

Análisis e interpretación de los resultadosTOP

En la Tabla 2 se muestran los resultados de edades aparentes de las aguas subterráneas obtenidos bajo diferentes modelos, tras la aplicación del programa del USGS (USGS, 2013). La Fig. 5 es una representación gráfica de los mismos resultados.

En primer lugar se debe valorar si las concentraciones de los gases medidos corresponden a los aportes atmosféricos originales o por el contrario han sufrido alteraciones que deban tenerse en cuenta en las interpretaciones. En este sentido señalamos los siguientes puntos:

En relación a la interpretación de los datos, en términos generales las muestras con menor concentración en CFCs y SF₆ se pueden considerar las más antiguas (muestras A y E, ver Tabla 1), teniendo en cuenta las curvas históricas de concentración de estos elementos en la atmósfera (Fig. 1), y aunque con ellas se puede llegar a algunas interpretaciones en relación al flujo en los acuíferos, es preferible realizar una interpretación bajo un modelo de mezcla que se ajuste a las condiciones del muestreo y a las características hidrodinámicas del acuífero (Darling et al, 2012).

Una forma de conocer el ajuste de las muestras a los modelos de mezcla posibles es representando los datos en gráficas binarias de concentración de un CFC respecto de otro o respecto del SF₆ (Darling et al., 2012). Teniendo en cuenta la tendencia decreciente de la concentración en la atmósfera de los CFCs, las combinaciones CFC-113 versus CFC-12 y especialmente la CFC-11 versus CFC-12, presentan ambigüedades. El único gráfico que no ofrece problemas es el CFC-12 versus SF₆, ya que la curva que representa el flujo pistón no se cruza con la línea que marca la mezcla binaria. Gooddy et al. (2007) diferencian la existencia de hasta 3 modelos de mezcla a partir de la representación en un gráfico de este tipo de muestras procedentes de una única trayectoria de flujo en un acuífero del Sur de Inglaterra. En nuestro caso, al distorsionarse las concentraciones de SF₆ por aportes del sustrato, se invalida su uso en las interpretaciones. Utilizamos en su lugar la gráfica CFC-113 versus CFC-12, que mostramos en la Fig. 6, que aunque presenta cierta ambigüedad a partir de los últimos años de la década de los 80, puede ofrecernos resultados más completos.

Vemos en la Fig. 6, que todas las muestras se encuentran siempre sobre alguna línea de mezcla, por encima de la línea que marca el flujo pistón.

Descartado el flujo pistón, los resultados pueden interpretarse como mezclas exponenciales o como mezclas binarias de edad. La mezcla exponencial aporta un único valor de tiempo de residencia (Tabla 2) en el que tiene más peso la señal del agua que se encuentra en mayor proporción, mientras que la interpretación según el modelo binario aporta información en relación a una mezcla simple de agua de dos diferentes edades. Así, a la muestra E, Navalpotro, se le asigna un valor de 100 años según el modelo exponencial o se interpreta, según el modelo de mezcla binaria, como una mezcla entre un 30% de agua post-1950 y un 70% de agua pre 1950. Con el modelo de mezcla binaria se cuantifica la proporción de esa mezcla y se aporta un valor de edad de la fracción “joven”. Ni en la mezcla exponencial ni en la binaria las edades aportadas son absolutas (Tabla 2).

Por la naturaleza carbonatada de los acuíferos estudiados, estos deben previsiblemente albergar al menos dos porosidades que derivan en flujos del agua subterránea que transcurren por vías diferentes; una lenta, a través de fracturas de pequeño tamaño y de la matriz de los bloques en los que se encuentra dividido el sustrato (especialmente en las unidades cretácicas) y otra rápida por los conductos kársticos mayores. Consideramos más ajustada por tanto la interpretación de los resultados a través de una mezcla binaria entendiendo que coexisten en las muestras al menos dos fracciones de edad, una de las cuales correspondería a aguas libres de CFCs que llegan al punto de muestreo por la vía lenta y la otra a agua “joven”, que viaja por los conductos kársticos que aunque proceda de múltiples entradas, se homogeniza en agua de infiltración relativamente reciente.

Teniendo en cuenta estas apreciaciones, llegamos a las siguientes interpretaciones sobre las muestras analizadas, que exponemos de forma separada por unidades acuíferas:

Discusión y conclusionesTOP

Al no disponerse de datos de los trazadores medidos en las zonas de descarga de las unidades jurásicas, ni en las zonas de recarga en las cretácicas no es posible llegar a valoraciones completas en relación al flujo del agua subterránea. Los datos sí apoyan el conocimiento previo deducido por otras vías.

A excepción del agua recogida en Jodra del Pinar, las muestras de las unidades cretácicas poseen mayor proporción de aguas jóvenes que las recogidas en las jurásicas, lo que puede indicarnos que los macizos cretácicos poseen mayor permeabilidad y presencia de vías de circulación del agua subterránea más rápidas, lo que está de acuerdo con el modelo conceptual propuesto. Esta suposición no se encuentra refrendada por los datos de conductividad hidráulica medidos en los sondeos, pero hay que tener en cuenta que estos no son representativos del comportamiento general de unos acuíferos de gran heterogeneidad y anisotropía. En este sentido los trazadores constituyen una herramienta discriminante.

Aun considerando que los datos son escasos y se debe seguir investigando, ha quedado claro que las interpretaciones de los trazadores no derivan en valores de edad que se ajusten a la complejidad natural de estos sistemas. Esta conclusión es clara en el caso de La Cabrera ya que es en este punto donde mayor información hidrodinámica se posee y ha podido ser comparada, pero se daría con toda probabilidad en el resto de los puntos de agua.

La elevada edad de las fracciones jóvenes de todos los puntos (siempre por encima de los 17 años) indica que se recogen sin un volumen significativo del agua que descarga de forma rápida tras las precipitaciones importantes y que circula por los conductos kársticos mayores. Esta fracción representa una parte pequeña del almacenamiento en este tipo de acuíferos (menos del 4%, según Worthington, 2000) y se debe estudiar a través del uso de técnicas de trazado con químicos solubles colorantes o con trazadores sólidos en suspensión. A su vez este dato muestra que debe ser importante el flujo que circula a través de la matriz y de las fracturas más pequeñas y medias.

En términos generales los CFCs han sido eficientes para confirmar y cuantificar mezclas de aguas en sistemas binarios, lo que ha supuesto la base de su utilidad en estos acuíferos, suficiente para apoyar un determinado sistema de flujo. Los resultados obtenidos, son coherentes con un modelo de escaso gradiente en la recarga y con la existencia de al menos una doble porosidad en todos los acuíferos. Los CFCs han servido además para poner de manifiesto la posible presencia de contaminación urbana en las muestras.

Además de la ausencia de puntos de agua susceptibles de muestrearse en algunas áreas de estos acuíferos, existen otros elementos que deben considerarse para realizar una mejor interpretación de los contenidos en CFCs y SF₆; caben destacarse un mejor conocimiento de la forma del flujo en la zona no saturada y la verificación de la ausencia de ambiente reductor en algunas zonas profundas de estos macizos aunque activas en el flujo.

Por otro lado, se ha puesto de manifiesto, una vez más, lo importante de acompañar los valores de CFCs y SF₆ con los de otros trazadores (Han, et al., 2006). El ³H y la relación ³H/³He, pueden ayudar a validar las edades aparentes obtenidas con los CFCs; el ¹⁴C puede ser útil para estimar las fracciones de mezcla y la edad de la porción pre-1950 en las muestras, donde se ha establecido la existencia de mezclas; para calcular la temperatura en la recarga es importante además contar con medidas de gases como el N₂ y el Ar (Lindsey, et al., 2003; Yager, et al., 2013).

Existen dos importantes puntualizaciones sobre el futuro del uso de los CFCs y el SF₆ en estos acuíferos carbonatados, ambas relacionadas; en primer lugar, si se confirma que el sustrato carbonatado de estos macizos aporta SF₆ al agua subterránea, lo que puede realizarse con mayores determinaciones, se invalidaría su uso como trazador. En segundo lugar y como consecuencia, cada vez será más difícil distinguir el modelo de mezcla al que puede ajustarse una muestra a partir de su contenido en CFC, elemento muy importante para su interpretación, ya que las curvas de concentración de CFCs en la atmósfera se encuentran decreciendo.

AGRADECIMIENTOSTOP

El muestreo de CFCs y SF₆ se realizó con la ayuda del proyecto de Investigación CGL2004-02179/BTE de la DGYCT y los trabajos posteriores recibieron el apoyo de la Junta de Comunidades de Castilla-La Mancha a través de la subvención GUGE002609. Los autores quisieran expresar su agradecimiento a Juan Antonio Calvo, Javier Temiño Vela y José Francisco García-Hidalgo Pallarés, del Departamento de Geología, Geografía y Medio Ambiente de la Universidad de Alcalá por su apoyo durante los trabajos realizados y a E. Busenberg, del Servicio Geológico de los Estados Unidos y a H. Oster del laboratorio de CFC de Wachenheim (Alemania), por sus comentarios y aportaciones en la interpretación de los datos. Agradecemos a los revisores de este trabajo sus importantes aportaciones.

ReferenciasTOP