高分三号卫星对海浪的首次定量遥感

高分三号(GF-3)是我国首颗C频段多极化高分辨率微波遥感卫星,于2016年8月10日在太原卫星发射中心成功发射。GF-3 SAR卫星入射角范围约为20°—50°,具备单极化、双极化和全极化等多极化工作能力,还是世界上成像模式最多的SAR卫星,具有12种成像模式。不仅涵盖了传统的条带、扫描成像模式,而且可在聚束、条带、扫描、波浪、全球观测、高低入射角等多种成像模式下实现自由切换,既可以探地,又可以观海,达到"一星多用"的效果。近日,国家海洋局第二海洋研究所卫星海洋环境动力学国家重点实验室利用首批GF-3合成孔径雷达(SAR)遥感数据(图1)对夏威夷西北部附近太平洋海域的海浪进行了首次定量分析和反演研究(图2)。图1为GF-3 SAR的灰度图像,成像时间为2016年9月2日8:30(GMT),卫星此时飞行速度约为7.6km,极化方式为VV极化,飞行方向为降轨,空间分辨率为8m×8m,中心入射角为28.32°。由图1可以看出,SAR图像上存在明显的黑白相间的海浪条纹,说明海浪在图像上能够顺利成像。通过提取灰度图像上的调制信息,并作傅里叶变换分析,可得到包含海浪信息的图像谱。进一步,基于经典的Hasselmann SAR海浪成像模型的准线性形式,同时估计倾斜调制、水动力调制和聚束调制等三种海浪调制函数(MTF),可以利用图像谱反演得到海浪谱,此时的海浪谱主要为较长波长的涌浪信息,至于较短波长的海浪信息提取,由于受到方位向截断效应的影响,则需要引入初猜谱加以补偿实现。图2为图1反演的海面涌浪参数。从图2可以看出,该海域海浪由西北向东南传播(即由外海向近岸传播),平均波长约200m,有效波高从2.5m到4.0m不等,能够反映浪场的分布差异。由于没有同步的现场观测资料和其他卫星遥感资料,本文将这些结果与欧洲中期天气预报中心(ECMWF)提供的ERA-Interim再分析数据进行了比对。初步反演与比对结果表明,两者有较好的一致性,但本文的反演结果反映了更多的细节,显示GF-3 SAR有能力对海面涌浪信息进行高分辨率的观测;同时,再次表明ERA-Interim再分析数据低估了有效波高,因此GF-3卫星的发射将有利于提高全球海浪的遥感观测水平。     

基于SVM回归模型的SAR图像有效波高估测

A new method for estimating significant wave height(SWH) from advanced synthetic aperture radar(ASAR) wave mode data based on a support vector machine(SVM) regression model is presented. The model is established based on a nonlinear relationship between σ0, the variance of the normalized SAR image, SAR image spectrum spectral decomposition parameters and ocean wave SWH. The feature parameters of the SAR images are the input parameters of the SVM regression model, and the SWH provided by the European Centre for Medium-range Weather Forecasts(ECMWF) is the output parameter. On the basis of ASAR matching data set, a particle swarm optimization(PSO) algorithm is used to optimize the input kernel parameters of the SVM regression model and to establish the SVM model. The SWH estimation results yielded by this model are compared with the ECMWF reanalysis data and the buoy data. The RMSE values of the SWH are 0.34 and 0.48 m, and the correlation coefficient is 0.94 and 0.81, respectively. The results show that the SVM regression model is an effective method for estimating the SWH from the SAR data. The advantage of this model is that SAR data may serve as an independent data source for retrieving the SWH, which can avoid the complicated solution process associated with wave spectra.     

一种利用浮标站资料改进海浪模式有效波高预报的方法

基于浮标站海浪历史数据,利用回归分析方法建立了海浪数值模式有效波高预报产品的一元二次回归方程订正统计模型。通过2017年7月1日-2018年10月10日期间业务试运行结果发现:订正方程能有效改善有效波高数值预报产品的预报精度,且预报时效越短订正效果越显著。其中,第6~11 h预报时效内的订正前后平均绝对误差值减小0.17~0. 241 m,第6~18 h预报时效内订正前后均方根误差减小幅度为0.103~0. 28 m。这说明应用订正统计模型对海浪模式输出产品进行订正,也是改进海浪模式预报准确率的一种有效途径。     

利用C波段RADARSAT-2单极化SAR图像的方位向截断波长估计有效波高

This paper proposes two simple models, look-up table(LUT) model and empirical model, to directly retrieve significant wave height(Hs) using synthetic aperture radar(SAR) azimuth cutoff(λc). Both models aim at C-band VV, HH, VH, and HV single-polarization SAR images. The LUT model relates Hs to λc, while the empirical model relates Hs to both λc and SAR range-to-velocity(β). The LUT model coefficients are derived by simulation under different sea states and observation conditions, which depend on incidence angle(θ), wave direction(dw), and βbut are independent of polarization. The empirical model coefficients are obtained by fitting the collocated data,which only depend on polarization. To fit empirical model coefficients and validate the two models, C-band RADARSAT-2 fine quad-polarization(VV+HH+VH+HV) single-look complex(SLC) SAR images and collocated buoy data are collected. Retrieved Hs, using Yang model and the two models proposed in this paper from four kinds of polarization SAR data, are compared with buoy Hs. Results show that both LUT and empirical models have the capacity of retrieving Hs from C-band RADARSAT-2 co-polarization SAR data, while Yang model is not suitable for these kinds of SAR data. Moreover, the empirical model is also valid for cross-polarization SAR data showing clear ocean wave stripes.     

利用X波段雷达图像估计有效波高

海浪有效波高与雷达的海杂波强度有关,但是无法直接由雷达图像得到.借鉴运用SAR图像计算有效波高的方法,即假设有效波高与雷达回波强度信噪比的平方根成线性关系,可以由X波段雷达图像计算得到海浪的有效波高.将在小麦岛和南海分别进行的岸基试验和船基试验获得的浮标资料和雷达资料结合起来分析,结果表明用X波段航海雷达测量有效波高的最大误差不超过9%.     

基于准线性近似方法的SAR图像海浪参数反演研究

本文根据相干斑噪声的时间快变特征和非海浪纹理现象的时间缓变特征，基于交叉谱提出了一种对相干斑噪声和大尺度非海浪纹理的抑制的方法，进而结合SAR图像谱和海浪谱之间的准线性映射关系，基于SAR数据对海浪参数进行了反演。在反演过程中，首先仿真分析了不同海况下准线性近似法的海浪反演能力，结果表明：风浪引起的方位向截断效应会显著影响反演精度，因此该方法在低风速时的涌浪反演精度更高。通过将基于Sentinel-1卫星2020年的波模式SAR数据的反演结果与欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts, ECMWF)提供的再分析数据进行对比，发现高海况海浪有效波高反演结果明显偏低，而且该反演误差与风速、方位向截断波长之间存在显著相关性。为了提高有效波高的反演精度，本文进一步给出了海浪有效波高反演误差与风速、方位向截断波长之间的经验校正函数模型，结果显示，通过该模型修正后的海浪有效波高反演结果与ECMWF数据和浮标测量数据具有良好一致性。     

“海洋二号”有效波高数据在多源卫星高度

利用反距离加权法根据不同的方案融合得到2011年10月1日至2012年3月8日“HY-2”和“Jason-1/2”, “ENVISAT”卫星高度计的10 d平均波高, 借助NDBC提供的浮标数据对融合结果进行验证和比较, 结果表明:不同的融合结果均可以很好地反应海洋表面的真实情况;“HY-2”卫星高度计数据在多源卫星高度计数据融合中可以与其他卫星高度计数据一样发挥作用, 甚至可以替代已经停止运行的“ENVISAT”参与多源卫星高度计数据融合。     

基于SVR和LSTM算法反演散射计下有效波高研究

有效波高反演对于海洋工程及海洋环境安全具有重要意义.我国海洋二号(HY-2A)卫星载有散射计和高度计等获取海洋要素的仪器.散射计可获取海洋风场数据但无法直接获取有效波高数据,高度计可获取海洋有效波高数据但覆盖区域狭小.本文将散射计与高度计各自优势结合,利用支持向量回归(SVR)和长短期记忆(LSTM)智能算法反演散射计...     

神经网络方法对海浪有效波高数值模拟的改进

结合东北太平洋浮标资料,使用神经网络模型对WAVEWATCHⅢ海浪模式模拟的有效波高进行训练模拟,并与增加风场作为输入项的神经网络模型作了对比分析。通过分析浮标观测资料、WAVEWATCHⅢ数值模式和神经网络模拟的海浪有效波高大小,可以看出使用神经网络结合数值模式能够较好地提高有效波高的模拟精度。     

SWAN模式对陆架浅水区有效波高的模拟研究及改进

验证了QSCAT/NCEP混合风场,并将其作为SWAN模式的驱动风场。以南黄海海域作为目标区域,对SWAN模式在陆架浅水区有效波高的模拟能力进行了研究。研究表明,默认参数下SWAN模式计算的有效波高较JASON-1卫星高度计数据偏小,最大偏差达0.6 m。通过对SWAN模式中各物理过程的分析,确定模式计算值偏小的原因是白浪耗散过大。采用参数修正法对白浪耗散项进行改进,将SWAN模式计算有效波高的均方根误差降低到0.16 m以下,相关系数提高到0.85以上。选择2002年中具有代表性的4个月对改进后SWAN模式进行验证,结果显示SWAN模式在研究区域具有良好的稳定性和适用性。     

使用参数化初猜测谱算法提取VV极化Sentinel-1 SAR图像的海浪要素并与ECMWF再分析数据以及浮标数据进行验证对比

The purpose is to study the accuracy of ocean wave parameters retrieved from C-band VV-polarization Sentinel-1Synthetic Aperture Radar(SAR) images, including both significant wave height(SWH) and mean wave period(MWP), which are both calculated from a SAR-derived wave spectrum. The wind direction from in situ buoys is used and then the wind speed is retrieved by using a new C-band geophysical model function(GMF) model,denoted as C-SARMOD. Continuously, an algorithm parameterized first-guess spectra method(PFSM) is employed to retrieve the SWH and the MWP by using the SAR-derived wind speed. Forty–five VV-polarization Sentinel-1 SAR images are collected, which cover the in situ buoys around US coastal waters. A total of 52 subscenes are selected from those images. The retrieval results are compared with the measurements from in situ buoys. The comparison performs good for a wind retrieval, showing a 1.6 m/s standard deviation(STD) of the wind speed, while a 0.54 m STD of the SWH and a 2.14 s STD of the MWP are exhibited with an acceptable error.Additional 50 images taken in China's seas were also implemented by using the algorithm PFSM, showing a 0.67 m STD of the SWH and a 2.21 s STD of the MWP compared with European Centre for Medium-range Weather Forecasts(ECMWF) reanalysis grids wave data. The results indicate that the algorithm PFSM works for the wave retrieval from VV-polarization Sentinel-1 SAR image through SAR-derived wind speed by using the new GMF C-SARMOD.     

台风中降雨对中国高分3号合成孔径雷达观测信号的影响

Gaofen-3(GF-3), a Chinese civil synthetic aperture radar(SAR) at C-band, has operated since August 2016.Remarkably, several typhoons have been captured by GF-3 around the China Seas over its last two-year mission.In this study, six images acquired in Global Observation(GLO) and Wide ScanSAR(WSC) modes at verticalvertical(VV) polarization channel are discussed. This work focuses on investigating the observation of rainfall using GF-3 SAR. These images were collocated with winds from the European Centre for Medium-Range Weather Forecasts(ECMWF), significant wave height simulated from the WAVEWATCH-III(WW3) model, sea surface currents from climate forecast system version 2(CFSv2) of the National Centers for Environmental Prediction(NCEP) and rain rate data from the Tropical Rainfall Measuring Mission(TRMM) satellite. Sea surface roughness,was compared with the normalized radar cross section(NRCS) from SAR observations, and indicated a 0.8 correlation(COR). We analyzed the dependences of the difference between model-simulated NRCS and SARmeasured NRCS on the TRMM rain rate and WW3-simulated significant wave height. It was found that the effects of rain on SAR damps the radar signal at incidence angles ranging from 15° to 30°, while it enhances the radar signal at incidence angles ranging from 30° to 45° and incidence angles smaller than 10°. This behavior is consistent with previous studies and an algorithm for rain rate retrieval is anticipated for GF-3 SAR.     

中国HY-2卫星雷达高度计有效波高真实性检验

Chinese Haiyang-2(HY-2) satellite is the first Chinese marine dynamic environment satellite. The dual-frequency(Ku and C band) radar altimeter onboard HY-2 has been working effective to provide operational significant wave height(SWH) for more than three years(October 1, 2011 to present).We validated along-track Ku-band SWH data of HY-2 satellite against National Data Buoy Center(NDBC) in-situ measurements over a time period of three years from October 1, 2011 to September 30, 2014, the root mean square error(RMSE) and mean bias of HY-2SWH is 0.38 m and(–0.13±0.35) m, respectively. We also did cross validation against Jason-2 altimeter SWH data,the RMSE and the mean bias is 0.36 m and(–0.22±0.28) m, respectively. In order to compare the statistical results between HY-2 and Jason-2 satellite SWH data, we validated the Jason-2 satellite radar altimeter along-track Ku-band SWH data against NDBC measurements using the same method. The results demonstrate the validation method in this study is scientific and the RMSE and mean bias of Jason-2 SWH data is 0.26 m and(0.00±0.26) m,respectively. We also validated both HY-2 and Jason-2 SWH data every month, the mean bias of Jason-2 SWH data almost equaled to zero all the time, while the mean bias of HY-2 SWH data was no less than –0.31 m before April2013 and dropped to zero after that time. These results indicate that the statistical results for HY-2 altimeter SWH are reliable and HY-2 altimeter along-track SWH data were steady and of high quality in the last three years. The results also indicate that HY-2 SWH data have greatly been improved and have the same accuracy with Jason-2SWH data after April, 2013. SWH data provided by HY-2 satellite radar altimeter are useful and acceptable for ocean operational applications.     

一种基于Envisat ASAR波模式数据的有效波高反演算法

The main objective of this paper is to propose a newly developed ocean Significant Wave Height(SWH) retrieval method from Envisat Advanced Synthetic Aperture Radar(ASAR) imagery. A series of wave mode imagery from January, April and May of 2011 are collocated with ERA-Interim reanalysis SWH data. Based on the matched datasets, a simplified empirical relationship between 22 types of SAR imagery parameters and SWH products is developed with the Genetic Algorithms Partial Least-Squares(GA-PLS) model. Two major features of the backscattering coefficient σ_0 and the frequency parameter S_(10) are chosen as the optimal training feature subset of SWH retrieval by using cross validation. In addition, we also present a comparison of the retrieval results of the simplified empirical relationship with the collocated ERA-Interim data. The results show that the assessment index of the correlation coefficient, the bias, the root-mean-square error of cross validation(RMSECV) and the scattering index(SI) are 0.78, 0.07 m, 0.76 m and 0.5, respectively. In addition, the comparison of the retrieved SWH data between our simplifying model and the Jason-2 radar altimeter data is proposed in our study.Moreover, we also make a comparison of the retrieval of SWH data between our developed model and the wellknown CWAVE_ENV model. The results show that satisfying retrieval results are acquired in the low-moderate sea state, but major bias appears in the high sea state, especially for SWH5 m.     

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.     

预测有效波高的深度学习模型研究

研究基于RNN、LSTM、GRU深度学习模型,针对NOAA浮标数据集中的44013、44014、44017浮标的数据,通过斯皮尔曼相关性分析提高模型预测效果。实验结果表明,在进行相关性分析后,S-RNN、S-LSTM、 S-GRU的预测效果均比原始RNN、LSTM、GRU模型预测效果好。此外,提出一种基于LSTM的LSTM-Attention 波高预测模型,并进行相关实验,量化LSTM-Attention模型的预测效果,实验结果表明LSTM-Attention模型有更好的预测效果。为评估模型的泛化能力,研究还提出了一种采用邻近浮标数据进行学习,预测浮标缺失数据的方 法。实验结果表明,该方法的预测精度可以达到97.93%。本研究为海浪预测提供了新的方法和思路,也为未来深 度学习模型在海浪预测中的应用提供了参考。     

Shenzhou-4 spaceborne altimeter waveform processing and significant wave height retrieval

1IntroductionSZ-4spaceborne,launched on30December2002,is the fourth unmanned spaceborne of China.Multimode microwave remote sensor is its primaryload,of which the altimeter mode is the most impor-tant one.During about five months,a large numberofSZ-4altim…