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.

The Role of Sea Water Viscidity in Modeling the Vertical Movement of Internal Tides

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

High-fidelity Simulation of Regular Waves Based on Multi-Moment Finite Volume Formulation and THINC Method

A CIP-based Numerical Simulation of Wave Interaction with a Fluid-Filled Membrane Submerged Breakwater

Modeling Vertical Heat Transport in the Wave Affected Surface Layer of the Ocean

Numerical simulation study on the residual currents in the Xiangshan Bay

Wave-Modified Ekman Current Solutions for the Time-Dependent Vertical Eddy Viscosity

Numerical Study of the Tidal Refraction Effect on the Open Coast in Southeastern Zhejiang

An Intercomparison and Validation of a Range of Turbulence Closure Schemes Used in Three Dimensional Tidal Models

Modelling the effect of sea ice in an ocean tide model

Formulation of a three-dimensional shelf edge model and its application to internal tide generation

3D Modal Solution for Tidally Induced Lagrangian Residual Velocity with Variations in Eddy Viscosity and Bathymetry in a Narrow Model Bay

INTERNAL LEE WAVES AND TURBULENCE MIXING OVER AN ISOLATED SEAMOUNT: RESULTS FROM TURBULENCE ENERGY MODELS

Internal Tide Generation and Dissipation by Small Periodic Topography in Deep Ocean

A Three‐dimensional Model of Internal Tides on the Malin‐Hebrides Shelf and Shelf Edge

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

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

The Influence of Topographic Features and Density Variations Upon the Internal Tides in Shelf Edge Regions

Ocean Turbulence. Part I: One-Point Closure Model—Momentum and Heat Vertical Diffusivities

A Model Study of Internal Tides in Coastal Frontal Zone