Abstract:The Xiangshangang Bay is a typical narrow semi-enclosed bay in the middle of the Zhejiang Province,located at the south edge of the Changjiang River Delta.Along its axis from southwest to northeast,the Xiangshangang Bay extends about 60 km with a water area of 563 km2,in which 171 km2 is tidal flat(above the base level of theoretical depth).The tidal range in the bay is large.It increases gradually from the entrance inwards with a mean value of above 3 m.At the innermost part it reaches a maximum of 6 m with a mean range of 4 m or so.The water area of the mudflats changes markedly with tide.When the water level is high,the whole bay is merged by water,and when the water level falls,the shallow areas are exposed.Therefore,the influence of drying and flooding of the mudflats should be considered when a numerical model is applied for such a region with wide mudflats and large tidal range.Various drying and of flooding algorithms have been developed very well for 2D models.However,if those algorithms used in 2D models are directly introduced into 3D models,the numerical oscillations or instability could occur,and thus it is necessary to find an easy and feasible mobile boundary technique for 3D models.The original POM(Princeton Ocean Model) with a fixed boundary cannot deal with the overbank flow.Base on the work of Burchard et al in 2004,the fixed boundary in original POM was replaced with a naturally variable boundary following a strategy to let the physics do the best to meet the requirement of flow reduction.Some simplifications of the model physics will be applied for very shallow water in order to describe the physical process in a proper way.In this method,a dimensionless damping factor,which varies linearly between unity(deep water) and zero(very shallow water),was introduced to the horizontal advection,horizontal diffusion,Coriolis force and baroclinicity terms of the momentum equations.When the whole water depth is larger than the set critical water depth,usual momentum equations are resulted,and when the whole water depth trends approaching the minimum value,the momentum equations are simplified to be the balance among the tendency,friction and external pressure gradient.The friction dominates in such a way that for water depths approaching zero,the current velocities converge to zero as well,and as a result,the unphysical negative water depths can be avoided.In addition to the simplification of the momentum equations at very shallow water depths,a correction to the barotropic pressure gradient should be applied where the sea level cell is below the sea bed at the adjacent grid cell,so as to avoid unreasonable pressure gradient.Furthermore,the expression of the frictional coefficient at the bottom was modified so that the model becomes more stable to deal with drying and flooding of mudflats.The improved POM was applied to simulate 3D tidal current field of the Xiangshangang Bay.The computed results matched very well the observed data from several survey stations including 2 tidal level stations,1 mooring site and 2 anchor stations.The model has reasonably reproduced the 3D structure and variation of the tidal current in the whole Xiangshangang Bay.Horizontal and vertical distributions of the tidal velocity in the innermost part of the bay at the beginning of the ebb tide,the peak ebb tidal currents,the beginning of the flood tide and the peak flood tidal currents show that the model is able to present the correct processes of drying and flooding of mudflats.It is hence indicated that the modification of the origin POM is very reasonable and feasible,and will enable POM to be used more extensively.

Numerical Study on Tidally Induced Cross-Frontal Mean Circulation

Numerical Study of Turbidity Current over a Three-Dimensional Seafloor

Numerical Simulation of Three-Dimensional Characteristics of Tidal Current and Residual Current in Yangshan Harbor

Numerical Study of the Tides and Residual Currents in the Qiongzhou Strait

Numerical Modeling of Lock-Exchange Gravity/turbidity Currents by a High-Order Upwinding Combined Compact Difference Scheme

A Numerical Study on the Effect of Tidal Flat’s Slope on Tidal Dynamics in the Xiangshan Bay, China

Numerical simulation study on the residual currents in the Xiangshan Bay

Numerical Study on the Velocity Structure Around Tidal Fronts in the Yellow Sea

Momentum Flux Balance At The Air-Sea Interface

Numerical Solution of Unsteady Ekman Equation Modified by Wave and Linear Friction Term

High-resolution Observations of Longitudinal Fronts in a Well-Mixed Tidal Channel

Numerical Modelling of Hydrodynamics and Sand Transport in the Tide-Dominated Coastal-to-estuarine Region

Effects of Wind Input and Wave Dissipation Formulations on the Steady Ekman Current Solution

A Numerical Study of Residual Flow Induced by Eddy Viscosity-Shear Covariance in a Tidally Energetic Estuary

The Effects of Random Surface Waves on the Steady Ekman Current Solutions

The Turbulent Structure Of The Marine Atmospheric Boundary Layer During And Before A Cold Front

Subtidal water flux through a multiple-inlet system: Observations before and during a cold front event and numerical experiments

Tide-Induced Suspended Sediment Transport: Depth-Averaged Concentrations and Horizontal Residual Fluxes

Observations of a Tidal Intrusion Front in a Tidal Channel

Modeling tides and their influence on the circulation in Prince William Sound, Alaska

3D Variable-Boundary Numerical Modeling of Tidal Current in the Xiangshangang Bay