¿Alguna vez te has preguntado por qué algunos minerales se forman en las profundidades volcánicas y otros en el fondo del océano? ¿O cómo es posible predecir si una reacción química ocurrirá espontáneamente en un acuífero subterráneo? La respuesta no está solo en la química, sino en la termodinámica geoquímica: el puente que conecta la energía, la temperatura, la presión y la composición química de la Tierra.
En los primeros 90 segundos de lectura comprenderás la idea central: la termodinámica geoquímica aplica las leyes de la termodinámica a sistemas naturales complejos (rocas, fluidos, magma) para predecir qué fases minerales son estables, qué reacciones ocurren y cómo se distribuyen los elementos químicos en nuestro planeta. Si dominas estos principios, podrás interpretar desde la génesis de un yacimiento de cobre hasta el comportamiento del carbono en el manto terrestre.
Ahora profundicemos paso a paso, con ejemplos concretos y aplicaciones reales que todo estudiante de geología, geoquímica o ingeniería geoambiental debería manejar.
¿Por qué la termodinámica y no solo la cinética?
Un error común es confundir lo que puede ocurrir (termodinámica) con lo rápido que ocurre (cinética). La termodinámica nos dice si una reacción es espontánea bajo ciertas condiciones de presión (P) y temperatura (T). La cinética nos dirá si esa reacción ocurre en segundos, años o millones de años.
En geoquímica, muchos procesos son extremadamente lentos (como la cristalización del granito), pero termodinámicamente están “permitidos”. Por eso, antes de estudiar velocidades, debemos dominar los elementos de equilibrio.
¿Qué es la hidrodinámica? Características, principios y ejemplos
Ejemplo clave: El diamante es metaestable en la superficie terrestre. Termodinámicamente debería convertirse en grafito, pero la barrera cinética es tan alta que nunca lo vemos cambiar. Sin embargo, en el manto profundo (>150 km), el diamante es la fase estable.
Conceptos fundamentales: sistema, entorno y variables
Para aplicar termodinámica a la Tierra, debemos definir nuestro sistema (porción del universo que estudiamos) y el entorno (todo lo demás). En geoquímica usamos principalmente tres tipos de sistemas:
- Abierto: intercambia masa y energía (ej. un géiser).
- Cerrado: intercambia energía pero no masa (ej. una inclusión fluida sellada).
- Aislado: no intercambia nada (ideal, raro en naturaleza).
Las variables clave en geotermodinámica son:
- T (temperatura): desde casi 0 K en el espacio hasta >6000 K en el núcleo interno.
- P (presión): desde 1 bar en superficie hasta ~360 GPa en el centro terrestre.
- V (volumen), composición química (actividades, fracciones molares).
Y las funciones de estado (energía interna U, entalpía H, entropía S, energía libre de Gibbs G) que no dependen del camino, solo del estado actual.
Las tres leyes aplicadas a procesos geológicos
Primera Ley: La energía se conserva
En cualquier proceso geológico (fusión de roca, evaporación de agua marina, metamorfismo), la variación de energía interna (ΔU) es igual al calor intercambiado (Q) menos el trabajo realizado (W). En geoquímica, el trabajo suele ser de expansión/compresión: W = P·ΔV.
Modelos de Ian Barbour: El Mapa para Entender Ciencia y Religión
Aplicación práctica: Cuando una placa oceánica se hunde (subducción), la presión aumenta drásticamente. El trabajo de compresión se convierte en calor, contribuyendo a la fusión parcial que genera magma andesítico.
Segunda Ley: La entropía del universo aumenta
La entropía (S) mide el desorden. En un sistema aislado, los procesos espontáneos son aquellos que aumentan la entropía total. Pero en geoquímica rara vez trabajamos con sistemas aislados; por eso usamos la energía libre de Gibbs (G = H – T·S).
Regla de oro: A P y T constantes (lo más común en la corteza y manto), una reacción es espontánea si ΔG < 0. Si ΔG = 0, hay equilibrio.
Tercera Ley: La entropía de un cristal perfecto es cero en el cero absoluto
Aunque en geología no alcanzamos 0 K, esta ley permite calcular entropías absolutas de minerales (valores tabulados en bases de datos como THERMOCALC o SUPCRTBL).
La energía libre de Gibbs: el predictor definitivo
La ecuación fundamental de la termodinámica geoquímica es:
Donde:
- ΔH = variación de entalpía (calor de reacción a P constante).
- ΔS = variación de entropía.
- T = temperatura absoluta (Kelvin).
Para reacciones que implican sólidos, líquidos y gases, debemos considerar también la presión. Para un gas ideal (como el CO₂ en una cámara magmática):
Para sólidos y líquidos puros, la actividad es 1 (fase pura). Para soluciones (ej. una aleación de olivino (Mg,Fe)₂SiO₄), usamos actividades (a_i) que dependen de la fracción molar y modelos de solución no ideal.
Caso real: Predecir si en una fuente hidrotermal a 250°C y 500 bar precipita calcita (CaCO₃) o se disuelve. Calculando ΔG de la reacción Ca²⁺ + CO₃²⁻ → CaCO₃(s), si ΔG < 0 → precipita.
Equilibrio de fases: diagramas P-T y regla de las fases de Gibbs
La regla de las fases de Gibbs para sistemas geológicos (sin considerar campos gravitatorios ni eléctricos) es:
Donde:
- F = grados de libertad (variables intensivas que podemos cambiar sin alterar el número de fases).
- C = número de componentes químicos independientes.
- F = número de fases (minerales + fluidos + fundido).
Aplicación: En un sistema metamórfico de pelitas (rocas ricas en Al, Si, O, H, K, Fe, Mg), con 5 componentes y 4 fases, los grados de libertad son 5-4+2=3 (podemos variar P, T y una composición sin desaparecer fases).
Los diagramas de fase (P-T) muestran campos de estabilidad de asociaciones minerales. Por ejemplo, el límite entre andalucita, silimanita y cianita (polimorfos del Al₂SiO₅) es clásico en geotermobarometría.
Geotermobarometría: midiendo P y T pasadas
Una de las aplicaciones más poderosas: usando reacciones de intercambio catiónico entre minerales, podemos calcular la presión y temperatura a la que se formó una roca.
Ejemplo: El par granate-biotita. El hierro y el magnesio se distribuyen entre ambos minerales según:
Esta constante de equilibrio varía sistemáticamente con T y P, mediante ecuaciones calibradas experimentalmente. Así, analizando un granito metamórfico con microsonda electrónica, obtenemos que se formó a 650°C y 8 kbar.
Actividad y fugacidad: cuando la realidad se aleja de la idealidad
En fluidos geológicos a alta presión, el comportamiento de gases no es ideal. Usamos fugacidad (f) como una “presión corregida”. Para agua pura a 1000 bar y 500°C, la fugacidad puede ser el doble de la presión real.
En soluciones sólidas (ej. plagioclasa: albita + anortita), la actividad (a) corrige la concentración real. Modelos comunes: solución ideal (a = X), regular (parámetros de interacción), o subregular (Van Laar).
Diagramas de actividad (log fO₂ – pH) y el mundo de los fluidos
En sistemas acuosos superficiales y subsuperficiales (acuíferos, lagos ácidos, minas), la termodinámica geoquímica se expresa con diagramas de Pourbaix (Eh-pH) o log fO₂ – pH. Estos muestran qué especie domina (Fe²⁺, Fe³⁺, hematita, pirita) según el oxígeno disponible y la acidez.
Aplicación ambiental: Predecir si el arsénico en un agua subterránea está como arseniato (movil) o como orpimente sólido (inmovilizado), ajustando el pH o añadiente oxidante.
Bases de datos termodinámicas y software
Ningún geoquímico moderno calcula ΔG a mano para cada mineral. Usamos bases de datos como:
- SUPCRT92/SUPCRTBL (fluidos y minerales hasta 1000°C, 5 kbar).
- THERMOCALC (metamorfismo de rocas).
- HCh (modelado de interacción agua-roca).
- PhreeqC (transporte reactivo en aguas subterráneas).
Estos programas contienen valores de ΔG°f, S°, H°, Cp para miles de especies. Con ellos simulamos desde la evaporación de agua marina hasta la alteración hidrotermal de un yacimiento de cobre porfídico.
Limitaciones y futuro de la termodinámica geoquímica
La termodinámica de equilibrio asume que el sistema alcanzó el mínimo de energía libre. Pero muchas rocas muestran texturas de desequilibrio (zonación, reemplazos parciales). Ahí debemos complementar con cinética o modelos de fraccionamiento de isótopos (que también se basa en ΔG de reacciones isotópicas).
Las fronteras actuales incluyen:
- Termodinámica de minerales nanoestructurados (diferentes propiedades por tamaño de grano).
- Efectos de la deformación mecánica sobre el equilibrio.
- Nuevos modelos para el núcleo terrestre a > 300 GPa y 6000 K, donde la teoría de funcionales de la densidad (DFT) complementa los datos experimentales.
Resultados de aprendizaje
Después de leer este artículo, el estudiante debería ser capaz de:
- Definir un sistema geoquímico y diferenciar entre sistema abierto, cerrado y aislado aplicado a procesos terrestres.
- Explicar la diferencia entre termodinámica y cinética, reconociendo que la espontaneidad (ΔG < 0) no implica velocidad rápida.
- Aplicar la primera y segunda ley de la termodinámica a procesos como subducción, fusión parcial y metamorfismo.
- Calcular la variación de energía libre de Gibbs (ΔG) para una reacción mineral simple a partir de datos tabulados de ΔH y S.
- Interpretar un diagrama de fase P-T (como el del Al₂SiO₅) para determinar condiciones de formación de rocas metamórficas.
- Utilizar la regla de las fases de Gibbs para predecir cuántas variables pueden cambiar independientemente en una asociación mineral.
- Distinguir entre fugacidad (para gases reales) y actividad (para soluciones sólidas y fluidos iónicos).
- Leer un diagrama de Pourbaix (Eh-pH) para predecir la movilidad de metales pesados en aguas naturales.
- Nombrar al menos tres bases de datos termodinámicas y un software de modelado geoquímico, describiendo su aplicación típica.
- Reconocer las limitaciones del equilibrio termodinámico en rocas naturales y cuándo es necesario incorporar cinética o modelos de no-equilibrio.
Explora más sobre este tema
Selecciona un tema y sigue aprendiendo...
