The effect of wave breaking on surface wave imaging by Synthetic Aperture Radar

The effect of ocean wave breaking as a non-Bragg mechanism on backscattering cross-section and modulation transfer functions （MTF） of radar was investigated based on Bragg resonance theory and parametric method. The result showed that the additional effect of wave breaking on backscattering cross-section is not more than 20% except for the small incident angle of VV polarized electromagnetic （e.m.） wave but is significant for HH polarized e.m. wave. Breaking waves lead to increase in the modulus of tilt modulation MTF and the larger the wind speed, the faster the increase. For large incident angle, the modulus of tilt modulation MTF with wave breaking decreases quickly with incident angle for HH polarization and approach to that without wave breaking for VV polarization. The hydrodynamic MTF increases 30%-60% when considering wave breaking and the increase is larger for HH polarization than for VV polarization.     

A wave modulation model of ripples over long surface waves

A study is presented on the modulation of ripples induced by a long surface wave (LW) and a new theoretical modulation model is proposed. In this model, the wind surface stress modulation is related to the modulation of ripple spectrum. The model results show that in the case of LW propagating in the wind direction with the wave age parameter of LW increasing, the area with enhanced shear stress shifts from the region near the LW crest on the upwind slope to the LW trough. With a smaller wave age parameter of LW, the ripple modulation has the maximum on the upwind slope in the vicinity of LW crest, while with a larger parameter the enhancement of ripple spectrum does not occur in that region. At low winds the amplitude of ripple modulation transfer function (MTF) is larger in the gravity wave range, while at moderate or high winds it changes little in the range from short gravity waves to capillary waves.     

Measurements of ocean wave and current field using dual polarized X-band radar

A new ocean wave and sea surface current monitoring system with horizontally-(HH) and vertically-(VV) polarized X-band radar was developed. Two experiments into the use of the radar system were carried out at two sites, respectively, for calibration process in Zhangzi Island of the Yellow Sea, and for validation in the Yellow Sea and South China Sea. Ocean wave parameters and sea surface current velocities were retrieved from the dual polarized radar image sequences based on an inverse method. The results obtained from dual-polarized radar data sets acquired in Zhangzi Island are compared with those from an ocean directional buoy. The results show that ocean wave parameters and sea surface current velocities retrieved from radar image sets are in a good agreement with those observed by the buoy. In particular, it has been found that the vertically-polarized radar is better than the horizontally-polarized radar in retrieving ocean wave parameters, especially in detecting the significant wave height below 1.0 m.     

基于GF-3和Landsat 8遥感数据的土壤水分反演研究

基于我国首颗全极化雷达卫星高分三号（GF-3）和Landsat8数据,研究浓密植被覆盖地表土壤水分反演方法。为了提高浓密植被覆盖地表土壤水分反演精度,首先利用PROSAIL模型、实测植被参数及Landsat8光学数据分析了8种植被指数与植被冠层含水量的相关性,从中优选出归一化差异水指数（NDWI₅）用于反演植被冠层含水量,并通过分析植被含水量和植被冠层含水量的关系,构建植被含水量模型;然后结合植被含水量反演模型和简化MIMICS模型校正了植被对雷达后向散射系数的影响,最后基于AIEM建立裸土后向散射系数模拟数据集,发展一种主动微波和光学数据协同反演浓密植被覆盖地表土壤水分模型,并以山东省禹城市为研究区,实现了玉米覆盖下HH、VV和HH+VV 3种模式土壤水分反演。实验结果表明： ① NDWI₅为最佳植被指数,对于去除植被影响有较好效果;② 基于此方法,利用GF-3和Landsat8卫星数据反演得到的土壤水分具有较高的精度;③ 相比HH和VV两种极化模式,HH+VV双通道模式对土壤水分反演结果更好,决定系数（R2）为0.4037,均方根误差（RMSE）为0.0667 m ³m ^-3。     

A new modulation transfer function for ocean wave spectra retrieval from X-band marine radar imagery

When imaging ocean surface waves by X-band marine radar,the radar backscatter from the sea surface is modulated by the long surface gravity waves. The modulation transfer function(MTF) comprises tilt,hydrodynamic,and shadowing modulations. A conventional linear MTF was derived using HH-polarized radar observations under conditions of deep water. In this study,we propose a new quadratic polynomial MTF based on VV-polarized radar measurements taken from heterogeneous nearshore wave fields. This new MTF is obtained using a radar-observed image spectrum and in situ buoy-measured wave frequency spectrum. We validate the MTF by comparing peak and mean wave periods retrieved from X-band marine radar image sequences with those measured by the buoy. It is shown that the retrieval accuracies of peak and mean wave periods of the new MTF are better than the conventional MTF. The results also show that the bias and root mean square errors of the peak and mean wave periods of the new MTF are 0.05 and 0.88 s,and 0.32 and 0.53 s,respectively,while those of the conventional MTF are 0.61 and 0.98 s,and 1.39 and 1.48 s,respectively. Moreover,it is also shown that the retrieval results are insensitive to the coefficients in the proposed MTF.     

ENVISAT ASAR数据与CBERS-02数据的融合与应用

本文首次利用ENVISAT ASAR交叉极化模式数据与中巴地球资源卫星02星(CBERS 02)CCD数据进行融合,对融合后的图像进行分析和讨论。通过对IHS变换融合、主成分分析融合和小波变换融合比较分析发现,采用IHS变换融合具有更好的目视判读效果,但是由于ASAR数据受到相干斑噪声的影响,在滤波去噪后,使融合图的空间分辨率变低。另外,该融合图包含了地物的VV极化和HH极化后向散射特征,使地物结构形态、复介电常数和地物的光谱信息等重要的地物参数信息都包含在融合图中,增强了区别地物的能力,就农作物的目视解译和分类而言,精度得到了较大的提高。     

Retrieval of Significant Wave Height Under Typhoon Conditions from Gaofen-3 SAR Imagery

The objective of this paper is to propose an empirical method to inverse significant wave height(SWH)under typhoon conditions from collected dual-polarization Gaofen(GF)-3 synthetic aperture radar(SAR)imagery.The typhoon scenes were cap-tured from narrow scan(NSC)and wide scan(WSC)images,and collocated with European Center for Medium-Range Weather Fore-casts reanalysis data of(ECMWF).To improve the quality of GF-3 SAR images,the recalibration over rainforest and de-scalloping were carried out.To establish the empirical relationship between SAR-derived parameters and collocated SWH,the sensitivity analysis of typical parameters about the normalized radar cross section(Nrcs)and imagery variance(Cvar)were performed to both VV and VH polarized images.Four scenes from GF-3 SAR imagery under typhoon conditions were used for training the model by the multivari-ate least square regression,and one scene was used for preliminary validation.It was found that the joint retrieval model based on VV and VH polarized SAR imagery performed better than any single polarized model.These results,verified by using ECMWF data,revealed the soundness of this approach,with a correlation of 0.95,bias of 0 m,RMSE of 0.44 and SI of 0.01 when VV polarization and VH polarization data were both used.     

A new model to estimate significant wave heights with ERS-1/2 scatterometer data

A new model is proposed to estimate the significant wave heights with ERS-1/2 scatterometer data. The results show that the relationship between wave parameters and radar backscattering cross section is similar to that between wind and the radar backscattering cross section. Therefore, the relationship between significant wave height and the radar backscattering cross section is established with a neural network algorithm, which is, if the average wave period is ≤7s, the root mean square of significant wave...     

Evaluation of Wind Retrieval from Co-Polarization Gaofen-3 SAR Imagery Around China Seas

Gaofen-3(GF-3) is the first Chinese space-borne satellite to carry the C-band multi-polarization synthetic aperture radar(SAR). Marine applications, i.e., winds and waves retrieved from GF-3 SAR images, have been operational since January 2017. In this study, we have collected more than 1000 quad-polarization(vertical-vertical(VV); horizontal-horizontal(HH); vertical-horizontal(VH); horizontal-vertical(HV)) GF-3 SAR images, which were acquired around the China Seas from September 2016 to September 2017. Wind streaks were visible in these images in co-polarization(VV and HH) channel. Geophysical model functions(GMFs), including the CMOD5N together with polarization ratio(PR) model and C-SARMOD, were used to retrieve winds from the collected co-polarization GF-3 SAR images. Wind directions were directly obtained from GF-3 SAR images. Then, the SAR-derived wind speeds were compared with the measurements at a 0.25? grid from the Advanced Scatterometer on board the Metop-A/B and microwave radiometer WindSAT. Based on the analysis, empirical corrections are proposed to improve the performance of the two GMFs. Results of this study show that the standard deviation of wind speed is 1.63 m s~(-1) with a 0.19 m s~(-1) bias and 1.71 m s~(-1) with a 0.26 m s~(-1) bias for VV-and HH-polarization GF-3 SAR, respectively. Our work not only systematically evaluates wind retrieval by using the two advanced GMFs and PR models but also proposes empirical corrections to improve the accuracy of wind retrievals from GF-3 SAR images around the China Seas and thus enhance the accuracy of near real-time operational SAR-derived wind products.     

Extracting the subsonic anti-symmetric lamb wave from a submerged thin spherical shell backscattering through iterative time reversal

The extraction of the weakly excited anti-symmetric Lamb wave from a submerged thin spherical shell backscattering is very difficult if the carrier frequency of the incident short tone burst is not at its frequency of greatest enhancement. Based on a single channel iterative time reversal technique, a method for isolating the subsonic anti-symmetric Lamb wave is proposed in this paper. The approach does not depend on the form function of a thin shell and any other priori knowledge, and it is also robust in the presence of some stochastic noise. Both theoretical and numerical results show that the subsonic anti-symmetric Lamb wave can be identified, even when the carrier frequency of the incident short tone burst is away from the frequency of greatest enhancement. The phenomenon may also be observed even in the case that the subsonic anti-symmetric Lamb wave is submerged in the noise, other than the case with the Signal to Noise Ratio being less than 10 dB, when the amplitude of the noise is comparable with the specular wave. In this paper, each iteration process contains a traditional transmission and time reversal transmission steps. The two steps can automatically compensate the time delay of the subsonic anti-symmetric Lamb wave relative to the specular wave and within-mode dispersion in the forward wave propagation.     

Underwater topography detection of Taiwan Shoal with SAR images

Under suitable conditions of tidal current and wind, underwater topography can be detected by synthetic aperture radar (SAR) indirectly. Underwater topography SAR imaging includes three physical processes: radar ocean surface backscattering, the modulation of sea surface short wave spectrum by the variations in sea surface currents, and the modulation of sea surface currents by the underwater topography. The first process is described usually by Bragg scattering theory because the incident angle of SAR is always between 20°–70°. The second process is described by the action balance equation. The third process is described by an ocean hydrodynamic model. Based on the SAR imaging mechanism for underwater topography, an underwater topography SAR detection model and a simplified method for its calculation are introduced. In the detection model, a two-dimensional hydrodynamic model — the shallow water model is used to describe the motion of tidal current. Due to the difficulty of determining the expression of SAR backscattering cross section in which some terms can not be determined, the backscattering cross section of SAR image used in the underwater topography SAR detection is pro-processed by the simulated SAR image of the coarse-grid water depth to simplify the calculation. Taiwan Shoal, located at the southwest outlet of Taiwan Strait, is selected as an evaluation area for this technique due to the occurrence of hundreds of sand waves. The underwater topography of Taiwan Shoal was detected by two scenes of ERS-2 SAR images which were acquired on 9 January 2000 and 6 June 2004. The detection results are compared with in situ measured water depths for three profiles. The average absolute and relative errors of the best detection result are 2.23 m and 7.5 %, respectively. These show that the detection model and the simplified method introduced in the paper is feasible.     

微波与光学遥感协同反演旱区地表土壤水分研究

土壤水分是水文循环中的关键因素,尤其对旱区的生态环境具有十分重要的意义。微波遥感是反演土壤水分的有效手段,而植被是影响土壤水分反演精度的重要因素。因此,对土壤水分的反演需要考虑植被的影响。本文以内蒙古乌审旗为研究区,利用Radarsat-2雷达数据与TM光学数据,对旱区稀疏植被覆盖地表土壤水分反演进行研究。利用TM数据,分别选取NDVI和NDWI指数对植被含水量进行反演,通过水云模型消除植被层对土壤后向散射系数的影响;在此基础上,根据研究区地表植被特性,提出一种基于AIEM 模型的反演土壤水分的改进算法,反演了不同粗糙度参数、不同极化（VV极化和HH极化）条件下的研究区土壤水分。反演结果与野外实测数据的对比结果表明,本文提出的基于地表植被特性的土壤水分改进算法,具有更好的适应性;土壤水分反演模式\begin{array}{c}M{v}_{\mathit{\sigma vv}\left(1\right)}^{\mathit{lh}}\end{array}（VV极化方式下采用NDVI去除植被影响的反演模式）更适合于旱区考虑稀疏植被覆盖影响的地表土壤水分的反演。     

CALCULATIONS OF WAVE BREAKING PROBABILITY AND DISTRIBUTION OF BREAKING WAVE HEIGHTS IN DEEP WATER

lareODUcrI0NThephenomenonofbreakingwavesintheoasnoocurswheneveramomentahiyhighcrestmecheSanunstablecondition.It0ccursintendtimhy,andtheoccurrencefre-qUencydePendsonthescastate.ManystudAshaddricaniedoutontheoccurmcefre-qUencyofbrmkingwavesindeepweter.0chiandTsai(l983),LongueHiggins(l975)andVanD0mandPamn(l975)nadstheptalictonofbrmkingwavesoccurrenceindeePwater,usingabaskingcriterionforindividualwavesbasedonthewavesmpness.ThecriterioncanbeexpmeedintennSofthewaVehdghtHandperiodTasH>PgT'(…     

基于SVR的旱区稀疏植被覆盖下土壤水分遥感反演

中国西北半干旱区降水稀少、蒸散强烈,土壤水分作为重要的生态因子,影响着土壤-大气界面的能量平衡。支持向量回归模型具有估算精度高、可处理非线性问题、泛化能力强等优点,近年来被应用于土壤水分反演研究中,但已有模型极少考虑地表粗糙度因素的影响,导致反演精度受到一定限制。因此,本文以内蒙古乌审旗为研究区,采用水云模型去除地表稀疏植被覆盖的影响,提取全极化Radarsat-2 SAR影像裸土后向散射系数（ {\sigma }_{\mathrm{soil}}^0 ）,并利用AIEM模型和Oh模型建立后向散射系数数据库,采用LUT法模拟地表有效粗糙度参数,构建基于支持向量回归的土壤水分反演模型,并系统地对比分析了不同极化方式的后向散射系数作为数据源的土壤水分反演结果。研究结果表明：不考虑粗糙度参数的单数据源作为模型参数时,同极化数据反演结果比交叉极化具有更高的反演精度;当模型参数为考虑粗糙度的多源数据时,不同极化数据的反演精度均有所提高,其中数据源为 {\sigma }_{\mathit{vv}}^0 和粗糙度参数时,反演结果最好（R ²=0.917,MAE=3.980%,RMSE=5.187%）。研究结果可为旱区稀疏植被覆盖地表土壤水分的遥感监测提供技术支持。     

Determination of fractional energy loss of waves in nearshore waters using an improved high-order Boussinesq-type model

Fractional energy losses of waves due to wave breaking when passing over a submerged bar are studied systematically using a modified numerical code that is based on the high-order Boussinesq-type equations. The model is first tested by the additional experimental data, and the model’s capability of simulating the wave transformation over both gentle slope and steep slope is demonstrated. Then, the model’s breaking index is replaced and tested. The new breaking index, which is optimized from the several breaking indices, is not sensitive to the spatial grid length and includes the bottom slopes. Numerical tests show that the modified model with the new breaking index is more stable and efficient for the shallow-water wave breaking. Finally, the modified model is used to study the fractional energy losses for the regular waves propagating and breaking over a submerged bar. Our results have revealed that how the nonlinearity and the dispersion of the incident waves as well as the dimensionless bar height (normalized by water depth) dominate the fractional energy losses. It is also found that the bar slope (limited to gentle slopes that less than 1:10) and the dimensionless bar length (normalized by incident wave length) have negligible effects on the fractional energy losses.     

Dynamic response analysis of a floating mooring system

An innovative floating mooring system with two or more independent floating mooring platforms in the middle and one rigid platform on each side is proposed for improving efficiency and safety in shallow water. For this new system, most of collision energy is absorbed through the displacement of floating platforms. In order to illustrate the validity of the system, a series of model tests were conducted at a scale of 1:40. The coupled motion characteristics of the floating mooring platforms were discussed under regular and irregular waves, and the influences of wave direction and other characteristics on dynamic response of the system were analyzed. The results show that the mooring system is safest at 0° of wave incident angle, whereas the most dangerous mooring state occurs at 90° of wave incident angle. Motion responses increase with the increase of wave height, but are not linearly related to changes in wave height.     

Wave Pressure Acting on V-Shaped Floating Breakwater in Random Seas

Wave pressure on the wet surface of a V-shaped floating breakwater in random seas is investigated. Considering the diffraction effect, the unit velocity potential caused by the single regular waves around the breakwater is solved using the finite-depth Green function and boundary element method, in which the Green function is solved by integral method. The Response-Amplitude Operator(RAO) of wave pressure is acquired according to the Longuet-Higgins' wave model and the linear Bernoulli equation. Furthermore, the wave pressure's response spectrum is calculated according to the wave spectrum by discretizing the frequency domain. The wave pressure's characteristic value corresponding to certain cumulative probability is determined according to the Rayleigh distribution of wave heights. The numerical results and field test results are compared, which indicates that the wave pressure calculated in random seas agrees with that of field measurements. It is found that the bigger angle between legs will cause the bigger pressure response, while the increase in leg length does not influence the pressure significantly. The pressure at the side of head sea is larger than that of back waves. When the incident wave angle changes from 0? to 90?, the pressure at the side of back waves decreases clearly, while at the side of head sea, the situation is more complicated and there seems no obvious tendency. The concentration of wave energy around low frequency(long wavelength) will induce bigger wave pressure, and more attention should be paid to this situation for the structure safety.     

Wave steepness retrieved from scatterometer data in a genetic algorithm

Wave steepness is an important characteristic of a high sea state, and is widely applied on wave propagations at ports, ships, offshore platforms, and CO2 circulation in the ocean. Obtaining wave steepness is a difficult task that depends heavily on theoretical research on wavelength distribution and direct observations. Development of remote-sensing techniques provides new opportunities to study wave steepness. At present, two formulas are proposed to estimate wave steepness from QuikSCAT and ERS-1/2 scatterometer data. We found that wave steepness retrieving is not affected by radar band, and polarization method, and that relationship of wave steepness with radar backscattering cross section is similar to that with wind. Therefore, we adopted and modified a genetic algorithm for relating wave steepness with radar backscattering cross section. Results show that the root-mean-square error of the wave steepness retrieved is 0.005 in two cases from ERS-1/2 scatterometer data and from QuikSCAT scatterometer data.     

A mechanism of the effect of non-uniform current on the spectrum of short wind waves

A mechanism is suggested in this paper concerning the effect of non-uniform current on the spectrum of short wind waves. According to this mechanism, a non-uniform current brings changes to the breaking criteria of short wind waves through modulating the surface drift, and hence enhances or weakens wave breaking. Some modification is proposed to the source term, which represents the spectral rate of wave energy dissipation due to wave breaking so that the source term can incorporate this mechanism. In order to illustrate whether this mechanism is significant, a real case is studied, in which the wind waves propagate on a tidal current flowing over the sea bottom covered with sand waves. Finally, the effect of the new mechanism on the equilibrium spectrum of small scale gravity waves is discussed. Numerical estimates suggest that, for water depths less than 50 m and wavelengths less than 1 m, this current field may result in distinct spatial variations of the wave breaking criteria, the spectral rate of wave energy dissipation and the equilibrium spectrum of short gravity waves.     

Laboratory Experiments on an Internal Solitary Wave over a Triangular Barrier

Laboratory experiments were conducted to investigate the evolution of interfacial internal solitary waves(ISWs) incident on a triangular barrier. ISWs with different amplitudes were generated by gravitational collapse. The ISW energy dissipation and turbulence processes were calculated as waves passed over the triangular barrier. Experimental results showed that ISWs were reflecting back off the triangular barrier, and shoaling ISWs led to wave breaking and mixing when waves propagated over the obstacle. Wave instability created the dissipation of energy as it was transmitted from waves to turbulence. The rate of ISW energy dissipation, the maximum turbulent dissipation, and the buoyancy diffusivity linearly increased with the increase in the incident wave energy.