Earth
Delta fluvial
Envía un río de 200 m³/s que lleva una centésima de metro cúbico de arena por segundo sobre una plataforma y mira crecer un siglo de delta: el agua se reparte en proporción a la pendiente hacia cada vecino más bajo, cada celda marina deja caer una fracción 1 − exp(−profundidad / 4 m) de lo que pasa y sube hasta quedar sobre el mar con pendiente suave, la costa avanza y el río encuentra varias bocas. Triplica el sedimento y el lóbulo cubre toda la plataforma; profundiza la plataforma y la misma arena construye menos tierra; sube el mar un centímetro al año y la ganancia es aún menor. El sedimento aportado iguala el depositado más el exportado salvo redondeo.
Qué ocurre
Un delta sobre una malla de 32 por 24 celdas de 100 m: el río entra abajo en el centro de una llanura costera inclinada, y su agua y sedimento se enrutan celda a celda en orden de potencial decreciente (el lecho en tierra, el nivel del mar menos un quinto de la profundidad en el mar para que la corriente se extienda hacia el agua profunda) en proporción a la pendiente hacia todos los vecinos más bajos. Cada celda deposita una fracción 1 − exp(−profundidad/λ) de lo que pasa, al menos un dos por ciento, hasta un techo que baja suavemente desde la boca para que la llanura conserve pendiente hacia el mar; lo que llega al borde o a una depresión se exporta. Actualización síncrona del lecho cada 0,05 años; aportado igual a depositado más exportado salvo redondeo. Un segundo de experimento es un año.
Qué probar
Mira avanzar el lóbulo y dividirse las bocas; triplica el sedimento, profundiza la plataforma, sube el mar y compara la tierra construida en un siglo.
El modelo y sus límites
Enrutamiento multidireccional determinista sin inercia, olas ni mareas, una regla de depósito por profundidad, sin incisión de cauce ni dinámica de avulsión. Las superficies de tierra son las de esta regla celular.
Ecuaciones del módulo
Ecuaciones η = 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
