Underwater Target Localization Using Opportunistic Ship Noise Recorded on a Compact Hydrophone Array

: In this research, a new application using broadband ship noise as a source-of-opportunity to estimate the scattering ﬁeld from the underwater targets is reported. For this purpose, a ﬁeld trial was conducted in collaboration with JASCO Applied Sciences at Duncan’s Cove, Canada in September 2020. A hydrophone array was deployed in the outbound shipping lane at a depth of approximately 71 m to collect broadband noise data from different ship types and effectively localize the underwater targets. In this experiment, a target was installed at a distance (93 m) from the hydrophone array at a depth of 25 m. In this study, a matched ﬁeld processing (MFP) algorithm is utilized for localization. Different propagation models are presented using Green’s function to generate the replica signal; this includes normal modes in a shallow water waveguide, the Lloyd-mirror pattern for deep water, as well as the image model. We use the MFP algorithm with different types of underwater environment models and a proposed estimator to ﬁnd the best match between the received signal and the replica signal. Finally, by applying the scatter function on the proposed multi-channel cross correlation coefﬁcient time-frequency localization algorithm, the location of target is detected.


Introduction
Nowadays, ocean acoustic scientists pay significant attention towards to the use of passive approaches instead of active techniques for ocean sensing due to the environmental regulations prohibiting the transmission of invasive noise source; passive sensing also reduces the cost and energy consumption for remote nodes. Furthermore, there has been a paradigm shift in ocean acoustics utilizing ambient noise as an opportunistic source. Opportunistic sources refer to acoustic energy existing in a maritime medium that has not been generated for the clear aim of ocean measurement. Wave noise and noise generated by sea ice and ships are known examples of acoustic sources of opportunity [1]. Additionally, ship noise can be considered as an active acoustic source for monitoring the ocean, despite its undesirable effects on the life of marine mammals. To date, a few studies have applied ship noise as an opportunistic active source to estimate the geoacoustic properties of the seabed using commercial or small ships [2][3][4][5]. Ships produce broadband noise at frequencies ranging from 20 Hz to 100 kHz, due to the hydrodynamics, on-board machinery, and more importantly, from propeller cavitation. The rotation of the blades around the propeller axis causes the cavitation phenomenon [6,7]. Specifically, the cavitation introduces a rapid change of pressure in a liquid, which leads to small vaporfilled cavity formations in places where the pressure is relatively low. At higher pressures, these cavities, called "bubbles" or "voids", collapse and can generate a high intensity shock wave. Cavitation depends on several factors, including the vessel size, its speed and load, the ocean conditions, the age of the vessel, the engine type, and the propeller size. As such, ships are powerful and underutilized acoustic sources of opportunity that can be used to obtain information about the marine environment and target localization [8].
The propagation of underwater acoustic is mathematically described by the wave equation and must take into consideration the appropriate boundary conditions which are unique to different environmental conditions. There are five methods to solve the propagation wave equation: the fast field program (FFP), normal modes, ray theory, the parabolic equation (PE), and finite element (FE) solution [9,10]. The choice of method depends on the application, the deployment geometry, and the frequency of operation. For example, matched field processing (MFP) when considering ray-theory as a replica signal creates misleading results at low frequencies [11].
A variety of techniques have been studied to localized the objects in the ocean. In [12], common approaches are reviewed to localize the sources with wireless sensor networks that utilize the intensity at the receiver, the time of arrival (ToA), the time difference of arrival (TDoA) [13], or the direction of arrival (DoA) which are obtained by combining the multiple sensor signals. However, other algorithms have been introduced recently. For example, in [14], the distributed sub-band valley fusion (DSVF) method is proposed for a multi-array system. The target position is estimated using the intersection of the bearing lines of the receiver arrays with respect to the target. For each array, the bearing line is extracted from the two-dimensional (2-D) x − y steered response power (SRP) spectrum of the received signals of the array by implementing a 2-D sub-band valley energy detection (SVED) process. Additionally, in [15], a cluster-based direct source localization approach for large-aperture horizontal line arrays is proposed. Nonetheless, to localize the sources in complex media, MFP has been recognized as the most accurate technique and obtains an estimate of the source by matching the receive data at receivers and replica signals. In [16], an MFP method with measured replica fields is described, which provides instantaneous horizontal range estimates of a surface vessel in shallow water using a single hydrophone located above the sea bottom.
Popular techniques for localizing the underwater targets combine underwater wireless sensor networks (UWSNs) with active sonar systems [17]. Despite their effectiveness, these methods have many limitations. More importantly, the active source produces a negative impact on the hearing of marine mammals, particularly for those that use echolocalization to find their hunt. UWSN systems also have limitations, including energy consumption, channel capacity limitation, and low data rates. Additionally, the motion of the multiple sensors deployed sparsely reduces the accuracy of the target location estimation. The method that will be presented in this study is an opportunistic use of the ambient noise signals to identify the underwater targets; it does not suffer from these limitations, and the algorithm is completed using a single compact array equipped with a signal processing engine.
The focus of this research is to localize the target echo passively using ship noise. For this purpose, the MFP algorithm is utilized and combined with the scattering function according to a bistatic geometry in an isovelocity waveguide [18]. The isovelocity environment is justified for many deployment scenarios, and particularly in shallow water conditions. In these shallow conditions, the propagation can be represented using a waveguide, and the spatial complexities of the acoustic field are utilized in a MFP algorithm, which is a generalized beam-forming method to localize the active sources as a function of range, depth. Specifically, MFP is a spatial filter that compares the received signal and a replica signal modeled using Green's function in a waveguide [19]. It can be implemented using conventional or adaptive formulations for both narrowband and wideband signal models. Several of Green's formulations to construct the replica signal can be utilized, and the most common are normal modes in a shallow water waveguide, the Lloyd-mirror pattern, and the image model [11]. So, for localization purposes, using an accurate propagation model with an applicable localization method as an estimator is very important. The main purpose of this paper is to measure the return echo of the ship noise signal due to the reflection of the target surface, by utilizing a simple normal mode acoustic model to represent underwater propagation in shallow environments.
The main contributions of this work are two-fold. First, the opportunistic use of ship noise is described for the localization of a target; this technique avoids the generation of additional man-made acoustic signals that can threaten marine life. Second, a novel estimator using MFP is developed in this paper. The proposed estimator is based on the multi-channel cross-correlation coefficient (MCCC) localization algorithm. In the proposed algorithm, a time-domain cross-correlation matrix that contains the information of a specific frequency band is calculated, and the algorithm is optimized to reduce the number of mathematical operations on a real-time processor. The algorithm also takes into consideration the features of the propagation model. Furthermore, the algorithm is validated in a high-traffic shipping lane on the littoral of North America.
The structure of this paper is as follows: Section 2 describes theoretical background, such as the underwater acoustic propagation models and the scattered field from a target; Section 3 presents the methods and material including the proposed time-frequency MCCC algorithm to localize the underwater active source and experimental setup and dataset; Additionally, Section 4 demonstrates the performance of the proposed algorithm; finally, Section 5 highlights the conclusions of this work.

Theoretical Background
In this section, different approaches to model the acoustic propagation with different levels of complexity and scattering field from a target in a waveguide will be described. Specifically, in Section 2.1, a review of the Lloyd-mirror pattern will be provided; then, in Sections 2.2 and 2.3 will be described the image model, and the normal mode model, respectively; finally, in Section 2.4, will be explained the scattering field from the targets in a waveguide.

The Lloyd-Mirror Pattern (LMP)
In this section, the Lloyd-mirror pattern presented in [9] is reviewed. Let us assume that S(ω) is the value of a given acoustic source at an angular frequency ω and located at p 0 (x 0 , y 0 , z 0 ) on a Cartesian coordinate system. For a hydrophone array with N elements, the received signal H n (ω) at the n th element located at p n (x n , y n , z n ) has a relationship with S(ω) depending on the propagation model.
A LMP occurs when the acoustic source is located near the ocean surface and the ocean surface acts as a perfect reflector. Using LMP, the interference pattern in the acoustic field results from the constructive and destructive interference between the direct and source-reflected acoustic waves. Thus, the multipath is modeled using a direct and a reflected path from the ocean surface. As can be seen from Figure 1, there are two paths with lengths R 1 and R 2 in this model. Given a transmit signal with intensity S(ω) at frequency ω, the received signal intensity is expressed as [9] H n (ω) = S(ω) · P(p 0 , p n , where f and k are frequency and wave number, respectively. Additionally, the direct path length is R 1 = (x 0 − x n ) 2 + (y 0 − y n ) 2 + (z 0 − z n ) 2 , while the length of the path with a surface bounce is R 2 = (x 0 − x n ) 2 + (y 0 − y n ) 2 + (z 0 + z n ) 2 .

The Image Model
The image model is a generalized case of the LMP and takes into consideration multiple surface and bottom reflections. The image model contains multiple modes, and in this work each mode contains four paths. Figure 2 shows the paths for the mode m = 0 in the image model. Using the image model, the relationship between the source signal S(ω) and the received signal H n (ω) at the n th hydrophone at angular frequency ω is defined as where m and D are the mode number and depth of ocean, respectively. Additionally, for the mth mode, the range for path j = 1, . . . , 4 between the source and hydrophone n th is The depths z mj are for the m th mode and for paths j = 1, . . . , 4.
note that the location of the receiver and source are determined by z and z s , respectively.

The Normal-Mode Model in an Isovelocity Environment
The normal-modes model is a solution of the wave equation, defined as [9].
where p, ρ, , and ∂ are the acoustic pressure, the density of environment, the divergence operator, and the derivative operator, respectively. The following boundary conditions are assumed in this work:
A perfectly rigid bottom with lossless reflection at the sea bottom, i.e., when z = D.
Note that, on a first order approximation, this assumption is valid for the test conditions in this work since the bottom is rocky [9].
The general solution of the wave equation for solving the normal-mode problem is equal to where m, ρ(z 0 ), Ψ m and k rm are the mode number, the density of water at the depth of source, the mode function or eigen-function, and the horizontal wave number or eigenvalues, respectively.
In an isovelocity problem, the sound speed profile and density are constant. As such, the mode function is obtained according to where k zm , r and D are the vertical wave number, the range between the array and source, and the depth of sea, respectively. The maximum number of propagation modes can be bound to a finite value [9], such that only modes which cause k rm to be a real number are considered. Note that the m th mode has m zeroes. Thus, the maximum number of modes for the frequency f is determined as Table 1 shows the conditions for the proposed deployment scenario in Duncan's Cove, on the littoral of Nova Scotia, Canada. In this case study, there are 140 propagating modes. We also obtained the transmission loss for this environment which is defined by Using (7), Figure 3a shows the transmission loss as a function of range for a receiver when the depth of propeller, and the receiver depth are 1.3 m and 70 m, respectively, using 140 modes. As can be observed, since the deployment scenario can be considered a waveguide, the modes form a multipath channel, which causes the transmission loss to vary significantly as a function of distance.  To assess the feasibility of detecting a target that is deployed at different depths, the transmission loss between a potential target is also assessed as a function of depth and range. The simulations are done using the Legendre-Galerkin algorithm [20] to simulate the normal mode. Figure 3b shows the signal attenuation. As can be observed, the blue color indicates the locations with the lowest transmission loss, and as will be demonstrated in Section 4, the target is more likely to be identified if located the blue area, near the ship.

Scattering Field from Targets in a Waveguide
In this section, Ingenito's solution [21,22] for the scattering problem from a spherical target based on the mode coupling method is reviewed. Additionally, using the planewave scattering function with specific boundary conditions in an isovelocity environment, the scattered field is obtained. The boundary conditions take into consideration the rigid or pressure-release boundary conditions, the continuity of pressure, and continuity for the normal particle velocity conditions. The geometry of the scattering problem is shown in Figure 4, where the surface and bottom of this isovelocity waveguide are at z = −D, and z = H − D (H is the depth of water), respectively. The origin of this system is assumed at the center of the target. The source is located at (x 0 , y 0 , z 0 ) and the received field is assumed at (x, y, z) while a point of the surface target at (x , y , z ) is encountered with the input field. The channel frequency response for the received scattering field using the mode-coupling method is modeled by where P st (p 0 , p n , ω), P tr (p 0 , p n , ω), and P sr (p 0 , p n , ω) are the acoustic channel response between the target and the source, the acoustic channel response between the target and receiver, and the acoustic channel response between the source and receiver, respectively. Additionally, the points p 0 = (r 0 , θ 0 , φ 0 ), and p t = (r t , θ t , φ t ) are the source point, and the target point, respectively in the spherical coordinates system. According to the Ingenito's solution and the near-horizontal propagation approximation, B n (r, r n , φ 0 , φ, z t ) = P n st (p 0 , p n , ω)P n tr (p 0 , p n , ω) is equal to where S(θ n = 0, φ n ; θ 0 = 0, φ 0 ), z t , r, r n are the far-field scattering function, the depth of the target, the range between the source and target, and the range between the target and receiver, respectively. The far-field scattering function for a soft radius a is expressed as while for the rigid sphere, the far-field scattering function is expressed as where cos(β) = sin(θ) sin(θ 0 ) + cos(θ) cos(θ 0 ) cos(φ − φ 0 ), and J n and H n are the spherical Bessel function, and Hankel function, respectively.
Using the normal mode model, the acoustic channel gain as a function of frequency between the source and receiver P sr (p 0 , p n , ω) is calculated using (4). Figure 5 shows the far-field scattering function for both the rigid and soft target with a radius of 1 m.
The frequency of source, and the sound speed are 2.5 kHz, and 1490 m/s, respectively. Additionally, the source is located at a coordinate of (θ 0 = 315.01, φ 0 = 0). As can be observed in Figure 5a, the scattering field of the soft target is smoother than for the rigid target. Additionally, the maximum field is aligned with the direction of the input source for the rigid and soft target.  A simulation for a target with a radius of 10 m is run to evaluate the effect of the target size on the scattering field. As shown in Figure 5b, it is obvious that the scattering field intensity is improved by about 10 dB by increasing the radius of the target by a factor of 10. Additionally, the directivity of the scattered field is increased in comparison to a smaller target. The effect of the frequency on the scattering field is also validated by simulation. Figure 5c shows the scattered field for a low frequency source at 1 Hz. According to the simulation results, the low-frequency source generates a stronger scattering field than the high-frequency source, and the field is broadcast in all directions uniformly.

Methods and Material
The proposed estimator for MFP and experimental dataset have been presented in Sections 3.1 and 3.2, respectively.

Proposed Localization Algorithm for MFP
The proposed algorithm is based on the multi-channel cross correlation coefficient (MCCC) localization algorithm [23]. MCCC is a time-domain approach to estimate the time delay of arrival (TDOA) in environment when the hydrophone (or microphone) array has more than two elements. Assume that s(t k ) is the time value of a given narrowband acoustic source at time t k , and is located at the point p 0 = (x 0 , y 0 , z 0 ) in the Cartesian coordinate system. This value is received at the n th hydrophone with the location p n = (x n , y n , z n ) at the time t k + τ n . Thus, for a hydrophone array with N elements in a noiseless free space, the received signal from single path at the n th element at time t is calculated as where α(p 0 , p n ) is the loss, α(p 0 , p n ) ∈ {R|0 ≤ α(p 0 , p n ) ≤ 1}. Additionally, where c is the sound speed. One can estimate the transmit signal s n (t k , p 0 ) from hydrophone n as The measured signals are modeled by r n (t k ), n = 1, 2, . . . , N at discrete time t k = t min + k ∆t, where ∆t is the sampling time k = 0, 1, . . . , L − 1, and s n (t k ,p 0 ) is the discrete-time transmit signal at positionp 0 . If the pointp 0 = (x 0 ,ỹ 0 ,z 0 ) is the correct location of the source then all N values s n (t k ,P 0 ), n = 1, 2, . . . , N are almost completely correlated with each other at the time t k . Considering the different propagation models defined, one can represent the transmit signal in the frequency domain as Note that R n (ω) is the l th frequency component of the FFT of the measured time signal r n (t k ), and P(p 0 , p n , ω) is the channel impulse response that is obtained according to the underwater acoustic propagation model (for example using the image model or the Lloyd-mirror pattern) between the source position p 0 , and the sensor at position p n . The discretized frequency samples are equal to ω l = l ∆ω, l = 0, 1, . . . , where . denotes the floor value of its argument, also, L is the length of signal. Because the signal is a real value, only half of the signal is considered.
is the FFT ofs(t k , p), Note that | . | denotes the absolute value of its argument. From [9], the IFFT of is defined asS(ω, p) is defined as where Re{.} denotes the real part of its argument. Additionally, 0 = 1, and l = 2 for l > 0. Then, one can represent (16) as where Im{.} denotes the imaginary part of its argument, and By calculating (17) in a specific frequency range for specific values of l ∈ {0, 1, . . . , (L/2) − 1}, that are corresponding to prefixed frequencies. The frequency limited time-domain crosscorrelation matrix can be simplified to where C c (p) = 1 L ∑ L−1 k=0 s (t k , p)s(t k , p) , is the cross-correlation matrix. Finally, the location of the source is estimated using [9] p 0 = arg min In the proposed algorithm, a time-domain cross-correlation matrix which contains the information of specific frequency band is calculated using an approximation on the inverse fast Fourier transform (IFFT). The proposed method, in addition to taking advantages of both time and frequency-domain methods, also considers the features of the propagation model of the environment. This is a very important feature particularly in a shallow water environment, in which the received signal arrives at the receiver through multiple paths. Indeed, using a propagation model that considers other paths in addition to the direct path improves the underwater localization accuracy. The proposed method can be used with different propagation models, such as the Lloyd-mirror, image, and normal-modes.

Experimental Setup and Data Collection
The acoustic data analyzed in this paper were produced by ships and were collected over 21 days. The experiment was conducted in collaboration with JASCO Applied Sciences in September 2020 in Duncan's Cove, Nova Scotia, Canada. A map showing the location of the hydrophone array and target, as well as the boundary for the shipping lanes are shown in Figure 6. A 4-element hydrophone array was placed at approximately 900 m away from the outbound traffic lanes, and the target was at a range of 93 m from the hydrophone array, between the hydrophone and the shipping lanes. The exact location of target and hydrophone array is shown in Table 2. The geometry of the hydrophone array was designed, such that three hydrophones formed a planar array, while the fourth hydrophone was placed above one of the hydrophones on the planar array, such that they formed together a 2-element vertical array. During the trials, the depth of water was 71.6 m but the tidal conditions caused changes in the water depth on different days as can be observed in Figure 7 where the variation of depth over time is offset from the location of the depth logger. Note that on 23 September 2020, Hurricane Teddy affected the depth log and acoustic data, which was discarded during this period. Because the depth is not constant over time, it should be noted that the propagation environment is not entirely static. The target is made of four 17-inch glass spheres from Vitrovex. They are mounted as pairs in two SUBS C-2 streamlined floats, that were joined together and installed at a depth 25 m. This experiment was run for 21 days from 8 September to 29 September 2020. Several vessels passed during this period and their noise was recorded by the hydrophone array.
Automatic Identification System (AIS) data was purchased from Marine Traffic (https://www.marinetraffic.com/) within the limited experimental area to obtain the exact position of each passing ship over the experimental area for the proposed experiment.

Simulation Using the Synthetic Data
Prior to the deployment, different scenarios were modeled to assess the localization algorithm in realistic environment conditions. The purpose of doing the pre-tests is to evaluate the ambiguity function in different areas if the target is in those areas. Figure 8 presents the simulated ambiguity functions for the different scenarios. In these simulations, the parameters, such as the position of the source (the ship propeller), the position of the receiver, the size of the target, the frequency of the emitted signal from source, and the number of modes, are assumed to be known. Then, the most likely situations in which the target can be detected are identified by simulating (9) for each point in space. The simulation results for six different scenarios demonstrate the best position for the target location, i.e., where the ambiguity surface has the most energy. It can be observed that the target should be located on the transmitter-receiver line to be most likely detected. Additionally, other parameters such as the target radius and depth of the target in water can affect the maximum energy of the ambiguity function. Specifically, as the target size increases, and as it is deployed closer to the surface, the far-field scattered acoustic pressure increases.
The main purpose of this paper is to localize an underwater target due to propeller noise scattered from the target surface. As such, the propeller ship noise acts as the active sonar probing sequence. To estimate the position of the target, the received scattering field at the hydrophone array is processed using the algorithms described in Section 2.4. The target localization procedure (by MATLAB) is summarized as shown in Figure 9. First, the received signals at the hydrophones are converted to the frequency domain using a discrete-time Fourier transform. Then, the replica signal and the suggested underwater channel defined by (8) are modeled using different models of the acoustic propagation media. Finally, by sweeping the space, the position of the target is determined. ((0,0,70m)  In this paper, we use the ship noise that was recorded on September 8 and 10 to evaluate the accuracy of the localization algorithm. First, JASCO's vessel present during the deployment on September 8 is used as a source of opportunity. The vessel is also tracked using the AIS data to provide a reference location.

Scattered value at the hydrophone
To localize the underwater target in a pseudo-active manner using existing ambient signals as a source, the acoustic signal produced by the vessel, combined with the scattering signal at the hydrophone array, is recorded. By post-processing the data, the received filtered signal is applied to the algorithm solving (8) to evaluate the estimator accuracy based on using different acoustic propagation models. The parameters used for the simulation are given in Table 3. According to the simulation results presented in Section 4.1, the best scenario for identifying the target position is to place the array, the target and the source (propeller) along a straight line. So, when the ship is located at (x = 442.7 m, y = 22.7 m, z = 1.2 m ), the target localization algorithm is run using the received signals at the recorder. Figure 10a shows the result simulating the proposed time-frequency MCCC estimator using the normal-mode model as the underwater acoustic propagation model. The white square in this figure shows the maximum value of energy on this ambiguity plot. This localization value for the target is estimated at a range of 93.6 m from the array and at a depth of 24.8 m below sea surface. This demonstrates the accuracy of the target localization algorithm. Additionally, Figure 10b indicates the accuracy of the target localization algorithm for the proposed time-frequency MCCC estimator using the image model with four paths. The maximum in this case is estimated at a distance of 91.4 m of the array. In comparison, the estimated value calculated using the Lloyd-mirror model produces an estimate of 116.1 m, as shown Figure 10c. The reason for this difference in the estimated values can be explained by the inaccuracy of Lloyd-mirror's model as an underwater acoustic propagation model, and specifically because the Lloyd-mirror pattern models multipath using only two paths. In comparison, the normal mode model considers the characteristics of the underwater environment such as water velocity (albeit assumed to be constant in this shallow environment) and density in addition to the multipath effect. Table 3. Parameters needed to localize the target using JASCO's ship noise.

Parameters Values
Target diameter 1.  Table 4, z represents the estimated depth of the target. Additionally, the horizontal range and the radial range (the radial range is the distance between the hydrophone array and the estimated target location) are, respectively, calculated as x 2 + y 2 and where z h is the depth of the hydrophone array (which is 70 m). Furthermore, the relative error is the absolute error divided by the magnitude of the true value.  Secondly, to confirm the validity of the obtained results, the algorithm is validated on 10 September 2020 using the CMA CGM BRAZIL shown in Figure 11. The CMA CGM BRAZIL is one of the world's biggest cargo ships which berthed in Halifax's harbor on 10 September 2020. This ship can carry 15,000 containers and is 366 m in length and 51 m wide. Since it passed in the predefined zone of the AIS data collected, it was possible to track its position. This ship was at approximately 1 km from the hydrophone array, as it was navigating towards the port of Halifax. We recorded its position to be at (x = 907 m, y = 130 m, z = 1.9 m), with respect to the array location. Using normal mode and the scattering model, the target location is estimated. In this test, the desired frequency band is between 500 Hz and 1500 Hz. Other simulation parameters are summarized in Table 3. As shown in Figure 12a, the maximum energy is obtained for a location of (x = 88 m, y = 4 m), so the estimated range of target from the hydrophone array is approximately 88 m.
The reason for the slight difference in the value of actual range from the estimated one can be due to the imprecision in the actual distance of the ship propeller from the array, as well as the limited accuracy that can be obtained with only two vertical elements in the array configuration.
Next, an attempt is done to estimate the location of the target using the image model and Lloyd-mirror pattern with the data collected. Figure 12b,c show the ambiguity function for the image model and Lloyd-mirror pattern, respectively. In this case, the maximum is estimated at a location near the hydrophone array instead of the target, but a yellow area (at a range of about 91 m) is seen at the location of the target on the ambiguity surface. This confirms the importance of an accurate propagation model to estimate the location of the target. Table 5 shows the estimated location for the target by considering two maximum values, the global or the first local maximum point (shown by the superscript (F)), and the second local maximum point (shown by the superscript (S)). The purpose of this Table is

Conclusions
In this paper, a new technique was presented based on the use of ambient noise as a type of secure signaling. The proposed technique prevents the negative impact of active sonar on the hearing of underwater animals, and as such it is an environmentally friendly solution. Additionally, using the propeller noise is more secure than active sonar, and consumes less energy. The signal generated by the rotation of the propellers of moving ships can act as a sonar signal source and as an opportunistic noise source that detects the position of the target by receiving its scattering from the surface of objects using the hydrophone arrays. Based on our knowledge, the present work is the first application using the scattering of ship noise towards underwater localization of targets. In this paper, the return echo from the target surface was processed using the proposed MCCC, which is a low complexity time-frequency method based on MFP. A comparison of the algorithm accuracy for different acoustic propagation models was done. The algorithm was validated in a shallow environment on the littoral of North America, and it was demonstrated that a target with radius of 1.6 m installed at a depth of 25 m could be detected using the propeller noise from the ship, as well as using the proposed estimator. The simulation results also showed that the proposed MCCC-time-frequency estimator and normal mode model are resistant to aliasing effect and ambient noise.
Although the proposed method does not require an active source, it still has its limitations, because the underwater channel is a non-stationary environment and has a random signature that varies with time. Therefore, a more realistic model of acoustic propagation will improve the localization accuracy. Additionally, an important challenge for the proposed method occurs when more than one ship is in the test area, in which case the noise generated by each ship has a detrimental effect on the signal generated by the propeller of the desired ship. The impact is most severe when the noise signal generated by other ships is in the same frequency band as the original signal, causing interference on the array, and causing positioning errors.
In future work, the proposed estimator will be simplified to a new equation for implementation on the field-programmable gate array (FPGA) board to localize the underwater active source.