Earth
I ghiacciai modellano le valli
Riempi la valle a V di un fiume con un ghiacciaio spesso seicento metri e guarda passare centomila anni. Il ghiaccio scorre più veloce al centro della sua superficie e per nulla contro il letto; lo sforzo che trascina sul letto è massimo sui fianchi della V e si annulla nel suo stretto fondo, così scivolamento e abrasione arretrano i fianchi e la valle si allarga verso una U - la vecchia V resta disegnata per confronto. Il ghiaccio pulito taglia molto meno di quello carico di roccia; senza scivolamento non taglia nulla; fai avanzare e ritirare il ghiacciaio a cicli. Ogni metro cubo perso dal letto va alla morena. In un canale semicircolare la soluzione è quella chiusa di Poisson.

Cosa succede
Una sezione di valle larga 2400 m su una griglia di 61 per 40. La velocità a valle del ghiaccio newtoniano (10¹³ Pa·s, pendenza 0,05) risolve eta grad²u = −rho g sin alpha con u = 0 sul letto e du/dz = 0 in superficie; il letto irregolare è trattato con differenze di Shortley-Weller e la soluzione per sovrarilassamento dalla precedente. Lo sforzo sul letto guida uno scivolamento proporzionale - una semplificazione dichiarata, lo scivolamento non entra nella soluzione - e lo scivolamento abrade il letto (Hallet), scalato dai detriti nel ghiaccio e dalla resistenza della roccia. Ogni metro cubo va alla morena. Passi di 1 millennio; un secondo è mille anni a scala temporale 1.
Da provare
Guarda la V allargarsi verso una U e leggi l’esponente del profilo; confronta ghiaccio pulito e sporco; ferma lo scivolamento; fai avanzare e ritirare il ghiacciaio a cicli; verifica il semicerchio contro la formula chiusa.
Il modello e i suoi limiti
Una sezione senza flusso lungo la valle, ghiaccio newtoniano, scivolamento proporzionale allo sforzo calcolato, una sola legge d’abrasione; estrazione e fratturazione della roccia non sono modellate. I profili appartengono a questa valle idealizzata.
Equazioni del modulo
Equazioni η·∇²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


