Experimental Investigation on the Vertical Structure Characteristics of Internal Solitary Waves

Abstract

An experimental investigation of the vertical structure characteristics of internal solitary waves (ISWs) was systematically carried out in a large gravitationally stratified fluid flume. Four different stratifications were established, and basic elements of ISWs were measured by a conductivity probe array. The vertical distributions of the amplitude, characteristic frequency and waveform of two types of ISWs under different stratifications were obtained, and the experimental results were compared with the theoretical model. The study shows that most vertical structures of the amplitude under different stratifications agree with those of the theoretical model, while there are some deviations for ISWs with large amplitudes. Neither the two-layer model nor the continuously stratified model can effectively describe the variation in the characteristic frequency at different depths with amplitude. For a single small-amplitude ISW, the characteristic frequency first increases and then decreases with increasing depth. The characteristic frequency is largest at the depth of the maximum buoyancy frequency. For an ISW with a relatively large amplitude, there is likely to be a local minimum of the characteristic frequency near the depth where the maximum buoyancy frequency lies. In different stratifications, the ${\mathrm{sech}}^2$ function of KdV theory can describe the waveforms of ISWs at different depths well.

1. Introduction

Internal solitary waves (ISWs), which are a kind of nonlinear internal wave, are ubiquitous in marginal sea areas [1,2,3]. Based on the balance between nonlinearity and dispersion, ISWs can maintain a stable shape for long distances in the horizontal direction [4]. Due to the opposite flow direction induced by the upper and lower layers of the wave, a strong shear flow can be vertically induced, which can have a significant impact on the safety of marine structures and underwater platforms [5,6,7,8,9,10]. In addition, ISWs with large amplitudes, which are often accompanied by relatively large vertical velocities, can have an important impact on seawater mixing [11], nutrient transport [12] and sediment resuspension [13] in the ocean. Therefore, studying the vertical structure characteristics of ISWs is greatly significant.

Vertical structures of ISWs usually refer to the distribution of amplitudes at different depths [14]. A large amount of information can be obtained from the vertical structure of ISWs, including the depth of the maximum amplitude and the vertical scale of ISWs. Additionally, vertical structures can reflect the evolution of ISWs to a certain extent [15]. Many theoretical studies on vertical structures of ISWs have been carried out thus far [16,17,18,19,20,21,22]. Fliegel and Hunkins [18] first calculated the dispersion relation of internal waves by using the Thomson–Haskell method and obtained vertical structures of internal waves in different modes. Based on the Thomson–Haskell method, Cai and Gan [19] calculated the internal wave dispersion relation without using the Boussinesq approximation and theoretically analyzed the vertical structure characteristics of different modes at two stations in the South China Sea. Zhang and Gao [20] obtained vertical structures and the dispersion relation of internal waves by using the transformation method of Russian scholars. Based on a numerical method, Deng et al. [21] analyzed the vertical structures of different modes of internal waves in background currents. Kurkina et al. [22] solved the theoretical vertical structures of first-mode and second-mode ISWs with nonlinear corrections in two-layer and three-layer fluids.

The propagation and evolution characteristics of ISWs have been investigated in several studies by using in situ observational data [23,24,25,26,27]. Vlasenko et al. [28] obtained measured data of ISWs in the Strait of Messina and found that the observed vertical structures of ISWs are consistent with numerical simulation results. Geng et al. [29] obtained the vertical structures of ISWs near the Dongsha Atoll of the South China Sea using seismic technology and found that the observational vertical structures of ISWs are generally consistent with the theoretical results except when ISWs interact with the topography. Similarly, Gong et al. [14] compared the vertical structures of ISWs obtained by the seismic method with two theories. The results show that ISWs with low nonlinearity can be described by a linear vertical structure function, while ISWs with high nonlinearity are in agreement with the first-order nonlinear vertical structure function. Because obtaining refined vertical structure data from field observations is difficult, the verification of the theoretical vertical structure of ISWs based on observational data is relatively limited. In terms of laboratory experiments, many experimental methods have been used to study the characteristics of ISWs [30,31,32,33,34,35]. At present, experimental research objects of ISWs are not limited to the first mode. Many experimental studies on the generation mechanism and evolution characteristics of second-mode ISWs have been carried out [36,37,38]. However, most experimental studies focus on the propagation and evolution characteristics of ISWs at the interface (position of maximum amplitude), while there are few studies on the vertical structure of ISWs in stratified fluids.

Korteweg de Vries’ (KdV) theory is the most commonly used theoretical model to describe the propagation and evolution characteristics of ISWs in laboratory experiments [39,40,41]. Since weak nonlinearity is an important condition for validating the KdV model, KdV theory is not applicable to describe ISWs with large amplitudes [42,43,44]. Michallet and Barthelemy [45] compared the theoretical and experimental results of the waveform, phase velocity and frequency–amplitude relationship of ISWs and found that KdV theory can describe small-amplitude solitary waves in different stratified fluids well. Grue et al. [46] analyzed the propagation characteristics of ISWs in two-layer fluids and found that KdV theory can provide predictions that are consistent with the experimental results of small-amplitude ISWs. Du et al. [47] compared the experimental results of the elevation- and depression-type ISWs with the results from the following three different theoretical models: KdV, extended KdV (eKdV) and Miyata–Choi–Camassa (MCC). The results show that the experimental results conform to those of KdV theory when the ratio of the amplitude of ISWs to the total water depth is less than 0.04. However, KdV theory is presently used to describe the characteristics of ISWs at the interface in laboratory experiments, and the validity of KdV theory describing the properties of ISWs at different depths in stratified fluids is not clear.

In view of the above discussion, an experimental method is used in this paper to study the vertical structure of the amplitude of ISWs in a large, gravitationally stratified fluid flume and to analyze the domain of validity of KdV theory by comparing the characteristics of ISWs at different depths with theoretical models. The influence of environmental factors on the vertical structure characteristics of ISWs is further analyzed by setting different stratifications and amplitudes. This paper is organized in the following manner: The KdV model is introduced in Section 2, including the two-layer model and the continuously stratified model. The experimental techniques and methods are presented in Section 3. The results and discussion are provided in Section 4, and the conclusion is presented in Section 5.

2. Theoretical Models

Considering a stratified fluid of a finite depth with rigid boundaries where buoyancy frequency is a continuous function of the depth, the eigenfunction $W$ is determined by the solution of the boundary value problem:

(1)$\frac{d^{2}W}{d{Z}^2}+\frac{N^2(z)}{c^2}W=0,W(0)=W(-H)=0.$

where $N(z)$ is the buoyancy frequency, which is determined by the formula $N(z)={(-gd\rho /\rho dz)}^{1/2}$; $c$ is the linear phase speed of ISWs; and $H$ is the total water depth. Equation (1) can be solved by using the Thomson–Haskell method [18], and $c$ can be obtained. Then, $W$ can be calculated under the normalization condition of $W_{\max }=1$.

For ISWs with a maximum amplitude ${\eta }_{\max }$ in a stratified fluid, their amplitude at different depths ${\eta }_0(z)$ can be calculated by the following expression (with the nonlinear terms removed):

(2)${\eta }_0(z)={\eta }_{\max }W(z).$

Under the Boussinesq approximation, the vertical isopycnal displacement $\eta (x,z,t)$ of ISWs propagating in the positive direction of axis $x$ can be described by the classical KdV equations:

(3)$\frac{\partial \eta }{\partial t}+c\frac{\partial \eta }{\partial x}+\alpha \eta \frac{\partial \eta }{\partial x}+\beta \frac{{\partial }^3\eta }{\partial x^3}=0.$

where $t$ is the time, and $x$ is the horizontal coordinate along the propagation direction of ISWs. The parameters $\alpha$ and $\beta$ are called the nonlinear coefficient and the dispersion coefficient, respectively.

The well-known solitary wave solution to Equation (3) is

(4)$\eta (x,t)={\eta }_0 {\mathrm{sech}}^2(\frac{x-Vt}{\lambda }),$

where ${\eta }_0$ is the amplitude of ISWs, and the nonlinear velocity $V$ and the characteristic width $\lambda$ of this soliton can be calculated by:

(5)$V=c+\frac{\alpha {\eta }_0}{3},\lambda ={(\frac{12\beta }{\alpha {\eta }_0})}^{\frac{1}{2}}.$

The characteristic frequency of the soliton can be obtained by $\omega =V/\lambda$. Since $V$ and $\lambda$ are both relative to ${\eta }_0$, the characteristic frequency $\omega$ will change with ${\eta }_0$. Moreover, the obtained frequency–amplitude relationship varies for the two-layer fluid model and the continuously stratified fluid model due to the different calculation methods of $\alpha$, $\beta$ and $c$ in these two models.

In the continuously stratified model, $c_1$ can also be obtained when the eigenfunction $W$ is calculated by the Thomson–Haskell method. The terms ${\alpha }_1$, ${\beta }_1$ can be determined by the following formula:

(6)${\alpha }_1=\frac{3c_1}{2}{\displaystyle {\int }_{-H}^0{(\frac{dW}{dZ})}^{3}dz/}{\displaystyle {\int }_{-H}^0{(\frac{dW}{dZ})}^{2}dz},$

(7)${\beta }_1=\frac{c_1}{2}{\displaystyle {\int }_{-H}^0{W}^{2}dz/}{\displaystyle {\int }_{-H}^0{(\frac{dW}{dZ})}^{2}dz}.$

In the two-layer model, $c_2$, ${\alpha }_2$ and ${\beta }_2$ can be determined by the following expressions:

(8)$c_2={(\frac{2g({\rho }_2-{\rho }_1)h_1{h}_2}{({\rho }_1+{\rho }_2)(h_1+h_2)})}^{\frac{1}{2}},$

(9)${\alpha }_2=\frac{3c_2(h_1-h_2)}{2h_1{h}_2},$

(10)${\beta }_2=\frac{c_2{h}_1{h}_2}{6}.$

where $h_1$ and $h_2$ are the thicknesses of the upper and lower layers, respectively, and ${\rho }_1$ and ${\rho }_2$ are the densities of the upper and lower layers, respectively.

3. Experimental Device and Method

The experiments were performed in a gravitationally stratified fluid flume that is 1200 cm long, 120 cm wide and 100 cm high. At the right end of the flume, the wave maker aimed to achieve gravity collapse in the stratified fluid, which is presented in Figure 1. The step length and depth are $L_0$ and $D_0$, respectively, in the experimental generation principle of ISWs. This device can produce high-quality elevation- or depression-type ISWs at arbitrary pycnocline depths. At the left end of the flume, a wedge-shaped wave absorber was installed. According to the amplitude of ISWs and interface depth, the wedge angle can be adjusted to effectively eliminate the reflected waves. A stratified fluid was created by using the classic “double-tank” method. At first, the flume was filled with fresh water of density ${\rho }_1$ to the depth $h_1$. Then a brine solution of density ${\rho }_2$ was slowly injected beneath the fresh water through two small openings along the bottom of the flume until the brine depth reached $h_2$, and the total depth was $H=h_1+h_2$. Due to the mixing of fresh water in the upper layer and brine solution in the lower layer, a continuous density distribution structure with a certain thickness was formed at the interface, namely the pycnocline, which can simulate the actual vertical density distribution of the ocean with the characteristics of a typical pycnocline.

The vertical distribution of density in stratified fluids and vertical structures of ISWs was measured by an internal wave dynamic measurement system. This system consists of a conductivity probe array, signal control circuit and data processing system. A vertically displaceable probe was used to measure the vertical density distribution of the fluid. A conductivity probe array consisting of seven probes vertically arranged at equal intervals was used to measure the vertical structure of ISWs. The horizontal distance between the probes and the wave maker is 6 m, so the fully developed ISWs can be detected (Figure 1). Adjusting the position of the conductivity probe array to cover the whole pycnocline according to the stratification conditions and amplitude of ISWs, the fluctuation of ISWs at different depths can be completely captured. The origin of the axis $z$ is set at the surface with the positive axis $z$ running downward. Figure 2a shows the vertical distribution of isopycnic lines of ISWs in a certain stratification. The back portion of the wave was contaminated by the small-amplitude disturbances in the tail. Therefore, as shown in Figure 2b, only the forward portion of the wave was used to calculate the amplitude and characteristic length of ISWs [40,41]. The characteristic length and frequency of ISWs can be experimentally calculated by the following expressions (Equations (10) and (11), respectively):

(11)$\lambda =\frac{V}{{\eta }_0}{\displaystyle {\int }_{-\infty }^0\eta (t,x_0)dt}$

(12)$\omega =\frac{V}{\lambda }=\frac{{\eta }_0}{{\displaystyle {\int }_{-\infty }^0\eta (t,x_0)dt}}$

4. Results and Discussion

To fully study the vertical structure characteristics of ISWs, four stratification conditions were set in the experiment, and two types of ISWs (elevation- and depression-type) were generated during operation. The specific experimental conditions are summarized in Table 1. S1 is a representative of the stratification in the deep sea. S2–S4 are used to simulate the different stratifications in the continental shelves. Both elevation-type ($\alpha >0$) and depression-type ($\alpha <0$) ISWs exist in real oceans. Compared with the depression-type ISWs, the amplitudes of the elevation-type ISWs are generally smaller. During the experiment, the total water depth $H=h_1+h_2=80\mathrm{cm}$ and the densities of upper and lower layers (${\rho }_1$ and ${\rho }_2$, respectively) remained unchanged. Additionally, the dimensionless maximum amplitude $\eta \prime =\left|{\eta }_{\max }/H\right|$ ranged from 0.008 to 0.120.

4.1. Vertical Structure of the Amplitude

According to the theoretical model, the distribution of the buoyancy frequency $N(z)$ can be obtained from the experimental density profile $\rho (z)$, and the vertical modal function of the internal wave $W(z)$ can be subsequently calculated. Figure 3 illustrates the vertical distribution of ISWs in S3, with Figure 3a showing the density profile and buoyancy frequency distribution and Figure 3b presenting the vertical modal function of ISWs obtained by the Thomson–Haskell method. Figure 3 shows that under the condition of this density stratification, the variation in the buoyancy frequency and internal wave vertical modal function with depth are consistent, which show a trend of first increasing and then decreasing with an increasing depth. The maximum value of the buoyancy frequency profile is at a depth of $z/H=0.255$, and the maximum value of the corresponding vertical modal function is at a depth of approximately $z/H=0.26$ (red line in Figure 3b). There is a small deviation between the depth of the maximum value of the buoyancy frequency profile and the depth of the maximum value of the vertical modal function that is obtained by the Thomson–Haskell method, which can also be found from the theoretical results of Kurkina et al. [22].

Figure 4 demonstrates a comparison between the theoretical and experimental results of the vertical structure of the amplitude of ISWs in four stratification conditions. Under different stratification conditions, whether elevation-type or depression-type, the theoretical vertical structure of the amplitude is basically consistent with most of the experimental results. In addition, the amplitude initially increases and then decreases with an increasing depth both theoretically and experimentally. For S1, the depths of the maximum amplitudes of ISWs under different cases are basically the same, while the depth of the maximum amplitude of ISW with large amplitude (${\eta }^{\prime }=0.120$) moves up. For S2, most of the vertical distributions of the amplitude conform to the theory, while the vertical amplitude structure of ISWs with relatively large amplitudes (${\eta }^{\prime }=0.108$ and ${\eta }^{\prime }=0.110$) moves up. For S3, the depths of the maximum amplitudes under different cases and the vertical structure of the amplitude are consistent with the theory. For S4, similar to S3, the theory can describe the vertical structure of ISWs under different cases well. Figure 5 shows the variation in root mean square error $E$ with increasing ${\eta }^{\prime }$ under different stratification conditions. Root mean square errors in different stratifications are relatively small when ${\eta }^{\prime }<0.1$. However, once ${\eta }^{\prime }>0.1$, the root mean square errors start to increase significantly and exceed the allowable range.

For S1 and S2, the vertical structures of ISWs basically conform to the theory in most cases, while for ISWs with large amplitudes (${\eta }^{\prime }>0.1$), the experimental vertical structures deviate from the theoretical model. These results occur because the effect of the nonlinear terms will be more significant with an increasing amplitude. In addition, with an increasing $h_1/H$, the difference between the theoretical and experimental vertical structures of ISWs with large amplitudes is more obvious. However, for S3 and S4, the vertical structures of ISWs are consistent with the theoretical value due to the relatively small amplitudes (${\eta }^{\prime }<0.1$).

4.2. Vertical Distribution of the Characteristic Frequency

As shown in Figure 6, the variation in the theoretical characteristic frequency of ISWs with amplitudes in continuously stratified fluid and two-layer fluid can be obtained according to the expression of the characteristic frequency $\omega =V/\lambda$. The results show that the characteristic frequency of ISWs increases with an increasing amplitude whether it is for the continuously stratified fluid or the two-layer fluid. The theoretical value of the characteristic frequency in the two-layer fluid is larger than that in the continuously stratified fluid at the same amplitude. Generally, this is mainly due to smaller absolute values of $\alpha$ and larger values of $\beta$ in the continuously stratified model for different stratifications. In addition, the variation in characteristic frequency is closely related to the stratification conditions. For the depression-type ISWs, both the characteristic frequency under the same amplitude and the difference in the theoretical values between the two-layer model and the continuously stratified model will correspondingly decrease with an increasing $h_1/H$. The variation in the characteristic frequency of the elevation-type ISWs with the amplitude in S4 is basically consistent with that in S3. These trends are also closely linked to values of $\alpha$ and $\beta$ in different stratification conditions.

The theoretical and experimental results of the frequency–amplitude relationship of ISWs in S1–S4 are presented in Figure 7. For S1, the experimental characteristic frequencies at small amplitudes (${\eta }_0/H<0.02$) are slightly larger than their theoretical values. The characteristic frequencies of ISWs are close to the theoretical values of the two-layer model under the condition of $0.02\le {\eta }_0/H\le 0.04$. When $0.04<{\eta }_0/H\le 0.08$ is satisfied, the experimental characteristic frequencies approach the theoretical values of the continuously stratified model. The characteristic frequencies of ISWs are smaller than the theoretical frequencies for relatively large amplitudes (${\eta }_0/H>0.08$). For S2, which is similar to S1, the experimental characteristic frequencies at small amplitudes (${\eta }_0/H<0.04$) are larger than their theoretical values. Regarding the amplitude that satisfies $0.04\le {\eta }_0/H\le 0.09$, the characteristic frequencies of ISWs are basically between the theoretical values of the two models but closer to the continuously stratified model. The characteristic frequencies of ISWs are still smaller than the theoretical frequencies for relatively large amplitudes (${\eta }_0/H>0.09$). For S3, the theoretical values of the two models are coinciding, but the characteristic frequencies of most ISWs obtained by experiments are larger. For S4, the theoretical values of the two models are also close, and the characteristic frequencies obtained experimentally are dispersed with the amplitude. Different from other stratification conditions, the experimental values of the characteristic frequencies can be larger or smaller than the theoretical ones under the same amplitude.

Comparing different stratification conditions, the characteristic frequency of ISWs increases with an increasing amplitude for the four stratifications, which is consistent with the trend of the theoretical model. Among these stratifications, S1 and S2 can describe the variation in the characteristic frequency with the amplitude to some extent by the two-layer model or the continuously stratified model. For S3 and S4, these two models have difficulty describing the frequency–amplitude relationship. For the depression-type ISWs (S1–S3), the characteristic frequency distribution is more concentrated with the amplitude, while the characteristic frequency distribution of elevation-type ISWs (S4) is relatively dispersed with the amplitude. One possible reason may be that the characteristic frequencies do not change much with the amplitude due to the interactive relationship of the waveforms at different depths. Therefore, the experimental characteristic frequency does not increase with an increasing amplitude as fast as that in models. In addition, we also consider that the waveforms of elevation-type ISWs can be easily affected by some disturbances due to the small amplitudes, which results in the dispersive amplitude–frequency relationship in S4. Another explanation for the discrepancy is that the density profile of the fluid may change with the extension of time, leading to the difference between the experimental characteristic frequency and the theoretical values.

Figure 8 shows the vertical variation in the characteristic frequency under different stratifications and cases, where the solid line is the fitted curve reflecting the variation in the characteristic frequency with depth. Figure 8a–c show that the variation in the characteristic frequency with depth can be divided into two stages. In the first stage, the characteristic frequency increases with an increasing depth and then begins to decrease when it reaches a certain value. Then, there is a local minimum of the characteristic frequency. In the second stage, the characteristic frequency initially increases with an increasing depth and then decreases when it reaches a certain value. Some studies have shown that the characteristic frequency initially increases with an increasing amplitude. When the amplitude reaches a certain value, the characteristic frequency will stop increasing and begin to decrease [46]. Since the amplitude first increases and then decreases with an increasing depth, such variation in the characteristic frequency with the amplitude is formed. The red line in Figure 7 represents the depth where the maximum amplitude of ISWs lies (position of the maximum buoyancy frequency in the stratified fluid), which corresponds to the depth where the local minimum of the characteristic frequency lies. Figure 8d shows that the characteristic frequency first increases and then decreases with an increasing depth under different cases in S4, and the characteristic frequency reaches the maximum near the depth of the maximum buoyancy frequency. This is mainly because the amplitudes in this stratification are relatively small. Therefore, the characteristic frequency only increases with an increasing amplitude in S4. In summary, for a single small-amplitude ISW, the characteristic frequency initially increases and then decreases with an increasing depth. The characteristic frequency is the largest at the depth where the maximum buoyancy frequency lies. For an ISW with relatively large amplitude, there is likely to be a local minimum of the characteristic frequency near the depth where the maximum buoyancy frequency is located.

4.3. Vertical Distribution of the Waveform

The ISW solution obtained by KdV theory is $\eta (x,t)={\eta }_0 {\mathrm{sech}}^2[(x-Vt)/\lambda ]$. This ISW solution contains three basic elements (amplitude, characteristic frequency and ${\mathrm{sech}}^2$ function). The amplitude and characteristic frequency of ISWs have been discussed above. To investigate whether the ${\mathrm{sech}}^2$ function in KdV theory can describe waveforms of ISWs at different vertical depths, the amplitude and characteristic frequency obtained by experiments are used to calculate the theoretical waveform distribution of ISWs. Then, the theoretical waveform distribution is compared with the experimental results to verify the validity of the ${\mathrm{sech}}^2$ function. Figure 9 illustrates the theoretical and experimental results of the vertical distribution of waveforms in four different stratification conditions. Figure 9 shows that in four different stratification conditions, the experimental distribution of solitary wave waveforms at different depths is consistent with the theory. Since some errors may occur when the conductivity probes are used to measure the density in the flume, the experimental waveforms of the ISWs at some depths deviate from the theory to some extent. Generally, the deviations are within the allowable range. Thus, the ${\mathrm{sech}}^2$ function in KdV theory can describe the waveforms of ISWs at different depths well.

5. Conclusions

In this paper, four stratifications are established by using a gravitationally stratified fluid flume. The experimental results of the vertical structure characteristics of ISWs in each stratification are compared with the theoretical model, and the influence of different stratifications and different maximum amplitudes on the agreement between the theoretical and experimental results is analyzed. The conclusions are as follows:

Under different stratification conditions, the vertical amplitude structures of most ISWs with relatively small amplitudes are in good agreement with the theoretical model, whether it is elevation-type or depression-type. Additionally, with the increase in the maximum amplitude of ISWs, the influence of the nonlinear terms of the vertical structure function will be more significant, which leads to deviations between the experimental results and the theoretical model.

For the same stratification and the same wave amplitude, the theoretical value of the characteristic frequency in a two-layer fluid is larger than that in a continuously stratified fluid. For the depression-type ISWs, with the increase in $h_1/H$, the theoretical values of the characteristic frequency in the two models are both larger, and the difference between the theoretical values of the two models is more significant under the same amplitude. The frequency–amplitude relationship in experiments indicates that the frequency–amplitude relationship can be described to a certain extent by the two-layer model or the continuously stratified model only when $h_1/H$ is relatively small.

The vertical distribution of the characteristic frequency of a single ISW in the experiment indicates that the characteristic frequency of a single small-amplitude ISW initially increases and then decreases with an increasing depth. The characteristic frequency is largest at the depth where the maximum buoyancy frequency lies. For an ISW with a relatively large amplitude, there is likely to be a local minimum of the characteristic frequency near the depth where the maximum buoyancy frequency is located.

Under different stratification conditions, the ${\mathrm{sech}}^2$ function of KdV theory can be in good agreement with the waveforms of ISWs at different depths regardless of whether the ISWs are depression-type or elevation-type.

Author Contributions

Conceptualization, P.P. and H.D.; methodology, P.P. and H.D.; software, H.D. and G.W.; validation, P.P., H.D., S.W. and P.X.; formal analysis, P.P.; investigation, P.P., H.D. and P.X.; resources, H.D., S.C. and J.X.; data curation, H.D. and S.W.; writing—original draft preparation, P.P.; writing—review and editing, H.D., S.C. and J.X.; visualization, P.P.; supervision, H.D., G.W., S.C. and J.X.; project administration, H.D. and G.W.; funding acquisition, S.W., H.D. and G.W. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank the reviewers for their valuable comments, which helped to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Table

Figure 1. Schematic diagram of the experimental apparatus.

Figure 2. Vertical distribution of isopycnic lines of ISWs in a certain stratification (a) and experimental method of calculating the amplitude and characteristic frequency of ISWs (b).

Figure 3. Profile of the density and buoyancy frequency in S3 (a) and the vertical modal function calculated by the Thomson–Haskell method in S3 (b).

Figure 4. Experimental vertical structure of the amplitude and theoretical vertical mode function: (a) S1; (b) S2; (c) S3; (d) S4.

Figure 5. Variation in [Forumla omitted. See PDF.] with [Forumla omitted. See PDF.] in different stratifications.

Figure 6. The theoretical frequency–amplitude relationship in S1–S4, where CS refers to continuously stratified.

Figure 7. Experimental and theoretical frequency–amplitude relationship of ISWs: (a) S1; (b) S2; (c) S3; (d) S4.

Figure 8. Experimental vertical distributions of the characteristic frequencies: (a) S1; (b) S2; (c) S3; (d) S4. The red line represents the depth of the maximum amplitude of ISWs.

Figure 9. Experimental and theoretical waveforms of ISWs at different depths: (a) S1; (b) S2; (c) S3; (d) S4.

Table 1

Summary of experimental conditions.

Stratifications	Type 1	${\mathit{h}}_1$ (cm)	${\mathit{h}}_2$ (cm)	${\mathit{\rho }}_1$ (g/cm3)	${\mathit{\rho }}_2$ (g/cm3)	Probes 2 (cm)	${\mathit{\eta }}^{\prime }$
S1	D	10	70	0.998	1.018	5–23	0.075–0.120
S2	D	15	65	0.998	1.018	10–28	0.017–0.110
S3	D	20	60	0.998	1.018	15–33	0.021–0.095
S4	E	60	20	0.998	1.018	52–64	0.008–0.068

¹ Type refers to the type of ISWs. D refers to depression-type, E refers to elevation-type. ² Probes refers to the range of the depth of probes.

References

1. Cai, S.; Xie, J.; He, J. An overview of internal solitary waves in the South China sea. Surv. Geophys.; 2012; 33, pp. 927-943. [DOI: https://dx.doi.org/10.1007/s10712-012-9176-0]

2. Cai, S.; Xie, J.; Xu, J.; Wang, D.; Chen, Z.; Deng, X.; Long, X. Monthly variation of some parameters about internal solitary waves in the South China sea. Deep-Sea. Res. Part I; 2014; 84, pp. 73-85. [DOI: https://dx.doi.org/10.1016/j.dsr.2013.10.008]

3. Alford, M.H.; Peacock, T.; MacKinnon, J.A.; Nash, J.D.; Buijsman, M.C.; Centurioni, L.R.; Chao, S.-Y.; Chang, M.-H.; Farmer, D.M.; Tang, T.Y. The formation and fate of internal waves in the South China Sea. Nature; 2015; 521, pp. 65-69. [DOI: https://dx.doi.org/10.1038/nature14399] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25951285]

4. Grimshaw, R.; Pelinovsky, E.; Talipova, T.; Kurkin, A. Simulation of the Transformation of Internal Solitary Waves on Oceanic Shelves. J. Phys. Oceanogr.; 2004; 34, pp. 2774-2791. [DOI: https://dx.doi.org/10.1175/JPO2652.1]

5. Cai, S.; Wang, S.; Long, X. A simple estimation of the force exerted by internal solitons on cylindrical piles. Ocean Eng.; 2006; 33, pp. 974-980. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2005.05.012]

6. Cai, S.; Long, X.; Wang, S. Forces and torques exerted by internal solitons in shear flows on cylindrical piles. Appl. Ocean Res.; 2008; 30, pp. 72-77. [DOI: https://dx.doi.org/10.1016/j.apor.2008.03.001]

7. Guo, H.; Zhang, L.; Li, X.; Lou, M. Dynamic responses of top tensioned riser under combined excitation of internal solitary wave, surface wave and vessel motion. J. Ocean Univ. China; 2013; 12, pp. 6-12. [DOI: https://dx.doi.org/10.1007/s11802-013-2079-y]

8. Fan, H.; Li, C.; Wang, Z.; Xu, L.; Wang, Y.; Feng, X. Dynamic analysis of a hang-off drilling riser considering internal solitary wave and vessel motion. J. Nat. Gas. Sci. Eng.; 2017; 37, pp. 512-522. [DOI: https://dx.doi.org/10.1016/j.jngse.2016.12.003]

9. Gong, Y.; Xie, J.; Xu, J.; Chen, Z.; He, Y.; Cai, S. Oceanic internal solitary waves at the Indonesian submarine wreckage site. Acta Oceanol. Sin.; 2022; 41, pp. 109-113. [DOI: https://dx.doi.org/10.1007/s13131-021-1893-0]

10. Wang, T.; Huang, X.; Zhao, W.; Zheng, S.; Yang, Y.; Tian, J. Internal Solitary Wave Activities near the Indonesian Submarine Wreck Site Inferred from Satellite Images. J. Mar. Sci. Eng.; 2022; 10, 197. [DOI: https://dx.doi.org/10.3390/jmse10020197]

11. La Forgia, G.; Cavaliere, D.; Adduce, C.; Falcini, F. Mixing Efficiency for Breaking Internal Solitary Waves. J. Geophys. Res. Oceans; 2021; 126, e2021JC017275. [DOI: https://dx.doi.org/10.1029/2021JC017275]

12. Woodson, C.B. The Fate and Impact of Internal Waves in Nearshore Ecosystems. Ann. Rev. Mar. Sci.; 2018; 10, pp. 421-441. [DOI: https://dx.doi.org/10.1146/annurev-marine-121916-063619]

13. Boegman, L.; Stastna, M. Sediment Resuspension and Transport by Internal Solitary Waves. Annu. Rev. Fluid Mech.; 2019; 51, pp. 129-154. [DOI: https://dx.doi.org/10.1146/annurev-fluid-122316-045049]

14. Gong, Y.; Song, H.; Zhao, Z.; Guan, Y.; Kuang, Y. On the vertical structure of internal solitary waves in the northeastern South China Sea. Deep-Sea. Res. Part I; 2021; 173, 103550. [DOI: https://dx.doi.org/10.1016/j.dsr.2021.103550]

15. Vlasenko, V.; Hutter, K. Numerical Experiments on the Breaking of Solitary Internal Waves over a Slope–Shelf Topography. J. Phys. Oceanogr.; 2002; 32, pp. 1779-1793. [DOI: https://dx.doi.org/10.1175/1520-0485(2002)032<1779:NEOTBO>2.0.CO;2]

16. Pelinovksy, E.; Talipova, T.; Ivanov, V. Estimations of the nonlinear properties of the internal wave field off the Israel coast. Nonlinear Proc. Geoph.; 1995; 2, pp. 80-88. [DOI: https://dx.doi.org/10.5194/npg-2-80-1995]

17. Shi, X.; Fan, Z.; Liu, H. A Numerical Calculation Method for Eigenvalue Problems of Nonlinear Internal Waves. J. Hydrodyn.; 2009; 21, pp. 373-378. [DOI: https://dx.doi.org/10.1016/S1001-6058(08)60159-6]

18. Fliegel, M.; Hunkins, K. Internal Wave Dispersion Calculated Using the Thomson-Haskell Method. J. Phys. Oceanogr.; 1975; 5, pp. 541-548. [DOI: https://dx.doi.org/10.1175/1520-0485(1975)005<0541:IWDCUT>2.0.CO;2]

19. Cai, S.; Gan, Z. A numerical method of internal waves dispersion relation. J. Trop. Oceanogr.; 1995; 1, pp. 22-30. (In Chinese)

20. Keben, Z.; Hushan, G. A numerical method for vertical structure of internal waves in ocean. J. Trop. Oceanogr.; 1997; 4, pp. 62-67. (In Chinese)

21. Bing, D.; Xiang, Z.; Ming, Z. Calculation and Analysis of Vertical Structure of Internal Wave in Background Current. Adv. Mar. Sci.; 2014; 2, pp. 121-129. (In Chinese)

22. Kurkina, O.; Rouvinskaya, E.; Kurkin, A.; Giniyatullin, A.; Pelinovsky, E. Vertical structure of the velocity field induced by mode-I and mode-II solitary waves in a stratified fluid. Eur. Phys. J. E Soft Matter.; 2018; 41, 47. [DOI: https://dx.doi.org/10.1140/epje/i2018-11654-3]

23. Zhao, Z.; Alford, M.H. Source and propagation of internal solitary waves in the northeastern South China Sea. J. Geophys. Res.; 2006; 111, C11012. [DOI: https://dx.doi.org/10.1029/2006JC003644]

24. Grimshaw, R.; Pelinovsky, E.; Talipova, T.; Kurkina, O. Internal solitary waves: Propagation, deformation and disintegration. Nonlinear Proc. Geoph.; 2010; 17, pp. 633-649. [DOI: https://dx.doi.org/10.5194/npg-17-633-2010]

25. Xudong, Z.; Jing, W.; Lina, S.; Junmin, M. Study on the amplitude inversion of internal waves at Wenchang area of the South China Sea. Acta Oceanol. Sin.; 2016; 7, pp. 14-19.

26. Chang, M.H.; Lien, R.C.; Lamb, K.G.; Diamessis, P.J. Long-Term Observations of Shoaling Internal Solitary Waves in the Northern South China Sea. J. Geophys. Res.; 2021; 126, pp. 2169-9275. [DOI: https://dx.doi.org/10.1029/2020JC017129]

27. Song, H.; Gong, Y.; Yang, S.; Guan, Y. Observations of Internal Structure Changes in Shoaling Internal Solitary Waves Based on Seismic Oceanography Method. Front. Mar. Sci.; 2021; 8, pp. 2296-7745. [DOI: https://dx.doi.org/10.3389/fmars.2021.733959]

28. Vlasenko, V.; Brandt, P.; Rubino, A. Structure of Large-Amplitude Internal Solitary Waves. J. Phys. Oceanogr.; 2000; 30, pp. 2172-2185. [DOI: https://dx.doi.org/10.1175/1520-0485(2000)030<2172:SOLAIS>2.0.CO;2]

29. Geng, M.; Song, H.; Guan, Y.; Bai, Y. Analyzing amplitudes of internal solitary waves in the northern South China Sea by use of seismic oceanography data. Deep-Sea. Res. Part I; 2019; 146, pp. 1-10.

30. Wang, S.-D.; Wei, G.; Du, H.; Wu, J.-L.; Wang, X.-L. Experimental investigation of the wave-flow structure of an oblique internal solitary wave and its force exerted on a slender body. Ocean Eng.; 2020; 201, 107057. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2020.107057]

31. Zhi, C.; Wang, H.; Chen, K.; You, Y. Theoretical and experimental investigation on strongly nonlinear internal solitary waves moving over slope-shelf topography. Ocean Eng.; 2021; 223, 108645. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2021.108645]

32. Du, H.; Wang, S.-D.; Wang, X.-L.; Xu, J.-N.; Guo, H.-L.; Wei, G. Experimental investigation of elevation internal solitary wave propagation over a ridge. Phys. Fluids; 2021; 33, pp. 1-9. [DOI: https://dx.doi.org/10.1063/5.0046407]

33. Tian, Z.; Jia, Y.; Du, Q.; Zhang, S.; Guo, X.; Tian, W.; Zhang, M.; Song, L. Shearing stress of shoaling internal solitary waves over the slope. Ocean Eng.; 2021; 241, pp. 0029-8018. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2021.110046]

34. Wang, J.; Zhang, M.; Mei, Y.; Lu, K.; Chen, X. Study on Inversion Amplitude of Internal Solitary Waves Applied to Shallow Sea in the Laboratory. IEEE. Geosci. Remote. Sens. Lett.; 2021; 18, pp. 577-581. [DOI: https://dx.doi.org/10.1109/LGRS.2020.2985123]

35. Zhang, M.; Wang, J.; Li, Z.; Liang, K.; Chen, X. Laboratory Study of the Impact of the Surface Solitary Waves Created by the Internal Solitary Waves on Optical Imaging. J. Geophys. Res.; 2022; 127, pp. 2169-9275. [DOI: https://dx.doi.org/10.1029/2021JC017800]

36. Vlasenko, V.; Hutter, K. Generation of second mode solitary waves by the interaction of a sill. Nonlinear Proc. Geoph.; 2001; 8, pp. 223-239. [DOI: https://dx.doi.org/10.5194/npg-8-223-2001]

37. Brandt, A.; Shipley, K.R. Laboratory experiments on mass transport by large amplitude mode-2 internal solitary waves. Phys. Fluids; 2014; 26, 046601. [DOI: https://dx.doi.org/10.1063/1.4869101]

38. Deepwell, D.; Stastna, M.; Carr, M.; Davies, P.A. Wave generation through the interaction of a mode-2 internal solitary wave and a broad, isolated ridge. Phys. Rev. Fluids; 2019; 4, 094802. [DOI: https://dx.doi.org/10.1103/PhysRevFluids.4.094802]

39. Walker, L.R. Interfacial solitary waves in a two-fluid medium. Phys. Fluids; 1974; 16, 1796. [DOI: https://dx.doi.org/10.1063/1.1694215]

40. Koop, C.G.; Butler, G. An investigation of internal solitary waves in a two-fluid system. J. Fluid Mech.; 1981; 112, pp. 225-251. [DOI: https://dx.doi.org/10.1017/S0022112081000372]

41. Kao, T.W.; Pan, F.-S.; Renouard, D. Internal solitons on the pycnocline: Generation, propagation, and shoaling and breaking over a slope. J. Fluid Mech.; 1985; 159, pp. 19-53. [DOI: https://dx.doi.org/10.1017/S0022112085003081]

42. Ostrovsky, L.A.; Stepanyants, Y.A. Internal solitons in laboratory experiments: Comparison with theoretical models. Chaos; 2005; 15, 37111. [DOI: https://dx.doi.org/10.1063/1.2107087]

43. Chen, M.; Chen, K.; You, Y.-X. Experimental investigation of internal solitary wave forces on a semi-submersible. Ocean Eng.; 2017; 141, pp. 205-214. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2017.06.027]

44. Barros, R.; Choi, W.; Milewski, P.A. Strongly nonlinear effects on internal solitary waves in three-layer flows. J. Fluid Mech.; 2020; 883, A16. [DOI: https://dx.doi.org/10.1017/jfm.2019.795]

45. Michallet, H.; Barthelemy, E. Experimental study of interfacial solitary waves. J. Fluid Mech.; 1998; 366, pp. 159-177. [DOI: https://dx.doi.org/10.1017/S002211209800127X]

46. Grue, J.; Jensen, A.; Rusas, P.-O.; Sveen, J.K. Properties of large-amplitude internal waves. J. Fluid Mech.; 1999; 380, pp. 257-278. [DOI: https://dx.doi.org/10.1017/S0022112098003528]

47. Du, H.; Wei, G.; Wang, S.-D.; Wang, X.-L. Experimental study of elevation- and depression-type internal solitary waves generated by gravity collapse. Phys. Fluids; 2019; 31, 102104.