Earth
Metamorfizm
Pogrzeb mułowiec z prędkością kilometra na milion lat wzdłuż geotermy 25 K/km, przytrzymaj go i wynieś z powrotem: jego punkt wędruje po siatce ciśnienie–temperatura facji, chloryt ustępuje biotytowi, granatowi i staurolitowi, gdy udziały równowagowe się przesuwają, a ponieważ szybkość reakcji spada tysiąckrotnie między 600 a 400 stopniami, paragenezą szczytową wraca na powierzchnię niemal nietknięta, trochę bardziej retrogresywna przy wolnej ekshumacji. Przystrom pogrzebanie pod zimnym gradientem, a ścieżka wejdzie w łupki niebieskie; pozwól plutonowi ogrzać płytką skałę, a wypiecze się na hornfels; trzymaj próbkę izotermicznie, a jej udziały relaksują na zamkniętej eksponencie.
Co się dzieje
Pakiet skał na ścieżce pogrzebania, postoju i ekshumacji: ciśnienie z nadkładu, temperatura z liniowej geotermy z opcjonalnym impulsem kontaktowym zanikającym w pół miliona lat, schematyczna siatka facji w temperaturze i ciśnieniu oraz pięć pelitowych minerałów wskaźnikowych, których udziały równowagowe są gładkimi oknami temperatury o sumie jeden, relaksowane do równowagi z szybkością Arrheniusa 10 na mln lat przy 600 °C i 200 kJ/mol. Relaksacja jest dokładna w każdym podkroku, tak że paragenezą szczytowa przeżywa szybką ekshumację, a wolna trochę ją retrogresuje, a krok izotermiczny trafia w postać zamkniętą. Sekunda doświadczenia to milion lat.
Co wypróbować
Śledź punkt po siatce facji i pasek minerałów przy przesuwaniu; porównaj szybki i wolny powrót; przejdź zimną głęboką ścieżkę w łupki niebieskie i płytką gorącą w hornfels.
Model i jego ograniczenia
Schematyczna siatka facji i pięć minerałów o gładkich oknach, jedna szybkość Arrheniusa dla każdej reakcji, liniowa geoterma bez przewodzenia nieustalonego, brak fluidów i deformacji. Ścieżka dydaktyczna, nie siatka petrogenetyczna.
Równania modułu
Równania 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)
