中国海洋湖沼学会主办。
文章信息
- 王娟娟, 高志一, 王久珂, 李本霞, 宋金宝. 2016.
- WANG Juan-Juan, GAO Zhi-Yi, WANG Jiu-Ke, LI Ben-Xia, SONG Jin-Bao. 2016.
- Jason-2有效波高产品在中国海域的准确性评估
- VALIDATION ON JASON-2 SIGNIFICANT WAVE HEIGHT PRODUCT FOR CHINA SEAS
- 海洋与湖沼, 47(3): 509-517
- Oceanologia et Limnologia Sinica, 47(3): 509-517.
- http://dx.doi.org/10.11693/hyhz20150700180
-
文章历史
- 收稿日期:2015-07-08
- 改回日期:2015-12-01
2. 浙江大学海洋学院 杭州 310058
2. Ocean College, Zhejiang University, Hangzhou 310058, China
20世纪70年代以来,遥感观测技术和数据反演算法不断发展进步。在卫星高度计方面,自1985年美国海军的Geosat卫星(Geodetic satellite)成功发射,已经有多个卫星高度计任务顺利实施,例如: ERS-1/2(European Remote Sensing-1/2)、TOPEX/Poseidon、Jason-1/2、Envisat(Environmental Satellite)、Cryosat-2和HY-2等(Queffeulou et al,2007; Durrant et al,2009; Zieger et al,2009; Abdalla et al,2011; Jiang et al,2012; Zhang et al,2015)。目前,卫星高度计对全球海面高度的测量精度已经达 到了厘米级(Ablain et al,2010; Desai et al,2010;Ray et al,2012)。除了海面高度以外,卫星高度计还可以获取海面的风速、有效波高和潮位等信息。与传统的海洋站、浮标和船舶观测相比,卫星高度计观测克服了时间和空间上的局限性,具有全天候、高分辨率和高空间覆盖率等特点。
随着ERS-1/2、TOPEX/Poseidon和Jason-1/2等卫星数据产品的大量应用(程永存等,2008; 于振涛等,2006; 范陈清等,2014),卫星高度计已经逐渐成为物理海洋学和海洋测绘学领域的重要研究手段。特别是Jason-2卫星高度计的数据产品,发展较为成熟,在大量研究中得到应用,是高度计研究及应用技术的典范。
关于Jason-2卫星高度计产品的有效性和准确性,已经有大量研究通过与浮标观测数据(例如: 美国数据浮标中心的NDBC浮标数据)的对比分析得以证明。Abdalla等(2010)利用与浮标观测匹配的24908组数据,对Jason-2的风速和有效波高产品进行评估验证,其有效波高的覆盖范围达11m,评估得到平均偏差Bias=0.0290m,标准偏差STD=0.4119m。此外,他们还认为在Hs(有效波高)>6m时,Jason-2有效波高产品存在偏大现象。Queffeulou等(2011)利用与浮标匹配的3919组数据对Jason-2有效波高产品进行验证,得到Bias=-0.0569m,STD=0.2169m。叶小敏等(2014)利用NDBC浮标数据对2011年的Jason-2海面风速和有效波高产品进行了检验,匹配数据1635组,评估得到相关系数R=0.97,均方根误差RMSE=0.28m,STD=0.28m,Bias=-0.02m。Zhang等(2015)比较了27个月的NDBC浮标和Jason-2有效波高产品,得到RMSE=0.292m,Bias=0.016m。此外,Shanas等(2014)使用印度海域四个浮标的观测数据,评价Jason-1、Jason-2、Envisat和ERS-2四个卫星高度计的有效波高产品,分析认为它们都高估了有效波高,且在Hs<0.5m的情况下相关性极差,其中Jason-2产品与浮标观测的相关性最好。
Jason-2有效波高产品在中国海域的评估验证相对较少。Chen等(2013)用南海现场观测和NDBC浮标的数据,评估Jason-1/2和HY-2的有效波高产品,获得Jason-2有效波高产品的RMSE=0.36m。但是在南海的匹配数据仅有10组,且均为有效波高小于3m的低海况数据。Park等(2013)使用韩国气象局的海洋气象浮标数据,验证Topex/Poseidon,Jason-1,Envisat和Jason-2的有效波高产品,浮标观测有效波高的上限为4m,其位于中国黄海的浮标(34.78°N,125.77°E)观测与Jason-2有效波高产品的相关性很差,对比组114个,RMSE=1.1310m。而另一位于中国黄海的浮标(36.25°N,125.75°E),相对远离大陆和岛屿,与Jason-2匹配的对比组为88个,评估得到RMS= 0.3304m,Bias=0.1413m,R=0.9161,准确性较好。
综合上述对Jason-2有效波高产品的评估研究,将其主要结论汇总在表 1中。
参考文献 | RMSE | Bias | STD | R | N | Hs上限及结论 |
Abdalla等(2010) | 0.4119 | 0.0290 | 24908 | Hs<11m; Hs>6m时,偏大 | ||
Queffeulou等(2011) | -0.0569 | 0.2169 | 3919 | Hs<10m | ||
叶小敏等(2014) | 0.2800 | -0.0200 | 0.2800 | 0.9700 | 1635 | NDBC,Hs<10m |
Zhang等(2015) | 0.2920 | 0.0160 | NDBC,Hs<10m | |||
Chen等(2013) | 0.3600 | NDBC,Hs<7m 南海,Hs<3m | ||||
Park等(2013) | 0.3304 | 0.1413 | 0.9161 | 88 | 黄海,Hs<4m | |
1.1310 | 114 | 黄海,Hs<4m | ||||
(1)RMSE为均方根误差,Bias为平均偏差,STD为标准偏差,R为相关系数,N为匹配数据量 | ||||||
RMSE: root mean square error; Bias: the mean bias; STD: standard deviation; R: correlation coefficient; N: data amount |
从表 1中可以看出,在中国海域以外的Jason-2有效波高产品验证,匹配数据量较多,有效波高覆盖范围较广,评估得到的RMSE在0.28m至0.4119之间,准确性良好。但是,在中国海域的评估研究却非常少,只有少数几个浮标短时间观测的评估,可用的匹配数据量严重不足,且有效波高普遍在4m以下,覆盖范围较窄。此外,不同浮标获得的评估参数也存在较大差异。
为了进一步了解Jason-2有效波高产品在中国海域的准确性,我们认为可以参考TOPEX和Jason-1卫星在中国海域的评估结果。Jason-2是TOPEX和Jason-1的后续卫星,它们在多数方面基本一致(Quartly,2010),关于它们在中国海域的准确性评估也有一些研究。Ebuchi等(1994)利用JMA(Japan Meteorological Agency)在中国东海(浮标位置: 28°10′N,126°20′E。)的浮标数据,比较了TOPEX的有效波高产品,评估参数Bias=0.31m(Ku-band)和0.32m(C-band),RMSE=0.50m(Ku-band)和0.54m(C-band)。刘花等(2013)利用南海北部近岸海域3个浮标的观测数据(观测有效波高低于3m),评估分析了Jason-1的有效波高产品,评估得到的RMSE小于0.35m,相关系数大于0.95,并认为Jason-1的有效波高产品整体偏大,且与浮标波高的差异随着波高的增大而逐渐增加。
TOPEX在东海的评估RMSE在0.5m左右,Jason-2在南海的评估RMSE小于0.35m,相比于Jason-2略偏大。但是,它们所使用的匹配数据量不足,有效波高的范围较窄,整体的参考意义有限。
由于缺乏大范围、长时间和多海况覆盖的海浪观测,导致在中国海域缺乏Jason-2有效波高产品的综合评估,这使得很多研究在中国海域使用Jason-2有效波高产品时,无法详细了解其准确程度和适用范围。
作者对国家海洋局在中国海域布放的28个水文气象浮标的历史数据进行了多渠道的收集和整理,完成了繁琐的解码和质量控制工作,获得了从2002年至2014年之间的中国海域水文气象观测数据。本文采用的是从2008年至2014年之间较为完整的海浪浮标观测数据,对渤黄海、东海和南海三个分海区,以及对不同海况下的Jason-2有效波高产品准确性进行评估。本文的结果可以为Jason-2有效波高产品在中国不同海域和多种海况下的应用研究提供参考。
1 数据 1.1 Jason-2卫星高度计有效波高产品作为OSTM(Ocean Surface Topography Mission)任务TOPEX/Poseidon和Jason-1的后继卫星,Jason-2卫星于2008年6月20日在美国加利福尼亚州的范登堡空军基地成功发射。它是由CNES(Centre National d’EtudesSpatiales)、NASA、EUMETSAT(the European Organisation for the Exploitation of Meteorological Satellites)和NOAA合作完成,其上搭载了双频高度计、三频微波辐射计和高精度定位系统。海面测高精度为2.5—3.4cm,有效波高精度为0.4m或10%。该卫星的轨道高度为1336km,轨道倾角66.039°,重复周期约为10天,与TOPEX/Poseidon和Jason-1卫星基本相同。
本文选用Jason-2经过完全校正的GDR(Geo-graphical Data Record)数据,其轨道精度高,并且进行了波形重定等校正,数据准确度较高。Jason-2有效波高产品有Ku波段和C波段观测两种,由于2011年OSTM/Jason-2产品说明书中认为Ku波段观测优于C波段,因此,本文采用Ku波段数据,其数据精度为0.001m。
使用该数据之前,还需要对其进行筛选和质量控制,例如: 剔除陆地、冰面、降雨条件下以及数据可信度标志为假的数据。本文采用2011年OSTM/ Jason-2产品说明书中推荐的判据进行质量控制。
1.2 国家海洋局水文气象浮标为了满足业务化海洋预报和海洋灾害预警的需要,国家海洋局在渤黄海、东海和南海总计布放了28个常年运行的浮标,用于观测海面上的气象和水文要素。最早的浮标观测数据始于2002年,但完整性较差。从2008年以后,所有浮标的观测数据才趋于完整和连续。
浮标观测的数据格式有FUB报文和可扩展标记语言Xml两种。渤黄海数据主要为FUB格式,东海和南海早期也为FUB格式,从2011年之后陆续更改为Xml格式。FUB报文按照国际统一的报文编码规范,有效波高的数据精度为0.5m,风速精度1m/s,有效波周期精度1s,浮标位置精度0.1°,风向精度10°,气温和水温精度0.1°C。观测要素包括: 风速、风向、气温、大气压、水温、有效波高、有效波周期、最大波高,最大波周期等。
Xml格式数据精度高于FUB报文。其有效波高的精度为0.1m,风速精度0.1m/s,有效波周期精度0.1s,浮标位置精度0.01′,风向精度1°,气温精度0.1°C,水温精度0.01°C。观测要素除FUB格式包含的所有要素外,还包括: 相对湿度、平均波高、平均波周期、平均波向、十分之一波高、十分之一波周期等。
为了保证浮标业务化运行,国家海洋局会定期进行维修和养护,可能会拖动或移动浮标,此外,还由于风吹和海流的动力拖曳作用,浮标位置从长期来看有小幅变化。
2 评估分析 2.1 数据匹配Jason-2有效波高产品和浮标数据的匹配方法为: 浮标在t时刻的位置经纬度(x,y),筛选Jason-2产品中同时满足下列两个条件的数据:(a)落在以浮标位置(x,y)为圆心,以匹配距离ds为半径的圆形范围内,(b)与t时刻的绝对时间差不超过dt(小时)。根据上述的距离和时间窗口,匹配出对应浮标t时刻的所有卫星数据点,并对存在多个卫星数据点的情况下剔除超过3倍标准差的异常点。
Jason-2卫星沿着固定轨道运行,并有一定的重复周期。而浮标位置长期来看是略有变化的,其观测频率一般为1小时1次。在实际的对比分析中不可能得到与浮标观测随时同步且位置绝对重合的Jason-2数据。所以,时间和距离窗口的选择不可避免成为两者差异的原因之一。在以往类似的评估研究中,一般采用距离50km和时间30min的窗口。
相比于Jason-2数据,本文使用的浮标数据在时间上较为连续,观测间隔为1小时,因此采用30min作为时间匹配窗口比较合适。在距离窗口的选取方面,本文对不同匹配距离下的误差和数据量情况进行了综合比较,见图 1,据此来分析不同匹配距离对评估结果的影响,以得到合适的距离窗口。
从图 1中可以看出,在匹配距离为20、30和40km时,四个误差统计参数(RMSE、Bias、AE、STD)随匹配距离的变化都非常小,而数据量则逐渐增加。在匹配距离为50km时,RMSE增大非常明显,而且大误差区域内的数据量明显增多。研究中为了保证尽量多的数据,又为了兼顾匹配距离的合理性和较小的误差,选用了40km作为匹配距离,比传统的50km略小。
按照上述窗口和匹配方法,对2008年至2014年期间的浮标数据和Jason-2产品进行匹配。由于Jason-2的轨道位置固定,而浮标位置的变化不大,得到20个可用的浮标,包括: 渤黄海10个、东海6个、南海4个。最终获得Jason-2和浮标的有效波高匹配数据26060组,其中在渤黄海、东海和南海分别为11976组、8678组和5406组。
此外,本文用于验证Jason-2有效波高产品准确性的评估参数包括: 假设Jason-2产品为XJason-2,浮标观测为Xbuoy,Jason-2的测量偏差x=XJason-2-Xbuoy,则平均偏差
基于上述匹配数据,计算Jason-2有效波高产品在各个浮标位置上的评估参数。为了能够清晰查看各项评估参数的空间分布特性,图 2中展示了各个浮标位置上的四项评估参数及数据量。由于渤海的浮标布放非常密集,将其放大后在左侧子图中进行展示。
从图 2可以看出,RMSE普遍不超过0.7m,但是在东海浮标位置上的RMSE相对较大,特别是台湾海峡北端。在渤黄海的RMSE则相对较小,南海次之,都在0.3m至0.5m之间。从相对误差RE的空间分布上来看,渤黄海和台湾海峡北端浮标的RE普遍在30%左右,相对较大。相关系数R方面,渤黄海浮标的相关系数R在0.4至0.7左右,明显小于南海和东海。综合来看,南海相比于渤黄海和东海,其RMSE、相对误差RE和绝对误差AE均相对较小,相关系数R较大,这说明Jason-2有效波高产品在南海比较准确。
Jason-2在渤黄海浮标位置上的RMSE相对较小,并不能说明在渤黄海的准确性好,因为其相对误差RE和相关系数R都表现不佳。RMSE较小可能是由于该海域的匹配数据集中在较小的有效波高范围内。为此,本文在图 3中展示了三个海区和整个中国海的Jason-2、浮标有效波高匹配数据对比。从图 3中可以发现,渤黄海的浮标观测最大有效波高仅为3.5m,且大部分都集中2.5m以下。这与渤黄海特殊的地理和海洋环境有关: 水深较浅,没有强盛的大风过程影响,大浪一般是由冷空气和温带气旋引起,这导致该海域的极值浪高普遍不大,而且全年出现大浪过程的次数相比东海和南海较少。
浮标观测有效波高在4m以上的数据量较少,且主要集中在东海和南海,有效波高大于6m的数据仅有在南海的13组。东海和南海每年都会受到强台风过程的影响,因此浮标能够捕捉到较大的浪高。但是,由于浮标数据需要与Jason-2数据进行匹配,所以有效波高4m以上的匹配数据量较少。
上述对Jason-2有效波高产品在三个海区的准确性给出了常规定性分析,下面将通过评估参数值的大小进行定量分析。表 2中列出了三个海区和整个中国海的5个评估参数和匹配数据量的值,矩形的长度代表数值的大小,更易于比较。
从偏差Bias上来看,三个海区均为正值,说明Jason-2的有效波高产品偏大。东海的偏差最大,渤黄海次之,南海接近于0,这说明Jason-2有效波高在东海的偏大程度最高,在南海则基本上无偏。从图 3中Jason-2与东海浮标的对比数据来看,Jason-2与浮标差别较大的数据主要分布在有效波高大于4m的情况。而且,由于浪高较大,此区域的偏差值对整体偏差的影响也较为明显。
从绝对误差AE和RMSE上来看,东海比渤黄海和南海都要大,初步猜测也是由于有效波高4m以上的数据所引起。为了进一步分析三个海区的准确性差异,我们参考相对误差RE和相关系数R。可以发现,东海的相对误差和相关系数均要优于渤黄海。这说明东海由于有效波高4m以上的大值点较多,导致Jason-2与浮标之间的绝对型误差参数AE、RMSE、Bias较大,但是由于相对误差和相关系数的值均较优,其准确性要高于渤黄海。
三个海区中,南海的5项评估参数均为最优,充分说明Jason-2有效波高在南海的准确性最高。只是南海的数据量相比于渤黄海要少得多,且4m以上有效波高的数量也少于东海。
综上,在整个中国海域,Jason-2有效波高产品与浮标观测的综合评估情况如下: 平均偏差Bias为正值0.091m,绝对偏差AE为0.358m,RMSE为0.445m,相对误差RE为25.21%,相关系数R为0.826。参考表 1中Jason-2评估研究的误差参数值,RMSE在0.28m至1.131m之间,Bias在-0.0569m至0.1413m之间,R在0.9161至0.97之间。因此,本文研究认为Jason-2有效波高产品在中国海域的准确性较好,具有良好的代表性; 在南海的准确性最好,在渤黄海的准确性相对较差。
2.3 不同海况下的评估上述3.2节对Jason-2有效波高产品在三个海区的准确性进行了评估,但是由于各个海区的有效波高最大值存在明显差异,导致评估结果存在一定的局限性。为了评价Jason-2有效波高产品在各种海况下的准确性,本文在图 4中展示了三个评估参数随有效波高的变化情况,并区分了各个浪级。轻浪: 0m<Hs<1.3m,中浪: 1.3m≤Hs<2.5m,大浪: 2.5m≤Hs<4m,巨浪: 4m≤Hs<6m,狂浪: 6m≤Hs<9m。其中,观测有效波高最大为6.2m。
从图 4中可以看出,RMSE随着有效波高增大而增加,在Hs>5m时开始迅速增大,说明Jason-2与浮标之间的偏差程度随波高而不断增大。偏差Bias的绝对值在Hs≤5m时均小于0.2m,十分接近于0,证明Jason-2与浮标观测有效波高之间是基本无偏的。而在Hs>5m时Bias的数值开始增大,当Hs>6m时Bias已经超过1.0m且为正值,说明在Hs>6m的情况下Jason-2有效波高产品明显偏大。这与Abdalla等(2010)研究的结论类似,他们将Jason-2有效波高产品与浮标数据对比后,认为Jason-2在Hs>6m时明显偏大。相对误差RE在2m≤Hs≤5m时普遍低于25%,而在Hs=1m和6m时均大于30%。
因此,综合来看,在2m≤Hs≤5m 的情况下,RMSE、偏差Bias和相对误差RE均较小,说明Jason-2有效波高产品在此范围内准确性最好。在Hs=1m时,相对误差RE略偏大,RMSE和Bias较小,Jason-2的准确度一般。在Hs>5m,特别是Hs≥6m时,RMSE和Bias都较大,Jason-2的准确度略差。
此外,表 3展示了从轻浪到狂浪级别下的偏差Bias、平均绝对误差AE、RMSE、相对误差RE、相关系数R以及该浪级下的数据量N。与上述结论基本类似。
在轻浪条件下,RMSE、Bias和绝对误差AE虽然较小,但是相对误差RE却比其余浪级明显偏大。在中浪和大浪条件下,各项误差统计参数均较小,RMSE在0.4m至0.6m之间,相对误差RE在16%至23%之间,相关系数大于0.68,数据量超过2000组,评估参数整体表现都较优。巨浪下的各项评估参数比中浪和大浪略差,相关系数仅为0.553,数据量也仅有125组。狂浪下的各项评估参数表现最差,特别是相关系数,但是由于数据量仅为13组,统计学上的
代表意义不足。因此,我们认为Jason-2有效波高产品在中浪和大浪下的准确性最好,在轻浪和巨浪下表现略差。
3 讨论 3.1 渤黄海准确性较差的原因分析对Jason-2有效波高产品在三个分海区的评估分析,认为其在渤黄海的准确性较差,究其原因,我们认为主要有以下4个方面:
(1)浮标观测有效波高的数据精度
渤黄海浮标数据格式70%以上为FUB报文,其有效波高的精度仅为0.5m,这会导致与Jason-2有效波高的相关性变差。
(2)浮标多数靠近近岸
渤黄海的浮标距离岸线普遍较近,而卫星高度计的信号越靠近陆地,越容易受到各种杂波信号的干扰。虽然我们对Jason-2数据进行了质量控制,剔除了陆地上和异常的轨道点,但是仍然不能完全排除陆地对卫星高度计数据的影响。
(3)水深较浅
渤黄海的水深较浅,特别是渤海水深仅有十几米,这导致波浪破碎、白冠、浪-流相互作用等非线性动力效应较强。这可能会导致海面反射回卫星高度计的回波信号发生异常,从而影响有效波高产品的准确度。
(4)有效波高以轻浪为主
Jason-2有效波高在不同海况下的评估分析表明其在轻浪条件下准确性稍差。而渤黄海浮标的观测有效波高50%以上都为轻浪,这也会导致其整体的准确性变差。
3.2 离岸距离对Jason-2准确性的影响为了研究远离陆地的距离和水深对Jason-2有效波高产品准确性的影响,我们利用与Jason-2匹配的浮标所在的离岸距离和水深,分析Jason-2的准确性随离岸距离和水深的变化情况。考虑到各个浮标的有效波高分布范围差异较大,采用绝对型的评估参数RMSE、Bias和AE等可能会对分析结果带来错误的判断,因此选择了相对型的评估参数: 相关系数和相对误差。评估参数随离岸距离的变化情况见图 5。
从图 5中可以看出,相关系数R随着离岸距离而逐渐增大,在最远的离岸距离279km位置上,相关系数达到了0.94,而在离岸距离低于100km时,相关系数可以低至0.66。相对误差随着离岸距离的增大而逐渐减小,在最远的离岸距离位置上,相对误差为17.8%,在离岸距离低于100km时,相对误差可以高达51.5%。因此,我们认为Jason-2有效波高产品的准确性随着离岸距离增加会显著提高。而且,从卫星高度计测高的原理上来说,陆地对卫星信号的污染效应是十分明显的,远离陆地的海域卫星高度计信号的可信度会更高。
Jason-2评估参数随水深的变化曲线与图 5大致类似,但是数据点较为分散,规律性不强,准确性随水深增加的特征不是很明显,因此本文没有进行展示和详细分析。
此外,在对Jason-2的评估分析中,匹配数据中巨浪和狂浪条件下的数据量分别为125组和13组,且有效波高的分布无法覆盖狂浪的全部范围,考虑到样本数太少带来的偶然误差影响,在一定程度上会降低巨浪和狂浪条件下的分析结果可靠性。
4 结论基于上述分析,总结得到Jason-2有效波高产品在中国海域的准确性评估结论如下:
(1)Jason-2有效波高产品在中国海域具有较好的准确性。偏差Bias=0.091m,绝对误差AE=0.358m,均方根误差RMSE=0.445m,相对误差RE=25.21%,相关系数R=0.826。
(2)Jason-2有效波高产品在南海的准确性最好,而在渤黄海的准确性较差,主要原因可能是数据精度较低、靠近近岸和水深较浅等。
(3)在不同海况下,Jason-2有效波高产品在2m≤Hs≤5m 的情况下准确性较好,在Hs>5m时准确性略差。
(4)Jason-2有效波高产品的准确性随着离岸距离的增加会显著提高。
此外需要注意,所有的匹配数据组中,有效波高观测值绝大多数都在4m以下,虽然符合近海的一般海况条件,但是缺乏高海况下的对比数据。因此,本文的评估分析结果有一定的适用局限性。