Computed and Observed Currents, Elevations, and Salinity in a Branching Estuary

Abstract

A one-dimensional, hydrodynamical model of the Tamar Estuary shows good agreement with measured tidal elevations and currents. Computed currents are used to drive a one-dimensional moving-element model of the salt balance. The moving-element model overcomes the numerical difficulties associated with strong tidal advection. Axial distributions of salinity at high water, computed using the moving-element model, compare well with measurements. The modelled and observed high water salinity distributions in this macrotidal estuary show little dependence on tidal range. The major variability in salinity is due to runoff. This strong and rapid dependence on runoff is a consequence of short residence (or flushing) times. Typically, residence times are less than one day throughout the year in the upper 10 km of estuary. The residence times maximize in summer, reaching 14 d for the whole estuary. During high runoff winter periods residence times are less than 5 d. Mixing coefficients for the moving-element salinity model are deduced from salinity measurements. Dispersion coefficients at fixed locations along the estuary are deduced from solutions of the salinity model. The spatially-averaged coefficients at mean spring and neap tides are 180 and 240 m2 s−1, respectively, for average runoff. Therefore, spring-neap variations in dispersion are fairly small and show a negative correlation with tidal range. The spatially-averaged dispersion coefficients at mean tides vary from 150 to 300 m2 s−1 for typical summer and winter runoff, respectively. The increase in dispersion with runoff and the decrease with tidal range implies that buoyancy-driven currents generate an important component of the shear dispersion in this estuary.