Earth
Delta fluviale
Manda un fiume da 200 m³/s che porta un centesimo di metro cubo di sabbia al secondo su una piattaforma e guarda crescere un secolo di delta: l’acqua si divide in proporzione alla pendenza verso ogni vicino più basso, ogni cella marina lascia cadere una frazione 1 − exp(−profondità / 4 m) di ciò che passa e sale finché sta sopra il mare con pendenza dolce, la costa avanza e il fiume trova più bocche. Triplica il sedimento e il lobo copre tutta la piattaforma; approfondisci la piattaforma e la stessa sabbia costruisce meno terra; alza il mare di un centimetro all’anno e il guadagno è ancora minore. Il sedimento apportato eguaglia il deposto più l’esportato a meno dell’arrotondamento.

Cosa succede
Un delta su una griglia di 32 per 24 celle di 100 m: il fiume entra in basso al centro di una piana costiera inclinata, e la sua acqua e il suo sedimento sono instradati cella per cella in ordine di potenziale decrescente (il fondo a terra, il livello del mare meno un quinto della profondità in mare perché la corrente si spanda verso l’acqua profonda) in proporzione alla pendenza verso tutti i vicini più bassi. Ogni cella deposita una frazione 1 − exp(−profondità/λ) di ciò che passa, almeno il due per cento, fino a un tetto che scende dolcemente dalla bocca perché la piana conservi pendenza verso il mare; ciò che raggiunge il bordo o una depressione è esportato. Aggiornamento sincrono del fondo ogni 0,05 anni; apportato uguale a deposto più esportato a meno dell’arrotondamento. Un secondo di esperimento è un anno.
Da provare
Guarda il lobo avanzare e le bocche dividersi; triplica il sedimento, approfondisci la piattaforma, alza il mare e confronta la terra costruita in un secolo.
Il modello e i suoi limiti
Instradamento multidirezionale deterministico senza inerzia, onde né maree, una regola di deposito per profondità, niente incisione di alveo né dinamica di avulsione. Le aree di terra sono quelle di questa regola cellulare.
Equazioni del modulo
Equazioni η = z (z ≥ S), η = S − r·(S − z) (z < S) Q_ab = Q_a·w_ab/Σw, w = (η_a − η_b)/d_ab f = max(f_land, 1 − e^(−depth/λ)), Δz = min(f·q, cap)/A Q_s·t = V_deposited + V_exported


