Earth
Metamorfismo
Enterra una lutita a un kilómetro por millón de años a lo largo de una geoterma de 25 K/km, mantenla y devuélvela: su punto recorre la malla presión-temperatura de las facies, la clorita cede a biotita, granate y estaurolita a medida que las fracciones de equilibrio se desplazan, y como la velocidad de reacción cae mil veces entre 600 y 400 grados, la paragénesis del pico vuelve a la superficie casi intacta, algo más retrogradada si la exhumación es lenta. Empina el enterramiento bajo un gradiente frío y el camino entra en los esquistos azules; deja que un plutón caliente una roca somera y se cuece en corneana; mantén una muestra isoterma y sus fracciones relajan sobre la exponencial cerrada.
Qué ocurre
Un paquete de rocas en un camino de enterramiento, parada y exhumación: presión de la sobrecarga, temperatura de una geoterma lineal con un pulso de contacto opcional que decae en medio millón de años, una malla esquemática de facies en temperatura y presión y cinco minerales índice pelíticos cuyas fracciones de equilibrio son ventanas de temperatura suaves de suma uno, relajadas hacia el equilibrio a una velocidad de Arrhenius de 10 por Ma a 600 °C y 200 kJ/mol. La relajación es exacta en cada subpaso, de modo que la paragénesis del pico sobrevive a una exhumación rápida mientras una lenta la retrograda un poco, y un paso isotermo cae en la forma cerrada. Un segundo de experimento es un millón de años.
Qué probar
Sigue el punto por la malla de facies y la barra de minerales que se desplaza; compara un retorno rápido y uno lento; recorre un camino frío y profundo hacia los esquistos azules y uno somero y caliente hacia la corneana.
El modelo y sus límites
Una malla de facies esquemática y cinco minerales con ventanas suaves, una velocidad de Arrhenius para toda reacción, una geoterma lineal sin conducción transitoria, sin fluidos ni deformación. Un camino didáctico, no una malla petrogenética.
Ecuaciones del módulo
Ecuaciones P = ρ·g·z, T = T_s + G·z + ΔT_0·e^(−(t − t_c)/τ) df_i/dt = r(T)·(f_i^eq − f_i), r = r_600·exp(−(E/R)·(1/T − 1/873 K)) f_i^eq ∝ ¼·(1 + tanh((T − T_lo)/w))·(1 − tanh((T − T_hi)/w)), Σ f_i^eq = 1 f(t) = f^eq + (f_0 − f^eq)·e^(−r·t)
