Earth
Marées dans un estuaire
Regardez la marée du soir remonter un estuaire en entonnoir devant le phare : le chenal se rétrécit vers la rivière, si bien que l’onde monte en avançant, que le marnage en tête dépasse de loin celui de l’embouchure, que le flot est plus court que le jusant et que l’eau de mer salée pousse des dizaines de kilomètres dans les terres contre la rivière. Augmentez le débit de la rivière et le sel est repoussé vers la mer ; accentuez l’entonnoir ou redressez-le ; fermez la tête et le chenal résonne en quart d’onde, la marée en tête a / cos(kL) de celle de la mer. Une eau calme fermée reste exactement calme, un estuaire salé sans rivière exactement salé, et l’eau et le sel s’équilibrent avec ce qui a franchi l’embouchure et la tête.

Ce qui se passe
Un chenal de 60 km et 8 m de profondeur dont la largeur décroît depuis 3 km comme exp(−x/25 km), en 120 cellules sur une grille décalée : la surface au centre des cellules, la vitesse aux faces. Des pas avant-arrière de 20 s mettent à jour la surface à partir des débits aux faces, puis la vitesse à partir de la nouvelle pente, avec advection décentrée amont et frottement de Manning implicite. La marée M2 de 1,5 m de la mer entre à l’embouchure sur deux périodes ; la rivière apporte 100 m³/s en tête. Le sel suit les mêmes débits en décentré amont avec une dispersion de marée de 150 m²/s, si bien qu’un estuaire salé sans rivière reste salé et qu’une eau fermée calme reste calme. Un chenal prismatique fermé vérifie le quart d’onde a/cos(kL).
À essayer
Regardez le marnage croître vers le haut de l’entonnoir et le flot durer moins que le jusant ; montez la rivière et chassez le sel ; redressez ou accentuez l’entonnoir ; fermez la tête et comparez au quart d’onde.
Le modèle et ses limites
Un chenal unidimensionnel bien mélangé de profondeur uniforme, sans assèchement des estrans, une salinité sans effet de densité, une seule composante de marée. Les marnages appartiennent à cet estuaire idéalisé.
Équations du module
Équations b·∂η/∂t + ∂(b·H·u)/∂x = 0, b = b0·e^(−x/L_b) ∂u/∂t + u·∂u/∂x = −g·∂η/∂x − g·n²·|u|·u/H^(4/3) ∂(b·H·S)/∂t + ∂(b·H·u·S − b·H·K·∂S/∂x)/∂x = 0 closed prismatic channel: a_head = a/cos(kL), k = ω/√(g·h)


