Spread Spectrum Modulation with Grassmannian Constellations for Mobile Multiple Access Underwater Acoustic Channels

The objective of this study is to evaluate Grassmannian constellations combined with a spread spectrum multiple access scheme for underwater acoustic mobile multiple access communication systems. These communication systems enable the coordination of a fleet of Autonomous Underwater Vehicles (AUVs) from a surface or bottom control unit, e.g., a boat. Due to its robustness against phase rotation, the demodulator of Grassmannian constellations uses non-coherent detection, and the main advantage of such modulation lies in the spectrum efficiency gain with respect to conventional differential modulation. The communication system under study in this paper consists of (i), at the transmitter side, a Grassmannian modulation used in an orthogonal spread spectrum multiple access scheme called Multiuser Hyperbolic Frequency Modulation (MU-HFM) and (ii), at the receiver side, a non-coherent array decoder. The modulation and demodulation are presented as well as the considered spreading sequences. Finally, performances of the proposed transmission scheme are evaluated over replayed underwater acoustic channel responses collected at sea by a multi-sensor acoustic acquisition system.


Introduction
Digital communication systems for Underwater Acoustic (UWA) channels are traditionally affected by long multipath channels (e.g., caused by the reverberation of the acoustic wave on the bottom or the surface of the water), a large Doppler spread, and movement-induced Doppler shifts. The speed of sound in water (c = 1500 m·s −1 ) is five orders of magnitude lower than the speed of light in air, and this causes long propagation delays and multipath channels. The communication system can be considered as ultra-wideband, since the center frequency is comparable with the available bandwidth (as an example, in this paper, we consider a center frequency of 27 kHz with a bandwidth equal to 4 kHz). Sound attenuation in water is frequency dependent and time-varying, and together with background acoustic noise, they limit the achievable data rate considerably as indicated in [1,2]. Because of these drawbacks, realizing a communication system that makes a fleet of Autonomous Underwater Vehicles (AUVs) a mobile connected network is a challenging task. An example of the considered scenario consists of multiple AUVs that communicate with a receiver array deployed over a surface control unit, which could be a boat or a buoy. In the following, we consider the access channel between a fleet of N u AUVs transmitting data and a receiver equipped with N r hydrophones situated at the sea surface.
An extensive multipath effect causes inter-symbol interference (ISI) in the received signal that can be addressed either by non-coherent modulation such as Frequency Shift Keying (FSK) [2] or coherent modulation associated with an advanced equalization scheme, like the Decision Feedback Equalizer (DFE) as proposed in [3][4][5]. Another approach to deal with Inter Symbol Interference (ISI) is the use of spread spectrum communication combined with Differential Phase Shift Keying (DPSK) modulation at the price of data rate reduction [6]. This loss in the data rate is accentuated by the insertion of a preamble sequence often required to estimate time-varying Doppler shifts or to detect the start of a data frame. In [7,8], a preamble called a Dual-HFM signal, which consists of up-chirp Hyperbolically Frequency Modulation (HFM) and down-chirp HFM, is used for Doppler shift estimation. In [7], Doppler shift estimation is done by matched filtering at the receiver. This makes it possible to find the positions of the 2 HFM signals. A method based on speed spectrum scanning, which consists of starting from a set of candidate speeds, is used in [8] in order to improve the accuracy of the Doppler shift estimation. Several variants exist on the combination of the preamble, as, for example, in [9], where the preamble is composed of chirp Linear Frequency Modulation (LFM) and chirp HFM. However, the methods cited above do not apply in the case of a multi-user scenario. In the rest of the paper, we will use the method given in [10], which consists in using several filter banks to find the Doppler shift estimation.
To be able to detect and correctly decode multiple users at the receiver, traditional orthogonal multiple access schemes are Time-Division Multiple Access (TDMA), where a time slot is assigned to each user with a guard period, or Code-Division Multiple Access (CDMA), where a specific spreading code is assigned to each user [11]. As an alternative to CDMA, we introduced in [12] a novel multiple access scheme, called MultiUser Hyperbolically Frequency Modulation (MU-HFM), which is based on a set of mutually orthogonal chirp-based waveforms that are robust with respect to Doppler and multipath effects. The objective of this proposed scheme is, on the one hand, to benefit from the robustness of chirps against UWA channel impairments and, on the other hand, to use orthogonality to separate multiple users at the receiver side, using a simple matched filter. However, the MU-HFM scheme proposed in [12] employs DPSK modulation that exhibits low spectral efficiency. In order to increase the data rate of the communication system, we propose in this paper to consider non-coherent Grassmannian modulation whose codeword is generated by the Cube-Split quantizer that was originally introduced in [13,14]. Several other Grassmannian modulations exist in the literature, as in [15,16]. However, this type of modulation requires that the modulation symbols at the receiver are stored. In this paper, we will consider a variant of this modulation scheme that has a lower demodulation complexity at the receiver. In the following, we compare the performance obtained with the Cube-Split modulation with that obtained by a differential phase shift keying modulation that has the same spectral efficiency. We then increase the order of the Grassmannian modulation in order to achieve a higher spectral efficiency, and therefore a larger amount of bits per symbol, to then compensate for the low spectral content typical of spread spectrum communication schemes. Higher order Grassmannian modulations require higher robustness against inter-symbol interference. This is achieved by exploiting the channel spatial diversity processing gain provided by the hydrophone array at the receiver. To the best of the authors' knowledge, there is no previously published work presenting Grassmannian modulations applied to UWA communication systems, so the main contribution of the paper lies in the use of Grassmannian modulations as an alternative of differential modulation for non-coherent UWA communication in the context of an underwater multiple-user mobile network.
The paper is organized as follows: in Section 2, we introduce the system model and detail the MU-HFM multiple access scheme. The construction of the Grassmannian modulation based on the Cube-Split quantizer and the associated decoding process is summarized in Section 3. In Section 4, we provide the performance obtained by using the underWaterAcousTic channEl Replay benchMARK (Watermark) channel simulator [17], where channel impulse responses collected during shallow water experiments conducted in the roadstead of Brest, France, are used. Finally, conclusions are drawn in Section 5.
This article is an extension of the conference paper [12] with the following additional features: • A study of a new Grassmanian modulation [13,14] combined with a spread spectrum communication in an UWA channel and a comparison with differential modulations of the DPSK type. • The implementation of Doppler shift estimation and frame synchronization processes at the receiver side for more realistic communication and a comparison of the MU-HFM spreading sequence against Pseudo-Noise (PN) sequence as preamble for Doppler shift estimation.

Mathematical Notations
In the following, j denotes the unit imaginary number, || · || 2 the Euclidean norm, · the scalar product, (·) * the complex conjugate, (·) T the transpose, and u * v the convolution product between u and v.

Transmitter
We assume a multi-user spread-spectrum transmission scheme with a baseband transmit signal per user i ∈ [1, N u ] expressed by with N s as the number of symbols, and d i,k as the k-th modulated symbol. each modulated symbol is assumed to be a preliminary bit interleaved and encoded with a Forward Error Correction (FEC) code of rate R C in order to increase robustness against channel impairments. For user i , T s is the waveform duration, which, in the case of a spread spectrum scheme, also corresponds to the symbol duration, and g i (t) is the spreading signal. In order to provide robustness against the distortions within the UWA channel, the spreading signals are constructed from the MU-HFM scheme, originally presented in [12], that consists of HFM-based signals combined in a way to ensure mutually orthogonality. The starting point of the construction is a family of narrow-band HFM signals noted by {c i (t)} N u i=1 and defined as follows: is the signal slope, f l and f h are, respectively, the start and stop frequencies, where B denotes the signal bandwidth, and T s the waveform duration. The parameter ζ ∈ R influences the bandwidth of the different c i (t) and is determined by using Simpson's method [18] in order to ensure mutual orthogonality between c i (t) values [12].
Finally, the different c i (t) values are combined with an HFM signal over the entire bandwidth. This makes it possible to have a larger BT s time-frequency product and to have a better correlation peak for synchronization [9]. To keep the mutual orthogonality across c i (t), we use the Gram-Schmidt process [19] given by where and e 0 (t) is a full bandwidth HFM signal defined as Figure 1 shows an example of auto and cross-correlation functions for the signals g i (t), whereas Figure 2 provides a diagram of the transmitter processing.

Receiver Architecture
A global diagram of the receiver is shown in Figure 3. The acoustic signals sent by the N u users, received by an array of N r hydrophone sensors such as the signal recorded in baseband at the p-th sensor, can be expressed as [21] where h i,p (τ, t) is the time-varying channel impulse response between the i-th user and the p-th receiver, f c denotes the center frequency of the passband transmitted signal, a i,p = v i,p c is the motion-induced Doppler scale factor for the i-th user, with v i as the relative speed of the i-th user with respect to the receiver p, and c 1500 m/s as the speed of sound in water. Finally, n p (t) is the assumed additive white Gaussian noise for the p-th receiver.
If the Doppler shift can be estimated, the Doppler effect is usually removed prior to decoding by resampling the received baseband signal and compensating the phase rotation as follows [22]: whereâ i,p denotes an estimation of the Doppler shift. By assuming perfect time symbol synchronization, the k-th information data k ∈ [1, N s ] of the i-th user for the p-th receiver can be estimated by matched filtering z i,p (t) with the transmit waveform of user i as follows: where γ i,k,p denotes the bias of the decoder, η i,k,p is the interference terms, and w i,k,p is the additive noise term for the p-th receiver; the exact expression of these three terms is provided in Appendix A. The final symbol estimation and decision is obtained by combining the N r received signals as follows: As explained in [23], if the modulation duration T s is greater than the UWA channel delay spread, the inter-symbol interference is negligible. On the other hand, due to the mutual orthogonality property of g i (t), the multi-user interference is also limited to obtain η i,k,p terms that decay towards 0 as T s tends to +∞. In [12], data symbols d i,k belong to a DPSK modulation constellation in order to avoid the use of channel-estimation-based equalizers at the receiver side to compensate for the channel coefficients γ i,k,p . A UWA communication channel may have rapid time variations, and differential modulation schemes have been demonstrated to provide good reliability vs. spectral efficiency trade-offs and, in some conditions, even outperform coherent modulation [24]. However, the main drawback of a DPSK modulation lies in its lower robustness against additive noise with respect to a Phase Shift Keying (PSK) modulation that limits its spectral efficiency at the same Bit Error Rate (BER) [25]. In the following, we will consider an alternative to the DPSK constellation by considering a Grassmannian constellation.

Grassmannian Constellations
In the literature, many Grassmannian constellation designs have been proposed to define the different d i,k [26]. Each Grassmannian constellation has its own method of symbol generation and decoding. In general, a good constellation should consist of codewords maximally distant from each other, while also ensuring sphere-packing properties, which translates in evenly spread symbols over a Grassmannian manifold. For example, in [27], the distance between codewords is maximized by a numerical optimization. This optimization generates a lack of structure in the modulation that requires that the codeword be stored on reception. To avoid this storage, codewords should be generated to have a given structure. For example, in [13,14,26], a new Grassmannian-structured constellation is generated through Cube-Split algorithms, whereby the Grassmannian space is partitioned into cells.
The main advantage of this new constellation lies in its robustness against single tap channel fading that causes such modulation to be decoded without Channel State Information (CSI) at the receiver and as a consequence without a complex equalization stage as in, for example, Quadrature Amplitude Modulation (QAM). Moreover, the advantage of Grassmannian constellation with respect to conventional Differential PSK modulation lies in its higher spectral efficiency that is achieved by a multi-level constellation format. To begin with, the encoded bit stream will be split into two parts. One part represents the number of the different cells of the Grassmannian space, while the other part will provide a set of coordinates corresponding to each cell of the Grassmannian space. The final constellation symbol is given in the form of a vector with M components.
In the following, we will detail the Cube-Split modulation design and the decoding process for a spread-spectrum communication scheme. In [26], the communication channel is a flat Rayleigh fading channel. In our study, we consider the underwater acoustic channel that is considered to be a non-flat fading channel [28], and the different users will be moving at a speed unknown by the receiver.

Transmitter Modulation Scheme Design
Let G(C M , 1) be the Grassmannian manifold with M = 2 n and n ∈ N * . A Grassmannian constellation symbol will be defined by a vector to remain consistent with Equation (1), where i ∈ [1, N u ] is the user number and k ∈ [1, N s ] is the processed constellation symbol. Each user will have the same modulation scheme and constellation. The general idea of the Cube-Split constellation is to partition the G(C M , 1) space into M cells and to define, in each of these cells, a local coordinate system using a bijective mapping. A Grassmannian constellation symbol is then defined by N = log 2 (M) + 2(M − 1)L 0 bits, where log 2 (M) bits indicate the cell index, and L 0 ∈ N * represents the number of local coordinates. The spectral efficiency of the Cube-split modulation is given by In the following, we give an overview of the Cube-Split modulation. For more details, the reader can refer to [13,26]. The cell l ∈ [1, M] of G(C M , 1) is defined by where is the chordal distance between x and ξ l .
To simplify Equation (13), we define a center for each cell using the canonical basis vectors. Let {e 1 , . . . , e M } be the centers of the cells with e l as the M × 1 canonical basis vector, with 1 at position l and 0 elsewhere. With this choice of vectors for the centers of the cells, Relation (13) can be simplified as [26] The codewords of the constellation are mapped in the different cells C l with a bijective mapping whose definition is given in Appendix B. ∀i ∈ [1, N u ], ∀k ∈ [1, N s ] : where is the Cartesian product, and A j is a subset of the interval [0, 1] defined ∀j ∈ [1, 2(M − 1)] by [26]: with L 0 ∈ N * representing the number of local coordinates. The constellation can then be defined ∀i ∈ [1, In spread-spectrum communication, with d i ∈ C N s M×1 and ∀q ∈ [1, N s M], d i,q = d i,k,l .
The expression of the signal transmitted in baseband for the Cube-Split modulation is then given by We can see that the difference comes from the upper bound of the sum, which is, in the case of the Cube-Split modulation, equal to N s M.

Receiver Demodulation Scheme Design
For the demodulation process, Equation (8) represents the output after matched filter. The receiver knows the parameters M and L 0 , which represent the modulation order and the number of local coordinates of the Cube-Split modulation. We can then define for the p-th receiver and ∀i ∈ [1, N u ] to remain consistent with Equation (19) andd i,p ∈ C N s M×1 . To find the symbol transmitted d i,k,p with k ∈ [1, N s ], we proceed in two steps: we decode first the cell number and then the local coordinates.
We begin by calculating the left singular vector u i,k,p = [u i,k,1,p , . . . , u i,k,M,p ] T corresponding to the largest singular value of the received signald i,k,p on the p-th hydrophone, which is also the solution for k ∈ [1, N s ] of arg max The cell index estimate is given by Equations (13) and (14): The local coordinates are estimated by first applying the inverse mapping φ −1 l to havẽ The decoded symbold i,k,p is then identified from the estimated parameters {l i,k,p ,q i,k,p } asd i,k,p = φ −1 l i,k,p (q i,k,p ). The final estimate symbol is obtained by combining the N r receive

Results and Discussion
The channel sounding experiments took place during summer 2019 in the roadstead of Brest, France, by using the IROMI-LMAIR platform [29]. At the transmitter side, one transducer immersed at a 2 m depth was used. At the receiver side, 5 hydrophone sensors were vertically deployed with a spacing of 1 m. Five channel soundings were performed at different transmission ranges D i between 65 and 540 m by using a 511-Maximal Length Sequence (MLS) as a probe signal [30] centered on f c = 27 kHz over a 4 kHz bandwidth. The different channel soundings are not calibrated with respect to propagation loss as a near-far effect. Time synchronization among the hydrophone channels is preserved. Figure 4 provides an example of the delay-Doppler spread extracted from the successive estimated CIR. Estimated channel delay and Doppler spreads are reported in Table 1.

Watermark Replay Channel
To emulate a realistic transmission experiment, we take as a basis the Watermark simulator [17] that is a replay channel simulator driven by measurements of time varying CIR. The principle of the simulator consists in distorting input waveforms by convolving them with measured channels. To simulate a realistic mobile multiuser communication, we considered the following procedure: (i) replay sounded channels on the pass-band signal transmitted by each user s i (t) by using the Watermark simulator, (ii) simulate the Doppler effect linked on the relative motion of each user, (iii) delay and scale the signal according to the relative range of each user in order to represent attenuation and delays of each user signal, and (iv) sum all signals obtained in order to represent a multi-user transmission. As a result, the final signal obtained on the p-th hydrophone can be expressed in baseband as where the Doppler shift is artificially added in the receive signal by resampling and phase rotating the transmitted signal s i (t),ĥ i,p (τ, t) denotes the recorded CIR of the i-th user to the p-th hydrophone,τ i is the communication delay between the i-th user and the receiver, and n p (t) is the additive noise on the p-th hydrophone assumed to be white, zero mean, and Gaussian distributed. γ i represents the propagation losses that generate the near-far effect in the case of multi-user communication. The value is assumed to be the same for each receiving hydrophone and is given by [2] where A(D i , f c ) denotes the acoustic power attenuation in decibels at range D i and frequency f c : with m as the spatial dispersion factor that is assumed to be cylindrical (m = 1), and α( f c ) is given by Thorp's formula [31]. The Signal-to-Noise Ratio (SNR) in dB is then defined ∀p ∈ [1, N r ]: In the following, we will take into account the problems of frame synchronization and Doppler shift estimation at the receive side. The estimation of both the user frame starting time and the Doppler shift is performed owing to a preamble known by the receiver. More especially, Doppler estimation is obtained with a Doppler-bank matched filter, cycling through velocities from [−2, 2] m/s in steps of 0.1 m/s [10]. The filter is matched to the signal preambles that are chosen as a PN sequence or MU-HFM waveform. Below, we recall the expression of PN sequence-based waveforms: with [c i,1 , c i,2 , . . . , c i,N SFpr ] as the spreading code of length N SF pr , where T c is the chip duration, N SF pr is the spreading factor, and φ(t) is the pulse shaping filter chosen as a Square Root Raised Cosine (SRRC) filter [25].
Transmission system parameters are provided in Table 2, where parameter ζ of MU-HFM is computed according to the signal bandwidth by using Simpson's method. For B = 4 kHz, we choose [23]

Performance Metrics
In order to measure the efficiency of the proposed transmission scheme, we choose as a performance metric the average effective data rate defined as follows [23,25]: where η denotes the spectral efficiency of the considered modulation defined as where R C denotes the channel coding rate, and FER is the Frame Error Rate. In the following, a frame is considered erroneous when at least one bit per frame after channel decoding is erroneous. As a benchmark protocol, we consider a conventional deterministic multiple access scheme for the UWA channel, i.e., TDMA and CDMA. In order to provide a fair comparison, the TDMA protocol is combined with Direct Sequence Spread Spectrum (DSSS) modulation with a spreading factor N SF equal to the CDMA one and equal to the time-bandwidth product of MU-HFM waveforms, such as the performance of all protocols that are comparable in the single user scenario. In order to prevent intersymbol interference, the symbol duration T s and the TDMA guard interval time T g are both chosen to be greater than the maximum duration of the various channel delays reported in Table 1. For each protocol, we consider the conventional Differential Quaternary Phase Shift Keying (DQPSK) constellation and the proposed Grassmannian modulation. The average effective data rate per user for CDMA and TDMA is given by Finally, to evaluate the accuracy of the Doppler shift estimationâ i,p , we use the Root Mean Square Error (RMSE) metric defined ∀i ∈ [1, N u ], ∀p ∈ [1, N r ] as

Performance Results
Compared to [12], in this part, we provide the performance results of the 3 studied protocols with and without Grassmannian modulation by using the UWA replay channel described before and by considering mobile AUVs, whose motion is emulated by adding the motion-induced Doppler scale at the output of the Watermark channel. For each frame, the speed value of each AUV is randomly selected in the interval [−2, 2] m/s and assumed constant along a frame. Moreover, for all simulations, frame synchronization and Doppler shift estimation are performed at the receive side in order to be as close as possible to a real experiment. Figure 5 shows the RMSE of the Doppler shift estimation as a function of the number of users. As a reference, we plot, in a green line, an estimation error of 0.1 m/s, leading to a Doppler shift error of 0.6 × 10 −4 , which represents the search step of the Doppler bank algorithm. We can see for a simultaneous communication that a preamble based on a PN sequence is not efficient since, from 2 users, the estimation is strongly degraded. This is due to the low resistance of the PN sequence against channel impairments. However, by using a MU-HFM-based preamble, the accuracy is satisfactory up to 4 users. At 5 users, the Doppler estimation is substantially erroneous, yielding a strong degradation of the decoding performance. As a reference, we can see that the TDMA-PN curve representing a PN preamble-based Doppler estimation in the TDMA protocol case leads to perfect Doppler estimation, regardless of the number of users. In Figure 6, FER and effective data rates results are presented by assuming a single stream receiver. For each protocol, we can see that DQPSK and Grassmannian modulation with M = 4 provide the same FER. Logically, by increasing the parameter M of the Grassmannian modulation to M = 8, the effective data rate is increased, since the number of local coordinates is also higher. However, the Grassmannian modulation requires precise decoding based on the index of the cell and on the local coordinates associated with the cell, so the increase in M makes the modulation less robust to noise and interference leading to higher FER. The CDMA-PN curve that represents the performance of CDMA with a PN sequence as a preamble are provided in order to show the ineffectiveness of this type of preamble in the case of simultaneous communication. As shown in Figure 5, this is due to the very low accuracy in the estimation of the Doppler shift. By comparing protocols associated with MU-HFM preambles, we can see that the maximum number of simultaneous users is 4. Beyond that, the Doppler shift estimate is too erroneous, as shown in Figure 5. The CDMA approach appears less efficient beyond 3 active users, whereas the traditional TDMA provides the best performance at 5 users. This poor performance of CDMA compared to that of MU-HFM is explained by the fact that, in MU-HFM, we have achieved an orthogonal combination using the Gram-Schmidt method between an HFM signal and a set of narrowband chirps. From 4 users, the MU-HFM has an effective data rate that is decreasing due to the higher multiuser interference terms yielding a higher FER. For the Grassmannian modulation with M = 8, we can see the MU-HFM still has the best performance compared to CDMA for the same reasons as found previously. In Figure 7, results are presented with 5 hydrophones at the receiver side and multichannel processing, as shown in (11). As for the system with N r = 1, the maximum number of users is 4 for protocols with a MU-HFM type preamble. For a preamble with a PN sequence, CDMA is only possible for 1 user. However, the preamble with a PN sequence allows TDMA communication up to 5 users because the estimation of the Doppler shift is perfect. In the case of the CDMA protocol with an MU-HFM preamble, the FER of the Grassmannian modulation with M = 8 is slightly improved for 2 users. In the case of MU-HFM protocol, the FER performance of both differential and Grassmannian modulations are substantially improved, yielding error free transmission up to 4 users with DQPSK or Grassmannian modulation with M = 4. With M = 8, error free transmission is reached for up to 4 users and again at 5 users. The transmission is ineffective due to the poor estimation of the Doppler shift. By comparison with Figure 6, this demonstrates that the MU-HFM protocol takes full advantage of the spatial processing gain provided by the multi-sensor architecture. On the other hand, the MU-HFM combined with Grassmannian modulation with M = 8 offers a constant user data rate of 10 bps for up to 4 mobile users, which demonstrates the superiority of the proposed MU-HFM protocol in comparison to conventional protocols. These performance results also show the spectral efficiency gain of Grassmannian modulation against differential modulation. In fact, in order to reach a similar data rate, the MU-HFM should be combined with a 16-state DPSK modulation that is clearly not feasible in the context of UWA communication.

Conclusions
In this paper, we have considered the use of a Grassmannian modulation, named Cube-Split modulation, as an alternative to differential modulation for multiuser communication in a scenario where an AUV fleet transmits acoustic data to a surface receiver within the same bandwidth. To reduce the detrimental effects of the UWA channel, we have considered spread-spectrum-based communications-more specifically, MU-HFM waveforms designed for multiple access in the UWA channel. Experimental results with the Watermark channel fed by sea channel sounding demonstrate that Grassmannian modulation outperforms conventional quaternary differential modulation when multi-hydrophone reception is considered. In particular, the MU-HFM offers a quasi-constant effective data rate as the number of active users increases, up to a maximum of 4 users. Beyond that, the estimation of the Doppler shift is no longer possible for a preamble composed of the MU-HFM waveform. Consequently, the use of Grassmannian modulation is demonstrated to be an efficient alternative of differential modulation that cannot exceed 8 states in practice. Indeed, to reach an equivalent data rate, coherent QAM modulation would be required, which implies an equalization step that eliminates the interference between symbols. However, such an approach would require precise estimates of multiuser channels based on a pilot sequence designed with respect to the channel coherence time that would limit the number of simultaneous users. Funding: This research was partly funded by the Brittany region council under the grant ARED-2017-161. The experimental platform used for this study was partly funded by the European Regional Development Fund (ERDF) under the grant EU000897. This work is also part of the SFI Smart Ocean (a Centre for Research-based Innovation). The Centre is funded by the partners in the Centre and the Research Council of Norway (project no. 309612).

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Conflicts of Interest:
The authors declare that they have no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: In Section 2, the received baseband signal after Doppler pre-processing at the p-th receiver can be expressed as (A8)