# Reliable Approach for EM Scattering Calculation from Sea Surface Covered with Foams Based on MFBM and VRT Theory

Authors: Tao Song; Chuang-Ming Tong; Li-Li Cong; Peng Peng
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

A reliable approach based on modified facet-based model and Vector Radiative Transfer theory is presented to calculate electromagnetic scattering from a particular electrically large sea surface superimposed with foams, which can handily give both monostatic and bistatic scattering results and could be applied to synthetic aperture radar imagery simulation. The facet model is derived from Fuks’ first-order small perturbation method function, and then the Kirchhoff model is introduced to revise the results in view of the dependency on cut-off wave number at near vertical incidence angles. Additionally, the contributions of foams are taken into consideration on the basis of Vector Radiative Transfer theory. The accuracy and superiority of this proposed approach are demonstrated in comparison with traditional facet model, which illustrates that the results of this approach agree better with experimental results. Moreover, several examples are given to verify that the proposed approach is of more significance at large incidence angles and high wind speed.


---

## Body

## 1. Introduction

Investigation on electromagnetic (EM) scattering from maritime scene is of great value in both civil and military applications, such as maritime environment monitoring and clutter rejection in target detection. Extensive endeavors have been devoted to evaluate the scattering characteristic in either theoretical or experimental terms. Usually, approaches for EM scattering calculation from rough sea surface in theoretical term can be classified as numerical category and analytical category. Among the numerical efforts, some recent attempts could be noted, for instance, the Multiple Sweep Method of Moments [1], the Accelerated Forward-Backward Method [2], and the Ordered Multiple Method [3]. Although these methods prevail over traditional numerical ones in terms of accelerated skills, they would still be too tricky and time-consuming for general applications, particularly for electrically large maritime scene. On the other hand, analytical methods have received more and more attention in view of their high efficiency, relative accuracy, and clear physical interpretation. Crombie firstly discovered the Bragg resonant phenomena that the Bragg waves resonant with EM waves dominate in the echo [4]. And then two-scale theory has been always a famous approach since it was presented in the 1970s [5]. Further, the modified two-scale method was presented and validated with better accuracy [6]. However, it is still a statistic approach, which obtains an average of the scattering coefficients under different conditions without a particular sea height map, and nothing is said about local information.In past decades, synthetic aperture radar (SAR) has been greatly developed in consideration of its independency on daylight and weather, and SAR imaging has become one of the most important means for remote sensing of maritime environment. At the same time, local information is useful and desirable in SAR imagery simulation of marine scene, especially the one with oil film and ship wake. This demand has encouraged the development of facet-based approaches, which try to break the surface into facets and then the scattering contributions from individual facets could be obtained. A facet backscattering model was presented by Franceschetti et al. by means of Kirchhoff approximation [7], but it loses sight of the Bragg phenomena. Chen et al. came up with a slope-determined facet model that could reflect both the specular and diffuse configurations in monostatic or bistatic system [8].In reality, with the increase of wind speed above sea surface, foams would come into being. What is more, one-third of sea surface would be covered with foams when the velocity of wind comes up to 25 m/s. So it is of great significance to take the influences of foams on sea surface scattering into account at high wind speed condition. Jin [9] endeavored to investigate the scattering characteristics of sea surface with foams by introducing Vector Radiative Transfer (VRT) theory into two-scale model (TSM). Liang and Guo [10] combined VRT theory with modified TSM. However, the two aforementioned statistical methods cannot be applied to SAR imagery simulation because of the deficiency of local information.In this work, a facet-based model derived from Fuks’ first-order small perturbation method function [11, 12] is proposed to calculate the scattering contribution from individual sea surface facets in both monostatic and bistatic configurations. Moreover, the Kirchhoff model is introduced to revise the results because the small perturbation method is not suitable and the results are dependent on cut-off wave number at near vertical incidence angles. Meanwhile, the scattering effects of foams are considered and calculated according to the VRT theory, which is of significance at high wind speed and grazing incidence consideration. The proposed approach would be attractive on the ground that the formulation is tractable and time-saving, and the mechanism is more comprehensive including both specular and diffuse features. What is more, as this scheme can calculate the local scattering contributions, it could be readily applied to SAR imagery simulation and the analysis of scattering contributions from complex sea surface with ship wake and oil film.


## 2. Model Description and Formulation

### 2.1. Modified Facet-Based Model

Sea surface can be seen as a two-scale model where small-scale waves are superimposed on the large-scale waves. The ELH spectrum [13] can involve both long and short wind-driven waves, so it is used to generate sea surface here. Then the sea surface is decomposed into facets with different slopes and the total scattering can be calculated by summing up the contributions from all the facets. Therefore, an adequate formula for scattering from arbitrary slope facet needs to be derived.As shown in Figure1, it is assumed that the mean of the microscopic rough surface is ζ=0 and the distribution satisfies statistically spatial homogeneity. When considering a unit incident plane wave Ei in xoz plane, we can get the scattering coefficient according to the first-order small perturbation method function presented by Fuks [11, 12]:(1)σpq0k^i,k^s=πk4ε-12Fpq2Sζql,where k^i is a unit vector directed from the transmitter, while k^s is a unit vector directed to the reception point; ql is the projection vector ofq on the plane z=0; q=k^s-k^i;k denotes the electromagnetic wave number; ε indicates the relative permittivity of the dielectric surface; p, q=h, v stand for the polarization direction of scattering and incident wave vector, respectively; Fpq represents the polarization factor and detailed expressions in Appendix; Sζql is the spatial power spectrum of microscopic rough surface.Figure 1
Geometry of microscopic rough scattering surface.Next, we need to calculate the scattering coefficient of arbitrary slope facet.As shown in Figure2, x^g,y^g,z^g is the global coordinate, while the local coordinate is created at the pointo′ of facet as follows:(2)z^l=n^,y^l=n^×k^in^×k^i,x^l=y^l×z^l,where n^=-Zxx^g-Zyy^g-z^g/1+Zx2+Zy2 indicates the normal vector of the facet; θi,θs,∅i,∅s and θil,θsl,∅il,∅sl are the global and local angles of incidence and scattering direction in different coordinates.Figure 2
Geometry of global coordinate and local coordinate.Referring to (1), we can write the scattering coefficient of arbitrary slope rough facet:(3)σPQfacetk^i,k^s=πk4ε-12F~PQ2Sζql,where F~PQ is defined as the polarization factor in global coordinate and can be obtained from the following relation:(4)F~VVF~VHF~HVF~HH=V^s·v^sH^s·v^sV^s·h^sH^s·h^sFvvFvhFhvFhhV^i·v^iV^i·h^iH^i·v^iH^i·h^i,where H^i,V^i,H^s,V^s and h^i,v^i,h^s,v^s are the unit polarization vectors in global and local coordinate, respectively.According to the Bragg resonance hypothesis, both the Bragg waves spreading along and against the radar sight direction contribute to the radar receiver; thus (3) can be rewritten as(5)σPQfacetk^i,k^s=πk4ε-12F~PQ2Scapiql,where Scapiql=1/2SEcapi-ql+SEcapiql is the Bragg waves component; SEcapi is the small-scale capillary ripple part of one-sided ELH spectrum.Thus, cut-off wave numberkcut is introduced to divide ELH spectrum into large-scale gravity wave part SEgrav and small-scale capillary ripple part SEcapi:(6)SEgravk=0,k≥π6SEk,k≤π6,SEcapik=SEk,k≥π60,k≤π6,where SEk is the two-dimensional ELH spectrum and the kcut here is equal to π/6.Assume that the lengths of the two-dimensional simulation sea surface sample inx and y directions are Lx and Ly, respectively. The area is A=LxLy, the numbers of discrete points are M and N, and the distances between two adjacent points are Δx and Δy. Then, according to (5), we can get the total scattering coefficient by summing up the contributions from all of the facets:(7)σPQtotalk^i,k^s=1A∑m=1M∑n=1NσPQ,mnTSPMk^i,k^sΔxΔy,where σPQ,mnTSPMk^i,k^s is the normalized Bragg scattering coefficient of the number mn facet.Unfortunately, there is no unified standard to choose the advisable cut-off wave number. As shown in Figure3, the small perturbation method is not suitable and the results are dependent on kcut at near vertical incident angles.Figure 3
Impact of cut-off wave number on scattering coefficients.The parameters in the simulation are in order as follows: the incident wave is at Ku-band (14.0 GHz), the wind speed at 10 m above the mean sea level is selected byu10=5.0 m/s, the angle of wind direction ∅w is defined as the included angle of wind direction and x^g axis, here ∅w=0°, the numbers of discrete points are M=N=256, the distances between two adjacent points are Δx=Δy=1.0 m, and the relative permittivity calculated by Klein model [14] at 20°C and 35‰ salinity is ε=55.845,37.731.In view of the dependency on cut-off wave number at near vertical incident angles, Kirchhoff model can be introduced to revise the slope-modulated perturbation coefficient [15]. Then (7) can be rewritten as(8)σPQtotalk^i,k^s=1A∑m=1M∑n=1NσPQ,mnKAMk^i,k^s+σPQ,mnTSPMk^i,k^sΔxΔy,where σPQ,mnKAMk^i,k^s is the normalized specular scattering coefficient of the number mn facet and can be obtained by(9)σPQ,mnKAMk^i,k^s=πk2q2qz4F~PQKAMPzxtan,zytan,where P• is the slope probability density function (refer to [16]). The detailed expressions of polarization factor F~PQKAM are given in Appendix.By now, we can calculate the total scattering coefficient of a particular sea surface sample by (8). At the same time, the local contribution from any facet can be readily obtained. So this modified facet-based scattering model could be applied to SAR imagery simulation.


### 2.2. Foams Scattering Model

In reality, with the increase of wind speed above the sea surface, foams come into being. So it is inaccurate to calculate the scattering from sea surface with forenamed modified facet-based model and the influences of foams on sea surface scattering should be taken into account at high-wind-speed condition.As shown in Figure4, the Vector Radiative Transfer (VRT) functions of the foam layer can be expressed by [9](10)cos⁡θsddzI¯θs,φs,z=-Ke¯¯θs,φs·I¯θs,φs,z+∫02πdφi∫0π/2dθi×sin⁡θs·P¯¯θs,φs;θi,φiI¯θi,φi,z+Q¯z,θs,φs,-cos⁡θsddzI¯π-θs,φs,z=-Ke¯¯θs,φs·I¯π-θs,φs,z+∫02πdφi∫0π/2dθisin⁡θs·P¯¯π-θs,φs;θi,φi·I¯θi,φi,z+∫02πdφi∫0π/2dθisin⁡θs·P¯¯π-θs,φs;π-θi,φi·I¯π-θi,φi,z+Q¯z,θs,φs,where θi, ∅i  θs, and ∅s represent incident angle, incident azimuth angle, scattering angle, and scattering azimuth angle, respectively. I-θ,φ,z is a 4×1 intensity vector. Ke̿ stands for the extinction matrix. Q- indicates the heat emitter undersea which is equal to zero when illuminated by EM wave above the sea surface. P̿ denotes the phase matrix including the multiple scattering and coupling relationship among particles which can be obtained by Rayleigh approximation. The boundary condition is given in [10]. VRT functions could be solved utilizing boundary condition together with constant variation function, and then each order result can be obtained by iterative method.Figure 4
Geometry of foams scattering.The bistatic scattering coefficient of sea surface with foams can be defined as [17](11)σpqθs,φs;θi,φi=4πcos⁡θsIpθs,φsIqθi,φiwhere Iq and Ip are the incident and scattering intensities, respectively.We can further get the zeroth-order and first-order scattering coefficientsσpqa0, σpqb1, σpqc1, σpqd1, and σpqe1 of sea surface with foams:(12)σpqa0θs,ϕs,θi,ϕi=σpqe-ked1/cos⁡θs+1/cos⁡θiσpqb1θs,ϕs,θi,ϕi=4πkecos⁡θsP¯¯pqθs,ϕs,π-θi,ϕi1-e-ked1/cos⁡θs+1/cos⁡θi,σpqc1θs,ϕs,θi,ϕi=4πkecos⁡θs∫02πdθ′sin⁡θ′∫02πdϕ′∑l=v,hP¯¯plθs,ϕs,θ′,ϕ′·R¯¯pqθ′,ϕ′,π-θi,ϕi×cos⁡θ′cos⁡θs-cos⁡θ′e-ked1/cos⁡θ′-1/cos⁡θs,σpqd1θs,ϕs,θi,ϕi=4πkecos⁡θs∫0π/2dθ′sin⁡θ′∫02πdϕ′∑m=v,hR¯¯pmθs,ϕs,π-θ′,ϕ′×P¯¯mqπ-θ′,ϕ′,π-θi,ϕi×cos⁡θicos⁡θi-cos⁡θ′e-ked1/cos⁡θ′-1/cos⁡θi,σpqe1θs,ϕs,θi,ϕi=4πkecos⁡θs∫0π/2dθ′sin⁡θ′∫02πdϕ′∑l=v,hR¯¯plθs,ϕs,π-θ′,ϕ′×∫0π/2dθ′′sin⁡θ′′∫02πdϕ′′×∑m=v,hP¯¯lmπ-θ′,ϕ′,θ′′,ϕ′′R¯¯mqθ′′,ϕ′′,θi,ϕi, where ke is the extinction coefficient and σpq is the scattering coefficient calculated by the presented modified facet-based model above. The physical scattering processes of different coefficients are demonstrated in Figure 5.Figure 5
Scattering processes of different coefficients.In consequence, the scattering coefficient of the area covered with foams can be calculated by(13)σpqfoamθs,φs;θi,φi=σpqa0θs,φs;θi,φi+σpqb1θs,φs;θi,φi+σpqc1θs,φs;θi,φi+σpqd1θs,φs;θi,φi+σpqe1θs,φs;θi,φi.In fact, it could be assumed that the foams are made up of many water-air particles which could be represented by the model in Figure6. Referring to [18], a=250μm is the outer radius, b is the inner radius, and a-b=20 μm is the thickness of water film. The thickness of foams here is presumed to be 2 cm. Then the extinction coefficient ke of this foams model can be calculated by [19](14)ke=N02πk2∑n=1∞2n+1Re-TnN-TnM,where N0 is the number of foams in unit volume and could be obtained according to [20].Figure 6
Geometry of foams.The coverage percentage of foams in wind-driven sea surface is intimately related to wind speed and seawater temperature. While the temperature difference between air and seawater is moderate, it can be calculated by(15)C=2.32×10-6U103.4988.Then the number of facets covered with foams isNf=M×N×C. According to the slope criterion presented by Longuet-Higgins and Fox [21], we can add foams to Nf facets with greater slopes. Thus, we can get the distribution of foams instead of the coverage percentage alone, which is of great significance to SAR imagery simulation of complex sea surface.As a result, we can calculate the total scattering coefficient of a particular wind-driven sea surface covered with foams by(16)σpqtotal=∑m=1MN-Nfσpq,mKAM+σpq,mTSPMΔxΔy+∑n=1Nfσpq,nfoamΔxΔyA.


## 2.1. Modified Facet-Based Model

Sea surface can be seen as a two-scale model where small-scale waves are superimposed on the large-scale waves. The ELH spectrum [13] can involve both long and short wind-driven waves, so it is used to generate sea surface here. Then the sea surface is decomposed into facets with different slopes and the total scattering can be calculated by summing up the contributions from all the facets. Therefore, an adequate formula for scattering from arbitrary slope facet needs to be derived.As shown in Figure1, it is assumed that the mean of the microscopic rough surface is ζ=0 and the distribution satisfies statistically spatial homogeneity. When considering a unit incident plane wave Ei in xoz plane, we can get the scattering coefficient according to the first-order small perturbation method function presented by Fuks [11, 12]:(1)σpq0k^i,k^s=πk4ε-12Fpq2Sζql,where k^i is a unit vector directed from the transmitter, while k^s is a unit vector directed to the reception point; ql is the projection vector ofq on the plane z=0; q=k^s-k^i;k denotes the electromagnetic wave number; ε indicates the relative permittivity of the dielectric surface; p, q=h, v stand for the polarization direction of scattering and incident wave vector, respectively; Fpq represents the polarization factor and detailed expressions in Appendix; Sζql is the spatial power spectrum of microscopic rough surface.Figure 1
Geometry of microscopic rough scattering surface.Next, we need to calculate the scattering coefficient of arbitrary slope facet.As shown in Figure2, x^g,y^g,z^g is the global coordinate, while the local coordinate is created at the pointo′ of facet as follows:(2)z^l=n^,y^l=n^×k^in^×k^i,x^l=y^l×z^l,where n^=-Zxx^g-Zyy^g-z^g/1+Zx2+Zy2 indicates the normal vector of the facet; θi,θs,∅i,∅s and θil,θsl,∅il,∅sl are the global and local angles of incidence and scattering direction in different coordinates.Figure 2
Geometry of global coordinate and local coordinate.Referring to (1), we can write the scattering coefficient of arbitrary slope rough facet:(3)σPQfacetk^i,k^s=πk4ε-12F~PQ2Sζql,where F~PQ is defined as the polarization factor in global coordinate and can be obtained from the following relation:(4)F~VVF~VHF~HVF~HH=V^s·v^sH^s·v^sV^s·h^sH^s·h^sFvvFvhFhvFhhV^i·v^iV^i·h^iH^i·v^iH^i·h^i,where H^i,V^i,H^s,V^s and h^i,v^i,h^s,v^s are the unit polarization vectors in global and local coordinate, respectively.According to the Bragg resonance hypothesis, both the Bragg waves spreading along and against the radar sight direction contribute to the radar receiver; thus (3) can be rewritten as(5)σPQfacetk^i,k^s=πk4ε-12F~PQ2Scapiql,where Scapiql=1/2SEcapi-ql+SEcapiql is the Bragg waves component; SEcapi is the small-scale capillary ripple part of one-sided ELH spectrum.Thus, cut-off wave numberkcut is introduced to divide ELH spectrum into large-scale gravity wave part SEgrav and small-scale capillary ripple part SEcapi:(6)SEgravk=0,k≥π6SEk,k≤π6,SEcapik=SEk,k≥π60,k≤π6,where SEk is the two-dimensional ELH spectrum and the kcut here is equal to π/6.Assume that the lengths of the two-dimensional simulation sea surface sample inx and y directions are Lx and Ly, respectively. The area is A=LxLy, the numbers of discrete points are M and N, and the distances between two adjacent points are Δx and Δy. Then, according to (5), we can get the total scattering coefficient by summing up the contributions from all of the facets:(7)σPQtotalk^i,k^s=1A∑m=1M∑n=1NσPQ,mnTSPMk^i,k^sΔxΔy,where σPQ,mnTSPMk^i,k^s is the normalized Bragg scattering coefficient of the number mn facet.Unfortunately, there is no unified standard to choose the advisable cut-off wave number. As shown in Figure3, the small perturbation method is not suitable and the results are dependent on kcut at near vertical incident angles.Figure 3
Impact of cut-off wave number on scattering coefficients.The parameters in the simulation are in order as follows: the incident wave is at Ku-band (14.0 GHz), the wind speed at 10 m above the mean sea level is selected byu10=5.0 m/s, the angle of wind direction ∅w is defined as the included angle of wind direction and x^g axis, here ∅w=0°, the numbers of discrete points are M=N=256, the distances between two adjacent points are Δx=Δy=1.0 m, and the relative permittivity calculated by Klein model [14] at 20°C and 35‰ salinity is ε=55.845,37.731.In view of the dependency on cut-off wave number at near vertical incident angles, Kirchhoff model can be introduced to revise the slope-modulated perturbation coefficient [15]. Then (7) can be rewritten as(8)σPQtotalk^i,k^s=1A∑m=1M∑n=1NσPQ,mnKAMk^i,k^s+σPQ,mnTSPMk^i,k^sΔxΔy,where σPQ,mnKAMk^i,k^s is the normalized specular scattering coefficient of the number mn facet and can be obtained by(9)σPQ,mnKAMk^i,k^s=πk2q2qz4F~PQKAMPzxtan,zytan,where P• is the slope probability density function (refer to [16]). The detailed expressions of polarization factor F~PQKAM are given in Appendix.By now, we can calculate the total scattering coefficient of a particular sea surface sample by (8). At the same time, the local contribution from any facet can be readily obtained. So this modified facet-based scattering model could be applied to SAR imagery simulation.


## 2.2. Foams Scattering Model

In reality, with the increase of wind speed above the sea surface, foams come into being. So it is inaccurate to calculate the scattering from sea surface with forenamed modified facet-based model and the influences of foams on sea surface scattering should be taken into account at high-wind-speed condition.As shown in Figure4, the Vector Radiative Transfer (VRT) functions of the foam layer can be expressed by [9](10)cos⁡θsddzI¯θs,φs,z=-Ke¯¯θs,φs·I¯θs,φs,z+∫02πdφi∫0π/2dθi×sin⁡θs·P¯¯θs,φs;θi,φiI¯θi,φi,z+Q¯z,θs,φs,-cos⁡θsddzI¯π-θs,φs,z=-Ke¯¯θs,φs·I¯π-θs,φs,z+∫02πdφi∫0π/2dθisin⁡θs·P¯¯π-θs,φs;θi,φi·I¯θi,φi,z+∫02πdφi∫0π/2dθisin⁡θs·P¯¯π-θs,φs;π-θi,φi·I¯π-θi,φi,z+Q¯z,θs,φs,where θi, ∅i  θs, and ∅s represent incident angle, incident azimuth angle, scattering angle, and scattering azimuth angle, respectively. I-θ,φ,z is a 4×1 intensity vector. Ke̿ stands for the extinction matrix. Q- indicates the heat emitter undersea which is equal to zero when illuminated by EM wave above the sea surface. P̿ denotes the phase matrix including the multiple scattering and coupling relationship among particles which can be obtained by Rayleigh approximation. The boundary condition is given in [10]. VRT functions could be solved utilizing boundary condition together with constant variation function, and then each order result can be obtained by iterative method.Figure 4
Geometry of foams scattering.The bistatic scattering coefficient of sea surface with foams can be defined as [17](11)σpqθs,φs;θi,φi=4πcos⁡θsIpθs,φsIqθi,φiwhere Iq and Ip are the incident and scattering intensities, respectively.We can further get the zeroth-order and first-order scattering coefficientsσpqa0, σpqb1, σpqc1, σpqd1, and σpqe1 of sea surface with foams:(12)σpqa0θs,ϕs,θi,ϕi=σpqe-ked1/cos⁡θs+1/cos⁡θiσpqb1θs,ϕs,θi,ϕi=4πkecos⁡θsP¯¯pqθs,ϕs,π-θi,ϕi1-e-ked1/cos⁡θs+1/cos⁡θi,σpqc1θs,ϕs,θi,ϕi=4πkecos⁡θs∫02πdθ′sin⁡θ′∫02πdϕ′∑l=v,hP¯¯plθs,ϕs,θ′,ϕ′·R¯¯pqθ′,ϕ′,π-θi,ϕi×cos⁡θ′cos⁡θs-cos⁡θ′e-ked1/cos⁡θ′-1/cos⁡θs,σpqd1θs,ϕs,θi,ϕi=4πkecos⁡θs∫0π/2dθ′sin⁡θ′∫02πdϕ′∑m=v,hR¯¯pmθs,ϕs,π-θ′,ϕ′×P¯¯mqπ-θ′,ϕ′,π-θi,ϕi×cos⁡θicos⁡θi-cos⁡θ′e-ked1/cos⁡θ′-1/cos⁡θi,σpqe1θs,ϕs,θi,ϕi=4πkecos⁡θs∫0π/2dθ′sin⁡θ′∫02πdϕ′∑l=v,hR¯¯plθs,ϕs,π-θ′,ϕ′×∫0π/2dθ′′sin⁡θ′′∫02πdϕ′′×∑m=v,hP¯¯lmπ-θ′,ϕ′,θ′′,ϕ′′R¯¯mqθ′′,ϕ′′,θi,ϕi, where ke is the extinction coefficient and σpq is the scattering coefficient calculated by the presented modified facet-based model above. The physical scattering processes of different coefficients are demonstrated in Figure 5.Figure 5
Scattering processes of different coefficients.In consequence, the scattering coefficient of the area covered with foams can be calculated by(13)σpqfoamθs,φs;θi,φi=σpqa0θs,φs;θi,φi+σpqb1θs,φs;θi,φi+σpqc1θs,φs;θi,φi+σpqd1θs,φs;θi,φi+σpqe1θs,φs;θi,φi.In fact, it could be assumed that the foams are made up of many water-air particles which could be represented by the model in Figure6. Referring to [18], a=250μm is the outer radius, b is the inner radius, and a-b=20 μm is the thickness of water film. The thickness of foams here is presumed to be 2 cm. Then the extinction coefficient ke of this foams model can be calculated by [19](14)ke=N02πk2∑n=1∞2n+1Re-TnN-TnM,where N0 is the number of foams in unit volume and could be obtained according to [20].Figure 6
Geometry of foams.The coverage percentage of foams in wind-driven sea surface is intimately related to wind speed and seawater temperature. While the temperature difference between air and seawater is moderate, it can be calculated by(15)C=2.32×10-6U103.4988.Then the number of facets covered with foams isNf=M×N×C. According to the slope criterion presented by Longuet-Higgins and Fox [21], we can add foams to Nf facets with greater slopes. Thus, we can get the distribution of foams instead of the coverage percentage alone, which is of great significance to SAR imagery simulation of complex sea surface.As a result, we can calculate the total scattering coefficient of a particular wind-driven sea surface covered with foams by(16)σpqtotal=∑m=1MN-Nfσpq,mKAM+σpq,mTSPMΔxΔy+∑n=1Nfσpq,nfoamΔxΔyA.


## 3. Model Validation and Analysis

The accuracy and validity of the scheme proposed in this paper are verified by comparing the backward-scattering results of the presented approach with those of traditional facet model and experimental results in [22]. In simulation, the incident frequency is at 13.9 GHz, the temperature on sea surface is 17.4°C, the salinity is 32.54‰, ε=44.43,39.55, and the wind speed at 10 m height is 10 m/s. As shown in Figure 7, the results of teh proposed approach in horizontal polarization are generally higher than those of the traditional facet model without foams, which match better with the experimental results, especially for large incidence angles. The reason for the difference is that at large incident angles the multiple scattering matters and the SPM and KA theory is unserviceable. Therefore the traditional facet model based on SPM and KA theory is inaccurate to some extent. While the foams scattering mechanism proposed in this paper includes the multiple scattering contributions which is in accordance with the physical situation, in vertical polarization, the effects of foams at different incidence angles are faint.Figure 7
Comparison of backward-scattering results between and among the presented approach, traditional facet model, and experiment. (a) HH polarization. (b) VV polarization.
(a)
(b)For further validation, Figure8 depicts the backward-scattering results with and without foams at different wind speeds of 4 m/s, 10 m/s, and 15 m/s for vertical polarization. The incident frequency is at Ku-band. It is clearly shown that the backward-scattering results increase as the wind speed increases. This is because the sea surface becomes rougher and the incoherent scattering grows with the wind breezing up. Similarly, the influences of foams at small angles are unnoticeable and increase with the angles growing up. The differences between the backward-scattering results of sea surface with and without foams are minimum at the speed 4 m/s and maximum at the speed 10 m/s. The reason is that the coverage rate of foams increases with the increase of wind speed and foams contribute more to the backward-scattering results.Figure 8
Comparison of backward-scattering results at different wind speeds.In order to investigate the impact of azimuth angles on backward-scattering coefficients, the results at different azimuth angles are demonstrated in Figure9. The parameters are the same as Figure 3. It is clearly observed that the backward-scattering coefficients for vertical polarization are larger than those for horizontal polarization, while the results in horizontal polarization are more sensitive to wind speeds and azimuth angles. Furthermore, the results in downwind and backwind directions are larger than those in side-wind directions. This is because the sea surface is rougher in downwind and backwind directions than in side-wind directions, which is more noticeable at higher wind speed. In practical situation, the backward-scattering coefficients in backwind directions are always larger than those in downwind directions because the wave crest is tilted down the wind. All the same, the results are of no differences in simulations owing to the oversight of this asymmetry.Figure 9
Impact of azimuth angles on backward-scattering coefficients.


## 4. Conclusion

In this paper, a reliable approach based on MFBM and VRT theory is developed for the calculation of scattering from a particular electrically large sea surface covered with foams. This proposed scheme is of more comprehensive significance because the contributions from not only the specular and diffuse configurations but also the foams are taken into consideration. Moreover, it can handily give both the monostatic and bistatic scattering results and could overwhelm the application to SAR imagery simulation. The comparisons show better agreement with experimental results than traditional facet model, which demonstrates the accuracy and superiority of this scheme. Additionally, numerical results show that the contributions of foams are assignable and of great significance at large incidence angles and high wind speeds.



---

*Source: 2035381-2018-02-19.xml*