Processing math: 100%
 注册 登录 English Version
综述 | 浏览量 : 3164 下载量: 1700 CSCD: 0 更多指标
  • 导出

  • 分享

  • 收藏

  • 专辑

    • 卫星高度计融合产品的研发综述

    • Review on the research and development of merged satellite altimeter products

    • 刘磊

      123

      蒋星亮

      24

      李志锦

      2

      费建芳

      1

      吴新荣

      5

      王辉赞

      1

      叶芳

      1

      张福颖

      4

      史文丽

      6
    • 2022年26卷第12期 页码:2473-2485   

      纸质出版日期: 2022-12-07

    • DOI: 10.11834/jrs.20210241     

    扫 描 看 全 文

  • 引用

    阅读全文PDF

  • 刘磊,蒋星亮,李志锦,费建芳,吴新荣,王辉赞,叶芳,张福颖,史文丽.2022.卫星高度计融合产品的研发综述.遥感学报,26(12): 2473-2485 DOI: 10.11834/jrs.20210241.
    Liu L,Jiang X L,Li Z J,Fei J F,Wu X R,Wang H Z,Ye F,Zhang F Y and Shi W L. 2022. Review on the research and development of merged satellite altimeter products. National Remote Sensing Bulletin, 26(12):2473-2485 DOI: 10.11834/jrs.20210241.
  •  
  •  
    论文导航

    摘要

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

    Abstract

    Merged satellite altimeter products have high scientific and social values and are widely used in many ocean related fields, such as marine environmental monitoring, ocean meso-scale structure research, and ocean numerical forecast. At present, 5—6 altimeter satellites are in orbit and simultaneously publishing along-track data. We focused on the main question of extracting effective information as much as possible from the existing satellite altimeters of constant number and improving the effective resolution and product quality. The development status of the existing merged satellite altimeter products were introduced in detail.

    We introduced DT series merged products made by DUACS in detail based on the widely used AVISO products. Moreover, we analyzed the influence of different merged methods on the product quality starting from the merged method and expounded the decisive role of the selection of background field in the merged method on the effective resolution of the merged product. When combined with sentinel-3a (S3a) and J3 satellites that have recently released data and CryoSat, Saral/Altika, and HY-2a satellites, the observation density of satellite altimeters along the orbit is significantly improved. Meanwhile, the effective resolution of DT2018 produced by DUACS is not significantly improved compared with DT2014 and DT2010. The main reason is that the 7, 20, and 25 year average MSS statistical fields were selected for the background field, which includes more large-scale signals in the background error, resulting in the background error. The scale of error covariance correlation coefficient is longer. The filtering effect of the observation data is more significant, and the effective resolution cannot be effectively improved. HY-2B and HY-2C satellites were successfully launched in October 2018 and September 2020, respectively, with sentinel-3b (S3b) already in operation. The accuracy and density of satellite altimeter observation data will reach a new stage together with the SWOT altimeter satellite, which will be launched in 2021. How to consider the correlation of the same orbit observation error, accurately estimate the background error, establish a more reasonable expression function of the correlation coefficient of the background error, and provide an effective method for further improving the merged product quality and effective resolution are the directions that scientists continue to explore.

    Previous studies showed that the quality of merged satellite altimeter products is mainly affected by three aspects(1) Type of satellite selected for the production of products: the space-time accuracy of sea surface height also varies due to the difference of orbit height and period of different types of satellites; (2) Selection of altimeter standards: instrument parameters, geophysical parameters, environmental correction, construction of mean sea level, etc.; (3) Data merged methods: selection of different merged methods in the background field, difference of background error correlation coefficient scale, and processing of observation error.

    关键词

    卫星高度计; 融合产品; 有效分辨率; 变分方法; 最优估计

    Keywords

    satellite altimeter; merged data product; effective resolution; variational method; optimal estimation

    1 引言

    海面高度是气候变化最重要的指标之一,反映了气候系统对人为强迫以及自然变化的响应(

    Abraham等,2022)。20世纪初开始,通过对实况潮汐观测发现全球海平面高度GMSL(Global Mean Sea Level)持续以每年1.7±0.3 mm的速度抬升(Jevrejeva等,2008Woppelmann等,2009Church等,2011Stocker等,2013);从90年代开始,高精度的卫星高度计探测手段在海表高度场测量中的使用,为海洋分析和预报提供了全天候、高分辨率的海表高度场观测数据(Lellouche 等,2013Oke 等,2015)。

    过去二十多年,全球范围内很多研究团队都在关注于卫星高度计资料的处理并提供全球或区域尺度的海表高度时间序列产品(

    Greenslade 等,1997;Le Traon和Dibarboure,1999)。与卫星高度计的沿轨数据相比,在时间和空间尺度上具有连续性的多卫星融合产品能够直观的表现海表地转流和海洋动力过程,具有很高的实际应用价值(Koblinsky等,1992Ubelmann等,2016)。卫星高度计融合产品的质量主要受以下3个方面的影响:(1)制作产品所选择的卫星种类:因不同种类卫星的轨道高度、周期差异,导致获取的海表高度时空精度也存在不同(Pascual等,2006Dibarboure等,2011);(2)高度计标准的选择:包括高度计L2P资料处理中的仪器参数、地球物理参数、环境校正、平均海面的构造等(Dibarboure 等,2011Ablain等,2015);(3)融合处理资料的方法:不同融合方法在背景场的选择、背景误差相关系数尺度的差异以及观测误差的处理等(Ducet等,2000Dussurget等,2011Griffin和Cahill,2012Escudier等,2013)。

    针对现有卫星高度计融合产品的研发现状,本文系统地介绍了当前卫星高度计融合产品的种类,对比分析产品质量,并针对卫星高度计的融合方法,着重介绍卫星高度计融合产品质量的提升途径,进一步强调融合方法与有效分辨率的关系,为提升卫星高度计融合产品质量,提高融合产品的有效分辨率,扩大卫星高度计融合产品的应用范围提供参考。

    2 卫星高度计融合产品现状

    1991年和1992年发射的ERS-1和TOPEX/Poseidon(TP)开启了卫星高度计高精度测量的历史。随着卫星自身高度准确定位技术取得突破,TP卫星的定位准确度从米量级提高到厘米量级,最终对海面高度测量的精度达到了几厘米,从而使得卫星高度计资料在海洋中尺度涡等研究中得到广泛应用(

    Ubelmann等,2016Ablain等,2017)。表1给出了国内外高度计卫星的发展情况(http://marine.‍copernicus.‍eu/documents/QUID/CMEMS-SL-QUID-008-032-062.pdf[2020-07-03])。

    表1  高度计卫星的概况
    Table 1  Basic information about altimeter satellites
    时间高度计卫星轨道高度/km周期/d全球范围平均高程精度/cm国家或地区备注
    1973-05—1974-02 Skylab 840 23 100 美国航天局(NASA) 原理性试验卫星
    1975-04—1978-12 GEOS-3 840 23 25—50 NASA 开启卫星高度计测高时代
    1978-06—1978-10 Seasat 800 17 20—30 NASA
    1985-03—1990-01 Geosat 800 17 10—20 美国海军
    2000-01—2008-09 GFO 800 17 3.2 美国海军 Geosat的后继卫星
    1991-07—2000-03 ERS-1 785 3/35/168 3.5 欧空局(ESA)
    1995-04—2011-07 ERS-2 785 35 3.8 ESA ERS-1的后继卫星
    2002-03—2012-04 ENVISAT 782.4—799.8 30/35 2.5 ESA ERS-1/2的后继卫星
    1992-08—2006-01 TOPEX/POSEIDON(T/P) 1336 10 2.9 NASA与法国国家空间研究中心(CNES) 开启高度计高精度测量时代
    2001-12—2013-07 Jason-1(J1) 1336 10 2.9 NASA/CNES T/P的后继卫星
    2008-06—今 Jason-2(J2) 1336 10 2.9 NASA/CNES T/P高度计的第二颗后续卫星
    2016-02-17—今 Jason-3(J3) 1336 10 2.9 NASA/CNES T/P高度计的第三颗后续卫星
    2010-04—今 Cryosat-2 717 369 2.5 ESA
    2011-08—今 Haiyang-2A(HY-2A) 973 14 3.1 中国
    2013-02—今 Saral/AltiKa 800 35 2.1 印度 ENVISAT的后续卫星
    2016-06—今 Sentinel-3A(S3A) 815 27 2.4 ESA
    2018-10—今 Haiyang-2B(HY-2B) 973 14 2.5 中国
    2018-11—今 Sentinel-3B(S3B) 815 27 2.4 ESA
    icon 下载:  导出CSV icon 下载:  下载表格图片

    Le Traon等(2003)最早融合了GEOSAT后续卫星、T/P和ERS-2这3颗卫星高度计资料。2001年12月,替代T/P高度计卫星的Jason-1(简称J1;Jason-2简称J2;Jason-3简称J3)高度计卫星升空,出现了4颗卫星同时在轨运行时期。众多科学家(

    Fu等,2003;Le Traon和Dibarboure,2004;Pascual等,2006)从研究角度提出多颗卫星高度计同时在轨时,得到的融合产品质量显著提高。但上述研究处于试验分析阶段,并没有形成多源卫星高度计融合产品。

    1997年AVISO(Archiving,Validation,and Interpretation of Satellite Oceanographic data)的高度计数据融合系统DUACS(Data Unification and Altimeter Combination System)开始制作发布在国际上影响最为广泛的卫星高度计产品,包括延迟几小时到一天的准实时NRT(Near Real Time)产品和延迟若干月发布的后处理产品DT(Delayed Time)。为了尽快获取多源高度计数据,NRT产品采用非中央处理时间窗口,仅处理窗口时间前六周的高度计数据,主要应用在业务预报领域;DT产品利用所有在轨运行的卫星高度计数据,得到经交叉校准的高精度长时间序列的卫星高度计融合产品,主要应用在科学研究领域。

    2015年5月起该系统由欧洲哥白尼计划(European Copernicus Program)负责,DT2018融合产品分别由哥白尼海事局CMEMS(Copernicus Marine Service)和哥白尼气候服务机构C3S(Copernicus Climate Change Service)生产和发布。其中,CMEMS产品主要针对全球或区域尺度范围内海洋模式和海洋环流分析中的中尺度信号反演,采用所有高度计卫星数据(图1),具有更为精确的海表高度精度(

    Taburet 等,2019)。C3S产品主要用来监测海表高度的长期变化,分析气候响应因子,因此选用了具有稳定性能的两个高度计卫星数据(T/P,J1-3系列与CMEMS产品中的另一颗卫星融合),上述产品的具体差异可以参阅Taburet等(2019)表2给出了卫星高度计融合产品(全球)的概况。

    fig

    图1  DT2018采用的卫星高度计数据时间线(

    Taburet 等,2019

    Fig. 1  Timeline of modern radar altimetry missions mapped by DT2018 (

    Taburet et al., 2019

    icon 下载:  原图 | 高精图 | 低精图
    表 2  卫星高度计融合产品(全球)概况,当前在轨运行卫星标粗
    Table 2  Introduction about merged satellite altimeter products (global), and the bold are running satellites at present
    资料类型后处理产品(DT)准实时(NRT)
    产品版本 DT2010 DT2014 DT2018 当前版本
    产品制作 DUACS
    发布机构 AVISO AVISO C3S CMEMS AVISO/CMEMS
    卫星种类 T/P, J1, J2,ERS1, ERS2, ENVISAT, GFO,融合所有在轨卫星数据,标记为:all_sat_merged;融合两颗卫星数据,T/P(后续J1,J2)与ERS1(后续ERS2、ENVISAT),标记为:two_sat_merged T/P, J1 J2 , ERS-1, ERS-2, ENVISAT, GFO, Cryosat-2, Saral/AltiKa, Y-2A,融合所有在轨卫星数据,标记为:all_sat_merged;融合两颗卫星数据,T/P(后续J1,J2)与ERS1(后续ERS2、ENVISAT、Saral/AltiKa或者没有ENVISAT和Saral/AltiKa时用Cryosat-2替代),标记为:two_sat_merged T/P, J1, J2, J3中的一颗卫星与CMEMS产品中的另一颗卫星融合 T/P, Jason-1, ERS-1, ERS-2, ENVISAT, GFO, J3S3A/B, J2, SARAL/AltiKaCryosat-2HY-2A
    参考面计算周期 1993—1999 1993—2012 1993—2017
    网格分辨率 1/3°×1/3° 麦卡托网格 1/4°×1/4° 笛卡尔网格
    时间分辨率 1 d
    产品时长 1993—2010 1993—2013 1993—2020.3 1993—今
    icon 下载:  导出CSV icon 下载:  下载表格图片

    DUACS每隔4年(会适当调整时间)会依据最新的处理方法更新DT产品,本文仅介绍DT2010之后的版本。当前最新的DT2018产品由CMEMS发布,采用了累计76个使命年(commission years)共12种不同的高度计数据(

    Taburet 等,2019),制作了从1993年1月1日开始超过了26年时间序列的高度计融合产品(图1)。产品要素包括绝对动力高度ADT(Absolute Dynamic Topography)融合场和沿轨ADT数据,同时也发布了地转流的绝对融合场和异常融合场。ADT场由平均动力高度MDT(Mean Dynamic Topography)和海表高度异常SLA(Sea Level Anomaly)组成,我们将在第3节中介绍ADT场的划分原因。为提高产品质量,DT2018主要从以下3个方面进行了改进:(1)增加了DT2014未考虑的卫星种类,如在部分地区增加了2016年3月—2017年2月的HY-2A卫星资料(林明森 等,2019);(2)采用了新的高度计标准和地理修正,原始数据的选择进一步优化,最优插值参数区分了全球和区域范围;(3)相对DT2010和DT2014的SLA计算中MSS统计场分别采用7年和20年平均,DT2018进一步提升到了25年平均,提高了SLA的计算准确性(Taburet 等,2019)。

    经上述改进,采用独立沿轨数据评估后发现,DT2018产品中尺度部分的均方根误差(RMSE)相对DT2014减小了3%—4%,尤其在沿岸区域更为显著,达到了10%,地转流(Geostrophic current)的均方根误差减小了5%到10%。DT2018产品相对DT2014有了显著提升,但仍然存在以下不足:沿轨资料处理采用的低通滤波以及格点产品制作过程中使用的最优插值方法,使海表高度场上波长在200 km到65 km范围内的波动能量损失60%,即3/5的中尺度信号缺失(

    Pujol 等,2016)。中纬度地区的水平分辨尺度为1.7°,平均200 km,使得产品无法有效分辨更小尺度的中尺度信号(Chelton 等,2011和2014)。

    针对现有卫星高度计融合产品中对海表高度场中尺度信号缺失的现状,本文着重从融合方法上分析原因。

    3 高度计数据获取与处理

    在DUACS制作DT系列高度计产品的流程中,首先选定计算SLA的高度计标准,包括仪器订正,地理参数、环境校正、平均海面MSS(Mean Sea Surface,)等。随着研究深入,标准不断进行调整(见表1),之后继续完成(

    Taburet 等,2019):高度计数据获取与均质化处理、输入数据质量控制、多星交叉校准、沿轨SLA数据生成,最后生成融合产品并进行最终质量控制等多个步骤(图2)。

    fig

    图2  DUACS系统处理流程(

    Taburet 等,2019

    Fig. 2  DUACS system processing(

    Taburet et al., 2019

    icon 下载:  原图 | 高精图 | 低精图

    在高度计数据获取环节,DT产品获取自1993年起若干年份的不同卫星种类的可用高度计数据,随着新的高度计卫星升空,DT也相应的补充了部分相关数据(详见表2),然后参照各研究机构和研究团队选定高度计标准,对数据进行均质化处理(

    Taburet 等,2019)。

    在输入数据质量控制中,为了检测出高度计观测中的错误虚假值,通过改变阈值进行质量控制(

    Ablain 等,2010)。如在DT2014产品的数据处理中发现,高纬地区以及离岸20 km区域MSS分辨率不高,导致非重复轨道卫星和新到达该地区的卫星的沿轨SLA出现异常值(Pujol 等,2016)。为避免该问题出现在DT2014产品中,在这些关键区域严格筛选了错误检测标准(Pujol 等,2016)。DT2018中,MSS数据质量在沿岸减小的问题得到解决,因此,在海岸附近尽可能地保留了有效测量数据(Taburet 等,2019)。

    在数据获取环节中,从不同高度计卫星获取的高度场已经均质化处理,但由于仪器误差、轨道残差以及处理误差等存在,仍需要多星交叉校准(

    Ablain 等,2015Fu和Haines,2013)。该环节对于气候信号具有很高的重要性,需要确保TP、J1和J2卫星数据中MSS一致,之后进行多卫星间的轨道误差校正。郭靖等(2013)对Haiyang-2A(HY-2A)轨道进行校正,定轨精度可达1—2 cm。DT系列(包括DT2010、DT2014和DT2018)产品制作中的校正方法保持一致,具体可参阅Le Traon和Ogor(1998)。为了减小相邻轨道间地理位置相近引起的观测误差相关性,采用长波误差减小算法,引入长波相关误差ELWPujol 等,2016

    海表高度SSH(Sea Surface Height)为卫星相对地球重力中心轨道高度减去高度计距离海面的高度,由大地水准面高度(Geoid)和ADT组成(图3)。其中,Geoid是等重力位势面,如果海洋处于静止状态—即仅在重力场下没有洋流,它将与海洋表面相对应(

    Nerem 等,1994Rapp 等,1996)。ADT是与地转流等动力结构密切相关的海洋高度场。当海洋受到风、热交换、降水以及其他动力因素的影响时,海洋表面会偏离Geoid面(Kelly和Gille,1990;Qiu,1994)。因此,Geoid面的变化提供了海洋动力学的信息。过去,由于准确获取地球Geoid面存在较大困难,为了得到具有动力意义的高度场信息,AVISO利用高度计卫星轨道的重复特性,通过SSH减去经多年时间平均计算的MSS后(如DT2018采用25年平均),得到SLA(Taburet 等,2019)。当高度计卫星为非重复轨道时,无法采用同轨道的时间平均计算MSS,SLA的计算采用格点化的MSS作为参考面(Ophaug等,2021)。近年来,随着GRACE和GOCE等重力卫星的升空(郑伟 等,2010冯伟 等,2012),获取较高测量精度的Geoid面成为可能,MSS减去Geoid后可获得MDT,SLA与MDT相加,即可得到ADT。

    fig

    图3  高度计测高原理(引用自AVISO手册:http://marine.copernicus.eu/documents/PUM/CMEMS-SL-PUM-008-032-062.pdf)

    Fig. 3  Altimetry principle (Cited from AVISO:http://marine.copernicus.eu/documents/PUM/CMEMS-SL-PUM-008-032-062.pdf)

    icon 下载:  原图 | 高精图 | 低精图

    SLA受到观测噪音的影响,在沿轨SLA数据的制作中,采用Lanczos低通沿轨滤波减小观测噪音,并随纬度不同而调整截断波长(Le Traon和Dibarboure,1999)。DT2018沿轨产品的生成方法与DT2014保持一致,只是采用了新的高度计标准和新的卫星高度计数据等。通过信噪比分析(

    Dufau 等,2016),沿轨滤波数据的截断尺度平均为65 km,表明沿轨滤波数据可有效分辨65 km以上波长尺度的海洋结构。滤波方法降低了沿轨产品的分辨率,为了在模式和预报系统里的应用而不降低分辨率,从DT2014开始增加了沿轨非滤波产品,保留了沿轨1 Hz的短波信号。

    生成沿轨滤波数据和沿轨非滤波数据后,即进入数据融合环节。采用融合方法对多源卫星高度计沿轨数据进行插值,在规则网格上产生多源卫星高度计融合产品(

    刘磊 等,2020)。不同的融合方法对产品质量是否有影响,在具体的资料融合中该如何选择?本文给出具体分析和介绍。

    4 融合方法与背景场

    DUACS发布的DT系列融合数据采用基于高斯—马尔科夫理论的最优插值(

    Bretherton 等,1976):

    hax=hbx+Mr=1Cxr[Mq=1A-1rq(hoq-hbq)]
    (1)

    相关误差e2可表示为:

    e2=Cxx-Mr=1Mq=1CxrCxqA-1rq
    (2)

    式中,hax是规则网格x上SSH估计值,由有限个观测点rr=1,…,M)和规则网格x上的hbx融合得到,hbx取MSS。观测点上的hoq由观测真值和误差组成:hoq=htq+εoq,其中εoq为观测误差;规则网格点上的背景场hbx,由背景真值和误差组成,hbx=htx+εbx,其中εbx为背景误差。

    Arq为观测点rq点之间的协方差矩阵,Arq=<εorεoq>+<εbrεbq>,尖括号为统计平均。DUACS的融合方法中观测误差考虑了两部分,第一部分是无关项,<εorεoq>=EδrqE为观测误差方差,δrq为克罗内克函数;第二部分为相关项,考虑了卫星观测同一周期内的沿轨相关误差(

    Ducet 等,2000;Le Traon等,2003),加入了长波相关误差ELW,这样<εbrεbq>=ELW

    Cxr=<εbxεbr>为观测点和网格点之间的协方差矩阵,空间相关尺度和时间相关尺度以及信号传播速度特征采用了Arhan和de Verdière(1985)的定义:

    C(m, t)=(1+am+16(am)2-16(am)3)e-amet2/T2  
    (3)

    式中,a=3.337,t为两点间信号传播时间,T为时间相关尺度,m=(dx-CpxdtLx)2+(dy-CpydtLy)2dxdy为纬向和经向方向的空间格距,dt格点时间,LxLy为空间相关尺度半径,CpxCpy为纬向和经向传播速度。DUACS的DT2018产品融合中,背景场hbx取25年平均的MSS,并以MSS作为参考面,更多的年际变化、季节变化以及大尺度涡等信号包含在距平场SLA中作为背景误差(

    Li 等,2016)。Daley(1991)指出,更多的大尺度信号包含在背景误差中时,背景误差相关系数尺度越长。经计算得到背景误差协方差相关系数尺度LxLy在80—400 km(图4),在南海及附近区域的背景误差相关系数尺度为120 km。

    fig

    图4  DT2018产品中纬向相关长度尺度(

    Pujol 等,2016

    Fig. 4  Zonal correlation scales in DT2018(

    Pujol et al.,2016

    icon 下载:  原图 | 高精图 | 低精图

    Bretherton等(1976)的最优插值估计广泛的应用在海表温度和海表高度的资料融合中(He 等,2003Guan和Kawamura,2004Barron和Kara,2006Kawai 等,2006)。在国内,徐广珺等(2013)利用最优插值反距离加权法融合了HY-2A、J1/2以及ENVISAT卫星高度计数据,验证了HY-2A在资料融合中可替代ENVISAT数据;崔伟(2016)采用Bretherton等(1976)的最优插值方法融合了ENVISAT、J1、J2以及HY-2A等4颗同时在轨运行的卫星高度计测高数据,得到融合产品质量与DUACS产品基本一致。

    高度计卫星获取的SSH场为准地转平衡动力场,可以由线性的最优插值方法将大部分的中尺度动力信息融合网格化(

    Ducet等,2000),但非线性的小尺度信息,不能被最优插值方法识别并融合。为减小插值误差,Ubelmann等(2015和2016)提出了动力插值方法,加入准地转位涡方程,并应用在高度计沿轨资料SLA的网格化中,估计场由最优插值方法表示为:

    xa=Ns=1BHTs(HsBHTs+Rs)-1yo
    (4)
    y0=Hsx+εN
    (5)

    式中,xa是估计场,yo是观测场,x为分析场,εN为测量误差,N为观测数量,B为背景误差协方差矩阵,Rs为观测误差协方差矩阵,Hs为观测算子,将网格点投射到观测点,T为矩阵转置。由于应用在SLA融合中,因此式(4)没有列出背景场。为了考虑SSH的非线性传播,引入准地转位涡方程:

    qt+J(ψ, q)=0
    (6)
    q=2ψ-z(f20N2ψz)ψ
    (7)
    ψ=gfSSH
    (8)

    式中,q为位涡,J为雅各比算子,ψ为流函数(可以分别由SSH、SLA或ADT表示),f为科氏参数,g为重力加速度,N2为浮频率。传播算子Μ描述SSH随时间演变,可以向前或向后积分。由分析场变量xtref及不可预报场变量xtnr得到在t时次的真实状态场:

    xt=Μ(xtref)+xtnr
    (9)

    t=trefM=Ixtnr=0。当Μ对SSH演变描述的越好,xtnr越小。将式(9)代入式(5),可以得到:

    y0=Hs[M(xtref)+xtnr]+εN
    (10)

    由于M为非线性算子,无法由y0反向线性拟合得到xtref,所以给出一个初猜场,初猜误差为准线性:

    M(xtref)M(xtrefg)+M[xtref-xtrefg]
    (11)

    其中M为线性算子,将式(11)代入式(10)可得:

    y0-Hs[M(xtrefg)-Mxtrefg]=HsMxtref+Hsxtnr+εN  
    (12)

    等号左侧为已知项,右侧为线性项,估计值由最优插值方法表示为

    xtrefa=Ns=1BMTHTs(HsMBMTHTs+Cnr+Rs)-1(yo-Hs[M(xtrefg)-Mxtrefg])
    (13)

    Cnr为观测空间不可分辨信号协方差矩阵,得到与正文(1)式相似的表达。

    xtrefa=Ns=1BMTHTs(HsMBMTHTs+Cnr+Rs)-1(yo-Hs[M(xtrefg)-Mxtrefg])
    (14)

    式中,xtrefa是估计场,yo是观测场,N为观测数量,B为背景误差协方差矩阵,Rs为观测误差协方差矩阵,Hs为观测算子,将网格点投射到观测点,T为矩阵转置,Cnr为观测空间不可分辨信号协方差矩阵。传播项Μ描述SSH随时间演变,可以向前或向后积分,M为M的线性算子,xtrefg为初猜场。假设初猜场误差为准线性,可以得到与式(1)相似的表达,该方法没有实际应用在多源高度计融合产品中。此外,该方法实施中背景场仍然取多年平均的MSS,并以MSS作为参考面,由此的得到的背景误差协方差尺度阻止了融合数据分辨率的提高(

    刘磊 等,2020)。

    为进一步优化产品质量,依据观测密度选取背景场,从而确定背景误差协方差相关系数尺度的二维变分方法(2DVAR)被提出(

    刘磊 等,2020),该方法与最优插值具有一致的数理基础。根据变分原理,代价函数可以写为:

    J(hx)=12(hx-hbx)TB-1(hx-hbx)+12Ns=1(Hshx-hos)TR-1s(Hshx-hos)
    (15)

    式中,hx为分析场,hbx是背景场,hos是观测场,其他符号与前文一致。为使代价函数最小化,使δJ(hx)/δhx=0,代价函数表示为:

    δJ(hx)/δhx=B-1(hx-hbx)+Ns=1HTsR-1s(Hshx-hos)
    (16)
    0=B-1(hx-hbx)+Ns=1HTsR-1s(Hshbx-hos)+Ns=1HTsR-1s[Hs(hx-hbx)]
    (17)
    (B-1+Ns=1HTsR-1sHs)(hx-hbx)=Ns=1HTsR-1s(hos-Hshbx)
    (18)
    hx-hbx=(B-1+Ns=1HTsR-1sHs)-1Ns=1HTsR-1s(hos-Hshbx)
    (19)
    hx=hbx+(B-1+Ns=1HTsR-1sHs)-1Ns=1ΗTsR-1s(hos-Hshbx)
    (20)

    另外,由于:

    Ns=1HTsR-1s(HsBHTs+Rs)=Ns=1HTsR-1sHsBHTs+Ns=1HTs
    (21)
    Ns=1HTsR-1s(HsBHTs+Rs)=Ns=1HTsR-1sHsBHTs+B-1BNs=1HTs
    (22)
    (Ns=1HTsR-1sHs+B-1)-1Ns=1HTsR-1s=Ns=1BHTs(HsBHTs+Rs)-1
    (23)

    最终可得:

    hx=hbx+Ns=1BHTs(HsBHTs+Rs)-1(hos-Hshbx)  
    (24)

    用到的矩阵运算法则包括:δhTxAhx/δhx=(A+AT)hxA为矩阵,当R-1为对称矩阵时,有δ12hTxR-1hx/δhx=12(R-1+(R-1)T)hx=R-1hx

    基于相同的数理基础,变分方法推导的式(24)与最优插值式(1)具有相同的形式。但是,二维变分方法相对最优插值具有以下计算优点:(1)在整个网格点上进行最小化,计算效率高,可以解决高分辨率条件下格点数巨大的求解问题;(2)背景和观测误差协方差矩阵灵活多变,可更方便的简化处理,同时有利于引进动力约束;(3)依据不同的观测算子Hs,方便融合不同性质的观测资料。

    在卫星高度计的融合方法中,依据选定的背景场,构建与该背景场对应的背景误差协方差矩阵。二维变分方法计算背景误差相关系数时,以高斯函数构建背景误差协方差矩阵,用最小二乘法进行高斯分布拟合得到的背景误差相关系数尺度为约50 km(图5),远小于DT2018。融合方法中背景误差相关系数尺度是影响融合产品质量的重要因素,而背景误差相关系数尺度又决定于背景场的选择;背景场的选取决定了背景误差协方差,从而影响了背景误差协方差相关系数尺度(

    Li 等,2015和2016)。因此,DT2018的背景误差相关系数远大于选用前一天的融合产品结果作为背景场的二维变分方法(图5),其中,背景误差以相应的融合场与沿轨滤波资料相减得到。红实线为2DVAR计算得到的相关系数分布;红虚线为用最小二乘法拟合高斯函数分布(相关系数尺度为50 km);蓝实线为AVISO的DT2018融合产品25年气候平均背景场计算得到的相关系数分布(刘磊 等,2020)。

    fig

    图5  背景误差相关系数

    Fig. 5  Background error correlation coefficients

    icon 下载:  原图 | 高精图 | 低精图

    5 有效分辨率的提升途径

    融合产品的有效分辨率是融合产品质量的重要指标,即为可分辨海面涡旋的最小尺度,定义为融合产品与沿轨资料能量谱密度的比值为1/2的尺度(

    Chelton 等,2007)。在多源卫星高度计资料的融合中,如何在现有在轨运行卫星高度计数量维持不变的前提下,进一步提升有效分辨率,是众多学者关注的热点问题。

    Chelton和Schlax(2003)Chelton等(2011和2014)通过计算SSH的融合误差或融合产品与沿轨数据的能量谱密度的比值,估计了早期SLA融合产品的有效分辨率。分析结果认为中纬度的有效分辨率在2°—6°之间,主要受融合产品中的高度计数量以及采样轨道的影响。

    Ballarotta等(2019)详细介绍了DUACS的DT系列融合产品有效分辨率的变化,并得出DT2014产品在强变化区域的有效分辨率相对DT2010产品提升2%,如湾流、黑潮和南极绕流海域。有效分辨率的提升主要在于改进了高度计卫星探测技术、大气指数修正、潮汐指数修正、校正方法以及融合方法中采用更小的相关系数尺度等Pujol 等,2016Taburet 等,2019。在赤道附近有效分辨率降低可能与该区域更长的相关系数尺度有关(Pujol 等,2016)。DT2018的有效分辨率总体上与DT2014一致,但在强变化区域有2%—10%的提升,改进的原因主要在于融合方法中观测值的优化和选用更长统计周期的背景场,由此得到更为准确的背景误差(Ballarotta 等,2019)。另外,融合产品中高度计卫星选用的数量也对有效分辨率有一定影响。在DUACS产品中,融合高度计卫星的数量由两颗提升到3颗时,有效分辨率平均增加5%(Ballarotta 等,2019)。Ballarotta等(2019)认为,DUACS产品有效分辨率随卫星数量的增加而增幅不明显的原因在于采用的最优插值方法中,误差协方差参数没有变化,对融合场施加了较强的滤波效果。

    Li等(2015)指出:融合产品的有效分辨率决定于背景误差相关系数尺度,与观测信息由观测点插值到规则网格点的信息量密切相关。当背景误差相关系数尺度给定时,无论观测资料密度如何变化,融合产品的有效分辨率都不高于两倍相关系数尺度。因此,可以从以下几个方面提升融合产品的有效分辨率:(1)提升高度计观测密度;(2)考虑同轨观测误差相关性;(3)依据选取的背景场准确估算背景误差;(4)构建更为合理的背景误差相关系数表示函数,考虑不同区域的背景误差相关系数尺度变化。最终使资料融合中背景误差相关系数尺度与观测密度相匹配。在一定的观测密度分布下,背景误差相关系数尺度小,融合产品的有效分辨率高,相反,背景误差相关系数尺度大,融合产品有效分辨率低(Daley,1991)。从理论上讲,如果选用的背景场合理,包含更多的大尺度信号,从而将更少的大尺度信号归于背景误差中,可显著降低背景误差协方差相关系数尺度,有效提升多源卫星高度计融合产品有效分辨率。依据上述原则选用前一日融合场作为背景场的2DVAR方法,可显著提高融合产品的有效分辨率,见图6:AVISO(蓝)和2DVAR(红)融合产品被插值到图6(a)3(Jason-3)和图6(b)S3A 轨道(黑)后的平均沿轨能量谱密度对比;AVISO(蓝)、2DVAR(红)与图6(c)J3以及图6(d)S3A平均沿轨能量谱密度的比值。由于卫星沿轨滤波数据滤除了65 km以下波动,没有画出小于65 km尺度部分(刘磊 等,2020)。

    fig
    icon 下载:  | 高精图 | 低精图

    图6  平均沿轨能量谱密度比值

    Fig. 6  Average energy spectral density and its ratio

    6 结论

    卫星高度计融合产品广泛应用在海洋环境监测、海洋中尺度系统研究以及海洋数值预报等业务领域,具有极高的科学和社会价值,为人类更好认识海洋、研究海洋提供了新的工具。卫星高度计融合产品的研发涉及数据采集、数据均一化处理、数据质量控制、多卫星交叉校准、沿轨SLA数据生成、多源卫星融合以及最终的产品质量控制等多个流程(

    Pujol 等,2016),是众多科学家跨专业通力合作的成果。如何立足现有观测密度不变的条件下,最大程度地提取有效观测信息,以改进融合产品的有效分辨率和进一步提升产品质量,是海洋学家们研究的热点和当前亟需解决的难点问题,需对当前卫星高度计产品的研究现状进行梳理。

    本文以国际上广泛应用的AVISO产品为主线,详细介绍了DUACS制作的DT系列融合产品。从融合方法入手,分析了不同融合方法对产品质量的影响,进一步阐述了融合方法中背景场的选择对融合产品有效分辨率的决定作用。结合近期已发布数据的Sentinel-3A(S3A)和J3卫星,以及CryoSat,SARAL/AltiKa和HY-2A等卫星,卫星高度计沿轨观测密度显著提升,而DUACS制作的DT2018产品的有效分辨率相较DT2014和DT2010,并没有显著提升,主要原因在于背景场依次选用了7年、20年和25年平均的MSS统计场,将更多的大尺度信号包含于背景误差中,导致由此统计得到的背景误差协方差相关系数尺度更长,对观测资料的滤波效果更为显著,有效分辨率无法有效提升。

    随着Sentinel-3B(S3B)进入运行阶段,HY-2B和HY-2C卫星也分别于2018年10月和2020年9月成功发射(

    Guo,2018),结合即将于2021年发射的SWOT高度计卫星(徐永生 等,2017),卫星高度计观测数据的精度和密度会跨上一个新的台阶。如何在观测中考虑同轨观测误差相关性,如何准确估算背景误差,如何建立更为合理的背景误差相关系数表示函数,为进一步提升卫星高度计融合产品质量和有效分辨尺度提供有效手段及方法改进,是科学家们继续探索的方向。

    参考文献(References)

    Ablain M, Cazenave A, Larnicol G, Balmaseda M, Cipollini P, Faugère Y, Fernandes M J, Henry O, Johannessen J A, Knudsen P, Andersen O, Legeais J, Meyssignac B, Picot N, Roca M, Rudenko S, Scharffenberg M G, Stammer D, Timms G and Benveniste J. 2015. Improved sea level record over the satellite altimetry era (19932010) from the climate change initiative project. Ocean Science, 11(1): 67-82 [DOI: 10.5194/os-11-67-2015] [百度学术] 

    Ablain M, Legeais J F, Prandi P, Marcos M, Fenoglio-Marc L, Dieng H B, Benveniste J and Cazenave A. 2017. Satellite altimetry-based sea level at global and regional scales. Surveys in Geophysics, 38(1): 7-31 [DOI: 10.1007/s10712-016-9389-8] [百度学术] 

    Ablain M, Philipps S, Picot N and Bronner E. 2010. Jason-2 global statistical assessment and cross-calibration with Jason-1. Marine Geodesy, 33(S1): 162-185 [DOI: 10.1080/01490419.2010.487805] [百度学术] 

    Abraham J, Cheng L, Mann M, Trenberth K, Schuckmann K. 2022. The ocean response to climate change guides both adaptation and mitigation efforts. Atmospheric and Oceanic Science Letters. 15: 100221 [DOI: 10.1016/j.aosl.2022.100221] [百度学术] 

    Arhan M and de Verdiere A C. 1985. Dynamics of eddy motions in the eastern North Atlantic. Journal of Physical Oceanography, 15(2): 153-170 [DOI: 10.1175/1520-0485(1985)015<0153:DOEMIT>2.0.CO;2] [百度学术] 

    Ballarotta M, Ubelmann C, Pujol M I, Taburet G, Fournier F, Legeais J F, Faugère Y, Delepoulle A, Chelton D, Dibarboure G and Picot N. 2019. On the resolutions of ocean altimetry maps. Ocean Science, 15(4): 1091-1109 [DOI: 10.5194/os-15-1091-2019] [百度学术] 

    Barron C N and Kara A B. 2006. Satellite-based daily SSTs over the global ocean. Geophysical Research Letters, 33(15): L15603 [DOI: 10.1029/2006GL026356] [百度学术] 

    Bretherton F P, Davis R E and Fandry C B. 1976. A technique for objective analysis and design of oceanographic experiments applied to MODE-73. Deep Sea Research and Oceanographic Abstracts, 23(7): 559-582 [DOI: 10.1016/0011-7471(76)90001-2] [百度学术] 

    Chelton D B, Dibarboure G, Pujol I, Taburet G and Schlax M G. 2014. The spatial resolution of AVISO gridded sea surface height fields. OSTST Lake Constance, Germany, 2831 October 2014, available at: http://meetings.aviso.altimetry.fr/fileadmin/user_upload/tx_ausyclsseminar/files/29Red0900-1_OSTST_Chelton.pdf [百度学术] 

    Chelton D B and Schlax M G. 2003. The accuracies of smoothed sea surface height fields constructed from tandem satellite altimeter datasets. Journal of Atmospheric and Oceanic Technology, 20(9): 1276-1302 [DOI: 10.1175/1520-0426(2003)020<1276:TAOSSS>2.0.CO;2] [百度学术] 

    Chelton D B, Schlax M G and Samelson R M. 2011. Global observations of nonlinear mesoscale eddies. Progress in Oceanography, 91(2): 167-216 [DOI: 10.1016/j.pocean.2011.01.002] [百度学术] 

    Chelton D B, Schlax M G, Samelson R M and de Szoeke R A. 2007. Global observations of large oceanic eddies. Geophysical Research Letters, 34(15): L15606 [DOI: 10.1029/2007GL030812] [百度学术] 

    Church J A, White N J, Konikow L F, Domingues C M, Cogley J G, Rignot E, Gregory J M, Van Den broeke M R, Monaghan A J and Velicogna I. 2011. Revisiting the Earth’s sea-level and energy budgets from 1961 to 2008. Geophysical Research Letters, 38(18): L18601 [DOI: 10.1029/2011gl048794] [百度学术] 

    Cui W. 2016. The Study of Merging Sea Level Anomalies Data Derived from Multi-satellite Altimeter (Doctoral dissertation). Qingdao: First Institute of Oceanography,State Oceanic Administration,People’s Republic of China (崔伟. 2016. 多源卫星高度计海面高度异常数据融合研究(博士论文). 青岛: 国家海洋局第一海洋研究所) [百度学术] 

    Daley R. 1991. Atmospheric Data Analysis. Cambridge, UK and New York, NY, USA: Cambridge University Press471 [百度学术] 

    Dibarboure G, Pujol M I, Briol F, Le Traon P Y, Larnicol G, Picot N, Mertz F and Ablain M. 2011. Jason-2 in DUACS: updated system description, first tandem results and impact on processing and products. Marine Geodesy, 34(3/4): 214-241 [DOI: 10.1080/01490419.2011.584826] [百度学术] 

    Ducet N, Le Traon P Y and Reverdin G. 2000. Global high-resolution mapping of ocean circulation from TOPEX/Poseidon and ERS-1 and -2. Journal of Geophysical Research: Oceans, 105(C8): 19477-19498 [DOI: 10.1029/2000JC900063] [百度学术] 

    Dufau C, Orsztynowicz M, Dibarboure G, Morrow R and Le Traon P Y. 2016. Mesoscale resolution capability of altimetry: present and future. Journal of Geophysical Research: Oceans, 121(7): 4910-4927 [DOI: 10.1002/2015JC010904] [百度学术] 

    Dussurget R, Birol F, Morrow R and Mey P D. 2011. Fine resolution altimetry data for a regional application in the Bay of Biscay. Marine Geodesy, 34(3/4): 447-476 [DOI: 10.1080/01490419.2011.584835] [百度学术] 

    Escudier R, Bouffard J, Pascual A, Poulain P M and Pujol M I. 2013. Improvement of coastal and mesoscale observation from space: application to the northwestern Mediterranean Sea. Geophysical Research Letters, 40(10): 2148-2153 [DOI: 10.1002/grl.50324] [百度学术] 

    Feng W, Lemoine J M, Zhong M and Hsu T T. 2012. Terrestrial water storage changes in the Amazon basin measured by GRACE during 20022010. Chinese Journal of Geophysics, 55(3): 814-821 [百度学术] 

    冯伟, Lemoine J M, 钟敏, 许厚泽. 2012. 利用重力卫星GRACE监测亚马逊流域2002—2010年的陆地水变化. 地球物理学报, 55(3): 814-821 [DOI: 10.6038/j.issn.0001-5733.2012.03.011] [百度学术] 

    Fu L L, Stammer D, Leben R R and Chelton D B. 2003. Improved spatial resolution of ocean surface topography from the T/P-Jason-1 altimeter mission. EOS, 84(26): 241-248 [DOI: 10.1029/2003EO260002] [百度学术] 

    Fu L L and Haines B J. 2013. The challenges in long-term altimetry calibration for addressing the problem of global sea level change. Advances in Space Research, 51:1284-1300[DOI: 10.1016/j.asr.2012.06.005] [百度学术] 

    Greenslade D J M, Chelton D B and Schlax M G. 1997. The midlatitude resolution capability of sea level fields constructed from single and multiple satellite altimeter datasets. Journal of Atmospheric and Oceanic Technology, 14(4): 849-870[DOI: 10.1175/1520-0426(1997)014<0849:TMRCOS>2.0.CO;2] [百度学术] 

    Griffin D and Cahill M. 2012. Assessment of Cryosat sea level anomaly using HF radar and SST imagery. Oral presentation, OSTST Venise, Italy 2012, available at: http://www.aviso.altimetry.fr/fileadmin/documents/OSTST/2012/oral/01_thursday_27/06_NRT_applications/07_NRT_Griffin.pdf [百度学术] 

    Guan L and Kawamura H. 2004. Merging satellite infrared and microwave SSTs: methodology and evaluation of the new SST. Journal of Oceanography, 60(5): 905-912[DOI: 10.1007/s10872-005-5782-5] [百度学术] 

    Guo J, Zhao Q L, Li M and Hu Z G. 2013. Centimeter level orbit determination for HY2A using GPS data. Geomatics and Information Science of Wuhan University, 38(1): 52-55 [百度学术] 

    郭靖, 赵齐乐, 李敏, 胡志刚. 2013. 利用星载GPS观测数据确定海洋2A卫星cm级精密轨道. 武汉大学学报(信息科学版), 38(1): 52-55 [DOI: 10.13203/j.whugis2013.01.021] [百度学术] 

    Guo Q. 2018. China launches HY-2B satellite atop a LM-4B. Aerospace China, (3): 57 [百度学术] 

    He R Y, Weisberg R H, Zhang H Y, Muller-Karger F E and Helber R W. 2003. A cloud-free, satellite-derived, sea surface temperature analysis for the West Florida Shelf. Geophysical Research Letters, 30(15): 1811 [DOI: 10.1029/2003GL017673] [百度学术] 

    Jevrejeva S, Moore J C, Grinsted A and Woodworth P L. 2008. Recent global sea level acceleration started over 200 years ago?. Geophysical Research Letters, 35(8): L08715 [DOI: 10.1029/2008GL033611] [百度学术] 

    Kawai Y, Kawamura H, Takahashi S, Hosoda K, Murakami H, Kachi M and Guan L. 2006. Satellite-based high-resolution global optimum interpolation sea surface temperature data. Journal of Geophysical Research, 111(C6): C06016 [DOI: 10.1029/2005JC003313] [百度学术] 

    Kelly K A and Gille S T. 1990, Gulf Stream surface transport and statistics at 69°W from the Geosat altimeter. Journal of Geophysical Research, 95(C3): 3149-3161 [DOI: 10.1029/JC095iC03p03149] [百度学术] 

    Koblinsky C J, Gaspar P and Lagerloef G. 1992. The future of spaceborne altimetry//Oceans and Climate Change: A Long-Term Strategy(Technical Report). National Aeronautics and Space Administration (NASA), USA. N-92-26121 [百度学术] 

    Lellouche J M, Le Galloudec O, Drévillon M, Régnier C, Greiner E, Garric G, Ferry N, Desportes C, Testut C E, Bricaud C, Bourdallé-Badie R, Tranchant B, Benkiran M, Drillet Y, Daudin A and De Nicola C. 2013. Evaluation of global monitoring and forecasting systems at Mercator Océan. Ocean Science, 9(1): 57-81 [DOI: 10.5194/os-9-57-2013] [百度学术] 

    Le Traon P Y and Dibarboure G. 2004. An illustration of the contribution of the TOPEX/Poseidon-Jason-1 tandem mission to mesoscale variability studies. Marine Geodesy, 27(1/2): 3-13 [DOI: 10.1080/01490410490489313] [百度学术] 

    Le Traon P Y and Dibarboure G. 1999. Mesoscale mapping capabilities of multiple-satellite altimeter missions. Journal of Atmospheric and Oceanic Technology, 16(9): 1208-1223 [DOI: 10.1175/1520-0426(1999)016<1208:MMCOMS>2.0.CO;2] [百度学术] 

    Le Traon P Y, Faugère Y, Hernandez F, Dorandeu J, Mertz F and Ablain M. 2003. Can we merge GEOSAT Follow-On with TOPEX/Poseidon and ERS-2 for an improved description of the ocean circulation?. Journal of Atmospheric and Oceanic Technology, 20(6): 889-895 [DOI: 10.1175/1520-0426(2003)020<0889:CWMGFW>2.0.CO;2] [百度学术] 

    Le Traon P Y and Ogor F. 1998. ERS-1/2 orbit improvement using TOPEX/POSEIDON: the 2 cm challenge. Journal of Geophysical Research: Oceans, 103(C4): 8045-8057 [DOI: 10.1029/97JC01917] [百度学术] 

    Li Z J, Cheng X P, Gustafson W I Jr and Vogelmann A M. 2016. Spectral characteristics of background error covariance and multiscale data assimilation. International Journal for Numerical Methods in Fluids, 82(12): 1035-1048 [DOI: 10.1002/fld.4253] [百度学术] 

    Li Z J, McWilliams J C, Ide K and Farrara J D. 2015. A multi-scale data assimilation scheme: formulation and illustration. Monthly Weather Review, 143(9): 3804-3822 [DOI: 10.1175/MWR-D-14-00384.1] [百度学术] 

    Lin M S, He X Q, Jia Y J, Bai Y, Ye X M and Gong F. 2019. Advances in marine satellite remote sensing technology in China. Acta Oceanologica Sinica, 41(10): 99-112 [百度学术] 

    林明森, 何贤强, 贾永君, 白雁, 叶小敏, 龚芳. 2019. 中国海洋卫星遥感技术进展. 海洋学报, 41(10): 99-112 [DOI: 10.3969/j.issn.0253-4193.2019.10.007] [百度学术] 

    Liu L, Jiang X L, Fei J F and Li Z J. 2020. Development and evaluation of a new merged sea surface height product from multi-satellite altimeters. Chinese Science Bulletin, 65(18): 1888-1897 [百度学术] 

    刘磊, 蒋星亮, 费建芳, Li Zhijin. 2020. 多源卫星高度计新融合产品的研发与评估. 科学通报, 65(18): 1888-1897 [DOI: 10.1360/TB-2020-0097] [百度学术] 

    Nerem R S, Schrama E J O, Koblinsky C J, and Beckley B D. 1994. A preliminary evaluation of ocean topography from the TOPEX/Poseidon mission. Journal of Geophysical Research, 99(C12): 24565- 24583 [DOI: 10.1029/94JC01431] [百度学术] 

    Oke P R, Larnicol G, Fujii Y, Smith G C, Lea D J, Guinehut S, Remy E, Balmaseda M A, Rykova T, Surcel-Colan D, Martin M J, Sellar A A, Mulet S and Turpin V. 2015. Assessing the impact of observations on ocean forecasts and reanalyses: Part 1, Global studies. Journal of Operational Oceanography, 8(S1): s49-s62 [DOI: 10.1080/1755876X.2015.1022067] [百度学术] 

    Ophaug V, Breili K, Andersen O B. 2021. A coastal mean sea surface with associated errors in Norway based on new-generation altimetry. Advances in Space Research, 68:1103-1115[DOI: 10.1016/j.asr.2019.08.010] [百度学术] 

    Pascual A, Faugère Y, Larnicol G and Le Traon P Y. 2006. Improved description of the ocean mesoscale variability by combining four satellite altimeters. Geophysical Research Letters, 33(2): L02611 [DOI: 10.1029/2005GL024633] [百度学术] 

    Pujol M I, Faugère Y, Taburet G, Dupuy S, Pelloquin C, Ablain M and Picot N. 2016. DUACS DT2014: the new multi-mission altimeter data set reprocessed over 20 years. Ocean Science, 12(5): 1067-1090 [DOI: 10.5194/os-12-1067-2016] [百度学术] 

    Qiu B. 1994. Determining the mean Gulf Stream and its recirculations through combining hydrographic and altimetric data. Journal of Geophysical Research, 99(C1): 951-962. [DOI: 10.1029/93JC03033] [百度学术] 

    Rapp R H, Zhang C, and Yi Y. 1996. Analysis of dynamic ocean topography using T/P and orthonormal functions. Journal of Geophysical Research, 101(C10): 22583-22598. [DOI: 10.1029/96JC01987] [百度学术] 

    Stocker, T F, Qin D, Plattner G K, and Midgley P M. 2013. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK and New York, NY, USA: Cambridge University Press [百度学术] 

    Taburet G, Sanchez-Roman A, Ballarotta M, Pujol M I, Legeais J F, Fournier F, Faugere Y and Dibarboure G. 2019. DUACS DT2018: 25 years of reprocessed sea level altimetry products. Ocean Science, 15(5): 1207-1224 [DOI: 10.5194/os-15-1207-2019] [百度学术] 

    Ubelmann C, Cornuelle B and Fu L L. 2016. Dynamic mapping of along-track ocean altimetry: method and performance from observing system simulation experiments. Journal of Atmospheric and Oceanic Technology, 33(8): 1691-1699 [DOI: 10.1175/JTECHD-15-0163.1] [百度学术] 

    Ubelmann C, Klein P and Fu L L. 2015. Dynamic interpolation of sea surface height and potential applications for future high-resolution altimetry mapping. Journal of Atmospheric and Oceanic Technology, 32(1): 177-184 [DOI: 10.1175/JTECH-D-14-00152.1] [百度学术] 

    Woppelmann G, Letetrel C, Santamaria A, Bouin M N, Collilieux X, Altamimi Z, Williams S D P and Miguez B M. 2009. Rates of sea-level change over the past century in a geocentric reference frame. Geophysical Research Letters, 36(12): L12607 [DOI: 10.1029/2009gl038720] [百度学术] 

    Xu G J, Yang J S, Xu Y, Pan Y F and Chen X Y. 2013. The data fusion of the significant wave height measured with “HY-2” satellite altimeter. Acta Oceanologica Sinica, 35(4): 208-213 [百度学术] 

    徐广珺, 杨劲松, 徐圆, 潘玉方, 陈小燕. 2013. “海洋二号”有效波高数据在多源卫星高度计数据融合中的应用. 海洋学报, 35(4): 208-213 [DOI: 10.3969/j.issn.0253-4193.2013.04.024] [百度学术] 

    Xu Y S, Gao L and Zhang Y H. 2017. New generation altimetry satellite SWOT and its reference to China’s swath altimetry satellite. Remote Sensing Technology and Application, 32(1): 84-94 [百度学术] 

    徐永生, 高乐, 张云华. 2017. 美国新一代测高卫星SWOT——评述我国宽刈幅干涉卫星的发展借鉴. 遥感技术与应用, 32(1): 84-94 [DOI: 10.11873/j.issn.1004-0323.2017.1.0084] [百度学术] 

    Zheng W, Xu H Z, Zhong M and Yun M J. 2010. Research progress in international gravity satellites and future satellite gravity measurement program in China. Science of Surveying and Mapping, 35(1): 5-9 [百度学术] 

    郑伟, 许厚泽, 钟敏, 员美娟. 2010. 国际重力卫星研究进展和我国将来卫星重力测量计划. 测绘科学, 35(1): 5-9 [DOI: 10.16251/j.cnki.1009-2307.2010.01.001] [百度学术] 

    文章被引用时,请邮件提醒。
    提交

    相关作者

    暂无数据

    相关机构

    潍坊学院 物理与光电工程学院
    青岛海洋科学与技术试点国家实验室
    中国海洋大学 信息科学与工程学部海洋技术学院
    中国人民解放军海洋环境专项办公室
    中国科学技术大学
    0