A model for predicting the temporal evolution of dissolved oxygen concentration in shallow estuaries

The design of sewage discharge systems in estuaries needs to consider the dissolved oxygen concentration among other water quality indicators. Due to the great number of factors affecting the dissolved oxygen, the prediction of the temporal evolution of this element requires the use of mathematical tools. In the case of shallow estuaries with extensive intertidal zones, the complexity of this task increases since the water domain varies continuously. This work describes a numerical model which solves the vertical integrated transport equation including the effect of extensive tidal flats. The model makes use of the tidal velocities computed by a previous run of a hydrodynamic model. This procedure allows the dissolved oxygen modelling to be performed using a larger time step than that used for the velocity field calculation. The numerical scheme developed for the model guarantees that even in the drying and wetting of intertidal areas the water column dissolved oxygen concentration is not modified by numerical errors. The representations of the processes included in the model (advection, dispersion and reaction) have been validated successfully in several theoretical cases. An application to the Urdaibai Estuary, a shallow estuary in the Basque Country, Northern Spain, is also presented.