The Fast Simulation of Scattering Characteristics from a Simplified Time Varying Sea Surface
This paper aims at applying a simplified sea surface model into the physical optics (PO) method to accelerate the scattering calculation from 1D time varying sea surface. To reduce the number of the segments and make further improvement on the efficiency of PO method, a simplified sea surface is proposed. In this simplified sea surface, the geometry of long waves is locally approximated by tilted facets that are much longer than the electromagnetic wavelength. The capillary waves are considered to be sinusoidal line superimposing on the long waves. The wavenumber of the sinusoidal waves is supposed to satisfy the resonant condition of Bragg waves which is dominant in all the scattered short wave components. Since the capillary wave is periodical within one facet, an analytical integration of the PO term can be performed. The backscattering coefficient obtained from a simplified sea surface model agrees well with that obtained from a realistic sea surface. The Doppler shifts and width also agree well with the realistic model since the capillary waves are taken into consideration. The good agreements indicate that the simplified model is reasonable and valid in predicting both the scattering coefficients and the Doppler spectra.
The calculation of electromagnetic (EM) scattering from a time varying surface is important in many fields such as radar surveillance, target tracking, and ocean remote sensing . Useful techniques have already been developed to provide realistic results. They can be based on exact numerical methods (MoM, FEM, FDTD, and so on [2–5]) or approximate approaches . Because numerical methods are unfortunately not efficient, the approximate approaches are widely used for the moment to calculate the scattered field from a large time varying surface. Among them, PO [7, 8] is most employed because it is simple and easy to implement.
However, the PO method is still limited by the number of unknowns when dealing with large time varying sea surface. When dealing with the scattering problem form the sea, it should be noted that the sea surface should be divided into small segments whose length of each segment on the realistic sea surface should be 1/8~1/10 wavelength of the incident wave to accurately reflect the geometry characteristics of the sea. Since each segment is that small, the number of the segments form the realistic sea surface will be very large. To reduce the number of segments, one can divide the sea surface into larger segments. However, the larger the segment is, the more inaccuracy will be shown in the scattering coefficient and Doppler spectrum. The inaccurate results may be caused by two reasons. Firstly, the phase difference on one segment will be neglected in this case. This will lead to the inaccuracy of the integration on each segment. Secondly, the capillary waves superimposed on each segment are not taken into consideration. This will cause the fact that the Doppler value in some spectrum regions cannot be detected by the radar. Since the segment cannot be larger than 1/8~1/10 wavelength, the number of the unknowns will be tremendous especially when dealing with time varying large sea surface, which will limit the efficiency of PO method.
In order to get accurate results with larger segments, as well as promote the efficiency of conventional PO method, we add capillary waves . The capillary wave considered to be sinusoidal line superimposing on each segment. The wavenumbers of the sinusoidal waves are supposed to satisfy the resonant condition of Bragg waves which are predominant in the scattered short wave component. In this way, we substitute the simplified sea surface for the realistic sea surface approximately. Since the capillary wave is periodical within one facet, the analytical expression of the induced currents can be given after we get the currents on the first period of the sinusoidal wave.
This paper is organized as follows. In Section 2, the PO method is introduced and applied on the simplified sea surface model. All the formulas are derived and all the expressions are given in this part; both horizontal (HH) polarization and vertical (VV) polarizations are considered. Several numerical simulations are exhibited in Section 3 to show the validity and efficiency of the new model compared with the realistic sea surface. Then, this model is used to investigate the characteristic of the Doppler spectrum of time varying sea surface. Section 4 ends with a summary of the new model and a proposition for further pertinent investigation.
2.1. The Physical Optics Formulation
The initial point of physical optics is the surface currents produced by an incoming electromagnetic wave (). Considering a 1D sea surface, the induced electric currents and magnetic currents on each segment (the length of each segment is set as 1/8~1/10 wavelength of the incident wave to meet the division criterion) are given by where is the unit normal vector of the surface and indicates the position of each segment. and are, respectively, the total electric and magnetic fields at the surface. For TE case, the horizontal polarization vector of the incident wave is along , is the wavenumber of the incident wave, and the vertical polarization vector is . The scattering problems in the global coordinate system and the local coordinate system are shown in Figures 1(a) and 1(b), respectively. Given the slop of one segment, the local coordinate system can be built as , and each coordinate component can be expressed as
Let a plane wave illuminate on the rough sea surface; then the incident field on the surface can be given aswhere is the position of each surface segment and is the wave impedance in the free space.
Inserting the TE reflection coefficient , the electric and magnetic reflected fields arewhere . Adding the contributions of incident and reflected field, the total electric and magnetic fields can be written as a function of the incident electric field, the reflection coefficients, and the intrinsic impedance of the first medium. Then, going back to the boundary conditions (1), the equivalent currents can be reckoned from the total fields. Consider
The field scattered by the dielectric sea surface can readily be determined as the radiation of and . The scattered field can be given by the well-known Stratton-Chu formula; that is,where is the unit normal vector of each segment. The resulting expressions are written in terms of Green’s function and . The far field approximation allows expressing these functions as . Considerwhere is the scattering wavenumber vector . Consequently, the scattered fields are
Similarly, when the sea surface is illuminated by the incident field (VV polarization), the induced electric and magnetic currents can be written asthe scattered fields are
The analytical expression of scattering coefficient can be defined aswherein is the length of the sea surface.
The Doppler spectrum is a power spectral density of the random time varying complex amplitude of the sea surface scattering field; to evaluate it, we used a standard spectral estimation technique. Thus, the expression of the Doppler spectrum is defined aswhere the angular bracket stands for the averaging over all the surface realizations. is the evolution time of the time varying sea surface.
2.2. The Simplified Sea Surface Model
In this simplified sea surface model, each segment on the sea surface is much longer than the wavelength of the incident wave. All segments are locally approximated by tilted facet. (In the following part of the paper, the “large segments” are all named as “facet.”) These facets are centered on the grid points that are geometrically described by a discrete set of . can be regarded as the gravity wave. The sinusoidal capillary waves are added on the planar facet. The surface profile within a planar facet can be represented aswhere is the height of the simplified sea surface, represents the amplitude of the superimposed capillary wave, and is the spatially short wavenumber, which is supposed to satisfy the resonant condition of Bragg waves. The definitions of and will be discussed in the following in detail; represents the circle frequency, , and m/s, rad/m; is the random phase of the ripple component. is the height of the gravity wave in the middle point of the facet, and is the slope of the facet.
The configuration of the sea surface is depicted in Figure 2. Figure 2(a) shows the large-scale surface. Each facet is flat in this model. In Figure 2(b), the capillary waves are added and each facet is modified by sinusoidal waves.
Since the sinusoidal ripple is superimposed at microscopic level and is relevant for the Bragg resonant scattering mechanism, we assume that the echo waves are reflected by the wave crest of the resonant sinusoidal ripples. As Figure 3(a) shows, , and is its projection on the plane facet. is the wavelength of the resonant sinusoidal ripples on the plane facet. is the corresponding wavelength on the -axis. To meet the resonant condition, is set as
The wavenumber of the sinusoidal ripple is , which can be written as . According to (15), can be determined asand the amplitude of the ripple in (14) could be expressed by ; is the capillary spectrum located in the higher part of the sea spectrum, and is the length of the facet.
Figure 3(b) shows that the facet is divided into periods, . The currents in the first period can be expressed as (5), (6), and (10). Each period is divided into small segments. The lengths of the small segments are set as 1/8~1/10 wavelength of the incident wave to meet the division criterion. For the case of HH polarization, by moving the term from (5) and (6) to (8) and moving the term from (8) to (5) and (6), in (5) is removed, since it will counteract with the in (7). Then the currents can be rewritten aswhere the subscripts and indicate the th segment in one period and the th period on one facet, respectively. Consider . , which is the slope of each segment, is periodical on one facet. is the differential of . On each facet, can be expressed asthe scattered field from one facet can be expressed aswhere is the number of the segments in one period, , and Since and vary periodically on one facet, (19) can be rewritten aswhere is fixed on one facet, and it can be expressed as . Thus the scattered field from one facet can be obtained by
Summarizing the scattered field from all the facets, we can obtain the total field from the sea surface.
3. Numerical Simulations and Discussions
To validate the proposed simplified sea surface model, simulation results obtained from the simplified sea surface, the realistic sea surface, and the simplified sea without capillary wave will be shown and compared in the following.
Figure 4 shows the comparison of the presented simplified model with the realistic sea surface for the angular distribution of the backscattering coefficient (BSC). The geometry of observation is monostatic. Both HH (a) and VV (b) polarization are taken into consideration. The results are obtained by averaging 50 surface realizations. The surfaces are generated using Monte Carlo  methods with 50 sets of Gaussian random numbers whose mean value is 0 and variance is 1. The random numbers are called using the IMSL function library with 50 different seeds. The parameters in the simulation are as follows: the frequency of the incident wave is at the X-band (10 GHz), the incident angle varies from 0° to 70°, the length of the sea surface is 245.76 m, the length of each facet is set to 1.92 m, the wind speed at 10 m above the mean sea level is m/s, and the relative permittivity of the sea water is calculated as according to the Klein model  at 20°C and 32.5% of salinity. It is obvious that the BSC obtained from the proposed sea surface model and the realistic model are in fairly good agreement for the whole backscattering region for both polarizations. The BSC obtained by a simplified sea surface without capillary wave have an obvious discrepancy with the other two curves. Thus, the presented model is reliable for predicting the backscattering coefficient. We can also conclude that the resonant capillary wave plays an important role in the scattering result. The simulation time obtained from the simplified sea surface and realistic sea surface is also presented in Table 1. The superiority of the simplified sea surface in terms of computation time is obvious. All results are obtained on a computer with a 2.93 GHz processor (Intel Core i3 CPU), 3.45 GB memory, and Visual Fortran 6.5 compiler. All the comparisons of scattering results and the simulating time indicate that the proposed sea surface model is an accurate and efficient model.
Subsequently, the Doppler behaviors are observed with incident angle and frequency GHz. The results obtained from the proposed simplified sea surface are compared with the result from the realistic sea surface and from the simplified sea without capillary wave (Figure 5). These results are averaged over 50 samples and normalized by their respective maximum values. For the time-evolving simulations, the time step is selected as 0.005 s, and 1024 steps are selected for the sea surface realization. Thus, the evolution time is 5.12 s. The other parameters are the same as those given in Figure 4. The two curves obtained from the proposed sea surface model and the realistic model well agree for both polarizations. Doppler spectrum obtained from the simplified sea without capillary waves cannot reflect the Doppler characteristics of the sea surface. It should be noticed that the Doppler shifts are almost the same for both methods, whereas the Doppler width of the proposed method is a little smaller than the realistic sea surface. It is mainly caused by the fact that the other components of the capillary waves are supposed to be absent. After all, we conclude that the simplified sea surface model is accurate enough and convincible to describe the scattering characteristics of the sea.
The simulation times of Doppler spectrum of both sea surface models are presented in Table 1. The simplified sea surface model takes much less time than the realistic sea surface, especially in the Doppler spectrum simulation.
In Figure 6, the Doppler spectra are illustrated for wind speed of 5 m/s at different incident angles. The other parameters in calculation are the same as those given in Figure 5. With the incident angle increasing, the widths of the Doppler spectra first become broader and then shrink. Figure 7 shows the Doppler spectrum at different wind speed. The incident angle is set as . We can find that with the increasing wind speeds, the velocities of the surface water and the orbital motion for the large-scale sea wave increase, and there has been a gradual increase for the Doppler spectrum frequency shift, which is coincident with the measuring result by Rozenberg et al. . On the other hand, the Doppler spectrum broadens with the wind speed increasing in the roughness of the sea surface. Thus the backscattering energy will be distributed over a wide region of the frequency domain, which has been discussed in [12, 13]. All these results can demonstrate that the simplified sea surface model is valid to analyze the Doppler characteristic of the time varying sea surface.
The PO method is combined with the simplified sea surface model to calculate the backscattering coefficient and Doppler spectrum from time varying sea surface. With the simplified model, the efficiency of PO method is promoted. The characteristic of the Doppler spectrum of time varying sea surface is also investigated by the proposed model. Some shortcomings still exist in this method. Since the coupling between different segments are not taken into consideration, this method is limited to deal with the sea whose sea state is low or moderate. And this method will be inaccurate if the incident angle is large. The 2D problem is not considered in this paper. We will pay more attention to these remaining problems in the future study.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by the National Natural Science Foundation for Distinguished Young Scholars of China (Grant no. 61225002), the Specialized Research Fund for the Doctoral Program of Higher Education (Grant no. 20100203110016), and the Fundamental Research Funds for the Central Universities (Grant no. K5051007001).
F. T. Ulaby, R. K. Moore, and A. K. Fung, Microwave Remote Sesing (Active and Passive), Addison Wesley, New York, NY, USA, 1982.
L. Tsang and J. A. Kong, Scattering of Electromagnetic Waves, vol. 2, John Wiley & Sons, 2001.
Y. W. Wei, L. X. Guo, A. Q. Wang, and Z. S. Wu, “Application of multiregion model to EM scattering from a dielectric rough surface with or without a target above it,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 11, pp. 5607–5620, 2013.View at: Publisher Site | Google Scholar | MathSciNet
J. T. Johnson, “On the geometrical optics and physical optics approximations for scattering from exponentially correlated surfaces,” IEEE Transactions on Antennas and Propagation, vol. 45, pp. 2619–2629, 1994.View at: Google Scholar
L. A. Klein and C. T. Swift, “An improved model for the dielectric constant of sea water at microwave frequencies,” IEEE Transactions on Antennas and Propagation, vol. AP-25, no. 1, pp. 104–111, 1977.View at: Google Scholar
J. V. Toporkov and G. S. Brown, “Numerical study of the extended Kirchhoff approach and the lowest order small slope approximation for scattering from ocean-like surfaces: doppler analysis,” IEEE Transactions on Antennas and Propagation, vol. 50, no. 4, pp. 417–425, 2002.View at: Publisher Site | Google Scholar