Earth
Pergélisol et dégel
Regardez le sol sous une maison de toundra au fil des saisons : chaque été le haut dégèle jusqu’à la profondeur de la couche active, chaque hiver il regèle, et en dessous le sol reste gelé au point de fusion de sa glace. Une épaisse couche de neige retient le froid de l’hiver et garde le sol cinq degrés plus chaud ; raclez-la et le pergélisol se durcit. Réchauffez le climat d’un degré par décennie et la couche active s’approfondit d’année en année jusqu’à ce que la couche riche en glace dégèle et que la maison s’affaisse de sa glace en excès. La chaleur latente de la glace retient le front de dégel exactement comme le dit la solution de Neumann, et la chaleur passée par la surface et la base est l’enthalpie gagnée par la colonne.

Ce qui se passe
Une colonne de sol de 15 m en 75 cellules, 35 % d’eau avec une couche riche en glace à 70 % entre 1,5 et 3 m. Chaque cellule tient son enthalpie relative au sol gelé à 0 °C : en dessous de zéro gelée, entre zéro et la chaleur latente de son eau un mélange à 0 °C, au-delà dégelée ; la conduction utilise des moyennes harmoniques des conductivités gelée et dégelée. L’air (−6 °C de moyenne, ±18 °C, une tendance au réchauffement) pilote le haut à travers 0,3 m de neige en hiver ; 60 mW/m² de chaleur terrestre entrent par la base. La glace en excès tasse au premier dégel. Pas explicites de 3 h ; la chaleur par le haut et la base égale l’enthalpie gagnée, et la solution de Neumann vérifie le front de dégel.
À essayer
Regardez une année de dégel et de gel ; raclez la neige ou épaississez-la ; réchauffez le climat d’un degré par décennie et attendez que la maison s’affaisse ; coupez la chaleur latente et voyez à quelle vitesse court le front.
Le modèle et ses limites
Conduction unidimensionnelle, un point de congélation net, des propriétés de sol fixes, sans migration d’eau, gonflement au gel ni végétation ; le tassement est une comptabilité déclarée de la glace en excès. Les profondeurs sont celles de cette colonne idéalisée.
Équations du module
Équations 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


