Improve the Accuracy in Numerical Modeling of Suspended Sediment Concentrations in the Hangzhou Bay by Assimilating Remote Sensing Data Utilizing Combined Techniques of Adjoint Data Assimilation and the Penalty Function Method

Suspended sediment dynamics play an important role in controlling nearshore and estuarine geomorphology and the associated ecological environments. Modeling the transport of suspended sediment is a complicated and challenging research topic. The goal of this study is to improve the accuracy of modeling the suspended sediment concentrations (SSCs) with newly developed techniques. Based on a three-dimensional suspended cohesive sediment transport model, the transport of suspended sediment and SSCs are simulated by assimilating SSCs retrieved from the Geostationary Ocean Color Imager (GOCI) with the adjoint data assimilation in the Hangzhou Bay, a typical strong tidal estuary along the coast of the East China Sea. To improve the effect of the data assimilation, the penalty function method, in which the reasonable constraints of the estimated model parameters are added to the cost function as penalty terms, will be introduced for the first time into the adjoint data assimilation in the SSCs modeling. In twin experiments, the prescribed spatially varying settling velocity is estimated by assimilating the synthetic SSC observations, and the results show that the penalty function method can further improve the effect of data assimilation and parameter estimation, regardless of synthetic SSC observations being contaminated by random artificial errors. In practical experiments, the spatially varying settling velocity is firstly estimated by assimilating the actual GOCI-retrieved SSCs. The results demonstrate that the simulated results can be improved by the adjoint data assimilation, and the penalty function method can additionally reduce the mean absolute error (MAE) between the independent check observations and the corresponding simulated SSCs from 1.44 &times: 10&minus:1 kg/m3 to 1.30 &times: 10&minus:1 kg/m3. To pursue greater simulation accuracy, the spatially varying settling velocity, resuspension rate, critical shear stress and initial conditions are simultaneously estimated by assimilating the actual GOCI-retrieved SSCs to simulate the SSCs in the Hangzhou Bay. When the adjoint data assimilation and the penalty function method are simultaneously used, the MAE between the independent check observations and the corresponding simulated SSCs is just 9.90 &times: 10&minus:2 kg/m3, which is substantially less than that when only the settling velocity is estimated. The MAE is also considerably less than that when the four model parameters are estimated to be without using the penalty function method. This study indicates that the adjoint data assimilation can effectively improve the SSC simulation accuracy, and the penalty function method can limit the variation range of the estimated model parameters to further improve the effect of data assimilation and parameter estimation.