An integrated geochemical and hydrodynamic model for tidal coastal environments

In this paper, the design, calibration and application of an integrated geochemical–hydrodynamic model are described. The model comprises three parts: A hydrodynamic submodel that was adopted from a depth-averaged, semi-implicit hydrodynamic model, a geochemical submodel based on equilibrium partitioning of chemicals between aqueous and particulate phases, and a particle dynamic submodel that simulates resuspension, transport and settling of suspended particulate matter (SPM). The integrated model was implemented in San Diego Bay (SDB), a heavily urbanized, semi-closed mesotidal embayment. A series of model calibrations were carried out based on observations on salinity, polychlorinated biphenyls (PCBs) and SPM. Salinity calibrations indicated that only 15% of precipitation in the drainage area of SDB could reach the bay, presumably due to the dams on the tributary rivers. Steady-state calibrations of PCBs based on fixed concentrations at known ‘hot spots’ have reproduced observed PCB concentrations in both the dissolved and particulate phases. SPM calibrations showed that shipping-induced resuspension produce more SPM than natural processes. Based on the calibrated model, the annual transport of PCBs out of SDB was estimated to be 3.85 kg (3.5 kg and 0.35 kg in the dissolved and particulate phases, respectively), much higher than the previous estimates based on steady-state assumptions. It was also found out that only a small portion of the fine sediment exported from SDB was derived from riverine input. This model can be applied to the studies of the transport and fate of other chemical species. It can be transplanted to other coastal areas as well. The integrated model represents a novel framework in which geochemical processes in coastal environments can be investigated on a truly dynamic basis.