Electromagnetic Scattering and Doppler Spectrum Simulation of Land–Sea Junction Composite Rough Surface

: In this paper, a weighted arctangent function is used in conjunction with the spectral method to generate a land–sea junction composite rough surface under the spatially homogeneous and time-stationary hypotheses. The exponential correlation function and the Joint North Sea Wave Project (JONSWAP) spectrum, combined with an experiment-veriﬁed shoaling coefﬁcient, are applied to model the land surfaces and the time-varying sea surfaces separately. The second-order small slope approximation (SSA-II) with tapered wave incidence is utilized for evaluating the electromagnetic scattering characteristics and Doppler characteristics of the generated composite rough surface. The inﬂuence of land–sea interface factors on radar cross-section (RCS) and Doppler shift of radar echoes is investigated in detail by comparing the RCS and Doppler spectra of the land–sea junction composite rough surfaces with those of ﬁnite-depth sea surfaces. It can be found that the Doppler spectra of the land–sea junction composite rough surface is narrower than that of the ﬁnite-depth sea surface under upwind directions and wider than that of the ﬁnite-depth sea surface under crosswind directions.


Introduction
With the development of computational electromagnetic (CEM), the electromagnetic (EM) scattering characteristics and Doppler characteristics of random rough surface have received a great deal of attention in recent years from a range of applications, such as radar imaging, surface remote sensing, and ocean wave spectra estimation. Thus far, several methods have been proposed to evaluate the EM scattering from a rough surface. For example, the method of moments (MoM) combined with Rao-Wilton-Glisson (RWG) basis functions and Poggio-Miller-Chang-Harrington-Wu (PMCHW) integral equations is proposed and utilized to solve EM scattering from a rough surface in [1]. The multilevel fast multipole algorithm (MLFMA), together with the impedance boundary conditions (IBC) and the self-dual integral equation, is applied to calculate the backscattering from a rough sea surface under low grazing incidence [2]. A stochastic solution of the EM scattering from a penetrable, randomly rough surface is derived by using the vector-based-finite-element method (FEM) in [3]. Li et al. combined the Finite-Difference Time-Domain (FDTD) approach with the uniaxial perfectly matched layer (UPML) absorbing boundary to solve the EM scattering problems of 1-D rough sea surfaces [4], 2-D rough sea surfaces [5], and two-layered rough surfaces [6]. The time domain integral equation (TDIE) is applied to simulate the transient scattering from the ocean surface with ship wake in [7]. In [8], a second stochastic degree iterative algorithm with sparse matrix (SM) and Chebyshev approximation is proposed to calculate the fully polarimetric bistatic scattering behavior that resulted from the interaction of an EM wave with all scales of ocean waves. Saddek Afifi uses the first-order small slope approximation (SSA-I) and the first-order small perturbation method (SPM-I) to study the EM scattering from a 2-D random rough surface that separates the vacuum from a perfect electromagnetic conductor (PEMC) [9]. The second-order small slope approximation (SSA-II) is applied to solve the bis-and mono-static scattering from linear and non-linear rough sea surfaces in [10][11][12]. The second-order small perturbation method (SPM-II) is applied to solve the EM scattering from layered mediums with rough interfaces by Hasan Zamani [13,14] and R. J. Burkholder [15]. A semi-deterministic facet model (SDFM) is proposed in [16] to simulate the EM scattering from an ocean-like surface.
In [17], the advanced integral equation model (AIEM) is applied to predict the backscattering, bistatic scattering, and emission of rough soil surfaces. Wang and Tong proposed an improved facet-based model derived from a two-scale model (TSM) to evaluate the scattering strength of an electrically large ocean surface [18]. In [19][20][21][22][23], the first-order SSA and TSM models are utilized to solve monostatic and bistatic scattering from randomly rough surfaces (sea surface covered with or without oil film, ground surface) in different frequency bands.
However, in all of the literature mentioned above, the processed rough surface is the sea surface or ground surface, which contains only one type of background surface and one power spectrum. It is possible that both the sea surface and ground surface exist in one beam in the real radar detection. Due to the difference in the dielectric parameters and the surface profile between the land-sea junction composite rough surface and the sea (ground) surface, the EM scattering characteristics of the land-sea junction composite rough surface are correspondingly different from those of the sea surface and those of the ground surface. Motivated by the reasons illustrated above, our focus in this paper is on the influence of land-sea interface factors on radar cross-section (RCS) and Doppler spectra of EM scattered echoes of the land-sea junction composite rough surface, which has a frozen land surface and a time-evolving sea surface. The research in this paper is potentially valuable for remote sensing, target detection in land-sea junction areas, etc. Due to desirable properties such as analytical tractability and numerical efficiency, in this paper the nonlinear hydrodynamic model of the JONSWAP spectrum combined with an experiment-verified shoaling coefficient is used to model a finite-depth sea surface, and the exponential correlation function is adopted to model the ground surface. A weighted arctangent function is introduced to smoothly connect the sea surface and the ground surface, and to model the land-sea junction composite rough surface. Compared with the classical models such as Kirchhoff approximation (KA), TSM, and SPM, SSA-II takes into account facets' tilt modulation and second-order Bragg scattering and is thus capable of reasonably predicting the RCS and Doppler spectra. Therefore, the SSA-II is employed to evaluate the RCS and Doppler spectra of radar echoes from sea surfaces and land-sea junction composite rough surfaces in this study.
In [24][25][26][27][28], several methods have been proposed to solve electromagnetic scattering from a land-sea junction composite rough surface or a target located above a coastal environment. This paper has advantages over the literature cited previously. First, the SSA-II is applied as the electromagnetic scattering kernel since it is more accurate than the techniques in the aforementioned literature and can account for second-order Bragg scattering. Second, this paper focuses on the analysis of Doppler characteristics of the land-sea junction composite rough surface, which is not involved in the prior literature. The weighted arctangent function was originally introduced in [29] to generate a land-sea junction composite rough surface, but there are several differences between [29] and this paper. First, although the weighted arctangent function has been introduced both in [29] and the manuscript under review, the two functions have different forms. Second, the Modified Equivalent Current Approximation (MECA) is used to calculate the EM scattering fields from the composite rough surface in [29], which can only take into account the specular reflection effect. In this paper, the SSA-II is applied as the electromagnetic scattering kernel, which is more accurate than MECA and can take into account second-order Bragg scattering. Third, the EM scattering from a target located in the land-sea junction area is the primary issue addressed in [29]. In this paper, the RCS and Doppler spectra of the land-sea junction composite rough surface are calculated and discussed in detail. The remainder of this paper is organized as follows: in Section 2, the weighted arctangent function is introduced, and the modeling of the land-sea junction composite rough surface is presented. Section 3 presents the SSA-II EM scattering model with tapered wave incidence for evaluating the RCS of rough surfaces (sea surface and land-sea junction composite rough surface). A comparison of the numerical results of the RCS and Doppler spectra for finite-depth sea surface and land-sea junction composite rough surfaces is presented and discussed in Section 4. Section 5 is devoted to the conclusions of this paper.

Sea Spectrum of Finite-Depth Sea Water
With the change of the water depth and seabed topography, the profile and the statistical characteristics of the sea surface in coastal areas are different from those in deep sea areas. In [30], a shoaling coefficient is proposed by McCormick to study the effect of water depth on the wave height, which can be written as in which d and k represent the depth and the wave number of the sea surface, respectively. According to McCormick's theory, the spectrum of actively growing wind waves in finite-depth seas can be analytically expressed as where S finite (ω) and S finite (k) are the spectrums of finite-depth seas. S JONSWAP (ω) denotes the JONSWAP spectrum, ω is the angular frequency of sea waves, which can be related according to the shallow water gravity capillarity dispersion relation as where g 0 is the acceleration of gravity and k m = 363.2rad/m represents the wave number with minimum phase speed. Figure 1 depicts the influence of water depth on the finite-depth sea spectrum. In Figure 1, the wind fetch involved in the finite-depth sea spectrum is set to 50 km, and the wind speed at 10 m height is set to u 10 = 10 m/s. It can be seen from Figure 1a that the peak value of the curves first decreases and then increases as the water depth decreases. Furthermore, it is evident that the peaks of angular frequency ω remain unchanged with the decrease of water depth, whereas the peaks of wave number k shift toward the high-frequency part as the water depth decreases.

Rough Surface Realizations
To analyze the radar echoes (RCS, Doppler spectra), the geometric properties of the land-sea junction composite rough surface must be captured. In this paper, the spectral

Rough Surface Realizations
To analyze the radar echoes (RCS, Doppler spectra), the geometric properties of the land-sea junction composite rough surface must be captured. In this paper, the spectral method under the spatially homogeneous and time-stationary hypothesis is applied to yield the rough finite-depth sea surface and ground surface, which is highly efficient due to the introduction of the fast Fourier transform. To establish a rough surface with the spectral method, a linear superposition of harmonic waves whose amplitudes are independent, normally distributed random values times the square root of the sea/ground surface spatial spectrum is generated, and the Fourier amplitude at any time t can be expressed as where k is the two-dimensional spatial wave vector, S(k, ϕ) is the two-dimensional sea/ground spectrum (for the finite-depth sea surface, Elfouhaily wind direction function is applied), and δk x = 2π/L x and δk y = 2π/L y represent the sampling interval along x and y directions, respectively. γ(k) is a complex Gaussian series with zero mean and unit standard deviation.
Thus, the sea/ground surface elevation h(r, t) at position r = (x, y) and time t can be expressed by Equation (6) can be efficiently accomplished by introducing the inverse fast Fourier transform. To ensure that h(r, t) is real, A(k, t) is required to satisfy the conjugate symmetry (Hermitian form) as follows: Figure 2 shows the simulated finite-depth sea surface under different water depths at a wind speed of u 10 = 10 m/s. It can be found that as the water depth decreases, the major sea wavelength decreases and, as a result, the major spatial wave number shifts to a high frequency, which is similar to that in Figure 1.
The exponential correlation function is utilized in the ground surface modeling for the realization of the land-sea junction composite rough surface.

Land-Sea Junction Composite Rough Surface Realizations
Once the finite-depth sea surface and the ground surface are generated, the only problem in realizing land-sea junction composite rough surface is how to smoothly connect the two rough surfaces with different undulations. The problem is solved in this study by introducing a weighted arctangent function, which can be described as: in which f com (x, y) represents the land-sea junction composite rough surface, f gro (x, y) and f sea (x, y) are the ground surface and finite-depth sea surface, respectively. The weighted factors ω 1 and ω 2 can be written as ω 1 = [π/2 − arctan(y − y b )]/π and ω 2 = [π/2 + arctan(y − y b )]/π, where y b is the boundary between the ground surface and the sea surface. Figure 3 is a simple example of how the weighted arctangent function works. In Figure 3, the weighted arctangent function shown in Equation (7) is applied to smoothly connect the two flat plates (z = 0.1 for y < 0 area, z = 0.5 for y > 0 area). The area surrounding the boundary (y = 0) of the two flat plates can be found to have continuous values and slopes.  The exponential correlation function is utilized in the ground surface modeling for the realization of the land-sea junction composite rough surface.

Land-Sea Junction Composite Rough Surface Realizations
Once the finite-depth sea surface and the ground surface are generated, the only problem in realizing land-sea junction composite rough surface is how to smoothly connect the two rough surfaces with different undulations. The problem is solved in this study by introducing a weighted arctangent function, which can be described as:      Figures 4 and 5, y < 0 area is the finite-depth sea surface and the y > 0 area is the ground surface. In the simulation, the length of the composite rough surface along the xand y-axes is 256 m, and the wind speed and wind fetch are set as 5 m/s and 50 km, respectively. The water depth is 1 m. In Figure 4, the correlation length and root mean square height are l x = l y = 3.0 m and δ = 0.1 m. In Figure 5, the correlation length and root mean square height are l x = l y = 10.0 m and δ = 1.0 m. The wind direction in Figures 4a and 5a is set as ϕ = 90 • , which means that the wind blew towards the shore. In Figures 4b and 5b, the wind direction is set as ϕ = 0 • , which means that the wind blew along the coastline. The surface height and slope of the coastline area are continuous, as seen in Figures 4 and 5.
is 1 m. In Figure 4, the correlation length and root mean square height are 3.0m x y l l = = and 0.1m δ = . In Figure 5, the correlation length and root mean square height are 10.0m which means that the wind blew towards the shore. In Figures 4b and 5b, the wind direction is set as 0 ϕ =  , which means that the wind blew along the coastline. The surface height and slope of the coastline area are continuous, as seen in Figures 4 and 5. (a) Figure 5. Simulated land-sea junction composite rough surface 10.0m x y l l = = 1.0m δ = .

SSA-II Model for EM Scattering from Land-Sea Junction Composite Rough Surface
The SSA proposed by Voronovich consists of a basic approximation of theory (SSA-I) and a second-order correction to it (SSA-II). The SSA-I only contains first-order Bragg scattering, which results in the fact that the SSA-I cannot predict the depolarization of wave scattering from a rough surface in the plane of incidence. Compared with other is 1 m. In Figure 4, the correlation length and root mean square height are 3.0m x y l l = = and 0.1m δ = . In Figure 5, the correlation length and root mean square height are 10.0m which means that the wind blew towards the shore. In Figures 4b and 5b, the wind direction is set as 0 ϕ =  , which means that the wind blew along the coastline. The surface height and slope of the coastline area are continuous, as seen in Figures 4 and 5. (a) (a) Figure 5. Simulated land-sea junction composite rough surface 10.0m x y l l = = 1.0m δ = .

SSA-II Model for EM Scattering from Land-Sea Junction Composite Rough Surface
The SSA proposed by Voronovich consists of a basic approximation of theory (SSA-I) and a second-order correction to it (SSA-II). The SSA-I only contains first-order Bragg scattering, which results in the fact that the SSA-I cannot predict the depolarization of wave scattering from a rough surface in the plane of incidence. Compared with other

SSA-II Model for EM Scattering from Land-Sea Junction Composite Rough Surface
The SSA proposed by Voronovich consists of a basic approximation of theory (SSA-I) and a second-order correction to it (SSA-II). The SSA-I only contains first-order Bragg scattering, which results in the fact that the SSA-I cannot predict the depolarization of wave scattering from a rough surface in the plane of incidence. Compared with other classical models (such as KA, TSM, and SPM), the SSA-II has its own advantages: (1) it takes into account the mutual transformation of the two linear polarization states caused by facets' tilts; (2) it takes into account the second-order Bragg scattering and is thus capable of reasonably predicting the depolarized scattering from rough surfaces both in and outside the plane of incidence; and (3) a continuous spectrum of roughness is required, and it is not necessary to divide roughness into large-and small-scale components (just like it is in TSM). Due to the advantages mentioned above, the SSA-II model is used to evaluate the bistatic RCS, monostatic RCS, and Doppler spectra of the land-sea junction composite rough surface.
The SSA is based on the transformation properties of scattering amplitude (SA) with respect to vertical shifts. It claims that the slope of roughness is generally the only small parameter underlying this theory. Consider a tapered plane wave illuminated on a land-sea junction composite rough surface to reduce the edge effect caused by the limited surface size of L x × L y , which can be expressed as where t x = (x cos θ inc cos ϕ inc + y cos θ inc sin ϕ inc + z sin θ inc ) 2 g 2 cos 2 θ inc (10)  (12) in which g represents the tapered wave factor, which is set as max L x , L y /6 in the following simulations. θ inc and ϕ inc denote the incident angle and incident azimuth angle, respectively. The geometry of the scattering problem and corresponding notation are shown in Figure 6. θ sca and ϕ sca denote the scattering angle and scattering azimuth angle of the receiver, respectively. k inc = k 0 + q 0ẑ and k sca = k 1 + qẑ are the incident wave vector and scattering wave vector, respectively; k 0 and k 1 denote the horizontal components of incident wave vector and scattering wave vector; and q 0 and q are the vertical components of incident wave vector and scattering wave vector.ĥ andv represent the polarizations of incident and scattering waves.  (11) in which g represents the tapered wave factor, which is set as max{ , } / 6 x y L L in the following simulations. inc θ and inc ϕ denote the incident angle and incident azimuth angle, respectively. The geometry of the scattering problem and corresponding notation are shown in Figure 6. sca θ and sca ϕ denote the scattering angle and scattering azimuth angle of the receiver, respectively. 0 0înc q z = + k k and 1ŝca qz = + k k are the incident wave vector and scattering wave vector, respectively; 0 k and 1 k denote the horizontal components of incident wave vector and scattering wave vector; and 0 q and q are the vertical components of incident wave vector and scattering wave vector. ĥ and v represent the polarizations of incident and scattering waves. Figure 6. Geometry of the land-sea junction composite rough surface scattering problem and notation.
Thus, after introducing the tapered plane incident wave, the SA of the SSA-II model can be expressed as where represents the Fourier transform of the surface elevation. P inc denotes the incident wave power captured by the rough surface. The kernel functions B and M are 2 × 2 matrices describing transitions of the EM waves of various polarizations into each other, which are mainly dependent on the configuration angles, polarizations, boundary conditions, and the permittivity of the lower medium. In fact, Bragg's scattering is described by function B in SSA, and the term relates to the function M proportional to the slopes of a rough surface rather than to the elevations themselves. The corresponding details of these two kernel functions can be found in [31]. Thus, the average scattering coefficient can be expressed as where the angle brackets denote the ensemble average over random surface realizations. Furthermore, the Doppler spectrum can be defined as the power spectral density of the random time-varying scattering amplitude and can be evaluated by utilizing a standard spectral estimation technique by the following equation where S(k 1 , k 0 ; t) is the backscattering amplitude from rough surface at time t. It should be noted that when calculating the Doppler spectra of the land-sea junction composite rough surface, the sea surface is time-evolving, whereas the ground surface is constant and time-independent.
To quantitatively measure the Doppler spectrum, the Doppler shift f c , that is, the spectral centroid, and the bandwidth of the Doppler spectrum f w can be defined as where f c represents the horizontal motion between ocean waves and observation radar, which can be used to invert the propulsion velocity of ocean waves. f w is mainly caused by the up-and-down motion of the points on the sea surface profile.

Numerical Results and Discussion
To validate SSA-II, a comparison of monostatic scattering coefficients of the sea surface derived from SSA-II with those from 3 months of Seasat microwave scatterometer (SASS) measurements [32] is performed and shown in Figure 7 under wind speed u 10 = 5 m/s and u 10 = 10 m/s for Ku-band (14.6 GHz) to validate SSA-II. According to the Debye expression [33], the relative permittivity of sea water is ε sea r = 44.78 + i39.12 at sea water temperature of 20 • C and a salinity of 30 parts per thousand. The sampling interval is λ/8, where λ is electromagnetic wavelength. The size of the sea surface is L x = L y = 10000λ. The wind fetch is 100 km and 300 km. The wind direction is set as ϕ = 0 • , and the tapered incident wave is illuminated upwind to the sea surface with tapering parameter g = λ/6. The final monostatic scattering coefficient is an ensemble average of 50 realizations of sea surface.

Numerical Results and Discussion
To validate SSA-II, a comparison of monostatic scattering coefficients of the sea surface derived from SSA-II with those from 3 months of Seasat microwave scatterometer (SASS) measurements [32] is performed and shown in Figure 7 under wind speed   It can be seen in Figure 7 that SSA-II is very consistent with the experimental data for both HH and VV polarizations. Additionally, the VV polarization has more accurate results than that of the HH polarization, especially at 10 10m/s u = . It is because the nonlinear effect of the sea surface is more obvious under HH polarization.
In the following, numerical simulations are performed at a frequency of 1.2GHz f =  It can be seen in Figure 7 that SSA-II is very consistent with the experimental data for both HH and VV polarizations. Additionally, the VV polarization has more accurate results than that of the HH polarization, especially at u 10 = 10 m/s. It is because the nonlinear effect of the sea surface is more obvious under HH polarization.
In the following, numerical simulations are performed at a frequency of f = 1.2GHz. The relative permittivity of sea surface and ground surface are ε sea r = 73.2 + i67.2 and ε gro r = 4.54 + i0.36, respectively. The length of three types of rough surfaces (sea surface, ground surface, and land-sea junction composite rough surface) is L x = L y = 128λ, with sampling interval of λ/8. The tapering parament is set to be g = λ/6. The water depth and wind fetch are fixed at 5 m and 30 km, respectively. For the ground surface, the correlation length and root mean square height are l x = l y = 3.0 m and δ = 0.1 m, respectively. All numerical results are derived by averaging 50 surface realizations. Figure 8 is a comparison of the bistatic and monostatic scattering coefficients of the land-sea junction composite rough surface derived by SSA-II and MoM. In the simulation, the wind speed is u 10 = 5 m/s in an upwind direction. The incident angle is θ i = 30 • , and the incident azimuth angle and scattering azimuth angle are both set to 0 • . It can be seen from Figure 8 Figure 9 shows the comparison of the bistatic scattering coefficients of the finitedepth sea surface, the ground surface, and the land-sea junction composite rough surface. The wind speed, incident angle, and incident (scattering) azimuth angle are the same as those in Figure 8. Figure 9 shows that under HH and VV polarizations, the bistatic scattering coefficients of three different types of rough surfaces exhibit spikes in the mirror direction of the incident wave. Regardless of polarization, the bistatic scattering coefficient of a land-sea junction composite rough surface is similar to but smaller than that of a finite-depth sea surface. The bistatic scattering coefficient of the ground surface is much smaller than that of the sea surface and that of the land-sea junction composite rough surface's overall scattering angles. This is because the ground surface is much flatter than the sea surface in this simulation.   Figure 9 shows the comparison of the bistatic scattering coefficients of the finite-depth sea surface, the ground surface, and the land-sea junction composite rough surface. The wind speed, incident angle, and incident (scattering) azimuth angle are the same as those in Figure 8. Figure 9 shows that under HH and VV polarizations, the bistatic scattering coefficients of three different types of rough surfaces exhibit spikes in the mirror direction of the incident wave. Regardless of polarization, the bistatic scattering coefficient of a land-sea junction composite rough surface is similar to but smaller than that of a finitedepth sea surface. The bistatic scattering coefficient of the ground surface is much smaller than that of the sea surface and that of the land-sea junction composite rough surface's overall scattering angles. This is because the ground surface is much flatter than the sea surface in this simulation. Figure 10 depicts a comparison of the monostatic scattering coefficients of the finite-depth sea surface, the ground surface, and the land-sea junction composite rough surface. The parameters used in the simulation are the same as those in Figure 9. A similar conclusion can be derived from Figures 9 and 10.
The comparison of the bistatic scattering coefficients of a land-sea junction composite rough surface with different root mean square heights is shown in Figure 11. The simulating parameters are the same as those in Figure 9, except that the root mean square of the ground region is set at 0.1 m, 0.3 m, and 0.5 m. From Figure 11, it can be seen that the peak value of the bistatic scattering coefficient decreases with the increase of the root mean square of ground for HH and VV polarizations. However, for HH and VV polarizations, the bistatic scattering coefficient increases with the increase of the root mean square of the ground at the moderate and large scattering angles (−40 • ≤ θ s ≤ 15 • and 40 • ≤ θ s ≤ 90 • ). This is because as the root mean square increases, so does the degree of ground fluctuation, and more energy is scattered rather than reflected. It can also be found that the scattering coefficient increases with the increase of the root mean square of the ground when the scattering angle θ s ≥ −60 • for cross-polarization.
The monostatic scattering coefficients of the land-sea junction composite rough surface versus different wind speeds and wind directions are illustrated in Figure 12. The simulating parameters are the same as in Figure 9, except that the wind speed is set at 5m/s and 10 m/s in upwind and crosswind directions, respectively. As shown in Figure 12, the monostatic scattering coefficient of the land-sea junction composite rough surface decreases with the increase of the wind speed at moderate and small incident angles (θ i ≤ 60 • ) and increases with the increase of the wind speed at large incident angles (θ i > 60 • ) under HH and VV polarizations. However, for cross-polarization, the monostatic scattering coefficient increases with the increase in wind speed over all incident angles. It can also be found that the monostatic scattering coefficient under cross wind is smaller than that under upwind for co-polarization and is greater than that under upwind for cross-polarization. Figures 13 and 14 show the average Doppler spectra of backscattered echoes from the finite-depth sea surface and the land-sea junction composite rough surface. The surface area, the corresponding sample interval, and the incident wave frequency are the same as those in the aforementioned scattering coefficient simulation. The wind speeds are u 10 = 5 m/s and u 10 = 10 m/s in Figures 13 and 14, respectively. The radar is looking upwind, with the incident angle changing from 5 • to 80 • . All of the results were averaged over 50 realizations of rough surfaces. Both co-polarization and cross-polarization results are presented. Due to the reciprocity of HV polarization and VH polarization in backscattering, only HV polarized Doppler spectra are presented. The corresponding Bragg frequencies f B = g sin θ i /(πλ) are also plotted in Figures 13 and 14 with vertical short-dotted lines.  Figure 9 shows the comparison of the bistatic scattering coefficients of the finitedepth sea surface, the ground surface, and the land-sea junction composite rough surface. The wind speed, incident angle, and incident (scattering) azimuth angle are the same as those in Figure 8. Figure 9 shows that under HH and VV polarizations, the bistatic scattering coefficients of three different types of rough surfaces exhibit spikes in the mirror direction of the incident wave. Regardless of polarization, the bistatic scattering coefficient of a land-sea junction composite rough surface is similar to but smaller than that of a finite-depth sea surface. The bistatic scattering coefficient of the ground surface is much smaller than that of the sea surface and that of the land-sea junction composite rough surface's overall scattering angles. This is because the ground surface is much flatter than the sea surface in this simulation.   Figure 10 depicts a comparison of the monostatic scattering coefficients of the finitedepth sea surface, the ground surface, and the land-sea junction composite rough surface. The parameters used in the simulation are the same as those in Figure 9. A similar conclusion can be derived from Figures 9 and 10. The comparison of the bistatic scattering coefficients of a land-sea junction composite rough surface with different root mean square heights is shown in Figure 11. The simulating parameters are the same as those in Figure 9, except that the root mean square of the ground region is set at 0.1 m, 0.3 m, and 0.5 m. From Figure 11, it can be seen that the peak value of the bistatic scattering coefficient decreases with the increase of the root mean square of ground for HH and VV polarizations. However, for HH and VV polarizations, the bistatic scattering coefficient increases with the increase of the root mean square of the . This is because as the root mean square increases, so does the degree of ground fluctuation, and more energy is scattered rather than reflected. It can also be found that the scattering coefficient increases with the increase of the root mean square of the ground when the scattering angle 60 s θ ≥ −  for cross-polarization.    The monostatic scattering coefficients of the land-sea junction composite rough surface versus different wind speeds and wind directions are illustrated in Figure 12. The simulating parameters are the same as in Figure 9, except that the wind speed is set at 5m/s and 10 m/s in upwind and crosswind directions, respectively. As shown in Figure  12, the monostatic scattering coefficient of the land-sea junction composite rough surface decreases with the increase of the wind speed at moderate and small incident angles ( 60 i θ ≤  ) and increases with the increase of the wind speed at large incident angles ( 60 i θ >  ) under HH and VV polarizations. However, for cross-polarization, the monostatic scattering coefficient increases with the increase in wind speed over all incident angles. It can also be found that the monostatic scattering coefficient under cross wind is smaller than that under upwind for co-polarization and is greater than that under upwind for cross-polarization.  The left column, the middle column, and the right column is for VV, HH, and HV Figure 13. Average Doppler spectra of backscattered echoes from the finite-depth sea surface and the land-sea junction composite rough surface at various incident angles with a wind speed of u 10 = 5 m/s. The left column, the middle column, and the right column is for VV, HH, and HV polarizations, respectively. The vertical short-dotted lines represent Bragg frequencies. The radar is looking upwind.   Figures 13 and 14 show that as the incident angle increases, the bandwidth of the Doppler spectrum widens first and then narrows, and the spectral peak moves closer to the corresponding Bragg frequency for both the finite-depth sea surface and the land-sea junction composite rough surface. Compared with co-polarization, the cross-polarized  Figures 13 and 14 show that as the incident angle increases, the bandwidth of the Doppler spectrum widens first and then narrows, and the spectral peak moves closer to the corresponding Bragg frequency for both the finite-depth sea surface and the land-sea junction composite rough surface. Compared with co-polarization, the cross-polarized Doppler spectra are much closer to the corresponding Bragg frequency. This attributes to the fact that the specular reflection plays a leading role in small and moderate incident angles, whereas the Bragg scattering is the primary scattering mechanism in large incident angles. It can also be observed that the Doppler spectra of the land-sea junction composite rough surface are narrower than those of the finite-depth sea surface for both co-polarization and cross-polarization. In addition, it is readily observed that the movement tendency of the Doppler spectral peak is different between co-polarization and cross-polarization. For the co-polarization case, the corresponding frequency of the Doppler spectral peak of the land-sea junction composite rough surface is smaller than that of the finite-depth sea surface, while for the cross-polarization case, the corresponding frequency of the Doppler spectral peak of the land-sea junction composite rough surface is greater than that of the finite-depth sea surface.
A quantitative comparison of the Doppler shifts (calculated by Equation (16)) and Doppler spectral bandwidths (calculated by Equation (16)) between the finite-depth sea surface and the land-sea junction composite rough surface is shown in Figures 15 and 16 Figures 13 and 14 show that as the incident angle increases, the bandwidth of the Doppler spectrum widens first and then narrows, and the spectral peak moves closer to the corresponding Bragg frequency for both the finite-depth sea surface and the land-sea junction composite rough surface. Compared with co-polarization, the cross-polarized Doppler spectra are much closer to the corresponding Bragg frequency. This attributes to the fact that the specular reflection plays a leading role in small and moderate incident angles, whereas the Bragg scattering is the primary scattering mechanism in large incident angles. It can also be observed that the Doppler spectra of the land-sea junction composite rough surface are narrower than those of the finite-depth sea surface for both co-polarization and cross-polarization. In addition, it is readily observed that the movement tendency of the Doppler spectral peak is different between co-polarization and cross-polarization. For the co-polarization case, the corresponding frequency of the Doppler spectral peak of the land-sea junction composite rough surface is smaller than that of the finite-depth sea surface, while for the cross-polarization case, the corresponding frequency of the Doppler spectral peak of the land-sea junction composite rough surface is greater than that of the finite-depth sea surface.
A quantitative comparison of the Doppler shifts (calculated by Equation (16)) and Doppler spectral bandwidths (calculated by Equation (16)) between the finite-depth sea surface and the land-sea junction composite rough surface is shown in Figures 15 and 16  A comparison of average Doppler spectra of backscattered echoes from finite-depth sea surface and land-sea junction composite rough surface under cross wind is shown in Figure 17. The simulating parameters are the same as those in Figure 13, except that the wind direction is set to 90 w ϕ =  , which indicates that the radar is looking crosswind. It can be readily observed that only one peak appears when the incident angle is smaller A comparison of average Doppler spectra of backscattered echoes from finite-depth sea surface and land-sea junction composite rough surface under cross wind is shown in Figure 17. The simulating parameters are the same as those in Figure 13, except that the wind direction is set to ϕ w = 90 • , which indicates that the radar is looking crosswind. It can be readily observed that only one peak appears when the incident angle is smaller than 30 • under co-polarization and smaller than 60 • under cross-polarization. However, with the increase in incident angle, two distinct peaks appear for both co-and cross-polarization. These two peaks correspond to the fact that the ocean waves are moving toward and away from the radar. Different from Figures 13-16, it can be obviously observed from Figure 17 that the Doppler spectra of the land-sea junction composite rough surface is wider than that of the finite-depth sea surface. In addition, the Doppler spectral peak of the land-sea junction composite rough surface shifts to anhigher frequency relative to that of the finite-depth sea surface.  Figure 13, except that the wind direction is ϕ w = 90 • .

Conclusions
In this paper, a land-sea junction composite rough surface is generated by introducing a weighted arctangent function. The SSA-II, together with a tapered plane incident wave, is applied to obtain the scattering efficient and Doppler spectra of the finite sea surface, the ground surface, and the land-sea junction composite rough surface. The simulation results show that both the bis-and mono-static scattering coefficients of the land-sea junction composite rough surface are smaller than those of the finite sea surface and greater than those of the ground surface for co-and cross-polarizations. By comparing the Doppler spectra of the land-sea junction composite rough surface with that of the finite-depth sea surface, it can be found that a narrowing of the Doppler spectra is observed under upwind direction and a widening of the Doppler spectra is observed under crosswind direction in the presence of the ground surface for both co-and cross-polarizations. Finally, it should be pointed out that only straight shorelines can be modeled with the approach proposed in this paper.