A Barotropic Tide Model for Global Ocean Based on Rotated Spherical Longitude-Latitude Grids

Ocean modeling and simulation are important for understanding the dynamic processes in the geophysical system, and the simulation of tidal dynamics is of great significance for understanding the dynamic evolution of the ocean. However, there are some problems in existing simulations, including lack of specific standards to produce a desirable discrete spherical mesh for global ocean modelling. Many global ocean numerical models based on conventional longitude-latitude (LL) coordinates suffer from the “pole problem” in regions adjacent to the North Pole due to the convergence of meridians, which seriously hinders global ocean simulations. In this paper, a new longitude-latitude spherical grid coupled with rotated coordinate mapping is proposed to overcome the problem. In the design of the numerical model, for spatial approximation, the finite volume method on staggered C grid is proposed to solve the two-dimensional tidal wave equations for the global ocean. For temporal integration, the third-order Adams-Bashforth method is used to explicitly extrapolate the value on the next time interval half layer, and then the fourth-order implicit Adams-Moulton method is used to correct the water level. Finally, the constructed model is used to simulate the dynamics of two-dimensional tidal waves in the global ocean, and the co-tidal maps of two major diurnal tide and semidiurnal tide components are shown. The results demonstrate that the proposed model can support the simulation of tidal dynamics in the global ocean, especially for the Arctic Ocean.