Tide simulation in a global eddy-resolving ocean model

Tide simulation in a global eddy-resolving ocean model

PDF

Zhiwei Tian¹, Caixia Wang¹, Zipeng Yu²^,^*, Hailong Liu²^,³^,^*, Pengfei Lin²^,⁴, Zhuhua Li⁵

Acta Oceanologica Sinica | 2024, 43(9) : 1 - 10

Less

Acta Oceanologica Sinica | 2024, 43(9): 1-10

• Physical Oceanography, Marine Meteorology and Marine Physics •

Tide simulation in a global eddy-resolving ocean model

Full

Zhiwei Tian¹, Caixia Wang¹, Zipeng Yu²^,^*, Hailong Liu²^,³^,^*, Pengfei Lin²^,⁴, Zhuhua Li⁵

Affiliations

¹ Physical Oceanography Laboratory, Ocean University of China, Qingdao 266100, China

² State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

³ Laoshan Laboratory, Qingdao 266237, China

⁴ College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China

⁵ Hohai University, Nanjing 210098, China

Published: 2024-09-25 doi: 10.1007/s13131-024-2352-5

Outline

Abstract

Less

The tide plays a pivotal role in the ocean, affecting the global ocean circulation and supplying the bulk of the energy for the global meridional overturning circulation. To further investigate internal tides and their impacts on circulation, it is imperative to incorporate tidal forcing into the eddy-resolving global ocean circulation model. In this study, we successfully incorporated explicit tides (eight major constituents) into a global eddy-resolving general ocean circulation model and evaluated its tidal simulation ability. We obtained harmonic constants by analyzing sea surface height through tidal harmonic analysis and compared them with the analysis data Topex Poseidon Cross-Overs v9 (TPXO9), the open ocean tide dataset from 102 open-ocean tide observations, and tide gauge stations from World Ocean Circulation Experiment. The results demonstrated that the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics (LASG/IAP) Climate System Ocean Model 3.0 (LICOM3.0) effectively simulated tides, with errors predominantly occurring in nearshore regions. The tidal amplitude simulated in LICOM3.0 was greater than that of TPXO9, and these high-amplitude areas exhibited greater errors. The amplitude error of the M₂ constituent was larger, while the phase error of the K₁ constituent was more significant. Furthermore, we further compared our results with those from other models.

Key words

tide / eddy resolution / ocean general circulation models / harmonic analysis / LICOM3.0

Cite this Article

Zhiwei Tian, Caixia Wang, Zipeng Yu, Hailong Liu, Pengfei Lin, Zhuhua Li. Tide simulation in a global eddy-resolving ocean model[J]. Acta Oceanologica Sinica, 2024 , 43 (9) : 1 -10 . DOI: 10.1007/s13131-024-2352-5

Full Text

Less

1 Introduction

Less

Diapycnal mixing is of great significance for global ocean stratification and driving meridional overturning circulation (Marshall and Speer, 2012), which can significantly affect local and global climate change through the diapycnal transport of matter and heat (Kostov et al., 2014). This process can have a profound impact on regional and global climate systems. The wind-driven upwelling of the Southern Ocean and the diapycnal mixing of the Pacific and Indian Oceans are considered two crucial driving processes (Talley, 2013), while the tidal process in the open ocean is the main source of mixing energy (Munk and Wunsch, 1998). Munk and Wunsch (1998) analyzed the tidal energy budget based on the altimeter results of Egbert (1997). The results showed that the Sun-Moon-Earth system delivered a total of 3.5 TW of energy to the ocean, of which 2.6 TW dissipated in the shallow water in the form of barotropic tides, and the remaining 0.9 TW was dissipated in the deep ocean. Deep ocean dissipation could be divided into local dissipation and remote dissipation. Nevertheless, local and remote dissipation could provide considerable energy for maintaining stratification and driving thermohaline circulation in the deep sea. Egbert and Ray (2000), Jayne and St. Laurent (2001), and Niwa and Hibiya (2014) estimated that the magnitude of dissipation due to tides in the deep ocean was approximately 1 TW. Therefore, studying tides and their dissipation is crucial for understanding ocean circulation.

Due to limitations in computing power and theoretical understanding, early tidal and global ocean circulation simulations were conducted separately. Tidal numerical models are mainly based on a two-dimensional numerical model of the tidal wave dynamic equation. In contrast, ocean general circulation models (OGCMs) are based on simplified Navier‒Stokes equations without considering tidal processes. However, with advancements in the development of satellite altimeter technology, two-dimensional tidal models have undergone significant improvements, resulting in relatively mature prediction and simulation capabilities (Stammer et al., 2014). Tidal data that closely match observations can be obtained by integrating two-dimensional tidal models with data assimilation (Shum et al., 1997). The two-dimensional tidal model can be refined by incorporating altimeter data, resulting in an increasingly precise tidal map. Nevertheless, several issues remain unresolved, including the correlation between tides and climate change, the interplay between tides and large-scale ocean circulation, the mechanism behind internal tidal dissipation in the open ocean, and the prediction of internal tides. Addressing these questions necessitates the concurrent simulation of global ocean circulation and tides (Arbic et al., 2018).

There are usually two methods for incorporating tides into OGCMs. The first approach is implicit, where the impacts of tides on the ocean, such as mixing induced by internal tides, are parameterized (St. Laurent et al., 2002). This parameterization enhances mixing in the open ocean and makes the coefficient of diapycnal mixing not uniform in the horizontal direction. Simmons et al. (2004) added the parameterization scheme of St. Laurent et al. (2002) to a low resolution OGCM. The scholars found that adding the parameterization of tidal mixing reduces the deviation of ocean simulation. Yu et al. (2017) showed that this parameterization significantly enhances the Atlantic meridional overturning circulation (AMOC). Although several works have shown that the parameterization of tidal mixing has effects on improving ocean circulation simulation (Jayne, 2009; Melet et al., 2013; Saenko and Merryfield, 2005; Simmons et al., 2004; Song et al., 2023; Yu et al., 2017), uncertainty remains a problem in the parameterization scheme. Changes in some coefficients in this scheme can significantly impact the results of ocean model simulations (Jayne, 2009). The second approach is explicit, where the tidal potential is directly introduced into the sea surface height formula of the ocean model. Thomas et al. (2001) and Schiller (2004) incorporated the tide explicitly into an OGCM. Nevertheless, due to some numerical problems with the model, Schiller (2004) could only strictly limit the tidal potential in the Indonesian sea. Schiller and Fiedler (2007) successfully incorporated the tidal potential of the eight constituents into an OGCM, and the results showed that the inclusion of tides enhanced the simulation of temperature and salinity in the Indonesian throughflow. Müller et al. (2010) discussed the impact of incorporating tidal potential into a climate model and demonstrated that the temperature simulation in the North Atlantic was improved. It is worth mentioning that these studies were conducted with low resolution OGCMs (Müller et al., 2010; Schiller, 2004; Schiller and Fiedler, 2007; Thomas et al., 2001).

Arbic et al. (2010) were the earliest researchers to incorporate explicit tides into an eddy-resolving global OGCM. The scholars incorporated the tidal gravitational potential of eight major tidal constituents in the hybrid coordinated ocean model (HYCOM) with realistic ocean stratification. The parameterized topographic wave drag was considered to simulate the tides accurately. The simulated barotropic tides were found to be as accurate as previous tidal models without data assimilation (Stammer et al., 2014). Since then, global eddy-resolving OGCMs that incorporate atmospheric and tidal forcing have emerged successively, such as STORMTIDE (Müller et al., 2012), Massachusetts Institute of Technology General Circulation Model (MITgcm) (Rocha et al., 2016a, b), and Icosahedral Nonhydrostatic Weather and Climate Model (ICON-O) (Logemann et al., 2021; von Storch et al., 2023). In addition, some researchers have considered tides using boundary conditions in regional high resolution ocean models (Peng et al., 2021; Siyanbola et al., 2023; Wang et al., 2016). In future studies of eddy-resolving ocean models, concurrent atmospheric and tidal forcings (wind plus tide) should be increasingly utilized (Arbic et al., 2018).

Global eddy-resolving OGCMs prioritize high-frequency movements in the ocean, the reciprocal impacts of internal tides and ocean circulation, and the prediction of internal tides (Arbic et al., 2012, 2018; Li and von Storch, 2020; Müller et al., 2012, 2014, 2015; Ansong et al., 2015, 2017; Buijsman et al., 2016; Savage et al., 2017a, b). However, evaluating the simulation of barotropic tides in a global eddy-resolving OGCM is essential before studying related physical processes. Arbic et al. (2010) reported an error of 8.26 cm for the M₂ constituent, which they found through a comparison of the 102 open-ocean tide observations made by Shum et al. (1997); this value is equivalent to the tides simulated by Arbic et al. (2004). Müller et al. (2014) conducted the same evaluation, and the error of the M₂ constituent reached 8.19 cm. Additionally, seasonal variability in the amplitude of the M₂ constituent was found. von Storch et al. (2023) found slightly larger errors in simulated barotropic tides than Arbic et al. (2010) and Müller et al. (2014). The researchers noted that improving the horizontal resolution of OGCM did not necessarily reduce open-ocean tide simulation errors.

We have performed some related studies with the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics (LASG/IAP) Climate System Ocean Model (LICOM). Yu et al. (2016) incorporated tidal forcing into LICOM2.0 with a 1° horizontal resolution, analyzed its impact on ocean circulation simulation, and investigated its effect on climate via coupled models (Yu et al., 2020). Jin et al. (2022) adopted a new formulation of an explicit tidal scheme and successfully tested it in LICOM2.0 with the same resolution. Furthermore, the global eddy-resolving version of LICOM3.0 has undergone preliminary evaluation without tidal forcing (Li et al., 2020). In addition, incorporation and long-term tests of tidal forcing have been recently completed. In the present study, we aim to incorporate tidal forcing into the eddy-resolving model LICOM3.0 and systematically evaluate the simulated barotropic tides, laying a foundation for future research. Section 2 details the model used, the model configuration, the tidal forcing incorporation method, and the error analysis method. Evaluations of the errors of barotropic tides are described in Sections 3 and 4, providing the conclusion of this article and the future outlook.

2 Model, experiment, and methods

Less

2.1 Astronomical tidal forcing

Based on the theory of equilibrium tides (Pugh, 1996), Yu et al. (2016) incorporated the tidal forcing of eight major constituents into LICOM2.0, including four major diurnal constituents (K₁, P₁, Q₁, and O₁) and four major semidiurnal constituents (M₂, S₂, N₂, and K₂). In this article, based on Yu et al. (2016) and Arbic et al. (2018), we incorporated the tidal forcing of eight constituents into LICOM3.0. Ocean tides were generated by the gravitational forces of the moon and the sun on seawater. The diurnal tide

$ {\eta }_{\mathrm{EQ}1} $

and the semidiurnal tide

$ {\eta }_{\mathrm{EQ}2} $

could be calculated by Eqs (1) and (2), respectively, as follows:

(1)

$ \begin{split}\eta_{\mathrm{EQ}1}=\, & Af\left(t_{\mathrm{ref}}\right)\left(1+k_2-h_2\right)\ \sin(2\phi)\times \\ &\cos\left(\omega(t-t_{\mathrm{ref}})+V_0(t_{\mathrm{ref}})+\mu(t_{\mathrm{ref}})+\mathit{\lambda}\right),\end{split} $

(2)

$ \begin{split}\eta_{\mathrm{EQ}2}=\, & Af\left(t_{\mathrm{ref}}\right)\left(1+k_2-h_2\right)\ \cos^2\phi\times \\ & \cos\left(\omega(t-t_{\mathrm{ref}})+V_0(t_{\mathrm{ref}})+\mu(t_{\mathrm{ref}})+2\mathit{\lambda}\right),\end{split} $

where the equilibrium tide

$ {\eta }_{\mathrm{E}\mathrm{Q}} $

is the sum of

$ {\eta }_{\mathrm{E}\mathrm{Q}1} $

and

$ {\eta }_{\mathrm{E}\mathrm{Q}2} $

$ \phi $

is the latitude,

$ \lambda $

is the longitude,

$ t $

is the time,

$ {t}_{\mathrm{r}\mathrm{e}\mathrm{f}} $

is the reference time of the model,

$ {V}_{0}({t}_{\mathrm{r}\mathrm{e}\mathrm{f}}) $

is the initial phase of the constituents at the reference time,

$ A $

and

$ \omega $

are the amplitude and frequency of the equilibrium tide for each constituent, respectively,

$ f({t}_{\mathrm{r}\mathrm{e}\mathrm{f}}) $

and

$ \mu ({t}_{\mathrm{r}\mathrm{e}\mathrm{f}}) $

are the nodal factor and nodal angle, respectively, representing the modulations of the amplitude and delay angle of the tides over 18.61 years (Pugh, 1996), and

$ {h}_{2} $

and

$ {k}_{2} $

are Love numbers, which explain the deformation of the solid Earth caused by astronomical forcing and the influence of the redistribution of Earth’s mass on the gravitational potential, respectively. In this study,

$ 1+{k}_{2}-{h}_{2} $

is taken as a constant value of 0.7.

The nodal factor

$ f\left({t}_{\mathrm{r}\mathrm{e}\mathrm{f}}\right) $

and the phase-dependent values

$ {V}_{0}({t}_{\mathrm{r}\mathrm{e}\mathrm{f}})+\mu ({t}_{\mathrm{r}\mathrm{e}\mathrm{f}}) $

are only related to astronomical variables. The formula for

$ {V}_{0}({t}_{\mathrm{r}\mathrm{e}\mathrm{f}}) $

(3)

$ V_0=n_1\tau_0+n_2s_0+n_3h_0+n_4p_0-n_5N_0+n_6p_0'+n_0\times90^{\circ}, $

(4)

$ {\tau }_{0}=180^{\circ}-{s}_{0}+{h}_{0}, $

where n ₀−n ₆ are Doodson numbers, found in Appendix Table 4 in Zuo et al. (2018).

$ {\tau }_{0} $

$ {s}_{0} $

$ {h}_{0} $

$ {p}_{0} $

$ {N}_{0} $

, and

$ {p}_{0}' $

are mean lunar time, mean longitude of moon, mean longitude of sun, mean longitude of lunar perigee, mean longitude of lunar ascending node, and mean longitude of solar perigee, respectively. The formula for

$ f\left({t}_{\mathrm{ref}}\right) $

and

$ \mu ({t}_{\mathrm{ref}}) $

are shown in Table 1. The astronomical parameters are calculated by

(5)

$ {s}_{0}=277.02^{\circ}+481\,267.89^{\circ}T+0.001\,1^{\circ}{T}^{\,2}, $

(6)

$ {h}_{0}=280.19^{\circ}+36\,000.77^{\circ}T+0.000\,3^{\circ}{T}^{\,2}, $

(7)

$ {p}_{0}=334.39^{\circ}+4\,069.04^{\circ}T-0.010\,3^{\circ}{T}^{\,2}, $

(8)

$ {N}_{0}=259.16^{\circ}-1\,934.14^{\circ}T+0.002\,1^{\circ}{T}^{\,2}, $

(9)

$ {p}_{0}'=281.22^{\circ}+1.72^{\circ}T+0.000\,5^{\circ}{T}^{\,2}, $

where T is Julian century since January 1, 1900, T =

$ [365(Y\ - 1\,900)+(D-1)+i]/{36\,525} $

$ Y $

is the year of model time,

$ D $

is the day of the model time in year

$ Y $

, and

$ i $

is the number of leap years from 1900 to year

$ Y $

. In the subsequent analysis, we compared the nodal factor

$ f\left({t}_{\mathrm{r}\mathrm{e}\mathrm{f}}\right) $

and the phase-dependent

$ {V}_{0}({t}_{\mathrm{r}\mathrm{e}\mathrm{f}})+ \mu ({t}_{\mathrm{r}\mathrm{e}\mathrm{f}}) $

values of the theoretical and LICOM3.0 model results; moreover, the equilibrium tide

$ {\eta }_{\mathrm{E}\mathrm{Q}} $

is shown.

By incorporating the sea surface height caused by tides into the sea surface height

$ \eta $

of the barotropic equation in the model, the gradient scheme for the sea surface height became the following formula:

(10)

$ \nabla \left[\left(1-\alpha \right)\eta -{\eta }_{\mathrm{E}\mathrm{Q}}\right], $

where the coefficient

$ \alpha $

represents the self-attraction and loading term (Hendershott, 1972). In this article,

$ \alpha $

is taken as 0.948.

The tide can change with the relative position of the Sun, Moon, and Earth. The tide is the strongest when the three celestial bodies are in a straight line, where it is called a spring tide. When the angle between the straight line connecting the Sun and the Moon to the Earth is a right angle, the tide is relatively weak, and it is called a neap tide.

To ensure model stability, we added the parameterized topographic wave drag term of Jayne and St. Laurent (2001), while in the work of Yu et al. (2016), only the bottom friction term was implemented. The formula for the wave drag term used was as follows:

(11)

$ {\tau }_{0}=c{N}_{0{\mathrm{b}}}{U}_{0}\kappa {h}^{2}, $

where c is a coefficient and can be adjusted,

$ {N}_{0{\mathrm{b}}} $

represents the buoyancy frequency at the bottom boundary,

$ {U}_{0} $

denotes the barotropic velocity in the direction of the reference drag layer,

$ {N}_{0{\mathrm{b}}} $

and

$ {U}_{0} $

are computed in LICOM3.0 with model density and velocities,

$ \kappa $

denotes the topographic wavenumber and is set to 2π/10 km (Jayne and St. Laurent, 2001), and

$ {h}^{2} $

represents the topographic roughness and is computed using the (1/30)° Earth Topography v2 (ETOPO2; National Centers for Environmental Information, 2006) grid as the average squared depth anomaly over a 1° square for a plane fit. Arbic et al. (2010) and Buijsman et al. (2020) included a parameterized topographic wave drag term and found that the wave drag term had a certain impact on the accuracy and energy distribution of the simulated tides.

2.2 Model and experiment

The model used in this study was LICOM (LASG/IAP Climate System Ocean Model). This model is an OGCM developed by the National Key Laboratory for Numerical Simulation of Atmospheric Science and Earth Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, China. The version used in this paper was the eddy-resolving LICOM3.0 (Li et al., 2020). The horizontal grid of LICOM3.0 was a tripolar grid (Murray, 1996), which could effectively address the singularity issues of latitude and longitude grids. The eddy-resolving version had a resolution of (1/10)°, with 55 layers of

$ \eta $

coordinates vertically. The sea ice model used was Community Ice Code version 4 (CICE4).

To study the influence of wave drag coefficient c on the accuracy of barotropic tides simulation, c is set to 0.5, 1.0, and 1.5 in the wave drag scheme; the three experiments are called wd0.5, wd1, and wd1.5. The astronomical tidal forcing is added to these three experiments, starting from January 1, 2015, to March 31, 2015, using the results of the phase 2 of the Ocean Model Intercomparison Project (OMIP2) experiment (Li et al., 2020) as initial values. The tide has gone through dozens of cycles in three months, so the barotropic tides are approximately stable. The sea surface height of the three experiments is conducted for harmonic analysis of tides. The results can be seen in Table 2. As c increases, the amplitude and phase errors increase of the four major constituents, as well as the total error. Due to the high cost of conducting high-resolution global experiments with hourly output, c = 0.5 is considered the relatively optimal parameter so far. So, in the following experiment, c is set to 0.5.

Based on the OMIP2 experiment (Li et al., 2020) using JRA55-do (Tsujino et al., 2018), the astronomical tidal forcing in the form of tidal potential was further incorporated, and the wind plus tide experiment, as defined by Arbic et al. (2018), was conducted. The standard OMIP2 experiment covered the period from 1958 to 2018, while the tidal experiment was integrated starting from January 1, 2015, to December 31, 2016, using the results of the OMIP2 experiment as initial values. The data from 2016 were analyzed in this paper. For tidal harmonic analysis, the data output interval of the tidal experiment was one hour.

2.3 Harmonic analysis and error analysis

To evaluate the simulated tides, we used the products of Topex Poseidon Cross-Overs v9 (TPXO9; Egbert and Erofeeva, 2002), which is the latest product of the global barotropic ocean tidal model TPXO, in this paper. We interpolated these results to a grid with a horizontal resolution of (1/10)°. In addition, we used 102 open-ocean tide observations (hereafter referred to as st102) from Shum et al. (1997) and data from two tide gauge stations, Yap (9.508 3°N, 138.128 3°E) and Kodiak (57.731 7°N, 152.511 7°W), selected from World Ocean Circulation Experiment (WOCE; Ponchaut et al., 2001). Moreover, we used the Python version of UTIDE (Codiga, 2011) to perform the harmonic analysis of the sea surface height of the tide experiment and obtained the average amplitude and phase values of the eight constituents.

Based on the definition of Arbic et al. (2004), the square of the total error of tides

$ {d}^{2} $

between LICOM3.0 and the observation data at a certain grid point could be defined as the sum of the square of the amplitude error

$ {d}_{\mathrm{a}}^{2} $

and the square of the phase error

$ {d}_{\mathrm{p}}^{2} $

(12)

$ \begin{split}{d}^{2}=\,&{d}_{\mathrm{a}}^{2}+{d}_{\mathrm{p}}^{2}=\frac{1}{2}{\left({A}_{\mathrm{model}}-{A}_{\mathrm{obs}}\right)}^{2}+\\&{A}_{\text{model}}{A}_{\text{obs}}\left(1-\mathrm{cos}\left({\phi }_{\mathrm{model}}-{\phi }_{\mathrm{obs}}\right)\right), \end{split}$

where

$ {A}_{\text{model}} $

and

$ {A}_{\text{obs}} $

represent the amplitude of the model and observation of a certain constituent, respectively, and

$ {\phi }_{\text{model}} $

and

$ {\phi }_{\mathrm{o}\mathrm{b}\mathrm{s}} $

represent the phase of the model and observation of a certain constituent, respectively. Obviously, the corresponding square of the regional or grid point average errors are as follows:

(13)

$ {D}^{2}=\overline{\left({d}^{2}\right)},$

(14)

${D}_{\mathrm{a}}^{2}=\overline{\left({d}_{\mathrm{a}}^{2}\right)},$

(15)

${D}_{\mathrm{p}}^{2}=\overline{\left({d}_{\mathrm{p}}^{2}\right)},$

where

$ \overline{(\cdot )} $

indicates an average over grid points or area.

To evaluate the simulations of all constituents, the Root Square Sum (RSS) values of all constituent errors were considered. The percentage of sea surface height variance captured was defined as follows, where RSS² is the sum of the square of average error D ² of all constituents and S is the signal of tides:

(16)

$ {V}_{\mathrm{c}\mathrm{a}\mathrm{p}\mathrm{t}}=100\%\times \left(1-\frac{\mathrm{R}\mathrm{S}{\mathrm{S}}^{2}}{{S}^{2}}\right), $

(17)

$ {S}^{2}=\frac{\overline{A}_{\mathrm{obs}}^{2}}{2}, $

when the percentage of sea surface height variance captured (V _capt) was 100%, the model perfectly simulated the barotropic tides. Therefore, the higher this value was, the more accurate the model simulation of the tides.

3 Results

Less

3.1 Tidal forcing

First, we verified whether the phase-dependent

$ {V}_{0}+\mu $

and nodal factor f in LICOM3.0 are correct. We output these two items from the module of tidal forcing and compare them with the results of Fisheries and Oceans Canada (2017 , 2018) (DFO) and the tidal harmonic analysis program UTIDE. We recorded the results on January 1 and July 1, as shown in Table 3. We found few differences between phase-dependent

$ {V}_{0}+\mu $

and nodal factor f in the DFO and UTIDE. Excluding the Q₁ constituent, the phase-dependent

$ {V}_{0}+\mu $

simulated by LICOM3.0 is within 0.31° of DFO and UTIDE, and the difference in the Q₁ constituent is slightly larger at approximately 1.632°. The relative error of nodal factor f simulated by LICOM3.0 compared to DFO and UTIDE is less than 1.22%. Therefore, the phase-dependent

$ {V}_{0}+\mu $

and nodal factor f calculated by the tidal module in LICOM3.0 are basically correct when compared with the references, with no significant differences.

The patterns of the daily changes of the spring and neap tides calculated by LICOM3.0 are shown in Fig. 1, where four time points on January 10 and January 18, 2016, are selected, respectively, at 04:00, 10:00, 16:00, 22:00, and 00:00, 06:00, 12:00, 18:00 (UTC+8). There are obvious differences in the amplitude of sea surface height between spring tides and neap tides. The maximum tide height during a spring tide is above 0.4 m, while the maximum tide height during a neap tide is only half of that during a spring tide, i.e., approximately 0.2 m. Both tides exhibit distinct characteristics of semidiurnal and diurnal cycles, with a phase shift from east to west. The change in tides in high-latitude areas is minimal. From west to east, the sequence of spring tides and neap tides follows a positive-negative-positive-negative pattern in latitude, which is consistent with the findings of the EXP1 experiment by Jin et al. (2022) and is aligned with the actual situation regarding ocean tides.

3.2 Tides simulated by LICOM3.0 and TPXO

We mainly analyzed the largest semidiurnal constituent, M₂, and the largest diurnal constituent, K₁. Figure 2 shows the amplitudes and phase of M₂ and K₁ constituents from LICOM3.0 and TPXO. Overall, the geographical distribution of the amplitude extreme values is approximately the same, and the amphidromic point locations can be mapped. For M₂, the amplitude of the tide simulated by LICOM3.0 is larger than that of TPXO. The global average M₂ amplitudes of LICOM3.0 and TPXO are 37.60 cm and 30.70 cm, respectively. The large amplitude values for TPXO displays are mostly distributed on the continent or the periphery of islands, such as the Bay of Alaska, the Weddell Sea, and the North Atlantic. In addition, in these regions, the amplitude of the M₂ constituent simulated by LICOM3.0 is relatively large, which accounts for the high global average amplitude. The phase of the M₂ constituent differs from TPXO in some areas, such as the Tasman Sea and the Southern Ocean. The gap is obvious in the Southern Ocean, which may be related to the absence of LICOM3.0 data at high latitudes, resulting in a northward shift of the amphidromic points. In contrast to M₂, large amplitude values for K₁ are mainly distributed along the entire coast of Antarctica and around Indonesia.

Similarly, the amplitude of the K₁ constituent simulated by LICOM3.0 is larger than that of TPXO, with global average amplitudes of 11.11 cm and 10.2 cm, respectively. From a spatial distribution perspective, the simulated tides surrounding the Antarctic mainland are particularly strong in LICOM3.0. The position of the amphidromic points simulated by LICOM3.0 is consistent with those simulated by TPXO, with significant differences in the North Atlantic and Indian oceans.

Notably, in both the M₂ and K₁ constituents, the phase contour lines simulated by LICOM3.0 display obvious perturbations in regions such as the North Pacific and the tropical Pacific. These perturbations are more obvious in areas where the internal tides are active. However, there is no such phenomenon in TPXO. These perturbations represent baroclinic tides, which signal internal tides. As TPXO is a two-dimensional barotropic model without internal tide processes, no such signals exist. The horizontal resolution of the LICOM3.0 eddy-resolving OGCM is approximately 10 km, which is much smaller than the shortest wavelength of mode-1 M₂ and K₁ internal tides (Carrere et al., 2021). LICOM3.0 is sufficient to partially resolve the M₂ and K₁ internal tides, resulting in such perturbations in the phase. A similar phenomenon appeared in the results reported by Arbic et al. (2012).

Figure 3 displays the amplitude error

$ {d}_{\mathrm{a}} $

, phase error

$ {d}_{\mathrm{p}} $

, and total error

$ d $

values of the M₂ and K₁ constituents. Table 4 presents the global average results. The results in Table 4 only comprise an ocean between 66°S and 66°N and a depth greater than 1 000 m. For the M₂ constituent, the global average total error of LICOM3.0 is 10.19 cm relative to TPXO9, with an amplitude error of 6.83 cm and a phase error of 6.09 cm. The amplitude error accounts for a relatively large proportion. From a spatial distribution perspective, the maximum value of the M₂ constituent total error appears near Greenland in the North Atlantic and the Ross Sea, and the error is large near other semi-enclosed bays. In areas where the total error is large, the amplitude and phase errors are both large, but the amplitude error is even larger. In such areas, the amplitude of the M₂ constituent is the largest in both TPXO9 and LICOM3.0, and the error in LICOM3.0 is the largest. For the K₁ constituent, the global average total error of LICOM3.0 is 2.85 cm relative to TPXO9, with an amplitude error of 1.66 cm and a phase error of 1.92 cm. The phase error accounts for a relatively large proportion. From a spatial distribution perspective, the total error is large near the Southern Ocean. However, the phase error contributes more to the total error in the K₁ constituent than in the M₂ constituent. ICON-O (von Storch et al., 2023) is a recently developed OGCM with tidal forcing, and we compared the results of LICOM3.0 and ICON-O. Table 4 shows the results of ICON-O (von Storch et al., 2023). A comparison shows that the results of LICOM3.0 are significantly better than those of ICON-O, both in terms of the total error and the errors of the two components. The total error of the two constituents is more than 50% lower in LICOM3.0. In addition, it is worth noting that, unlike LICOM3.0, the amplitude error is larger than the phase error of the M₂ constituent in ICON-O results, while amplitude and phase errors contribute almost equally in the K₁ constituent.

We have also investigated other six constituents (not shown). We found that the results for other constituents are similar to the most dominant semidiurnal (M₂) and diurnal (K₁) constituents.

In summary, we found that the errors of simulated M₂ and K₁ constituents vary in different sea areas in LICOM3.0, and the amplitude and phase contribute differently to the total error. The results from different OGCMs are different, and it is impossible to identify any apparent commonalities. Further analysis of results from additional models is needed.

3.3 LICOM3.0 simulated tides and observed tides

We further utilized the st102 data from Shum et al. (1997) and the data from two tidal gauge stations of World Ocean Circulation Experiment (WOCE) for validation. These observational data are often used for tide simulation verification. Table 5 shows the bias of the eight major constituents in the LICOM3.0 and st102 (Shum et al., 1997) dataset. Unlike TPXO, the data from these 102 stations are observation data, which, although they do not cover as wide a range as TPXO, are representative of the true tides. The comparison shows that the simulation error is still the largest for the M₂ constituent, with a total error of 13.85 cm and a relative error of 41.7%, where the amplitude contributes more to the total error. The total error for the K₁ constituent is 3.44 cm, and the relative error is 30.6%, with the phase contributing more to the total error. Overall, the relative sizes of the total errors for the M₂ and K₁ constituents and the contributions of amplitude and phase to the total errors are the same as those found in the previous section by LICOM3.0 and TPXO9. For comparison, we evaluated the results of other models, with the total errors for the M₂ and K₁ constituents being 8.26 cm and 2.48 cm for HYCOM (Arbic et al., 2010) and 8.91 cm and 2.52 cm for STORMTIDE (Müller et al., 2014), respectively, which are smaller than those from LICOM3.0.

Regarding the percentage of sea surface height variance that we captured, the simulated M₂ and K₁ constituents of LICOM3.0 were found to reach 82.62% and 90.67%, respectively. Among the eight constituents, the accuracy of the simulated N₂ constituent is the best, with a variance of 91.07%. The poorest accuracy constituent was O₁, with a variance of 35.89%. The RSS of the eight constituents in LICOM3.0 is approximately 80.76% relative to the st102 data, which indicates that LICOM3.0 can generally simulate the eight major constituents well. In addition, this RSS was compared with the results of Müller et al. (2014) and Arbic et al. (2010), whose variances reached 92.8% and 92.6%, respectively. In addition, von Storch et al. (2023) increased the variance explained by adjusting the model time step and altering the coefficients of the bottom friction term and the topographic wave drag term; however, this variance was still below 80%. The tide simulation of LICOM3.0 is slightly inferior to that of Müller et al. (2014) and Arbic et al. (2010) but superior to that of von Storch et al. (2023).

Figure 4 displays the spectral analysis of the sea level data from the Yap and Kodiak tide gauge stations (Ponchaut et al., 2001) and the sea surface height data from LICOM3.0. The results from LICOM3.0 at the frequency bands of the eight major constituents are similar to the tide gauge station results; however, the amplitude at the frequency band of the M₂ constituent is significantly lower in LICOM3.0 than in the tide gauge station results. The amplitudes of the major constituents are similar in both the LICOM3.0 and tide gauge station results. Furthermore, there are some errors in low-frequency bands between LICOM3.0 and the tide gauge station results, but tides do not cause these errors and are outside the scope of this discussion. These low-frequency errors have been observed in the results of other articles, such as Jin et al. (2022).

4 Conclusions and discussion

Less

In this study, we incorporated the astronomical tidal potentials of eight major constituents into a global eddy-resolving OGCM, LICOM3.0. Herein, we comprehensively evaluated the tidal simulation ability of the model. The resulting sea surface height was analyzed through tidal harmonic analysis to obtain the harmonic constants. The harmonic constants were compared with the data from TPXO9 and st102 and with previous results. In addition, the sea surface height was compared with WOCE tide gauge stations through spectral analysis. The main conclusions are as follows:

(1) LICOM3.0 can effectively simulate barotropic tides. The amplitude and phase of the M₂ and K₁ constituents simulated by LICOM3.0 are approximately the same as the TPXO data in spatial distribution, and the position of the amphidromic points can also be correspondingly mapped. The simulation error is smaller in the open ocean and larger in the nearshore regions.

(2) Relative to the TPXO data, the total errors for the M₂ and K₁ constituents are 10.19 cm and 2.85 cm, respectively. Specifically, the amplitude error is more significant in the M₂ constituent, whereas the K₁ constituent exhibits a larger phase error. This situation is different from those in other models.

(3) In LICOM3.0, the simulated tidal amplitude is greater than that of TPXO, and the errors in these high-amplitude regions are relatively great.

(4) LICOM performs comparably to other models. Relative to the st102 data, the tidal error in LICOM3.0 is smaller than that of ICON-O (von Storch et al., 2023) but slightly larger than those of STORMTIDE (Müller et al., 2014) and HYCOM (Arbic et al., 2010).

In the discussion of this article, the tidal results simulated by LICOM3.0 were found to be suboptimal, indicating a need for further adjustment of the coefficients of the topographic wave drag term. Although this paper featured only an initial evaluation of the tidal simulation in the eddy-resolving version of LICOM3.0, the simulated tides were sufficiently accurate to capture the signal of internal tides, allowing for the study of some statistical characteristics. To investigate the influences of simulated tides on ocean energy distribution and circulation, we should further explore the synergistic impact of internal tides and circulation, building on existing barotropic tide simulation results.

Funding

Less

The National Natural Science Foundation of China under contract Nos 41931182, 42090040, 42176024, and 42206006; theNational Key Program for Developing Basic Sciences under contract No. 2022YFC3104802.

References

Less

Ansong J K, Arbic B K, Alford M H, et al. 2017. Semidiurnal internal tide energy fluxes and their variability in a global ocean model and moored observations. Journal of Geophysical Research: Oceans, 122(3): 1882–1900, doi: 10.1002/2016JC012184

Ansong J K, Arbic B K, Buijsman M C, et al. 2015. Indirect evidence for substantial damping of low-mode internal tides in the open ocean. Journal of Geophysical Research: Oceans, 120(9): 6057–6071, doi: 10.1002/2015JC010998

Arbic B K, Alford M H, Ansong J K, et al. 2018. A primer on global internal tide and internal gravity wave continuum modeling in HYCOM and MITgcm. In: Chassignet E P, Pascual A, Tintoré J, et al, eds. New Frontiers in Operational Oceanography. GODAE OceanView, 308–334

Arbic B K, Garner S T, Hallberg R W, et al. 2004. The accuracy of surface elevations in forward global barotropic and baroclinic tide models. Deep-Sea Research Part II: Topical Studies in Oceanography, 51(25/26): 3069–3101, doi: 10.1016/j.dsr2.2004.09.014

Arbic B K, Richman J G, Shriver J F, et al. 2012. Global modeling of internal tides within an eddying ocean general circulation model. Oceanography, 25(2): 20–29, doi: 10.5670/oceanog.2012.38

Arbic B K, Wallcraft A J, Metzger E J. 2010. Concurrent simulation of the eddying general circulation and tides in a global ocean model. Ocean Modelling, 32(3/4): 175–187, doi: 10.1016/j.ocemod.2010.01.007

Buijsman M C, Ansong J K, Arbic B K, et al. 2016. Impact of parameterized internal wave drag on the semidiurnal energy balance in a global ocean circulation model. Journal of Physical Oceanography, 46(5): 1399–1419., doi: 10.1175/JPO-D-15-0074.1

Buijsman M C, Stephenson G R, Ansong J K, et al. 2020. On the interplay between horizontal resolution and wave drag and their effect on tidal baroclinic mode waves in realistic global ocean simulations. Ocean Modelling, 152: 101656, doi: 10.1016/j.ocemod.2020.101656

Carrere L, Arbic B K, Dushaw B, et al. 2021. Accuracy assessment of global internal-tide models using satellite altimetry. Ocean Science, 17(1): 147–180, doi: 10.5194/os-17-147-2021

Codiga D L. 2011. Unified tidal analysis and prediction using the UTide matlab functions. Technical Report 2011-01. Narragansett: Graduate School of Oceanography, University of Rhode Island, 59

Egbert G D. 1997. Tidal data inversion: interpolation and inference. Progress in Oceanography, 40(1–4): 53–80, doi: 10.1016/S0079-6611(97)00023-2

Egbert G D, Erofeeva S Y. 2002. Efficient inverse modeling of Barotropic ocean tides. Journal of Atmospheric and Oceanic Technology, 19(2): 183–204, doi: 10.1175/1520-0426(2002)019<0183:EIMOBO>2.0.CO;2

Egbert G D, Ray R D. 2000. Significant dissipation of tidal energy in the deep ocean inferred from satellite altimeter data. Nature, 405(6788): 775–778, doi: 10.1038/35015531

Fisheries and Oceans Canada. 2017. Table for the astronomical argument V+U. https://www.dfo-mpo.gc.ca/science/data-donnees/tidal-marees/argument-u-v-eng.html [2017-01-26/2023-09-15]

Fisheries and Oceans Canada. 2018. Table of values for the node factor f. https://www.dfo-mpo.gc.ca/science/data-donnees/tidal-marees/facteur-node-factor-eng.html[2018-02-27/2023-09-15]

Hendershott M C. 1972. The effects of solid earth deformation on global ocean tides. Geophysical Journal of the Royal Astronomical Society, 29(4): 389–402, doi: 10.1111/j.1365-246X.1972.tb06167.x

Jayne S R. 2009. The impact of abyssal mixing parameterizations in an ocean general circulation model. Journal of Physical Oceanography, 39(7): 1756–1775, doi: 10.1175/2009JPO4085.1

Jayne S R, St. Laurent L C. 2001. Parameterizing tidal dissipation over rough topography. Geophysical Research Letters, 28(5): 811–814, doi: 10.1029/2000GL012044

Jin Jiangbo, Guo Run, Zhang Minghua, et al. 2022. Formulation of a new explicit tidal scheme in revised LICOM2.0. Geoscientific Model Development, 15(10): 4259–4273, doi: 10.5194/gmd-15-4259-2022

Kostov Y, Armour K C, Marshall J. 2014. Impact of the Atlantic meridional overturning circulation on ocean heat storage and transient climate change. Geophysical Research Letters, 41(6): 2108–2116, doi: 10.1002/2013GL058998

Li Yiwen, Liu, Hailong, Ding Mengrong, et al. 2020. Eddy-resolving simulation of CAS-LICOM3 for phase 2 of the ocean model intercomparison project. Advances in Atmospheric Sciences, 37(10): 1067–1080, doi: 10.1007/s00376-020-0057-z

Li Zhuhua, von Storch J S. 2020. M₂ internal-tide generation in STORMTIDE2. Journal of Geophysical Research: Oceans, 125(8): e2019JC015453, doi: 10.1029/2019JC015453

Logemann K, Linardakis L, Korn P, et al. 2021. Global tide simulations with ICON-O: testing the model performance on highly irregular meshes. Ocean Dynamics, 71(1): 43–57, doi: 10.1007/s10236-020-01428-7

Marshall J, Speer K. 2012. Closure of the meridional overturning circulation through Southern Ocean upwelling. Nature Geoscience, 5(3): 171–180, doi: 10.1038/ngeo1391

Melet A, Hallberg R, Legg S, et al. 2013. Sensitivity of the ocean state to the vertical distribution of internal-tide-driven mixing. Journal of Physical Oceanography, 43(3): 602–615, doi: 10.1175/JPO-D-12-055.1

Müller M, Arbic B K, Richman J G, et al. 2015. Toward an internal gravity wave spectrum in global ocean models. Geophysical Research Letters, 42(9): 3474–3481, doi: 10.1002/2015GL063365

Müller M, Cherniawsky J Y, Foreman M G G, et al. 2012. Global M ₂ internal tide and its seasonal variability from high resolution ocean circulation and tide modeling. Geophysical Research Letters, 39(19): l19607

Müller M, Cherniawsky J Y, Foreman M G G, et al. 2014. Seasonal variation of the M ₂ tide. Ocean Dynamics, 64(2): 159–177, doi: 10.1007/s10236-013-0679-0

Müller M, Haak H, Jungclaus J H, et al. 2010. The effect of ocean tides on a climate model simulation. Ocean Modelling, 35(4): 304–313, doi: 10.1016/j.ocemod.2010.09.001

Munk W, Wunsch C. 1998. Abyssal recipes II: energetics of tidal and wind mixing. Deep-Sea Research Part I: Oceanographic Research Papers, 45(12): 1977–2010, doi: 10.1016/S0967-0637(98)00070-3

Murray R J. 1996. Explicit generation of orthogonal grids for ocean models. Journal of Computational Physics, 126(2): 251–273, doi: 10.1006/jcph.1996.0136

National Centers for Environmental Information. 2006. 2-Minute Gridded Global Relief Data (ETOPO2) v2. Beijing, China: National Geophysical Data Center, National Centers for Environmental Information, doi: 10.7289/V5J1012Q

Niwa Y, Hibiya T. 2014. Generation of baroclinic tide energy in a global three-dimensional numerical model with different spatial grid resolutions. Ocean Modelling, 80: 59–73, doi: 10.1016/j.ocemod.2014.05.003

Peng Shiqiu, Liao Jiawen, Wang Xiaowei, et al. 2021. Energetics-based estimation of the diapycnal mixing induced by internal tides in the Andaman Sea. Journal of Geophysical Research: Oceans, 126(4): e2020JC016521, doi: 10.1029/2020JC016521

Ponchaut F, Lyard F, Le Provost C. 2001. An analysis of the tidal signal in the WOCE Sea level dataset. Journal of Atmospheric and Oceanic Technology, 18(1): 77–91, doi: 10.1175/1520-0426(2001)018<0077:AAOTTS>2.0.CO;2

Pugh D T. 1996. Tides, Surges and Mean Sea-Level (Reprinted with Corrections). Chichester: John Wiley & Sons Ltd, 59–140

Rocha C B, Chereskin T K, Gille S T, et al. 2016a. Mesoscale to submesoscale wavenumber spectra in drake passage. Journal of Physical Oceanography, 46(2): 601–620, doi: 10.1175/JPO-D-15-0087.1

Rocha C B, Gille S T, Chereskin T K, et al. 2016b. Seasonality of submesoscale dynamics in the Kuroshio Extension. Geophysical Research Letters, 43(21): 11304–11311

Saenko O A, Merryfield W J. 2005. On the effect of topographically enhanced mixing on the global ocean circulation. Journal of Physical Oceanography, 35(5): 826–834, doi: 10.1175/JPO2722.1

Savage A C, Arbic B K, Alford M H, et al. 2017a. Spectral decomposition of internal gravity wave sea surface height in global models. Journal of Geophysical Research: Oceans, 122(10): 7803–7821, doi: 10.1002/2017JC013009

Savage A C, Arbic B K, Richman J G, et al. 2017b. Frequency content of sea surface height variability from internal gravity waves to mesoscale eddies. Journal of Geophysical Research: Oceans, 122(3): 2519–2538, doi: 10.1002/2016JC012331

Schiller A. 2004. Effects of explicit tidal forcing in an OGCM on the water-mass structure and circulation in the Indonesian throughflow region. Ocean Modelling, 6(1): 31–49, doi: 10.1016/S1463-5003(02)00057-4

Schiller A, Fiedler R. 2007. Explicit tidal forcing in an ocean general circulation model. Geophysical Research Letters, 34(3): L03611

Shum C K, Woodworth P L, Andersen O B, et al. 1997. Accuracy assessment of recent ocean tide models. Journal of Geophysical Research: Oceans, 102(C11): 25173–25194, doi: 10.1029/97JC00445

Simmons H L, Jayne S R, St. Laurent L C, et al. 2004. Tidally driven mixing in a numerical model of the ocean general circulation. Ocean Modelling, 6(3/4): 245–263, doi: 10.1016/S1463-5003(03)00011-8

Siyanbola O Q, Buijsman M C, Delpech A, et al. 2023. Remote internal wave forcing of regional ocean simulations near the U.S. West Coast. Ocean Modelling, 181: 102154, doi: 10.1016/j.ocemod.2022.102154

Song Pengyang, Sidorenko D, Scholz P, et al. 2023. The tidal effects in the finite-volumE sea ice–ocean model (FESOM2.1): a comparison between parameterised tidal mixing and explicit tidal forcing. Geoscientific Model Development, 16(1): 383–405, doi: 10.5194/gmd-16-383-2023

St. Laurent L C, Simmons H L Jayne S R. 2002. Estimating tidally driven mixing in the deep ocean. Geophysical Research Letters, 29(23): 21

Stammer D, Ray R D, Andersen O B, et al. 2014. Accuracy assessment of global barotropic ocean tide models. Reviews of Geophysics, 52(3): 243–282, doi: 10.1002/2014RG000450

Talley L D. 2013. Closure of the global overturning circulation through the Indian, Pacific, and southern oceans: schematics and transports. Oceanography, 26(1): 80–97, doi: 10.5670/oceanog.2013.07

Thomas M, Sündermann J, Maier-Reimer E. 2001. Consideration of ocean tides in an OGCM and impacts on subseasonal to decadal polar motion excitation. Geophysical Research Letters, 28(12): 2457–2460, doi: 10.1029/2000GL012234

Tsujino H, Urakawa S, Nakano H, et al. 2018. JRA-55 based surface dataset for driving ocean-sea-ice models (JRA55-do). Ocean Modelling, 130: 79–139, doi: 10.1016/j.ocemod.2018.07.002

von Storch J S, Hertwig E, Lüschow V, et al. 2023. Open-ocean tides simulated by ICON-O, version icon-2.6. 6. Geoscience Model Development, 16(17): 5179–5196, doi: 10.5194/gmd-16-5179-2023

Wang Xiaowei, Peng Shiqiu, Liu Zhiyu, et al. 2016. Tidal mixing in the South China Sea: an estimate based on the internal tide energetics. Journal of Physical Oceanography, 46(1): 107–124, doi: 10.1175/JPO-D-15-0082.1

Yu Yi, Liu Hailong, Lan Jian. 2016. The influence of explicit tidal forcing in a climate ocean circulation model. Acta Oceanologica Sinica, 35(9): 42–50, doi: 10.1007/s13131-016-0931-9

Yu Zipeng, Liu Hailong, Lin Pengfei. 2017. A numerical study of the influence of tidal mixing on Atlantic Meridional Overturning Circulation (AMOC) simulation. Chinese Journal of Atmospheric Sciences (in Chinese), 41(5): 1087–1100

Yu Yi, Liu Hailong, Lin Pengfei, et al. 2020. The impact of oceanic processes on the transient climate response: a tidal forcing experiment. Acta Oceanologica Sinica, 39(1): 52–62, doi: 10.1007/s13131-019-1466-0

Zuo Juncheng, Du Ling, Chen Meixiang, et al. 2018. Analysis Method of Ocean Hydrological Environmental Factors (in Chinese). Beijing, China: Science Press, 249–252

Appendix

Less

Year 2024 volume 43 Issue 9

PDF

Cite this Article

BibTeX

Article Info

doi: 10.1007/s13131-024-2352-5

Receive Date：2023-11-26
Online Date：2025-11-19
Published：2024-09-25

Article Data

Affiliations

History

Received：2023-11-26
Accepted：2024-04-27

Funding

The National Natural Science Foundation of China under contract Nos 41931182, 42090040, 42176024, and 42206006; theNational Key Program for Developing Basic Sciences under contract No. 2022YFC3104802.

Affiliations

¹ Physical Oceanography Laboratory, Ocean University of China, Qingdao 266100, China

² State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

³ Laoshan Laboratory, Qingdao 266237, China

⁴ College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China

⁵ Hohai University, Nanjing 210098, China

Corresponding:

* Email: yzp@lasg.iap.ac.cn

lhl@lasg.iap.ac.cn

References

Share

https://castjournals.cast.org.cn/joweb/aos/EN/10.1007/s13131-024-2352-5

Share to

Scan QR to access full text

Cite this article

BibTeX

Citations

表12种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

Table 1. Calculation formulas for the nodal factor f and the nodal angle μ

Constituent	f	μ/(°)
K₁	1.0060 + 0.1150 cos $ {N}_{0} $ – 0.0088 cos(2 $ {N}_{0} $) + 0.0006 cos(3 $ {N}_{0} $)	–8.86° sin $ {N}_{0} $ + 0.68° sin(2 $ {N}_{0} $) − 0.07° sin(3 $ {N}_{0} $)
O₁	1.0089 + 0.1871 cos $ {N}_{0} $ − 0.0147 cos(2 $ {N}_{0} $) + 0.0014 cos(3 $ {N}_{0} $)	10.80° sin $ {N}_{0} $ − 1.34° sin(2 $ {N}_{0} $) + 0.19° sin(3 $ {N}_{0} $)
P₁	1	0
Q₁	same as O₁	same as O₁
M₂	1.0004 − 0.0373 cos $ {N}_{0} $ + 0.0003 cos(2 $ {N}_{0} $)	–2.14° sin $ {N}_{0} $
S₂	1	0
N₂	same as M₂	same as M₂
K₂	1.0241 + 0.2863 cos $ {N}_{0} $ + 0.0083 cos(2 $ {N}_{0} $) − 0.0015 cos(3 $ {N}_{0} $)	–17.74° sin $ {{N}}_{0} $ + 0.68° sin(2 $ {N}_{0} $) − 0.04° sin(3 $ {N}_{0} $)

Table 2. Total (d), amplitude (d _a), and phase (d _p) errors of the four major constituents of LICOM3.0 relative to TPXO

Experiment	Error/cm
	M₂			S₂			K₁			O₁
	d	d _a	d _p	d	d _a	d _p	d	d _a	d _p	d	d _a	d _p
wd0.5	9.86	7.21	5.37	5.12	4.00	2.47	3.58	2.36	2.25	2.59	1.26	1.99
wd1	10.70	8.16	5.45	5.38	4.25	2.57	3.77	2.49	2.36	2.57	1.25	1.94
wd1.5	11.39	8.88	5.59	5.59	4.43	2.67	3.92	2.59	2.45	2.57	1.26	1.91

Table 3. Phases-dependent $ {V}_{0}+\mu $ and nodal factor f values of each constituent calculated by DFO, tidal harmonic analysis program UTIDE and LICOM3.0 tidal module in 2016, corresponding to the dates of January 1 and July 1 of that year

Constituent	( $ {V}_{0}+\mu $) (Jan. 1)/(°)			f (Jul. 1)
Constituent	DFO	UTIDE	LICOM3.0	DFO	UTIDE	LICOM3.0
M₂	210.745	210.729	210.857	1.037	1.037	1.037
K₁	9.242	9.227	9.299	0.886	0.886	0.886
O₁	202.050	202.106	202.000	0.809	0.809	0.813
S₂	0.001	0	0	0.998	0.998	1.000
P₁	349.847	349.852	349.903	1.011	1.011	1.000
Q₁	43.050	43.384	41.752	0.821	0.823	0.813
N₂	50.690	50.690	50.610	1.037	1.037	1.037
K₂	198.451	198.412	198.722	0.752	0.752	0.755

Table 4. Errors relative to TPXO in LICOM3.0 and ICON-O (von Storch et al., 2023), including total error $ d $, amplitude error $ {d}_{\mathrm{a}} $, and phase error $ {d}_{\mathrm{p}} $

Constituent	Error/cm
	LICOM3.0			ICON-O
	$ d $	$ {d}_{\mathrm{a}} $	$ {d}_{\mathrm{p}} $	$ d $	$ {d}_{\mathrm{a}} $	$ {d}_{\mathrm{p}} $
M₂	10.19	6.83	6.09	14.25	10.63	9.39
K₁	2.85	1.66	1.92	4.60	3.26	3.25
O₁	4.30	2.55	2.98	5.19	3.89	3.43
S₂	5.37	3.12	3.84	7.90	4.71	1.50
P₁	1.12	0.52	0.86	1.49	1.08	1.02
Q₁	0.77	0.47	0.56	0.86	0.48	0.72
N₂	1.51	0.87	1.06	2.30	1.50	1.74
K₂	1.19	0.78	0.73	2.30	1.19	1.96

Table 5. Signal $ S $, total error $ d $, amplitude error $ {d}_{\mathrm{a}} $, and phase error $ {d}_{\mathrm{p}} $ of the eight major constituents of LICOM3.0 relative to st102 data, the percentages of captured sea surface height variance $ {V}_{\mathrm{c}\mathrm{a}\mathrm{p}\mathrm{t}} $ for each constituent of LICOM3.0, STORMTIDE (Müller et al., 2014), and HYCOM (Arbic et al., 2010), and the root square sum (RSS) values of the eight major constituents

Constituent and RSS	Error/cm				$ {V}_{\mathrm{capt}} $/%
Constituent and RSS	$ S $	$ d $	$ {d}_{\mathrm{a}} $	$ {d}_{\mathrm{p}} $	LICOM3.0	STORMTIDE	HYCOM
M₂	33.22	13.85	9.56	8.12	82.62	93.9	93.8
K₁	11.26	3.44	1.89	2.39	90.67	95.0	95.1
O₁	7.76	6.21	3.31	4.62	35.89	83.2	89.7
S₂	12.62	6.35	4.33	4.03	74.66	86.9	83.2
P₁	3.62	1.61	0.73	1.26	80.22	94.7	95.2
Q₁	1.62	1.06	0.66	0.73	57.19	64.7	82.1
N₂	6.86	2.02	1.36	1.21	91.07	96.0	95.9
K₂	3.43	1.58	1.23	0.74	78.78	89.7	76.9
RSS	39.04	17.11	11.36	10.64	80.76	92.8	92.6

Fig. 1. Changes in the sea surface caused by the equilibrium tide calculated by the LICOM3.0 tidal module: the results during a spring tide at 04:00, 10:00, 16:00, and 22:00 (a–d) and the results during a neap tide at 00:00, 06:00, 12:00, and 18:00 (e–h). Time in UTC+0.

Fig. 2. Spatial patterns of the amplitude (colored, cm) and phase (contour, °) characteristics of the M₂ constituent for TPXO (a) and LICOM3.0 (b) and of the K₁ constituent for TPXO (c) and LICOM3.0 (d). The black solid lines are separated by 30°.

Fig. 3. Errors of the simulated M₂ constituent (a–c) and K₁ constituent (d–f) relative to TPXO in LICOM3.0: amplitude error (a, d), phase error (b, e), and total error (c, f).

Fig. 4. Spectral analysis of sea levels observed (obs) by Yap (9.5083°N, 138.1283°E) and Kodiak (57.7317°N, 152.511 7°W) in WOCE (Ponchaut et al., 2001) and simulated by LICOM3.0.

Articles: Latest Articles; Most Read; Collections

Updates: Events; News; Multimedia

About: About Us

Contact

No. 86 Xueyuan South Road, Haidian District, Beijing

100081

010-62199257

qkjq@cast.org.cn

Copyright © 2025 China Association for Science and Technology. All rights reserved. For all open access content, the relevant licensing terms apply.
Sponsored by the Office of the Leading Group for Cybersecurity and Informatization of CAST, and supported by Science and Technology Review Publishing House