Composite Backscatter Characteristics of Conductive/Dielectric Ships and Sea Surfaces with Breaking Waves under High Sea Conditions

When a radar detects marine targets, the radar echo is influenced by the shape, size and dielectric properties of the targets, as well as the sea surface under different sea conditions and the coupling scattering between them. This paper presents a composite backscattering model of the sea surface and conductive and dielectric ships under different sea conditions. The ship scattering is calculated using the equivalent edge electromagnetic current (EEC) theory. The scattering of the sea surface with wedge-like breaking waves is calculated using the capillary wave phase perturbation method combined with the multi-path scattering method. The coupling scattering between ship and sea surface is obtained using the modified four-path model. The results reveal that the backscattering RCS of the dielectric target is significantly reduced compared with the conducting target. Furthermore, the composite backscattering of the sea surface and ship increases significantly in both HH and VV polarizations when considering the effect of breaking waves under high sea conditions at low grazing angles in the upwind direction, especially for HH polarization. This research offers valuable insights into optimizing radar detection of marine targets in varying sea conditions.


Introduction
Composite scattering modeling for sea surfaces and targets is of great significance for monitoring sea environment, identifying and intercepting targets [1,2]. The difficulties of the coupling scattering mechanism between the sea surface and the target are primarily twofold. Firstly, with the development of stealth technology, radar stealth can be achieved by using coated absorbing materials to reduce the signals of the target within a certain range. In this case, the target is no longer a perfect conductor, necessitating a study of the scattering characteristics of dielectric targets. For complex dielectric targets with wedge structures, the scattering problem becomes more complicated due to the diffraction coefficient determined by coupling the electric and magnetic fields for the dielectric wedge. Secondly, with advances in shipbuilding and navigation technology, some types of ships can sail under high sea conditions, making it necessary to fully consider the impact of breaking waves on the scattered echo. Therefore, it is very urgent and necessary to study the composite scattering characteristics of the sea surface and the target, particularly the composite scattering characteristics of the sea surface and the dielectric target under high sea conditions.
There have been many studies on sea surface scattering [3][4][5]. Research on the composite scattering of sea and target has primarily involved numerical acceleration algorithms, such as the generalized forward and backward iteration algorithm (GFBM) [6], the double iteration model of locally coupled field technology [7], the multi-layer fast multi-pole method [8], high-low frequency hybrid algorithms such as KA and MoM (KA-MoM) [9,10], etc. The high-frequency approximation method has been widely used for calculating the composite scattering of an electrically large sea surface and targets with high computational efficiency, mainly including the reciprocity theorem, four-path method, bouncing ray method, etc. Chen and Zhang simulated scattering echo and SAR imaging of a ship above a rough sea surface using a modified four-path method [11]. Yan simulated the imaging of a ship above the sea with single-frequency and single-view using the bouncing ray method [12]. He used the physical optics and iterative physical optics (PO-IPO) hybrid method to calculate the electromagnetic scattering characteristics of a random rough sea and objects above it [13]. Li analyzed the electromagnetic scattering from electrically large ship targets in a marine environment via the combining of the factbased asymptotic method (FBAM) and the geometrical optics and physical optics (GO/PO) hybrid method [14].
Previous studies were dedicated to analyzing the scattering from mental targets or dielectric targets, but without consideration of the diffracted field for a dielectric wedge. On the other hand, sea-only scattering is not competent to fully describe the couple scattering mechanism between the sea surface and breaking waves under high sea conditions. There have been few studies on the composite scattering of a sea surface with breaking waves and the dielectric target above it under high sea conditions. In this paper, the deterministic distribution of breaking waves on the sea surface under high sea conditions is established, and the scattering of the wedge-like breaking waves and the multiple scattering between the breaking wave and the surrounding sea surface are calculated. The numerical matching technique is introduced to quickly calculate the scattering field of electrically large dielectric targets with wedge structures. Combined with the modified four-path method, the composite scattering characteristics of conducting and dielectric ships under different sea conditions are calculated. The influence of breaking waves and the dielectric property of the target on composite scattering of a sea surface and ship under high sea conditions are analyzed. The results provide valuable insights for improving radar efficiency in detecting sea surface targets.

Calculation Method of Electromagnetic Scattering from Conducting/Dielectric Targets EEC Method for Conducting Targets
The total field of the conducting target can be obtained by adding a physical optical field and diffraction field, namely GTDEEC=POEEC+PTDEEC. POEEC is calculated according to the integration direction selected by Cui and Wu [15], and PTDEEC is obtained based on Michaeli's work [16]. The far scattering field E s can then be obtained according to EEC theory and the Gordon integral formula: where k is the wavenumber of the incident wave and η 0 is the wave impedance. The vector w is related to the incident wave vectork i and scattered wave vectork s and satisfies w =k s −k i . r 0 = (r 2 + r 1 )/2 is the central position of the wedge, with which r 1 and r 2 are the two vertices of the wedge. l is the wedge length.t is the tangential vector of the wedge. I(r ) and M(r ) are the equivalent edge electromagnetic currents. The total scattering field of the impedance wedge includes the physical optical field and edge diffraction field. When the equivalent edge electromagnetic current is used to solve the edge diffraction field, the incremental length diffraction coefficient D f of the edge wave needs to be solved first. The solution needs to deduct the contribution of the PO diffraction coefficient, so the PO diffraction coefficient needs to be derived from the physical optical field first. The local coordinate system {x,ŷ,ẑ} is established on the upper surface of the impedance wedge with an exterior angle nπ, as shown in Figure 1. fi is the angle between the incident vector, and the tangential direction of the wedge and the angles of incidence and of diffraction, measured with respect to the 0 face, are denoted by φ i and φ, respectively.
The total scattering field of the impedance wedge includes the physical optical field and edge diffraction field. When the equivalent edge electromagnetic current is used to solve the edge diffraction field, the incremental length diffraction coefficient f D of the edge wave needs to be solved first. The solution needs to deduct the contribution of the PO diffraction coefficient, so the PO diffraction coefficient needs to be derived from the physical optical field first. The local coordinate system {ˆˆx, y,z } is established on the upper surface of the impedance wedge with an exterior angle π n , as shown in Figure 1. β is the angle between the incident vector, and the tangential direction of the wedge and the angles of incidence and of diffraction, measured with respect to the 0 face, are denoted by φ i and φ , respectively.
(2 ) n − π The PO diffraction coefficient matrix where η + is the impedance of the upper face of the wedge and p is the introduced coefficient to simplify the formula, which has the form where η + is the impedance of the upper face of the wedge and p is the introduced coefficient to simplify the formula, which has the form The PO diffraction coefficient of the lower wedge can be obtained using the diffraction coefficient of the upper wedge with transforming η When calculating the diffraction coefficient of the impedance wedge, we use the numerical matching technique to expand the spectrum function of the impedance wedge with oblique incidence and arbitrary wedge angle according to [17]: where Ψ e,h (α) is the traditional Maliuzhinets' function, and σ φ i (α) = 1 n sin φ i n /(sin α n − cos φ i n ). The auxiliary function Λ e,h (α) is as follows: where B e,h 1,2 (α) is the matching coefficient [18]. Then, the impedance wedge diffraction coefficient matrix D(φ, φ i , ρ, α 1 , α 2 ) can be written as where nπ denotes the exterior wedge face, is the special function mentioned in Equation (3) and B(α) is the matching coefficient mentioned in Equation (4).
, the equivalent edge currents I and M can be derived as where η 0 is the impedance of the wedge. E i z and H i z are the longitudinal components of incident electric and magnetic fields, respectively.
The diffraction field of the wedge is obtained by applying the electromagnetic current radiation integral formula, and the total scattering field can be obtained by superposing the physical optical field of the target and diffraction field of the wedge. The scattering field of a single surface facet can be expressed as [11] E f acet where S PQ (k i ,k s ) is the scattering amplitude. Since the main contribution to the radar echo is the Bragg wave propagating near to and far away from the radar, the contribution of the Bragg wave components in the positive and negative directions is considered, respectively, and the scattering amplitude of the small facet can be expressed as: where ∆S is the area of each single facet, n z is the z component of the normal direction of the small facet, F PQ is polarization factor, J n (·) represents the Bessel function of the first kind of order and B(k c ) is the amplitude of the capillary wave. It can be defined as the corresponding energy amplitude in the wave spectrum of its resonant wave number k c , that is, where S capi E (k c ) is the Elfouhaily capillary spectrum. The integral term I 0 (k c ) in Equation (8) reflects the phase modulation of the capillary wave, and its expression is as follows: The length of the simulated two-dimensional sea surface in x g and y g directions is L x and L y , the area is A = L x L y , the number of sample points is M and N and the distance between adjacent two points is ∆x g and ∆y g , respectively. Ignoring the multiple scattering effect between each small facet, the total scattering field of the sea surface can be written as

Wedge-like Breaking Waves Scattering
In this paper, the deterministic distribution of breaking waves on the sea surface is obtained by using the whitecap coverage (F) model proposed by Monahan [19], combined with the slope of the simulated sea surface. Specifically, assuming that the number of sea surface facets is M × N, we selected the 100F × M × N% facets with the largest sea surface slope and altered them with the wedge-like breaking waves as shown in Figure 2. The wedge-like breaking wave was constructed of rectangular plates with dimensions BC for the downwind side and DE for the upwind side (AB and CD represent the sea surface surrounding the breaking wave). Taking downwind observations as an example, there are six scattering path combinations from the radar to the wedge-like breaking waves, each with different induced currents: (i) single scattering from AB; (ii) double scattering from AB to BC; (iii) single scattering from BC; (iv) single scattering from CD; (v) double scattering from CD to DE; (vi) single scattering from DE. We still treat AB and DE as rough surfaces, and the single scattering from rough AB and DE is obtained using the capillary wave phase perturbation method mentioned in Section 2.2.1.

Wedge-Like Breaking Waves Scattering
In this paper, the deterministic distribution of breaking waves on the sea surfac obtained by using the whitecap coverage ( F ) model proposed by Monahan [19], co bined with the slope of the simulated sea surface. Specifically, assuming that the num of sea surface facets is facets with the largest surface slope and altered them with the wedge-like breaking waves as shown in Figur The wedge-like breaking wave was constructed of rectangular plates with dimensions for the downwind side and DE for the upwind side (AB and CD represent the sea surf surrounding the breaking wave). Taking downwind observations as an example, there six scattering path combinations from the radar to the wedge-like breaking waves, e with different induced currents: (i) single scattering from AB; (ii) double scattering fr AB to BC; (iii) single scattering from BC; (iv) single scattering from CD; (v) double scatt ing from CD to DE; (vi) single scattering from DE. We still treat AB and DE as rough s faces, and the single scattering from rough AB and DE is obtained using the capillary w phase perturbation method mentioned in Section 2.2.1. The multiple scattering of the breaking waves and the surrounding sea surface w calculated using multi-path scattering. The specific calculation formulas for multi-pat and 2 are given below: Path 1: The total induced electromagnetic current of the reflected field at BC w position r 2 is The multiple scattering of the breaking waves and the surrounding sea surface was calculated using multi-path scattering. The specific calculation formulas for multi-path 1 and 2 are given below: Path 1: The total induced electromagnetic current of the reflected field at BC with position r 2 is Path 2: The total induced electromagnetic current of the reflected field at DE with position r 4 is The far scattering field of the breaking waves can be obtained by substituting the induced electromagnetic current at BC and DE into the electromagnetic current radiation integral equation: For the sea surface with breaking waves under high sea conditions, the total scattering field can be written as where F is the whitecap coverage.

Modified Four-Path Model
The coupling scattering between the ship and sea surface can be represented as the coherent superposition of the scattering fields of four paths, as shown in Figure 3.
For the sea surface with breaking waves under high sea conditions, the total scattering field can be written as where F is the whitecap coverage.

Modified Four-Path Model
The coupling scattering between the ship and sea surface can be represented as the coherent superposition of the scattering fields of four paths, as shown in Figure 3. path 1. single scattering of the ship; path 2. forward scattering of the ship and then reflected by the sea surface; path 3. the reflection field of the sea surface and then scattered by the ship; path 4. the field reflected from the sea surface and scattered by the ship and then reflected from the sea surface again.
According to the mirror principle, the latter three equivalent paths are as follows: path 2 is equivalent to I→II→IV, path 3 is equivalent to IV→II→III and path 4 is equivalent to V →II→III→IV. When radar detects the sea surface, the spatial distribution of reflected energy is related to the roughness of the sea surface. For slightly rough surfaces, the reflected energy is mainly concentrated in the mirror direction, while for high sea conditions, the diffuse path 1. single scattering of the ship; path 2. forward scattering of the ship and then reflected by the sea surface; path 3. the reflection field of the sea surface and then scattered by the ship; path 4. the field reflected from the sea surface and scattered by the ship and then reflected from the sea surface again.
According to the mirror principle, the latter three equivalent paths are as follows: path 2 is equivalent to I→II→IV, path 3 is equivalent to IV→II→III and path 4 is equivalent to V→II→III→IV.
When radar detects the sea surface, the spatial distribution of reflected energy is related to the roughness of the sea surface. For slightly rough surfaces, the reflected energy is mainly concentrated in the mirror direction, while for high sea conditions, the diffuse reflection component is dominant. Thus, the reflection coefficient needs to be modified. The modified complex reflection coefficient is where R v,h is the reflection coefficient of each sea surface facet, and ρ s is the sea surface reflection factor as follows: where τ = σ h cos θ i /λ and σ h is the standard deviation of the height of the rough surface. Then, the total scattering field can be obtained by summing up the scattering field of the sea surface and the four paths as follows: where ρ is the modified complex reflection coefficient. E sea is the scattering field of the sea surface, E target is the scattering filed of the target. E path2 , E path3 and E path4 are the scattering fields depicted in Figure 3. The incident and scattering vectors in Equation (18) are depicted in Figure 3.

Numerical Results
To verify the correctness of the model, we presented the backscatter coefficient of sea surfaces generated using the Monte Carlo method (with the size of 256 × 256 m 2 and a windspeed of 5 m/s, averaging 50 sea surface samples) along with a comparison to SASS-II measured data [20], as shown in Figure 4. Additionally, we presented the backscatter coefficient of the sea surface with breaking waves (with the size of 256 × 256 m 2 and a windspeed of 13 m/s, averaging 50 sea surface samples) along with the comparison to measured data [21], as shown in Figure 5. We observe that even after employing Monte Carlo methods for statistical averaging under high sea conditions, the scattering curve is still not entirely smooth due to the presence of breaking wave structures. In future work, we will explore the effects of different breaking wave structures on the scattering characteristics. coefficient of the sea surface with breaking waves (with the size of 256 × 256 m and a windspeed of 13 m/s, averaging 50 sea surface samples) along with the comparison to measured data [21], as shown in Figure 5. We observe that even after employing Monte Carlo methods for statistical averaging under high sea conditions, the scattering curve is still not entirely smooth due to the presence of breaking wave structures. In future work, we will explore the effects of different breaking wave structures on the scattering characteristics.   As shown in Figure 6, the ship model has a length of 154 m, a width of 16 m and a draft of 6 m. Figure 6 shows the backscattering RCS of the ship model under different polarizations using the conducting and dielectric equivalent edge electromagnetic current method. The incidence angle is −90º~ 90º, corresponding to the half-space incidence from the stern to the bow. The incidence frequency is 10 GHz (X-band), the dielectric parameter of the dielectric ship is µ ε + + =4 1.5i, =2 i r r and the coating thickness is λ 0.04 . It can be seen from Figure 6 that the ship scattering peaks appear at horizontal incidence to the bow and stern and perpendicular incidence to the ship deck. There are also scattering peaks at several angles due to the complex structure of the tower, and the difference between HH and VV polarization is not obvious. The backscattering RCS of the dielectric ship decreases due to the existence of the coated dielectric layer, which shows that this material can achieve stealth in a certain range of angles at X-band for radar detection.  As shown in Figure 6, the ship model has a length of 154 m, a width of 16 m and a draft of 6 m. Figure 6 shows the backscattering RCS of the ship model under different polarizations using the conducting and dielectric equivalent edge electromagnetic current method. The incidence angle is −90 •~9 0 • , corresponding to the half-space incidence from the stern to the bow. The incidence frequency is 10 GHz (X-band), the dielectric parameter of the dielectric ship is µ r = 4 + 1.5i, ε r = 2 + i and the coating thickness is 0.04λ. It can be seen from Figure 6 that the ship scattering peaks appear at horizontal incidence to the bow and stern and perpendicular incidence to the ship deck. There are also scattering peaks at several angles due to the complex structure of the tower, and the difference between HH and VV polarization is not obvious. The backscattering RCS of the dielectric ship decreases due to the existence of the coated dielectric layer, which shows that this material can achieve stealth in a certain range of angles at X-band for radar detection. Figure 7 shows the composite scattering results of the ship and sea surface calculated with the equivalent edge electromagnetic current method combined with the sea surface scattering theory mentioned in Section 2. The simulated sea surface area is 512 m × 512 m, with a 1 m × 1 m interval. The sea surface wind speed is 5 m/s, and the wind direction is 180. It can be seen from Figure 7 that the total backscattering RCS increases with the consideration of the sea surface, compared with that in Figure 6. For near-vertical incidence, the sea surface scattering is dominant, and the target is completely submerged in sea background. For small grazing angle incidences, ship scattering is dominant, and the tower scattering peak (±75 incidence) is visible. seen from Figure 6 that the ship scattering peaks appear at horizontal incidence to the bow and stern and perpendicular incidence to the ship deck. There are also scattering peaks at several angles due to the complex structure of the tower, and the difference between HH and VV polarization is not obvious. The backscattering RCS of the dielectric ship decreases due to the existence of the coated dielectric layer, which shows that this material can achieve stealth in a certain range of angles at X-band for radar detection.   Figure 7 that the total backscattering RCS increases with the consideration of the sea surface, compared with that in Figure 6. For near-vertical incidence, the sea surface scattering is dominant, and the target is completely submerged in sea background. For small grazing angle incidences, ship scattering is dominant, and the tower scattering peak (±75 incidence) is visible.  Figure 8 shows the composite scattering results of the ship and sea surface with a wind speed of 10 m/s. It can be seen that when the wind speed increases, the sea surface backscatter RCS increases. The sea surface scattering still dominates for near-vertical incidence. At small grazing angles (incidence angle less than −60), the backscattering RCS of the sea surface (especially for HH polarization) is enhanced due to the breaking wave scattering, but the tower scattering peak (−75 incidence at the stern) is still visible.  Figure 9 shows the composite scattering results of the ship and sea surface with a wind speed of 15 m/s. It can be seen that with the increase in the wind speed, the sea surface scattering is enhanced. The backscattering RCS of the sea surface with small grazing angle incidence enhances especially for HH polarization due to the impact of the scattering of the breaking waves. However, the scattering peak (−75 incidence at the stern) of  Figure 8 shows the composite scattering results of the ship and sea surface with a wind speed of 10 m/s. It can be seen that when the wind speed increases, the sea surface backscatter RCS increases. The sea surface scattering still dominates for near-vertical incidence. At small grazing angles (incidence angle less than −60), the backscattering RCS of the sea surface (especially for HH polarization) is enhanced due to the breaking wave scattering, but the tower scattering peak (−75 incidence at the stern) is still visible.  Figure 8 shows the composite scattering results of the ship and sea surface with a wind speed of 10 m/s. It can be seen that when the wind speed increases, the sea surface backscatter RCS increases. The sea surface scattering still dominates for near-vertical incidence. At small grazing angles (incidence angle less than −60), the backscattering RCS of the sea surface (especially for HH polarization) is enhanced due to the breaking wave scattering, but the tower scattering peak (−75 incidence at the stern) is still visible.  Figure 9 shows the composite scattering results of the ship and sea surface with a wind speed of 15 m/s. It can be seen that with the increase in the wind speed, the sea surface scattering is enhanced. The backscattering RCS of the sea surface with small grazing angle incidence enhances especially for HH polarization due to the impact of the scattering of the breaking waves. However, the scattering peak (−75 incidence at the stern) of the tower is still visible.  Figure 9 shows the composite scattering results of the ship and sea surface with a wind speed of 15 m/s. It can be seen that with the increase in the wind speed, the sea surface scattering is enhanced. The backscattering RCS of the sea surface with small grazing angle incidence enhances especially for HH polarization due to the impact of the scattering of the breaking waves. However, the scattering peak (−75 incidence at the stern) of the tower is still visible.  Figures 10 and 11 show the composite backscattering of a dielectric ship on sea surfaces with different wind speeds. It can be seen that the composite backscattering RCS of the target is reduced compared with the conducting target due to the coated dielectric layer. With the increase in wind speed, the scattering enhancement of the breaking waves with small grazing angle makes the tower scattering peak submerged in the sea background. Figure 12 shows the azimuthal variation of the scattering of a dielectric ship and sea surfaces at different wind speeds with 60° incidence. It is noticeable that for ship-only scattering, four peaks emerge at the azimuth angle of 90° (bow), 180° (hull side), 270° (stern) and 360° or 0° (hull side), which corresponding to the four sides of the ship, respectively. When the scattering of the sea surface is taken into consideration, the RCS value significantly increases and the sea surface scattering becomes dominant, resulting in the disappearance of peaks at 90° and 270°. As the wind speed increases to 10 m/s, the sea surface scattering increases and the scattering peaks of the target are completely submerged.   Figures 10 and 11 show the composite backscattering of a dielectric ship on sea surfaces with different wind speeds. It can be seen that the composite backscattering RCS of the target is reduced compared with the conducting target due to the coated dielectric layer. With the increase in wind speed, the scattering enhancement of the breaking waves with small grazing angle makes the tower scattering peak submerged in the sea background. Figure 12 shows the azimuthal variation of the scattering of a dielectric ship and sea surfaces at different wind speeds with 60° incidence. It is noticeable that for ship-only scattering, four peaks emerge at the azimuth angle of 90° (bow), 180° (hull side), 270° (stern) and 360° or 0° (hull side), which corresponding to the four sides of the ship, respectively. When the scattering of the sea surface is taken into consideration, the RCS value significantly increases and the sea surface scattering becomes dominant, resulting in the disappearance of peaks at 90° and 270°. As the wind speed increases to 10 m/s, the sea surface scattering increases and the scattering peaks of the target are completely submerged.

Discussion
This study aimed to provide a composite scattering model of the sea surface with a conducting and dielectric ship above it under different sea conditions. According to the simulated results in Section 3, the following statements can be made.
(1) The backscattering RCS of the ship decreases due to the existence of the coated dielectric layer, indicating that radar stealth can be achieved by using coated absorbing materials to reduce the signal intensity of the target. (2) The sea surface scattering dominates for near-vertical incidence, while for small grazing angle incidences, ship scattering is dominant, and the tower scattering peak is visible.

Discussion
This study aimed to provide a composite scattering model of the sea surface with a conducting and dielectric ship above it under different sea conditions. According to the simulated results in Section 3, the following statements can be made.
(1) The backscattering RCS of the ship decreases due to the existence of the coated dielectric layer, indicating that radar stealth can be achieved by using coated absorbing materials to reduce the signal intensity of the target. (2) The sea surface scattering dominates for near-vertical incidence, while for small grazing angle incidences, ship scattering is dominant, and the tower scattering peak is visible.

Discussion
This study aimed to provide a composite scattering model of the sea surface with a conducting and dielectric ship above it under different sea conditions. According to the simulated results in Section 3, the following statements can be made.
(1) The backscattering RCS of the ship decreases due to the existence of the coated dielectric layer, indicating that radar stealth can be achieved by using coated absorbing materials to reduce the signal intensity of the target. (2) The sea surface scattering dominates for near-vertical incidence, while for small grazing angle incidences, ship scattering is dominant, and the tower scattering peak is visible.
(3) Under high sea conditions, it is crucial to consider the impact of the breaking waves on the scattered echo. The mirror scattering decreases, but the incoherent scattering increases with the increase in wind speed due to the roughness of the sea surface. The backscattering RCS of the sea surface (especially for HH polarization) is enhanced due to breaking wave scattering at small grazing angles. As wind speed increases, the scattering enhancement of breaking waves with small grazing angles causes the tower scattering peak to submerge in the sea background. (4) The azimuthal variation of the ship scattering shows several peaks at the angle that is perpendicular to the four sides of the ship. When the scattering of the sea surfaces is taken into consideration, the sea surface scattering results in the disappearance of scattering peaks at the bow and stern. When considering the scattering of the sea surface, the scattering peaks at the bow and stern disappear. Additionally, for high sea conditions, the scattering peaks of the target are completely submerged.
To sum up, it is necessary to consider the impact of the breaking waves and the dielectric materiel of the target on the scattered echo. We believe that the model can provide a theoretical basis for radar to detect marine targets with better efficiency.

Conclusions
In this paper, a composite scattering model from the sea surface with breaking waves and a conducting and dielectric ship above it is presented. The EEC theory is used to calculate the ship scattering, while the scattering of the sea surface with wedge-like breaking waves is calculated using the capillary wave phase perturbation method and the multi-path scattering model. The coupling scattering between them is obtained using the modified fourpath model. Using the proposed model, the RCS results of the conducting and dielectric ships on sea surfaces at low, medium and high wind speeds under different polarizations are simulated numerically. The numerical results show that radar stealth can be achieved in a certain range by using coated absorbing materials to reduce the signal intensity of the target. When the wind speed increases, the composite backscattering from the sea surface and target for both polarizations increases with the consideration of the effect of breaking waves, especially for upwind under high sea conditions and small grazing angle incidence, and the scattering enhancement of HH polarization is more obvious. The tower scattering peak of the conductor ship is visible at low and high wind speeds. However, it is difficult to separate the target from the sea background under high sea conditions due to the scattering weakness of the dielectric target itself, the enhancement of the scattering effect of the sea surface and the spike scattering of the breaking waves.
Author Contributions: Conceptualization, Z.W.; formal analysis, X.Z.; funding acquisition, X.Z.; investigation, X.Z.; methodology, X.Z. and X.S.; writing-original draft, X.Z.; writing-review and editing, X.Z. and X.S. All authors have read and agreed to the published version of the manuscript.