Earth
Permafrost und Tauen
Sehen Sie den Boden unter einem Tundrahaus durch die Jahreszeiten: jeden Sommer taut die Oberschicht bis zur Tiefe der Auftauschicht, jeden Winter friert sie wieder, und darunter bleibt der Boden am Schmelzpunkt seines Eises gefroren. Eine dicke Schneedecke hält die Winterkälte ab und den Boden fünf Grad wärmer; schieben Sie ihn weg, und der Permafrost wird härter. Erwärmen Sie das Klima um ein Grad je Jahrzehnt, und die Auftauschicht vertieft sich Jahr für Jahr, bis die eisreiche Schicht taut und das Haus um sein Überschusseis sackt. Die latente Wärme des Eises hält die Tauffront genau so zurück, wie Neumanns Lösung sagt, und die Wärme durch Oberfläche und Basis ist die Enthalpie, die die Säule gewann.

Was passiert
Eine Bodensäule von 15 m in 75 Zellen, 35 % Wasser mit einer eisreichen Schicht von 70 % zwischen 1,5 und 3 m. Jede Zelle hält ihre Enthalpie relativ zu gefrorenem Boden bei 0 °C: darunter gefroren, zwischen null und der latenten Wärme ihres Wassers ein Gemisch bei 0 °C, darüber getaut; die Wärmeleitung nutzt harmonische Mittel der gefrorenen und getauten Leitfähigkeit. Die Luft (−6 °C Mittel, ±18 °C, ein Erwärmungstrend) treibt die Oberseite im Winter durch 0,3 m Schnee; 60 mW/m² Erdwärme treten an der Basis ein. Das Überschusseis setzt sich beim ersten Tauen. Explizite Schritte von 3 h; die Wärme durch Ober- und Unterseite gleicht der gewonnenen Enthalpie, Neumanns Lösung prüft die Tauffront.
Zum Ausprobieren
Sehen Sie ein Jahr Tauen und Frieren; schieben Sie den Schnee weg oder vertiefen Sie ihn; erwärmen Sie das Klima um ein Grad je Jahrzehnt und warten Sie, bis das Haus sackt; schalten Sie die latente Wärme ab und sehen Sie, wie schnell die Front läuft.
Das Modell und seine Grenzen
Eindimensionale Wärmeleitung, ein scharfer Gefrierpunkt, feste Bodeneigenschaften, keine Wanderung des Wassers, kein Frosthub, keine Vegetation; die Setzung ist eine erklärte Buchführung des Überschusseises. Tiefen sind die dieser idealisierten Säule.
Gleichungen des Moduls
Gleichungen 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


