Next Article in Journal
Design of a Tri-Band Wearable Antenna for Millimeter-Wave 5G Applications
Next Article in Special Issue
Communication System Based on Magnetic Coils for Underwater Vehicles
Previous Article in Journal
Colorful Image Colorization with Classification and Asymmetric Feature Fusion
Previous Article in Special Issue
Energy Harvesting for TDS-OFDM in NOMA-Based Underwater Communication Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Long-Range Biomimetic Covert Communication Method Mimicking Large Whale

Agency of Defense Development, Changwon-si 516852, Korea
*
Author to whom correspondence should be addressed.
Submission received: 25 August 2022 / Revised: 17 October 2022 / Accepted: 19 October 2022 / Published: 20 October 2022
(This article belongs to the Special Issue Underwater Wireless Communications)

Abstract

:
Short-range biomimetic covert communications have been developed using dolphin whistles for underwater acoustic covert communications. Due to a channel characteristics difference by range, the conventional short-range methods cannot be directly applied to long-range communications. To enable long-range biomimicking communication, overcoming the large multipath delay and a high degree of mimic (DoM) in the low-frequency band is required. This paper proposes a novel biomimetic communication method that preserves a low bit-error rate (BER) with a large DoM in the low-frequency band. For the transmission, the proposed method utilizes the time-dependent frequency change of the whistle, and its receiver obtains additional SNR gain from the multipath delay. Computer simulations and practical ocean experiments were executed to demonstrate that the BER performance of the proposed method is better than the conventional methods. For the DoM assessment, the novel machine learning-based method was utilized, and the result shows that the whistles generated by the proposed method were recognized as the actual whistle of the right humpback whale.

1. Introduction

Recently, a biomimetic communication method has been studied for military purposes. Biomimetic communication methods make the enemy confuse our communication signals with dolphin or whale sounds to achieve covertness [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. Traditional covert communication systems decrease the signal strength and hide the existence of the signal using various techniques (e.g., spread spectrum, chaotic modulation, etc.) [15,16,17,18]. Low-power spectral density of the transmitted signal causes a low SNR at the receiver, which results in low bit-error-rate (BER) performance. Especially for long-distance communication, low frequency and narrow bandwidth are inevitable, and it is difficult to conceal the signal using these methods [7,8,9,10,11,12,13,14]. Biomimetic underwater acoustic communication has been developed to overcome the disadvantage of traditional covert communication. This method allows high transmit power without considering a low detection probability, and large SNR at the receiver can be preserved.
Conventional biomimetic communication methods utilize various modulation methods such as chirp spread spectrum (CSS), frequency shift keying (FSK), phase shift keying (PSK), and time-frequency shift keying (TFSK) to modulate bits into whistles [7,8,9,10,11,12,13,14]. Conventional biomimetic communication studies have demonstrated their communication performance through ocean experiments. In addition, the conventional methods have evaluated the degree of mimic (DoM) by comparing the shape of mimicking the signal with the real whistle. Especially in studies of CV-CFM (based on PSK) and TFSK, the authors performed the mean-of-score (MOS) test and showed that humans could hardly distinguish a mimicking signal from a real dolphin whistle [11,12,13,14]. All of the conventional studies mimicked the whistle of small dolphins (such as white-sided dolphins utilizing tens of kilohertz sounds) and conducted ocean experiments at a short range within 5 km [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. However, considering the operation area of military-purpose underwater vehicles from tens to hundreds of km, biomimetic communication should support long ranges. For long-distance communication, using a low-frequency signal (i.e., less than several hundred hertz) is advantageous over a high-frequency signal (i.e., over tens of kilohertz) due to absorption loss, and large whales use these low-frequency whistles. Thus, it is necessary to mimic a large whale for long-range biomimicking communication, and long-range underwater acoustic communication is difficult with conventional biomimetic methods.
The first problem is an increase in BER due to a large multipath delay. As the range between transmitter and receiver increases to tens of kilometers, the multipath delay also increases up to a few seconds [19,20,21,22,23,24,25]. However, conventional biomimetic communication methods do not consider a few seconds of multipath delay because these methods are designed for a short communication range within 5 km. Thus, when a multipath delay occurs for a few seconds, the BER of the conventional biomimetic methods increases. The second problem is a degradation of the DoM performance. Human psychoacoustic ability is the most sensitive to sound in hundreds of hertz to kilohertz [26]. On the other hand, small dolphins utilize high-frequency whistles exceeding the human auditory sense. When the conventional methods mimic small dolphins, these facts have an advantage in the DoM performance. As a result of the DoM evaluation in Ref. [14], it was determined that the DoM of killer whales (2 kHz~4 kHz) is lower than that of white-sided dolphins (15 kHz~25 kHz). Therefore, two things are required to enable long-range communication with biomimetic modulation. The first one is the robustness against large multipath delay, and the second is the performance improvement of the DoM for whistles with hundreds of hertz to kilohertz. In this paper, we propose a novel long-range biomimetic modulation method that has robustness against large multipath delay and higher DoM performance than the conventional methods.
Since the whistle is a non-linear chirp signal, it has strong robustness against the distortion of large multipath delay. Thus, if the bit can be modulated while maintaining the shape of the whistle, the modulated signal has tolerance to the large multipath delay and the high DoM performance. Whales generate whistles with various frequency contours, but they have a habit when generating the sounds of the whistle. The whistle-generating habit is that the following frequency values are limited by the frequency value of the previous time. The proposed biomimicking modulation method utilizes this frequency-change habit of the whistle.
The main contributions of the paper are summarized as follows:
  • We modeled the time-dependent frequency change of the whistle of the right humpback whale (RHW) using the Markov chain (MC). The modeled results can be employed to select the frequency values used by the whistle. Note that we considered the time-frequency orthogonality of the symbol to preserve low BER when modeling the frequency change of the whistle.
  • The proposed modulation method preserves the signal shape of the whistle because it modulates bits using frequency values that whales use to generate a whistle. Thus, the whistle-mimicking communication signal generated by the proposed modulation method shows improved performance of the DoM and tolerance of large multipath delay.
  • The proposed demodulation method estimates the multipath delay profile using a preamble and increases received signal gain using the estimated profile. Thus, the proposed method achieves low BER in a long-range underwater environment with a large multipath delay.
  • Computer simulations and practical ocean experiments were conducted and demonstrated that the proposed method had a lower BER than conventional covert communication methods.
  • Since the proposed modulation method makes a whistle using modeling results, a machine learning-based DoM assessment was conducted. The assessment results show that the trained machine learning classifier recognized the whistle-mimicking signal generated by the proposed modulation method as the whistle of the RHW. This result shows that machine learning can be used as an effective evaluation method for the DoM performance of biomimetic communication.
This paper is organized as follows. Section 2 describes how to model the time-dependent frequency change of the whistle of the RHW using MC. In Section 3, the proposed biomimicking modulation and demodulation methods are described. Section 4 shows the result of a machine learning-based DoM assessment. In Section 5, through computer simulations and ocean experiments, it is shown that the proposed biomimicking communication method has a lower BER than conventional methods. Section 6 concludes the paper.

2. Modeling Whistles of the Right Humpback Whale

Whales communicate with each other using whistles. For the long-range biomimicking communication method, we selected the RHW as a mimicking model. The RHW is about 15 m long and makes sounds from tens to hundreds of hertz. Figure 1 shows the spectrogram of the whistles of the RHW.
Figure 1 shows three frequency contours of the whale whistle starting at 150 Hz. In Figure 1a,d, the frequency of the whistle is 150 Hz until halfway through the whistle. Then it increases to 250~300 Hz. In Figure 1b,e, the whistle frequency increases linearly from 150 Hz to 300 Hz. The whistles in Figure 1c,f also start at 150 Hz but end at 200 Hz. As shown in Figure 1, whales generate various whistles, but we can see that there is a specific habit when whales make whistle sounds. The frequency values of the whistle are limited according to their habits [27,28,29,30]. For example, if the RHW starts the whistle at 150 Hz, we can infer that whales are used to generating sounds under 300 Hz halfway through the whistle. In this paper, we propose a biomimicking communication method using the whistle-generating habit of whales.
One of the simple examples of the proposed method is shown in Figure 1d,e. The whistles in Figure 1d,e start at 150 Hz. Then, at 0.55 s of whistle time, bit 0 is modulated when the frequency value of the whistle is 150 Hz (red box). If the frequency value of the whistle is 200 Hz, it indicates bit 1 (white box). Likewise, if the whistle-generating habit of the RHW is modeled, we can transmit binary information according to the modeling results. As shown in Figure 1, pre-recorded whistles are already contaminated with noise, and communication performance cannot be guaranteed when the pre-recorded whistle is utilized directly. Therefore, in this paper, a whale whistle is modeled.
The whistle-generating habit of the whale is that the subsequent frequency values are restricted by the frequency value of the previous time [27,28,29,30]. This habit is suitable for Markov chain (MC)-based modeling, and some studies have analyzed the time-dependent frequency change of whistles using MC [31,32,33,34,35,36,37]. However, applying conventional modeling methods directly to the proposed methods is difficult due to the degradation of BER performance. The time-frequency orthogonality between symbols must be maintained for the low BER. However, conventional modeling methods do not consider this orthogonality. Therefore, the time-dependent frequency change of the whistles is modeled in a time-frequency unit that satisfies the orthogonality. The proposed method represents information bits using the frequency change with time. The symbol can be distinguished when the time and frequency change is greater than twice the reciprocal of the symbol duration, such as linear frequency modulation (LFM).
For the modeling, the frequency values of a whistle over time were extracted by short-time Fourier transform (STFT). When the time window length of STFT is Δ T , the minimum frequency unit that guarantees the orthogonality between two different LFMs is 2 / Δ T . Thus, the frequency resolution of STFT Δ F is 2 / Δ T . As a result of STFT, the whistle is represented in the 2D domain with K -length of the time sample and L -length of the frequency sample. Figure 2a shows the ( L × K ) size of the 2D whistle. In Figure 2,   T w is the time length of the whistle. We can obtain one frequency value for each K column from the 2D whistle, and the time-dependent frequency change is obtained as F = [ f 1 , , f k , , f K ] . If we obtain the F matrix from N whistles, the probability distribution of the frequency change over time can be obtained using Markov-chain modeling.
It is assumed that there are M whistles with a frequency value of Δ F × l at an arbitrary time Δ T × k . Among these M whistles, only M whistles have the increased frequency value Δ F × ( l + 1 ) at the next time Δ T × ( k + 1 ) . When the frequency value of the whistle is Δ F × l at Δ T × k , the probability that the frequency value is Δ F × ( l + 1 ) in the next time interval Δ T × ( k + 1 ) becomes M / M . This conditional probability can be expressed as P (   f k = Δ F × l |   f k + 1 = Δ F × ( l + 1 ) ) . If this conditional probability is 0, the whale does not use the whistle pattern that changes from Δ F × l to Δ F × ( l + 1 ) . By applying the Markov chain to the calculated conditional probabilities, we can model the time-dependent frequency-change rule of the whistle of the RHW. As a result of modeling, when the frequency value is   f k 1 at Δ T × k , the set of the frequency values generated by the whale at Δ T × ( k + 1 ) is expressed as Equation (1):
F k + 1 = {   Δ F × l   | Pr (   f k |   f k + 1 = Δ F × l ) > Thr ,   l = 1 ~ L ,   f k F k   } ,   k = 1 ~ K .
In Equation (1), F k + 1 is the set of frequency values used by a whale at Δ T × ( k + 1 ) . Thr is a threshold value for selecting only frequency values usually generated by whales, and we set it as 0.5. Note that frequency values in F k + 1 are quantized by Δ F , which guarantees symbol orthogonality. On the other hand, the actual whale whistle has time-frequency ambiguity, as shown in Figure 1. The time-frequency ambiguity of the whistle is also considered by modeling the time-dependent frequency change of the whistle with various lengths of the STFT window. Figure 2b shows the STFT results with another window length.
In Figure 2, Δ T is smaller than Δ T ^ . For example, let us assume that Δ T ^ (e.g., Δ T ^ : 0.2) is 2 Δ T (e.g., Δ T : 0.1). The time resolution of Figure 2a is higher than that of Figure 2b. However, the frequency resolution of Figure 2a is lower than that of Figure 2b. As shown in Figure 2, the time-frequency ambiguity of the whistle is considered by conducting STFT with various lengths of the window. In this paper, we modeled the whistle-generating habit using 6700 whistles of the RHW [38] and six lengths of the STFT window (i.e., 0.075 s, 0.1 s, 0.125 s, 0.15 s, 0.175 s, 0.2 s). Figure 3 shows some examples of the modeling results.
In Figure 3, blue lines show the possible variation of whistle frequency. As shown in Figure 3, whistle lengths are quantized to the mean value of each time length of T w . Since Δ T of Figure 3a is shorter than that of Figure 3b, whistles are divided into three parts in (a) and two parts in (b). In Figure 3, the RHWs are used to start within 130~165 Hz when generating a whistle. Thus, in Figure 3b, starting frequency values ( F 1 ) are quantized as (135 Hz, 145 Hz, 155 Hz) and frequency values at the half-time of the whistle ( F 2 ) are quantized as (145 Hz, 155 Hz, 165 Hz, 175 Hz, 185 Hz). As a result of the model in Figure 3b, the RHW does not generate exceeding 165 Hz at half of the whistle when the start frequency of the whistle is 135 Hz. The proposed modulation method utilizes this whistle-generating habit of the RHW to modulate bits. A simple example is represented in Figure 3b. When the whistle starts at 135 Hz, there are two available frequency values (145 Hz and 155 Hz) at 0.2 s. In this case, 145 Hz represents 1 bit and 155 Hz is 0 bits. In Figure 3b, the given parameters satisfy the time-frequency orthogonality condition ( Δ F 2 / Δ T ) and low BER can be achieved. Section 3 describes the proposed method in detail.

3. The Proposed Biomimicking Communication Method

The proposed method utilizes the model of the whistle-generating habit of the RHW when it modulates/de-modulates bits. The transmitter and receiver share these modeling results. The proposed transmitter uses a pseudo-random number generator and randomly selects one model among the models with various lengths of the window. It is assumed that the receiver also selects the same model using the same pseudo-random number generator. Firstly, the proposed biomimicking modulation method is described.

3.1. The Proposed Biomimicking Modulation Method

The time-frequency units of the model ( Δ T and Δ F ) satisfy symbol orthogonality and become the symbol time frequency. The whistle is divided by the symbol length Δ T . Let us assume that whistle is divided into K symbols. Information bits are mapped to the available frequency values at each symbol period. Figure 4 shows an example. In the ( k 1 ) th symbol period, there are four available frequency values and two bits are modulated. If 10 is modulated to the ( k 1 ) th symbol, there are two available frequencies in the k th symbol period. Thus, the   k th symbol modulates one bit, whereas the ( k 1 ) th symbol modulates two bits. That is, the number of modulation bits in the k th symbol depends on the frequency value selected in the ( k 1 ) th symbol. M k is the modulation order in the k th symbol.
In Figure 4, the proposed method selects the available whistle frequency according to the bits using gray coding and generates the entire biomimicking whistle using the selected frequency. The modulation process is described in detail. When k is 1, available frequency values are elements of the set F 1 in Equation (1) and the modulation order M 1 is calculated as l o g 2 ( | F 1 | ) . Thus, M 1 bits are transmitted with the first symbol. Assume that B M 1 is the set of all M 1 bit sequences. Each element of B M 1 is mapped one by one to an element in a set F 1 using gray coding. When all the bit sequences in B M 1 are converted to decimals, the converted result is a set of integers from 0 to 2 M 1 1 . This set of integers becomes the symbol set ( S 1 ) for the first symbol and each element of S 1 is an index of F 1 . Let us assume that b 1 is the Tx bits sequence and it is converted to the symbol s 1 ( s 1 S 1 ). The whistle frequency of the first symbol is selected with F 1 ( s 1 ) . When k is 2 or more, the available frequency in the k -th symbol depends on the selected frequency in the ( k 1 ) -th symbol. Therefore, all frequencies in the set F k cannot be used, as shown in Figure 4. When G k is the set of frequency values available in the k -th symbol and G k is expressed by Equation (2), then
{ k 2 , G k = { Δ F × l | Pr ( G k 1 ( s k 1 ) = f k 1 | f k = Δ F × l ) > T h r , l = 1 ~ L } k = 1 , G 1 = F 1 .
In Equation (2), when k is 1, any frequency in the set F 1 is available. Thus, G 1 is the same as F 1 . If k is greater than or equal to 2, a frequency value is selected from the set G k according to the Tx bits b k and the modulation order M k is calculated as l o g 2 ( | G k | ) . Therefore, likewise, when k is 1, the whistle frequency for the k -th symbol is selected as G k ( s k ) and the biomimicking communication signal in the k -th symbol period is expressed by Equation (3):
W k ( t ) = cos ( 2 π ( ( G k 1 ( s k 1 ) ) G k ( s k ) ) 2 Δ T t 2 + G k 1 ( s k 1 ) t ) ) .
In Equation (3), the biomimicking communication signal in the k -th symbol period is LFM, as shown in Figure 4. The entire biomimicking communication signal is the sum of all the signals from 1 to the K -th symbol period, and it is obtained by Equation (4):
W ( t ) = k = 2 K W k ( t Δ T k ) .
The block diagram of the proposed biomimicking modulation is shown in Figure 5.
A preamble is inserted in front of the biomimicking communication signal to find the received signal and estimate the channel delay at the receiver. This preamble is the original whistle of the whale. Let us assume that the receiver also has the same preamble. The frame structure of the biomimicking communication signal is shown in Figure 6.
The proposed biomimicking whistle is the summation of several LFMs. Therefore, like LFM, the auto-correlation function becomes a Kronecker delta function. If the channel delay can be estimated, an additional signal gain can be obtained using signals received through multiple paths, reducing bit error. The following section describes the reception method.

3.2. The Proposed Biomimicking Demodulation Method

The transmitted signals generated using Equation (4) pass through the underwater channel and are received at the receiver. The received signal y ( t ) is modeled as
y ( t ) = i = 1 I H i × W ( t τ i ) + n ( t ) .
In Equation (5), H i is the channel gain of the i -th path and n ( t ) denotes an additive white Gaussian noise (AWGN). The number of the path is I . The receiver generally finds the received signal on the path with the greatest channel gain using the preamble. Thus, we assume that τ 1 is zero and H i has the maximum value when i is 1. For the precise detection of received bits, a maximum likelihood sequence detection (MLSD)-based receiver is proposed. In Equation (2), the k -th symbol is modulated depending on the ( k 1 ) -th symbol. Therefore, ML detecting for the sequence with K symbols shows a lower error rate than detecting each of the K symbols as a single one. When the sequence of the transmitted K symbols is S Tx * = { s 1 * , , s k * , , s K * } , detecting the sequence of symbols with the maximum likelihood probability for S Tx * is expressed as
S ^ Rx = max S 1 , , S K P (   G 1 ( s 1 * ) =   G 1 ( s 1 ) , , G k ( s k * ) = G k ( s K ) , , G K ( s K * ) =   G K ( s K ) ) = max S 1 , , S K k = 2 K P (   G k 1 ( s k 1 * ) =   G k 1 ( s k 1 ) |   G k ( s k * ) =   G k ( s k ) ) .
To find a symbol sequence that satisfies Equation (6), the proposed demodulation method calculates a correlation between the transmitted whistle and all whistles that the model can generate. When the generated whistle is W ^ ( t ) with the given model and K symbols { s 1 , , s K }, the correlation between the generated whistle W ^ ( t ) and transmitted whistle W ( t ) is expressed as
R ( τ ) =   W ( t τ ) × W ^ ( t ) d t = k = 2 K   W k ( t ) × W ^ k ( t ) d t   = k = 2 K R k ( τ ) .
As derived from Equation (7), the correlation between W ^ ( t ) and W ( t ) is equal to the summation of K symbol correlation values R k ( τ ) that are calculated with the transmitted and generated whistle for each symbol. Therefore, if the symbol sequence satisfies Equation (6), the signal generated by the corresponding sequence has the maximum correlation with the received signal. The correlation between y ( t ) and W ^ ( t ) is expressed as Equation (8).
y ( t ) W ^ ( t ) = i = 1 I H i × R ( τ τ i ) + n ( τ ) ,
where   denotes a correlation operation. As mentioned above, when W ^ k ( t ) is equal to W k ( t ) , R ( τ ) is Kronecker delta function δ ( τ )   and the result of Equation (8) is the same as the delay profile. Note that H i has the maximum value when i is 1. Thus, the maximum value of Equation (8) is the gain of the first path H 1 , whereas the proposed demodulation method obtains an additional signal gain from the received signal through the other paths. Let us assume that the preamble has gone through the same multipath as the received signal. The estimated channel using the preamble is i = 1 I   H i × δ ( τ τ i ) . To obtain the additional signal gain, integrate and dump the estimated channel and Equation (8). The result of the integration and dump is expressed as follows:
{ H 0 ( W k ( t ) = W ^ k ( t ) ) : P = i = 1 I H l 2 + N H 1 ( W k ( t ) W ^ k ( t ) ) : P = i = 1 I H l 2 ×   R ( τ ) d τ + N .
If W ^ k ( t ) is equal to W k ( t ) , Equation (9) is the summation of I path gain, whereas when W ^ k ( t ) and W k ( t ) are different, R ( τ ) d τ has a very small value because R ( τ ) becomes the cross-correlation of different LFM signals. Therefore, Equation (9) becomes the maximum value when W ^ k ( t ) is equal to W k ( t ) and the received symbol sequence is detected as Equation (10):
S ^ Rx = max S 1 , , S K P .
The block diagram of the proposed biomimicking demodulation method is depicted in Figure 7.

4. Evaluation of Degree-of-Mimic Performance

This section evaluates the degree-of-mimic (DoM) performance of the proposed method. Quantitative (i.e., spectral correlation) and qualitative methods (i.e., MOS-test) have been proposed to measure the degree of mimic in conventional biomimetic studies [8,9,10,11,12,13,14]. Spectral correlation measures the similarity of the time-frequency change between the actual whale whistle, and the MOS test evaluates whether a person recognizes a mimicking whistle as an actual whale whistle. The results of conventional studies show that the biomimicking whistle has the best score in both DOM evaluation methods (i.e., spectral correlation and MOS test) when it is not distorted during the modulation process [10,11]. Since the proposed method does not distort the whale whistle, it has the highest score of the two methods and is undetected by conventional communication signal-recognition techniques, which find the artificial and periodical signal characteristics. Therefore, this paper proposes a novel DoM evaluation method to more precisely evaluate the DoM of biomimicking signals.
The proposed DoM evaluation method discriminates even if the biomimicking whistle is the sound of a specific whale species. The proposed method not only judges whether a biomimicking whistle sounds like a whale, but also determines whether it simply sounds like a specific species (e.g., right humpback, killer, etc.). Determining a species of whale with only whistles is a professional task and requires much experience. Thus, to complete this task effectively, whistle classifiers based on machine learning have been studied [39,40]. In this paper, the machine-learning classifier was utilized to evaluate the DoM performance of the proposed method. The classifier that won the Marine Explore and Cornell University Whale Detection Challenge competition was used for evaluation [38]. To evaluate the DoM, 1000 biomimicking whistles were generated using the model in Section 3. These 1000 whistles were modulated with conventional methods (i.e., CV-CFM and TFSK [10,13] and FSK [6,7,8]) to compare the DoM performance of the proposed method with that of conventional methods. TFSK-modulated biomimicking whistles are mathematically the same signal as CV-CFM-modulated whistles [13]. Thus, these two modulation schemes shared whistles for evaluation. Note that all models with six window lengths were utilized to generate 1000 whistles. Figure 8 shows examples of the whistles generated for the DoM evaluation.
The whistle in Figure 8 was generated by the model with a 0.2 s window length. The results of the DoM evaluation are shown in Table 1.
In Table 1, 95% of whistles generated by the proposed method were recognized as the whistle of the RHW. Note that the classifier did not simply identify the given sound as an actual whistle of a whale but rather as the whistle of the RHW. Thus, this result demonstrates that the proposed model generated a whistle like the real RHW whistle. However, in cases of FSK- and CV-CFM-modulated whistles, only 10% and 60% were recognized as the whistle of the RHW, respectively. Note that these three types of whistles have the same frequency contour; the only difference is the modulation scheme. Therefore, this result shows that the conventional methods distorted the whistle when modulating bits. This distortion decreased the DoM performance. In addition, as seen in Figure 8b,c, the existing method distorted the signal so humans could recognize this difference and the MOS score decreased. The value of spectral correlation also decreased. As a result of the DoM evaluation, the DoM performance of the proposed method is superior to that of the conventional methods. In the next section, the computer simulations and ocean-experiment results are shown for communication-performance comparisons of the proposed method with the conventional methods.

5. Evaluation of Communication Performance

5.1. Simulation

The BER performance of the proposed biomimicking method was compared with that of conventional methods (i.e., FSK, CV-CFM, and TFSK) in the long-range underwater environment. The modulation order and parameter were selected for the fair BER comparison to yield a similar data rate. These values are given in Table 2. M is the modulation order.
In Table 2, since all models with six window lengths were utilized for the proposed method, the data rate is an average value. In Table 2, let N be the number of transmitted whistles for simulation. b n is the number of bits transmitted by the n -th whistle and T n is the time length of this whistle. The data rate is obtained as ( n = 1 N b n / T n ) / N . The data rate of the conventional methods is about 12~15 bps, except for TFSK. The simulations were performed with the bandwidth used by the whale, and the data rate of the TFSK is limited.
There are many simulators for modeling underwater acoustic sound propagation. However, we referred to the measured channel from the experiments of long-range (i.e., longer than 50 km) underwater acoustic communication for realistic simulation. Table 3 summarizes the cases of the long-range underwater acoustic communication experiment conducted in the bandwidth of the RHW. The depiction of the channel measured in the experiments is also described in Table 3.
In Table 3, the maximum channel delay is 2.5 s and the maximum number of the multipath is 17. Based on these facts, three channels were generated by the bellhop according to maximum delay (i.e., 30 ms, 1 s, and 2 s), and these three channels were utilized for simulation. The three channels used for simulation and the BER results are displayed in Figure 9.
In Figure 9, the orange line denotes the BERs of the proposed method. The blue line is the BER of the proposed method without multipath gain. The yellow, purple, and green indicate the BERs of the conventional biomimicking communication methods (i.e., CV-CFM, FSK, TFSK, respectively). In Figure 9, the proposed method showed better BER performances than the conventional algorithms. The conventional methods had error floors at 10−2~10−1. When SNR was lower than −10 dB, TFSK showed the lowest BER. However, the data rate of the TFSK was also half that of the others. When proposed without multipath gain, the BER performance of the proposed method was similar to that of TFSK. Please note that, as shown in Table 2, the data rate of the proposed method was higher than that of TFSK. Therefore, the BER performance of the proposed method was better even without multipath combining. In addition, increasing SNR did not decrease the BER of the TFSK. The proposed method decreased BER by taking multipath gain. These results show that it is difficult for conventional methods to overcome the large multipath delay in the long-range underwater communication environment. However, the proposed method obtained the additional SNR gain using signals received from multipath and showed the lowest BER with 1~2 s multipath delay, as depicted in Figure 9b,c. In the following subsection, the results of the practical ocean experiments are described.

5.2. Ocean Experiment

Practical ocean experiments with the proposed and conventional methods were executed 60 km away, eastbound from Pohang city, on 18 May 2022. The range between the transmitter and the receiver was over 60 km. The average depth of the ocean was about 1100 m. The transmitter was located at 295 m with a depth of 850 m. The bandwidth of the transmitter was 2.6 kHz to 4.6 kHz, and we upconverted the center frequency of the whistle to 3 kHz. The three-channel receiver was utilized, and each sensor was positioned at 205 m, 230 m, and 255 m below the ocean surface. The depth of the receiver was about 1200 m. Figure 10 shows the experiment location and the configuration of the experiments.
The parameters of the proposed and conventional methods were the same as those of computer simulations in Table 2. When bits are utilized as orthogonal code, CV-CFM-modulated whistles become TFSK-modulated biomimicking whistles [14]. Thus, these two modulation schemes shared whistles for evaluation. The DoM performance of the FSK biomimicking signal showed about 10%, which is too low to use practically. Thus, the FSK biomimicking method was not transmitted during the ocean experiment. For the BER performance calculation, 945 whistles were transmitted for each method. The proposed and conventional methods (i.e., CV-CFM and TFSK) were alternatively sent and went through the same underwater acoustic channel. Figure 11a,b depict the multipath delay profile of the proposed method and the conventional methods.
In Figure 11, the same multipath delays were measured for each method. It can be seen that more than two paths occurred between 30 ms and 40 ms. The three sensors were 25 m apart from each other, and the depths of the three sensors were different. Therefore, the received signal of each sensor went through a different underwater communication channel, as shown in Figure 11. The uncoded BERs of each sensor are shown in Table 4.
The average BER of the proposed method was 0.0057, whereas the conventional method exhibited 0.079 and 0.1165. The BER of the proposed method without the multipath was 0.009, which is also lower than conventional methods. This result shows that the proposed method obtained gains from the multipath. In the ocean experiments, the SNR of the received signal was estimated at about 5 dB, and the multipath delays of the second and third sensors were similar to those of the simulation experiment in Figure 9a. Therefore, the BERs of those two sensors were similar to the BER in Figure 9a. As shown in Figure 11, the multipath of the first sensor was the largest, so the BER of the first sensor increased like the simulation results in Figure 9. Table 4 demonstrates that the proposed method showed the lowest BER. This result shows that the proposed method can practically overcome the multipath delay.

6. Conclusions

This paper proposes a biomimicking modulation method for long-range underwater acoustic communication. The proposed method mimicked the whistle of a large-size whale (i.e., a right humpback whale) and modeled the time-dependent frequency change of the whistles. This model was utilized to overcome the large multipath delay in long-range underwater acoustic communication. Results of the computer simulations and the ocean experiment demonstrate that the BER performance of the proposed method is better than that of conventional studies. In this paper, the novel DoM assessment based on a machine-learning classifier was also utilized to demonstrate that whistles generated by the proposed method are the same as those of the RHW. The assessments showed that the whistles generated with the proposed method were determined as a natural RHW sound, whereas the conventional methods were not. Thus, the proposed method shows better BER and DoM performance than that of conventional methods.

Author Contributions

Conceptualization, J.A. and W.K.; methodology, J.A. and W.K.; software, J.A.; validation, J.A.; formal analysis, J.A.; investigation, J.A. and W.K; resources, J.A. and W.K; data curation, J.A., D.D. and W.K.; writing—original draft preparation, J.A.; writing—review and editing, J.A., D.D. and W.K.; visualization, J.A.; supervision, J.A., D.D. and W.K.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, S.; Qiao, G.; Ismail, A. Covert underwater acoustic communication using dolphin sounds. J. Acoust. Soc. Am. 2013, 133, EL300–EL306. [Google Scholar] [CrossRef] [Green Version]
  2. Han, X.; Yin, J.; Du, P.; Zhang, X. Experimental demonstration of underwater acoustic communication using bionic signals. J. Appl. Acoust. 2013, 78, 7–10. [Google Scholar] [CrossRef]
  3. ElMoslimany, A.; Zhou, M.; Duman, T.M.; Papandreou-Suppappola, A. An underwater acoustic communication scheme exploiting biological sounds. J. Wirel. Commun. Mob. Comput. 2016, 16, 2194–2211. [Google Scholar] [CrossRef]
  4. Liu, S.; Wang, M.; Ma, T.; Qiao, G.; Bilal, M. Covert underwater communication by camouflaging sea piling sounds. J. Appl. Acoust. 2018, 142, 29–35. [Google Scholar] [CrossRef]
  5. Bilal, M.; Liu, S.; Qiao, G.; Wan, L.; Tao, Y. Bionic Morse coding mimicking humpback whale song for covert underwater communication. J. Appl. Sci. 2019, 10, 186. [Google Scholar] [CrossRef] [Green Version]
  6. Qiao, G.; Ma, T.; Liu, S.; Bilal, M. A frequency hopping pattern inspired bionic underwater acoustic communication. J. Phys. Commun. 2021, 46, 101288. [Google Scholar] [CrossRef]
  7. Ma, T.; Liu, S.; Qiao, G.; Pu, W. Bionic underwater acoustic communication by mimicking dolphin whistle based on frequency shift keying. J. Electron. Inf. Tech. 2022, 44, 720. [Google Scholar]
  8. Ahn, J.; Lee, H.; Kim, Y.; Kim, W.; Chung, J. Multipath combining method for frequency shift keying underwater communications mimicking dolphin whistle. J. Acoust. Soc. Korea 2018, 37, 404–411. [Google Scholar]
  9. Lee, H.; Ahn, J.; Kim, Y.; Lee, S.; Chung, J. A biomimetic communication method based on time shift using dolphin whistle. J. Acoust. Soc. Korea 2019, 38, 580–586. [Google Scholar]
  10. Lee, H.; Ahn, J.; Kim, Y.; Lee, S.; Chung, J. Time-frequency modulation based mimicking dolphin whistle for covert underwater acoustic communication. J. Jpn. Appl. Phys. 2020, 59, SKKF03. [Google Scholar] [CrossRef]
  11. Ahn, J.; Lee, H.; Kim, Y.; Kim, W.; Chung, J. Machine learning based biomimetic underwater covert acoustic communication method using dolphin whistle contours. J. Sens. 2020, 20, 6166. [Google Scholar] [CrossRef]
  12. Lee, H.; Kim, Y.; Seol, S.; Kim, W.; Ahn, J.; Chung, J. Mimicking Multiple Whale Whistles-Based Underwater Covert Communication. J. IEEE Access 2022, 10, 43934–43942. [Google Scholar] [CrossRef]
  13. Ahn, J.; Lee, H.; Kim, Y.; Lee, S.; Chung, J. Mimicking dolphin whistles with continuously varying carrier frequency modulation for covert underwater acoustic communication. J. Jpn. Appl. Phys. 2019, 58, SGGF05. [Google Scholar] [CrossRef]
  14. Ahn, J.; Lee, H.; Kim, Y.; Kim, W.; Chung, J. Biomimicking covert communication by time-frequency shift modulation for increasing mimicking and BER performances. J. Sens. 2021, 6, 2184. [Google Scholar] [CrossRef] [PubMed]
  15. Yang, T.C.; Wenbin, Y. Low probability of detection underwater acoustic communications using direct-sequence spread spectrum. J. Acoust. Soc. Am. 2008, 124, 3632–3647. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Jun, J.; Hao, H.; William, R. Covert underwater acoustic communications. J. Acoust. Soc. Am. 2010, 148, 2898–2909. [Google Scholar]
  17. Xiujun, S.; Jun, W.; Haibin, W.; Xiaoxia, Y. Chaotic direct sequence spread spectrum for secure underwater acoustic communication. J. Appl. Sci. 2016, 104, 57–66. [Google Scholar]
  18. Fengzhong, Q.; Xiangzhao, Q.; Liuqing, Y.; Yang, T.C. Spread-spectrum method using multiple sequences for underwater acoustic communications. J. IEEE Ocean Eng. 2018, 43, 1215–1226. [Google Scholar]
  19. Freitag, L.; Stojanovic, M. Basin-scale acoustic communication: A feasibility study using tomography m-sequences. In Proceedings of the MTS/IEEE Oceans, Honolulu, HI, USA, 5–8 November 2011. [Google Scholar]
  20. Song, H.; Cho, S.; Kang, T.; Hodgkiss, W.; Preston, J. Long range acoustic communication in deep water using a towed array. J. Acoust. Soc. Am. 2011, 129, 71–75. [Google Scholar] [CrossRef]
  21. Kang, T.; Song, H.; Hodgkiss, W. Long-range multicarrier acoustic communication in deep water using a towed horizontal array. J. Acoust. Soc. Am. 2012, 131, 4665–4671. [Google Scholar] [CrossRef] [Green Version]
  22. Shimura, T.; Ochi, H.; Watanabe, Y.; Hattori, T. Results of basic at-sea experiments on time-reversal communication in the deep ocean. In Proceedings of the MTS/IEEE Oceans, Kobe, Japan, 8–11 April 2008. [Google Scholar]
  23. Shimura, T.; Ochi, H.; Watanabe, Y.; Hattori, T. Experiment results of time-reversal communication at the range of 300 km. J. Jpn. Appl. Phys. 2011, 49, 07HG11. [Google Scholar] [CrossRef]
  24. Shimura, T.; Watanabe, Y.; Ochi, H.; Song, H. Long-range time reversal communication in deep water Experimental results. J. Acoust. Soc. Am. 2012, 132, 49–53. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Shimura, T.; Ochi, H.; Watanabe, Y.; Hitori, T. Time reversal communication in deep ocean results of recent experiments. In Proceedings of the 2011 IEEE Workshop on Scientific Use of Submarine Cables and Related Technologies, Tokyo, Japan, 5–8 April 2011. [Google Scholar]
  26. Moore, B. An Introduction to the Psychology of Hearing, 6th ed.; Emerald Group Publishing: Binglry, UK, 2012. [Google Scholar]
  27. Munger, L.; Wiggins, S.; Hidebrand, J. North Pacific right whale up-call source levels and propagation range on the southeastern Bering Sea Shelf. J. Acoust. Soc. Am. 2011, 129, 4047–4054. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Crance, J.; Berchok, C.; Wright, D.; Brewer, A.; Woodrich, D. Song production by the North Pacific right whale, Eubalaena japonica. J. Acoust. Soc. Am. 2011, 145, 3467. [Google Scholar] [CrossRef] [PubMed]
  29. Parks, S.; Tyack, P. Sound production by North Atlantic right whales (Eubalaena glacialis) in surface active groups. J. Acoust. Soc. Am. 2005, 117, 3297. [Google Scholar] [CrossRef]
  30. Parks, S.; Tyack, P. Short- and long-term changes in right whale calling behavior: The potential effects of noise on acoustic communication. J. Acoust. Soc. Am. 2007, 6, 3725. [Google Scholar] [CrossRef]
  31. Lafferty, J.; Zhai, C. Document language models, query models, and risk minimization for information retrieval. In Proceedings of the 24th Annual International ACM SIGIR Conference on Research and Development in Information Retrieval, New Orleans, LA, USA, 9–13 September 2001. [Google Scholar]
  32. Geyer, J. Practical markov chain monte carlo. J. Stat. Sci. 1992, 7, 473–483. [Google Scholar] [CrossRef]
  33. Stephen, P. Markov chain Monte Carlo method and its application. J. R. Stat. Soc. 1998, 47, 69–100. [Google Scholar]
  34. Popov, V.; Langrock, R.; Deruiter, S.; Visser, F. An analysis of pilot whale vocalization activity using hidden Markov models. J. Acoust. Soc. Am. 2017, 141, 159. [Google Scholar] [CrossRef] [Green Version]
  35. Tennessen, J.; Holt, M.; Ward, E.; Hanson, M.; Emmons, C.; Giles, D.; Hogan, J. Hidden Markov models reveal temporal patterns and sex differences in killer whale behavior. J. Sci. Rep. 2019, 9, 14951. [Google Scholar] [CrossRef] [Green Version]
  36. Jiang, J.; Wang, X.; Duan, F.; Liu, W.; Bu, L.; Li, F.; Li, C.; Sun, Z.; Ma, S.; Deng, C. Study of the relationship between pilot whale (Globicephala mela) behavior and the ambiguity function of its sounds. J. Appl. Acoust. 2019, 146, 31–37. [Google Scholar] [CrossRef]
  37. Pace, F. Automated Classification of Humpback Whale (Megaptera novaeangliae) Songs Using Hidden Markov Models. Ph.D. Thesis, University of Southampton, Southampton, UK, May 2013. [Google Scholar]
  38. Available online: https://www.kaggle.com/competitions/whale-detection-challenge/data (accessed on 25 August 2022).
  39. Esfahanian, M.; Erdol, N.; Gerstein, E.; Zhuang, H. Two-stage detection of north Atlantic right whale upcalls using local binary patterns and machine learning algorithms. J. Appl. Acoust. 2017, 120, 158–166. [Google Scholar] [CrossRef]
  40. Available online: https://github.com/nmkridler/moby (accessed on 25 August 2022).
Figure 1. Spectrogram of whistles of the RHW: (a), (d) Pattern 1; (b), (e) Pattern 2; (c), (f) Pattern 3.
Figure 1. Spectrogram of whistles of the RHW: (a), (d) Pattern 1; (b), (e) Pattern 2; (c), (f) Pattern 3.
Sensors 22 08011 g001
Figure 2. Example of the STFT results of whistles with various Δ T and Δ F : (a) STFT with Δ T and (b) STFT with Δ T ^ .
Figure 2. Example of the STFT results of whistles with various Δ T and Δ F : (a) STFT with Δ T and (b) STFT with Δ T ^ .
Sensors 22 08011 g002
Figure 3. Modeling results of the time-dependent frequency change of the RHW whistle: (a) 0.4 s ≤ T w ≤ 0.5 s, Δ T : 0.15 s, Δ F : 13.3 Hz and (b) 0.3 s ≤ T w ≤ 0.5 s, Δ T : 0.2 s, Δ F : 10 Hz.
Figure 3. Modeling results of the time-dependent frequency change of the RHW whistle: (a) 0.4 s ≤ T w ≤ 0.5 s, Δ T : 0.15 s, Δ F : 13.3 Hz and (b) 0.3 s ≤ T w ≤ 0.5 s, Δ T : 0.2 s, Δ F : 10 Hz.
Sensors 22 08011 g003
Figure 4. Example of a whistle-mimicking modulation signal.
Figure 4. Example of a whistle-mimicking modulation signal.
Sensors 22 08011 g004
Figure 5. Block diagram of the proposed biomimicking modulation.
Figure 5. Block diagram of the proposed biomimicking modulation.
Sensors 22 08011 g005
Figure 6. The frame structure of the proposed biomimicking signal.
Figure 6. The frame structure of the proposed biomimicking signal.
Sensors 22 08011 g006
Figure 7. Block diagram of the proposed biomimicking demodulation.
Figure 7. Block diagram of the proposed biomimicking demodulation.
Sensors 22 08011 g007
Figure 8. Mimicking whistles for the DoM evaluation: (a) the proposed method, (b) FSK, and (c) CVC-FM and TFSK.
Figure 8. Mimicking whistles for the DoM evaluation: (a) the proposed method, (b) FSK, and (c) CVC-FM and TFSK.
Sensors 22 08011 g008
Figure 9. Simulation results of BER performance by delay profile: (a) 30 ms channel, (b) 1 s channel, and (c) 2 s channel.
Figure 9. Simulation results of BER performance by delay profile: (a) 30 ms channel, (b) 1 s channel, and (c) 2 s channel.
Sensors 22 08011 g009
Figure 10. Ocean experiment: (a) location, (b) configuration.
Figure 10. Ocean experiment: (a) location, (b) configuration.
Sensors 22 08011 g010
Figure 11. Multipath delay profile of the ocean experiment: (a) the proposed method and (b) CV-CFM and TFSK.
Figure 11. Multipath delay profile of the ocean experiment: (a) the proposed method and (b) CV-CFM and TFSK.
Sensors 22 08011 g011
Table 1. The result of the DoM evaluation.
Table 1. The result of the DoM evaluation.
Modulation SchemeProposedFSKCV-CFM and TFSK
Recognition probability95%10%60%
Table 2. Modulation parameters and data rate.
Table 2. Modulation parameters and data rate.
Proposed CV-CFMFSKTFSK
Δ T (s)
symbol duration
0.0750.10.1250.150.1750.20.0670.06670.067
M = 4
Δ F (Hz)
bandwidth
26.6201613.411.51015
M = 2
15
M = 2
30
M = 4
Data rate12 bps12 bps15 bps7 bps
Table 3. Summary of long-range acoustic communication experiment and channel.
Table 3. Summary of long-range acoustic communication experiment and channel.
Ref.BandwidthRangeNumber of PathsMaximum Delay
[16]Fc: 250 Hz
BW: 100 Hz
550 km230 ms
[17]Fc: 250 Hz
BW: 50 Hz
500 km
700 km
1~2600 ms
[18]Fc: 500 Hz
BW: 100 Hz
100 km3~41 s
[19]Fc: 400 Hz
BW: 100 Hz
300 km51 s
[20]Fc: 500 Hz
BW: 100 Hz
600 km15~172 s
[21]Fc: 500 Hz
BW: 100 Hz
1000 km5~81.5~2.5 s
Table 4. BER results of the ocean experiment.
Table 4. BER results of the ocean experiment.
Mod. Scheme and SensorProposedProposed w.o. Ch.-Comb.CV-CFMTFSK
1st Sensor0.00940.01130.09060.1394
2nd Sensor0.00300.00570.08300.1210
3rd Sensor0.00450.00990.06350.0890
Avg.0.00560.00900.07900.1165
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ahn, J.; Do, D.; Kim, W. The Long-Range Biomimetic Covert Communication Method Mimicking Large Whale. Sensors 2022, 22, 8011. https://0-doi-org.brum.beds.ac.uk/10.3390/s22208011

AMA Style

Ahn J, Do D, Kim W. The Long-Range Biomimetic Covert Communication Method Mimicking Large Whale. Sensors. 2022; 22(20):8011. https://0-doi-org.brum.beds.ac.uk/10.3390/s22208011

Chicago/Turabian Style

Ahn, Jongmin, Deawon Do, and Wanjin Kim. 2022. "The Long-Range Biomimetic Covert Communication Method Mimicking Large Whale" Sensors 22, no. 20: 8011. https://0-doi-org.brum.beds.ac.uk/10.3390/s22208011

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop