# Impact of Ship Motions on Maritime Radio Links

Authors: William Hubert; Yvon-Marie Le Roux; Michel Ney; Anne Flamand
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

The improvement of maritime radio links often requires an increase of emitted power or receiver sensitivity. Another way is to replace the poor antenna gains of traditional surface ship whips by novel antenna structures with directive properties. However, ship motions, especially when launches are involved, may affect the radio link quality. Therefore, a tool for modelling the impact of ship motions on such antenna structures was developed. It helps the specification of antenna radiating parameters and improves the radio link performance evaluation. This tool includes a deterministic two-ray model for radio-wave propagation over the sea surface. The geometrical 3D engine for ship motions has the advantage to be compatible with any propagation model and antenna structure.


---

## Body

## 1. Introduction

This paper presents the development of a deterministic transmission channel model for the improvement of a high data rate (>10 Mbits·s−1) radio-link in the S-band (2 to 4 GHz). This link is supposed to reach the largest possible range between small surface ships. First, the link budget at the approach of the radio horizon by means of the ITU P.526 model [1] was computed. It shows that establishing such a link between ships requires a huge increase in emitted power if poor antenna gains (about 2 dBi) are used. Hence, the purpose of the transmission model is to help the specification of novel antenna structures with gains above 10 dBi. To reproduce the real situation as accurately as possible, especially when ship motions are relevant, a 3D geometrical model was developed. The computational engine can accommodate any antenna and propagation model. However, a simple two-ray model and dipole antennas are used in this paper for illustrative example.After a description of the propagation model, a section is then devoted to the presentation of the geometrical context for the integration of ship motions. The next section is dedicated to the computation of the maximum antenna angular deviation for given sea-wave data. Then, the integration of antenna radiating parameters is explained. Finally, an expression of the received power is derived using channel matrix formalism, followed by a presentation of simulated results.


## 2. Propagation Model

The propagation is described by a two-ray model that includes a direct path between the transmitter (A) and the receiver (B) and a path coherently reflected (in P) by the sea surface (Figure 1). This model has recently been confirmed in [2] by measurements. It is applied here on a spherical earth of effective radius Re which corresponds to a mean atmospheric refractive index. The coherently reflected path is obtained according to the Snell-Descartes law and its amplitude is primarily determined by the Fresnel reflection coefficients:
(1)ρ∥=sinφ-(1/η)η-cos2φsinφ+(1/η)η-cos2φ
for the vertical (or electric field parallel to the incidence plane) polarization,
(2)ρ⊥=sinφ-η-cos2φsinφ+η-cos2φ
for the horizontal (or orthogonal to the incidence plane) polarization.Figure 1
2D geometry of the two-ray model.In the above expressions,η stands for the relative complex permittivity of the sea water surface. It can be defined by the Debye formula applied by Stogryn in [3].The reflected energy is also affected by a divergence factor which produces a larger aperture of the beam after reflection. Its expression is given by Beckmann and Spizzichino in [4] and depends on geometric parameters previously defined in Figure 1:
(3)FD=[1+2r1r2Re(r1+r2)sinφ]-1/2[1+2r1r2Re(r1+r2)]-1/2.
The model validity is primarily bounded by the variations of FD. This factor must not be used after a limiting distance which corresponds to the transition between the two-ray model and the ITU P.526 model. An attempt of integrating the latter does not give enough satisfaction with regard to the transition range. Moreover, the program objective is to assist the antenna specification. So, the design of the link is done at long range by means of the ITU model. Comparison between these different models is shown in Figure 2, obtained at 2.4 GHz, with an effective isotropic radiated power (EIRP) of 1 Watt, with both antennas at 6 meters above the sea surface.Figure 2
Field strength (dB V/m) versus distance (log scale).Expected results are obtained in Figure2 since good vertical and poor horizontal spatial diversities are observed. Besides, it appears worthless to take benefits from the polarization diversity observed at short ranges because of the unfavorable two-ray recombination in horizontal polarization.Finally, a second factor is applied to the specularly reflected energy in presence of sea waves. If a normal distribution hypothesis of sea surface elevations is made, the reduction factorρr depends on their standard deviation σ. According to Miller et al. [5], the reduction factor is given by the following expression:
(4)ρr=exp[-2(2πg)2]·I0[2(2πg)2],
in which(i)
g=(σsinφ)/λ is the Rayleigh criterion,(ii)
I0[·] is the modified Bessel function of the first kind and order zero, and(iii)
λ is the EM wavelength in meters.Figure3 is obtained with the same parameters as before but for a sea state 3 on Douglas sea scale (σ=0.5 meters). A slight swell is sufficient to significantly reduce the coherently reflected energy. The result is a partial vanishing of fading figures. The energy which is no longer coherently reflected leads to a diffuse scattering. However, diffusion is not yet integrated in this deterministic model. Further improvements could be made by introducing that diffuse component in the propagation model. A good starting point would be the analysis proposed in [6] by Karasawa and Shiokawa.Figure 3
Field strength (dB V/m) versus distance (log scale), sea state 3 on Douglas scale.


## 3. Integration of Ship Motions

The previously defined 2D geometry has to evolve to three dimensions to integrate ship motions. A global coordinate system(O,X→,Y→,Z→), in which O is the Earth center, is shown in Figure 4. Angle α is defined according to the Intership distance D (Figure 1) and β stands for their relative orientation.Figure 4
Global coordinate system for ship positions.All coordinates will be further expressed in this global system by means of translation vectors and rotation matrices. Then, a local coordinate system attached to each ship center of mass is defined wherein ship motions are described (Figure5).Figure 5
Local coordinate system for the ship inA, front view of the ship.Regarding the ship motions, one shall define below the dedicated terms related to the local coordinate system:(i)
heave is a translation along z→ axis (vertical motion),(ii)
roll is a rotation about x→ axis,(iii)
pitch is a rotation about y→ axis.The above motions can be simply described by sinus functions. Therefore, they are characterized here by their amplitude and period only. Further improvements could be made by introducing recorded data of real ship motions. Antenna locations are also described in this system by their onboard coordinates with a proper orientation. Finally, a coordinate system attached to the antenna center (A) is introduced according to antenna measurements or simulations (Figure 6).Figure 6
Antenna coordinate system with unit vectors for field components.By expressing all geometric objects in the global coordinate system, antenna position and velocity, rays arrival, and departure angles can be determined at any time.


## 4. Maximum Angular Deviation

Without real ship motion records, approximate geometrical relations can be derived for antenna maximum angular deviation from vertical direction. Knowing the maximum crest-to-trough wave heightHmax and the wavelength λs, one can evaluate the maximum antenna angular deviation:
(5)θmax=arcsin(πHmaxλs2+π2Hmax2).
Several simulations have shown that this angular deviation represents the maximum deviation from the position at rest for distances above 100 meters. For linearly polarized antennas, it also leads to a maximum polarization loss factor (PLF) given by
(6)PLFdB=20log10(cos2θmax).
Using the fundamental mode of the Pierson-Moskowitz spectrum for fully developed regular wind waves [7], a complete data set is defined for the studied case. It includes the period Ts and the standard deviation of sea-wave elevations σ:
(7)Hmax≃5.7m,Ts≃9.2s,λs≃131.4m,σ≃1m.
What is introduced in (5) leads to
(8)θmax=7.6°.
Note that the result on θmax remains approximately constant for any Hmax. This is due to the fully developed regular wind wave assumption. However, if both Hmax and λs are determined by any means (models or measurements), (5) still holds for small ships.


## 5. Antenna Integration

A radiation vector functionF is introduced to account for polarisation state and antenna gain along any direction. As in [8], F depends on the realized gain (which integrates impedance mismatch) for a given set of departure (emitting antenna) or arrival (receiving antenna) angles:
(9)F=G(χ,ψ)U=G(χ,ψ)(Uχ(χ,ψ)Uψ(χ,ψ)).
In the above expression, χ and ψ stand for the elevation angle and the azimuth angle, respectively. U is a normalized vector which corresponds to the complex ratio (amplitude and phase) emitted (or received) along u→χ and u→ψ components (see Figure 6).A database was created for the simulations. It consists of multiple tables of data extracted from electromagnetic (EM) simulators, in our case HFSS. These tables represent realised gain and antenna field components as expressed in (9). Finally, for wide band communications, the radiation vector function database F would have to be generated for several frequencies.


## 6. Expression of the Received Power

First of all, a channel matrix is built by using the unit vectors (see Figure6) expressed in the same coordinate frame. It expresses intrinsically the PLF. Then, the spherical-divergence and phase terms are applied to obtain the direct-path channel matrix:
(10)CDP=1r(u→χA·u→χBu→ψA·u→χBu→χA·u→ψBu→ψA·u→ψB)e-j2πr/λ.
Radiation vector functions for the transmitting antenna FA(χ,ψ) and for the receiving antenna FB(χ,ψ) are then applied to both sides of this matrix to obtain a transmission scalar sDP. Thus, with (9) and (10),
(11)sDP=1rGAGB(UχBUψB)×(u→χA·u→χBu→ψA·u→χBu→χA·u→ψBu→ψA·u→ψB)(UχAUψA)e-j2πr/λ.
For the coherently reflected path, unit vectors u→χ,ψA,B have to be decomposed into a parallel u→∥ and an orthogonal u→⊥ components to the sea surface as shown in Figure 7.Figure 7
Unit vector definition for reflection.Then, one can define the reflection dyadicR¯¯:
(12)R¯¯=FD·ρr·(ρ∥00ρ⊥)
in which ρ∥, ρ⊥, FD, and ρr are defined in (1), (2), (3), and (4), respectively. Thus, the propagation channel matrix for the reflected path becomes
(13)CRP=1r1+r2(u→χB·u→∥ru→χB·u→⊥u→ψB·u→∥ru→ψB·u→⊥)×R¯¯(e→χA·u→∥iu→ψA·u→∥ie→χA·u→⊥u→ψA·u→⊥)e-j2π(r1+r2)/λ.
The same approach used to obtain (11) is applied to derive the second scalar for the coherently reflected path:
(14)sRP=1r1+r2GAGB(UχBUψB)(u→χB·u→∥ru→χB·u→⊥u→ψB·u→∥ru→ψB·u→⊥)×R¯¯(u→χA·u→∥iu→ψA·u→∥iu→χA·u→⊥u→ψA·u→⊥)(UχAUψA)e-j2π(r1+r2)/λ.
Note that the unit vectors u→χ,ψA,B and radiation vector functions in (11) and (14) are not the same because antenna polarization and gain are subject to change with the propagation direction. Finally, using classical formulas, one can evaluate the received power Pr as a function of the supplied power at antenna terminals Pf:
(15)Pr=(λ4π)2Pf|sDP+sRP|2
in which sDP and sRP are given by (11) and (14), respectively. In the above expressions, time (related to ship motions) and frequency (for electromagnetic wave propagation) are involved.


## 7. Simulated Results

A simulation with a transmitter which undergoes the previously defined motions (sea state 5 on Douglas scale) and a receiver ashore (motionless) is proposed. Both antennas are at six meters above the sea surface with a distance of 700 meters between them. As shown in Figure2, this range naturally leads to a ray recombination equal to the free-space path loss. Both antennas are dipoles with maximum gain of 2.1 dBi. Variations of the received power versus time are shown in Figure 8. Different simulations have shown that those variations are mainly due to the two-ray recombination which changes with movements. On the other hand, significant losses are observed if antennas with relatively high directivity are used. To illustrate this, Figure 9 shows the gain at transmission for the direct path and for antenna arrays which have higher directivity. Note that similar results are obtained for the reflected path. The maximum gain of the linear vertical array of 8 dipoles is greater than 9 dBi, but its half-power beam aperture angle is less than 2θmax. As a result, instantaneous effective gain fluctuates to values below the ones obtained by linear vertical array of 4 dipoles. This confirms that the −3 dB elevation aperture angle has to be at least in the order of 2θmax (given by (5)) to preclude large gain reductions. Note that ship motions in the azimuth plane have negligible impact on the gain variation.Figure 8
Received power versus time, sea-state 5, ship motions in the incidence plane.Figure 9
Gains at transmission on the direct path with dipole and arrays of dipoles.Finally, several simulations show that ship motions lead to a polarization loss factor less (PLF) than 3 dB. This indicates that using circular polarization does not bring any advantage because of the deep fadings produced by horizontal polarization (see Figure2).


## 8. Conclusion

The description of a maritime radio link channel simulator was presented. It includes a two-ray propagation model, antenna parameters, and three-dimensional (3D) ship motion model. Each component of the simulator is independent as the ship motion calculation engine can be used with any propagation model or antennas. The objective of the simulator is to provide a reliable tool for antenna specification and channel performance evaluation. First results show that classical behaviour of the two-ray model remains the main contribution to the channel characteristics if the condition on the maximum deviation angleθmax is fulfilled. Antenna −3 dB elevation beam aperture has to be chosen according to ship motions, given a radio-link-operational sea state. In addition, the use of vertically polarized antenna seems to be the most advantageous option. The maximum deviation angle expression is provided and, for small-ship motions, both the polarisation loss factor and the decrease in antenna gain can be predicted. For larger ships, the only parameter that has to be known is the maximum angular deviation between roll and pitch.



---

*Source: 507094-2012-12-06.xml*