Tertiary Waves Measured during 2017 Pohang Earthquake Using an Underwater Glider

: An underwater glider equipped with a hydrophone observed the acoustic sounds of an earthquake that occurred on 15 November 2017 05:29:32 (UTC) in the Pohang area. The underwater glider observed the earthquake sounds after 19 s (05:29:51) at approximately 140 km from the Pohang epicenter. In order to distinguish the earthquake sound from the glider’s operation noise, the noise sources and Sound Pressure Level (SPL) of the underwater glider were analyzed and measured at laboratory tank and sea. The earthquake acoustic signal was distinguished from glider’s self-noises of ﬁn, pumped Conductivity-Temperature-Depth proﬁler (CTD) and altimeter which exist over 100 Hz. The dominant frequencies of the earthquake acoustic signals due to the earthquake were 10 Hz. Frequencies at which the spectra had dropped 60 dB were 50 Hz. By analysis of time correlation with seismic waves detected by ﬁve seismic land stations and the earthquake acoustic signal, it is clearly shown that the seismic waves converted to Tertiary waves and then detected by the underwater glider. The results allow constraining the acoustic sound level of the earthquake and suggest that the glider provides an e ﬀ ective platform for enhancing the earth seismic observation systems and monitoring natural and anthropogenic ocean sounds. noise level before the earthquake occurred. The earthquake acoustic signal can be distinguished from glider’s self-noises of ﬁn, pumped CTD and altimeter. The ﬁn noise exists at frequencies of 800–950 Hz every 3 to 4 s. The pumped CTD sensor noise exists at discrete frequencies of 357, 410, 710, 1450 Hz. The dominant CTD noise frequency is 710 Hz. The altimeter noise exists at frequencies above 100 Hz.


Introduction
Seismic energy from submarine earthquakes is converted into acoustic energy at the seafloor-water boundary. A Tertiary wave (or T-wave) is the acoustic signal from these earthquakes. A T-wave typically has frequencies ranging from 4 to 50 Hz. T-waves propagate efficiently in the ocean compared to seismic waves through the earth and can be detected at great distances. Underwater sound observation techniques can directly observe oceanic earthquakes, landslides, volcanic eruptions, and corresponding underwater explosions in the ocean. Although not many, there have been some reports of direct observations of hydro-acoustic waves due to underwater seismic activity from an offshore observatory [1][2][3][4]. These are the first direct measurements ever performed in a tsunami generation area.
It is becoming increasingly important to acquire low-cost, high-efficiency, small-scale, unmanned, direct observation technology to replace existing high-cost, low-efficiency, large-scale, manned exploration platforms such as marine research vessels and stations in [1][2][3][4][5][6]. Unmanned observation platforms such as Argo floats, Underwater Gliders, WaveGliders, Autonomous Underwater Vehicles and Autonomous Surface Vehicles have been widely used in marine surveys and military defense operations. Among these platforms, underwater gliders have been actively used for measurement of environmental underwater ambient noise. The addition of a single hydrophone sensor to an underwater glider provides a capability for measurement of underwater noise level [7][8][9][10][11][12][13][14] and the addition of a hydrophone array enables surveillance of the underwater acoustic environment [15,16].
Research programs using these technologies have been conducted to investigate geophysical causes of seafloor earthquakes, landslides, volcanic eruptions and corresponding underwater explosions using hydrophone recorders [17,18], an underwater glider [19]. Reference [19] observed underwater sounds at depths of 200 to 600 m using an underwater glider at a distance of 10 km from the Lau Basin submarine volcano. The sound frequency due to the volcano was 10-300 Hz, and the sound level was −25 dB higher than the ambient noise level. The prototype and concept as ocean surface gateway (OSG) using WaverGliders for retransmitting data of ocean bottom package to shore were demonstrated in [20].
In this paper we report the analysis of earthquake sounds recorded by an underwater glider equipped with a single hydrophone during the magnitude 5.5 earthquake of 15 November 2017, at Pohang in the Republic of Korea. The underwater glider measured the earthquake sounds at about 140 km from the epicenter, in a water depth of about 200 m during an underwater noise measurement survey from Uljin to Geoseong along the east coast of South KOREA. In order to distinguish the earthquake sound from the glider's operation noise, the noise sources and Sound Pressure Level (SPL) of the underwater glider were analyzed and measured. Next, the earthquake sounds measured by the gilder was analyzed in time and frequency domain. Finally, the time correlation of the seismic waves and earthquake sounds observed by the seismic land station and the glider were described.

Noise Sources of the Underwater Glider
Underwater gliders are propelled by adjusting the buoyancy, thereby using the resulting vertical motion to drive the glider. Underwater gliders have no propeller, so there is little operating noise. Its possible noise sources, such as the air bladder, fin steering, battery movement and volume piston are described in [21]. Chandrayadula et al. [11] showed that high-intensity sounds such as the sloshing of the glider at the sea surface, the bladder pump, the ballast pump, and the battery packs moving around are all broadband in nature. The bulk of the broadband noises were concentrated at the lower frequencies 0-3 kHz. However, these only describe possible noise sources and no specific frequency and noise-level analysis of each of the noise sources, including the driving and sensor component parts, has yet been studied. Recently, An analysis has been reported on the self-noise of only the driving part of the Petrel II underwater glider developed in China [22].
In this paper, the noise sources of the driving and sensor parts of the underwater glider observed the earthquake sounds are described in detail. A Slocum underwater glider (200 m class) manufactured by Teledyne-Webb Research is used in this study ( Figure 1). The glider has an overall length of 1.5 m and a mass of 60 kg. The buoyancy engine is an electrically powered piston drive, located in the nose section of the glider. The drive allows the glider to take in and expel water, thereby changing its overall buoyancy. The mechanism allows a nearly neutrally buoyant trimmed glider to change its displacement in water by ± 250 cc which corresponds to approximately ± 0.5% of the total volume displaced. This change in buoyancy generates a vertical force which is translated through two swept wings into a combined forward and up/downward motion. Due to the location of the piston drive, also called the buoyancy engine, the change in direction of the buoyant force also creates the main pitching moment for the glider. Besides the buoyancy engine the glider possesses two more control actuators, a 9.1 kg battery pack, referred to as the sliding mass, that can be linearly translated along the main axis of the glider and a rudder attached to the vehicle tail fin structure. The sliding mass is used for fine tuning the pitch angle. A pumped Conductivity-Temperature-Depth profiler (CTD) sensor is contained in the payload bay of the glider. When the underwater glider moves vertically and horizontally, the driving and sensor parts generate noise. The noise sources associated with the driving parts include a piston, pitch controller, fin, and air bladder. The sensor noise sources are a pumped CTD and an altimeter.

Driving-Part Noise Sources
To change its direction of vertical motion the underwater glider piston pulls sea water to change the effective hull volume. This produces a combined noise from the dc motor, reducer, ball screw (linear motion) and piston. When the glider's vertical motion changes, the glider moves the battery pack back and forth to adjust the angle of the bow to match the hull glide angle. This produces a combined noise of dc motor, reducer, lead screw and battery pack sliding. The underwater glider moves the tail fin to control the direction of movement every 3 to 4 seconds. This generates a noise from a small stepper-motor drive. An air bladder operating under the tail antenna prevents the tail antenna from sinking into the water when the underwater glider is surfacing for Iridium communication. An air pump drive starts as the glider ascends through the uppermost 10 m water depth; this pump generates noise.

Sensor-Part Noise Sources
For stable CTD measurement, water is pumped through the CTD sensor and this generates noise. The pumped CTD operates continuously. The rate of pumping can be changed according to whether the glider is descending or ascending. The altimeter is normally used only while diving to prevent the glider from hitting the sea floor. This is yet another source of noise.

Glider's Operation Noise Measurement at Laboratory Tank
Sound pressure level (SPL) is usually expressed as twenty times the logarithm to base 10 of the ratio of the root-mean-square sound pressure to the reference sound pressure: SPL = 20 log (/p0) where ̂ and p0 are the root-mean-square sound pressure and reference sound pressure, respectively. SPL has units of dB re 1 µ Pa 2 . SPL values of the driving-and sensor-part noise sources are discussed.
Noises from both the driving and sensor parts of the underwater glider were measured in a laboratory tank ( Figure 2). The size of the tank is 20 m (W) × 50 m (L) × 10 m (D) and water depth in the tank was about 8 m. Since the water tank contained fresh water, the buoyancy of the underwater glider was adjusted in advance to be neutrally buoyant in fresh water. In order to record the noise of the underwater glider, a hydrophone was mounted on the underwater glider body. The reference value was −168 dB re 1 µ Pa 2 . The gain and resolution of the recorder were 30 dB and 16 bits, respectively.

Driving-Part Noise Sources
To change its direction of vertical motion the underwater glider piston pulls sea water to change the effective hull volume. This produces a combined noise from the dc motor, reducer, ball screw (linear motion) and piston. When the glider's vertical motion changes, the glider moves the battery pack back and forth to adjust the angle of the bow to match the hull glide angle. This produces a combined noise of dc motor, reducer, lead screw and battery pack sliding. The underwater glider moves the tail fin to control the direction of movement every 3 to 4 s. This generates a noise from a small stepper-motor drive. An air bladder operating under the tail antenna prevents the tail antenna from sinking into the water when the underwater glider is surfacing for Iridium communication. An air pump drive starts as the glider ascends through the uppermost 10 m water depth; this pump generates noise.

Sensor-Part Noise Sources
For stable CTD measurement, water is pumped through the CTD sensor and this generates noise. The pumped CTD operates continuously. The rate of pumping can be changed according to whether the glider is descending or ascending. The altimeter is normally used only while diving to prevent the glider from hitting the sea floor. This is yet another source of noise.

Glider's Operation Noise Measurement at Laboratory Tank
Sound pressure level (SPL) is usually expressed as twenty times the logarithm to base 10 of the ratio of the root-mean-square sound pressure to the reference sound pressure: SPL = 20 log (p/p 0 ) wherep and p 0 are the root-mean-square sound pressure and reference sound pressure, respectively. SPL has units of dB re 1 µPa 2 . SPL values of the driving-and sensor-part noise sources are discussed.
Noises from both the driving and sensor parts of the underwater glider were measured in a laboratory tank ( Figure 2). The size of the tank is 20 m (W) × 50 m (L) × 10 m (D) and water depth in the tank was about 8 m. Since the water tank contained fresh water, the buoyancy of the underwater glider was adjusted in advance to be neutrally buoyant in fresh water. In order to record the noise of the underwater glider, a hydrophone was mounted on the underwater glider body. The reference value was −168 dB re 1 µPa 2 . The gain and resolution of the recorder were 30 dB and 16 bits, respectively. Except for the tail, which was used for radio frequency communication with the operator, the underwater glider was fully submerged. Noise was generated by inputting an operation command arbitrarily to measure the SPL of each driving and sensor component ( Figure 3). Fin noise measurement, however, was not possible because the tail was not submerged, to allow for radio frequency communication. The noise analysis of fin operation was obtained from the recorded data at sea, as described in the next section. The pumped CTD did not work due to the low electrical conductivity value of the fresh water used in the tank; so, the pumped CTD noise analysis is also described in the next section. The Piston noise level spectrum rises above the background level by more than 5 dB in all frequency bands up to 10 kHz ( Figure 4a). The Pitch noise level likewise rises by more than 5 dB in frequency bands above 1 kHz (Figure 4b). The Air bladder noise level rises more than 10 dB at 155, 310, 463, 522, 552, 618 Hz and to a lesser extent in a broad frequency band ( Figure  4c). The Altimeter noise level rises more than 5 dB at frequencies above 100 Hz (Figure 4d). To record noise from the flow between the underwater glider and the hydrophone, the glider performed a "yo" (descent and ascent) mission. By measuring the noise from repeated yo missions, it was confirmed that there was no flow noise caused by turbulence between the underwater glider body and the hydrophone.  Except for the tail, which was used for radio frequency communication with the operator, the underwater glider was fully submerged. Noise was generated by inputting an operation command arbitrarily to measure the SPL of each driving and sensor component ( Figure 3). Fin noise measurement, however, was not possible because the tail was not submerged, to allow for radio frequency communication. The noise analysis of fin operation was obtained from the recorded data at sea, as described in the next section. The pumped CTD did not work due to the low electrical conductivity value of the fresh water used in the tank; so, the pumped CTD noise analysis is also described in the next section. The Piston noise level spectrum rises above the background level by more than 5 dB in all frequency bands up to 10 kHz (Figure 4a). The Pitch noise level likewise rises by more than 5 dB in frequency bands above 1 kHz (Figure 4b). The Air bladder noise level rises more than 10 dB at 155, 310, 463, 522, 552, 618 Hz and to a lesser extent in a broad frequency band ( Figure 4c). The Altimeter noise level rises more than 5 dB at frequencies above 100 Hz (Figure 4d). To record noise from the flow between the underwater glider and the hydrophone, the glider performed a "yo" (descent and ascent) mission. By measuring the noise from repeated yo missions, it was confirmed that there was no flow noise caused by turbulence between the underwater glider body and the hydrophone. Except for the tail, which was used for radio frequency communication with the operator, the underwater glider was fully submerged. Noise was generated by inputting an operation command arbitrarily to measure the SPL of each driving and sensor component ( Figure 3). Fin noise measurement, however, was not possible because the tail was not submerged, to allow for radio frequency communication. The noise analysis of fin operation was obtained from the recorded data at sea, as described in the next section. The pumped CTD did not work due to the low electrical conductivity value of the fresh water used in the tank; so, the pumped CTD noise analysis is also described in the next section. The Piston noise level spectrum rises above the background level by more than 5 dB in all frequency bands up to 10 kHz (Figure 4a). The Pitch noise level likewise rises by more than 5 dB in frequency bands above 1 kHz (Figure 4b). The Air bladder noise level rises more than 10 dB at 155, 310, 463, 522, 552, 618 Hz and to a lesser extent in a broad frequency band ( Figure  4c). The Altimeter noise level rises more than 5 dB at frequencies above 100 Hz (Figure 4d). To record noise from the flow between the underwater glider and the hydrophone, the glider performed a "yo" (descent and ascent) mission. By measuring the noise from repeated yo missions, it was confirmed that there was no flow noise caused by turbulence between the underwater glider body and the hydrophone.

Glider's Operation Noise Measurement at Sea
The underwater glider dived to a depth of 10-200 m and repeated six yo missions during about 6 hours. It took about one hour for each yo. Figure 5 shows the command values of the driving part according to the water depth when the underwater glider is operated at sea. The driving-unit command value is −1 (low) to +1 (high).
When the underwater glider changes from descent to ascent, the values of the buoyancy control (Vol Piston) and the Pitch control (Pitch) change from −1 to +1; these values are negative during descent and positive during ascent. The Fin controller operated once every 3 to 4 seconds and the command value changed fast, regardless of descent or ascent. It can be seen that the Air Bladder command value is 0 except near the end when it is 1 as the bladder inflates raising the glider to the surface thereby allowing communication with the operator via satellite.
The Fin and pumped CTD noise levels were calculated from recorded data at sea. The same hydrophone was mounted on the same position of the underwater glider body as at tank experiment. Due to the ambient noise contained in the recordings at sea, the reference sound level of Fin and pumped CTD is greater than that in the laboratory tank. The Fin noise level rose 4-5 dB at frequencies of 800-950 Hz (Figure 6a). The pumped CTD noise level rose by more than 10 dB at discrete frequencies of 357, 410, 710, and 1450 Hz (Figure 6b).

Glider's Operation Noise Measurement at Sea
The underwater glider dived to a depth of 10-200 m and repeated six yo missions during about 6 h. It took about one hour for each yo. Figure 5 shows the command values of the driving part according to the water depth when the underwater glider is operated at sea. The driving-unit command value is −1 (low) to +1 (high).
Appl. Sci. 2019, 9, x FOR PEER REVIEW 6 of 13 Figure 5. Command values when the glider is operating at sea. Figure 5. Command values when the glider is operating at sea.
When the underwater glider changes from descent to ascent, the values of the buoyancy control (Vol Piston) and the Pitch control (Pitch) change from −1 to +1; these values are negative during descent and positive during ascent. The Fin controller operated once every 3 to 4 s and the command value changed fast, regardless of descent or ascent. It can be seen that the Air Bladder command value is 0 except near the end when it is 1 as the bladder inflates raising the glider to the surface thereby allowing communication with the operator via satellite.
The Fin and pumped CTD noise levels were calculated from recorded data at sea. The same hydrophone was mounted on the same position of the underwater glider body as at tank experiment. Due to the ambient noise contained in the recordings at sea, the reference sound level of Fin and pumped CTD is greater than that in the laboratory tank. The Fin noise level rose 4-5 dB at frequencies of 800-950 Hz (Figure 6a). The pumped CTD noise level rose by more than 10 dB at discrete frequencies of 357, 410, 710, and 1450 Hz (Figure 6b). A single omni-directional hydrophone (HTI92MIN from High Tech with a sensitivity of −168 dB re 1V/1μPa) was mounted on the underwater glider. We synchronized the time of the hydrophone with the underwater glider before launched the glider. The signal was digitized by the internal acoustic data logger at a 24 kHz sample rate with 16-bit resolution and 30 dB gain. Each file was approximately 255 MB in size, 60 minutes long; a total of 18.3 GB was stored on the 128 GB internal

In Situ Underwater Ambient Noise Measurement
The measurement survey of underwater ambient noise was conducted using an underwater glider with a single hydrophone from 13 to 17 November 2017 from Ulgin to the Goseong area, along the east coast of South Korea in Figure 7. The underwater glider was launched at Ulgin 13 November 05:04 (UTC) and recovered at Goesong 17 November 06:20 (UTC) using the R/V Eodo of Korea Institute Ocean Science and Technology. The underwater glider moved vertically between depths of about 10 m and 200 m. During the survey, the glider advanced at an average horizontal speed of 0.02 m/s, logging data while repeatedly taking profiles between 10 and 200 m depth. The glider drifted when it surfaced between dives to transmit CTD data and Global Positioning System (GPS) data to the shore station through a satellite connection. The average surface current calculated from the glider's GPS data was northward with a velocity of 0.35 m/s. A single omni-directional hydrophone (HTI92MIN from High Tech with a sensitivity of −168 dB re 1V/1µPa) was mounted on the underwater glider. We synchronized the time of the hydrophone with the underwater glider before launched the glider. The signal was digitized by the internal acoustic data logger at a 24 kHz sample rate with 16-bit resolution and 30 dB gain. Each file was approximately 255 MB in size, 60 min long; a total of 18.3 GB was stored on the 128 GB internal flash memory during the 4-day survey. No acoustic data were transmitted in real time due to the large size of the files. All recorded data were processed using MATLAB software after the recovery of the glider. Recordings from the survey produced a one-dimensional time series of sound level spanning approximately 150 km from Ulgin to Geseong. During the survey, on 15 November 05:29:32 (UTC), a magnitude 5.5 earthquake occurred at Pohang (epicenter 36.12 • N, 129.36 • E) [23,24]. Underwater sounds resulting from the momentary earthquake were observed in the glider's acoustic record.
flash memory during the 4-day survey. No acoustic data were transmitted in real time due to the large size of the files. All recorded data were processed using MATLAB software after the recovery of the glider. Recordings from the survey produced a one-dimensional time series of sound level spanning approximately 150 km from Ulgin to Geseong. During the survey, on 15 November 05:29:32 (UTC), a magnitude 5.5 earthquake occurred at Pohang (epicenter 36.12°N, 129.36°E) [23,24]. Underwater sounds resulting from the momentary earthquake were observed in the glider's acoustic record.

Time and Frequency Analysis of Earthquake Underwater Sound
It is most common for noise signals to be represented in the frequency domain as noise spectral levels plotted as a function of frequency [25,26]. The data for ambient noise is usually expressed as spectral density (SD) levels or power spectral density (PSD), where the data in each frequency band has been normalized by dividing by the bandwidth of the frequency band. The units of the levels in each band are then dB re 1 μPa 2 /Hz. We used the following analysis method to distinguish earthquake underwater sounds. Extraction of necessary timing data from the recording files; Calculation of Power Spectrum Density (PSD) using a time window of chosen length in each

Time and Frequency Analysis of Earthquake Underwater Sound
It is most common for noise signals to be represented in the frequency domain as noise spectral levels plotted as a function of frequency [25,26]. The data for ambient noise is usually expressed as spectral density (SD) levels or power spectral density (PSD), where the data in each frequency band has been normalized by dividing by the bandwidth of the frequency band. The units of the levels in each band are then dB re 1 µPa 2 /Hz. We used the following analysis method to distinguish earthquake underwater sounds. Extraction of necessary timing data from the recording files; Calculation of Power Spectrum Density (PSD) using a time window of chosen length in each recording. The length and overlap of FFT window for 1s and 50 %, respectively; Correction of PSD for recorder's sensitivity, gain and, when possible, frequency response to get the output data in sound pressure units (µPa). Figure 8 shows PSD of a 1 h recording started at 05:29:35 (UTC) collected by the glider when the earthquake occurred at 05:29:32 (UTC). It contains underwater earthquake sound as well as the operating noise of the glider's machinery mentioned in Section 3, e.g., piston and pitch motors, fin, airbladder, altimeter, pumped CTD. The glider descended with pumped CTD and altimeter turned on until it reached the turning depth (200 m), and it then turned on the piston and pitch motor for ascending and inflated the external air bladder to begin surfacing for communication. The pumped CTD and altimeter were turned off when the glider was ascending. The earthquake sounds were dominated below about 50 Hz.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 8 of 13 recording. The length and overlap of FFT window for 1s and 50 %, respectively; Correction of PSD for recorder's sensitivity, gain and, when possible, frequency response to get the output data in sound pressure units (μPa). Figure 8 shows PSD of a 1 hour recording started at 05:29:35 (UTC) collected by the glider when the earthquake occurred at 05:29:32 (UTC). It contains underwater earthquake sound as well as the operating noise of the glider's machinery mentioned in section III e.g., piston and pitch motors, fin, airbladder, altimeter, pumped CTD. The glider descended with pumped CTD and altimeter turned on until it reached the turning depth (200 m), and it then turned on the piston and pitch motor for ascending and inflated the external air bladder to begin surfacing for communication. The pumped CTD and altimeter were turned off when the glider was ascending. The earthquake sounds were dominated below about 50 Hz. Spectrum and spectrogram analyses of the earthquake sound have been carried out. By carefully examining the frequency spectrum of the record in Figure 9b, it can be seen that the dominant frequency of the earthquake sounds was about 10 Hz. Figure 9c shows that after the earthquake occurred, the average ambient noise level at frequencies below 50 Hz was approximately 60 dB higher than the ambient noise level before the earthquake occurred. The earthquake acoustic signal can be distinguished from glider's self-noises of fin, pumped CTD and altimeter. The fin noise exists at frequencies of 800-950 Hz every 3 to 4 seconds. The pumped CTD sensor noise exists at discrete frequencies of 357, 410, 710, 1450 Hz. The dominant CTD noise frequency is 710 Hz. The altimeter noise exists at frequencies above 100 Hz. Spectrum and spectrogram analyses of the earthquake sound have been carried out. By carefully examining the frequency spectrum of the record in Figure 9b, it can be seen that the dominant frequency of the earthquake sounds was about 10 Hz. Figure 9c shows that after the earthquake occurred, the average ambient noise level at frequencies below 50 Hz was approximately 60 dB higher than the ambient noise level before the earthquake occurred. The earthquake acoustic signal can be distinguished from glider's self-noises of fin, pumped CTD and altimeter.

Correlation with Situ Seismic Waves and Earthquake Sounds
Earthquake generated acoustic waves in the ocean are referred to as T-wave because in certain conditions, when they reach the shore, they may convert back to seismic waves and arrive third after the P-and S-seismic waves on near-shore seismological stations. The time series of the seismic and earthquake acoustic signal detected at the five stations and the underwater glider are presented in Figure 10. The time of arrival of the seismic waves observed at the station is proportional to the distance from the epicenter. Since the underwater glider was closer than Chungju station to the epicenter, the arrival time of the T-wave (i.e., earthquake acoustic signal) observed by the underwater glider was earlier than that of the seismic waves observed at the Chungju station. The seismic P-and S-waves travel at velocities from 2000 to 7000 m/s in the crust just up to the seafloor bottom (about 430 m) of detection point, then in the water column and finally detected by the underwater glider (about 200 m). It is clearly shown that the acoustic signal detected by the glider from the Pohang earthquake is a Tertiary wave (T-wave) in Figure 11.

Correlation with Situ Seismic Waves and Earthquake Sounds
Earthquake generated acoustic waves in the ocean are referred to as T-wave because in certain conditions, when they reach the shore, they may convert back to seismic waves and arrive third after the P-and S-seismic waves on near-shore seismological stations.
IRIS (Incorporated Research Institutions for Seismology) archives waveform (time-series) data from stations around the world. Seismic wave data of the Pohang Mw 5. The time series of the seismic and earthquake acoustic signal detected at the five stations and the underwater glider are presented in Figure 10. The time of arrival of the seismic waves observed at the station is proportional to the distance from the epicenter. Since the underwater glider was closer than Chungju station to the epicenter, the arrival time of the T-wave (i.e., earthquake acoustic signal) observed by the underwater glider was earlier than that of the seismic waves observed at the Chungju station. The seismic P-and S-waves travel at velocities from 2000 to 7000 m/s in the crust just up to the seafloor bottom (about 430 m) of detection point, then in the water column and finally detected by the underwater glider (about 200 m). It is clearly shown that the acoustic signal detected by the glider from the Pohang earthquake is a Tertiary wave (T-wave) in Figure 11. Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 13

Discussion
Early models of the Slocum underwater glider were equipped with a non-pumped CTD. The non-pumped CTD has a large salinity uncertainty because the flow rate in the conductivity cell is not

Discussion
Early models of the Slocum underwater glider were equipped with a non-pumped CTD. The non-pumped CTD has a large salinity uncertainty because the flow rate in the conductivity cell is not

Discussion
Early models of the Slocum underwater glider were equipped with a non-pumped CTD. The non-pumped CTD has a large salinity uncertainty because the flow rate in the conductivity cell is not constant. The salinity measurement uncertainty of a non-pumped CTD was about 30 times greater than that of a pumped CTD [27]. To obtain good quality ocean observation data, recently-produced Slocum underwater gliders are equipped only with a pumped CTD. In terms of the overall operating time of the underwater glider, the pumped CTD noise is higher and it occurs for a longer time than noise from other sources. It is the main noise source in underwater glider operation.
When an acoustic source emits in a test tank, acoustic waves are reflected by the walls, the surface, and the bottom. The multiple reflections give rise to a reverberated field. Measuring underwater sound power in reverberant environment such as a tank has advantages of speed but is limited in accuracy and frequency range, but it has been used to good effect in specific applications. For small sources (or small noisy components) like underwater glider, it is possible to make measurements of radiated noise in reverberant tanks [28].
For the seismic observation of seafloor and coastal earthquake disaster risk assessment, a comprehensive analysis of seismic sources is needed. This requires simultaneously observing underwater phenomena caused by seafloor earthquakes. Moreover, submarine earthquake underwater acoustical observation technology using underwater gliders is very synergistic because it can provide basic research data on natural, anthropogenic, and ecological underwater sources in the ocean. In addition, underwater sound observations using underwater gliders can be of great assistance in developing earth observation systems by providing submarine platforms ensuring stable operational technology and linking with existing earthquake observation networks.
To this end, improving the method for estimating a glider's underwater trajectory is important because it is a challenging task due to several errors in the linear interpolation method. It is not realistic because the position is estimated by dividing the distance between the real GPS fixed position by transit time, and the glider's direction and speed remain constant. In other words, the trajectory looks like a straight line. Our future work aims to provide a better method for estimating a glider's optimal underwater trajectory, focusing on the importance of the ocean currents and underwater wireless communication technology to transmit the acoustic data observed by the underwater glider to the shore in real time through the OSG such as WaveGliders in [20].

Conclusions
In this study, the acoustic sounds of earthquake occurred on 15 November 2017 in the Pohang area of South Korea observed by underwater glider were analyzed. By analysis of time correlation with seismic waves detected by five seismic land stations and the earthquake acoustic signal, it is clearly shown that the seismic waves converted to Tertiary waves (i.e., earthquake acoustic signal) and then it detected by the underwater glider. In order to distinguish the earthquake sound from the glider's operation noise, Sound Pressure Level (SPL) of the noise sources associated with the driving and sensor parts of the underwater glider were analyzed and measured at laboratory tank and sea. The earthquake acoustic signal in the range of 10-50 Hz was distinguished from glider's self-noises which are exist over 100 Hz. The results allow constraining the acoustic sound level of the earthquake and suggest that the glider provides an effective platform for enhancing the earth seismic observation systems and monitoring natural and anthropogenic ocean sounds.