Earth
Permafrost y deshielo
Mira el suelo bajo una casa de la tundra a lo largo de las estaciones: cada verano la parte superior se deshiela hasta la profundidad de la capa activa, cada invierno vuelve a helarse, y debajo el suelo sigue helado en el punto de fusión de su hielo. Una capa de nieve gruesa aísla del frío invernal y deja el suelo cinco grados más cálido; quítala y el permafrost se endurece. Calienta el clima un grado por década y la capa activa se hace más profunda año tras año hasta que la capa rica en hielo se deshiela y la casa se asienta lo que mide su hielo en exceso. El calor latente del hielo frena el frente de deshielo justo como dice la solución de Neumann, y el calor que cruzó la superficie y la base es la entalpía que ganó la columna.

Qué ocurre
Una columna de suelo de 15 m en 75 celdas, con 35 % de agua y una capa rica en hielo del 70 % entre 1,5 y 3 m. Cada celda guarda su entalpía respecto al suelo helado a 0 °C: bajo cero está helada, entre cero y el calor latente de su agua es una mezcla a 0 °C, por encima, deshelada; la conducción usa medias armónicas de las conductividades helada y deshelada. El aire (−6 °C de media, ±18 °C, una tendencia de calentamiento) actúa sobre la superficie a través de 0,3 m de nieve en invierno; 60 mW/m² del calor de la Tierra entran por la base. El hielo en exceso se asienta la primera vez que se deshiela. Pasos explícitos de 3 h; el calor por superficie y base iguala la entalpía ganada, y la solución de Neumann comprueba el frente.
Qué probar
Mira el deshielo y la helada de un año; quita la nieve o hazla más gruesa; calienta el clima un grado por década y espera a que la casa se asiente; apaga el calor latente y mira qué rápido corre el frente.
El modelo y sus límites
Conducción unidimensional, punto de congelación neto, propiedades del suelo fijas, sin migración de agua, levantamiento por helada ni vegetación; el asentamiento es una contabilidad declarada del hielo en exceso. Las profundidades son las de esta columna idealizada.
Ecuaciones del módulo
Ecuaciones H = C_f·T (T < 0), 0…L (T = 0), L + C_u·T (T > 0) dH/dt = ∂/∂z(k·∂T/∂z), q0 = (T_air − T1)/(R_snow + Δz/2k) X(t) = 2λ·√(α·t), λ·e^(λ²)·erf λ = St/√π ΔH_column = ∫q0 dt + ∫q_geo dt, s = Σ(w_ice − w)·Δz


