# Study on Numerical Simulation and Electromagnetic Scattering of Time-Invariant Freak Waves

Authors: Lichen Han; Gengkun Wu; Bin Liu
Journal: International Journal of Antennas and Propagation
 (2022)

---

## Abstract

The occurrence probability of freak waves is related to the sea wave spectrum. In this paper, different wave spectrums are used to simulate time-invariant three-dimensional freak waves. Freak waves that meet the international standards are generated at fixed time and location by adjusting the energy of the wavelets. We studied the occurrence probability of freak waves under the conditions of different wave spectrums, different wind speeds, and different modulation ratios and optimized the calculation speed of the model. Simulation data show that the difference in the shape of the wave spectrum affects the probability of freak waves occurrence. The model conforms to the Benjamin–Feir index (BFI), and the ratio of wave steepness to spectrum bandwidth is the key. In this paper, the Kirchhoff approximation theory is used to study the electromagnetic scattering (EM scattering) properties of freak waves on the large scale. We ideally calculate the Normalized Radar Cross-Section (NRCS) from the sea surface with freak waves, under different wind speeds and different grazing angles. The NRCS of freak waves is extremely low, and the increase of wind speed and the decrease of the grazing angle will make the detection of freak waves more difficult. The possibility of detection of freak waves is higher at high grazing angles (low incidence angles). The numerical simulations provide reference for engineering.


---

## Body

## 1. Introduction

Freak waves are abnormally large-amplitude waves that appear suddenly from nowhere in different regions of the world ocean. This dangerous natural phenomenon unusually lasts for tens of seconds, and it is difficult to obtain complete and long-term observation data [1–3]. It is practically impossible for a ship’s crew to take any precaution measures to prevent freak waves with wave heights of 10∼30 meters. These waves can pose devastating damage on ships and marine structures and have a serious impact on the offshore industry [4–6]. Such extreme ocean phenomena have been extensively studied. The main contents of the theoretical part include definition, genesis, observation and record, numerical analysis, and simulation. In terms of application, the prevention of the hazard of freak waves is an important topic. Wave force under freak waves is a conventional and at the same time very frontier problem [7]. The early warning and forecasting of freak waves are also very important [8, 9].The mechanisms of freak waves generation are still under discussion. Freak waves can appear in extreme sea conditions and can suddenly appear in normal seas. Many mechanisms have been proposed to explain these rogue waves [10]. Wave focusing: Wave energy may be concentrated in a small area because of wave trapping, refraction, and reflection, such as wave–current interaction; wind force; bottom topography influence, etc. [11–13]. The modulation instability, also known as Benjamin-Feir instability [14], has been accepted by many researchers to explain the occurrence of freak waves. Benjamin-Feir instability has achieved good results in experimental studies [15, 16]. New mechanisms are constantly being proposed to explain this extreme phenomenon. For example, Fedele [17] questioned the modulational instability, arguing that freak waves are likely to be rare occurrences of weakly nonlinear random seas. Various studies on freak waves have been developing. For example, Lavrenov studied freak waves with deep trough [18]; L. Zou studied the rapid variation of water depth, such as slope and seamount, contributing to the occurrence probability of freak waves [19]; Takuji Waseda studied the relationship between freak waves and spectral [20]; Babanin studied the modulational instability, steepness, and spectral bandwidth [21, 22]. Exploring the generation mechanism of freak waves depends on the progress of experimental equipment and further research.Numerical simulation is an important part of the study of freak waves and has practical application value. Numerical simulation methods: the linear wave superposition theory; the nonlinear Schrödinger (NLS) equation; computational fluid dynamics (CFD) based on viscous flow theory; Korteweg–de Vries (KdV) equation, and Kadomtsev–Petviashvili (KP) equation, etc. [23]. The key point of this paper is to reduce the amount of calculation, improve the calculation efficiency, and quickly simulate the freak waves. The linear wave focusing method is chosen to simulate the freak waves.More attention has been paid in the study of EM scattering from rough sea surface. Obtaining marine information and target feature information through microwave remote sensing technology is meaningful for marine remote sensing, marine target detection, and identification, etc. For random rough sea surface, the strength of EM scattering ability (usually EM backscattering) has commonly used the scattering coefficient (or normalized radar cross-section, NRCS). There are many methods to calculate EM scattering from random sea surfaces, which can be divided into numerical methods and approximate methods. The former includes FDTD (Finite Difference Time Domain) [24], MoM (Method of Moment) [25], etc., while the latter includes KA (Kirchhoff Approximation) [26], SPM (Small Perturbation Method) [27], SSA (Small Slope Approximation) [28], TSM (Two-Scale Model) [29], etc.The sea surface is a time-variant random rough surface, which will lead to many difficulties. In fact, it is difficult for current observation instruments to distinguish between noise and freak waves in the case of extreme sea conditions, due to the influence of droplets, foam, inherent noise of radar, and so on. Observational data show that most of the freak waves have noise characteristics, which makes the recording and research of the freak waves more difficult. This requires the development of wave observation instruments. Numerical simulations can ideally study time-invariant freak waves, ignoring complex environmental factors. The numerical simulation of electromagnetic scattering from sea surface is not affected by sea clutter, and the calculation speed is fast. The disadvantage is that the accuracy is not enough, which is different from the engineering project. Although it cannot be directly applied to engineering projects, the preliminary study of EM scattering of freak waves in theory can provide reference for solving practical problems.In this paper, an ideal sea surface model with freak waves is established to study the EM scattering of freak waves. Freak waves are larger gravity waves. We study the approximate simulation of freak waves without considering the influence of surface tension on the spectrum, which can improve the calculation efficiency while maintaining the characteristics of freak waves. PM spectrum, JONSWAP spectrum, WEN spectrum, etc., can describe fully developed wind-driven sea surface and are suitable for simulating gravity waves. Based on the comparison of different wave spectra, JONSWAP wave spectrum and SWOP directional function are selected to simulate the random sea surface, and the freak waves are generated at fixed time and location on the sea surface. The freak waves simulated by JONSWAP spectrum is a sea surface with obvious macroscopic features. Small-scale waves ranging from a few centimeters to tens of centimeters and millimeter droplets have little effect on the numerical simulation. For the above reasons, without discussing the microwave scattering from the sea surface, this paper uses the Kirchhoff approximation method in the study of EM scattering from random sea surface to establish the EM backscattering model. We study the EM scattering from the freak waves under the condition of analytical approximation and compare the backscattering coefficient (NRCS) of freak waves and background waves.


## 2. Numerical Simulation of Freak Waves Based on Different Wave Spectra

The notion of freak waves was first introduced by Drapper (1965) [30], and extremely large-amplitude waves is its obvious feature. Klinting and Sand (1987) [31] proposed a definition of freak waves including three conditions:(1)
D1:H0/Hs≥2.(2)
D2:η/H0≥0.65.(3)
D3:H0/H0−1≥2,D4:H0/H0+1≥2.The wave heightH0 of a freak wave is 2 times greater than the significant wave height Hs. The height of the freak wave crest η is greater than or equal to 0.65 times the wave height of this freak wave. H0 is greater than or equal to 2 times the wave heights of two adjacent waves H0−1 and H0+1. The definition of freak waves includes wave height, the relationship between adjacent wave heights, and the relationship between crests and troughs. Accurate definition is the key to distinguish freak waves from common extreme waves. But there is no widely accepted new definition so far. This paper will generate freak waves that meet the above general definition D1∼D4.


### 2.1. Three Kinds of Sea Wave Spectra

In this paper, the wave height is obtained by(1)ζx,y,t=∑m=1M∑n=1NAmncosωmt−kmxcosθn−kmsinθn+φmn#.Amn is the amplitude of the wavelet, ωm is the angular frequency, km is the wave number, θn is the direction angle −π/2,π/2, and φmn is the random initial phase 0,2π.Amn=2Sω,θΔωΔθ. Sω,θ=SωGω,θ. Sω and Gω,θ are the sea spectrum and the directional spectrum, respectively. The sea wave spectrum describes the distribution of the energy inside the wave with respect to frequency and direction. Common wave spectra include Elfouhaily’s sea spectrum, JONSWAP spectrum, PM spectrum, WEN spectrum, etc. [32, 33]. The last three wave spectra are adapted to generate sea surface.JONSWAP spectrum describes waves that continue to develop in limited wind field:(2)Sjonω=αg2ω5exp−54ωjω4γrr=exp−ω−ωj22σ2ωj2#,where γ is the peak elevation factor. In this paper, γ=7 is chosen to represent the narrow-banded spectrum. α is the energy scale factor, ωp is the spectral peak circular frequency, and σ is the peak shape parameters. α=0.076U102/Fg0.22, ωj=22g2/U10F1/3, σ=0.07ω≤ωj, and 0.09ω>ωj. Fetch F is the distance from a lee shore, in this paper, F = 130000.PM spectrum describes a fully developed sea wave:(3)Spmω=αg2ω5exp−βωpω4#.ω=2πf, α=8.1×10−3, β=0.74, ωp=g/U19.5, U10, U19.5 is the wind speed at a place 10m, 19.5m from the sea surface. For most of the air currents over the ocean, U19.5≈1.026×U10.WEN spectrum is the wave spectrum of China offshore. It can describe the different stages of wind wave growth and adapt to different water depths:(4)Sω=m0pωwexp−95lnp5.813−5.137η6.77−1.088p+0.013p21.307−1.426ηωω0−112/50⩽ω⩽1.15ωwSω=m0ωw6.77−1.088p+0.013p21.307−1.426η5.813−5.137η1.15ωwωmω⩾1.15ωw#,where ω0 is spectrum peak frequency, m0 is zeroth moment, p is peakedness, and η is depth parameter. ωw=5.72/T1/3, m0=1/16H1/32, p=95.3H1/31.35/T1/32.7, η=0.625H1/3/d0<η<0.5, and m=22−η. H1/3 and T1/3 represent the effective wave height and effective period, respectively.In order to study the influence of different spectra on the occurrence probability of freak waves, this paper uses the SWOP (stereo wave observation project) directional spectrum uniformly:(5)Gω,θ=1π1+0.5+0.82exp−w/w042cos2θ+0.32exp−w/w042#,where w0=0.855×g/U10.WhenU10=7m/s, the sea surface is beginning to froth, while more than 30% of the sea surface area is covered with droplet and foam as U10=25m/s. In this paper, freak waves are generated under the condition of U10=13.8m/s to reduce the influence of foam on EM scattering calculation.According to (1), we can generate a huge wave height at a specific location by focusing the energy. Specific ranges of ω and θ are chosen to focus their full energy on raising the wave height. We generate freak waves according to Figure 1 and Figure 2. This paper selects some high-energy areas to focus energy according to the distribution of energy (Figure 2). In order to optimize the calculation, we need to choose an appropriate range, which will be discussed in Section 2.2. Before simulation, we need to test the model. In this paper, the sea surface is simulated, and the correctness of the model is verified by spectral estimation. We use the background waves to estimate the spectrum. Test the model by comparing the fit of the simulated spectrum to the target spectrum. The following is the result of the spectrum estimation:Figure 1
Three kinds of wave spectra and SWOP directional spectrum. (a) JONSWAP spectrum, PM spectrum, WEN spectrum. (b) SWOP directional spectrum.
(a)(b)Figure 2
The amplitude of the wavelet. Use the same SWOP spectrum. (a) JONSWAP spectrum. (b) PM spectrum. (c) WEN spectrum.
(a)(b)(c)Figure3 shows that the fitting effect is very good. We simulate freak waves based on the reliable sea surface model. Numerous studies have shown that the freak wave occurrence may be related to wave energy focusing. For some reasons, the wave energy focuses on a certain location, resulting in giant waves. In this paper, we adjust the energy of a part of the wavelets to generate the freak waves at a fixed time and location. It is difficult to judge freak waves and extreme waves by their shapes. Through the above definitions D1∼D4 to determine whether the simulated wave is a freak wave.Figure 3
We choose the background wave of the spatial sequence for spectral estimation. (a) A wave waveform selected from the background wave. (b) Spectrum estimation, solid line represents ideal spectrum (target spectrum), circle represents simulation data.
(a)(b)We conducted1000 simulations with different types of sea surface and randomly selected a set of data (Table 1). It can be seen from Figure 4 that the freak wave and the extreme wave are similar in shape. The definitions of freak waves are important. The wave height of the freak wave is not a decisive factor; its value can range from 7∼15m,U10=13.8m/s. The relationship between adjacent waves is the key to find freak waves. While extremum waves often appear in groups for a long time, the freak wave often appear suddenly alone. The occurrence probability of freak waves is much lower than extreme waves. We strictly follow the definitions when we model freak waves.Table 1
Three different types of sea surface. Record the definition of freak waveD1∼D4. Judging freak waves by parameters.
TypeSea waveExtreme waveFreak waveD11.582.852.70D20.580.900.75D31.281.753.79D42.121.032.75××√Figure 4
Three different types of sea surface are observed from the sectional shape. (a) Normal seawave. (b) Extreme wave with abnormal wave height. (c) Freak wave.
(a)(b)(c)


### 2.2. Numerical Simulation of Freak Wave

In this paper, the sea surface of800×800m is selected to simulate freak waves. Reduce the amount of calculation and improve the speed of the model while maintaining the wave spectrum structure. With a controlled focusing in both time and space, we decided to use the wave focusing approach, one of the fastest and powerful methods to simulate ideal freak waves. The angular frequency ω0.3,1.3 is divided into 50 parts 1,50, and the direction angle θ−π/2,π/2 is divided into 35 parts 1,35. To generate freak waves, we choose the range of wavelets used to focus the energy. Different wave spectra and different wind speeds affect the selection range of wavelets, and a reasonable range directly affects the occurrence of freak waves. The results are as follows.Freak waves are often accompanied by strong winds. We do not know whether the formation of freak waves is caused by the wind, and how the wind affects the generation of freak waves. Wind speed affects three sea wave spectra, and we study the occurrence probability of freak waves under different winds. Select medium to high wind speed13.8,17.1,20.7. Use SWOP direction function.θ10,18 used to focus energy. We have designed some sets of angular frequencies; the frequencies ω in the set are used to focus energy. The selection of elements in the set is the key to fast simulate the freak waves.The complete freak wave process is an important topic. The duration of the freak wave is very short. We simulated the whole process based on the observed freak wave data. According to the data, we simulated the occurrence process of freak waves for15s. The disadvantage is that the sea surface does not change with time, which is not a continuous process, for example, Figure 5. The independent time-invariant sea surface at different times helps us to ideally observe each stage of the freak wave, which is also helpful for the subsequent EM scattering calculation. We discuss the range of ω1,50. The results are as follows.Figure 5
Simulated sea surface800×800. (a) Normal sea surface. (b) Freak wave.
(a)(b)It can be seen from Table2 that JONSWAP spectrum is more suitable for simulation of freak waves and can adapt to various wind speeds compared with other sea spectra. The JONSWAP spectrum of γ=7 is a typical narrow band spectrum. It is reasonable to use narrow band spectrum for the simulation of freak waves [20]. Observe the shape of the narrow band spectrum; it narrows in both the frequency bandwidth and the directional spreading. The energy is concentrated at the wave crest, and it is easy to focus energy suddenly. Many studies have shown that the occurrence probability of freak waves is closely related to the shape of the wave spectrum [34]. The occurrence probability of freak waves is inhomogeneous in space. Researchers try to link freak wave events to particular meteorological conditions. A high probability freak wave region can be formed under certain synoptic meteorological conditions [20]. Prediction of the freak wave by using wave spectra is a possible and interesting filed [35]. Considering that the wave forecast systems are based on spectral modelling. The formation of narrow band wave spectrum may be a prerequisite for certain meteorological states. The relationship between meteorological states, wave spectrum, and freak waves is the focus of the study.Table 2
The occurrence probability of freak waves under different wind speeds, wave spectra, and wavelet ranges.
SpectralU103–27 (%)5–29 (%)7–31 (%)9–33 (%)11–35 (%)13–37 (%)15–39 (%)17–41 (%)19–43 (%)JONSWAP13.82183458708484868417.116324664647268686220.734364864646464584013.8829094948884724834PM17.182806456482462120.7564828124121013.8889492908080623216WEN17.184846038261022120.7543816812210


## 2.1. Three Kinds of Sea Wave Spectra

In this paper, the wave height is obtained by(1)ζx,y,t=∑m=1M∑n=1NAmncosωmt−kmxcosθn−kmsinθn+φmn#.Amn is the amplitude of the wavelet, ωm is the angular frequency, km is the wave number, θn is the direction angle −π/2,π/2, and φmn is the random initial phase 0,2π.Amn=2Sω,θΔωΔθ. Sω,θ=SωGω,θ. Sω and Gω,θ are the sea spectrum and the directional spectrum, respectively. The sea wave spectrum describes the distribution of the energy inside the wave with respect to frequency and direction. Common wave spectra include Elfouhaily’s sea spectrum, JONSWAP spectrum, PM spectrum, WEN spectrum, etc. [32, 33]. The last three wave spectra are adapted to generate sea surface.JONSWAP spectrum describes waves that continue to develop in limited wind field:(2)Sjonω=αg2ω5exp−54ωjω4γrr=exp−ω−ωj22σ2ωj2#,where γ is the peak elevation factor. In this paper, γ=7 is chosen to represent the narrow-banded spectrum. α is the energy scale factor, ωp is the spectral peak circular frequency, and σ is the peak shape parameters. α=0.076U102/Fg0.22, ωj=22g2/U10F1/3, σ=0.07ω≤ωj, and 0.09ω>ωj. Fetch F is the distance from a lee shore, in this paper, F = 130000.PM spectrum describes a fully developed sea wave:(3)Spmω=αg2ω5exp−βωpω4#.ω=2πf, α=8.1×10−3, β=0.74, ωp=g/U19.5, U10, U19.5 is the wind speed at a place 10m, 19.5m from the sea surface. For most of the air currents over the ocean, U19.5≈1.026×U10.WEN spectrum is the wave spectrum of China offshore. It can describe the different stages of wind wave growth and adapt to different water depths:(4)Sω=m0pωwexp−95lnp5.813−5.137η6.77−1.088p+0.013p21.307−1.426ηωω0−112/50⩽ω⩽1.15ωwSω=m0ωw6.77−1.088p+0.013p21.307−1.426η5.813−5.137η1.15ωwωmω⩾1.15ωw#,where ω0 is spectrum peak frequency, m0 is zeroth moment, p is peakedness, and η is depth parameter. ωw=5.72/T1/3, m0=1/16H1/32, p=95.3H1/31.35/T1/32.7, η=0.625H1/3/d0<η<0.5, and m=22−η. H1/3 and T1/3 represent the effective wave height and effective period, respectively.In order to study the influence of different spectra on the occurrence probability of freak waves, this paper uses the SWOP (stereo wave observation project) directional spectrum uniformly:(5)Gω,θ=1π1+0.5+0.82exp−w/w042cos2θ+0.32exp−w/w042#,where w0=0.855×g/U10.WhenU10=7m/s, the sea surface is beginning to froth, while more than 30% of the sea surface area is covered with droplet and foam as U10=25m/s. In this paper, freak waves are generated under the condition of U10=13.8m/s to reduce the influence of foam on EM scattering calculation.According to (1), we can generate a huge wave height at a specific location by focusing the energy. Specific ranges of ω and θ are chosen to focus their full energy on raising the wave height. We generate freak waves according to Figure 1 and Figure 2. This paper selects some high-energy areas to focus energy according to the distribution of energy (Figure 2). In order to optimize the calculation, we need to choose an appropriate range, which will be discussed in Section 2.2. Before simulation, we need to test the model. In this paper, the sea surface is simulated, and the correctness of the model is verified by spectral estimation. We use the background waves to estimate the spectrum. Test the model by comparing the fit of the simulated spectrum to the target spectrum. The following is the result of the spectrum estimation:Figure 1
Three kinds of wave spectra and SWOP directional spectrum. (a) JONSWAP spectrum, PM spectrum, WEN spectrum. (b) SWOP directional spectrum.
(a)(b)Figure 2
The amplitude of the wavelet. Use the same SWOP spectrum. (a) JONSWAP spectrum. (b) PM spectrum. (c) WEN spectrum.
(a)(b)(c)Figure3 shows that the fitting effect is very good. We simulate freak waves based on the reliable sea surface model. Numerous studies have shown that the freak wave occurrence may be related to wave energy focusing. For some reasons, the wave energy focuses on a certain location, resulting in giant waves. In this paper, we adjust the energy of a part of the wavelets to generate the freak waves at a fixed time and location. It is difficult to judge freak waves and extreme waves by their shapes. Through the above definitions D1∼D4 to determine whether the simulated wave is a freak wave.Figure 3
We choose the background wave of the spatial sequence for spectral estimation. (a) A wave waveform selected from the background wave. (b) Spectrum estimation, solid line represents ideal spectrum (target spectrum), circle represents simulation data.
(a)(b)We conducted1000 simulations with different types of sea surface and randomly selected a set of data (Table 1). It can be seen from Figure 4 that the freak wave and the extreme wave are similar in shape. The definitions of freak waves are important. The wave height of the freak wave is not a decisive factor; its value can range from 7∼15m,U10=13.8m/s. The relationship between adjacent waves is the key to find freak waves. While extremum waves often appear in groups for a long time, the freak wave often appear suddenly alone. The occurrence probability of freak waves is much lower than extreme waves. We strictly follow the definitions when we model freak waves.Table 1
Three different types of sea surface. Record the definition of freak waveD1∼D4. Judging freak waves by parameters.
TypeSea waveExtreme waveFreak waveD11.582.852.70D20.580.900.75D31.281.753.79D42.121.032.75××√Figure 4
Three different types of sea surface are observed from the sectional shape. (a) Normal seawave. (b) Extreme wave with abnormal wave height. (c) Freak wave.
(a)(b)(c)


## 2.2. Numerical Simulation of Freak Wave

In this paper, the sea surface of800×800m is selected to simulate freak waves. Reduce the amount of calculation and improve the speed of the model while maintaining the wave spectrum structure. With a controlled focusing in both time and space, we decided to use the wave focusing approach, one of the fastest and powerful methods to simulate ideal freak waves. The angular frequency ω0.3,1.3 is divided into 50 parts 1,50, and the direction angle θ−π/2,π/2 is divided into 35 parts 1,35. To generate freak waves, we choose the range of wavelets used to focus the energy. Different wave spectra and different wind speeds affect the selection range of wavelets, and a reasonable range directly affects the occurrence of freak waves. The results are as follows.Freak waves are often accompanied by strong winds. We do not know whether the formation of freak waves is caused by the wind, and how the wind affects the generation of freak waves. Wind speed affects three sea wave spectra, and we study the occurrence probability of freak waves under different winds. Select medium to high wind speed13.8,17.1,20.7. Use SWOP direction function.θ10,18 used to focus energy. We have designed some sets of angular frequencies; the frequencies ω in the set are used to focus energy. The selection of elements in the set is the key to fast simulate the freak waves.The complete freak wave process is an important topic. The duration of the freak wave is very short. We simulated the whole process based on the observed freak wave data. According to the data, we simulated the occurrence process of freak waves for15s. The disadvantage is that the sea surface does not change with time, which is not a continuous process, for example, Figure 5. The independent time-invariant sea surface at different times helps us to ideally observe each stage of the freak wave, which is also helpful for the subsequent EM scattering calculation. We discuss the range of ω1,50. The results are as follows.Figure 5
Simulated sea surface800×800. (a) Normal sea surface. (b) Freak wave.
(a)(b)It can be seen from Table2 that JONSWAP spectrum is more suitable for simulation of freak waves and can adapt to various wind speeds compared with other sea spectra. The JONSWAP spectrum of γ=7 is a typical narrow band spectrum. It is reasonable to use narrow band spectrum for the simulation of freak waves [20]. Observe the shape of the narrow band spectrum; it narrows in both the frequency bandwidth and the directional spreading. The energy is concentrated at the wave crest, and it is easy to focus energy suddenly. Many studies have shown that the occurrence probability of freak waves is closely related to the shape of the wave spectrum [34]. The occurrence probability of freak waves is inhomogeneous in space. Researchers try to link freak wave events to particular meteorological conditions. A high probability freak wave region can be formed under certain synoptic meteorological conditions [20]. Prediction of the freak wave by using wave spectra is a possible and interesting filed [35]. Considering that the wave forecast systems are based on spectral modelling. The formation of narrow band wave spectrum may be a prerequisite for certain meteorological states. The relationship between meteorological states, wave spectrum, and freak waves is the focus of the study.Table 2
The occurrence probability of freak waves under different wind speeds, wave spectra, and wavelet ranges.
SpectralU103–27 (%)5–29 (%)7–31 (%)9–33 (%)11–35 (%)13–37 (%)15–39 (%)17–41 (%)19–43 (%)JONSWAP13.82183458708484868417.116324664647268686220.734364864646464584013.8829094948884724834PM17.182806456482462120.7564828124121013.8889492908080623216WEN17.184846038261022120.7543816812210


## 3. Electromagnetic Scattering Calculation Model of Time-Invariant Freak Wave

The EM scattering characteristic of freak waves is a frontier field. In this paper, we study the EM scattering characteristics of freak waves from time-invariant sea surface and make a preliminary study on this problem under ideal conditions. We discuss this content under ideal conditions.Radar cross-section (RCS) or normalized radar cross-section (NRCS) represents the scattering ability of the target. It shows the backscattering characteristics of the target and in far field conditions can be expressed as(6)σ=limr⟶∞4πr2Es2AEi2#,where r is the distance between scattering body center and electromagnetic wave source. A denotes the illuminated area.Es and Ei are the incident field and the scattered field, respectively.For an ideal sea surface with the freak wave, we calculate its specular scattering. KA can be used to simulate EM backscattering from the large-scale sea surface. KA is broadly employed for studying EM scattering from large-scale sea surface due to its advantage of simple and practical KA method; the scattered field at any point within a source-free region bounded by a closed surface can be represented by a tangential field on the surface. The applicable range of KA is(7)kl>6Rc>λ#,where k is the incidence wave number, l is the correlation length, Rc is the radius of curvature, and λ is the incidence wavelength.According to the KA, ignoring edge effects such as diffraction, the scattered field is expressed as(8)Es=jkejkR0Ei4πR0∫−XX∫−YYa∂fx,y∂x+b∂fx,y∂y−cejv¯⋅r¯dxdy#,where a=1−Rsinθi+1+Rsinθscosφ and b=1+Rsinθssinφ.v⟶=ksinθi−sinθscosφx^−sinθssinφy^−cosθi+cosθsz^, θi is the incident angle, θs is the scattering angle, and φ is the angle between the projection of scattered rays and the x−axis.The reflection coefficientR is the Fresnel reflection coefficient, and, in the case of backscattering, the cross-polarization coefficient is 0. The expression is as follows:(9)Rvv=Ycosθi−Ycosθi−sin2θiYcosθi+Ycosθi−sin2θiRhh=cosθi−Ycosθi−sin2θicosθi+Ycosθi−sin2θi#,where Y is the relative permittivity of seawater.Y=81.5 for normal sea water. Substituting (9) into (8), the EM backscattering coefficient is(10)σ=−14XYcosθic+avxvz∫−XX∫−YYejvxx+vzfx,ydxdy#KA is used to calculate the EM backscattering coefficient from the time-invariant 3D freak waves. C-band SAR is one of the most widely used in ocean observation. This article takes C-band as an example. The freak waves are easier to detect from the background sea surface with HH polarization than with VV polarization. Set simulation parameters:U10=13.8m/s,20.7m/s; the frequency of incident waves is 5.3 GHz; HH polarization; grazing angles θg are 40°,80°. The results are as follows.We can ideally observe the various stages of the freak wave, like Figure6, and calculate each NRCS. Different colors on Figure 7 represent electromagnetic scattering of different intensities. Blue indicates that the backscattered field in this area is relatively weak and is shielded by adjacent waves. Red indicates the backscattered field strength. We compared the NRCS from freak waves under different conditions. Wind speed U10 and grazing angle θg have great influence on the numerical simulation. We conducted simulations with different parameters, as follows:Figure 6
The process of freak waves. Characteristics of each stage. These stages are discontinuous. (a)∼(c) Different stages of time-invariant freak waves simulated by JONSWAP spectrum. (d)∼(f) Simulated by PM spectrum. (g)∼(i) Simulated by WEN spectrum.
(a)(b)(c)(d)(e)(f)(g)(h)(i)The above figures show the two stages before and after the occurrence of the freak wave. We use different grazing angles to observe. We can see the NRCS difference clearly from Figure8. T=2s is the stage when the freak wave is forming. At this stage, the NRCS difference between the freak wave and the background wave is little. It is difficult to capture the freak wave from its initial stage. T=15s is the stage when the freak wave disappears, and, at this time, there is no obvious difference from the normal sea surface. It is also difficult to capture freak waves from the top view Figure 9. The following is the NRCS under different grazing angles and different wind speeds:Figure 7
Backscattering (NRCS) from sea surface with the freak wave.U10=13.8m/s, θg=80°.Figure 8
X-Z view of NRCS from sea surface with different times andθ. We can clearly see the NRCS difference between the freak wave and the background wave. (a) Preliminary stage of freak wave θ=40°. (b) End stage of freak wave θ=40°. (c) Preliminary stage θ=80°. (d) End stage θ=80°.
(a)(b)(c)(d)Figure 9
X–Y view of NRCS distribution from sea surface with different times andθ. We can see NRCS from the whole sea surface clearly. (a) Preliminary stage of freak wave θ=40°. (b) End stage of freak wave θ=40°. (c) Preliminary stage θ=80°. (d) End stage θ=80°.
(a)(b)(c)(d)According to the simulation results, we can directly see the EM scattering from the freak waves. There are obvious differences between the Figures8–11. The EM backscattering (NRCS) from freak waves is significantly different from the background waves (Figure 10). We conducted numerical simulations with different grazing angles and different wind speeds. Observing the change of grazing angle, the NRCS difference between the freak wave and the background wave increases with the increase of the grazing angle. The blue part in the Figure 11 obviously increases, but the shape of the red range indicating the freak wave becomes smaller. With the increase of wind speed, the wave height of the background wave increases, and it becomes difficult to capture the freak wave. The red part under the low grazing angle increased significantly, and the freak wave can be clearly distinguished until the grazing angle is increased. High grazing angles (low incidence angles) make it more possible to detect freak waves.Figure 10
X-Z view of NRCS from freak wave sea surface with differentU10 and θ. (a) U10=13.8m/s, θ=40°. (b) U10=20.7m/s, θ=40°. (c) U10=13.8m/s, θ=80°. (d) U10=20.7m/s, θ=80°.
(a)(b)(c)(d)Figure 11
X–Y view of NRCS distribution from freak wave sea surface with differentU10 and θ. (a) U10=13.8m/s, θ=40°. (b) U10=20.7m/s, θ=40°. (c) U10=13.8m/s, θ=80°. (d) U10=20.7m/s, θ=80°.
(a)(b)(c)(d)Freak waves are difficult to capture in practical radar because of its noise characteristics. The real radar image is always very noisy due to the speckle noise. Freak waves are very localized phenomena, of the order of SAR image pixel size. On high resolution radar, the imaging of the freak wave is not obvious, and it is difficult to distinguish from noise. Thus, how this speckle noise (about10∼20dB in magnitude) can be distinguished from freak waves is a research topic. We can obtain a value according to the obvious difference between the freak wave and the background wave, like Figures 10 and 11, which is used to distinguish the background wave, noise, and freak wave. According to numerous simulations, we suggest that the NRCS threshold value of the freak wave is 25∼50dB. Many factors can affect the NRCS threshold value of freak waves, such as significant wave heights in the area, ocean currents, synoptic meteorological, etc. Even the radar itself can affect the threshold value. Dynamically calculating the NRCS threshold value of freak waves is a challenge. The study of EM scattering from sea surface with freak waves is a frontier subject.


## 4. Conclusions

The ‘freak waves’ means individual waves with a crest of an extremely high slope compared with background waves. According to the observed freak wave data, we simulated the whole process. This paper discusses the influence of different wave spectra on the simulation of freak waves under different wind speeds. On the validated numerical simulation model of sea surface, the time-invariant 3D freak wave is generated by adjusting the wavelet energy. The numerical simulation shows that the occurrence probability of simulated freak waves is significantly different with different wave spectra under different conditions. The shape of the spectrum affects the occurrence probability of freak waves. The narrow band spectrum is more suitable for the simulation of freak waves. We usually do wave forecast based on spectral modelling, which helps researchers to study the predictability of rogue events by means of probabilistic approaches.In this paper, we choose JONSWAP spectrum to simulate the freak wave for EM scattering calculation. The EM scattering characteristics of sea surface with freak waves are studied by KA method. The EM backscattering (NRCS) from freak wave and background wave is studied under different grazing angles and different wind speeds. The numerical simulations show that the backscattering coefficients of the freak wave and the background wave are significantly different. The wind speed and the grazing angle have great influence on NRCS. In real observation, it is difficult to find freak waves by definingD1∼D4. Radar echoes can catch them faster and more accurately. We try to use the difference between the NRCS of the background wave and the freak wave as the judgment condition. This theory is feasible under ideal conditions. However, in practical engineering, it is still difficult to distinguish the freak wave from the speckle noise, which makes the identification and capture of the freak wave more difficult. Accurately detecting freak waves on real radar and distinguishing noise is the next topic.



---

*Source: 1951162-2022-09-20.xml*