Earth
Permafrost e disgelo
Guarda il terreno sotto una casa della tundra lungo le stagioni: ogni estate la parte alta disgela fino alla profondità dello strato attivo, ogni inverno rigela, e sotto il terreno resta gelato al punto di fusione del suo ghiaccio. Uno spesso manto nevoso tiene fuori il freddo invernale e il terreno cinque gradi più caldo; toglilo e il permafrost si rafforza. Scalda il clima di un grado per decennio e lo strato attivo si approfondisce anno dopo anno finché lo strato ricco di ghiaccio disgela e la casa cede quanto il suo ghiaccio in eccesso. Il calore latente del ghiaccio trattiene il fronte di disgelo esattamente come dice la soluzione di Neumann, e il calore passato per la superficie e la base è l’entalpia guadagnata dalla colonna.

Cosa succede
Una colonna di terreno profonda 15 m in 75 celle, al 35 % d’acqua con uno strato ricco di ghiaccio al 70 % tra 1,5 e 3 m. Ogni cella tiene la sua entalpia rispetto al terreno gelato a 0 °C: sotto zero è gelata, tra zero e il calore latente della sua acqua una miscela a 0 °C, oltre disgelata; la conduzione usa medie armoniche delle conducibilità gelata e disgelata. L’aria (media −6 °C, ±18 °C, una tendenza al riscaldamento) agisce sulla superficie attraverso 0,3 m di neve d’inverno; 60 mW/m² del calore della Terra entrano dalla base. Il ghiaccio in eccesso cede la prima volta che disgela. Passi espliciti di 3 h; il calore per superficie e base eguaglia l’entalpia guadagnata, e la soluzione di Neumann verifica il fronte.
Da provare
Guarda il disgelo e il gelo di un anno; togli la neve o aumentala; scalda il clima di un grado per decennio e aspetta che la casa ceda; spegni il calore latente e guarda quanto corre il fronte.
Il modello e i suoi limiti
Conduzione unidimensionale, punto di congelamento netto, proprietà del terreno fisse, niente migrazione d’acqua, sollevamento da gelo né vegetazione; il cedimento è una contabilità dichiarata del ghiaccio in eccesso. Le profondità sono quelle di questa colonna idealizzata.
Equazioni del modulo
Equazioni 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


