Earth
Cycle de l’eau
Laissez une année de temps passer sur un bassin fait de forêt, de ville et de lac : l’humidité arrive par le vent, pleut au prorata des surfaces, la canopée en intercepte un quart et le rend par évaporation, le sol forestier boit cinquante millimètres par jour quand le pavé en prend cinq, les sols évaporent et percolent, la nappe suinte vers la rivière en trois cents jours et la rivière se vide à la mer en deux. Faites tomber un orage de 50 mm au jour 100 et comparez la pointe de crue du bassin boisé à celle du bassin pavé ; pompez fort la nappe et voyez le débit de base mourir ; coupez l’humidité advectée et les réservoirs se vident ensemble. Neuf réservoirs plus les exportations moins les importations gardent la somme initiale à l’arrondi près, et la décrue de la nappe est l’exponentielle fermée.

Ce qui se passe
Le cycle de l’eau d’un bassin en neuf réservoirs concentrés en millimètres sur toute la surface : l’air au-dessus, la canopée de la forêt, la surface et le sol de la forêt, la surface et le sol de la ville, la nappe, la rivière et le lac. L’humidité advectée et un orage alimentent l’air, la pluie tombe au prorata des surfaces, la canopée en intercepte un quart et l’évapore, les surfaces infiltrent à la vitesse de la forêt ou de la ville et ruissellent à leur propre rythme, les sols évaporent et percolent, la nappe alimente la rivière et est pompée, le lac reçoit la moitié du ruissellement forestier et évapore, la rivière se vide à la mer, l’air exporte de la vapeur. Pas explicites d’un centième de jour avec les sorties d’un réservoir réduites ensemble pour qu’aucun ne devienne négatif ; une seconde d’expérience vaut un jour.
À essayer
Regardez les niveaux s’établir, puis l’orage du jour 100 faire monter la rivière ; pavez le bassin et comparez la pointe ; pompez la nappe et voyez le débit de base mourir ; vérifiez la décrue contre l’exponentielle.
Le modèle et ses limites
Réservoirs concentrés linéaires et saturants, un seul orage, pas de neige, d’évaporation saisonnière ni de délais de propagation au-delà des constantes de temps. Un bassin pédagogique, pas un modèle hydrologique.
Équations du module
Équations dS_i/dt = Σ F_in − Σ F_out, R = S_atm/τ_rain I = K_inf·f·S_surf/(S_surf + S_½)·(1 − S_soil/S_max) Q_base = S_gw/τ_gw, S_gw(t) = S_gw(0)·e^(−t/τ_gw) (no recharge) Σ S_i + E_out − I_in = const


