A three‐dimensional numerical model of suspended sediment transport in Howe sound, British Columbia

Abstract

Numerical models of the hydrodynamics of Howe Sound and of the transport and deposition of suspended sediment from the Squamish River entering Howe Sound are described. Two different versions of a three‐dimensional hydrodynamic model of Howe Sound have been developed. One is a fine‐grid model employing a 390‐m grid size, the other a coarse‐grid model with a 1.95‐km grid size. In the present implementation, both models have 9 layers in the vertical, and calculate the three components of the velocity field as well as the time‐varying density field. The models are forced at their open boundaries by a larger three‐dimensional model of the Strait of Georgia, by surface wind stress and by the Squamish River flow. The coarse‐grid model is computationally less demanding and is used for long‐term simulations, whereas the fine‐grid model is used for detailed simulations for shorter time periods.

The sediment model derives its velocity field from the coarse grid three‐dimensional hydrodynamic model. Sediment transport and deposition are simulated using a Monte Carlo procedure, in which the transport, diffusion and settling of suspended sediment are reproduced by tracking the fate of a large number of pseudo‐particles, each having dynamics corresponding to the sediment grain size it represents. The effects of diffusion are simulated by prescribing a random velocity for each particle, in addition to the velocity derived from the three‐dimensional model and the intrinsic settling rate. Temporal variations in supply are simulated by corresponding variations in the rate of installation of pseudo‐particles. The grain size distribution of the source sediment is simulated by randomly selecting, for each pseudo‐particle,.a sinking velocity from a range of sinking velocities whose distribution reproduces the distribution of sinking velocities in the supply. The model was run for a 69‐day period, sufficient time for virtually all the sediment in the summer Squamish River freshet to settle to the floor of Howe Sound. The resulting annual rate of deposition is similar to calculations reported by Hicken (1989).

Résumé

On décrit les modéles numèriques de l'hydrodynamique de la baie Howe ainsi que du transport et du dépôt d'alluvions à l'embouchure de la rivière Squamish. On a développé deux versions différentes du modèle hydrodynamique à trois dimensions (3‐D) de la baie Howe : un modèle à maille fine utilisant une maille de 390 m et un modèle plus général de 1,95 km. Actuellement, les deux modèles utilisent 9 couches dans la verticale et calculent les trois composantes du champ de vélocité ainsi que le champ de densité temporalisé. Ils sont forcés à leurs limites ouvertes par un plus grand modèle 3‐D du détroit de Géorgie, par la force d'entraînement du vent et par l'écoulement de la Squamish. Le modèle général, utilisé pour des simulations à long terme, demande moins de travail alors que celui à maille fine sert à des simulations plus détaillées et sur de plus courtes périodes.

Le modèle des alluvions dérive son champ de vélocité du modèle hydrodynamique général 3‐D. Le transport et le dépôt des alluvions sont simulés en utilisant la méthode de Monte Carlo dans laquelle le transport, la diffusion et la décantation des alluvions suspendues sont reproduits en suivant le devenir d'un grand nombre de pseudo‐particules, chacune ayant une dynamique correspondant à la grosseur des grains qu ‘elle représente. On simule les effets de la diffusion en allouant une vélocité complexe à chaque particule en plus de la vélocité provenant du modèle 3‐D et du taux intrinsèque de décantation. Les variations temporelles de l'approvisionnement sont simulées par des variations correspondantes dans le taux d'installation des pseudo‐particules. La distribution de la grosseur des grains des sédiments à la source est simulée en choisissant au hazard, pour chacune des pseudoparticules, un temps de chute dans une fourchette dont la distribution reproduit la distribution des temps de chute dans l'approvisionnement. Le modèle a été exécuté sur une période de 69 jours, temps suffisant pour presque toutes les alluvions de la crue estivale de se déposer au fond de la baie. Les taux annuels de déposition sont semblables aux calculs de Hicken (1989).