GPS反射信号的海洋应用

本文描述了全球定位导航系统(GPS)的在海面的反射信号在海洋领域的应用，同时建立了利用GPS前向散射信号测量海面粗糙度和海面地形(Topography)的反演算法；海面粗糙度与风速和风向有直接的关系，回波功率的前缘形状及时延与海面地形相关；利用卫星主动雷达和国外机载数据的结果分析比较表明，GPS反射计作为遥感工具有两个优点：即比传统微波主动雷达高的空间分辨率和快速的时间分辨率。     

岸基GNSS反射信号延迟多普勒图变化特征分析

GNSS-R(Global Navigation Satellite System-Reflections)技术是一种可用于海面参数遥感探测的新兴技术,是当前国内外的研究热点之一。GNSS海面反射信号的延迟多普勒图像(DelayDoppler Map,DDM)是GNSS-R技术应用的重要理论基础。基于Z-V模型确定了岸基GNSS反射信号DDM的理论仿真方法和流程,并进行了实验验证,实验结果显示,仿真DDM与实测DDM之间具有很好的一致性。同时,全面系统地分析了风速、风向、接收天线高度和卫星仰角4个参数对DDM变化特征的影响规律,仿真结果显示:风速和接收天线高度的影响呈单调特性;风向的影响呈周期特性;卫星仰角的影响则比较特殊,在30°以下时呈现明显的单调性,大于30°时误差几乎相同,即具有窗函数选择特性,这种变化特征为利用DDM反演海面波浪、风场数和海面粗糙度等参数提供了理论依据和参考。     

利用机载GNSS反射信号反演海面风速的研究

全球卫星导航定位系统的反射信号(GNSS-R)遥感技术作为一种新型的、低成本的、高机动性的海面微波遥感测风技术,与其他测风手段优势互补,可以增加测风手段的多样性,弥补局部测风手段不足的状况。研究了接收机在机载高度时,GPS反射信号功率理论模型四部分函数的性质,在此基础之上,数值模拟了机载高度下理论相关功率波形,基于海面风速对波形峰值与后沿的影响,提出了一种能够兼顾所有理论波形信息的二维插值风速反演方法。利用该方法,结合实测机载数据对海面风速进行反演,反演的风速均值与附近测站风速均值相差为1.4 m/s,与浮标数据相一致。     

GNSS海面散射信号的海面风场遥感

首先介绍了GNSS-R海洋遥感的发展过程,简要叙述了利用GNSS信号进行海洋动力学遥感测量的基本原理;然后分析了接收机天线的部分参数对GNSS信号可见数量的影响;最后详细阐述了GNSS海面散射信号的理论模型,该模型是利用几何光学限制的基尔霍夫近似建立起来的,并根据理论模型进行了数值仿真试验。结果证实了通过散射信号功率波形来反演海面风场理论上的可行性。     

国外导航卫星的反射信号在海洋遥感中的应用分析

为了详细了解国外导航卫星反射信号(GNSS-R)新技术在海洋遥感方面的研究进展,在外文文献调研的基础上,详细阐述了GNSS-R的反射率和双基雷达两种方法在海洋监测应用的遥感技术原理,介绍了海面风场遥感、海冰遥感和海面测高等方面的的理论模型和方法。高空间分辨率和高时间分辨率是导航卫星反射信号在海洋中应用的最大优势。虽然GNSS-R的海洋气象要素反演技术只有十几年的发展历史,但是它已经在海洋天气预报和全球气候变化等领域显示出广阔的应用前景和巨大潜力,未来必将成为海洋大气探测的重要技术手段。     

GNSS-R信号测量海面有效波高的应用

全球导航卫星系统反射信号GNSS-R(Global Navigation Satellite System Reflection)可以用来探测海面状态,如海面有效波高。介绍了利用GNSS反射信号测量海面有效波高的观测设备和技术原理,对2009年博贺实验的数据进行了分析,证实了反射信号相关复数场的有效相关时间与有效波高之间存在相关,利用前期观测的数据进行拟合,得出反演公式,并与测波雷达数据进行了结果对比,平均误差0.136m,相对偏差12.02%,均方根偏差0.169m。     

基于MISR多角度光学遥感影像的海面风速估算实验

太阳耀光是来自粗糙海面的直接太阳反射光,其强度与海面粗糙度密切相关,而海面粗糙度主要受海面风场影响。因此,包含太阳耀光信息的光学遥感影像在海洋动力过程和海面风速探测中具有积极意义。本文利用2016年2月到2017年3月期间成像的25幅Terra卫星MISR(Multi-angle Imaging Spectro Radiometer)传感器的多角度遥感影像,分别提取了太阳的高度角和方位角、正视和后视影像的卫星观测角、方位角等信息,校正获得正视和后视影像的太阳耀光辐射强度,进一步反演海表面粗糙度信息,进而计算海面风速。最后利用ECMWF(European Centre for Medium-Range Weather Forecasts)的模式风速数据与反演获得的风速结果进行对比验证。结果表明,两者的相关系数较高(R=0.745),均方根误差和平均绝对偏差值分别为1.514 m·s^-1和1.319 m·s^-1。初步实验结果表明,利用MISR多角度光学遥感影像估算海表面风速是可行性的。     

基于SAR图像雨团足印的海面风向提取方法

利用地球物理模式函数进行SAR海面风速反演时,需以风向作为地球物理模式函数的输入。本文应用了一种利用SAR图像上雨团足印顺风一侧比逆风一侧明亮的图像特征的海面风向提取方法,以进行海面风速反演。4景RADARSAT-2卫星SAR示例数据风向提取结果相对于ASCAT散射计的风向均方根误差满足不大于16°。分别以本文方法提取的风向和ASCAT散射计风向作为输入,利用地球物理模式函数CMOD5进行海面风速的SAR反演,两者的风速反演结果基本一致,其均方根误差差值不超过0.3 m/s。本文利用SAR图像雨团足印信息的风向提取方法准确可靠,可应用于SAR海面风速反演。     

星载合成孔径雷达图像反演海面风场方法综述

合成孔径雷达作为微波雷达,具有全天候、全天时、高分辨率对海面成像的能力,能实现多波段、多极化、多视角的观测海面,可以提供大范围、高精度的实时动态海面信息,对其进行深入研究有着重要的意义。基于合成孔径雷达图像,分别以风向和风速反演为主线,对国内外海面风场反演的工作进行了详细的总结。     

星载SAR对雨团催生海面风场的观测研究

雨团或对流雨是热带与亚热带地区的主要降雨形式,较易被高分辨率星载合成孔径雷达（SAR）探测到。SAR图像上的雨团足印是由大气中雨滴的散射与吸收、下沉气流等共同导致形成的。本文以RADARSAT-2卫星100 m分辨率的SAR图像上雨团引起的海面风场及其结构反演与解译作为实例进行分析。使用CMOD4地球物理模式函数,分别以NCEP再分析数据、欧洲MetOp-A卫星先进散射计（ASCAT）和中国HY-2卫星微波散射计的风向为外部风向,进行了SAR图像的海面风场反演。反演的海面风速相对于NCEP、ASCAT和HY-2的均方根误差（RMSE）分别为1.48 m/s,1.64 m/s和2.14 m/s。SAR图像上一侧明亮另一侧昏暗的圆形信号图斑被解译为雨团携带的下沉气流对海面风场（海面粗糙度）的改变所致。平行于海面背景风场其通过雨团圆形足印中心的剖面上的风速变化可拟合为正弦或余弦曲线,其拟合线性相关系数均不低于0.80。背景风场的风速大小、雨团引起的风速大小以及雨团足印的直径可利用拟合曲线获得,雨团足印的直径大小一般为数千米或数十千米,本文的8例个例解译与分析均验证了该结论。     

Shallow water sonar conditions and radar backscatter

Sonar and radar signals are scattered by small-scale roughness elements of the sea surface. The characteristics of the roughness depend on wind, sea state, and other parameters. The dynamic behavior of this roughness is related to radar backscatter. This provides a basis for oceanographic radar remote sensing and links radar backscatter to sea surface sonar conditions. Radar, sonar, and roughness experiments have been conducted at the research platform NORDSEE. The results are combined and provide a method to estimate shallow water sonar conditions from the radar backscatter cross section     

Bistatic radar cross section of ship targets

Data describing simultaneous bistatic (BRCS) and monostatic (MRCS) radar cross sections of freighters atXband near grazing incidence are presented. Details of the measurement system, calibration procedures, and target-radar geometry are discussed. Measured values of BRCS were typically less than the MRCS and decreased with increasing bistatic angle.     

各向异性海面全极化微波双站散射的仿真与分析

为了解各向异性随机粗糙海面的微波双站散射机制及其特性,本文利用解析近似的积分方程模型以及一种改进的半经验海浪谱模型实现了对各向异性随机粗糙海面的全极化微波散射仿真模拟,并与卫星观测数据、经验的地球物理模式函数及已有的解析近似散射模型仿真结果进行了对比,验证了仿真结果的可行性和准确性。利用该模型分析了入射波频率、入射角、极化方式、海面风速及风向等参数对各向异性海面双站散射的影响。模拟结果表明,在不同的入射角、散射角及方位角等观测几何条件下,海面不同波段的双站散射表现出不同的空间散射特性,且对风速、风向等海面动力学参数表现出不同的敏感性,以L波段为例,海面向后半球双站散射在各个极化方式下都对风速较为敏感,而在同极化方式下,其对风向的响应在中低风速和高风速条件下相反,整体而言,低风速下海面双站散射对风向更为敏感。这表明对于海面动力参数的反演,双站散射可以提供比传统单站雷达后向散射更丰富的物理信息。本文探讨了各向异性海面微波双站散射特性,为基于主动式及分布式微波传感器的海洋动力参数遥感反演提供了理论分析基础。     

Extraction of Ocean Wave Spectra From Simulated Noisy Bistatic High-Frequency Radar Data

An algorithm is developed for the inversion of bistatic high-frequency (HF) radar sea echo to give the nondirectional wave spectrum. The bistatic HF radar second-order cross section of patch scattering, consisting of a combination of four Fredholm-type integral equations, contains a nonlinear product of ocean wave directional spectrum factors. The energy inside the first-order cross section is used to normalize this integrand. The unknown ocean wave spectrum is represented by a truncated Fourier series. The integral equation is then converted to a matrix equation and a singular value decomposition (SVD) method is invoked to pseudoinvert the kernel matrix. The new algorithm is verified with simulated radar Doppler spectrum for varying water depths, wind velocities, and radar operating frequencies. To make the simulation more realistic, zero-mean Gaussian noise from external sources is also taken into account     

Vandemark-Chapron算法与Young算法联合反演Jason-1雷达高度计海面风速方法研究

针对现行业务运行的宙达高度计风速反演算法从考虑0～20 m/s的缺陷,提出了vc算法(Vandemark-Chapron Algorithm)与Young算法联合反演高度计海曲风速的方法,通过对Jason-1资料的统计试验,确定了后向散射截面临界点,当高度计后向散射截面大于等于该临界点时采用VC算法反演风速,反之采用Y...     

Wind- and wave-field measurements using marine X-band radar-image sequences

This paper describes two algorithms for the retrieval of high-resolution wind and wave fields from radar-image sequences acquired by a marine X-band radar. The wind-field retrieval algorithm consists of two parts. In the first part, wind directions are extracted from wind-induced streaks, which are approximately in line with the mean surface wind direction. The methodology is based on the retrieval of local gradients from the mean radar backscatter image and assumes the surface wind direction to be oriented normal to the local gradient. In the second part, wind speeds are derived from the mean radar cross section. Therefore, the dependence of the radar backscatter on the wind vector and imaging geometry has to be determined. Such a relationship is developed by using neural networks (NNs). For the verification of the algorithm, wind directions and speeds from nearly 3300 radar-image sequences are compared to in situ data from a colocated wind sensor. The wave retrieval algorithm is based on a methodology that, for the first time, enables the inversion of marine radar-image sequences to an elevation-map time series of the ocean surface without prior calibration of the acquisition system, and therefore, independent of external sensors. The retrieved ocean-surface elevation maps are validated by comparison of the resulting radar-derived significant wave heights, with the significant wave heights acquired from three colocated in situ sensors. It is shown that the accuracy of the radar-retrieved significant wave height is consistent with the accuracy of the in situ sensors.     

On the Development of a Second-Order Bistatic Radar Cross Section of the Ocean Surface: A High-Frequency Result for a Finite Scattering Patch

The development of a model for the second-order bistatic high-frequency (HF) radar cross section on an ocean surface patch remote from the transmitter and receiver is addressed. A new approach is taken that allows a direct comparison with existing monostatic cross sections for finite regions of the ocean surface. The derivation starts with a general expression for the bistatically received second-order electric field in which the scattering surface is assumed to be of small height and slope. The source field is taken to be that of a vertically polarized dipole, and it is assumed that the ocean surface can be described, as is usually done, by a Fourier series in which the coefficients are zero-mean Gaussian random variables. Subsequently, a bistatic cross section of the surface, normalized to patch area, is derived. The result is verified by the following two means: 1) the complete form of the bistatic HF radar cross section in backscattering case is shown to contain an earlier monostatic result that has, itself, been used extensively in radio oceanography applications; and 2) the bistatic electromagnetic coupling coefficient is shown to reduce exactly to the monostatic result when backscattering geometry is imposed. The model is also depicted and discussed based on simulated data     

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.     

非线性随机表面Kirchhoff散射模式

自从Longuet-Higgins(1963)根据非线性作用导出较正态分布为准确的波面高度分布以来,高阶矩在军事、高科技等很多方面得到了应用。Huang等(1980)曾检验了高阶非线性分布,他们发现当波高概率分布中包含直至四阶项时与观测结果相符,但如果考虑更高阶不但效果不好,反而更差。尽管非线性随机过程在海浪理论中获得了广泛的应用,但对海面雷达散射研究仅讨论到三阶矩的影响(Fung et al.,1991;Chen et al.,1992)。电磁随机表面散射理论有适应大尺度随机起伏的粗糙面的Kirchhoff散射模式、轻度粗糙表面的微扰散射模式、大小尺度独立叠加的双尺度散射模式、全波散射模式(Bahar,1987)和积分方程散射模式(Chen et al.,1992)等。对随机粗糙Kirchhoff表面电磁波散射问题,尽管几十年来许多科学家已经进行了大量研究(Fung et al.,1991; Ulaby et al.,1982;Wu et al.,1988),Eom等(1983)曾对Gaussian面和非 Gaussian面的散射特性进行了比较研究,Fung等(1991)将Kirchhoff散射模式推广应用到三阶粗糙随机表面,然而更高阶矩对散射截面的影响还未见报道。本文在Fung等(1991)的基础上将Kirchhoff散射模式推广应用到四阶项,并对模式的应用进行了分析和讨论。