卫星高度计反演海面风速——模式函数与应用实例

对迄今为止国际上公开发表的１４ 种卫星高度计反演海面风速的模式函数进行了系统的分析和比较。在此基础上提出了一种“折中”模式函数,用于ＴＯＰＥＸ 资料海面风速反演效果良好。文中对高度计获取的全球海面风速分布特征进行了分析,最后探讨了改善海面风速模式函数反演精度的途径。     

利用GPS反射信号遥感全球海态

最近几年，迅速发展一种新的空基海洋遥感技术。它在低轨卫星上接收经海面或地面反射的GPS信号，通过分析C／A或P码的相关函数波形及其后沿，荻得全球海面平均高度、浪高、海面风场以及土壤湿度等，提供海洋中尺度涡旋发展、演变信息，监测和分析高纬地区冰层厚度。它采用无源探测方法，具有成本低、荻取数据量多和快速等突出特点，是传统卫星高度计探测的重要补充和发展。该文讨论了这项技术发展现状、反射信号的时间延迟以及海态参数计算等关键技术；也给出作者利用掩星观测资料分析得到的海面反射回波一些重要结果。     

卫星高度计融合产品的研发综述EI北大核心CSCD

卫星高度计融合产品广泛应用于海洋环境监测、海洋中尺度系统研究和海洋数值预报等业务领域,具有极高的科学和社会价值。当前同时在轨运行并发布沿轨数据的高度计卫星维持在5—6颗,为了立足现有数量的观测卫星,最大程度地提取有效观测信息,本文以改进融合产品的有效分辨率并提升高度计产品质量为主线,详细介绍了现有卫星高度计融合产品的研发现状。前人的研究结果表明:卫星高度计融合产品的质量主要3个方面的影响:(1)制作产品所选择的卫星种类:因不同种类卫星的轨道高度、周期差异,导致获取的海表高度时空精度也存在不同;(2)高度计标准的选择:资料处理中的仪器参数、地球物理参数、环境校正、平均海面的构造等;(3)融合处理资料的方法:不同融合方法在背景场的选择、背景误差相关系数尺度的差异以及观测误差的处理等。针对众多学者关注的热点和亟需解决的难点问题,本文系统地介绍了当前卫星高度计融合产品的种类,对比分析产品质量。同时针对卫星高度计的融合方法,着重介绍了卫星高度计融合产品质量的提升途径。文中进一步强调了融合方法与有效分辨率的关系,为提升卫星高度计融合产品质量,提高融合产品的有效分辨率,扩大卫星高度计融合产品的应用范围提供参考。     

“凤凰”台风期间海面风、浪的高频地波雷达OSMAR071遥感

介绍了OSMAR071高频地波雷达系统“凤凰”台风期间海面风、浪遥测结果。文中给出高频地波雷达OSMAR071 风、浪反演算法：利用SBM模型结合多波束采样法反演风向;利用改进的Barrick模型进行有效波高的反演;利用SMB关系式进行风速反演。将雷达定点观测结果与浮标数据进行对比,表现出很好的一致性,其中有效波高相关系数为0.72,均方根误差为0.48 m,风向相关系数为0.97,均方根误差为27.7°,风速的相关系数为0.6,均方根误差为 3.5 m/s。文中还对影响探测精度的因素进行了分析。     

GNSS浮标海面高测量中海浪的影响分析

以骑浪式浮标为对象,研究了海浪对于GNSS浮标测量海面高的影响。以单向规则波的形态简化模拟了海浪以及GNSS浮标随海浪的运动过程。在某固定海况条件下,比较了浮标天线参考点高不同时海面高测量误差;在浮标天线参考点高固定的条件下,比较了不同海况时海面高测量误差;给出了GNSS浮标在卫星高度计定标应用中的建议。结果表明:浮标测量海面高误差与天线参考点高成正比,海浪周期越长、波高越小时,测量误差越小;天线参考点高为5 cm,海浪振幅为1 m、周期为7 s、传播速度为5 m/s时,测量海面高平均比实际海面高低0.4 mm;相同浮标在小于3~4级的海况条件下,由海浪引起的海面高测量偏差约在1 mm以内,对于卫星高度计定标工作而言,该误差影响可以忽略。     

遥感和GIS技术在全球海面风速分析中的应用

将遥感技术和地理信息系统（GIS）技术相结合，建立了一个基于TOPEX卫星的全球海面风速分析的海洋地理信息系统（MGIS）。阐述了全球海面风速分析模型的建立方法，定义了描述风场特征的相关参数，介绍了全球海面风速海洋地理信息系统的结构、功能和工作流程，另外，对系统的3个应用实例进行了地学分析。     

海洋卫星测高及其反演全球海洋重力场和海底地形模型研究进展

海洋卫星测高在全球和区域大地水准面建模、全球海洋重力场反演、海底地形探测、海平面变化监测、构造板块运动研究等大地测量领域至关重要。本文概述了海洋微波测高卫星的简要发展历程,重点梳理了卫星测高在全球海洋重力场和全球海底地形建模方面取得的丰硕成果,对比分析了主流的海洋重力场和海底地形模型;介绍了合成孔径雷达高度计、Ka频段雷达高度计、合成孔径雷达干涉仪3种先进微波测高技术,并分析了其各自的优缺点,表明它们将在未来若干年呈并驱发展趋势;较为系统地阐述了海洋卫星测高的另一新型技术,即GNSS反射信号测量技术的研究动态,给出了GNSS-R(GNSS reflectometry)类(试验)卫星的发展脉络和发展前景。卫星测高的发展趋势之一是多颗测高卫星的组网观测,本文概括了曾经提出的和拟议中的若干组网测高计划,扼要介绍了由我国提出并正在实施的双星跟飞测高模式;最后指出了卫星测高发展的几个主要关注点,包括双星跟飞测高和SWOT(surface water ocean topography)任务的2维海面高(差)测量、卫星测高反演海底地形与高级地形激光高度计观测数据及遥感卫星图像的结合、星载GNSS-R厘米级海面高的载波相位测量、人工智能技术在卫星测高中的潜在应用等。     

NOAA AVHRR数据的高精度导航定位

在ＮＯＡＡＡＶＨＲＲ数据的导航定位中,一般只进行卫星轨道根数调整和卫星时间漂移校正,而未考虑卫星姿态系统偏差的影响。该文在定量分析影响导航定位精度的主要因素基础上,提出了一种新的导航算法。该导航算法不仅对卫星时间漂移、卫星轨道根数进行约束最优调整,而且通过估算卫星姿态偏差将其融入到导航定位中加以校正。同时,将导航算法和卫星轨道计算、地面扫描点地理经纬度的计算方法以及参考图像的高精度配准等结合起来考虑,利用具有足够定位精度的１∶１００ 万数字地图数据和高分辨率Ｌａｎｄｓａｔ ＴＭ图像,提取地面控制点,使ＡＶＨＲＲ数据导航定位精度在整条轨道大范围内提高到１ 个像元,在星下点附近小于１ 个像元。     

利用多代卫星测高资料监测1993～2011年全球海平面变化

联合多代测高卫星T/P、Jason-1/2海面高数据和验潮站数据,确定了各卫星高度计的长期低频漂移,建立了统一的卫星测高海面高观测值,研究了1993～2011年间全球平均海平面的变化。结果表明,近18a全球平均海平面以3.12±0.4mm·a-1的速率上升,其中,海平面年际信号与ENSO事件表现出较强的相关性。     

老铁山水道漏油事故卫星监测

本文通过对ＬａｎｄｓａｔＴＭ图像和ＮＯＡＡＡＶＨＲＲ图像进行特殊的处理，将石油污染区与周围海水区别开来，利用Ｌａｎｄｓａｔ卫星精度和ＮＯＡＡ卫星时相的优势，可以确定污染面积、扩散方向和扩散速度。     

Local crossover analysis of exactly repeating satellite altimeter data

The crossover adjustment plays a central role in the processing of satellite altimeter measurements. The usual procedure is to form sea surface height differences at crossover points, solve for the radial orbit error (with due attention to the singular nature of the estimation problem) and then to construct altimetric sea-level maps using the mean sea surface heights at the crossover points. Our approach is very different, to make direct use of measurements at crossover points without differencing and to estimate simultaneously orbit parameters, mean sea surface height and sea surface height variability in a single, unified adjustment. The technique is suited for repeat data over an area small enough that adjoining passes may be considered to be parallel and to permit the solution of a set of linear equations of dimension equal to the number of crossover points. The size of the numerical problem is almost independent of the number of repeat cycles of the altimeter mission. Explicit recognition is given to the rank defect of the least-squares estimation problem; we show that, for an orbit model with r parameters, the rank defect of the local crossover problem is exactly r ². The defect may be overcome by choosing an appropriate set of constraints – either giving a best fit of mean sea surface heights to a reference surface, or minimising orbit parameters, or a minimum norm solution in which both mean sea surface heights and orbit parameters are minimised. There is no need to choose a reference pass, all passes are treated equally and data gaps are easily accommodated. Numerical results are presented for the south-western Indian Ocean, based on the first 2 years of altimeter data from the Geosat Exact Repeat Mission. Received: 31 May 1996 / Accepted: 19 April 1997     

The determination of a one year mean sea surface height track from Geosat altimeter data and ocean variability implications

One year (November 1986 to October 1987) of Geosat altimeter data with improved orbits produced at The Ohio State University has been used to define sea surface heights for 22 ERM and one year averaged Geosat track. All sea surface heights were referenced to the single reference track through the application of geoid gradient corrections. The root mean square (RMS) gradient correction was on the order of ±1 cm although it could reach 20 cm with data points in trench areas. 10 values used to form the mean were considered. Although this study was initially driven by a need for a good reference sea surface for geodetic applications the formation of the reference track yields information on the variability of the ocean surface in the first year of the Geosat ERM. The RMS point variability was ± 12.6 cm with only a very small number of values exceeding 50 cm when a depth editing criteria was used. Global plots of the sea surface variability clearly reveal the major ocean currents and their variations in position in the year. Examination of the 1° × 1° averaged sea surface height variations show average and maximum variability values as follows: Gulf Stream (29 and 50 cm); Kurshio Current (24 and 49 cm); Agulhas Current (24 and 52 cm) and the Gulf of Mexico (18 and 31 cm). These magnitudes may be dependent on the radial orbit correction procedure. To investigate this effect sea slope variations were also computed. These results also showed clear current structures but also high frequency gravity field information despite efforts to average out such information. The data described in the paper is available from the authors for numerous other studies, some of which are suggested in the paper.     

MSS93A: A new stationary sea surface combining one year upgraded ERS-1 fast delivery data and 1987 GEOSAT altimeter data

Summary At the German Processing and Archiving Facility (D-PAF) an off-line altimeter data processing system has been developed with the capability to handle different satellite missions. For combining data from different missions an algorithm has been implemented that processes a reference sea surface from altimeter data of one mission while other mission data is taken to improve the spatial resolution of this reference model. In this way a stationary sea surface (MSS93A) was computed by compiling upgraded ERS-1 fast delivery data of the first year of the 35 day repeat cycle together with Geosat T2 GDRs of 1987.MSS93A was validated by means of comparisons to external models, gradient method and visualization techniques.     

Making the connection between Geosat and TOPEX/Poseidon

We can presently construct two independent time series of sea level, each at a precision of a few centimeters, from Geosat (1985–1988) and TOPEX/Poseidon (1992–1995) collinear altimetry. Both are based on precise satellite orbits computed using a common geopotential model, JGM-2 (Nerem et al. 1994). We have attempted to connect these series using Geosat-T/P crossover differences in order to assess long-term ocean changes between these missions. Unfortunately, the observed result are large-scale sea level differences which appear to be due to a combination of geodetic and geopotential error sources. The most significant geodetic component seems to be a coordinate system bias for Geosat sea level (relative to T/P) of −7 to −9 cm in the y-axis (towards the Eastern Pacific). The Geosat-T/P sea height differences at crossovers (with JGM-2 orbits) probably also contain stationary geopotential-orbit error of about the same magnitude which also distort any oceanographic interpretation of the observed changes. We also found JGM-3 Geosat orbits have not resolved the datum errors evident from the JGM-2 Geosat -T/P results. We conclude that the direct altimetric approach to accurate determination of sea level change using Geosat and T/P data still depends on further improvement in the Geosat orbits, including definition of the geocenter. Received: 11 March 1996; Accepted: 19 September 1996     

Observation of mesoscale ocean currents from GEOSAT altimeter data

This paper discusses an altimeter data processing technique designed to compute time series of the mesoscale dynamic sea surface and to produce mean sea surfaces and surface variability. The technique has been applied to Geosat data collected over the North and South Atlantic and the South Indian Ocean. The computed mean sea surfaces show a high correlation with ocean bottom topography, whereas the variability is found to be associated with mesoscale ocean currents. High variability levels are spotted near the Gulfstream Extension and the Agulhas Return Current.Detailed examination of the sea surface and related flow field time series made it possible to identify a large number of eddies and to keep track of them in both the nort-west and south-east Atlantic. Additionally, some of the eddy characteristics have been resolved such as translation and swirl velocity. It is found that the eddy motion is affected by ocean bottom slopes.     

Coastal Gravity Anomalies from Retracked Geosat/GM Altimetry: Improvement, Limitation and the Role of Airborne Gravity Data

We process geophysical and waveform data records of the Geosat/GM (geodetic mission) satellite altimeter mission for waveform retracking and applications. An improved threshold retracker is developed. The performances of the Beta-5, threshold and improved threshold retrackers are assessed over waters around Taiwan. The improved threshold retracker outperforms the other two. The improvement in the accuracy of sea surface height (SSH) is investigated according to marine zone and the distance of waters to the shore. The improvement rate increases closer to the land, with the largest improvement rate of about 20% in waters within 10 km of the shore. Over waters around islands and coasts, there are still retracked SSHs with large errors. Least-squares collocation is used to compute gravity anomalies from the Geosat/GM altimeter data. Use of retracked SSHs improves the accuracy of gravity anomalies by about 11%. Adding airborne gravity data further improves the accuracy, especially in the immediate vicinity of the coasts. Tide model errors over coastal waters remain a problem in altimetry applications, even if the waveforms are properly retracked.     

合成孔径高度计的海面回波仿真

本文首先采用经验海浪谱仿真了三维海面,进而在此基础上根据镜面反射理论计算了海面小面元散射系数。然后根据雷达方程和合成孔径高度计的合成孔径处理计算了海面回波信号的矩阵,以便给出最终的功率回波。在目前国内外缺少相应试验数据的条件下,验证了合成孔径高度计工作方式的有效性,经验证其回波模型基本正确。同时,为了与传统高度计比较,给出了传统高度计在相同参数下的回波仿真结果。文章还通过展示合成孔径高度计在同种波高、不同涌浪和风浪方向情况下的回波,说明波浪方向也会影响合成孔径高度计回波。     

Retracking considerations in spaceborne GNSS-R altimetry

The European Space Agency Passive Reflectometry and Interferometry System In-orbit Demonstrator (IoD) aims to perform mesoscale altimetric observations by measuring the Global Navigation Satellite System (GNSS) opportunity signals reflected over the sea surface. Altimetry based on GNSS reflectometry (GNSS-R) is significantly affected by satellite motion, since it requires relatively long integration times to reduce noise. We present the impact of the satellite motion on the GNSS-R observables and the need to retrack the waveforms. By using a detailed GNSS-R space mission simulator, the change of delay difference between the direct and the reflected signals during the incoherent averaging of the waveform has been investigated. Their effects on the waveform shape and the altimetric performance are presented comparing the aligned and non-aligned waveforms. Results show that the performance of spaceborne GNSS-R altimeter is seriously degraded without a proper alignment of the waveform samples.     

QuikSCAT在台风监测中的应用

利用成熟发展的台风自身在风向上的涡旋型特征和风眼处风速相对周围较低的特征, 分别采用搜索QuikSCAT反演风向的涡旋中心, 风速或后向散射系数的台风中心区域在局部最小值点的方法定位获得了台风中心位置信息。在风矢量反演过程中采用QSCAT-1模型、借助Holland的台风模式, 修正了反演过程中的风向误差, 提高了台风中心定位的准确性。同时, 将该方法应用到对台风路径和强度监测中, 利用QuikSCAT对台风的连续观测资料分析得出台风强度和路径信息, 其中台风路径结果与美国国家飓风中心(NHC)通过最佳路径分析得到的台风路径结果基本一致, 但在台风强度结果上存在较大误差。为提高台风强度监测精度, 在反演过程中采用大风地球物理模型NN-T-GMF代替QSCAT-1模型, 使台风强度监测结果精度提高。结果表明, QuikSCAT可以有效监测海上台风路径和强度发展, 为进一步推断台风的强度发展和移动趋势提供帮助。     

联合多代测高数据建立中国海域重力异常模型

本文联合T/P数据、T/P新轨道数据、ERS数据、GFO数据、GeosatGM数据和ERS-1/168数据,用测高卫星记录点的位置信息直接计算沿轨大地水准面的方向导数,结合测线轨迹方向的方位角在交叉点处推求垂线偏差,然后利用逆Vening-Meinesz公式计算了中国近海(0o～41oN,105o～132oN)2′×2′格网分辨率的海域重力异常模型。将其与CLS_SHOW99重力异常模型比较,统计结果表示与该模型差异的RMS为8.15mgal,在剔除差值大于20mgal的点(剔除3.3%)以后,RMS为4.72mgal;与某海区船测重力异常比较的RMS为8.91mgal。