GPS浮标数据反演海浪谱的理论仿真与试验验证

GPS浮标作为一种新型的海洋测量设备,近年来在海面高度现场测量和星载高度计定标方面取得了重要应用。通过仿真试验对反演海浪谱的方法和流程进行研究,旨在探索从GPS浮标测量的海面高度序列中提取海浪谱的方法。首先,使用Longuest模型生成了海浪波面位移时间序列,并通过Pierson-Moscowitz风浪谱对波面位移的统计特性进行约束,其随机性由相位引入。结合典型潮汐和GPS浮标仪器噪声的仿真时间序列,合成了仿真时间长度1h的1Hz(每秒1次)随机海面高度序列。然后,利用自相关函数法,进行高通滤波和数据压缩,得到了仿真的海浪谱。该仿真结果和理论海浪谱非常接近,可满足海浪谱反演的需求。最后,通过山东石岛外海的GPS浮标现场试验,验证了本文提出的反演方法的适用性。本文的研究解决了GPS浮标反演海浪谱的关键问题,丰富了海浪谱反演的手段,拓展了GPS浮标的应用领域,有望为未来我国的星载波谱仪定标服务。     

A semiempirical algorithm for SAR wave height retrieval and its validation using Envisat ASAR wave mode data

Many synthetic aperture radar(SAR) wave height retrieval algorithms have been developed.However,the wave height retrievals from most existing methods either depend on other input as the first guess or are restricted to the long wave regime.A semiempirical algorithm is presented,which has the objective to estimate the wave height from SAR imagery without any prior knowledge.The proposed novel algorithm was developed based on the theoretical SAR ocean wave imaging mechanism and the empirical relation between two types of wave period.The dependency of the proposed model on radar incident and wave direction was analyzed.For Envisat advanced synthetic aperture radar(ASAR) wave mode data,the model can be reduced to the simple form with two input parameters,i.e.,the cutoff wavelength and peak wavelength of ocean wave,which can be retrieved from SAR imagery without any prior knowledge of wind or wave.Using Envisat ASAR wave mode data and the collocated buoy measurements from NDBC,the semiempirical algorithm is validated and compared with the Envisat ASAR level 2 products.The root-mean-square-error(RMSE) and scatter index(SI) in respect to the in situ measurements are 0.52 m and 19% respectively.Validation results indicate that,for Envisat ASAR wave mode data,the proposed method works well.     

降雨条件下的导航X波段雷达海浪参数反演算法研究

X波段的电磁波受降雨影响容易产生衰减,这导致导航X波段雷达在降雨时无法用于海浪观测。本文提出了一种新的降低降雨影响的算法来反演海浪参数。首先,对X波段雷达图像做主成分分析,获得波浪变化的主成分,利用一维傅里叶变换得到波数谱,对其滤波减小降雨对雷达图像的影响;然后,选取JONSWAP(Joint North Sea Wave Project)谱作为理论谱,建立以观测谱与理论谱的最小化差异为目标函数的模型,求解该模型估算海浪的有效波高。与浮标测量的有效波高相比,该方法反演的有效波高的均方根误差是0.23 m,证明了该方法的可行性。     

基于相参X-波段海洋雷达的海表轮廓测量研究

X-波段海洋雷达测量所得海面散射单元的多普勒信息与散射单元的雷达视向速度密切相关。首先,基于符号多普勒估计方法,对X-波段雷达海面回波的多普勒频移信息进行了估计;在此基础上,应用各分辨单元回波的多普勒频移信息,建立了海浪表面轮廓的反演算法。该算法中,同时考虑了雷达入射角、方位角和雷达空间分辨率等因素对反演结果的影响。通过将反演结果与浮标测量数据相比较,发现雷达空间分辨率起到了类似低通滤波的作用,该作用对短重力波谱影响显著。同时,还应用加拿大麦克马斯特大学的IPIX雷达数据对海表轮廓进行了反演,并将反演所得有效波高、海浪周期与现场测量数据进行了比较,反演结果与现场测量结果吻合较好。     

Satellite-based RAR performance simulation for measuring directional ocean wave spectrum based on SAR inversion spectrum

Some missions have been carried out to measure wave directional spectrum by synthetic aperture radar (SAR) and airborne real aperture radar (RAR) at a low incidence. Both them have their own advantages and limitations. Scientists hope that SAR and satellite-based RAR can complement each other for the research on wave properties in the future. For this study, the authors aim to simulate the satellite-based RAR system to validate performance for measuring the directional wave spectrum. The principal measurements are introduced and the simulation methods based on the one developed by Hauser are adopted and slightly modified. To enhance the authenticity of input spectrum and the wave spectrum measuring consistency for SAR and satellite-based RAR, the wave height spectrum inversed from Envisat ASAR data by cross spectrum technology is used as the input spectrum of the simulation system. In the process of simulation, the sea surface, backscattering signal, modulation spectrum and the estimated wave height spectrum are simulated in each look direction. Directional wave spectrum are measured based on the simulated observations from 0° to 360°. From the estimated wave spectrum, it has an 180° ambiguity like SAR, but it has no special high wave number cut off in all the direction. Finally, the estimated spectrum is compared with the input one in terms of the dominant wave wavelength, direction and SWH and the results are promising. The simulation shows that satellite-based RAR should be capable of measuring the directional wave properties. Moreover, it indicates satellite-based RAR basically can measure waves that SAR can measure.     

Estimates of ocean wavelength and direction from X-and L-band synthetic aperture radar data collected during the Marineland experiment

Simultaneously obtainedX- andL-band synthetic aperture radar (SAR) data collected during the Marineland Experiment were spectrally analyzed by fast Fourier transform (FFT) techniques to estimate ocean wavelength and direction. An eight-sided flight pattern was flown over the same ocean area in order to study the sensitivity of the spectral estimate on radar look direction. These spectral estimates were compared with in situ wave measurements made by a pitch-and-roll buoy. The comparison revealed that theX-band SAR detected all gravity waves independent of radar look direction, while theL-band SAR detected all range-traveling gravity waves but failed to detect waves in three of four cases in which the waves were traveling within 25° of the azimuth direction. The analysis also indicates that azimuth-traveling waves appear longer and more range-traveling in the SAR imagery than observed by in situ instrumentation. It is postulated that degraded azimuth resolution due to scatterer motion is responsible for these observations.     

Estimation of wave directional spreading in shallow water

This paper describes wave directional spreading in shallow water. Waves were measured for a period of 2 months using the Datawell directional waverider buoy at 15 m water depth on the east coast of India in the Bay of Bengal. The study also showed that in shallow water wave directional spreading was narrowest at peak frequency and widened towards lower and higher frequencies. The wind direction was found to deviate from the wave direction during most of the time. The unidirectional spectrum was found to be satisfactorily represented by Scott spectra.     

机载波谱仪海浪谱反演方法及其验证

为了验证波谱仪反演二维海浪谱的功能,根据海浪波谱仪的信号形成机制,总结了机载波谱仪反演海浪的流程。利用机载波谱仪回波数据,通过自相关和互相关两种功率谱估计方法,反演了二维海浪谱。最后通过与浮标测量的二维海浪谱进行对比,验证了该机载波谱仪探测二维海浪谱的有效性。结果表明,无论采用自相关函数还是互相关函数进行功率谱估计,得到的主波波长和有效波高与实际二维海浪谱基本一致。互相关函数法得到的交叉谱能去除180°模糊现象,其在计算有效波高时相对于自相关函数会稍微偏小。在计算斜率方差时可以采用5°~12°入射角范围的后向散射系数进行公式拟合,因此定标与否并不影响最后的二维海浪谱结果,未来星载波谱仪只有靠多波束联合才能实现。     

Influence of the type of sea waves on the backscattered radar cross section at medium incidence angles

We consider the influence of the sea surface state on the backscattered radar cross section and the accuracy of the wind speed retrieval from the scatterometer data. We used a joint set of radars and buoys to determine the type of sea waves. Three types of sea waves were distinguished: developing wind waves, fully developed wind waves, and mixed sea. It is shown that the retrieval error of the near surface wind speed using a one-parameter algorithm is minimal in the case of fully developed wind waves. We compared these data with the results of radio-altimeter data analysis and showed that in both cases underestimation of the retrieval wind speed exists for developing wind waves and overestimation occurs for mixed sea. A variety of swell parameters (length of the dominating wave, swell height, swell age) significantly influence the backscattered radar cross section, leading to a growth in the mean square error of the retrieved wind speed during vertical sounding (radio-altimeter data), and only slightly influence the mean square error of the scatterometer data (medium incidence angles). It is necessary to include the information about the parameters of sea waves in the algorithms and take into account the regional wave properties to increase the accuracy of wind speed retrieval.     

对ENVISAT ASAR level2算法固有误差的分析

欧洲空间局的ENVISAT ASAR level 2算法是从合成孔径雷达(SAR)单视复图像反演涌浪方向谱的算法.该算法假设双峰海浪谱的SAR图像交叉谱是涌浪的图像交叉谱和风浪的图像交叉谱之和.实际上双峰海浪谱的SAR图像交叉谱中还有一个混合项,正是该混合项导致ENVI-SAT ASAR level 2算法有固有误差.利用遥感仿真的方法分析了不同海况条件下该算法的这一固有误差,结果表明,只有在有效波高较小、或风浪的成分较少、或双峰海浪的传播方向较靠近SAR距离向、或波长较长时固有误差才较小,ENVISAT ASAR level 2算法对海浪谱的反演才较为适用.     

Measurements of directional spectrum of wind waves using an array of wave detectors Part II. Field observation

The directional spectrum of wind waves was observed using seven wave gauges at the marine observation tower. Directional spectrum was calculated by the method described inFujinawa (1974 a) which assures the high directional resolving power. Under the nearly uniform condition of the wind the process of wave development was studied by the use of the data of the directional spectrum. The analysis revealed that, 1) the mean wave direction does not necessarily coincide with the wind direction, 2) the directional spreading is much narrower than hitherto reported, 3) the growth rate is closely proportional to the square of cosine of the angle between the wave direction and the wind direction, 4) the resonancetype wave-wave interaction plays only minor role in the process of wave development for any wave component in the earlier stage of development and for the most rapidly developing component in the main stage of development.     

A Spectral Approach for Determining Altimeter Wind Speed Model Functions

We propose a new analytical algorithm for the estimation of wind speeds from altimeter data using the mean square slope of the ocean surface, which is obtained by integration of a widely accepted wind-wave spectrum including the gravity-capillary wave range. It indicates that the normalized radar cross section depends not only on the wind speed but also on the wave age. The wave state effect on the altimeter radar return becomes remarkable with increasing wind speed and cannot be neglected at high wind speeds. A relationship between wave age and nondimensional wave height based on buoy observational data is applied to compute the wave age using the significant wave height of ocean waves, which could be simultaneously obtained from altimeter data. Comparison with actual data shows that this new algorithm produces more reliable wind speeds than do empirical algorithms. This revised version was published online in July 2006 with corrections to the Cover Date.     

基于海面多普勒谱模型的海浪反演精度分析

Microwave remote sensing is one of the most useful methods for observing the ocean parameters. The Doppler frequency or interferometric phase of the radar echoes can be used for an ocean surface current speed retrieval,which is widely used in spaceborne and airborne radars. While the effect of the ocean currents and waves is interactional. It is impossible to retrieve the ocean surface current speed from Doppler frequency shift directly. In order to study the relationship between the ocean surface current speed and the Doppler frequency shift, a numerical ocean surface Doppler spectrum model is established and validated with a reference. The input parameters of ocean Doppler spectrum include an ocean wave elevation model, a directional distribution function, and wind speed and direction. The suitable ocean wave elevation spectrum and the directional distribution function are selected by comparing the ocean Doppler spectrum in C band with an empirical geophysical model function(CDOP). What is more, the error sensitivities of ocean surface current speed to the wind speed and direction are analyzed. All these simulations are in Ku band. The simulation results show that the ocean surface current speed error is sensitive to the wind speed and direction errors. With VV polarization, the ocean surface current speed error is about 0.15 m/s when the wind speed error is 2 m/s, and the ocean surface current speed error is smaller than 0.3 m/s when the wind direction error is within 20° in the cross wind direction.     

海浪成长状态对TOPEX/Poseidon高度计风速反演的影响

以墨西哥湾同步高度计、浮标资料为例,研究了海浪成长状态对高度计风速反演的影响。同步的高度计风速和浮标风速比较显示,在墨西哥湾地区,海浪成长状态对高度计风速反演有较大影响。在考虑海浪成长状态影响的条件下,利用谱模型反演高度计风速,取得了较好的效果。与目前TOPEX/Poseidon高度计风速反演业务化算法相比,在海浪未充分成长条件下,考虑海浪成长状态影响后,根据谱模型反演获得的风速与浮标风速之间的均方根误差减小了30%,平均误差减小了83%。在利用谱模型算法反演高度计风速时,谱模型中的波龄因子(表示海浪成长状态)可以根据高度计测得的有效波高和风速获得,因此该方法具有广泛的适用性。     

Ocean wave parameters retrieved directly from compact polarimetric SAR data

We aim to directly invert wave parameters by using the data of a compact polarimetric synthetic aperture radar(CP SAR) and validate the effectiveness of ocean wave parameter retrieval from the circular transmit/linear receive mode and π/4 compact polarimetric mode. Relevant data from the RADARSAT-2 fully polarimetric SAR on the C-band were used to obtain the compact polarimetric SAR images, and a polarimetric SAR wave retrieval algorithm was used to verify the sea surface wave measurements. Usin...     

极低频海洋噪声与风的相关性研究

风通过影响海洋表面从而产生200 Hz以上的深海环境噪声,但有研究指出,通过风生表面波之间的非线性相互作用产生的驻波,能够与海床共振构成海底微震,从而产生10 Hz以下的噪声。针对这一新型风生噪声机制,本研究对威克岛海域10 Hz以下的极低频噪声进行了分析。比较了不同频率下海洋环境噪声功率谱级与风速的相关性,并讨论了风速和风向对设立在威克岛南北部二组水听器三联体信号的影响,结果表明2 Hz处的海洋环境噪声级与风速相关性最好,而风速和风向变化越剧烈海洋环境极低频噪声与风速风向的相关性越好。     

CODAR wave measurements from a North Sea semisubmersible

CODAR, a high-frequency (HF) compact radar system, was operated continuously over several weeks aboard the semisubmersible oil platform Treasure Saga for the purpose of wave-height directional measurement and comparison. During North Sea winter storm conditions, the system operated at two different frequencies, depending on the sea state. Wave data are extracted from the second-order backscatter Doppler spectrum produced by nonlinearities in the hydrodynamic wave/wave and electromagnetic wave/scatter interactions. Because the floating oil rig itself moves in response to long waves, a technique has been developed and successfully demonstrated to eliminate to second order the resulting phase-modulation contamination of the echo, using separate accelerometer measurement of the platform's lateral motions. CODAR wave height, mean direction, and period are compared with data from a Norwegian directional wave buoy; in storm seas with wave heights that exceeded 9 m, the two height measurements agreed to within 20 cm RMS, and the mean direction to better than 15° RMS     

分析海浪方向谱的扩展本征矢方法Ⅰ.方法的导出

基于交叉谱矩阵可以按本征值划分为信号和含有噪音部分的思想，提出了一种自直接测量数据估计海浪方向谱的方法。该方法称为扩展本征矢方法（EEV），可应用于单点测量系统、仪器阵列以及由二者构成的复合阵列。现有的某些估计方法（如最大似然方法及其扩展形式等）仅是此方法的某种特例。     

利用高频雷达海面回波中奇异峰估计浪高的方法

The popular methods to estimate wave height with high-frequency(HF) radar depend on the integration over the second-order spectral region and thus may come under from even not strong external interference. To improve the accuracy and increase the valid detection range of the wave height measurement, particularly by the smallaperture radar, it is turned to singular peaks which often exceed the power of other frequency components. The power of three kinds of singular peaks, i.e., those around ±1,±2~(1/2) and ±1(2~(1/2)) times the Bragg frequency, are retrieved from a one-month-long radar data set collected by an ocean state monitoring and analyzing radar,model S(OSMAR-S), and in situ buoy records are used to make some comparisons. The power response to a wave height is found to be described with a new model quite well, by which obvious improvement on the wave height estimation is achieved. With the buoy measurements as reference, a correlation coefficient is increased to 0.90 and a root mean square error(RMSE) is decreased to 0.35 m at the range of 7.5 km compared with the results by the second-order method. The further analysis of the fitting performance across range suggests that the peak has the best fit and maintains a good performance as far as 40 km. The correlation coefficient is 0.78 and the RMSE is 0.62 m at 40 km. These results show the effectiveness of the new empirical method, which opens a new way for the wave height estimation with the HF radar.