Espeleogenesis de los proto-conductos 1
INTRODUCCIÓN
La reacción química que es responsable de la disolución de la roca caliza por agua ácido, ya era conocida por los primeros investigadores de la espeleología. También era conocido que esta reacción solo es capaz de disolver pequeñas cantidades de roca, pero que a su vez ocurre muy rápida. El agua que se infiltra en una estrecha fisura llega a una saturación casi completa, a solo unos pocos metros de profundidad. Esto ponía a los investigadores ante un enigma, porque si el agua de infiltración ya estaba saturado (respecto a la caliza) a unos metros bajo la superficie: ¿cómo es posible que haya suficiente solución a decenas e incluso a cientos de metros de profundidad para formar cuevas? Para que el agua llegara a estas profundidades con bastante capacidad de disolución, se suponían la existencia de unas grietas relativamente anchas. A estas grietas las llamaron proto-conductos. Un proto-conducto es una fisura con un diámetro mínimo de entre 0,5 cm y 1 cm. Como se formaron no lo sabían, pero debieron de existir porque sin ellas no se podían explicar las cuevas.
Actualmente existe una teoría bastante sofisticada que explica la formación de los proto-conductos.
Este artículo trata de intentar de explicar algunos de los mecanismos que juegan un papel en el ensanchamiento de una fisura muy estrecha hasta un proto-conducto. Tratamos a los mecanismos según históricamente fueron descubiertos.
UN POCO DE QUÍMICA ELEMENTAL
La mayoría de las cuevas del mundo se encuentran en roca caliza. La caliza, compuesta principalmente por carbonato cálcico (CaCO3), es casi insoluble en agua pura. Sin
embargo, en realidad el agua se carga con el gas anhídrido carbónico (CO2),
presente en el aire y en el suelo, y una vez en solución se transforma
rápidamente en ácido carbónico (H2CO3).
En formula la reacción es la
siguiente:
(1) H2O +
CO2 ↔ H2CO3
(agua) +
(anhídrido carbónico) ↔ (ácido carbónico)
La presencia de ácido carbónico hace que la roca caliza si
se disuelve ligeramente en agua, formando bicarbonato de calcio.
En formula la reacción es la siguiente:
(2) H2CO3 +
CaCO3 ↔ Ca2+ + 2HCO3¯
(ácido carbónico) + (carbonato cálcico) ↔
(bicarbonato de calcio)
El bicarbonato de calcio también se suele escribir como
Ca(HCO3)2 . Aquí hemos elegido la anotación en iones porque
es la cantidad de calcio disuelto en agua (Ca2+) que condiciona
todos los mecanismos.
Podemos combinar ambas reacciones en una sola formula:
(3) H2O +
CO2 + CaCO3 ↔
Ca2+ + 2HCO3¯
(agua)
+ (anhídrido carbónico) + (carbonato cálcico) ↔ (bicarbonato de calcio)
Esta formula indica que cada
molécula de anhídrido
carbónico puede disolver una molécula de carbonato cálcico.
Hay que darse cuenta que se trata de una reacción de
equilibrio, es decir no todo el ácido carbónico se convierte en bicarbonato de
calcio. La reacción funciona hacia los dos lados y llega un momento que se
convierte la misma cantidad de bicarbonato de calcio en ácido carbónico que
ácido carbónico en bicarbonato de calcio, es decir que ya no se puede disolver
más caliza. Se dice que en este momento la solución de agua ha llegado a su
concentración de equilibrio y que está saturada.
En la formula 2 el ácido carbónico está disuelto en el agua
y el carbonato de calcio es la caliza que forma las paredes de una grieta o
cueva. El bicarbonato no es otra cosa que la cantidad de caliza disuelto por el
ácido carbónico. Cuando alguna vez se habla de la velocidad de disolución
(siempre respecto a la caliza), hay que darse cuenta que esto es lo mismo que
la velocidad con que se retroceden las paredes de la grieta o sea, lo rápido
que se ensancha la apertura (diámetro) de una fractura.
Para no complicarnos demasiado la vida con la química, a
partir de ahora llamamos al carbonato cálcico simplemente caliza e indicamos el
bicarbonato con Ca2+.
EL GRAN ENIGMA
Hemos visto que el agua cargada de H2CO3
es capaz de disolver cierta cantidad de caliza. ¿Que pasa cuando este agua
entra en las estrechas fisuras que se encuentran en la roca caliza por debajo
del suelo? Pues estas fisuras se ensancharán, pero la velocidad con que el agua
puede pasar por fisuras muy estrechas (alrededor de 0,2 mm de anchura) es tan
baja que el agua llega a saturarse, respecto al Ca2+, casi por
completo a solo unos metros de profundidad y por eso apenas puede disolver más
caliza..
¿Entonces como es posible que se han
formado cuevas, a menudo a cientos de metros de profundidad?
Las antiguas teorías sobre la formación de cuevas partieron
de la suposición de la existencia de fisuras lo bastante anchas como para dejar
pasar un flujo de agua considerable *1. De donde vinieron estas
fisuras era un enigma, pero debieron de existir porque sin ellas no se podían
explicar las cuevas. Estas fisuras se llamaban “proto conductos” (proto conduits).
*1 La
anchura mínima de una fisura que se necesita para que un flujo de agua puede
penetrar profundamente en una roca caliza sin perder la mayoría de su capacidad
de disolver caliza, son varios milímetros, mientras que la verdadera anchura de
las fisuras es desde 0,2 mm hacia abajo.
Hay que darse cuenta que se trata de una reacción de equilibrio, es decir no todo el ácido carbónico se convierte en bicarbonato de calcio. La reacción funciona hacia los dos lados y llega un momento que se convierte la misma cantidad de bicarbonato de calcio en ácido carbónico que ácido carbónico en bicarbonato de calcio, es decir que ya no se puede disolver más caliza. Se dice que en este momento la solución de agua ha llegado a su concentración de equilibrio y que está saturada.
CORROSIÓN DE
MEZCLAS
En el año 1964 el suizo Alfred Bögli publicó su famoso
principio de la “corrosión de mezclas” (mixing corrosion). Bögli demostró que
la mezcla de dos aguas completamente saturadas, pero con diferentes
concentraciones de Ca2+, posee otra vez una cierta capacidad de
disolución (figura 1) *2. Era la explicación perfecta para la
milagrosa reaparición de la capacidad de disolución a grandes profundidades.
Especialmente en los
años setenta casi todos los aspectos de la formación de
cuevas se explicaban con este principio. Sin embargo, incluso con la “corrosión
de mezclas” no se podía explicar satisfactoriamente el ensanchamiento de
estrechas fisuras hasta proto-conductos.
*2 Explicación de la figura 1: la línea negra indica le relación
entre la cantidad de carbonato cálcico (CaCO3) y la cantidad de anhídrido carbónico
(CO2) disueltas, cuando la solución está en equilibrio. Todo el área
por debajo de esta línea representa soluciones que no están saturadas. La línea
roja indica las posibles mezclas entre las soluciones saturadas A y B. En esta
caso la mezcla está indicado como C1 (se observa que C1
se
encuentra un poco más cerca de B, lo que indica que la solución B aporta un
poco más de agua que la solución A). El punto C1 se encuentra por debajo de
la línea negra lo que significa que no está saturada y por lo tanto puede
disolver una cantidad de CaCO3 adicional.
Cuando esta solución llega a saturarse otra ve,z su composición química está
indicado por C2 (situado sobre la línea
negra). La cantidad adicional que se ha disuelta equivale a la diferencia entre
X1 y X2.
El trayecto inclinado de C1-C2 tiene el siguiente
significado: Cuando se disuelve algo de
CaCO3 adicional, se consume un poco de CO2 (en cantidades
moleculares iguales como indica la formula 3). Esta disminución de CO2 se
manifiesta en un desplazamiento hacia la izquierda.
.
La
figura 1 muestra la cantidad de CaCO3 que se puede disolver
adicionalmente cuando dos soluciones saturadas y de diferentes concentraciones
se mezclan.
EL RÉGIMEN DE DISOLUCIÓN LENTA
En 1977 W.B. White publicó su trabajo acerca del comportamiento de una solución de agua cuando casi está saturada respecto al Ca2+. ¿Qué es exactamente lo que descubrió?
White descubrió que llega un momento que la rápida disolución de la roca caliza se ralentiza enormemente cuando la solución se acerca a su saturación. Es decir hay un régimen de rápida disolución y hay un régimen de disolución lenta. El cambio entre estos dos regímenes ocurre cuando la saturación llega a aproximadamente un 90 por ciento (figura 2). Este cambio abrupto en el comportamiento de la disolución es debido a ciertas impurezas, como fosfatos y silicatos, que se encuentran en las calizas naturales. Hasta este momento se había pensado que la velocidad con que una solución de agua se satura era rápida y que la saturación era completa. Es decir que el agua perdía toda su capacidad de disolución en los primeros metros de una estrecha fisura. Ahora se dieron cuenta que la velocidad con que se disuelve la caliza disminuye tanto (a partir de una saturación de 90 por ciento), que la solución no llega a saturarse por completo nunca.
La consecuencia es que una solución de agua sigue manteniendo una pequeña capacidad de disolución, incluso después de haber penetrado por una estrecha fisura hasta una distancia de varios kilómetros.
En resumen: White había descubierto que una solución nunca pierda toda su capacidad de disolución. Es verdad que esta capacidad se hace muy pequeña cuando la solución se acerca a una saturación casi completa, pero no se pierde totalmente, ni siquiera después de haber penetrado por una estrecha fisura hasta una distancia de varios kilómetros.
La figura 2 muestra la velocidad de disolución de Ca2+
respecto al porcentaje de saturación. Hasta un 90 por ciento la velocidad es alta y la caliza de las paredes de la fisura se disuelve rápidamente. Si la saturación supera el 90 por ciento, la velocidad de disolución baja enormemente y se disuelve muy poca caliza. Por eso la solución no llega nunca a saturarse por completa, manteniendo una pequeña capacidad de disolución. Un flujo de agua sigue disolviendo (un poco de) caliza, incluso después de haber penetrado por una estrecha fisura hasta una distancia de varios kilómetros. La unidad de velocidad de disolución es mol/cm2.
En números: si un flujo de agua se mantiene constante en el tiempo, el diámetro de una fisura se puede doblar en unos 100.000 años.
Al final se había encontrado un mecanismo capaz de explicar el ensanchamiento de una fisura muy estrecha, y por toda su longitud y además en un tiempo geológicamente aceptable.
LOS MODELOS NUMÉRICOS
Sin embargo, no fue hasta principios de los años noventa que aparecieron los primeros modelos numéricos acerca del desarrollo de una fisura estrecha hacia un conducto penetrable (Dreybrodt 1990, Palmer 1991). Con estos modelos, basados en simulaciones matemáticas hechos por ordenador, se pudieron cuantificar los procesos de disolución. Por ejemplo, se puede calcular con bastante exactitud el tiempo que necesita una aislada y estrecha grieta, para convertirse en una galería ancha. Ahora vamos a describir los procesos que actúan en el ensanchamiento de una fisura, partiendo de la situación mostrada en la figura 3.
La figura 3 muestra un macizo de caliza. Un flujo de agua atraviesa este macizo por una estrecha fisura de longitud L. El agua entra en punto A y sale en punto B, y la diferencia en altura entre los dos puntos es Dh.
Un arroyo llega a tener contacto con una formación de caliza. Cierta cantidad de agua entra en la caliza (en punto A) por una estrecha fisura y sale otra vez a una distancia L de 1000 metros (en punto B). La anchura (diámetro) de la fisura es de 0,2 mm y la diferencia en altura entre A y B es de 50 metros. La figura 3 represente esta situación.
MÁS CAUDAL, MÁS DISOLUCIÓN
La diferencia en altura (∆h) causa una presión hidrostática que “empuja” un flujo por la fisura. Cuando el agua fluye por una fisura, la caliza es disuelta muy lentamente de sus paredes y por eso la fisura se ensancha un poco. Como consecuencia puede pasar algo más de agua por esta fisura lo que hace que también se disuelve un poco más de caliza de sus paredes. De esta manera el proceso entra en un “bucle positivo”, es decir: cuanto más cantidad de agua pasa por la fisura, cada vez mayor es la cantidad de caliza que se disuelve de sus paredes por segundo.
En resumen: si un litro de agua disuelve una porción de caliza, pues 2 litros de agua disuelven 2 porciones de caliza.
Por si solo, este mecanismo puede doblar el diámetro de una fisura en bastante menos de 100.000 años, pero además existe un segundo mecanismo de igual importancia.
MÁS VELOCIDAD DEL FLUJO, MAS DISOLUCIÓN
Cuanto más ancha es una fisura, mayor es la velocidad con que fluye el agua (por razones de disminución de resistencia al flujo). Entonces, mientras sigue el proceso de ensanchamiento de la fisura, el agua permanece cada vez menos tiempo dentro de la grieta (porque fluye más rápido) y por eso sale con un menor porcentaje de saturación. En la figura 2 observamos que cuando una solución esta cerca de saturación total (concentración de equilibrio), una pequeña disminución en el porcentaje de saturación aumenta muchísimo su velocidad de disolución de caliza. Como consecuencia el flujo de agua es capaz de disolver la caliza más rápido, lo que aumenta la velocidad de ensanchamiento de la fisura *3.
En resumen: si una corriente de agua disuelve una porción de caliza por segundo, pues una corriente de agua que fluye algo más rápida (y por eso está un poco menos saturado) disuelve más de una porción por segundo.
*3 Si el agua sale de la fisura con un porcentaje menor de saturación, entonces cada litro de agua disuelve algo menos de caliza, lo que parece indicar una disminución de la velocidad de disolución. Sin embargo, en realidad cada litro solo disuelve una pequeña fracción menos, mientras el aumento en tanto la anchura de la fisura como la velocidad del flujo causan que pasa mucho más agua por la fisura por segundo. De este modo, el flujo en su totalidad puede disolver mucho más caliza (el efecto negativo es totalmente superado por los efectos positivos).
RÉGIMEN DE DISOLUCIÓN RÁPIDA POR TODA LA FRACTURA
Los dos mecanismos juntos hacen que el diámetro de la fisura se ensancha cada vez más rápido , que cada vez pasa una mayor cantidad de agua por la grieta y que la concentración de Ca2+sea cada vez menor.
Pues puede llegar el momento que en la salida de la fisura (punto B) el porcentaje de saturación baja hasta un 90 por ciento, imponiendo un régimen de disolución rápida por toda la fisura. Esto significa que a partir de ahora la disolución es máxima e igual en toda la longitud de la fisura. Desde este momento la fisura puede ensancharse con una velocidad de hasta 0,1 mm por año, lo cual significa que una fisura en estas condiciones se puede ensanchar hasta un conducto con un diámetro de 2 metros en tan solo 10.000 años.
COMBINACIÓN DEL RÉGIMEN DE DISOLUCIÓN LENTA Y RÁPIDA
Ahora es cuando realmente se puede apreciar el significado del descubrimiento de White. Un régimen de disolución rápida no es capaz de abrir una fisura estrecha, porque la solución llega enseguida a saturarse. Con un régimen de disolución lenta no se puede llegar a formar una cueva en un tiempo geológicamente aceptable. Pues ninguno de los dos es eficaz por si solo. Sin embargo, la combinación de los dos nos ofrece un modelo muy útil: la disolución lenta abre la fisura hasta el tamaño de un proto-conducto y la disolución rápida lo aumenta hacia el tamaño de una cueva. La duración de la fase de disolución lenta es de decenas de miles hasta cientos de miles de años y la duración de la fase de disolución rápida es de decenas de miles de años.
EVENTO DE RUPTURA
El rápido ensanchamiento de la fisura también aumenta enormemente la cantidad de agua (Q) que puede pasar por la fisura. La figura 4 muestra que el aumento es muy abrupto y este proceso es llamado “evento de ruptura” (breakthrough) y el tiempo que ha transcurrido desde el inicio del ensanchamiento de la fisura hasta el momento en que ocurre el evento de ruptura se llama “tiempo de ruptura” (breakthrough time).
El concepto de evento de ruptura es muy importante, porque indica, entre otros, que un conducto (fisura) está creciendo a máxima velocidad. También indica que se ha formado un proto-conducto (con un diámetro de alrededor de 0,5cm) y que empieza a fluir tanto agua por el conducto, que muy probablemente no se puede mantener el régimen de la presión hidrostática. Se puede decir que a partir del evento de ruptura la cueva realmente empieza a formarse.
En nuestro ejemplo (figura 3) el tiempo de ruptura es de unos 19.000 años (figura 4) y para formar una galería de 2 metros de diámetro solo pasarían 29.000 años desde el momento que el agua empezó a entrar en la estrecha fisura inicial.
La figura 4 muestra la cantidad de agua que pasa por una fisura respecto al tiempo. TB indica indica el tiempo de ruptura. Al principio el aumento es lento, lo que corresponde con un lento ensanchamiento de la fisura. Llega un momento que el flujo de agua aumenta muchísimo en muy poco tiempo, lo cual corresponde con el tiempo de ruptura. En este momento la velocidad de ensanchamiento es constante con un máximo de aproximadamente 0,1mm por año. Finalmente el aumento de la cantidad de agua que pasa por la fisura baja mucho, lo cual corresponde con el hundimiento de la presión hidrostática.
UN
POCO DE MATEMÁTICAS
En
1996 Dreybrodt era capaz de derivar una fórmula para determinar el tiempo de
ruptura. La formula parece bastante complicada y solo se menciona aquí porque
demuestra el comportamiento exacto del tiempo de ruptura (TB) en
relación con el diámetro original de la fisura (a0), la longitud de la fisura
(L) y la presión hidrostática (h).
1)
Resulta que TB ≈ a0-3 lo que significa que para
las fisuras realmente muy estrechas el tiempo de ruptura aumenta tanto que hace
imposible la formación de proto conductos. Por ejemplo, si la anchura de una
fisura es 10 veces más estrecha, el tiempo de ruptura se aumenta por 1000. En
nuestro ejemplo (figura 2) significaría que para una fisura con un diámetro de
0,02 mm el TB sería unos 19.000.000 años, un tiempo geológicamente
no muy probable.
2)
Resulta que TB ≈ L8/3 lo que significa que el tiempo de
ruptura se aumenta de un modo muy significante respecto a un aumento en
longitud. Por ejemplo, si la longitud es 3 veces mayor, el tiempo de ruptura se
aumentará con 18,7 veces. En nuestro ejemplo significaría que para una fisura
con una longitud de 3 kilómetros el TB sería unos 350.000 años, un
tiempo geológicamente todavía aceptable.
3)
Resulta que TB ≈ h-4/3 lo que significa que el tiempo de
ruptura se aumenta respecto a una disminución de la presión hidrostática. Por
ejemplo, si la presión es 3 veces menor, el tiempo de ruptura se aumentará con
4,3 veces. En nuestro ejemplo
significaría que para una diferencia de altura entre la entrada y la
salida (h) de unos 17 metros el TB sería unos 82.000 años, un tiempo
geológicamente muy aceptable.
Otro
aspecto es la cantidad de H2CO3 disuelto en el agua
(concentración de equilibrio). Si esta cantidad es alta, la solución puede
disolver más caliza antes de llegar a saturarse. En general, si el agua que
entra en la roca caliza ha pasado por un suelo bien desarrollado, el contenido
de H2CO3 es
bastante mayor que el contenido de un agua que infiltra directamente en una
roca caliza descubierta (karst desnudo). En un karst desnudo el tiempo de
ruptura se multiplica por lo menos con un factor 10.
RESUMEN
Agua que contiene el gas anhídrido carbónico es capaz de disolver el carbonato cálcico (caliza). Esto es un proceso rápido y enseguida la solución llega a estar casi saturada.
Entonces si las cuevas se han formado por la disolución de la roca caliza, ¿cómo es posible que existan cuevas a gran profundidad? Para explicarlo se suponían la existencia de unas fracturas relativamente muy anchas y lo llamaban proto-conductos. No sabían cómo se habían formado. Cuando en 1964 Alfred Bögli publicó su descubrimiento de la corrosión de las mezclas, solucionó numerosos problemas acerca de la génesis de cuevas. Sin embargo, todavía no se podían explicar la formación de los proto-conductos, hasta que en 1977 William White descubrió el régimen de la disolución lenta. Resulta que una solución no pierde nunca toda su capacidad de disolución, independiente de la distancia que ha penetrado en una fisura. Los mecanismos de Más Caudal-Más Disolución y Más Velocidad del Flujo- Más Disolución ya eran conocidos y por eso entendieron que incluso una muy pequeña capacidad de disolución es capaz de ensanchar una fisura estrecha hasta un diámetro que permite la disolución rápida por toda esta fisura. Finalmente se pudieron explicar la existencia de los proto-conductos. Además, la combinación del régimen de disolución rápida (la de toda la vida) y el régimen de disolución muy lenta (descubierto por White) era perfecta. La disolución lenta abre la fisura hasta el tamaño de un proto-conducto y a partir de aquí la solución rápida lo abre hasta formar una cueva. Con la mejora de los ordenadores en los años 90, los investigadores han sido capaz de crear modelos numéricos con lo cual han podido simular los procesos de disolución. En 1996 Dreybrodt ha derivado la formula para determinar el tiempo de ruptura, mostrando su relación con el diámetro original de la fisura (a0), la longitud de la fisura (L) y la presión hidrostática (h).
NOTA FINAL
En este artículo se ha intentado de explicar el concepto básico de cómo un proto-conducto se puede formar (con una aproximación en cuanto puede tardar) partiendo de una fisura muy estrecha. Se han omitido muchos aspectos que no son directamente necesarios para poder entender el concepto, pero que si son importantes. Algunos de estos aspectos son: el grado de saturación del agua que entra en la fisura, los parámetros químicos que determinan la disolución de la caliza, el régimen del flujo de agua (laminar o de turbulencia), la distribución del diámetro de la fisura a lo largo del conducto y la distribución de la presión hidrostática dentro del conducto.
La figura 2 muestra la situación geológica hace 3 millones de años. Los ríos que salieron de la montaña desembocaban en un grande lago interior.
La figura 5A muestra esquemáticamente la distribución del porcentaje de saturación en un momento dado. La figura 5B muestra la distribución de saturación en la misma fractura en un tiempo posterior. El punto S indica el cambio del régimen de disolución rápida hacia el régimen de disolución lenta. El porcentaje del cambio es el 90 por ciento. Las flechas negras indican la dirección y velocidad del flujo. Las líneas rojas indican puntos correspondientes. La fractura no está dibujada a escala y los porcentajes de saturación solo son orientativos.
La figura 2A muestra una galería con un arroyo vadoso, que desemboca en un valle al ras de un río. La figura 2B muestra que después de un tiempo el río se encaje y por el desnivel formado entre la galería y el río, empieza a fluir agua por las fracturas. Solo el flujo por las fracturas principales (1, 2 y 3) es indicado (flechas azules). La figura 2C muestra la formación de una galería inferior y el abandono de la galería principal desde el punto de pérdida (P1) hasta la salida. También muestra que el flujo que entra por la fractura 3 ha aumentado, consecuencia del aumento del gradiente hidráulico. El flujo en la galería inferior es freático. La figura 2D muestra el desarrollo de un nuevo punto de pérdida (P2) y el abandono del trayecto P1-P2. De este modo el abandono de la galería es progresivamente hacia dentro. El flujo en la galería inferior es una mezcla de vadoso y freático.
La figura 3 muestra el alcance de los tres regímenes de disolución. El intervalo de 30 a 36 por ciento muestra el cambio progresivo entre la disolución súper rápida y la disolución rápida.
