Earth
Les glaciers façonnent les vallées
Remplissez la vallée en V d’une rivière d’un glacier de six cents mètres d’épaisseur et regardez passer cent mille ans. La glace coule le plus vite au milieu de sa surface et pas du tout contre le lit ; le cisaillement qu’elle exerce sur le lit est maximal sur les flancs du V et s’annule dans son fond étroit, si bien que glissement et abrasion entaillent les flancs et que la vallée s’élargit en U - l’ancien V reste tracé pour comparer. La glace propre creuse bien moins que la glace chargée de roches ; sans glissement elle ne creuse rien ; faites avancer et reculer le glacier par cycles. Chaque mètre cube perdu par le lit va à la moraine. Dans un chenal semi-circulaire la solution est celle de Poisson en formule fermée.

Ce qui se passe
Une section de vallée de 2400 m de large sur une grille de 61 sur 40. La vitesse vers l’aval d’une glace newtonienne (10¹³ Pa·s, pente 0.05) résout eta grad²u = −rho g sin alpha avec u = 0 sur le lit et du/dz = 0 en surface ; le lit irrégulier est traité par différences de Shortley-Weller et la solution par surrelaxation depuis la précédente. Le cisaillement sur le lit entraîne un glissement proportionnel - une simplification déclarée, le glissement n’entre pas dans la résolution - et le glissement abrase le lit (Hallet), selon les débris dans la glace et la résistance de la roche. Chaque mètre cube va à la moraine. Pas de 1 millier d’années ; une seconde vaut mille ans à l’échelle de temps 1.
À essayer
Regardez le V s’élargir en U et lisez l’exposant du profil ; comparez glace propre et glace sale ; arrêtez le glissement ; faites avancer et reculer le glacier par cycles ; vérifiez le demi-cercle contre la formule.
Le modèle et ses limites
Une section sans écoulement le long de la vallée, une glace newtonienne, un glissement proportionnel au cisaillement calculé, une seule loi d’abrasion ; arrachement et fracturation de la roche ne sont pas modélisés. Les profils appartiennent à cette vallée idéalisée.
Équations du module
Équations η·∇²u = −ρ·g·sin α, u = 0 on the bed, ∂u/∂z = 0 at the surface τ_b = η·|∇u| at the bed, u_b = c_s·τ_b dz_b/dt = −K·D/S·u_b·√(1 + (∂z_b/∂y)²) semicircle: u = ρ·g·sin α·(R² − r²)/(4η), τ_b = ρ·g·sin α·R/2


