Underwater Sound Characteristics of a Ship with Controllable Pitch Propeller

: The time-dependent spectral characteristics of underwater sound radiated by an ocean vessel has complex dependencies on ship machinery, propeller dynamics, hydrodynamics of ship exhaust and motion, as well as ship board activities. Here the underwater sound radiated by a ship equipped with a controllable pitch propeller (CPP) is analyzed and quantiﬁed via its (i) power spectral density for signal energetics, (ii) temporal coherence for machinery tonal sound, and (iii) spectral coherence for propeller amplitude-modulated cavitation noise. Frequency-modulated (FM) tonal signals are also characterized in terms of their frequency variations. These characteristics are compared for different propeller pitch ratios ranging from 20% to 82% at ﬁxed propeller revolutions per minute (RPM). The efﬁcacy and robustness of ship parameter estimation at different pitches are discussed. Finally, analysis of one special measurement is provided, when ship changes speed, propeller pitch and RPM over the duration of the measurement. The 50% pitch is found to be a crucial point for this ship about which tonal characteristics of its underwater radiated sound attain their peak values, while broadband sound and associated spectral coherences are at a minimum. The ﬁndings here elucidate the effects of pitch variation on underwater sound radiated by ships with controllable pitch propellers and has applications in ship design and underwater noise mitigation.

The mechanisms for the sound generated by ships are described in [9,29,[34][35][36][37][38][39]. The sources of ship-radiated underwater sound include (i) persistent onboard machinery, such as engines and air-conditioning systems, (ii) propulsion dynamics, such as propeller and shaft rotation, (iii) hydrodynamical effects, such as cavitation, vortex shedding, intakes and exhausts, and (iv) other time-limited or transient activities, such as cargo handling and mooring and slamming motions [4,9]. Machinery and propeller noise dominate a ship's radiated underwater sound [9]. The overall underwater sound emanated by a ship comprises colored broadband noise over a wide frequency range embedded with several prominent narrowband tonal components from onboard machinery, along with time-varying and amplitude-modulated sound pressure levels due to propeller rotation.
Here, we analyze and characterize the underwater sound radiated by a specific ship with a controllable pitch propeller (CPP) and received on three vertically deployed hydrophones. The ship-radiated underwater sound characterization follows the approach of Refs. [40,41] and includes the calculation of (i) the power spectral density (PSD) for broadband sound energetics, (ii) the temporal coherence for machinery tonal sound, and (iii) the spectral coherence for amplitude-modulated propeller cavitation noise. The frequency variations of frequency-modulated (FM) tonal signals that are probably due to the non-uniform wake from the propeller blades caused by the ship's physical motions, such as yaw, roll or pitch, are also quantified. The calculations are initially made on measurements with ship at fixed propeller pitch ratio and revolutions per minute (RPM). These features are then compared for different propeller pitch ratios ranging from 20% to 82% at a fixed RPM. The 50% propeller pitch ratio was found to be a crucial point, at which multiple features of ship noise attain either their peak or minimum values. The efficacy and robustness of the ship parameter estimation at different pitches are discussed. Finally, an analysis of one special measurement is provided when the ship's speed, propeller pitch, and RPM are varied over the duration of the measurement.
Studies of the effects of propeller pitch variation on the underwater sound radiated by ships with controllable pitch propellers are limited since most previous works analyzed data from ships with fixed-pitch propellers.. The underwater sound radiated by small ships is analyzed in [39,42,43], while that from large commercial ships, including tanker, carrier, cargo, and container ships, is compared and quantified in terms of source level and spectral shape in [44]. Approaches for quantifying the cyclostationarity of ship propeller noise are developed in [45,46].
Underwater ship-radiated sound analyses in the published literature are often based on single hydrophone measurements. Other sensing modalities that deploy more than a single hydrophone have been implemented. In [4], ship cavitation noise was analyzed using the Stevens Passive Acoustic Detection System (SPADES) with two crossed pairs of hydrophones deployed near a harbour. A large aperture, densely-sampled coherent hydrophone array with 160-elements was employed to analyze the underwater sound radiated from multiple ships using the Passive Ocean Acoustic Waveguide Remote Sensing (POAWRS) technique over instantaneous wide areas, with >100 km diameters, in several continental shelf environments [40,41,47,48]. The analysis was focused on the detection and characterization of ship tonal signals, as well as passive ship localization in [47,48]. Three different approaches were applied to the POAWRS dataset in [40,41] for the simultaneous and automatic detection, bearing-time estimation, and acoustic signature characterization of multiple ships over long ranges, which included broadband energetics, temporal coherence for tonal extraction, and spectral coherence for cyclostationary propeller cavitation noise analysis.
Models of ship-radiated underwater sound spectra for a wide range of ships are useful in ship design and enhancement, sonar performance prediction, and modeling the impact of ship noise on communication between marine fauna. In [49], ship noise spectral levels predicted using empirical models were compared to measurements for several ships. In [50,51], a parametric analysis of the noise spectra radiated by commercial ships operating in different conditions was developed. A new spectral parametrization model was proposed that delivers a better performance than traditional methods due to greater model flexibility, based on calibration with empirical data. These analyses included ships with controllable pitch propellers, one of which is analyzed here for signature characterization in terms of propeller pitch ratios.

Data Collection
The acoustic data were acquired in an experimental campaign carried out within the framework of EU FP7 collaborative projects AQUO [52]. The measurement procedure followed the ANSI/ASA standard [53]. The ship-radiated underwater sound was received on buoy-mounted hydrophones (see Figure 1a) [49]. The ship under test passed the hydrophone array in a straight course at different operative conditions and different lateral distances (see Figure 1b). The ship's course and speed were kept constant from the COMEX to FINEX points. The underwater acoustic data was recorded on the hydrophones during the time window from the "Start Data" to "End Data" points identified in Figure 1, spanning a 30 degree angle on both sides of the closest point of approach. For the present study, only data acquired by the top hydrophone was analyzed. The sampling frequency was 13.11 kHz. The ship used for the experimental campaign was the research vessel (RV) Nawigator XXI, equipped with a controllable pitch propeller (CPP). The basic physical properties of this vessel are listed in Table 1, while the ship's propeller  and engine parameters are provided in Tables 2 and 3, respectively. The ship's machinery rotation parameters were obtained, including cylinder firing rate CFR, engine firing rate EFR, propeller shaft frequency SF, and propeller blade pass frequency BPF, as shown in Table 4. These parameters were calculated as follows: where n represents the nth harmonic order number, RPM e is the engine RPM, RPM p is the propeller RPM, RGR is the reduction gear rate, N c is the number of cylinders, N s is the number of strokes, and N b is the number of blades on the propeller.

Ship-Radiated Underwater Sound Analysis Methods
Underwater sound radiated by a ship is a complex superimposition of both broadband noise, mainly due to cavitation and the flow around the hull, and tonal components with specific frequencies due to the periodicity of on-board machinery and the propulsion system. Consequently, different techniques can be used to highlight different characteristics of the sound emitted by a ship. Here, we employed the approaches of Refs. [40,41] to analyze the underwater acoustic hydrophone data containing ship-radiated sound from RV Nawigator XXI with a controllable pitch propeller, as described below.

Power Spectral Density for Signal Energetics
Traditional power spectral density (PSD) is calculated to reveal signal's power distribution across a broad frequency band. The spectral curve variations reveal the features of ship-radiated broadband noise. Energetic narrowband tonal signal components can also be observed in signal PSD, especially at low frequencies below 2 kHz. The Welch's averaged modified periodogram method [54] was used to estimate PSD: f n = n L n = 0, . . . , L/2 (7) where A k (n) represents the finite Fourier transforms of the windowed sequence x k (j)w(j), I k ( f n ) represents the corresponding modified periodogram, and U represents the factor of window weight. The spectral estimationP( f n ) was then obtained from the average of the periodograms.

Temporal Coherence Analysis for Tonal Signals
The persistent narrowband tonal sounds at distinct frequencies that are radiated from ship machinery are typically present over long time durations and hence, are highly temporally coherent. These tonal sounds can be detected and quantified using temporal coherence calculations [40]. As the tonal signals have energies that are confined to relatively narrow frequency bands and centered at certain frequencies over long time durations, the signals in adjacent time periods are highly correlated at those tonal frequencies. Here, the magnitude-squared coherence MSC [55,56] between two adjacent time-windowed signals was first calculated as a function of signal frequency: where x w (t) and x w+1 (t) represent signals in two adjacent time windows w and w + 1. Then, the cross-power spectral density P x w ,x w+1 and the power spectral densities P x w ,x w and P x w+1 ,x w+1 were estimated using Welch's averaged modified periodogram method, as shown in Equations (5)-(9). Next, the temporal coherence in one data unit could be obtained as the averaged magnitude-squared coherence MSC: where N is the total window number in one data unit. The tonal components were obtained from the peaks of temporal coherence plotted as a function of frequency TC( f ). Additionally, variable window sizes could be chosen to provide variable frequency resolutions within different frequency ranges.

Frequency-Modulated Tonal Signal Analysis
At high frequencies, tonal signals have been found to have large frequency oscillations around their center frequencies [48] because they have been frequency-modulated (FM). To analyze these FM tonal signals, tone tracks of these signals in spectrograms are extracted. Tone tracks describe the time variation of the fundamental frequency of the tonal signal [57,58]. To estimate the modulation frequency, DC offset is first removed from the tone tracks of FM tonal signals. Then, by applying a Fourier transform on the resulting tone track, the dominant modulation frequency is obtained from the peak.

Spectral Coherence Analysis for Cyclostationary Cavitation Noise
A ship's rotating propellers generate cavitation noise that is amplitude-modulated [9] or cyclostationary and has hidden periodicities or repetitive patterns, for instance, periodic energy flow [59,60]. This kind of signal is spectrally coherent, inherently arising from interfering spectral components spaced apart by the fundamental and harmonic frequencies of the propeller rotation [40]. Spectral coherence analysis provides an approach for determining the cyclic frequencies. It is based on spectral correlation, defined as the Fourier transform of the cyclic auto-correlation function [61]: where R α x (τ) is the cyclic auto-correlation function of signal x(t) at cyclic frequency α and S α x ( f ) is the corresponding spectral correlation. The fast spectral correlation estimator S Fast x (α, f ) [60] was utilized to estimate the spectral correlation, which has been proven to have a great statistical performance along with the advantage of rapid computation, where α represents cyclic frequency and f c represents carrier frequency: where S x (α, f k ; p) represents the "scanning spectral correlation", which estimates spectral correlation from the discrete Fourier transform of the interactions between the STFT coefficients in frequency bins f k and f p−k , not necessarily spaced apart by exact α but by p∆ f , which is the closet frequency bin to α with some residue δ, in order to smoothly envelope the wave packet. R w (α) represents the window kernel from window function w with a central time index at N 0 . K is the number of windows, N w is the window length, and P is the maximum value of p. Based on this estimator, the spectral coherence γ x (α, f c ) could be obtained: The mean spectral coherence γ x (α) in the cyclic frequency α domain was obtained by averaging across the carrier frequency f c : The periodicities in propeller cavitation noise energy flow could be estimated from the peak components of the mean spectral coherence, corresponding to rotation frequencies including harmonics of ship's propellers.

Ship Energetics via Power Spectral Density Analysis
The power spectral density calculations revealed the ship's cylinder firing rate CFR. The spectrogram analyses based on the short-time Fourier transform of the recorded data and the corresponding power spectral density are shown in Figure 2a-d, respectively, for the fixed propeller pitch ratio of 65%. Ship noise was dominant in the measured data from 0 Hz to 3 kHz and could be significant up to 30 kHz at and near the closest point of approach. The ship-radiated sound included narrowband machinery tonal signals and broadband propeller cavitation noise. The sequence of click signals around 50 kHz was due to the ship's echo-sounder (see Figure 2a). At and near the closest point of approach, Lloyd's mirror effects were present in the measurements (Figure 2b) and the tonal signals were most dominant.
A series of harmonic tonal signals were present in the PSD distributions (see Figure 2d). The fundamental frequency could be estimated from the spacings between adjacent harmonic components, estimated to be around 6.3 Hz, and corresponded to the cylinder firing rate CFR. These CFR harmonics could only be detected from the seventh harmonic onwards, which was around 44.1 Hz. The strongest CFR component was located at 56.7 Hz, corresponding to the ninth harmonic. This component was slightly larger than the theoretical first harmonic component of the engine firing rate EFR that occurs at the eighth harmonic of CFR. As a result, EFR and hence, the cylinder number as the ratio of EFR and CFR could not be accurately estimated.

Ship Tonal Signals via Temporal Coherence Analysis
The ship's tonal signal components were determined via temporal coherence analysis in Figure 3. Here, three window sizes were used to provide different frequency resolutions in distinct frequency ranges. At higher frequency ranges, a wider frequency resolution was used since the tonal signals there have larger frequency oscillations. The temporal coherence calculations detrended the data so that the CFR tonal frequency components could be revealed more readily than in the PSD distribution (compare Figure 3a to Figure 2d).

Cyclostationary Cavitation Noise via Spectral Coherence Analysis
The spectral coherence analysis of the ship-radiated sound enabled both the propeller blade pass frequency BPF and shaft frequency SF to be determined, the ratio of which was the blade number. The two-dimensional spectral coherence distribution as a function of cyclic frequency (α) and carrier frequency ( f c ) from the analysis of the recorded data is shown in Figure 4a. The dominant spectral coherence (SC) values were located between 10 kHz and 40 kHz in the carrier frequency f c corresponding to the ship's propeller cavitation noise and around 50 kHz in carrier frequency f c corresponding to echo sounder signals. The mean spectral coherence obtained by averaging the entire carrier frequency f c range is shown in Figure 4b. By selecting different averaging intervals in carrier frequency f c , the spectral coherence means for the echo-sounder signals and the ship's propeller noise could be separated, as shown in Figure 4c,d, respectively. The harmonic group shown in Figure 4c had a fundamental frequency of 2.885 Hz, which corresponded to the ship's echo-sounder transmission rate. The harmonic group shown in Figure 4d had a fundamental frequency of 3.384 Hz, which corresponded to the propeller shaft frequency SF. Furthermore, the propeller blade pass frequency BPF given by the most dominant component shown in Figure 4d was 13.54 Hz and corresponded to the fourth harmonic of the shaft frequency SF. Hence, the propeller blade number could be estimated to be 4. These estimates were in good or exact agreement with the values provided in Tables 2 and 4 for the known ship's parameters.

Power at Different Propeller Pitches
Here, we compare the ship-radiated underwater sound power spectral density (PSD) levels at seven different propeller pitches. For signals between 100 Hz and 2000 Hz, the PSD level was the strongest at a 20% propeller pitch ratio (Figure 5a). The PSD then decreased with the increase in pitch until 51% pitch, where it was at a minimum, after which the PSD level increased with the increase in pitch. The PSD level below 250 Hz showed more complex variations with pitch (see Figure 5b). For instance, the 20% pitch led to the lowest PSD levels below 75 Hz. The differences were due to the different kinds of cavitation taking place at the propeller. For a detailed analysis of the different noise emissions of each cavitation phenomenon, see [50].

Tonal Signal Variations at Different Propeller Pitches
The power spectral density (PSD) level and temporal coherence values TC are plotted for significant tonal signals, i.e., the 7th to 13th harmonics of CFR, as a function of propeller pitch ratios in Figure 6a,b. The ninth harmonic of CFR had the highest PSD level and TC value at every propeller pitch. The PSD level of all harmonic tonals of CFR varied only slightly, between 2 dB to 4 dB with changes in pitch, with the exception of eighth harmonic where its PSD level changed by as much as 10 dB with variations in pitch, but still remained below the PSD level of the ninth harmonic. The temporal coherence values TC for most of the CFR tonals were high or at the maximum at 51% propeller pitch. This could be explained by the fact that the cavitation noise spectra was shifted toward higher frequencies above 500 Hz for the 51% propeller pitch, as can be noted from Figure 5b. Consequently, the tonal components below 100 Hz emerged more clearly in the spectrum.   Frequency-modulated (FM) tonal signal features are also compared at different propeller pitches in this paper. Figure 7 shows a spectrogram containing two narrowband FM tonal signals with frequencies wavering in the ranges of 1435-1460 Hz and 1145-1170 Hz, respectively. The mean and standard deviation of the frequency of these two FM tonal signals are plotted as functions of propeller pitch ratios in Figure 8a,b. The mean frequencies of both FM tonal signals could vary by 10 Hz with changes in propeller pitch and were consistently low between 60% and 80% pitch. The standard deviations of FM tonal frequencies were also the largest at high propeller pitch. The modulation frequencies of these FM tonal signals were the largest near 50% pitch, implying that the signal tone frequencies fluctuated the most at this pitch (see Figure 8c).

Propeller Cavitation Noise at Different Pitches
Next, propeller cavitation noise is compared at different propeller pitches via spectral coherence. The mean spectral coherence γ(α) distribution of propeller cavitation noise plotted as a function of cyclic frequency α is shown in Figure 9a-g for different propeller pitches. Table 5 summarizes the capacity for the estimation of the ship's propeller parameters from spectral coherence distributions at different propeller pitches, where SF1-SF4 are the first to fourth harmonics of shaft frequency and the fourth harmonic of shaft frequency SF4 is the propeller blade pass frequency BPF.   The mean spectral coherence γ(α) distributions were ideal at 82% and 79% pitches, as shown in Figure 9a, since the propeller shaft frequency SF, blade pass frequency BPF, and blade number could be readily and correctly estimated. These propeller parameters could also be estimated at 65% and 31% pitches, but the mean spectral coherence values were smaller by a factor of four to six. At 65% pitch, the second harmonic of shaft frequency SF2 could not be identified (Figure 9c). At other pitch ratios, the situation was far from ideal. At 51% pitch, the second harmonic of shaft frequency SF2 was much larger than the fundamental SF1 and the blade pass frequency BPF, which was SF4, implying that the correct blade number could not be estimated (Figure 9d). At 40% pitch, the magnitudes of the first four shaft frequency harmonics decreased with the increase in harmonic order (Figure 9e), so both blade pass frequency BPF and blade number were uncertain. The worst case scenario occurred at 20% pitch (Figure 9g), since only the fundamental shaft frequency SF1 could be identified. This could be due to pressure side cavitation taking place for negative angles of attack. Such an operative condition was very far from the designed one and the broadband cavitation noise was very high, as can also be seen from the black curve in Figure 5b for 20% pitch. The dependence of the spectral coherence on the carrier frequency f c at shaft frequencies provides further insights into ship propeller cavitation noise [62]. For several representative pitch ratios, the spectral coherences of the first four shaft frequency harmonics SF1-SF4 are plotted as functions of carrier frequency f c in Figure 10a-d. At 82% pitch, both the first and fourth harmonics of shaft frequency, SF1 and SF4, had high spectral coherence values in the carrier frequency f c range from 10 kHz to 30 kHz. The second and third harmonics of shaft frequency, SF2 and SF3, had comparatively low coherence values in the same carrier frequency range but were still significant. This scenario ensured that the propeller parameters could be estimated correctly. At the other pitches shown in Figure 10b-d, the fourth harmonic of shaft frequency SF4 was consistently weak or similar in magnitude to the other harmonics, so propeller blade pass frequency and blade number could be harder to infer.
The spectral coherence γ(α, f c ) as a function of carrier frequency f c at seven distinct propeller pitch ratios are plotted separately for the first four shaft frequency harmonics in Figure 11a-d, respectively. The spectral coherence behavior was similar at 82% and 79% pitches, with significantly higher values for the first, third, and fourth harmonics of shaft frequency, which led to the estimation of propeller parameters with high confidence. The behavior at 65% and 31% pitches were similar, but with significantly lower spectral coherence values. The spectral coherence values were low for all pitches at the fourth harmonic of shaft frequency, except for the idealized cases of 82% and 79% pitches.

Robustness of the Estimation of Ship Parameters
Here, we summarize the robustness of the estimation of the ship's parameters from the underwater radiated sound at different propeller pitches.

Engine-Related Parameters
Among the various engine-related parameters, only the cylinder firing rate CFR could be well estimated at different pitches from the differences in adjacent harmonic frequencies since the fundamental CFR component was not detected. Under normal conditions, the engine firing rate EFR can be estimated to be the strongest harmonic component of CFR, whose order number is the number of cylinders in the engine [9,63,64]. However, in this analysis, EFR could not be correctly estimated for this ship since the corresponding eighth harmonic was not the strongest tonal harmonic of CFR. However, the ninth harmonic of CFR was found to be the dominant tonal, both in terms of the highest power and the highest temporal coherence at every propeller pitch, as shown in Figure 6. As a result, the correct cylinder number could not be estimated, but was probably near to nine (the cylinder number provided for the ship under study is eight in Table 3).

Propeller-Related Parameters
Propeller-related parameters, such as shaft frequency SF, blade pass frequency BPF, and blade number, can be estimated from the spectral coherence distributions of shipradiated cavitation noise [4,9,45,46,62,65]. The magnitude of the spectral coherence at one specific frequency represents how dominant the corresponding time period is, at which propeller blades generate cavitation noise. In other words, it indicates how many propeller blades generate cavitation noise at approximately the same time for each rotation cycle and whether these blades generate cavitation noise with uniform time intervals between rotation cycles.
SF can be obtained as the fundamental frequency of the spectral coherence distribution, if at least one blade generates cavitation noise with a uniform time interval. Under perfect conditions, especially when all propeller blades generate cavitation noise at approximately the same time during each propeller rotation cycle, BPF can be estimated to be the second strongest harmonic component after SF. Then, the blade number can be obtained as the harmonic order of BPF with respect to SF. However, under other scenarios, especially when different propeller blades generate cavitation noise at different time instants during each rotation cycle or only few blades generate cavitation noise, then BPF and the blade number cannot be properly estimated, since the spectral coherence at BPF is no longer outstanding among the harmonics of SF [4,62]. (d) Figure 11. The spectral coherences (SC) γ(α, f c ) of the (a) first, α = SF1, (b) second, α = SF2, (c) third, α = SF3, and (d) fourth, α = SF4, harmonics of shaft frequency across the entire carrier frequency f c range at different propeller pitch ratios, ranging from 20% to 82%.
As shown in Figure 9 and Table 5, the shaft frequency SF could be obtained at all propeller pitches analyzed here, which means there was always at least one blade generating cavitation noise with the change of propeller pitch. However, the blade pass frequency BPF and blade number might not be correctly estimated. The most challenging scenario was at 20% pitch, where only the shaft frequency SF could be obtained. This indicates that, at this pitch, only one blade generated cavitation noise, but with a uniform time interval.

Gear-Related Parameters
Gear-related parameters can also be obtained, which is the reduction gear rate RGR. According to Equations (1)-(4), after the shaft frequency SF and cylinder firing rate CFR are estimated, the reduction gear rate RGR can be obtained from the ratio of those two [9,62]. However, there are multiple possibilities due to different engine stroke numbers. For a 2-stroke engine, the reduction gear rate is calculated as RGR = CFR/SF. For a 4-stroke engine, the reduction gear rate is obtained as RGR = CFR · 2/SF = 6.3 · 2/3.384 = 3.723, which matches the gear parameter for the ship under study in Table 3.

Changes in Both Propeller Pitch and RPM
Here, we analyze and discuss one special measurement where the ship changed speed, and both propeller pitch and RPM. The ship's speed changed continuously while propeller RPM underwent a jump over the measurement time window, as shown in Figure 12a,b. The ship continuously increased speed from 5.5 knots to 12.5 knots. The shaft frequency SF hovered around 3.37 Hz for the initial 3/4 of the measurement time window and then increased to 3.8 Hz for the remaining 1/4. The jump in shaft frequency was a result of the jump increase in propeller RPM. Detailed propeller pitch information is not available for this measurement.
The power spectral density (PSD), temporal coherence TC, and mean spectral coherence γ(α) for this special measurement are shown in Figure 13a-c, respectively. The dominant tonal at 56.8 Hz had the largest PSD and TC, similar to previous results, corresponding to the ninth harmonic of CFR associated with a shaft frequency of 3.37 Hz and a lower propeller RPM of 760. Due to the longer measurement time at low RPM (3/4 of the measurement duration), the tonal signals were dominant at frequencies associated with low RPM. The mean spectral coherence distribution, however, displayed harmonic groups corresponding to shaft frequency SF of 3.8 Hz instead, which was associated with a higher RPM of approximately 860. This result implies that the cavitation noise at high RPM was significantly stronger and more dominant than that at low RPM, despite the smaller measurement time at high RPM (only 1/4 of the measurement duration).    Figure 9, with respect to the first four harmonics of shaft frequency SF.
The mean spectral coherence variation with the increase in harmonic order, shown in Figure 13c, was similar to that of Figure 9d, corresponding to 51% pitch, following an "increase-decrease-increase" trend among the first four harmonics of shaft frequency SF. The correlation coefficient between the mean spectral coherence values for the first four harmonics of shaft frequency for the special measurement here, shown in Figure 13c, and at different pitches, shown in Figure 9a-g, were calculated and plotted in Figure 13d. The highest correlation for the special case here was with the 51% pitch. This indicates that the most dominant propeller pitch for the special measurement here could be around 51%. Besides the first four harmonics of shaft frequency, the mean spectral coherence distribution for the special measurement, shown in Figure 13c, also clearly revealed higher order shaft frequency harmonics, ranging from the fifth to the ninth, and twelfth.

Discussion
The underwater sound radiated by a ship with a controllable pitch propeller (CPP) was analyzed and quantified as a function of propeller pitch ratio, including acoustic feature extraction and ship parameter estimation. The ship's acoustic signature was characterized in terms of (i) signal energetics via power spectral density, (ii) machinery tonal sound via temporal coherence, (iii) propeller cavitation noise via spectral coherence, and (iv) frequencymodulated (FM) tonal signal mean and standard deviation in frequency variation. Each of these quantities were calculated from ship-radiated underwater sound measurement at different propeller pitch ratios, ranging from 20% to 80%, and then compared.
The 51% pitch was found to be a crucial point for this ship, because (a) the power spectral density of the radiated sound reached the lowest levels over a broad range of frequencies, (b) the power spectral density and temporal coherence of the ship's significant tonal sounds attained their highest values, and (c) the modulation frequency of the FM tonal signals were the largest. However, the spectral coherences of shaft frequency harmonics for the ship's cavitation noise were among the lowest at 51% pitch.
The efficacy and robustness of the estimation of the ship's parameters under different propeller pitch ratios were investigated. Our analysis indicates that the cylinder firing rate CFR could be well estimated at all propeller pitches from the frequency differences in adjacent CFR tonal harmonics. The engine firing rate EFR and cylinder number could not be accurately estimated since the strongest CFR harmonic was found to be one order higher than that provided for the ship under study. The propeller shaft frequency (SF) could be obtained from the spectral coherence analysis of the cavitation noise at every propeller pitch. The blade pass frequency (BPF) and blade number could be estimated for several pitches, but not all. The reduction gear rate (RGR) could be calculated from the CFR and SF frequencies with an appropriate engine stroke number assumption. The shaft frequency harmonics were most coherent at 79% and 82% pitches, providing the ideal distributions for estimating the blade pass frequency, blade number, and reduction gear rate of the ship.
An analysis of one special measurement where ship changed speed, propeller pitch and RPM over the course of the measurement indicates that the cavitation noise at high RPM was significantly stronger than at low RPM.
These findings elucidate the effects of pitch variation on underwater sound radiated by a specific ship with a controllable pitch propeller, which brings some valuable insights into ship propeller design and operation. However, in order to obtain a much deeper and more accurate understanding of underwater sound from ships with controllable pitch propellers, more analyses should be conducted on a sufficient number of ships. More factors should also be taken into account, including ship type, ship tonnage, ship speed, engine type, propeller shape, propeller size, propeller material, etc.
Author Contributions: Conceptualization, data curation, formal analysis, methodology, and investigation, C.Z., T.G. and P.R.; software and visualization, C.Z.; writing-original draft, C.Z.; writing-review and editing, C.Z., T.G. and P.R.; funding acquisition, P.R.; supervision, N.C.M. All authors have read and agreed to the published version of the manuscript.