Although the nonlinear interactions between a single short gravity wave and a long wave can be solved analytically, the solution is less tractable in more general cases involving multiple short waves. In this work we present a numerical method of studying nonlinear interactions between a long wave and multiple short harmonic waves in infinitely deep water. Specifically, this method is applied to the calculation of the temporal and spatial evolutions of the surface elevations in which a given long wave interacts with several short harmonic waves. Another important application of our method is to quantitatively analyse the nonlinear interactions between an arbitrary short wave train and another short wave train. From simulation results, we obtain that the mechanism for the nonlinear interactions between one short wave train and another short wave train （expressed as wave train 2） leads to the energy focusing of the other short wave train （expressed as wave train 3）. This mechanism occurs on wave components with a narrow frequency bandwidth, whose frequencies are near that of wave train 3.
To study the electromagnetic （EM） backscatter characteristics of freak waves at moderate incidence angles, we establish an EM backscattering model for freak waves in （1＋1）-dimensional deep water. The nonlinear interaction between freak waves and Bragg short waves is considered to be the basic hydrodynamic spectra modulation mechanism in the model. Numerical results suggest that the EM backscattering intensities of freak waves are less than those from the background sea surface at moderate incidence angles. The normalised radar cross sections （NRCSs） from freak waves are highly polarisation dependent, even at low incidence angles, which is different from the situation for normal sea waves; moreover, the NRCS of freak waves is more polarisation dependent than the background sea surface. NRCS discrepancies between freak waves and the background sea surface with using horizontal transmitting horizomtal （HH） polarisation are larger than those using vertical transmitting vertical （VV） polarisation, at moderate incident angles. NRCS discrepancies between freak waves and background sea surface decreases with the increase of incidence angle, in both HH and VV polarisation radars. As an application, in the synthetic-aperture radar （SAR） imaging of freak waves, we suggest that freak waves should have extremely low backscatter NRCSs for the freak wave facet with the strongest slope. Compared with the background sea surface, the freak waves should be darker in HH polarisation echo images than in VV echo images, in SAR images. Freak waves can be more easily detected from the background sea surface in HH polarisation images than in VV polarisation images. The possibility of detection of freak waves at low incidence angles is much higher than at high incidence angles.
Using the theory of nonlinear interactions between long and short waves, a nonlinear fractal sea surface model is presented for a one dimensional deep sea. Numerical simulation results show that spectra intensity changes at different locations （in both the wave number domain and temporal-frequency domain）, and the system obeys the energy conservation principle. Finally, a method to limit the fractal parameters is also presented to ensure that the model system does not become ill-posed,
