Ocean Wave Measurement Using Short-Range K-Band Narrow Beam Continuous Wave Radar

We describe a technique to measure ocean wave period, height and direction. The technique is based on the characteristics of transmission and backscattering of short-range K-band narrow beam continuous wave radar at the sea surface. The short-range K-band radar transmits and receives continuous signals close to the sea surface at a low-grazing angle. By sensing the motions of a dominant facet at the sea surface that strongly scatters signals back and is located directly in front of the radar, the wave orbital velocity can be measured from the Doppler shift of the received radar signal. The period, height and direction of ocean wave are determined from the relationships among wave orbital velocity, ocean wave characteristics and the Doppler shift. Numerical simulations were performed to validate that the dominant facet exists and ocean waves are measured by sensing its motion. Validation experiments were conducted in a wave tank to verify the feasibility of the proposed ocean wave measurement method. The results of simulations and experiments demonstrate the effectiveness of the short-range K-band narrow beam continuous wave radar for the measurement of ocean waves.


Introduction
There are two general approaches for the measurement of ocean waves: in situ and remote sensing measurements.The in situ measurements are usually made with the instrument located in the proximity of the site of interest and typically produces a more accurate local result.For example, motion and attitude sensors mounted on a wave rider buoy [1], moored at the sea surface, detect local wave motions.Acoustic and pressure sensors fixed below the sea surface measure distance to the surface and water pressure to detect wave elevation; radar and laser wave sensors mounted above the surface look downward observing wave elevation by measuring the distance to the surface [2].Remote sensing is conducted at a location distant from sea surface.HF radar [3][4][5][6] is employed to measure the sea wave spectrum, wave height, wind field, and surface current velocity.Both non-coherent and coherent microwave radars have been used to measure wave spectra and estimate wave characteristics at sea surface [7][8][9][10].ISR CORrad Digital Imaging Radar [11] and Sigma S6 WaMoS II Wave Radar [12] are two commercially available X-band marine radars for which the wave measurements are based on the analysis of various radar images.Synthetic-aperture radar employs L-band and C-band to measure directional and power spectra of wave fields with interferometric configuration [13] and measure near-surface ocean currents with the phase difference between Synthetic-aperture radar image where G T is the gain of radar transmitter, G R is the gain of radar receiver, σ is the sea surface RCS, λ is the wavelength of the continuous radar wave, and R is the range to the scattering area.P T , G T , G R , and λ are the inherent properties of the radar, σ is the property of the sea surface, and R is related to the relative position of the scattering area to the radar.
If the short-range K-band NBCWR is set close to the sea surface and transmits at low-grazing angle, a dominant facet F 0 will form and it dominates the backscattered power to the radar.The formation of the dominant facet is due to the narrow radar beam width in the azimuth and a sharp attenuation in the sea surface RCS for a low-grazing K-band radar.The narrow radar beam in the azimuth concentrates the radar's transmission and reception at a small angle, and the sharp attenuation of RCS ensures that the near surface scatters much higher power back to radar than the far surface.The radar's low height can reduce the size of the dominant facet.Its size is much smaller than the gravity waves so that it can be considered as a surface object moving with orbital motion of the gravity waves and, furthermore, it is much greater than the radar wavelength so that it can produce Bragg or breaking wave scattering.The orbital motions of the dominant facet cause the Doppler shift frequency f d which can be obtained through analyzing the spectrogram of the continuous IF outputs.The period and the amplitude of the Doppler shift frequency variation are related to the wave period, height, and direction.

Measurement of Wave Period, Height and Direction
Wave measurement by measuring the orbital motions of gravity waves using X-band radar (8~12 GHz) has already been conducted [7,11].Here, the relationship among the wave orbital velocity, ocean wave characteristics, and Doppler shift frequency is elaborated with the orbital motions of the dominant facet.If T, H, and α denote the wave period, height, and direction, respectively, and suppose v w is the wave orbital velocity, the three relationships are conceptually given as γ 1 : T, H, α v w , γ 2 : v w f d , and γ 3 : f d T, H, α, and illustrated in Figure 1.
Wave measurement by measuring the orbital motions of gravity waves using X-band radar (8~12 GHz) has already been conducted [7,11].Here, the relationship among the wave orbital velocity, ocean wave characteristics, and Doppler shift frequency is elaborated with the orbital motions of the dominant facet.If , , and  denote the wave period, height, and direction, respectively, and suppose  is the wave orbital velocity, the three relationships are conceptually given as  : , ,    ,  :    , and  :   , , , and illustrated in Figure 1.Theoretically, γ 1 is defined by linear wave theory and γ 2 through the Doppler effect.γ 3 is consequently deduced from γ 1 and γ 2 .These relationships suggest that ocean wave characteristics can be obtained from the deterministic relationship for γ 3 if the Doppler shift frequency f d is measured.
The relationship for γ 1 is discussed first.On the basis of linear wave theory [20], the orbit of the orbital motion of ocean wave is circular in deep water or elliptical in intermediate and shallow waters.The wave orbital velocity v w can be decomposed into a horizontal component v x and a vertical component v z : v z = πH T sinhk(d + z) sinhkd sin(kx − ωt), (6) where, k is the wave number, d is the water depth, ω is the wave angular frequency and equal to 2π/T, x is the horizontal coordinate, z is the vertical coordinate, and t is the time.v x points toward the wave direction α and v z is vertically upward.Equations ( 5) and ( 6) express the relationship γ 1 .The velocity components v x and v z can be approximately calculated by setting the coordinates x and z to a fixed point in the orbit [21].Therefore, given H, T, k, d, and ω, and assuming the coordinates x and z to be the mean position (x 0 , z 0 ) of the orbit, Equations ( 5) and ( 6) become periodic cosine and sine functions of time t with period of T.
Figure 2 illustrates a dominant facet F 0 (black dot) that is formed in the radiation of a short-range K-band NBCWR and traveling in a clockwise circle along with a long surface wave which propagates from left to right.The radar grazing angle is θ g .The radial velocity to the radar, v LOS , is the vector projection of wave orbital velocity v w on the radar line of sight (LOS) and can be calculated as: range K-band NBCWR and traveling in a clockwise circle along with a long surface wave which propagates from left to right.The radar grazing angle is  .The radial velocity to the radar,  , is the vector projection of wave orbital velocity  on the radar line of sight (LOS) and can be calculated as: Equation ( 7) is similar to Equation ( 5) except for the constant values  and  , and the phase difference  between  and  .Furthermore, supposing that  is a constant,  is also a periodic cosine function with period of .If a continuous radar signal with a stable frequency  is transmitted to the dominant facet  along the radar LOS, the Doppler effect modifies the continuous wave, and the Doppler shift frequency  is calculated with: Figure 2. Velocities of the dominant facet F 0 orbiting with a long surface wave.The radar is not shown here, but it transmits from right to left along the LOS and it is close to the sea surface.v x , v z and v LOS are three projections of v w to x-axis, z-axis, and LOS, respectively.Compared with v x , the phase of v LOS is delayed by θ g .Equation ( 7) is similar to Equation ( 5) except for the constant values x 0 and z 0 , and the phase difference θ g between v x and v LOS .Furthermore, supposing that θ g is a constant, v LOS is also a periodic cosine function with period of T. If a continuous radar signal with a stable frequency f T is transmitted to the dominant facet F 0 along the radar LOS, the Doppler effect modifies the continuous wave, and the Doppler shift frequency f d is calculated with: where c is the speed of light.f d is proportional to the radial velocity v LOS , so f d can also be expressed by a periodic cosine function with period of T and the term before the cos notation is the amplitude of the cosine function.Equation (8) expresses the relationship for γ 2 .
According to the relationships γ 1 and γ 2 , we know that the wave period coincides with the periods of radial velocity v LOS and Doppler shift frequency f d : where T LOS is the period of radial velocity v LOS , and T d is the period of the Doppler shift frequency variation.The Doppler shift frequency reaches its amplitude as the dominant facet F 0 moving along the LOS, i.e., cos kx 0 − ωt − θ g = 1.Therefore, the wave height can be deduced from Equation (8) and calculated from: where A d is the amplitude of the Doppler shift frequency variation.In Equation ( 10), T has the measured value from Equation ( 9), k = 2π/L, and the wavelength L is calculated with dispersion relation for ocean waves [22].Equations ( 9) and ( 10) express the relationship for γ 3 .Consequently, the wave height T and wave period H are obtained through measuring the period T d and the amplitude A d of the Doppler shift frequency variation.The Doppler shift frequency can be measured by sensing the orbital motion of the dominant facet F 0 at the sea surface with the short-range K-band NBCWR which functions as a super heterodyne radar [23].
The measured wave period is considered as the true value because the period of orbital motion does not change in the radar view even if the radar is not transmitting along the wave direction.However, the measured wave height is not necessarily the true value because the angle between the wave direction and radar looking direction (here referred to as angle to the wave, β) will affect the wave height estimation.In Equation (10), the measured amplitude A d of the Doppler shift frequency variation decreases with the coefficient |cos(β)|.In practice, the wave direction α and the angle to the wave β are unknown.This problem can be solved using the geometrical relationship between the true value and the several measured values.Figure 3 shows an example of the geometrical relationship for three measured wave heights.
by sensing the orbital motion of the dominant facet  at the sea surface with the short-range K-band NBCWR which functions as a super heterodyne radar [23].
The measured wave period is considered as the true value because the period of orbital motion does not change in the radar view even if the radar is not transmitting along the wave direction.However, the measured wave height is not necessarily the true value because the angle between the wave direction and radar looking direction (here referred to as angle to the wave, ) will affect the wave height estimation.In Equation (10), the measured amplitude  of the Doppler shift frequency variation decreases with the coefficient |cos ()|.In practice, the wave direction α and the angle to the wave  are unknown.This problem can be solved using the geometrical relationship between the true value and the several measured values.Figure 3 shows an example of the geometrical relationship for three measured wave heights.Three vectors can be illustrated for three different magnitudes  ,  , and  and measurement directions.The three created vectors  ⃗ ,  ⃗ , and  ⃗ are the true wave height projections in the measurement directions.By combining the three vectors at point A and constructing three perpendiculars at the end of each vector, the intersection point B is determined.The length of the line segment AB is the true value of wave height and the wave direction runs from H 3 are the true wave height projections in the measurement directions.By combining the three vectors at point A and constructing three perpendiculars at the end of each vector, the intersection point B is determined.The length of the line segment AB is the true value of wave height and the wave direction runs from point A to point B or the reverse.The angles β 1 , β 2 , and β 3 between line segment AB and the three vectors are also measured, respectively.The wave direction α is the angle between line segment AB and the north.In practice, a more precise result can be produced if more measured wave heights are applied.

Multiple Ocean Waves
In reality, many waves are simultaneously present at the sea surface traveling with different periods, heights, and directions.The orbital motion of the dominant facet F 0 depends on the propagation of all the waves.The wave orbital velocity v w is: where v n is a velocity component which originates from the nth ocean wave, and N is the number of ocean waves.Therefore, the radial velocity to the radar v LOS and the Doppler shift frequency f d are calculated as: where v LOS n is the radial velocity of component v n to radar, and f d n is the nth Doppler shift frequency component.Equation (13) implies that f d results from the sum of all the radial velocities and contains different frequency components, i.e., f d is the linear combination of a set of periodic cosine functions of different periods T d and amplitudes A d .It is known that the periods and the amplitudes can be demonstrated through analyzing the spectrum of f d .According to Equation (9), the resultant periods are the wave periods.Meanwhile, wave height can be calculated by Equation (10) for each wave.

Numerical Simulations
The measurement method is first validated via numerical simulations.The important thing is to prove that a dominant facet exists at the sea surface when the K-band NBCWR transmits and receives signals close to sea surface at low-grazing angle.Due to the short detection range, the length and width of the wave field are limited within a few hundred meters.In order to describe the orbital motion of the gravity wave and emphasize that the measurement method is based on sensing a small dominant facet at the sea surface, we use a regular wave in the wave field which is generated using the spectral peak in the Pierson-Moskowitz spectrum [24].Each facet in the wave field scatters energy back to the radar.The received power from each facet can be calculated using Equation (4) in which the antenna pattern gains and the sea surface RCS for the K-band radar are two important parameters.In the numerical simulations, we use the typical antenna pattern gains of a K-band radar transceiver module and the calculated K-band RCS for Bragg backscattering.Through calculating the IF power received from each facet in the wave field and comparing the powers of all facets, we find that the dominant facet exists and oscillates with wave propagations.
The time series of Doppler shift frequency caused by the motions of the dominant facet can be identified and extracted from spectrogram of IF signals.We can find significant peaks in the amplitude spectrum of the time series of Doppler shift frequency.The amplitude spectrum is calculated using Fourier Transforms.Each peak gives two values: T d and A d .Consequently, the wave period and height are obtained by substituting these two values into Equations ( 9) and (10).Different wave heights can be measured if different radar looking directions are used.The true wave height is retrieved according to the method shown in Figure 3, in which wave direction is also measured.All the details of the process of the numerical simulations will be described in the following sub-sections.

Wave Field
By collecting the data of surface waves under differing wind conditions, the ocean wave spectra were established to describe the wave energy distribution for different wave frequencies.The Pierson-Moskowitz spectrum, which defines a wave spectrum for a fully developed sea, is given by: where f is wave frequency, α 0 is equal to 8.1 × 10 −3 , g is the gravitational acceleration, and U is the wind speed at an elevation of 19.5 m above sea surface.The significant wave height H s and the wave frequency at the spectral peak f p can be deduced with Equation ( 14) and expressed by the following equations: Table 1 shows the characteristics of ocean waves at the spectral peaks under wind conditions from 8-18 m/s.The wavelength L p is calculated with dispersion relation for deep water (d = 1000 m).The wave at the spectral peak is dominant because it has more energy than the others.The dominant wave generated by the 8 m/s wind is used first.Figure 4 shows the wave field with the generated ocean waves.The modeled sea surface size is 200 m × 300 m, and the unit facet size is 1 m × 1 m.The wave is propagating towards 140 degrees, indicated with the upper arrow.The colored stripes reveal the waveform.The radar transceiver is set 2 m above the sea level and horizontally points to the north.The small black circle at (0, 0) is its location and the lower arrow indicates the direction of the radar transceiver.Each facet is much smaller than the waves and the orbital velocities of the points in one facet are nearly identical.In addition to the waves shown in the wave field, there are also wind-generated capillary waves which are not shown in the wave field but generate backscatter.Therefore, the facets at the surface are not only orbiting with the gravity waves but also scattering radar signals back to the radar.The impacts of capillary waves will be incorporated in the calculation of RCS of sea surface.

Radar IF Output Signal
The IF output signal is regarded as the sum of all the IF components that are received from all facets at the sea surface in Figure 4 and already processed through the radar transceiver.The IF component from each facet is given by Equation (3).The Doppler shift frequency f d can be calculated by Equation (8).It should be noted that the angle to wave β must be included in the calculation.The phases of incidence and scattered waves in a single facet have possible values ranging from 0 to 2π since the radar wavelength λ is much less than the facet size.Therefore, the phase ϕ of each IF component from one facet has a random value in [0, 2π].The power of the IF component is calculated using the radar power equation, expressed by Equation ( 4).
The radar transmitter gain G T is the antenna pattern gains, while the radar receiver gain G R is composed of the global antenna gain, the amplifier gain and the antennae pattern gains [25].They are expressed as: where G T (θ a , θ e ) and G R (θ a , θ e ) are respectively the gains of transmitter and receiver in a direction denoted by the azimuth angle θ a and elevation angle θ e .G pa (θ a ) and G pe (θ e ) are antenna pattern gains in the azimuth and elevation, respectively, G an and G am are the global antenna gain and total gain of the amplifiers, respectively.All the gains in Equations ( 16) and ( 17) have values in decibels.In the numerical simulation, the antenna pattern gains of VV (vertical transmit and vertical receive) polarized K-MC3 transceiver module made by RFbeam Microwave GmbH in St. Gallen, Switzerland will be used, which transmits 24 GHz continuous wave with narrow beam widths, seven degrees in azimuth and 25 degrees in elevation, as shown in Figure 5.The seven-degree beam width in the azimuth enables the radar transceiver to transmit and receive in a specified direction; and the 25 degrees beam width in elevation includes not only the far facets but also the close ones at the sea surface.Another important parameter in Equation ( 4) is the sea surface RCS () .The microwave backscattering from the sea surface has been investigated by many researchers.One of the dominant scattering mechanisms at microwave frequencies and low-to-medium grazing angle is Bragg scattering, and non-Bragg scattering from the breaking waves was also observed and studied at a low-grazing angle [26][27][28][29].According to the electromagnetic scattering perturbation theory, the firstorder backscatter cross-sections were established in [30].The radar cross-section for backscattering from capillary waves for VV and HH (horizontal transmit and horizontal receive) polarization and the results of measurements in the wind wave tank with an X-band radar was introduced in [31].A mathematical model of normalized radar cross-section (NRCS) for Bragg backscattering and breaking wave in equilibrium conditions was presented in [32].Microwave Doppler spectra models based on Bragg-scattering and composite-surface theory were developed and used to show backscattering from rough water surfaces under many wind speed and incidence angle conditions [33].The backscattering mechanism from rough sea surface for both Ku-and C-bands was studied with the small-slope approximation method and a theoretical model for numerical calculations of a radar backscattering cross section was presented in [34].
The transceiver module has VV polarization so that the backscattering is dominated by Bragg scattering.Therefore, we use the mathematical model of NRCS for Bragg backscattering established under equilibrium condition [32]: where  is the attenuation coefficient,  the direction of observation relative to wind,  the incidence angle,  * the friction velocity,  the radar wavenumber, and  () the first-order scattering coefficient for VV polarization [35]: where  is the relative complex dielectric constant of sea water.Although the NRCS model was established for equilibrium conditions, Phillips compared it with some radar experiments and concluded that it can be extended to a wide range of applications.Table 2 shows the parameters used for the NRCS calculation.Another important parameter in Equation ( 4) is the sea surface RCS (σ).The microwave backscattering from the sea surface has been investigated by many researchers.One of the dominant scattering mechanisms at microwave frequencies and low-to-medium grazing angle is Bragg scattering, and non-Bragg scattering from the breaking waves was also observed and studied at a low-grazing angle [26][27][28][29].According to the electromagnetic scattering perturbation theory, the first-order backscatter cross-sections were established in [30].The radar cross-section for backscattering from capillary waves for VV and HH (horizontal transmit and horizontal receive) polarization and the results of measurements in the wind wave tank with an X-band radar was introduced in [31].A mathematical model of normalized radar cross-section (NRCS) for Bragg backscattering and breaking wave in equilibrium conditions was presented in [32].Microwave Doppler spectra models based on Bragg-scattering and composite-surface theory were developed and used to show backscattering from rough water surfaces under many wind speed and incidence angle conditions [33].The backscattering mechanism from rough sea surface for both Ku-and C-bands was studied with the small-slope approximation method and a theoretical model for numerical calculations of a radar backscattering cross section was presented in [34].
The transceiver module has VV polarization so that the backscattering is dominated by Bragg scattering.Therefore, we use the mathematical model of NRCS for Bragg backscattering established under equilibrium condition [32]: where β 0 is the attenuation coefficient, Φ the direction of observation relative to wind, θ the incidence angle, u * the friction velocity, κ the radar wavenumber, and F vv (θ) the first-order scattering coefficient for VV polarization [35]: where ε r is the relative complex dielectric constant of sea water.Although the NRCS model was established for equilibrium conditions, Phillips compared it with some radar experiments and concluded that it can be extended to a wide range of applications.Table 2 shows the parameters used for the NRCS calculation.Here, ε r is the relative dielectric constant for sea water with salinity of 35 ppt and a temperature of 0 • C [36], u * is a measured value with 7.3 m/s wind speed [37].The calculated NRCS is shown in Figure 6.The value at 90-degree is not shown in the figure because of much lower than others.The backscattering from the simulated sea surface is weak because the radar transceiver transmits and receives at a low-grazing angle, i.e., most of the incidence angles to the facets are approaching 90 degrees.Therefore, depressing the radar transceiver slightly is recommended to increase the radar grazing angle.Here the grazing angle is five degrees.Here,  is the relative dielectric constant for sea water with salinity of 35 ppt and a temperature of 0 °C [36],  * is a measured value with 7.3 m/s wind speed [37].The calculated NRCS is shown in Figure 6.The value at 90-degree is not shown in the figure because of much lower than others.The backscattering from the simulated sea surface is weak because the radar transceiver transmits and receives at a low-grazing angle, i.e., most of the incidence angles to the facets are approaching 90 degrees.Therefore, depressing the radar transceiver slightly is recommended to increase the radar grazing angle.Here the grazing angle is five degrees.Table 3 shows the parameters of numerical simulation conditions about ocean wave, wave field, radar transceiver, and IF signal.

Parameters
Values Unit Period () 5.85 s Height () 1.37 m Figure 6.NRCS σ 0 for a VV polarized antenna (upwind).This is a calculated NRCS for a K-band radar under certain wind and sea water conditions as shown in Table 2.For the given simulation conditions, the IF component power from each facet can be determined from Equation ( 4).The IF component power distribution, as shown in Figure 7, is obtained if the IF component powers are laid out on the wave field.The blank areas are the geometric shadows in which the facets cannot be seen by the radar transceiver.For microwave backscattering, partial shadowing [38] may cause some signals from the geometrically shadowed areas.However, they are invisible and below the noise floor because of low incident power.The warm color areas, which stand for the relatively strong scattering, are confined within the narrow azimuthal beam width of the radar transceiver.For the given simulation conditions, the IF component power from each facet can be determined from Equation ( 4).The IF component power distribution, as shown in Figure 7, is obtained if the IF component powers are laid out on the wave field.The blank areas are the geometric shadows in which the facets cannot be seen by the radar transceiver.For microwave backscattering, partial shadowing [38] may cause some signals from the geometrically shadowed areas.However, they are invisible and below the noise floor because of low incident power.The warm color areas, which stand for the relatively strong scattering, are confined within the narrow azimuthal beam width of the radar transceiver.
Figure 7. IF component power distribution.The strongest power −53 dBW occurs at the dominant facet which is located at (0, 7) in the main lobe.Some directional spreading around the main lobe also presents because of the side lobes as shown in Figure 5, but all of them are over 100 dBW under the power from dominant facet so that they have hardly any impact on the measurement.
From the geometric viewpoint, the closer the wave approaches to the radar, the more surface of the wave the radar can see (see Figure 8).Thus, the white area at shorter range is narrower than those at longer ranges.Contrarily, the white areas are gradually getting wider at long range, as shown in Figure 7. Strong scattering appears at locations close to the radar transceiver, and Figure 9 shows the 3-D plot of the IF power distribution of the strong scattering area which is marked as the dashed rectangle with a size of 55 m × 20 m in Figure 7. Figure 10 shows the powers distribution with the Figure 7. IF component power distribution.The strongest power −53 dBW occurs at the dominant facet which is located at (0, 7) in the main lobe.Some directional spreading around the main lobe also presents because of the side lobes as shown in Figure 5, but all of them are over 100 dBW under the power from dominant facet so that they have hardly any impact on the measurement.
From the geometric viewpoint, the closer the wave approaches to the radar, the more surface of the wave the radar can see (see Figure 8).Thus, the white area at shorter range is narrower than those at longer ranges.Contrarily, the white areas are gradually getting wider at long range, as shown in Figure 7. Strong scattering appears at locations close to the radar transceiver, and Figure 9 shows the 3-D plot of the IF power distribution of the strong scattering area which is marked as the dashed rectangle with a size of 55 m × 20 m in Figure 7. Figure 10 shows the powers distribution with the west-east coordinates from −2 m to 2 m.The power level at dominant facet the F 0 is much higher than for the side facets whose coordinates are −2 m, −1 m, 1 m, and 2 m.          Figure 11 shows the power levels and Doppler shift frequencies of dominant facet F 0 and its four adjacent facets with the west-east coordinate being 0 m.The power of F 0 is not much higher than that of the two facets next to it, F 1 and F 2 , but the Doppler shift frequencies of these three facets are almost same.However, the power of F 0 is much higher than that of the two facets F 3 and F 4 .
Remote Sens. 2018, 10, x FOR PEER REVIEW 13 of 26 Figure 11 shows the power levels and Doppler shift frequencies of dominant facet  and its four adjacent facets with the west-east coordinate being 0 m.The power of  is not much higher than that of the two facets next to it,  and  , but the Doppler shift frequencies of these three facets are almost same.However, the power of  is much higher than that of the two facets  and  .Figure 12 shows the horizontal velocity and vertical velocity of the dominant facet F 0 in five wave periods (29.25 s).The curves are sinusoids and oscillate with same amplitude, and the phase difference is π/2.Therefore, the motions of the dominant facet F 0 represent the orbital motions of the generated ocean waves.Figure 12 shows the horizontal velocity and vertical velocity of the dominant facet  in five wave periods (29.25 s).The curves are sinusoids and oscillate with same amplitude, and the phase difference is π/2.Therefore, the motions of the dominant facet  represent the orbital motions of the generated ocean waves.5) and ( 6), in which the variation of spatial position (x, z) of the dominant facet is also taken into account.5) and ( 6), in which the variation of spatial position (x, z) of the dominant facet is also taken into account.

Signal Processing and Results
The  9) and ( 10), the measured wave period and wave height are 5.86 s and 1.05 m, respectively.

Signal Processing and Results
The calculations of the power  , frequency  , and phase  of the IF component originated from each facet of the simulated ocean waves have been introduced.The total IF signal is the sum of the IF components and noise components.The IF signal and its spectrogram are shown in Figure 13.In the spectrogram, a bright line appears and is almost synchronized with the horizontal velocity of the dominant facet  in Figure 12.This implies that the dominant facet  dominates the received power from the simulated sea surface, and the bright line illustrates the variation of the Doppler shift frequency caused by the motions of the dominant facet  .The bright line is extracted by finding the maximum value at each moment.The extracted bright line, i.e., the time series of the Doppler shift frequency  , and its amplitude spectrum are plotted in Figure 14.A significant peak can be seen in the amplitude spectrum and the corresponding amplitude and frequency are 90.68Hz and 0.1706 Hz.Therefore,  is 90.68 Hz and  is 5.86 s.Substituting  and  into Equations ( 9) and ( 10), the measured wave period and wave height are 5.86 s and 1.05 m, respectively.In addition to the significant peak at 0.1706 Hz, some other secondary peaks show up at the integer multiples of the frequency 0.1706 Hz because the velocities of the dominant facet are calculated according to Equations ( 5) and ( 6), which contain not only sinusoidal terms but also hyperbolic terms.According to the Taylor Series expansion of hyperbolic functions, the terms of cosh and sinh with variable z can be expressed as an infinite sum of terms which contain z n , n = 0, 1, 2, 3, ….Therefore the harmonics of the wave frequency exist in the velocities of the dominant facet as the result of the vertical coordinate z varying with same period of the generated wave.
Two other simulations were executed for two radar directions, 30 degrees and 330 degrees.Two spectrograms of the IF signals are shown in Figure 15.Two extracted Doppler shift frequencies and two amplitude spectra are shown in Figure 16.In addition to the significant peak at 0.1706 Hz, some other secondary peaks show up at the integer multiples of the frequency 0.1706 Hz because the velocities of the dominant facet are calculated according to Equations ( 5) and ( 6), which contain not only sinusoidal terms but also hyperbolic terms.According to the Taylor Series expansion of hyperbolic functions, the terms of cosh and sinh with variable z can be expressed as an infinite sum of terms which contain z n , n = 0, 1, 2, 3, . . . .Therefore the harmonics of the wave frequency exist in the velocities of the dominant facet as the result of the vertical coordinate z varying with same period of the generated wave.
Two other simulations were executed for two radar directions, 30 degrees and 330 degrees.Two spectrograms of the IF signals are shown in Figure 15.Two extracted Doppler shift frequencies and two amplitude spectra are shown in Figure 16.In addition to the significant peak at 0.1706 Hz, some other secondary peaks show up at the integer multiples of the frequency 0.1706 Hz because the velocities of the dominant facet are calculated according to Equations ( 5) and ( 6), which contain not only sinusoidal terms but also hyperbolic terms.According to the Taylor Series expansion of hyperbolic functions, the terms of cosh and sinh with variable z can be expressed as an infinite sum of terms which contain z n , n = 0, 1, 2, 3, ….Therefore the harmonics of the wave frequency exist in the velocities of the dominant facet as the result of the vertical coordinate z varying with same period of the generated wave.
Two other simulations were executed for two radar directions, 30 degrees and 330 degrees.Two spectrograms of the IF signals are shown in Figure 15.Two extracted Doppler shift frequencies and two amplitude spectra are shown in Figure 16.Similarly, two sets of results of wave period and height are obtained and shown in Table 4 with the first results.The three measured wave periods are identical, 5.84 Hz, which are very close to the period of the simulated ocean wave.The error mainly comes from the frequency resolution of the amplitude spectrum.Three measured wave heights are less than 1.37 m because of the angles to the wave direction.Similarly, the true wave height can be retrieved with three measured wave heights and three radar directions, as shown in Figure 17.Similarly, two sets of results of wave period and height are obtained and shown in Table 4 with the first results.The three measured wave periods are identical, 5.84 Hz, which are very close to the period of the simulated ocean wave.The error mainly comes from the frequency resolution of the amplitude spectrum.Three measured wave heights are less than 1.37 m because of the angles to the wave direction.Similarly, the true wave height can be retrieved with three measured wave heights and three radar directions, as shown in Figure 17.In this case, these perpendicular lines cannot cross because of the wave height measurement errors.Here, the centroid B of the intersecting points of the perpendiculars is used.The length of the line segment AB, 1.33 m, is the final result of wave height measurement and the wave direction is measured to be about 142.13 degrees.The resultant measurement approximates the wave height and direction in the simulation conditions but with errors of 3.7% and 1.5%, respectively.There are three main error sources: the frequency resolution in the spectrogram of IF signal; the noise components in In this case, these perpendicular lines cannot cross because of the wave height measurement errors.Here, the centroid B of the intersecting points of the perpendiculars is used.The length of the line segment AB, 1.33 m, is the final result of wave height measurement and the wave direction is measured to be about 142.13 degrees.The resultant measurement approximates the wave height and direction in the simulation conditions but with errors of 3.7% and respectively.There are three main error sources: the frequency resolution in the spectrogram of IF signal; the noise components in IF signal; the constructive or destructive scattering originated from the facets adjacent to the dominant facet F 0 .The errors can be greatly decreased if the frequency resolution and the extraction are improved.

Measurement of Multiple Ocean Waves
Considering the ocean waves generated by a distant storm and propagating into the wave field shown in Figure 4, two ocean waves present and the respective characteristics are shown in Table 5.The wave 1 is identical with the previous simulations.The period of wave 2 is 10.97 s which is consistent with the period T p generated by 15 m/s wind, but the wave height has reduced to 1 m.Wave 1 and wave 2 are superposed in the wave field (see Figure 18).Figure 19 includes the in-phase IF signal and the spectrogram of the IF signal in fifteen periods of wave 1 (87.75 s).The bright line in the spectrogram demonstrates the fluctuations of the Doppler shift frequency, and it obviously contains multiple frequency components.Figure 20 shows the extracted Doppler shift frequency and its amplitude spectrum.Two significant peaks stand out indicating the amplitudes and frequencies for the two waves.Table 6 shows the simulation results as well as some wave characteristics.The wave heights are calculated with Equation (10) and amended with the angles to wave.The measured wave periods and the amended wave heights are very close to the settings shown in Table 5.Furthermore, the true values of the wave heights can be obtained by measuring in several different directions when the wave directions and the angles to wave are unknown.Table 6 shows the simulation results as well as some wave characteristics.The wave heights are calculated with Equation (10) and amended with the angles to wave.The measured wave periods and the amended wave heights are very close to the settings shown in Table 5.Furthermore, the true values of the wave heights can be obtained by measuring in several different directions when the wave directions and the angles to wave are unknown.The validation experiment was conducted in the wave tank of the Fluids Laboratory of Memorial University of Newfoundland.The tank is 50 m long and 4.5 m wide, and 1.8 m deep.The short-range K-band NBCWR was mounted on the towing carriage about 2 m above the water surface.Figure 21 shows the radar installation in the wave tank.Additionally, two sensors were also used to profile the generated waves.The experimental configurations are shown in Table 7.The wave period and height were set by the wave generator of the wave tank.The wavelength was measured with the two wave sensors.Figure 22 shows the wave profiles output of one of the wave sensors.It can be seen that six waves passed in 30 s because of the 5-s wave period, and the troughs are wider than the crests because the troughs were slowed down by the friction from the bottom of the tank.The measurement results are shown in Table 8 and the bottom row gives the wave height obtained with the geometrical relationship of the two results, H 1 and H 2 .The measured period (4.81 s) and height (0.1552 m) are very close to the configured values of the generated wave, and errors are less than 3.8% and 3.5%, respectively.The measurement accuracy could be improved if more signals were sampled and more directions were applied.

Discussion
We performed the numerical simulation and validated the existence of the dominant facet, as shown in Figures 10 and 11.Given space limitations, we cannot present more cases with different simulation conditions.However, the waves were successfully measured in the simulation for multiple waves and the experiment in wave tank.It implies that dominant facets existed and oscillated with wave propagations.The dominant facet just has a size of a few square meters and is very close to the radar.Therefore, the technique provides a possibility of measuring ocean waves by sensing a small facet in a limited local area but not scanning a large area of sea surface to generate a radar image.The technique can be applied to surface craft and offshore platform to measure local ocean waves for enhancing operation safety.
In order to improve measurement accuracy, the size of dominant facet must be formed as small as possible because the orbital velocities of the points on a small dominant facet are much more uniform than those on a large one.The facet size can be reduced by using radar transceiver with a much narrower beam and setting the radar transceiver much closer to the sea surface.For this reason, the technique cannot be applied if the sensor is installed at a relatively high location, in which case a large dominant facet is formed, but inconsistent orbital velocities exist.
In the numerical simulations and the wave tank experiment, we also found that there were always errors in the measured heights, periods, and directions.There are three main error sources: the frequency resolution in the spectrogram of IF signal; the noise components in IF signal; and the constructive or destructive scattering originated from the facets adjacent to the dominant facet.As shown in Figure 17, the true wave height measurement has error because of the errors of three wave height measurements.This also causes wave direction measurement errors as well as a 180-degree ambiguity.However, the errors can be greatly decreased if the frequency resolution and the extraction method are improved.The frequency resolution can be improved by extending the total sampling time because the frequency resolution of amplitude spectrum is equal to the reciprocal of total sampling time.With a higher frequency resolution, we can select the significant peaks at the frequencies much closer to the true wave frequencies.In additions, the extracted time series of Doppler shift frequency have some noise components which might introduce errors even if the regular waves are used in numerical simulations and tank experiments.The error may be reduced by considering the weighted average of the values around the maximum value at a certain moment on bright line.Improved method and a study of the effect of sampling time on measurement error at complex sea sate will be investigated in the future.
Meanwhile, we considered only the wave propagations which contributed the periodical components of Doppler shift frequency.In practice, the phase velocity of Bragg scattering wave, velocity of the surface current, and wind-induced surface drift usually add non-periodical components to the Doppler shift frequency.However, these non-periodical components can be identified and removed in the phase of data processing.In addition, platform motions cause additional components which can be identified by recording the platform's relative velocity to the average position of the dominate facet.These operations cannot completely remove undesired components but the remainders are weak and negligible.
Regular waves were used in the numerical simulations and the wave tank experiments.Therefore, the orbital motions of the regular waves can be explained by linear wave theory, and the wave periods and heights were calculated using Equations ( 9) and (10), which are based on linear wave theory.The impacts of wind and wave age are not included.The simulated waves had the same period and height as the spectral peak wave in the Pierson-Moskowitz spectrum; the waves in the tank were generated by the back-and-forth motions of the wave generator but not by wind.However wind waves at the real sea surface are not regular waves.The horizontal and vertical components of the dominant facet's velocities caused by wind waves can be explained more precisely using nonlinear wave theory.According to the Stokes' second-order theory, the horizontal and vertical components are functions of wave period, height, and steepness [39].Wind wave steepness is a function of wave age and young age wind waves are much steeper, as described in [40], and a wind-wave spectrum model described in [41] shows how wave age affects the spectral peak and energy distribution, which is also used to retrieve the wind speed [42].Consequently, if wind waves show up in an amplitude spectrum, the wave age must be taken into consideration to improve the accuracy of the retrieval of wave period and height.
The amplitude spectrum measured at the real sea surface can be transformed to the wave spectrum in a specific direction.If many wave spectra for different directions are obtained, directional wave spectra can be established, which presents the wave energy distribution in the wave frequency and direction.The range of wave spectrum measurements is limited.It depends on the relative size of the dominant facet to the wave lengths.A capillary wave cannot be measured if its whole waveform is totally covered by the dominant facet, and a long period gravity wave cannot be measured if the component of Doppler shift frequency caused by it is less than the frequency resolution in the spectrogram of the IF signal.

Conclusions
In this paper, a measurement method of period, height and direction of ocean wave is presented.The short-range K-band NBCWR is employed as the wave sensor to measure orbital motions of ocean waves at low-grazing angle.The measurement is not based on the image generated by scanning sea surface but based on a dominant facet at the sea surface.
In the numerical simulations, we used the dominant waves from the Pierson-Moskowitz spectrum, the typical antenna pattern gains of the K-MC3 transceiver module and the K-band NRCS for Bragg backscattering.IF power distribution illustrates the existence of the dominant facet.It can be identified that it oscillates with wave propagations and its motions represent the orbital motions of the ocean waves.Wave period and height were obtained through measuring the period and amplitude of the variation of Doppler shift frequency caused by the dominant facet, and wave direction was deduced with the geometrical relationship of the measurement results in different directions.The multiple waves were measured as well.
In the wave tank validation experiments, we configured the radar with the same conditions as those for the numerical simulations.The generated wave is measured in two angles to the wave, 0 degrees and 50 degrees, respectively.The IF signals were strong enough to allow us to distinguish the Doppler shift frequencies and the strong signal in spectrograms shows the existence of the dominant facet.The measurement results were consistent with the characteristics of the generated wave.Therefore, the proposed method is reasonable and potentially feasible for ocean wave measurement.

Figure 1 .
Figure 1.Relationships between wave orbital velocity, ocean wave characteristics, and Doppler shift frequency.The wave measurement starts from measuring Doppler shift frequency with the radar.

Figure 1 .
Figure 1.Relationships between wave orbital velocity, ocean wave characteristics, and Doppler shift frequency.The wave measurement starts from measuring Doppler shift frequency with the radar.The variation of the measured Doppler shift frequency contains the information of the wave period, height, and direction.It is explained with the relationship γ 3 , which is deduced from γ 1 and γ 2 .

Figure 2 .
Figure 2. Velocities of the dominant facet  orbiting with a long surface wave.The radar is not shown here, but it transmits from right to left along the LOS and it is close to the sea surface. ,  and  are three projections of  to x-axis, z-axis, and LOS, respectively.Compared with  , the phase of  is delayed by  .

Figure 3 .
Figure 3. Geometrical relationship between the true wave height and three measured wave heights.All the directions are measured with respect to north.

Figure 3 .
Figure 3. Geometrical relationship between the true wave height and three measured wave heights.All the directions are measured with respect to north.

Figure 4 .
Figure 4. Top view of regular wave field, T = 5.85 s and L = 53.35m.The radar's location is (0, 0) and it transmits to the north.Wave direction is 140 degrees, measured clockwise from the north.

Figure 4 .
Figure 4. Top view of regular wave field, T = 5.85 s and L = 53.35m.The radar's location is (0, 0) and it transmits to the north.Wave direction is 140 degrees, measured clockwise from the north.

26 Figure 5 .
Figure 5. Antenna pattern of the K-MC3 transceiver module.The antenna's 3 dB beam widths in azimuth and elevation are seven degrees and 25 degrees, respectively.

Figure 5 .
Figure 5. Antenna pattern of the K-MC3 transceiver module.The antenna's 3 dB beam widths in azimuth and elevation are seven degrees and 25 degrees, respectively.

Figure 6 .
Figure 6.NRCS  for a VV polarized antenna (upwind).This is a calculated NRCS for a K-band radar under certain wind and sea water conditions as shown in Table2.
Remote Sens. 2018, 10, x FOR PEER REVIEW 12 of 26 west-east coordinates from −2 m to 2 m.The power level at dominant facet the  is much higher than for the side facets whose coordinates are −2 m, −1 m, 1 m, and 2 m.

Figure 8 .
Figure 8. Size variation of white area.The solid parts of the curve are visible to the radar but the dashed parts are invisible.

Figure 9 .
Figure 9. IF power distribution for the strong scattering area.The scattering power from the dominant facet located at (0, 7) is the strongest.

Figure 10 .
Figure 10.Power comparison around the dominant facet  .The absence of middle parts of the curves is due to the shadowed facets.The power levels at the dominant facet are much higher than the other four curves.Therefore, the scattering from side facets are negligible in the radar IF output.

Figure 8 .
Figure 8. Size variation of white area.The solid parts of the curve are visible to the radar but the dashed parts are invisible.

Figure 8 .
Figure 8. Size variation of white area.The solid parts of the curve are visible to the radar but the dashed parts are invisible.

Figure 9 .
Figure 9. IF power distribution for the strong scattering area.The scattering power from the dominant facet located at (0, 7) is the strongest.Figure 9. IF power distribution for the strong scattering area.The scattering power from the dominant facet located at (0, 7) is the strongest.

Figure 9 .
Figure 9. IF power distribution for the strong scattering area.The scattering power from the dominant facet located at (0, 7) is the strongest.Figure 9. IF power distribution for the strong scattering area.The scattering power from the dominant facet located at (0, 7) is the strongest.

Figure 10 .
Figure 10.Power comparison around the dominant facet  .The absence of middle parts of the curves is due to the shadowed facets.The power levels at the dominant facet are much higher than the other four curves.Therefore, the scattering from side facets are negligible in the radar IF output.

Figure 10 .
Figure 10.Power comparison around the dominant facet F 0 .The absence of middle parts of the curves is due to the shadowed facets.The power levels at the dominant facet are much higher than the other four curves.Therefore, the scattering from side facets are negligible in the radar IF output.

Figure 11 .
Figure 11.Comparison of Doppler shift frequency and signal power for dominant facet  and four adjacent facets.The power level of the adjacent facets decline rapidly so that the received power at the radar is dominated by the facet  .

Figure 12
Figure12shows the horizontal velocity and vertical velocity of the dominant facet  in five wave periods (29.25 s).The curves are sinusoids and oscillate with same amplitude, and the phase difference is π/2.Therefore, the motions of the dominant facet  represent the orbital motions of the generated ocean waves.

Figure 11 .
Figure 11.Comparison of Doppler shift frequency and signal power for dominant facet F 0 and four adjacent facets.The power level of the adjacent facets decline rapidly so that the received power at the radar is dominated by the facet F 0 .

Figure 11 .
Figure 11.Comparison of Doppler shift frequency and signal power for dominant facet  and four adjacent facets.The power level of the adjacent facets decline rapidly so that the received power at the radar is dominated by the facet  .

Figure 12 .
Figure 12.Horizontal velocity and vertical velocity of the dominant facet  .The velocities are calculated according to Equations (5) and (6), in which the variation of spatial position (x, z) of the dominant facet is also taken into account.

Figure 12 .
Figure 12.Horizontal velocity and vertical velocity of the dominant facet 0 .The velocities are calculated according to Equations (5) and (6), in which the variation of spatial position (x, z) of the dominant facet is also taken into account.
calculations of the power P IF , frequency f d , and phase ϕ of the IF component originated from each facet of the simulated ocean waves have been introduced.The total IF signal is the sum of the IF components and noise components.The IF signal and its spectrogram are shown in Figure 13.In the spectrogram, a bright line appears and is almost synchronized with the horizontal velocity of the dominant facet F 0 in Figure 12.This implies that the dominant facet F 0 dominates the received power from the simulated sea surface, and the bright line illustrates the variation of the Doppler shift frequency caused by the motions of the dominant facet F 0 .The bright line is extracted by finding the maximum value at each moment.The extracted bright line, i.e., the time series of the Doppler shift frequency f d , and its amplitude spectrum are plotted in Figure 14.A significant peak can be seen in the amplitude spectrum and the corresponding amplitude and frequency are 90.68Hz and 0.1706 Hz.Therefore, A d is 90.68 Hz and T d is 5.86 s.Substituting T d and A d into Equations (

Figure 13 .
Figure 13.IF signal and its spectrogram.The upper figure is the in-phase signal and the fluctuation of its amplitude reflects the wave propagation.The quadrature signal is not shown here.The in-phase signal and quadrature signal are the real part and the imaginary part of the complex IF signal, respectively.The strong signal in the lower figure originates from the dominant facet  .

Figure 13 .
Figure 13.IF signal and its spectrogram.The upper figure is the in-phase signal and the fluctuation of its amplitude reflects the wave propagation.The quadrature signal is not shown here.The in-phase signal and quadrature signal are the real part and the imaginary part of the complex IF signal, respectively.The strong signal in the lower figure originates from the dominant facet F 0 .

Figure 14 .
Figure 14.Extracted Doppler shift frequency and its amplitude spectrum.A harmonic of wave frequency can be found at the frequency 0.3412 Hz.

Figure 15 .
Figure 15.Spectrograms for two radar looking directions.Compared with the spectrogram in Figure 13, the Doppler shift frequency fluctuates with lower amplitude in the upper figure but higher amplitude in the lower figure.This is due to the coefficient |cos ()|.The wave direction is 140 degrees, and the radar looking directions are 0, 30, and 130 degrees.Therefore, the angles to wave:  = 40

Figure 14 .
Figure 14.Extracted Doppler shift frequency and its amplitude spectrum.A harmonic of wave frequency can be found at the frequency 0.3412 Hz.

26 Figure 14 .
Figure 14.Extracted Doppler shift frequency and its amplitude spectrum.A harmonic of wave frequency can be found at the frequency 0.3412 Hz.

Figure 15 .
Figure 15.Spectrograms for two radar looking directions.Compared with the spectrogram in Figure 13, the Doppler shift frequency fluctuates with lower amplitude in the upper figure but higher amplitude in the lower figure.This is due to the coefficient |cos ()|.The wave direction is 140 degrees, and the radar looking directions are 0, 30, and 130 degrees.Therefore, the angles to wave:  = 40

Figure 16 .
Figure 16.Two extracted Doppler shift frequencies and amplitude spectrums.The harmonics of wave frequency exist.

Figure 16 .
Figure 16.Two extracted Doppler shift frequencies and amplitude spectrums.The harmonics of wave frequency exist.

Figure 18 .
Figure 18.Two ocean waves.The highest points and the lowest points are 1.185 m and −1.185 m and located in the dark red areas and the dark blue areas, respectively.

Figure 18 .
Figure 18.Two ocean waves.The highest points and the lowest points are 1.185 m and −1.185 m and located in the dark red areas and the dark blue areas, respectively.

Figure 19 .
Figure 19.IF signal and its spectrogram for multiple ocean waves.In the lower figure, the Doppler shift frequency fluctuates with two frequencies corresponding to two wave frequencies 0.1709 Hz and 0.0912 Hz.

Figure 19 .
Figure 19.IF signal and its spectrogram for multiple ocean waves.In the lower figure, the Doppler shift frequency fluctuates with two frequencies corresponding to two wave frequencies 0.1709 Hz and 0.0912 Hz.
Figure 20 shows the extracted Doppler shift frequency and its amplitude spectrum.Two significant peaks stand out indicating the amplitudes and frequencies for the two waves.Remote Sens. 2018, 10, x FOR PEER REVIEW 19 of 26

Figure 20 .
Figure 20.Extracted Doppler shift frequency and its amplitude spectrum.The harmonics of the wave frequency exist.

Figure 20 .
Figure 20.Extracted Doppler shift frequency and its amplitude spectrum.The harmonics of the wave frequency exist.

Figure 21 .
Figure 21.Experimental set-up in wave tank.The radar is fixed on the stationary blue carriage.The radar transceiver is the white square in the center of the red circle.It can rotate to transmit in different directions to adjust the angle to wave, .

Figure 22 .
Figure 22.Wave profile in the wave tank.The horizontal straight line at 0 m is assumed to be still

Figure 21 .
Figure 21.Experimental set-up in wave tank.The radar is fixed on the stationary blue carriage.The radar transceiver is the white square in the center of the red circle.It can rotate to transmit in different directions to adjust the angle to wave, β.

Figure 23 .
Figure 23.Spectrograms for two measurement directions.The upper figure is the spectrogram which is obtained by pointing the radar transceiver to the wave coming direction, i.e.,  = 0 degree; the lower figure is obtained by rotating the radar transceiver counter-clockwise 50 degrees from the wave coming direction, i.e.,  = 50 degrees.

Figure 23 .
Figure 23.Spectrograms for two measurement directions.The upper figure is the spectrogram which is obtained by pointing the radar transceiver to the wave coming direction, i.e., β = 0 degree; the lower figure is obtained by rotating the radar transceiver counter-clockwise 50 degrees from the wave coming direction, i.e., β = 50 degrees.

Figure 23 .
Figure 23.Spectrograms for two measurement directions.The upper figure is the spectrogram which is obtained by pointing the radar transceiver to the wave coming direction, i.e.,  = 0 degree; the lower figure is obtained by rotating the radar transceiver counter-clockwise 50 degrees from the wave coming direction, i.e.,  = 50 degrees.

Table 1 .
Characteristics of waves at the spectral peak.

Table 1 .
Characteristics of waves at the spectral peak.

Table 3
shows the parameters of numerical simulation conditions about ocean wave, wave field, radar transceiver, and IF signal.

Table 4 .
Three simulation results.

Table 4 .
Three simulation results.

Table 5 .
Characteristics of two ocean waves.

Table 6 .
Simulation results for two waves.

Table 6 .
Simulation results for two waves.

Table 8 .
Wave tank measurement results.