Abstract

This paper is concerned with the sinking of the Korean naval warship (ROKS Cheonan) and the reported spectra of the seismic signals recorded at the time of the incident. The spectra of seismic signals show prominently amplitude peaks at around 8.5 Hz and its harmonics. These frequencies were explained with the vibrations of a water column due to an underwater explosion. This explanation is highly doubtful and concerns about its validity have already been raised in the scientific community. In this work an alternative explanation is presented: it is shown that the recorded seismic spectra are consistent with the natural frequencies of vibrations of a large submarine with a length of around 113 m. This finding raises the possibility that the ROKS Cheonan sunk because of the collision with a large submarine rather than the explosion of a torpedo or an underwater mine.

1. Introduction

Several years have passed since the incident of the Korean naval warship (ROKS Cheonan) sinking that occurred on the 26th of March 2010 [1, 2]. The Cheonan warship was split largely into two parts as shown in Figure 1 and sank off Korea’s west coast near Baekryeong Island in the Yellow Sea. The seismic signals were recorded at the time of the incident in many stations including Baekryeong Island. The analysis of the seismic signals followed soon after the incident as reported in [3]. Based on the official government report of a summary of its investigation on 20 May 2010, the warship was sunk by a North Korean torpedo fired by a midget submarine (JIG’s report [1] printed in September 2010). Nevertheless controversy over the scientific legitimacy of the report has been aroused among scientists [48]. The doubts raised by each of the reference are summarized in Appendix A.

This work was motivated from the report of [3, 9] showing the spectra of the seismic signals generated from the Cheonan incident: the spectral Figure 7 in [3] shows very interesting spectral peaks at around 8.5 Hz and its harmonics. The author in [3] has explained these characteristic frequencies as those of the hydroacoustic reverberation waves generated by an underwater explosion. This explanation is based on the theory in the report of [10] that underwater explosion generates reverberation waves in the column of seawater with characteristic odd harmonic series of a fundamental frequency by , where is the wave speed in water (~1.5 km/s) and is the water column depth (for  Hz,  m). However contrary to this prediction, Figure 7 in [3] clearly shows a series including  Hz and  Hz, which is not an odd harmonic series. Thus it is very questionable for the validity of the analysis result in [3]. Recently in another paper based on the same underwater explosion theory as in [9] where the explosive source was indicated as a South Korean land control mine rather than a North Korean torpedo this problem was recognized and the authors stated that “these unusual unexplained spectral peaks at 17 Hz and its second harmonic frequency at 34 Hz can be attributed to the nonlinear deformation of the hull structure by shock waves along with overlapping readings from later arriving T-phases.”

This statement in our opinion is a speculation and no proof or simulations to support it were shown in [9]; besides the corresponding amplitude matching of the recorded spectra had not been attempted in [3, 9].

We think that it is worthwhile to search an alternative theory to explain more properly the spectra including amplitudes of the seismic signals. For instance, it will be shown in this work that a possible explanation for the sinking of the ROKS Cheonan warship is a collision with a large submarine. The reason for this explanation relies on the fact the frequencies observed in the spectra of the seismic data are consistent with the natural frequencies of vibrations of a large submarine, as it will be shown later in the paper. The collision of the ship with the submarine would generate a large force transmitted to the submarine hull with the consequent structure born sound radiated in the sea-water. The sound waves then transmit into the earth crust as seismic waves and can be recorded at seismic detecting stations. The purpose of this paper is to explore this possibility to gain insight into the reality of the cause of the Cheonan sinking. For this aim both a simplified and more accurate model of a submarine hull will be presented in the following sections. A simplified model will be used first to find the main dimension of the submarine hull; a more realistic model will be then used to properly match the seismic signature.

2. Theoretical Formulation and Calculation

2.1. Natural Vibrations of a Submarine Hull

A model of a submarine is presented according to the approach of [11, 12], where the submarine was modeled as having a main cylindrical hull with internal bulkheads and ring stiffeners. The cylindrical shell is closed by truncated conical shells, which are closed by circular plates at each end as shown in Figure 2(a). Hereafter we call this the realistic model of a submarine; however we also consider a simplified version, as shown in Figure 2(b), where the hull is modeled as a cylindrical shell with ring stiffeners closed by circular plates at ends of the shell. This is because the effect of the end cones on the structural response is minimal [11] and the purpose of using the simplified model is to get a closed form solution for the main parameters of the hull and get a first matching with the recorded seismic frequencies. The detailed model as in [12] will be used later to calculate a more realistic structural and acoustic response of the submarine.

At first, let us consider a cylindrical shell of length with shear diaphragm boundary conditions to model the main part of the hull (Figure 2(b)). The equations of motion and the solutions and the boundary condition equation are given in Appendix B. The solution of the free vibration of the cylindrical shell satisfying the Flügge differential equations of motion can be written as given by Caresta et al. [12, 13]: Here , , and are the orthogonal components of shell displacements in the -axial, the -circumferential, and the -radial directions, respectively. is the circumferential mode number (, 1, 2), is the imaginary unit, is the angular frequency, and is the axial wave number. The shear diaphragm boundary conditions are given by where is the membrane force and the is the bending moment. With these boundary conditions, the axial wave number is given by Substituting (1) and (3) into the Flügge equations gives the eigenvalue equation in the matrix form In (4)    is the dimensionless frequency parameter given by . The elements of the matrix are given by where In the previous equations, is the mean radius of the cylindrical shell, is the shell thickness, is the cross-section area of the ring stiffener, is the stiffener spacing, and is the distance between the shell mid surface and the center of a ring. , , and are Yong’s modulus, density, and Poisson’s ratio of the shell, respectively. is the longitudinal wave speed in the shell and is the equivalent distributed mass on the shell to take into account the onboard equipment and ballast tanks. For the each set of (, ) the matrix equation of (4) gives three real eigenvalues of and the corresponding eigenvectors (, , and ) under the usual normalization condition .

2.2. Axisymmetric Modes of Vibrations

If , the wave modes are independent of and if we do not consider the torsional modes, we can set . The matrix equation (4) simplifies in a matrix equation: where Under the condition, , the off-diagonal elements become small values and then the solutions of (6a) approximately give for the nonzero of The first equation (7) with (3) and the first relation in (5j) and in (6b) give It is very important to note that (8) is a harmonic series with a fundamental frequency of , which can be used to match the characteristic frequencies of the seismic signals reported in [3, 9]. The validity of (8) in the approximation can be seen up to the wave number from the solution of (6a). The first 4 frequencies and the corresponding normalized eigenvectors of (6a) are listed in Table 1. Using (8) and (5j), we can obtain analytically the length of a submarine giving  Hz as to match the recorded seismic signature. The submarine hull is assumed to be made of steel with physical properties of 800 kg/,  N/, and ;  km/s. It is widely known that the typical dimension of submarines has the ratio of ~0.1 and   ~0.01. In , and can be used as in [12]. The stiffeners parameters and give minor contribution in so it is not necessary to know accurate values. is supposed to be less than the value estimated by the relation where and are the volume and the surface area of submarine, respectively. is the density of seawater, about 1025 kg/m3. If and , . Inserting these values in (8) gives . This value is the typical length of a large submarine. In summary we have , (radius of the cylindrical model), (thickness of the hull), and  kg/m2 giving . This model of submarine according to (8) has the natural frequencies as a harmonic series of 8.5 Hz. In this result the equivalent mass is found to be  kg/m2 and it seems a little too high. For a submarine with a length of 113 m, considering about 400 tons added at the ends to simulate the ballast tanks a more realistic value should be around ~2150 kg/m2. If this distributed mass is used, and becomes 9.6 Hz. However by considering the realistic model of the submarine with the conical end caps and the fluid loading with   ~2150 kg/m2, the match of the frequencies is still achieved, as can be seen in Table 1. Therefore we will keep using the simplified model of a submarine in vacuum with  kg/m2 for a quick guide. It has been noted that the adjusting dimensional parameters of length , radius , hull thickness , distributed mass , and so forth within 10% can still give harmonics of 8.5 Hz for the natural frequencies vibrations of the models. This result means the range of the uncertainty of the dimension of this model structure is less than 10%. The solution of the second equation in (7) for the second branch of the eigenvalues gives under the condition of small and small value of Here and we get  Hz. The accurate calculation of the corresponding frequencies gives 81.6 Hz (), 81.7 Hz (), 81.8 (), and 85.0 Hz. Since these frequencies are greater than the cut-off frequency  Hz of the seismic signals reported for this event, we will not consider this branch of frequencies any more.

2.3. Asymmetric Modes of Vibration

For the bending modes the solution of in (4) gives three branches of eigenvalues and the corresponding three eigenvectors. The frequency for each can be calculated from the eigenvalues of . The lowest frequency branch represents mainly radial motion and is of interest here; the first 6 frequencies and the corresponding normalized eigenvectors are listed in Table 1. The corresponding resonant frequencies calculated with the realistic model are also listed in Table 1. The frequencies of the second and the third branches of are greater than the frequency range investigated and they will not be considered in this study. The vibrations with are out of consideration in this work because those give little contribution to the response as it will be shown later. It should be noted that in Table 1 the difference between the frequencies calculated using the simplified model with  kg/ and the realistic model with  kg/m2 is small supporting the usefulness of the simplified model.

3. Structural and Acoustic Responses

3.1. Structural Response

Figures 3 and 4 show the structural frequency response function (FRF) of the complete submarine model for the axial and radial displacement at the driving point, as applying a point unit force at one end of the cylindrical hull in axial direction. For the reciprocity principle the FRFs shown can be also interpreted as the response for the radial and axial displacement in response to a radial force applied at the same point. The calculations in these plots were carried out using the modeling approach presented in [12]. The contribution of the and modes is also shown in the same figure; modes with were excluded by the simulation since it did not give a significant contribution to the response. It has to be noted that those results are for a point force, in case the excitation force is well distributed around the hull, such as in the case of a front collision, the contribution on the response would be given only by the modes. The natural frequencies corresponding to the sharp peaks in the FRF are summarized in Table 1.

The striking feature in Figure 3 is the resonant frequencies of the mode, being around 8.6 Hz, 17.0 Hz, 25.0 Hz, and 36.0 Hz. The contribution of the modes is also important and an increase in response is observed at around 35.0 Hz where the waves of Class II (the second branch of W2) cut on, as it was also observed in [12]. These results show a qualitative similarity between the seismic signature and the natural frequencies of a large submarine, validating the possibility of a collision discussed before. The result also suggest a collision with the front part of the submarines with a main excitation of the axisymmetric () modes of vibration. It can be seen in Figure 4 that the contribution of to the radial displacement is much lower than the axial displacements as shown in Figure 3. This aspect is consisting with the fact of low ratios of in the eigenvectors , as in Table 1. Also it can be noted that the modes have a mainly radial displacement. This aspect can be expected by noting the ratios of in the eigenvectors , , and in Table 1.

It has to be mentioned that the purpose of this calculation is to show the FRF of the submarine model and its natural frequencies, and a possible qualitative correlation with the seismic data recorded. In reality the exact impact point on the hull, the actual force transmitted to the submarine, and its frequency content would surely affect the quantitative results but an exact reconstruction of the impact it is surely a huge and uncertain task and is out of the scopes of this work.

3.2. Radiated Sound at Far Field

The sound pressure radiated at far field was also calculated using the model in [12] and the corresponding wave displacement is shown in Figure 5 normalized as a displacement per unit force and per km. The far field location is chosen to be in front of the submarine as a standard point for discussion. The trend of the figure is not much different from the seismic spectra reported in [3, 9] as discussed in Section 4. However this matching does not mean that the front of the submarine was toward to the BAR station, where Korea Meteorological Administration Station (37.9771N, 124.7142E) in Baekryong Island is [9] from where the seismic data analyzed in [3, 9] was obtained. The reason is that the waves from the vibration source of the conical shells of the submarine will propagate spherically in the far field through mediums. The path of this spherically propagated waves arrived at the BAR station was suggested in Figure 6 of [3] (see Figure 6 in this paper).

4. Discussion

Figure 7 of [3] showed the frequency spectrum of the -wave signals in -direction which was the Fourier transform of the first arrived -waves at the BAR recording system in one second time window. Also Figure 3(b) of [9] showed the normalized spectra of seismic signals in a time window of 15 seconds for each , EW (east-west direction) and NS (north-south direction), separately using the same BAR recordings of the seismic signals. In this case the analyzed signals contain not only the first arrived -wave, but the later arrived -wave, Rayleigh, Love waves, and so forth. The Fourier transforms of all these waves gives the spectrum which is the sum of the spectrum of each signal travelled via several independent paths. As expected the spectrum in each direction of , EW, and NS in Figure 3 of [9] showed similar patterns to each other. For the comparison of the spectrum in Figure 5 with these reported spectra, the data of the spectra peaks in Figure 7 of [3] and in Figure 3(b) of [9] are summarized in Tables 2 and 3, respectively. The frequencies of the peaks of Figure 5 are compared with the peaks of the reported spectra of Figure 7 of [3] in Table 2 and with the peaks of Figure 3(b) of [9] in Table 3.

In [9] the original unit scale was not specified in the spectra, but supposed to be in Pascal. Thus the difference between those in D column in Table 2 and those in the average in Table 3 may indicate uncertainty of amplitudes in the recorded seismic spectra. On the other hand, the peak ratios in Figure 5 are estimated to be 1.0 for 8.6 Hz, 0.23 for 17.0 Hz, 0.04 for 25.0 Hz, and 0.016 for 36.0 Hz as listed in Cal.-column in Table 2. These ratios are close to those in Table 2: C-column, implying good fitting for both frequencies and amplitudes ratios too with the reported spectra. We also noted that the frequencies of natural vibrations of the hull frame of ROKS Cheonan were calculated and listed in Table -6-2 (p160) in [1], as 2.32 Hz, 4.74 Hz, 7.71 Hz, 10.41 Hz, and 13.40 Hz in the bending motion only. (No axial vibration is expected for the hull frame of the warship.) Clearly these natural frequencies could not be correlated with the characteristic harmonic frequencies (8.5 Hz and the multiples) of the seismic signals recorded.

In summary the main peaks at 8.5 Hz and its multiples in the reported spectra of the seismic signals are reasonably consistent with the natural frequencies of the submarine with a length of around 113 m. This result suggests a possibility of the collision between the Cheonan warship and the submarine. Then one might wonder about the damage that the submarine would have following a collision. We believe that the submarine would be negligibly damaged by the collision since the thickness of the hull of a large submarine is supposed to be more than 6 cm with and made of high strength steels, while the hull thickness of the ROKS Cheonan is known to be about 1.2 cm with steels and aluminum alloys in the upper parts. Besides we have found that it is no difficult with this collision theory discussed in this paper to illustrate damage aspects including deformations on the recovered bow and the stern of the split Cheonan warship from the sinking.

5. Conclusion

This study shows that the characteristic frequencies at spectral peaks of the seismic signals recorded at the time of the Cheonan incident are consistent with the natural frequencies of vibrations of a large submarine with a length of around 113 m. The matching is particularly good with the axisymmetric () modes of vibrations suggesting a front collision with the submarine that might have caused the ROKS Cheonan sinking. The authors hope this work may be the starting point for new investigations to shine some light on the mysterious causes of the Cheonan sinking that caused the death of 46 people and for which a clear and definitive explanation has not been given yet.

Appendices

A. Doubts about the JIG’s Report

Several doubts have been raised about the JIG’s report and they are explained in what follows. The first issue of the controversy over the JIG’s report [1] is on the validity of composition analysis results of the Al-containing white powder absorbed. These materials were found from the following three locations:(A)the split ROKS Cheonan recovered in April of 2010;(B)the remnants of a torpedo retrieved near the incident site on 15 May 2010;(C)the explosion products from the small-scale explosion experiment using a highly aluminized explosive source of 15 g in a tank filled with 4.5 tons of seawater.

The composition analysis was performed with data obtained from SEM (scanning electron microscopy), EDS (energy dispersive spectrometer), and XRD (X-ray diffraction) for all three cases. The analysis was based on their knowledge or assumption: XRD data of an amorphous aluminum oxide will not show any noticeable diffraction peaks. When the aluminum is exposed to moisture, acids, and bases as in seawater environment for a long time, it forms naturally white corrosion products; the major components of these corrosion products are aluminum hydroxide (Al(OH)3, bayerite) along with boehmite (AlO(OH)) and Al2O3, all of which are known to be crystalline, rather than amorphous. In the JIG analysis the EDS data of A, B, and C samples showed commonly prominent intensity peaks of oxygen and aluminum with ratio in I(O) : I(Al) = 0.9 : 1. JIG regarded this result as evidence that the absorbed materials of A, B, and C samples were aluminum oxides. The XRD data of the A and B samples did not show any noticeable X-ray diffraction peaks of an aluminum oxide crystal and the XRD data of C sample showed prominent aluminum Bragg diffraction peaks and other weak peaks. These weak peaks were analyzed to be irrelevant to aluminum oxide crystals. However JIG did not give explicitly the reason why the strong Al Bragg peaks appeared in the XRD data and they thought that the XRD data of C sample also did not show any noticeable X-ray diffraction peaks relevant to the aluminum oxide. These results led JIG to the conclusion that the A, B, and C samples were almost all (~100%) an amorphous aluminum oxide by . Because the product was an amorphous material, these materials could not be natural corrosion products of aluminum which are known to be crystalline by . Then JIG concluded that all these results were clear evidence that the adsorbed materials of A and B were amorphous aluminum oxides, which are the explosion product of the torpedo whose remnants were found near the incident site. We think the steps stated above were JIG’s chain logic to draw the conclusion.

The second issue about JIG’s report is whether the torpedo remnants were genuine or not. Inside the rear section of the torpedo remnants, JIG found Korean handwrite marking “1bun (No. 1 in English in blue ink)”, similar to the marking of a North Korean test torpedo obtained in 2003, as seen in Figure Summary-4 and Figure Summary-5 of JIG’s report. JIG declared promptly “we found conclusive evidence” for the cause of the incident. JIG said these evidences confirmed that the recovered torpedo parts were manufactured in North Korea and concluded that ROKS Cheonan was split and sunk due to shockwave and bubble effects generated by the underwater explosion of a torpedo, manufactured by North Korea. However in JIG’s composition analysis data, Lee [4] indentified the XRD weak peaks for C sample, in fact, came from diffraction of a small amount of well crystallized -Al2O3, contrary to JIG’s claim that the weak peaks were irrelevant to aluminum oxide crystals. One can see clearly weak -Al2O3 Bragg diffraction peaks in Figure Appendix v-5-3, p280 in JIG’s report [1]. We think this finding broke JIG’s chain logic.

This result means that at least the C sample was quite different from the A and B samples, and in turn JIG’s conclusion above is very questionable, as Lee has concluded. Lee and Yang then reported their own experimental results in another study [5]. They presented the XRD data in Figure 2 of [5] for the sample of an Al-powder that underwent melting followed by rapid quenching. This figure showed prominent Al Bragg diffraction peaks and weak -Al2O3 Bragg diffraction peaks. They said that “this clearly indicates Al-powder oxides partially, not entirely, during the heating and rapid quenching.” This oxidation was supposed to occur on the surface thin layers of each grain of Al powder. The XRD data in Figure 2 of [5] showed the similar pattern as in Figure Appendix V-5-3 of JIG’s report for C sample. Nevertheless the EDS data in Figure 1 of [5] shows markedly different from Figure Appendix V-5-2 of JIG’s report for C sample. That is, the peak ratio of I(O)/I(Al) in Figure 1 of [5] was measured as 0.25, whereas Figure Appendix v-5-2 showed very different values of 0.81. On the other hand, the authors in [5] noticed that the high value of around 0.8-0.9 for I(O)/I(Al) in EDS data was expected for aluminum hydroxide, such as Al(OH)3. They also showed the simulation results of Al(OH)3 and Al2O3 in Figure 4 of [5], supporting their expectation. (It should be noted that the EDS intensity data comes from the surface thin layers of the sample. That is, the EDS data of Figure 1 of [5] came from the surface thin layers of oxidized grains with -Al2O3 of the heat treated Al-powder sample.) From these results they concluded that JIG’s adsorbed materials taken from the warship and the torpedo remnants were not associated with any explosion and that JIG’s EDS data of their test-explosion sample (Figure Appendix v-5-2) were likely fabricated. The South Korean government has adamantly denied any fabrication or any major problems with its interpretation of the data. Later on, independently the A and B samples were analyzed to be a basaluminite [Al4(SO4)(OH)10·4-5H2O] materials by Professor Jeong GY in Earth & Environmental Science from Andong National University in Korea, using several instruments including a TEM. He said that the materials appeared to be formed for a long time such as in seawater, implying they had nothing to do with explosions. Unfortunately he was not able to test the C sample.

Cyranoski [6] reported in Nature online the summary of the controversy over JIG’s report. Some important points raised in the article are as follows.(a)An expert investigator was placed on the JIG by the opposition party—Shin Sang-chul, a former officer in the South Korean navy who had also worked at a shipbuilding company—suggested, before the report was even released, that an accidental collision with a US warship, and not North Korea, was to blame. The United States and South Korea had been carrying out military exercises in the area at the time. Now he is insisting on the collision with a submarine with a length of around 60 m, not by a US warship.(b)The report’s claim that a torpedo-induced water column sank the Cheonan contradicted earlier testimony from survivors that they did not see that water column. They also testified neither seeing any explosion flash from the underwater nor smelling of explosives at the time of the incident.

Lee and Suh in [7] reported in Policy Forum 10-039 “Inconsistencies in South Korea’s Cheonan Report.” In the article we found particularly interesting the following point.

“If the bottom of the ship was hit by a bubble, it should show a spherical concave deformation resembling the shape of a bubble, as the JIG’s own simulation suggests in the Appendix of the report, but it does not. The bottom of the front part of the ship is pushed up in an angular shape, as the yellow line shows in Figure III-1-7 of the report, more consistent with a collision with a hard object. Suh has further argued that, contrary to common reports that the Cheonan was split in half as also seen in JIG’s simulation, it actually broke into two larger pieces as well as a third smaller piece (see Figure 1 in this paper).” Shin Sang-chul clearly pointed out if the warship was split by a bubble jet generated from the non-contact underwater explosion, the top hull of the recovered broken warship should show upward bending, but it does not seem the case as can be seen in Figure II-3-3 and 4 and Figure III-1-8 of the report.

Kim et al. in “Foreign Policy In Focus” [8] reported briefly skeptics on JIG’s judgment about the Korean handwrite marking “1bun” on the rear section of the torpedo remnants. JIG claimed this handwriting marking had been written by a North Korean in the process of manufacturing the torpedo and is the conclusive evidence that the torpedo was made in North Korea. Nevertheless Shin Sang-chul found the trace that the rusts on the panel of the rear section of the torpedo remnants had been removed by a sand-paper; after the rust was removed the marking of Korean words “1bun” was written on the panel. He said, therefore, this marking could not be “conclusive evidence” that the torpedo was made in North Korea. One can judge whether his judgment is correct or not from Figure Summary-4 of the report. We personally think Shin Sang-chul is right. It is reminded that the remnants were fragments following the torpedo explosion and had been deposited in the mud of seabed for 50 days and recovered by the special fishing net [1]. Considering all the several doubts explained above we believe the torpedo remnants they found near the incident site could not be a genuine conclusive evidence for the cause of the incident. Furthermore the conclusion was drawn after only five days of work (the remnant was found on 15 May 2010 and JIG’s report of a summary of its investigation was released on 20 May 2010).

B. The Free Vibrations of a Cylindrical Shell with Shear-Diaphragm Boundary Conditions

The Flügge equations of motion for the vibrations of a ring stiffened cylindrical shell are given by Caresta and Kessissoglou [12]: Here , , and are the orthogonal components of shell displacements in the -axial, the -circumferential, and the -radial directions in time , respectively. The meanings of the symbols in the equations are defined in Section 2. The general solutions for shear diaphragm boundary conditions can be written as in [13]: These solutions represent standing waves in both the axial and circumferential directions, with nodal lines in the direction and with nodal cross-sections spaced at a distance in the direction. is the axial wave number and is the circumferential mode number. is the imaginary unit. The shear diaphragm boundary conditions at , , are given by The membrane force and the bending moment are expressed as functions of the normal strains in the middle surface, and , as well as the mid-surface changes in curvature, , . The expressions for the strain and changes in curvature are given by Leissa [13] as As explained in [13], the shear diaphragm boundary conditions in (B.5) to (B.8) can be justified for a rigid thin, circular cover plate at each end, so that and displacements are restrained at both end boundaries. However the plates, by virtue of their thinness, would have very little stiffness in the direction transverse to their planes; consequently, they would generate negligible bending moment and longitudinal membrane force in the shell as the shell deforms. Under these boundary conditions, the axial wave number becomes , where is the axial mode number. The natural frequencies of the shell can be found substituting the solution (B.4) with into (B.3) and arrange in matrix form. For a nontrivial solution the determinant of the matrix must be zero. Expansion of the determinant leads to the characteristic equation of the shell. In vacuum, the dispersion equation is of sixth order in . From this equation three different natural frequencies can found for each value of .

B.1. Frequency Response Function (FRF)

Suppose that the shell is axially excited at one end by a point force of unity amplitude located at (, ). The force can be described in terms of a Dirac delta function in function of the arc length ; The equilibrium of the membrane force, , given by (B.5) evaluated at becomes The delta function can be expanded as a Fourier series around the circumference as With , (B.11) can be rewritten as For every circumferential mode number , the 8 boundary equations, as expressed in [12] for the 3 membrane force, 3 bending moments, transverse shearing, and the Kelvin-Kirchoff shear force-together with the equilibrium of the forces under point force excitation, can be arranged in the matrix form . is the vector of the 8 unknown displacement coefficients and is an force vector with only one nonzero term , where if and if . Structural damping can be introduced using a complex Young modulus where (about 0.02) is the structural loss factor. Solving the system for each circumferential mode number gives the steady state shell displacement response at a certain frequency (FRF).

C. The Path of the Signals Arriving at the BAR Recording System

The optimized path of the signals arriving at the BAR station from the submarine is drawn in Figure 6. In this figure the waves having the velocity 1.5 km/s in seawater with the incident angle 10.9° at o on the boundary of the seabed/crust, refracted into the crust with the angle 52.4° in case of -waves (the velocity 6.3 km/s in crust) by Snell’s law. Under this condition the refracted waves will have maximum amplitudes with the transmission coefficient 0.174. Then the waves reflected on the Moho boundary should continuously travel to the BAR recording system. This path is drawn as o-o′-d in Figure 6. For the case of -waves with velocity of 3.5 km/s in crust, the incident angle on the boundary must be 19.9° to have the same path of o-o′-d. For all and -waves the refraction angle is 90° at o′ with the zero transmission coefficients. To have this optimized o-o′-d with 14.22 km distance, the depth of the Moho boundary must be ~4.3 km. It is generally known that the depth of Moho boundary under sea is usually 5~7 km. As considering this fact, we think the estimated depth of the Moho boundary about 4.3 km is acceptable in this case.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.