Earth
Delta fluvial
Envoyez une rivière de 200 m³/s portant un centième de mètre cube de sable par seconde sur un plateau continental et regardez un siècle de delta croître : l’eau se partage au prorata de la pente vers chaque voisin plus bas, chaque cellule marine laisse tomber une fraction 1 − exp(−profondeur / 4 m) de ce qui passe et monte jusqu’à dominer la mer d’une pente douce, le rivage avance et la rivière trouve plusieurs bouches. Triplez le sédiment et le lobe recouvre tout le plateau ; approfondissez le plateau et le même sable construit moins de terre ; montez la mer d’un centimètre par an et le gain est encore plus faible. Le sédiment apporté égale le déposé plus l’exporté à l’arrondi près.

Ce qui se passe
Un delta sur une grille de 32 par 24 cellules de 100 m : la rivière entre en bas au milieu d’une plaine côtière inclinée, et son eau et son sédiment sont routés cellule par cellule par ordre de potentiel décroissant (le lit à terre, le niveau de la mer moins un cinquième de la profondeur en mer pour que le courant s’étale vers l’eau profonde) au prorata de la pente vers tous les voisins plus bas. Chaque cellule dépose une fraction 1 − exp(−profondeur/λ) de ce qui passe, au moins deux pour cent, jusqu’à un plafond qui descend doucement depuis la bouche pour que la plaine garde une pente vers la mer ; ce qui atteint le bord ou une dépression est exporté. Mise à jour synchrone du lit tous les 0,05 an ; apporté égale déposé plus exporté à l’arrondi près. Une seconde d’expérience vaut une année.
À essayer
Regardez le lobe avancer et les bouches se diviser ; triplez le sédiment, approfondissez le plateau, montez la mer et comparez la terre construite en un siècle.
Le modèle et ses limites
Routage multidirectionnel déterministe sans inertie, vagues ni marées, une règle de dépôt par la profondeur, pas d’incision de chenal ni de dynamique d’avulsion. Les surfaces de terre sont celles de cette règle cellulaire.
Équations du module
Équations η = 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


