Modulation of Internal Tides by Turbulent Mixing in the South China Sea

Modulations of internal tides (ITs) including the baroclinic tidal energy budget, the incoherency, and the nonlinear interactions among different tidal components by turbulent mixing in the South China Sea (SCS) are investigated through numerical simulations. The baroclinic tidal energy budget can hardly be affected by the structure of mixing. Meanwhile, change in the mixing intensity in a reasonable range also cannot obviously modulate the baroclinic tidal energy budget in the SCS. Compared to the baroclinic energy budget, the distributions of conversion and dissipation are more sensitive to the change of mixing. Turbulent mixing also modulates the incoherency of ITs by changing the horizontal density in the ocean. The horizontal variation of density adds incoherence to ITs largely by affecting the internal tidal amplitudes. Furthermore, nonlinear interactions among different components of ITs are generally modulated by the mixing intensity, whereas the variation of the mixing structure can hardly influence the nonlinear interactions. Therefore, the diapycnal diffusivity can set to be horizontally and vertically homogeneous in most of the internal tidal simulations, except for those in which the incoherency of ITs needs to be simulated. However, excessive strong mixing will destroy the stratification. Thus, the optimum range for IT simulations in the SCS is from O (10–5) to O (10–3) m2s–1.


INTRODUCTION
Internal tides (ITs) are internal waves generated as surface tides flow over rough topography such as seamounts and ridges (Nash et al., 2006;Xu et al., 2014;Cao et al., 2017). ITs play important roles in enhancing mixing (Egbert and Ray, 2000;Niwa and Hibiya, 2001;Tian et al., 2010;Zhao, 2014;Peng et al., 2021). Turbulent mixing provides energy for the upwelling of deep-water circulation, which is an essential process of global overturning circulation (Toggweiler and Samuels, 1995;Webb and Suginohara, 2001). However, the averaged diapycnal diffusivity, which is necessary to maintain the global meridional overturning circulation, is one order of magnitude larger than those observed in the open ocean (Munk, 1966). In other words, turbulent mixing in regional seas will be much stronger than that in the open ocean.
ITs represent intermediate steps, which transfer energy from surface tides to turbulence mixing (Foreman et al., 2004;Carter et al., 2006;Liu, 2010;Polzin et al., 2014). Breaking of ITs contributes apparently to regional turbulent mixing. The diapycnal diffusivities in regions with strong ITs can be more than three orders of magnitude greater (Lueck and Osborn, 1985;Rudnick et al., 2003;Carter et al., 2005). Haren et al. (2015) found that at the steepest slope, turbulence was 100 times more energetic than at the shallowest slope. Munk (1998) found that ITs provide approximately half of the energy necessary for the global overturning circulation. Carter et al. (2005) also attributed the enhancement of turbulence to the presence of ITs and nonlinear internal waves. Using observations and numerical simulations, Rudnick et al. (2003) investigated tides-to-turbulence cascade along the Hawaiian Ridge and suggested that tide-induced mixing is evidently stronger than that in the open ocean and isolated ridges are essential sites for enhanced mixing. Johnston et al. (2011) found that averaged turbulence mixing is 50% stronger along the internal tidal beams than outside the beams near the Monterey Bay and peak values in the beam can be an order of magnitude larger. Based on the observations, Zhao et al. (2012) examined the bottom intensification of turbulence in the Monterey Submarine Canyon and suggested that the temporal varying of mixing is modulated by ITs. Legg (2003) investigated the high-mode structure of ITs and the enhancement of mixing near the continental slope and pointed out that shear may lead to local mixing. MacKinnon and Gregg (2005) also attributed the elevation of turbulent mixing to the increment of shear, which is caused by the reinforcement of stratification. Strong near-surface stratification can apparently influence mixing by isolating the deeper water from mechanical penetration of atmospheric energy (Jinadasa et al., 2016).
The South China Sea (SCS), connected by the Luzon Strait (LS) to the West Pacific, is abundant in internal waves such as ITs, nonlinear internal waves, and near-inertial waves, which results in the enhancement of turbulent mixing. Liu and Lozovatsky (2012) investigated the turbulent kinetic energy dissipation rate in the upper pycnocline of the northern SCS. They found that internal waves originated in the LS are associated with the predominant mixing that appears in regions north of 20 • N. Wang et al. (2017) investigated the influence of ITs on the circulation and mass transformation in the SCS. They suggested that tidal mixing would strengthen the horizontal density gradient and the basin-scale cyclonic circulation, promoting overturning circulation in the SCS. Yang et al. (2016) observed spatial distribution of turbulent mixing in the SCS, based on the Gregg-Henyey-Polzin parameterization. They revealed that turbulent mixing generally increases with depth. In the horizontal direction, mixing is elevated in the south and east regions of the SCS.
Furthermore, temporal and spatial variations of mixing cause corresponding changes in stratification. Delorme et al. (2021) found that enhanced mixing results in the weakening of the abyssal stratification. Ivey (2004) found that during the period of strong mixing, the stratification is weakened or even absent. Drake et al. (2020) suggested that mixing, which generates flows in a bottom boundary layer, can effectively affect the global ocean stratification. Garrett (2003) also found that mixing affects ocean stratification and ocean circulation. Variations of stratification can effectively modulate ITs. Foreman and Thomson (1997) found that seasonal changes in stratification not only influence the intensity of ITs, enhancing bottom currents, but also eventually affect surface tides. Dushaw (2002) pointed out that variation of stratification can cause changes in phase speed of mode-1 ITs. Guo et al. (2012) found that the increased stratification in autumn is beneficial to propagating internal tidal energy and to preserving higher-mode ITs. They also suggested that changes of stratification play important roles in the dissipation of ITs. Nakamura et al. (2014) found that the phase relationship with seasonal variation in stratification affects the strength of tidal mixing. They also suggested that the strength of tide-induced mixing depends on the temporal variations of vertical diffusivity. In other words, diffusivity is envisioned as a key character, influencing the ITs through stratification. However, the modulation of ITs by turbulent mixing is rarely investigated, especially in the SCS, and the clear picture to a large extent remains unknown. Moreover, the diapycnal diffusivity is usually set to be spatial uniform when we simulate ITs. The spatial structure of mixing can effectively influence ITs is not clear. Therefore, investigations need to be carried to better understand the modulation of ITs by turbulent mixing. In this study, simulations of ITs based on the Massachusetts Institute of Technology General Circulation Model (MITgcm) are performed to explore modulations of ITs in the SCS by mixing. This study is organized as follows. Methods of model configuration and the baroclinic tidal energy budget analysis are introduced in section "Materials and Methods". Modulation of ITs by mixing is presented in section "Results". Summary and discussion are shown in section "Summary and Discussion."

Baroclinic Tide Model Configuration Method
The ocean model used in this study is the MITgcm, which has been widely used for simulating ITs around the LS and in the SCS (Buijsman et al., 2012;Wang et al., 2016). The simulation area includes the LS and part of the SCS (12 • -24 • N, 104 • -125 • E), which is shown in Figure 1A and has a horizontal resolution of 1/24 • × 1/24 • . The model topography is from the General Bathymetric Chart of the Oceans (GEBCO) bathymetry data with a high resolution of 30 arcs. The initial temperature and salinity profiles are derived from the monthly mean climatology of the Generalized Digital Environmental Model, version 3 (GDEMv3). The initial fields are set to be horizontally homogeneous by using temperature and salinity at sites 20.50 • N and 120.25 • E in January. In the vertical direction, there are 60 uneven vertical layers, varying from 0 at the top to 5,700 m at the bottom ( Figure 1B). The model is forced by barotropic tidal currents at the open boundaries. Principal constituents of diurnal and semidiurnal ITs, including K 1 , O 1 , M 2 , and S 2 , are taken into consideration. The amplitudes and phases of these constituents are extracted from the regional solution for the China Sea of the Oregon State University inverse barotropic tidal model (OTIS). A 0.5 • width sponge layer is applied.
The mixing parameterization is set as follows. The horizontal eddy viscosity and vertical eddy viscosity are 100 m 2 s −1 and 10 −4 m 2 s −1 , respectively. Configurations of the diapycnal diffusivity are estimated according to observations of Yang et al. (2016) and the Gregg-Henyey-Polzin parameterization is employed. Experiments are carried out to investigate modulations of generation, dissipation, and coherency of ITs and nonlinear interaction among different components of ITs. In the first group of experiments (IE1), IE1-1 is conducted with a horizontally and vertically homogeneous value of the diapycnal diffusivity of 1.7 × 10 −3 m 2 s −1 , which is the mean value of the observed diapycnal diffusivity in the SCS according to Yang et al. (2016). IE1-2 is conducted with a three-dimensional spatial structure of the diapycnal diffusivity. The distribution of the diapycnal diffusivity along a section of 20 • N is shown in Figure 2A. Because the observation is relatively sparse, the three-dimensional structure of the diapycnal diffusivity can not be accurately obtained by interpolation. The spatial structure, thus, can be estimated from the horizontal distribution and is calculated as: where, K hz is the horizontal structure of the diapycnal diffusivity, K 0 = 1.0 × 10 −5 m 2 s −1 is the background diffusivity, and F is the vertical structure function, which is calculated according to Wang et al. (2016).
Here, D = H +ξ is the total water depth, ξ is the time-mean sea level, H is the water depth, and β is the vertical decay scale. We vertically average observed data to calculate the horizontal structure. However, the area at the east ridge of the LS is not covered by observation, so the horizontal structure of the diapycnal diffusivity cannot be given directly. Considering that the LS is the generation site of ITs, which will effectively  , assuming that the distribution of topography and the distribution of the diapycnal diffusivity between east and west ridges are symmetrical. Therefore, horizontal distribution of the diapycnal diffusivity in the LS is obtained ( Figure 2B). Then, the three-dimensional structure of the diapycnal diffusivity can be calculated according to Equation 1.
In addition, moored current observation in the northern SCS is used for testifying the accuracy of the baroclinic tide simulations. The water depth in the mooring site is 2,480 m. As a part of the SCS Internal Wave Experiment, the subsurface mooring (21.05 • N, 118.45 • E) was deployed to observe the layers near the surface (from 80 to 355 m) ( Figure 1A). The velocity data were interpolated onto uniform levels with 5 m intervals and hourly averaged data from August 2010 to January 2011 were used for checking the accuracy of IT simulation.

Baroclinic Tidal Energy Budget Analysis Method
The internal tidal energy is approximately calculated according to Wang et al. (2016).
Where, C is the depth-integrated conversion of barotropic to baroclinic energy, FD is the divergence of depth-integrated baroclinic energy flux, and DIS is the dissipation of baroclinic energy. C and FD are given as: where, g is the acceleration of gravity, ρ is the density perturbation, wbt is the vertical velocity of barotropic tides, and u' and p' are the horizontal baroclinic velocity (u and v for the eastward and northward direction, respectively) and pressure perturbation, respectively. ρ and wbt can be calculated as: where, ρ is the instantaneous density,ρ is the density averaged in a tidal period, and σ is defined as σ = (z − ξ )/D. u and p are given as: where,ū is the velocity averaged in a tidal period. Figure 3 shows a comparison between simulated and observed amplitudes of eastward velocities of M 2 , S 2 , K 1 , and O 1 in IE1-1 (result of IE1-2 is similar and not shown). Table 1 displays the differences (vertical mean absolute errors of amplitudes) between observations and simulated results. Errors of semidiurnal ITs are smaller than those of diurnal counterparts. Vertical mean absolute errors of both the M 2 and S 2 are 0.6 cm/s. For diurnal constituents, larger differences in K 1 mostly appear beneath 200 m. The simulated O 1 is largely stronger than observation and makes the difference larger than other components. The conversion and dissipation of the baroclinic tidal energy budget in all the experiments are shown in Table 2. ITs are largely generated in the LS, but more energy is dissipated in the SCS. Differences of conversion between IE1-1 and IE1-2 are small in both the LS and the SCS. When structures of turbulent mixing alters from spatial uniform to spatial varying, baroclinic energy budgets of IE1-1 and IE1-2 are close. This suggests that the spatial structure of the diapycnal diffusivity can hardly modulate the baroclinic tidal energy budget in both the LS and the SCS.

Modulation on the Baroclinic Tidal Energy Budget
Group of idealized experiments 2 (IE2) is further conducted to investigate the influence of mixing intensity on ITs. IE2-1 and IE2-2 are conducted with the diapycnal diffusivity of 9.3 × 10 −6 m 2 s −1 and 1.0 × 10 −1 m 2 s −1 , which are the minimum and maximum of the observed diapycnal diffusivity, respectively, in the SCS according to Yang et al. (2016). Differences in the baroclinic tidal energy budget between IE2-1 and IE1-1 are still extraordinarily small even the intensity of turbulent mixing in IE2-1 becomes about two orders of magnitudes smaller. However, apparent differences emerge, when the intensity of the diapycnal diffusivity increases by about two orders of magnitudes and reaches the observed maximum. In both the LS and the SCS, barotropic-to-baroclinic conversions are decreased evidently. Meanwhile, modulations of energy dissipation also are evident. But, the dissipation shows opposite varying character between the LS and the SCS. The dissipation will be reinforced in the LS, but weakened in the SCS when mixing enhances. ITs in the SCS are largely generated in the LS. The decrease in conversion and the increase in dissipation in the LS reduce the energy, which propagated into the SCS and eventually degrade the dissipation rate in the SCS. In addition, the dissipation in the LS is more sensitive to the change of mixing. In the LS, dissipation in IE2-2 is 51% stronger than that in IE1-1. But, dissipation in the SCS is only 20% weaker.
Group of idealized experiments 3 (IE3) is carried out to testifying whether the spatial variation of mixing in vertical and horizontal directions both can hardly influence the total energy budget or they have the opposite effect, which cancels each other out. IE3-1 is conducted with a vertical varying structure of mixing. Vertical structure of the diapycnal diffusivity can be obtained by averaging the observation horizontally and values under 5,000 m depth are set to be the same as that at 5,000 m because of the lack of observation under this depth. IE3-2 is conducted with a horizontal varying structure of mixing ( Figure 2B). Results show that both the conversion and dissipation are close to each other. Differences between IE3 and IE1-1 are generally smaller than those between IE1-1 and IE1-2, which indicate that the spatial variation of mixing in both the vertical and horizontal directions can hardly influence the baroclinic tidal energy budget. Figure 4 shows simulated distributions of conversion rate during the period of spring tide in all the experiments. The black contours represent strong-conversion regions in which the conversion rate is equal to or larger than 0.5 Wm −2 . Patterns of conversion in IE1-2 are similar to that in IE1-1, which suggest that distribution of conversion can hardly be modulated when mixing alters from spatially uniform form to three-dimensional varying structure. Moreover, patterns of conversion between IE1-1 and IE3-1 are also very close. However, the distinction of distribution in the SCS between IE1-1 and IE3-2 is relatively evident. When mixing alters from spatially uniform to simple horizontal varying structure, regions of positive conversion enlarge, but areas of strong conversion (indicates by the black lines) reduce in the SCS. This indicated that compared with threedimensional changes of mixing structure, the simple variation in the horizontal direction has more effect on the distribution of conversion rate, especially in the SCS. The distinction of distribution is more evident when mixing intensity changes. Although the baroclinic tidal energy budget is close between IE1-1 and IE2-1, distributions show more obvious distinctions.
Regions of positive conversion also enlarge and areas of strong conversion reduce when mixing intensity weakens. However, when the mixing intensity strengthens, large parts of regions with positive conversion turn into negative conversion in the SCS, which eventually reduces the total energy converted from barotropic tides to ITs. Figure 5 shows simulated distributions of dissipation rate during the period of spring tide in all the experiments. Gray contours indicate the strong-dissipation regions in which the dissipation rate is equal to or larger than 0.5 Wm −2 . The distinction of dissipation distribution is largely concentrated near the LS. In IE1-2, the area of dissipation narrows, but the strength of dissipation increases in the LS. Differences of dissipation distribution also exist in IE3-1 and IE3-2, but not as evident as that in IE1-2. In addition, changes in the intensity of mixing (IE2) also cause corresponding variations in the distribution of dissipation. We also noticed that modulation of dissipation distribution is more apparent than that of conversion distribution.
It is interesting that except for IE2-2, the baroclinic tidal energy budget is close among each experiment. In IE2-2, the baroclinic tidal energy budget changes apparently. Previous studies found that variations of turbulent mixing can obviously affect stratification (Garrett, 2003;Drake et al., 2020;Delorme et al., 2021). Reinforcement of mixing can weaken the stratification or even destroy the structure of the pycnocline (Ivey, 2004). Changes of stratification may module the generation and dissipation of ITs. Figure 6A illustrates buoyancy frequency in the central of the LS (20.0 • N, 121.0 • E). Intensity of turbulent mixing greatly influences stratification. Stratification is weakened when mixing strengthened. When mixing raises to the magnitude of the observed maximum diffusivity, stratification will be destroyed or even absent, especially in the upper layer (above 400 m). Structure of the diapycnal diffusivity also plays an important role   Variations of stratification can result in the alteration of density perturbation, but only to a small extent for experiments, except IE2-2. The amplitude of density perturbation is influenced and the phase difference between density perturbation and vertical barotropic velocity is either enlarged or engendered slightly. However, in IE2-2, the amplitude of density perturbation is obviously weakened and nearly loses the temporally varying trend, which can be attributed to the absence of stratification. Therefore, the conversion in IE2-2 is impaired evidently. When turbulent mixing keeps increasing and becomes strong enough to destroy stratification, the variation trend of density perturbation will be destructed and eventually, the generation of ITs will be affected significantly. In other words, the variation of stratification only can influence the conversion of ITs to a small extent until the stratification is destroyed by excessive strong mixing.

Modulation on Incoherency of ITs
Turbulent mixing modulates the distribution of conversion and dissipation of ITs, which may influence local velocities. The baroclinic currents were calculated by removing the barotropic currents, which are represented by depth-averaged currents from the simulated currents. A fourth-order Butterworth filter was introduced to time series of baroclinic currents in diurnal (0.8-1.2 cpd) and semidiurnal (1.73-2.13 cpd) bands. Then, currents of diurnal and semidiurnal ITs can be obtained. Simulated results of the last 15 days are considered for analysis. To avoid errors from the band pass filter, results at the two ends (48 h for each end) are not included for investigation and only results in the 11 middle days are used. Figure 7 shows the eastward internal tidal currents in the northern SCS (20.0 • N, 115.0 • E) in diurnal and semidiurnal frequency bands. In the northern SCS, diurnal ITs are stronger than their semidiurnal counterparts. Distinctions in velocities exist among different experiments. For diurnal ITs, intensities and phases of velocity will change when different schemes of turbulent mixing are applied, especially for currents in the upper layer. Variations in phases of semidiurnal currents are also obvious and differences of semidiurnal velocities are evident in depths between 200 and 400 m. The variation of the current phase will make ITs in the SCS not totally phase locked with surface-tidal forcing at the LS, which may increase the incoherence.
Incoherency of ITs in the west of SCS (sites C and N in Figure 8) is calculated. Based on the least square method, harmonic analysis was applied to baroclinic currents to extract coherent ITs. Incoherent ITs were obtained by subtracting the coherent ITs from the band pass filtered baroclinic tidal currents. Vertically and temporally averaged incoherency of diurnal and semidiurnal ITs at sites C and N are shown in Table 3. In these sites, the incoherence varies among different experiments, indicating turbulent mixing can modulate the incoherency of ITs. Except IE2-2, incoherency is larger when the three-dimensional structure of the diapycnal diffusivity is employed for both the diurnal and semidiurnal ITs in the northern site (site C). The incoherency in IE1-2 can be as large as two times of those in other experiments. Meanwhile, for the southern site (site N), incoherency of ITs has no evident difference among idealized experiments, also except IE2-2. To investigate the difference in incoherency between sites C and N, density anomaly in each  experiment is calculated according to: where, ρ h is the mean density in each layer. Distributions of density anomaly at the end of simulation duration at 100 m depth are shown in Figure 8. The horizontal variations of density can lead to the incoherence of ITs (Li et al., 2020). Differences of density anomaly emerge when turbulent mixing changes, which further influence the incoherence of ITs. χ 2 test method is applied to investigate the differences in density anomaly between site C and site N. The χ 2 test results indicate that density anomaly has a significant difference (p < 0.01) at these two sites in IE1-2, IE2-1, and IE3-2. Meanwhile, in other experiments, differences of density anomaly are insignificant. Differences in incoherency between site C and site N are also more obvious in experiments with significant χ 2 test results, which indicate differences of density anomaly that lead to the corresponding variations in incoherency. We noticed that difference in incoherency of diurnal ITs between site C and site N in IE2-2 is large, but results of χ 2 test indicate that the difference is insignificance and the incoherency in IE2-2 is much larger than that in other experiments. In addition to the influence of horizontal variation of density, the evident incoherence of ITs in IE2-2 can also possibly be attributed to the absence of vertical stratification. The excessive strong mixing in IE2-2 destroys the stratification (Figure 6A), whereas damage of stratification could possibly influence the propagation of ITs, which may change the amplitude and phase of ITs and eventually add incoherence to ITs. Incoherency of ITs can be influenced by background currents and mesoscale eddies. Moreover, background currents and mesoscale eddies are thought as essential factors, which modulate the incoherent feature of ITs. The mesoscale eddies can increase the intensity of background currents and induce horizontal variations of density, which would eventually add incoherence to ITs (Li et al., 2020). But, the relationship between incoherency of ITs and turbulent mixing has rarely been discussed. Our findings suggest that turbulent mixing can also increase incoherence of ITs (site C in IE2-1) by changing horizontal density in the ocean. Figure 9 shows eastward velocities of diurnal ITs, coherent signals, and incoherent signals at the depth of 100 m at site C. Figure 10 demonstrates the results of semidiurnal ITs. The variation in horizontal density modulates incoherency of both the diurnal and semidiurnal ITs largely by changing the internal tidal amplitudes. In addition, phases of ITs are also slightly affected.

Modulation of Nonlinear Interactions Among Different Tidal Components
Previous studies found that both the terdiurnal and quarterdiurnal constituents can be generated from nonlinear interactions of ITs in the SCS (Xie et al., 2010). Xie et al. (2010) found that the intensity of M 2 plays important role in those nonlinear interactions. Currents of higher harmonic components of ITs are extracted by applying band pass filter (with cutoff frequency band between 2.8 and 4.2 cpd, which contains principle overtides of MK 3 , MO 3 , and M 4 ) to time series of simulated baroclinic currents. The kinetic energy (KE) of overtides for per unit volume is calculated according to: where, ρ 0 = 1,025 g·cm −3 is the mean density of the sea water and u t and v t are filtered eastward and northward components of overtides. Temporally averaged KE at sites A to E ( Figure 11A) is shown in Figures 11B-F. At all these sites, differences in KE among experiments in which the diapycnal diffusivity is set to be homogeneous are larger than those of experiments in which varying spatial structure is applied. This suggests that the variation in intensity of turbulent mixing has more effect on the interaction of ITs. KE in IE2-2 is nearly the smallest among those sites and KE in IE2-1 is larger at most depths, indicating that  Frontiers in Marine Science | www.frontiersin.org the elevation of turbulent mixing works against the nonlinear interactions among ITs and the reduction of mixing intensity is in favor of the nonlinear interactions. It is interesting to note that overtides in the LS are stronger than those in the SCS. During the propagation from sites A to C, energy of overtide decreases at most depths. This suggests that strength of overtide decays during propagation. Based on observations, Xu et al. (2013Xu et al. ( , 2014 also found that quarterdiurnal ITs on the continental slope in the northwestern SCS are weaker than those in the deep basin on the west of the LS. Empirical orthogonal function (EOF) analysis is applied to baroclinic currents of overtides in IE1-1. The leading EOF modes (explaining more than 70% of total variance) are shown in Figure 12. At site A, in the eastward direction, the first and the second EOF modes can explain 63 and 19% of the total variance, respectively, and the first and the second EOF modes correspond to the first and the third baroclinic modes, respectively. In the northward direction, mode I, which corresponds to the first baroclinic mode, contributes more than 70% of the total variance. At sites D and E, which are also located near the LS, more than 70% of the total variance of currents in the eastward direction need to be explained by two EOF modes, which correspond to the second and the first baroclinic modes, respectively. Modes in the northward direction are more complicated. At site D, the top two EOF modes correspond to the first and the second baroclinic modes, respectively. But, the first and the second EOF modes in site E are equal to the third and the second baroclinic modes, respectively. Meanwhile, at the sites in the SCS (B and C), the first EOF mode, which is equivalent to the first baroclinic mode, can contribute more than 70% of the total variance. Generally, sites near the LS have higher modes than those in the SCS. The overtides are generated at the LS and have higher vertical modes, which are more likely to dissipate during propagation. Therefore, the intensity of overtides may decrease significantly during the propagation from the LS to the SCS.

SUMMARY AND DISCUSSION
In this study, modulations of ITs on the baroclinic tidal energy budget, the incoherency, and the nonlinear interactions among different tidal components by turbulent mixing are investigated through numerical simulations. The spatial structure of turbulent mixing can hardly influence the baroclinic tidal energy budget. In addition, change in the mixing intensity in a reasonable range also cannot obviously affect the total baroclinic energy converted or dissipated in the LS and the SCS. However, excessively strong mixing will destroy the vertical stratification, the variation trend of density perturbation will be destructed, and eventually, the generation of ITs will be affected significantly. Meanwhile, the dissipation rate will also change apparently, when turbulent mixing keeps increasing and becomes strong enough to destroy stratification. Comparing to the baroclinic tidal energy budget, the influence on distributions of conversion and dissipation is more evident when turbulent mixing changes. Distinction of conversion distribution is largely concentrated near the SCS, while alteration of dissipation distribution is more apparent in the LS. Moreover, modulation on the distribution of dissipation is more apparent than that of conversion.
The relationship between incoherency of ITs and turbulent mixing is also explored. The findings suggest that in addition to background currents and mesoscale eddies, turbulent mixing can also increase incoherence of ITs. Variations of turbulent FIGURE 12 | Empirical orthogonal function (EOF) modes of eastward velocity of overtides at (A-E) sites A to E and corresponding results of (F-J) the northward velocity.
mixing change horizontal density in the ocean. The horizontal variations of density add incoherence to ITs largely by changing the internal tidal amplitudes. Moreover, variation of mixing intensity modulates the nonlinear interactions among different components of ITs. Elevation of mixing intensity works against the nonlinear interactions and the reduction of mixing intensity is in favor of the nonlinear interactions. However, the structure of turbulent mixing can hardly influence the interaction among ITs. We also find that overtides decrease evidently during the propagation from the LS to the SCS. For the overtides, they are generated at the LS and have higher vertical modes, which are more likely to dissipate during propagation.
It is notable that the baroclinic tidal energy budget can hardly be affected by mixing structure and the spatial structure of turbulent mixing also has little effect on the nonlinear interaction among ITs. The nonlinear interactions are modulated by the mixing intensity. Therefore, the influence of mixing structure can be neglectable in most of the internal tidal simulations. and the diapycnal diffusivity can set to be horizontally and vertically homogeneous. Vertical stratification will be destroyed when turbulent mixing is excessively strong. The damage of stratification greatly influences the conversion and dissipation of ITs. Therefore, excessive strong mixing needs to be avoided. Considering that the observed minimum and mean value of the diapycnal diffusivity in the SCS are 9.3 × 10 −6 m 2 s −1 and 1.7 × 10 −3 m 2 s −1 , which are close to O (10 −5 ) m 2 s −1 and O (10 −3 ) m 2 s −1 , respectively, the optimum range of the diapycnal diffusivity in the SCS is more appropriate to be from O (10 −5 ) to O (10 −3 ) m 2 s −1 for IT simulations. However, the spatial structure of turbulent mixing can modulate the incoherency of ITs. Consequently, the spatial structure of the diapycnal diffusivity should be considered when the incoherency of ITs needs to be simulated.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
LD and XL designed the research. BL and SP performed the research. BL wrote the manuscript. XM and YY analyzed the data. All authors contributed to the article and approved the submitted version.