• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    干涉雷達(dá)高度計(jì)定標(biāo)檢驗(yàn)進(jìn)展

    2020-11-12 09:44:02周興華楊磊徐永生朱
    海洋科學(xué)進(jìn)展 2020年4期
    關(guān)鍵詞:高度計(jì)中尺度波數(shù)

    周興華楊 磊徐永生朱 琳

    (1.自然資源部 第一海洋研究所,山東 青島266061;2.山東科技大學(xué) 海洋工程技術(shù)學(xué)院,山東 青島266590;3.中國科學(xué)院 海洋研究所,山東 青島266071;4.山東科技大學(xué) 測繪科學(xué)與技術(shù)學(xué)院,山東 青島266590)

    經(jīng)歷30余a的發(fā)展,衛(wèi)星測高獲得的海面動力學(xué)資料已成為物理海洋學(xué)、大地測量學(xué)等學(xué)科的重要基石[1-2]。至今全球已發(fā)射了數(shù)十顆不同系列高度計(jì)衛(wèi)星,從早期的SEASAT[3]、Geosat[4]到現(xiàn)在的Jason-3[5]、Sentinel3[6]等,數(shù)據(jù)觀測精度從最初的米級提高到現(xiàn)在的3 cm 左右[2],2011年和2017年我國分別發(fā)射的HY-2A 和HY-2B高度計(jì)衛(wèi)星使我國具備了衛(wèi)星海洋動力觀測的業(yè)務(wù)化能力,HY-2A/B 全球測高資料在國際測高領(lǐng)域發(fā)揮了重要的作用[7-12]。由于不同系列衛(wèi)星高度計(jì)的測量偏差各有差異,同一顆衛(wèi)星高度計(jì)的測量偏差也可能隨儀器老化等原因而產(chǎn)生漂移,這就需要開展持續(xù)性的定標(biāo)檢驗(yàn)工作[13-15]。

    定標(biāo)檢驗(yàn)(Calibration Validation)是監(jiān)測衛(wèi)星高度計(jì)測量數(shù)據(jù)偏差和系統(tǒng)漂移的基本手段,是保障業(yè)務(wù)化衛(wèi)星資料精度和質(zhì)量長期一致性和連續(xù)性的關(guān)鍵[16-18]。定標(biāo)檢驗(yàn)通常是利用地面實(shí)測真值評估衛(wèi)星測量的精度和質(zhì)量,地面實(shí)測真值通常又包含:1)具有參考橢球基準(zhǔn)的驗(yàn)潮站[19];2)GNSS浮標(biāo)[20-23];3)有源定標(biāo)器[24-26];4)全球驗(yàn)潮站網(wǎng)、海底壓力傳感器、Argo網(wǎng)[27-30]①LEGEAIS J F.Validation of altimeter data by comparison with in-situ T/S argo profiles,2012.等。上述地面實(shí)測數(shù)據(jù)中以統(tǒng)一基準(zhǔn)的驗(yàn)潮站和GNSS浮標(biāo)使用最為廣泛,而有源定標(biāo)器等方法的應(yīng)用研究偏少[13]。

    目前國外主要的定標(biāo)場:美國NASA 運(yùn)行的Harvest石油平臺定標(biāo)場[31-33]、法國CNES運(yùn)行的Corsica定標(biāo)場[6,34-38]、希臘克里特理工大學(xué)運(yùn)行的Crete定標(biāo)場[39-45]、澳大利亞塔斯馬尼亞大學(xué)負(fù)責(zé)的Bass海峽定標(biāo)場[46-48]。中國建設(shè)的衛(wèi)星高度計(jì)定標(biāo)站點(diǎn)有千里巖定標(biāo)場[49-52]和珠海萬山綜合定標(biāo)場等[13,53-55]。這些定標(biāo)場最長的觀測近30余a,為傳統(tǒng)雷達(dá)高度計(jì)提供服務(wù),對于探測衛(wèi)星測高偏差,統(tǒng)一衛(wèi)星測高數(shù)據(jù)基準(zhǔn)發(fā)揮了重要作用[16]。

    傳統(tǒng)高度計(jì)的噪聲水平被限定在100 km 以上的空間分辨率[56-57],為了以更高的分辨率觀測海洋現(xiàn)象,美國NASA 和法國CNES于2009年正式提出SWOT 干涉雷達(dá)高度計(jì)計(jì)劃②LEE L F,DOUG L A,ERENSTO R,et al.The SWOT(Surface Water and Ocean Topography)mission:spaceborne radar interferometry for oceanographic and hydrological applications,2009.,經(jīng)過10余a的論證和發(fā)展,SWOT 計(jì)劃于2022年發(fā)射[58]。該衛(wèi)星將利用干涉雷達(dá)測量技術(shù),獲得120 km 寬度刈幅的海面高度信息,使高度計(jì)觀測中尺度、亞中尺度海洋現(xiàn)象成為可能,對于進(jìn)一步研究海洋和全球氣候變化的互相作用具有重大意義[59]。同時該衛(wèi)星還將觀測高分辨率的陸地水表面高度,對研究氣候變化和全球陸地水分布、儲量變化的相互作用意義重大[60]。然而目前現(xiàn)有的定標(biāo)檢驗(yàn)技術(shù)均服務(wù)于傳統(tǒng)星下點(diǎn)高度計(jì),其點(diǎn)對點(diǎn)式的定標(biāo)方法很難適用于具有百余公里寬度刈幅的干涉雷達(dá)高度計(jì)[61]。如何利用地面實(shí)測資料實(shí)現(xiàn)干涉雷達(dá)高度計(jì)的定標(biāo)檢驗(yàn)已經(jīng)成為國內(nèi)外學(xué)者所關(guān)注的重點(diǎn)問題之一[59,62]。因此,NASA 研制了AirSWOT 機(jī)載干涉雷達(dá)高度計(jì)系統(tǒng),并進(jìn)行了相關(guān)定標(biāo)檢驗(yàn)試驗(yàn)[63-64],開發(fā)了SWOT simulator模擬程序以及SWOT OSSE定標(biāo)檢驗(yàn)?zāi)M系統(tǒng)[61,65-66]③GAULTIER L,UBELMANN C,FU L L.SWOT simulator documentation,2017.,進(jìn)行SWOT 短波分量的定標(biāo)檢驗(yàn)?zāi)M工作;我國也開展了相應(yīng)的預(yù)研工作,其中“天宮”二號搭載全球首個成像雷達(dá)高度計(jì)獲取了高分辨率的寬刈幅海面高度[67],我國也開展了機(jī)載干涉雷達(dá)高度計(jì)飛行試驗(yàn),對衛(wèi)星載荷原型進(jìn)行了性能檢驗(yàn)和評估[62,68]。

    1 SWOT定標(biāo)檢驗(yàn)技術(shù)

    圖1 SWOT 干涉衛(wèi)星高度計(jì)工作原理以及誤差波數(shù)譜需求④DANIEL E F.SWOT project mission performance and error budget revision A:JPL D-79084,2017.Fig.1 Conceptual illustration of the SWOT mission measurement and the sea surface height error spectrum requirement as a function of wavenumber④

    SWOT 主要載荷包括Ka波段的雷達(dá)干涉計(jì)(KaRIN),它可以通過干涉測量方式獲取海洋和陸地水的二維高度信息,此外還有一個傳統(tǒng)的剖面高度計(jì)載荷用于星下點(diǎn)測量,在這種測量模式下,SWOT 可以覆蓋寬20~120 km、平行于基線方向的范圍,其網(wǎng)格像素大小在10~70 m 范圍,在垂直于基線方向可以獲得高達(dá)2.5 m 的理論分辨率④DANIEL E F.SWOT project mission performance and error budget revision A:JPL D-79084,2017.(圖1)。在7.5 km×7.5 km的范圍內(nèi),通過濾波可以使SWOT噪聲水平降低到2cm2/(cycle·km-1),而Jason系列的噪聲水平約為100 cm2/(cycle·km-1)[69],因此SWOT 對海洋信號的分辨能力與傳統(tǒng)高度計(jì)比較提高一個數(shù)量級,其對海洋信號的目標(biāo)分辨率高達(dá)15 km[65]。針對SWOT二維海面高度特征以及全新的科學(xué)目標(biāo),其定標(biāo)檢驗(yàn)也有別于傳統(tǒng)高度計(jì),例如傳統(tǒng)高度計(jì)的定標(biāo)檢驗(yàn)?zāi)康闹饕獮楂@取測高偏差,并以均值和均方根誤差的形式進(jìn)行評估,而SWOT 的定標(biāo)檢驗(yàn)開創(chuàng)性地將高度計(jì)誤差分解為波數(shù)譜的形式[69],并設(shè)定2個目標(biāo):大地測量檢驗(yàn)和海洋學(xué)定標(biāo)檢驗(yàn)。其中大地測量檢驗(yàn)主要評估SSH(Sea Surface Height)是否滿足SWOT 波數(shù)譜的科學(xué)需求,主要可用的設(shè)備包括GNSS和LiDAR;海洋學(xué)檢驗(yàn)主要是評估SWOT 滿足既定科學(xué)目標(biāo)的能力,如亞中尺度海洋信號的探測等,可用的現(xiàn)場測量設(shè)備有溫鹽潛標(biāo)、Glider、海洋上層剖面儀等⑤⑥FU L L.Ocean in-situ Cal Val,2018.。此外SWOT 還觀測陸地水域,因此其定標(biāo)檢驗(yàn)也在河流湖泊開展,目前國外已經(jīng)利用機(jī)載干涉雷達(dá)高度計(jì)和GNSS浮毯、無人船、LiDAR 等開展了河流觀測和定標(biāo)檢驗(yàn)⑦⑧FROIDEVAL L,LAIGNEL B,PICOT N.Cal Val LiDAR,2018.。

    1.1 海洋長波信號定標(biāo)

    結(jié)合SWOT 的科學(xué)目標(biāo),其海洋領(lǐng)域的定標(biāo)檢驗(yàn)需在15~1000 km 波長范圍內(nèi)開展,其中短波部分(15~120 km)需使用現(xiàn)場觀測資料進(jìn)行定標(biāo)檢驗(yàn),而長波部分(120~1000 km)則使用SWOT 攜帶的Jason型傳統(tǒng)高度計(jì)[59,70]。研究表明,Jason級的傳統(tǒng)高度計(jì)在大于100 km 波長范圍海洋信號強(qiáng)于噪聲,在SWOT 刈幅寬度120 km 處的信噪比為5,具備在長波部分評估SWOT 的條件[59],但Jason系列高度計(jì)無法識別小于70 km 波長的海洋信號[71](圖2)。Wang等詳細(xì)研究了SWOT 傳統(tǒng)高度計(jì)評估干涉雷達(dá)高度計(jì)長波信號的能力,考慮了SWOT 誤差與沿軌距離的相關(guān)性,選擇刈幅的內(nèi)邊緣(距離星下點(diǎn)最近)和中間線(儀器噪聲最小)分別與傳統(tǒng)高度計(jì)比較波數(shù)譜,通過模擬實(shí)驗(yàn)證明傳統(tǒng)高度計(jì)可評估最短為120 km 波長的SWOT 數(shù)據(jù)[70],因此僅波長小于120 km 的海洋信號需要通過現(xiàn)場觀測進(jìn)行檢驗(yàn)。

    圖2 Jason-2和Jason-1串聯(lián)任務(wù)觀測的波數(shù)譜以及SWOT 傳統(tǒng)高度計(jì)對干涉高度計(jì)長波信號的定標(biāo)波長概率分布函數(shù)Fig.2 The sea surface height wavenumber spectra of Jason-1 and Jason-2 during the tandem mission and the probability function of the globally sampled nadir scale for SWOT

    1.2 海洋短波信號檢驗(yàn)

    SWOT 設(shè)計(jì)的主要科學(xué)目標(biāo)之一為通過海面地形觀測海洋亞中尺度現(xiàn)象,因此SWOT 短波部分的定標(biāo)檢驗(yàn)非常重要[61,66]。由于SWOT 攜帶的傳統(tǒng)高度計(jì)可在長波部分對SSH 進(jìn)行評估,故地面觀測主要用于SWOT 短波部分(15~120 km)的定標(biāo)檢驗(yàn)。為了開展短波信號的定標(biāo)檢驗(yàn),SWOT 設(shè)計(jì)了總計(jì)90d、重訪周期1 d的快速Cal/Val飛行軌道,和傳統(tǒng)高度計(jì)比較,該設(shè)計(jì)可大大提高時間分辨率,提高地面固定定標(biāo)場的定標(biāo)檢驗(yàn)效率[69]。短波信號的定標(biāo)檢驗(yàn)同樣采用波數(shù)譜的形式,可利用的定標(biāo)檢驗(yàn)數(shù)據(jù)包括機(jī)載LiDAR、機(jī)載干涉高度計(jì)和海上實(shí)測資料等。目前除了機(jī)載數(shù)據(jù),尚無海面現(xiàn)場觀測資料可以提供15~120 km 波長的可靠波數(shù)譜,通過現(xiàn)場觀測資料評估SWOT 短波信號是一項(xiàng)具有挑戰(zhàn)性的工作[69]。Wang等研究了利用固定站位Glider、溫鹽錨系定標(biāo)檢驗(yàn)SWOT 短波信號的可行性,為使現(xiàn)場觀測資料能滿足15 km波長的信號解析能力,站位空間間隔至少為7.5 km,為滿足150 km 的最大波長信號解析能力,需沿著SWOT 軌跡布設(shè)20個觀測站位[61]。Wang等結(jié)合傳統(tǒng)高度計(jì)解析長波信號的能力修改了短波信號的上限波長為120 km[70],因此現(xiàn)場站位數(shù)可降低為16個。此外,NASA 關(guān)于GNSS浮標(biāo)用于SWOT 短波定標(biāo)檢驗(yàn)的可行性正在評估之中[69]⑨HAINES B,DESAI S,MEINIG C,et al.Status,results and plans for development of GPS buoys:potential for SWOT in-situ Cal-Val,2018.。然而成像雷達(dá)高度計(jì)獲得的直接海面高度信息為相對于參考橢球的SSH,如與Glider/溫鹽錨系結(jié)果比較,需要轉(zhuǎn)化成動力高度SLA,因此定標(biāo)結(jié)果中將包含大地水準(zhǔn)面誤差。目前大地水準(zhǔn)面在短波分量的誤差尚難以確定[72-73],所以在實(shí)際的海洋定標(biāo)過程中可能產(chǎn)生由大地水準(zhǔn)面誤差引起的潛在問題。

    圖3 SWOT 和現(xiàn)場設(shè)備模擬的波數(shù)譜誤差[61]Fig.3 Simulated wavenumber spectra of SWOT and field facilities[61]

    因機(jī)載干涉雷達(dá)高度計(jì)存在顯著的wave-bunching效應(yīng),機(jī)載LiDAR 將替代AirSWOT 進(jìn)行SWOT短波信號的定標(biāo)檢驗(yàn)[74]。美國Scripps海洋研究所開發(fā)了一套MASS(Modular Aerial Sensing System)系統(tǒng),主要攜帶機(jī)載地形測量LiDAR、視頻、紅外高光譜成像儀等設(shè)備,可以直接觀測SSH、海浪譜和海流等參數(shù)[74]。因?yàn)楹Q蟊砻嬲鎸?shí)值難以確定,所以MASS選擇不變的地面目標(biāo)且與海洋波動類似的沙丘進(jìn)行了系統(tǒng)噪聲檢驗(yàn),結(jié)果如圖4所示,MASS的噪聲水平在傳統(tǒng)高度計(jì)之下,符合SWOT 需求。然而MASS尚未在海洋領(lǐng)域獲得可靠的結(jié)論,主要原因是缺少海上實(shí)測資料支持15~120 km 的準(zhǔn)確波數(shù)譜觀測,相關(guān)研究仍在進(jìn)行之中。

    圖4 MASS在沙漠中觀測的高度波數(shù)譜和SWOT 目標(biāo)誤差⑩Fig.4 MASS lidar height error spectra obtained over the dunes,with the SWOT ocean requirements and typical nadir altimetry signal levels overlaid⑩ LENAIN L,MELVILLE K,STATOM N,et al.Broad range airborne ocean topography measurements:Modular Aerial Sensing System(MASS)in support of SWOT Calval,2019.

    1.3 AirSWOT河流定標(biāo)檢驗(yàn)

    SWOT 陸地水觀測的科學(xué)目標(biāo)為直接觀測地表水面高度,并進(jìn)一步獲得水儲量、流量、洪水動力等參數(shù),陸地水觀測的誤差分配(error budget)包含高度誤差、坡度誤差和水分類誤差,SWOT設(shè)計(jì)的陸地水表面高度測量目標(biāo)為1 km2面積平均后的誤差不大于10 cm,坡度測量目標(biāo)為10 km 長度誤差不大于1 cm/km,對面積大于250 m2的水體或者寬度大于100 m 的河流識別誤差小于15%。和傳統(tǒng)高度計(jì)比較SWOT 可以觀測不低于70 m 分辨率的二維陸地水文信息,通過觀測水面坡度可估計(jì)得到水流速度、流量等信息,對于研究水文動力學(xué)意義重大[75-78]。

    圖5 AirSWOT 河流測量實(shí)驗(yàn)區(qū)域及WSE觀測的定標(biāo)檢驗(yàn)結(jié)果[63]Fig.5 Location of the Tanana River study by AirSWOT,and the Cal/Val results from the WSE measurement[63]

    為驗(yàn)證SWOT 所設(shè)定的科學(xué)目標(biāo)能否實(shí)現(xiàn),美國科學(xué)家利用AirSWOT 干涉雷達(dá)高度計(jì)進(jìn)行了大量的河流觀測定標(biāo)檢驗(yàn)試驗(yàn)[63-64]⑩LENAIN L,MELVILLE K,STATOM N,et al.Broad range airborne ocean topography measurements:Modular Aerial Sensing System(MASS)in support of SWOT Calval,2019.。AirSWOT 河流定標(biāo)檢驗(yàn)地點(diǎn)之一是美國阿拉斯加州的Tanana河,通過機(jī)載干涉雷達(dá)高度計(jì)得到了河流的水面高度(Water Surface Elevation,WSE)和水面坡度兩個參數(shù),地面驗(yàn)證數(shù)據(jù)采用船載測繪型GPS和河底壓力計(jì)。該試驗(yàn)結(jié)果表明,AirSWOT觀測到的WSE誤差均方根為9.0 cm(1 km 平均),坡度均方根誤差為1.0 cm/km(10 km 長度),符合SWOT 科學(xué)目標(biāo)的精度需求,因此未來可使用AirSWOT 開展SWOT 陸地水域的現(xiàn)場定標(biāo)檢驗(yàn)工作[63]。此外,AirSWOT 在美國俄勒岡州的Willamette河也開展了定標(biāo)檢驗(yàn)試驗(yàn),其WSE誤差均方根為11.6 cm,坡度誤差均方根為3.2 cm/km[64],與文獻(xiàn)[63]有所不同,高出SWOT 的目標(biāo)誤差值,這也表明SWOT 的地面定標(biāo)檢驗(yàn)仍舊存在挑戰(zhàn)和不確定性。

    2 中國干涉雷達(dá)高度計(jì)定標(biāo)檢驗(yàn)技術(shù)

    目前我國多家機(jī)構(gòu)也計(jì)劃研制干涉雷達(dá)高度計(jì),其中“觀瀾”衛(wèi)星是一顆由青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室牽頭研發(fā)的干涉雷達(dá)高度計(jì)衛(wèi)星,除干涉雷達(dá)高度計(jì)外,“觀瀾”還計(jì)劃搭載激光雷達(dá)進(jìn)行海洋的剖面觀測,“觀瀾”干涉雷達(dá)高度計(jì)的主要科學(xué)目標(biāo)與SWOT 類似,也是海洋中尺度和亞中尺度海洋信號的觀測[79]。但“觀瀾”干涉雷達(dá)高度計(jì)設(shè)計(jì)參數(shù)與SWOT 略有不同,其設(shè)計(jì)的刈幅寬度為166.4 km,觀測信息的像素大小為5 km,設(shè)定的目標(biāo)是識別10 km 尺度的亞中尺度海洋信息。為了驗(yàn)證“觀瀾”干涉雷達(dá)高度計(jì)的概念設(shè)計(jì),目前已開展機(jī)載海洋觀測進(jìn)行定標(biāo)試驗(yàn)[62]。

    2.1 機(jī)載系統(tǒng)定標(biāo)檢驗(yàn)

    機(jī)載干涉雷達(dá)高度計(jì)系統(tǒng)(Airborne Interferometric Radar Altimeter System,AIRAS)設(shè)計(jì)波段為Ka波段,刈幅覆蓋范圍為1 km(3000 m 高度,1°~15°入射角),預(yù)期測量精度為3.4 cm。2019年3月自然資源部第一海洋研究所等在青島小麥島海域進(jìn)行了機(jī)載觀測和現(xiàn)場定標(biāo)檢驗(yàn)試驗(yàn),現(xiàn)場觀測設(shè)備為2套同類型的GNSS測高浮標(biāo),同時搜集了小麥島海洋觀測站采集的驗(yàn)潮、氣象、波浪等資料,通過GNSS靜態(tài)測量和水準(zhǔn)測量,統(tǒng)一了小麥島驗(yàn)潮站和機(jī)載資料的參考基準(zhǔn)面。如何利用有限的GNSS浮標(biāo)對機(jī)載干涉雷達(dá)高度計(jì)的二維海面高度數(shù)據(jù)進(jìn)行評估是一項(xiàng)前人所未涉及的研究,由于GNSS僅提供固定點(diǎn)位的SSH 時間序列(1 Hz),而機(jī)載干涉雷達(dá)高度計(jì)則測量二維SSH(原始分辨率0.3 m),因此無法使用直接比較的辦法進(jìn)行誤差分析。對此,提出一種適用于GNSS浮標(biāo)評估二維海面高度的方法,對于時空域中波浪引起的海平面高程(Wind-induced Sea Surface Elevation,WSSE)的均勻場,假定以下3個陳述正確[62]:

    1)如果被測WSSE的時間序列足夠長,則在域內(nèi)任何位置測得的WSSE的方差都相同;

    2)如果被測WSSE的面積足夠大,則域內(nèi)隨時測量的WSSE的方差相同;

    3)如果條件1)和2)為真,則波浪引起的WSSE在空間和時間上是均勻的場,并且1)和2)的方差相等。

    進(jìn)而通過功率密度譜分析可以得到空間域的波數(shù)譜和時間域的頻率譜,通過積分可以獲得海洋波動信號方差,GNSS浮標(biāo)和機(jī)載數(shù)據(jù)的在時空域的積分方差差異可視為二者的總體差異。然而時空域的波數(shù)譜采用了不同的數(shù)學(xué)單位,因此無法直接分析誤差在單獨(dú)的空間域或者時間域的分布特征。為此,Yang等[62]提出一種時空譜統(tǒng)一方法,通過頻散關(guān)系和方差守恒關(guān)系可以得到:

    式中:f為頻率;k為波數(shù);S和Q分別為時間域和空間域的能量密度函數(shù)。

    通過上述轉(zhuǎn)換關(guān)系可統(tǒng)一GNSS浮標(biāo)一維時間頻率譜和機(jī)載干涉雷達(dá)高度計(jì)二維空間波數(shù)譜,使之在同樣的坐標(biāo)單位下進(jìn)行比較。結(jié)果表明,機(jī)載干涉雷達(dá)高度計(jì)和GNSS浮標(biāo)的方差差異小于8 cm2,對應(yīng)約2.8 cm 的標(biāo)準(zhǔn)偏差,二者波數(shù)譜的相關(guān)性大于0.9,表明二者在頻率域的方差分布具有很高的一致性。

    圖6 機(jī)機(jī)載干涉雷達(dá)高度計(jì)和GNSS浮標(biāo)統(tǒng)一后風(fēng)浪引起的海面高度波數(shù)譜比較[64]Fig.6 Unified VPSin frequency domain of Wind-induced Sea Surface Elevation measured by AIRAS and GNSS buoys[64]

    2.2 衛(wèi)星定標(biāo)檢驗(yàn)計(jì)劃

    我國“觀瀾”干涉雷達(dá)高度計(jì)設(shè)計(jì)的定標(biāo)方案同SWOT 相似,同樣采用波數(shù)譜的方式進(jìn)行誤差評估。目前我國在西太平洋、印度洋、南海等海域布設(shè)了溫鹽潛標(biāo),這些已有的設(shè)施可為我國相關(guān)衛(wèi)星的定標(biāo)檢驗(yàn)提供潛在的數(shù)據(jù)支持,衛(wèi)星發(fā)射后的定標(biāo)檢驗(yàn)計(jì)劃還包括布設(shè)GNSS 浮標(biāo)陣列、PIES(Pressure-sensorequipped Inverted Echo Sounders)陣列,以及在南海石油平臺安裝GNSS、驗(yàn)潮儀定標(biāo)設(shè)備等[79]。

    圖7 “觀瀾”任務(wù)的定標(biāo)示意圖[79]Fig.7 A schematic illustration of calibration/validation facility for the′Guanlan′mission[79]

    3 其他方法

    除現(xiàn)場觀測和機(jī)載激光測量方法,SWOT 二維海面高度的誤差評估還可以采用交叉譜分析方法[80]。由于SWOT 的觀測數(shù)據(jù)誤差源多樣(如地球物理誤差、儀器誤差等),并且測量誤差和距離向的長度相關(guān),如何定量分離各項(xiàng)誤差成為一項(xiàng)挑戰(zhàn)性的研究。Ubelmann等研究表明利用交叉譜分析可以探測雷達(dá)距離向的絕大多數(shù)誤差并給出能量密度譜,對于天線滾動誤差、相位誤差、基線膨脹誤差具有非常高的精度[80]。由于SWOT 的誤差評估采用的波數(shù)譜的策略,因此交叉譜誤差分析的方法對于SWOT 的定標(biāo)檢驗(yàn)具有一定的適用性。此外,該方法還可用于SWOT 數(shù)據(jù)去噪及計(jì)算誤差協(xié)方差等方面,對于SWOT 網(wǎng)格數(shù)據(jù)產(chǎn)品的制作有一定價值。交叉定標(biāo)也是SWOT 干涉雷達(dá)高度計(jì)數(shù)據(jù)定標(biāo)檢驗(yàn)的方法之一[81-82],Dibarboure等研究表明利用經(jīng)驗(yàn)交叉定標(biāo)方法可以探測并顯著降低SWOT 天線的姿態(tài)誤差,例如由天線基線長度誤差和翻滾誤差帶來的最大可為70 cm 的誤差通過交叉定標(biāo)可降低為2 cm[82]。

    4 挑戰(zhàn)和展望

    SWOT 及其同類型干涉雷達(dá)高度計(jì)的科學(xué)目標(biāo)之一是通過二維海面高度的觀測來解析海洋中尺度和亞中尺度信號,因此設(shè)定的目標(biāo)分辨率為10~15 km,較傳統(tǒng)高度計(jì)高出一個數(shù)量級。干涉雷達(dá)高度計(jì)通過InSAR 技術(shù)進(jìn)行海面高觀測,與傳統(tǒng)高度計(jì)測量方式不同,其測量方式和科學(xué)目標(biāo)使得定標(biāo)檢驗(yàn)具有一定的挑戰(zhàn)性,且尚無成熟經(jīng)驗(yàn)可以借鑒。其主要挑戰(zhàn)可簡述為以下幾方面:1)衛(wèi)星干涉雷達(dá)高度計(jì)觀測亞中尺度海洋信號本身即一項(xiàng)具有挑戰(zhàn)性的工作,實(shí)現(xiàn)此目標(biāo)應(yīng)首先將干涉雷達(dá)高度計(jì)儀器噪聲水平控制在約2 cm2/(cycle·km),較現(xiàn)在傳統(tǒng)雷達(dá)高度計(jì)的噪聲水平需提高2個數(shù)量級,因此對儀器硬件噪聲控制以及數(shù)據(jù)濾波算法都提出了很高要求,特別是在利用機(jī)載干涉雷達(dá)高度計(jì)標(biāo)定衛(wèi)星數(shù)據(jù)時,機(jī)載資料的噪聲水平不應(yīng)高于衛(wèi)星的噪聲水平,同理對浮標(biāo)、PIES等現(xiàn)場設(shè)備也存在噪聲水平的控制問題;2)干涉雷達(dá)高度計(jì)定標(biāo)的主要目標(biāo)之一是獲得噪聲誤差譜,分析其在不同尺度(特別是亞中尺度)的誤差水平,平均海平面模型(MSS)或者大地水準(zhǔn)面模型(Geoid)是海洋動力信號提取的重要參考基準(zhǔn),MSS包含了大地水準(zhǔn)面高度和穩(wěn)態(tài)海面地形(Mean Dynamic Topography,MDT)高度,MSS是海面高度異常(SLA)的參考基準(zhǔn),大地水準(zhǔn)面是絕對動力高度(Absolute Dynamic Topography,ADT)的參考基準(zhǔn)[83],目前常用的MSS模型為多年測高平均模型(如DTU MSS 2018),大地水準(zhǔn)面模型為理論模型EGM2008等,其中多年測高數(shù)據(jù)平均得到的MSS模型忽略了可能存在的亞中尺度的季節(jié)性變化,其引入的誤差可能對分離較弱的亞中尺度海洋信號帶來干擾,進(jìn)而影響誤差譜分析,因此解決MSS基準(zhǔn)模型在亞中尺度波段的精度評估和優(yōu)化格外重要[72]。

    本文對SWOT 發(fā)展至今10余年內(nèi)的國外定標(biāo)檢驗(yàn)工作和研究成果進(jìn)行了總結(jié)分析,并對我國相關(guān)的同類型干涉雷達(dá)高度計(jì)的機(jī)載定標(biāo)檢驗(yàn)工作做了闡述。由于SWOT 和我國的同類衛(wèi)星都尚未發(fā)射,目前絕大部分工作都是基于理論模擬數(shù)據(jù)開展,在真實(shí)海洋中的相關(guān)機(jī)載定標(biāo)檢驗(yàn)工作依舊缺乏。目前美國NASA/JPL證實(shí)AirSWOT 尚不具備海洋亞中尺度信號的定標(biāo)檢驗(yàn)?zāi)芰?而轉(zhuǎn)為采用機(jī)載激光測高系統(tǒng)MASS,這也給我國的相關(guān)機(jī)載試驗(yàn)帶來啟示,需盡早開展機(jī)載干涉雷達(dá)高度計(jì)的海洋觀測性能分析,并對海洋干涉雷達(dá)高度計(jì)數(shù)據(jù)處理算法進(jìn)行優(yōu)化設(shè)計(jì),論證是否可以滿足相關(guān)要求。SWOT 和我國同類衛(wèi)星的科學(xué)目標(biāo)之一還包括陸地水域測高,因此陸地水域的定標(biāo)檢驗(yàn)也是重要工作之一。陸地河、湖、濕地等資源是寶貴的自然資源,通過新型干涉雷達(dá)高度計(jì)觀測其水位、儲量和流量意義重大,我國在開展海洋領(lǐng)域定標(biāo)檢驗(yàn)研究同時,建議同步考慮陸地水域的測高定標(biāo)研究工作。

    參考文獻(xiàn)(References):

    [1]XU Y S,GAO L,ZHANG Y H.New generation altimetry satellite SWOT and its reference to China′s swath altimetry satellite[J].Remote Sensing Technology and Application,2017,32(1):84-94.徐永生,高樂,張?jiān)迫A.美國新一代測高衛(wèi)星SWOT 及對我國寬刈幅干涉衛(wèi)星的發(fā)展借鑒[J].遙感技術(shù)與應(yīng)用,2017,32(1):84-94.

    [2]CHELTON D,RIES J,HAINES B,et al.Satellite altimetry and earth sciences:a handbook of techniques and applications[M].California:Academic Press,2001.

    [3]TAPLEY B D,BORN G H,PARKE M E.The SEASAT altimeter data and its accuracy assessment[J].Journal of Geophysical Research:Oceans,1982,87(C5):3179-3188.DOI:10.1029/JC087iC05p03179.

    [4]CHENEY R E,DOUGLAS B C,MILLER L.Evaluation of Geosat altimeter data with application to tropical Pacific sea level variability[J].Journal of Geophysical Research:Oceans,1989,94(C4):4737-4747.DOI:10.1029/JC094iC04p04737.

    [5]BIANCAMARIA S,SCHAEDELE T,BLUMSTEIN D,et al.Validation of Jason-3 tracking modes over French rivers[J].Remote Sensing of Environment,2018(209):77-89.DOI:10.1016/j.rse.2018.02.037.

    [6]BONNEFOND P,LAURAIN O,EXERTIER P,et al.Calibrating the SAR SSH of Sentinel-3A and CryoSat-2 over the Corsica facilities[J/OL].Remote Sensing,2018,10(1):92.http:∥www.mdpi.com/2072-4292/10/1/92.DOI:10.3390/rs10010092.

    [7]CHEN N,HAN G,YANG J,et al.Hurricane Sandy storm surges observed by HY-2A satellite altimetry and tide gauges[J].Journal of Geophysical Research:Oceans,2014,119(7):4542-4548.DOI:10.1002/2013jc009782.

    [8]LI M,ZHAO C,ZHAO Y,et al.Polar sea ice monitoring using HY-2A scatterometer measurements[J].Remote Sensing,2016,8(8):688.DOI:10.3390/rs8080688.

    [9]JIANG X W,LIN M S,ZHANG Y G.An overview of HY-2 satellite ground application system[J].Strategic Study of CAE,2014,16(6):4-12.蔣興偉,林明森,張有廣.HY-2 衛(wèi)星地面應(yīng)用系統(tǒng)綜述[J].中國工程科學(xué),2014,16(6):4-12.

    [10]YANG L,ZHOU X H,WANG Z Y,et al.Absolute calibration of sea surface height for HY-2A satellite altimeter by GNSS buoy and tide gauge data[J].Haiyang Xuebao,2017,39(1):111-120.楊磊,周興華,王朝陽,等.基于GNSS浮標(biāo)和驗(yàn)潮資料的HY-2A 衛(wèi)星雷達(dá)高度計(jì)絕對定標(biāo)[J].海洋學(xué)報,2017,39(1):111-120.

    [11]PENG H L,LIN M S,MU B,et al.Global statistical evaluation and performance analysis of HY-2A satellite radar altimeter data[J].Haiyang Xuebao,2015,32(7):54-66.彭海龍,林明森,穆博,等.HY-2A 衛(wèi)星雷達(dá)高度計(jì)數(shù)據(jù)的全球統(tǒng)計(jì)評價及質(zhì)量分析[J].海洋學(xué)報,2015,32(7):54-66.

    [12]ZHU C,GUO J,HWANG C,et al.How HY-2A/GM altimeter performs in marine gravity derivation:assessment in the South China Sea[J].Geophysical Journal International,2019,219(2):1056-1064.DOI:10.1093/gji/ggz330.

    [13]YANG L,ZHOU X H,XU Q J,et al.Current status of satellite altimeter calibration[J].Journal of Remote Sensing,2019,23(3):392-407.楊磊,周興華,徐全軍,等.衛(wèi)星高度計(jì)定標(biāo)現(xiàn)狀[J].遙感學(xué)報,2019,23(3):392-407.DOI:10.11834/jrs.20198262.

    [14]ANDERSEN O B,CHENG Y.Long term changes of altimeter range and geophysical corrections at altimetry calibration sites[J].Advances in Space Research,2013,51(8):1468-1477.DOI:10.1016/j.asr.2012.11.027.

    [15]BROWN S,DESAI S,KEIHM S,et al.Microwave radiometer calibration on decadal time scales using on-earth brightness temperature references:application to the TOPEX microwave radiometer[J].Journal of Atmospheric and Oceanic Technology,2009,26(12):2579-2591.DOI:10.1175/2009jtecha1305.1.

    [16]FU L L,HAINES B J.The challenges in long-term altimetry calibration for addressing the problem of global sea level change[J].Advances in Space Research,2013,51(8):1284-1300.

    [17]BONNEFOND P,HAINES B J,WATSON C.In situ absolute calibration and validation:a link from coastal to open-ocean altimetry[M/OL].VIGNUDELLI S,KOSTIANOY A G,CIPOLLINI P.∥Coastal Altimetry.Springer Berlin Heidelberg,2011:259-296.http:∥dx.doi.org/10.1007/978-3-642-12796-0_11.DOI:10.1007/978-3-642-12796-0_11.

    [18]LEULIETTE E W,NEREM R S,MITCHUM G T.Calibration of TOPEX/Poseidon and Jason altimeter data to construct a continuous record of mean sea level change[J].Marine Geodesy,2004,27(1-2):79-94.

    [19]HAINES B J,DONG D,BORN G H,et al.The harvest experiment:monitoring Jason-1 and TOPEX/POSEIDON from a California offshore platform special issue:Jason-1 calibration/validation[J].Marine Geodesy,2003,26(3-4):239-259.DOI:10.1080/714044520.

    [20]XU X,XU K,SHEN H,et al.Sea surface height and significant wave height calibration methodology by a GNSS buoy campaign for HY-2A altimeter[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2016,9(11):5252-5261.

    [21]BORN G H,PARKE M E,AXELRAD P,et al.Calibration of the TOPEX altimeter using a GPS buoy[J].Journal of Geophysical Research:Oceans,1994,99(C12):24517-24526.DOI:10.1029/94jc00920.

    [22]CHADWELL C D,BOCK Y.Direct estimation of absolute precipitable water in oceanic regions by GPS tracking of a coastal buoy[J].Geophysical Research Letters,2001,28(19):3701-3704.

    [23]WATSON C S.Satellite altimeter calibration and validation using GPS buoy technology[D].Tasmania:University of Tasmania,2005.

    [24]WAN J,GUO W,ZHAO F,et al.HY-2A radar altimeter ultrastable oscillator drift estimation using reconstructive transponder with its validation by multimission cross calibration[J].IEEE Transactions on Geoscience and Remote Sensing,2015(9):5229-5236.

    [25]ZHEN G A O.Simulations of satellite-borne altimeter tracking system with transponder in calibration test[J].International Journal of Remote Sensing,1993,14(2):365-382.DOI:10.1080/01431169308904342.

    [26]WANG C,GUO W,ZHAO F,et al.Development of the reconstructive transponder for in-orbit calibration of HY-2A altimeter[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2015,9(6):2709-2719.

    [27]CHAMBERS D,RIES J,SHUM C,et al.On the use of tide gauges to determine altimeter drift[J].Journal of Geophysical Research:Oceans(1978-2012),1998,103(C6):12885-12890.

    [28]FENOGLIO-MARC L,SCH?NE T,ILLIGNER J,et al.Sea level change and vertical motion from satellite altimetry,tide gauges and GPS in the indonesian region[J].Marine Geodesy,2012,35(Suppl.1):137-150.DOI:10.1080/01490419.2012.718682.

    [29]VALLADEAU G,LEGEAIS J F,ABLAIN M,et al.Comparing altimetry with tide gauges and Argo profiling floats for data quality assessment and mean sea level studies[J].Marine Geodesy,2012,35(Suppl.1):42-60.DOI:10.1080/01490419.2012.718226.

    [30]GE L,XU Y S,YIN B S.Comparison of sea surface height anomalies derived by pressure-sensor-equipped inverted echo sounder[J].Marine Sciences,2016,40(2):128-135.葛磊,徐永生,尹寶樹.日本海西南海域現(xiàn)場觀測和衛(wèi)星高度計(jì)獲取的海面高度距平的比較研究[J].海洋科學(xué),2016,40(2):128-135.

    [31]CHRISTENSEN E J,HAINES B J,KEIHM SJ,et al.Calibration of TOPEX/POSEIDON at platform harvest[J].Journal of Geophysical Research:Oceans,1994,99(C12):24465-24485.DOI:10.1029/94jc01641.

    [32]HAINES B J,DESAI S D,BORN G H.The harvest experiment:calibration of the climate data record from TOPEX/Poseidon,Jason-1 and the Ocean Surface Topography Mission[J].Marine Geodesy,2010,33(Suppl.1):91-113.DOI:10.1080/01490419.2010.491028.

    [33]MNARD Y,JEANSOU E,VINCENT P.Calibration of the TOPEX/POSEIDON altimeters at Lampedusa:additional results at harvest[J].Journal of Geophysical Research:Oceans(1978-2012),1994,99(C12):24487-24504.

    [34]BONNEFOND P,EXERTIER P,MENARD Y,et al.Calibration of radar altimeters and validation of orbit determination in the Corsica-Capraia area[J].Molecular Human Reproduction,1997,15(414):205-213.

    [35]BONNEFOND P,EXERTIER P,LAURAIN O,et al.Absolute Calibration of Jason-1 and TOPEX/Poseidon altimeters in corsica special issue:Jason-1 calibration/validation[J].Marine Geodesy,2003,26(3-4):261-284.DOI:10.1080/714044521.

    [36]BONNEFOND P,EXERTIER P,LAURAIN O,et al.Absolute calibration of Jason-1 and Jason-2 altimeters in Corsica during the formation flight phase[J].Marine Geodesy,2010,33(Suppl.1):80-90.

    [37]BONNEFOND P,EXERTIER P,LAURAIN O,et al.SARAL/AltiKa absolute calibration from the multi-mission corsica facilities[J].Marine Geodesy,2015,38(Suppl.1):171-192.

    [38]GOURINE B,KAHLOUCHE S,EXERTIER P,et al.Corsica SLR positioning campaigns(2002 and 2005)for satellite altimeter calibration missions[J].Marine Geodesy,2008,31(2):103-116.DOI:10.1080/01490410802053658.

    [39]MERTIKAS S P,ZHOU X,QIAO F,et al.First preliminary results for the absolute calibration of the Chinese HY-2 altimetric mission using the CRS1 calibration facilities in West Crete,Greece[J].Advances in Space Research,2015,57(1):78-95.DOI:10.1016/j.asr.2015.10.016.

    [40]MERTIKAS S P,DASKALAKIS A,TZIAVOS I N,et al.Ascending and descending passes for the determination of the altimeter bias of Jason satellites using the gavdos facility[J].Marine Geodesy,2011,34(3-4):261-276.DOI:10.1080/01490419.2011.584837.

    [41]MERTIKAS S P,DASKALAKIS A,TZIAVOS I N,et al.First calibration results for the SARAL/AltiKa altimetric mission using the gavdos permanent facilities[J].Marine Geodesy,2015(Suppl.1):249-259.

    [42]MERTIKAS S P,DASKALAKIS A,TZIAVOS I N,et al.Altimetry,bathymetry and geoid variations at the Gavdos permanent Cal/Val facility[J].Advances in Space Research,2013,51(8):1418-1437.DOI:10.1016/j.asr.2012.10.021.

    [43]PAVLIS E C,MERTIKASS P.The GAVDOS mean sea level and altimeter calibration facility:results for Jason-1[J].Marine Geodesy,2004,27(3-4):631-655.DOI:10.1080/01490410490902106.

    [44]TZIAVOS I N,VERGOS G S,MERTIKAS S P,et al.The contribution of local gravimetric geoid models to the calibration of satellite altimetry data and an outlook of the latest GOCE GGM performance in Gavdos[J].Advances in Space Research,2013,51(8):1502-1522.DOI:10.1016/j.asr.2012.06.013.

    [45]MERTIKAS,DONLON,VUILLEUMIER,et al.An action plan towards fiducial reference measurements for satellite altimetry[J].Remote Sensing,2019,11(17):1993.DOI:10.3390/rs11171993.

    [46]WATSON C,WHITE N,COLEMAN R,et al.TOPEX/Poseidon and Jason-1:absolute calibration in Bass Strait,Australia[J].Marine Geodesy,2004,27(1-2):107-131.DOI:10.1080/01490410490465373.

    [47]WATSON C,WHITE N,CHURCH J,et al.Absolute calibration in bass strait,Australia:TOPEX,Jason-1 and OSTM/Jason-2[J].Marine Geodesy,2011,34(3-4):242-260.DOI:10.1080/01490419.2011.584834.

    [48]WATSON C,COLEMAN R,WHITE N,et al.Absolute calibration of TOPEX/Poseidon and Jason-1 using GPS buoys in Bass Strait,Australia special issue:Jason-1 calibration/validation[J].Marine Geodesy,2003,26(3-4):285-304.DOI:10.1080/714044522.

    [49]LEI Y,XINGHUA Z,P.M S,et al.First calibration results of Jason-2 and Saral/AltiKa satellite altimeters from the Qianliyan permanent facilities[J].Advances in Space Research,2017,59(12):2831-2842.

    [50]ZHOU X,YANG L,LIN M,et al.Absolute calibration of HY-2,Jason-2 and Saral/AltiKa from China in-situ calibration site:Qian Li Yan[C]∥2015 IEEE International Geoscience and Remote Sensing Symposium(IGARSS).DOI:10.1109/IGARSS.2015.7326616.

    [51]ZHOU X,YANG L,LEI N,et al.Calibration results of multiple satellite altimetry missions from QianliYan permanent CAL/VAL facilities[C]∥2017 IEEE International Geoscience and Remote Sensing Symposium(IGARSS).DOI:10.1109/IGARSS.2017.8127392.

    [52]ZHOU X,YANG L,WANG Y,et al.Research progress of satellite altimeter calibration in China[C]∥IGARSS 2019-2019 IEEE International Geoscience and Remote Sensing Symposium.DOI:10.1109/IGARSS.2019.8898571.

    [53]ZHANG Y F,JIANG X W,MA C F.The progress of the altimeter calibration based on the in-situ observation[J].Remote Sensing Technology and Application,2019,34(6):1190-1196.張宇飛,蔣興偉,馬超飛.基于現(xiàn)場觀測的星載雷達(dá)高度計(jì)定標(biāo)方法進(jìn)展[J].遙感技術(shù)與應(yīng)用,2019,34(6):1190-1196.

    [54]ZHAN Y G,LIN M S.Research progress of calibration site and calibration method of satellite altimeter[J].Marine Science Bulletin,2007,26(3):87-92.張有廣,林明森.衛(wèi)星高度計(jì)海上定標(biāo)場及定標(biāo)方法研究進(jìn)展[J].海洋通報,2007,26(3):87-92.

    [55]LIU Y L.Research on HY-2 radar altimeter sea height calibration technology[D].Qingdao:Ocean University of China,2014.劉亞龍.HY-2雷達(dá)高度計(jì)海面高度定標(biāo)技術(shù)研究[D].青島:中國海洋大學(xué),2014.

    [56]XU Y,FU L L.The effects of altimeter instrument noise on the estimation of the wave number spectrum of sea surface height[J].Journal of Physical Oceanography,2012,42:2229-2233.

    [57]XU Y,FU L L.Global variability of the wavenumber spectrum of oceanic mesoscale turbulence[J].Journal of Physical Oceanography,2011,41:802-809.

    [58]CHEVALIER L,DESROCHES D,LAIGNEL B,et al.High-resolution SWOT simulations of the macrotidal seine estuary in different hydrodynamic conditions[J].IEEE Geoscience and Remote Sensing Letters,2019,16(1):5-9.DOI:10.1109/LGRS.2018.2862470.

    [59]FU L L,UBELMANN C.On the transition from profile altimeter to Swath altimeter for observing global ocean surface topography[J].Journal of Atmospheric and Oceanic Technology,2014,31(2):560-568.DOI:10.1175/JTECH-D-13-00109.1.

    [60]PRIGENT C,LETTENMAIER D P,AIRES F,et al.Toward a high-resolution monitoring of continental surface water extent and dynamics,at global scale:from GIEMS(Global Inundation Extent from Multi-Satellites)to SWOT(Surface Water Ocean Topography)[M/OL]∥CAZENAVE A,CHAMPOLLION N,BENVENISTE J.Remote Sensing and Water Resources.Cham:Springer International Publishing,2016:149-165[2018-10-17].https:∥doi.org/10.1007/978-3-319-32449-4_7.DOI:10.1007/978-3-319-32449-4_7.

    [61]WANG J,FU L L,QIU B,et al.An observing system simulation experiment for the calibration and validation of the surface water ocean topography sea surface height measurement using in situ platforms[J].Journal of Atmospheric and Oceanic Technology,2017,35(2):281-297.DOI:10.1175/JTECH-D-17-0076.1.

    [62]YANG L,XU Y,ZHOU X H,et al.Calibration of an airborne interferometric radar altimeter over the Qingdao coast sea,China:10[J].Remote Sensing,2020,12(10):1651.DOI:10.3390/rs12101651.

    [63]ALTENAU E H,PAVELSKY T M,MOLLER D,et al.AirSWOT measurements of river water surface elevation and slope:Tanana river,AK[J].Geophysical Research Letters,2017,44(1):181-189.DOI:10.1002/2016GL071577.

    [64]TUOZZOLO S,LIND G,OVERSTREET B,et al.Estimating river discharge with Swath altimetry:a proof of concept using AirSWOT observations[J].Geophysical Research Letters,2019,46(3):1459-1466.DOI:10.1029/2018GL080771.

    [65]WANG J,FU L L,TORRES H S,et al.On the spatial scales to be resolved by the surface water and ocean topography Ka-band radar interferometer[J].Journal of Atmospheric and Oceanic Technology,2019,36(1):87-99.DOI:10.1175/JTECH-D-18-0119.1.

    [66]LI Z,WANG J,FU L L.An observing system simulation experiment for ocean state estimation to assess the performance of the SWOT mission:part 1 a twin experiment[J].Journal of Geophysical Research:Oceans,2019,124(7):4838-4855.DOI:10.1029/2018JC014869.

    [67]ZHANG Y,SHI X,WANG H,et al.Interferometric imaging radar altimeter on board Chinese Tiangong-2 space laboratory[C]∥2018 Asia-Pacific Microwave Conference(APMC).DOI:10.23919/APMC.2018.8617189.

    [68]CHEN H,GUO W,YANG S B,et al.A method of calibration of SAR altimeter[J].Journal of Electronics&Information Technology,2017,39(9):2232-2237.諶華,郭偉,楊雙寶,等.一種合成孔徑雷達(dá)高度計(jì)定標(biāo)方法[J].電子與信息學(xué)報,2017,39(9):2232-2237

    [69]MORROW R,FU L-L,ARDHUIN F,et al.Global observations of fine-scale ocean surface topography with the surface water and ocean topography(SWOT)mission[J/OL].Frontiers in Marine Science,2019,6[2020-04-04].https:∥www.frontiersin.org/articles/10.3389/fmars.2019.00232/full.DOI:10.3389/fmars.2019.00232.

    [70]WANG J,FU L L.On the long-wavelength validation of the SWOT KaRIn measurement[J].Journal of Atmospheric and Oceanic Technology,2019,36(5):843-848.DOI:10.1175/JTECH-D-18-0148.1.

    [71]DUFAU C,ORSZTYNOWICZ M,DIBARBOURE G,et al.Mesoscale resolution capability of altimetry:present and future[J].Journal of Geophysical Research:Oceans,2016,121(7):4910-4927.DOI:10.1002/2015JC010904.

    [72]PUJOL M I,SCHAEFFER P,FAUGRE Y,et al.Gauging the improvement of recent mean sea surface models:a new approach for identifying and quantifying their errors[J].Journal of Geophysical Research:Oceans,2018,123(8):5889-5911.DOI:10.1029/2017JC013503.

    [73]DIBARBOURE G,PUJOL M I.Improving the quality of sentinel-3A data with a Hybrid Mean Sea Surface Model,and implications for sentinel-3B and SWOT[J/OL].Advances in Space Research,2019[2020-04-12].http:∥www.sciencedirect.com/science/article/pii/S027311771930434X.DOI:10.1016/j.asr.2019.06.018.

    [74]MELVILLE W K,LENAIN L,CAYAN D R,et al.The modular aerial sensing system[J].Journal of Atmospheric and Oceanic Technology,2016,33(6):1169-1184.DOI:10.1175/JTECH-D-15-0067.1.

    [75]BARATELLI F,FLIPO N,RIVIèRE A,et al.Retrieving river baseflow from SWOT spaceborne mission[J].Remote Sensing of Environment,2018,218:44-54.DOI:10.1016/j.rse.2018.09.013.

    [76]GRIPPA M,ROUZIES C,BIANCAMARIA S,et al.Potential of SWOT for monitoring water volumes in Sahelian Ponds and Lakes[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2019:1-9.DOI:10.1109/JSTARS.2019.2901434.

    [77]DOMENEGHETTI A,SCHUMANN G J-P,FRASSON R P M,et al.Characterizing water surface elevation under different flow conditions for the upcoming SWOT mission[J].Journal of Hydrology,2018,561:848-861.DOI:10.1016/j.jhydrol.2018.04.046.

    [78]BIANCAMARIA S,DURAND M,ANDREADIS K M,et al.Assimilation of virtual wide swath altimetry to improve Arctic river modeling[J].Remote Sensing of Environment,2011,115(2):373-381.DOI:10.1016/j.rse.2010.09.008.

    [79]CHEN G,TANG J,ZHAO C,et al.Concept design of the"Guanlan"science mission:China′s novel contribution to space oceanography[J].Frontiers in Marine Science,2019,6.

    [80]UBELMANN C,DIBARBOURE G,DUBOIS P.A cross-spectral approach to measure the error budget of the SWOT altimetry mission over the ocean[J].Journal of Atmospheric and Oceanic Technology,2017,35(4):845-857.DOI:10.1175/JTECH-D-17-0061.1.

    [81]DIBARBOURE G,UBELMANN C.Investigating the performance of four empirical cross-calibration methods for the proposed SWOT mission[J].Remote Sensing,2014,6(6):4831-4869.

    [82]DIBARBOURE G,LABROUE S,ABLAIN M,et al.Empirical cross-calibration of coherent SWOT errors using external references and the altimetry constellation[J].IEEE Transactions on Geoscience and Remote Sensing,2012,50(6):2325-2344.DOI:10.1109/TGRS.2011.2171976.

    [83]XIE L,ZHENG Q,TIAN J,et al.Cruise observation of Rossby waves with finite wavelengths propagating from the Pacific to the South China Sea[J].Journal of Physical Oceanography,2016,46(10):2897-2913.DOI:10.1175/JPO-D-16-0071.1.

    猜你喜歡
    高度計(jì)中尺度波數(shù)
    聲場波數(shù)積分截斷波數(shù)自適應(yīng)選取方法
    一種基于SOM神經(jīng)網(wǎng)絡(luò)中藥材分類識別系統(tǒng)
    電子測試(2022年16期)2022-10-17 09:32:26
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計(jì)有效波高研究
    海洋通報(2021年3期)2021-08-14 02:20:46
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運(yùn)動研究
    海洋通報(2020年5期)2021-01-14 09:26:52
    基于深度學(xué)習(xí)的中尺度渦檢測技術(shù)及其在聲場中的應(yīng)用
    MIMU/GNSS/ODO/高度計(jì)/航姿儀組合導(dǎo)航微系統(tǒng)硬件設(shè)計(jì)
    航天控制(2020年5期)2020-03-29 02:10:36
    同化衛(wèi)星高度計(jì)觀測對CAS-ESM-C上層海洋溫度模擬的改進(jìn)
    2016年7月四川持續(xù)性強(qiáng)降水的中尺度濾波分析
    黃淮地區(qū)一次暖區(qū)大暴雨的中尺度特征分析
    合成孔徑雷達(dá)高度計(jì)與傳統(tǒng)高度計(jì)精度比對分析與機(jī)載試驗(yàn)驗(yàn)證
    老鸭窝网址在线观看| 三级毛片av免费| 色噜噜av男人的天堂激情| 精品午夜福利视频在线观看一区| 长腿黑丝高跟| 三级男女做爰猛烈吃奶摸视频| 18禁黄网站禁片午夜丰满| 国产精品久久久av美女十八| 国产单亲对白刺激| 免费一级毛片在线播放高清视频| 国产高清视频在线观看网站| 人妻久久中文字幕网| 亚洲乱码一区二区免费版| 中文字幕久久专区| 亚洲五月婷婷丁香| 亚洲国产精品999在线| 欧美久久黑人一区二区| 99精品久久久久人妻精品| 男女床上黄色一级片免费看| 久久久国产成人免费| 欧美绝顶高潮抽搐喷水| 别揉我奶头~嗯~啊~动态视频| 久久精品影院6| 色综合欧美亚洲国产小说| 88av欧美| 成人国产综合亚洲| 久久性视频一级片| 成在线人永久免费视频| 亚洲av成人不卡在线观看播放网| 久久久久久久久免费视频了| 男男h啪啪无遮挡| 制服丝袜大香蕉在线| 1024香蕉在线观看| 一级毛片精品| videosex国产| 婷婷精品国产亚洲av| 亚洲第一电影网av| 老汉色av国产亚洲站长工具| 国产欧美日韩一区二区精品| 三级国产精品欧美在线观看 | 国产精品 欧美亚洲| 一边摸一边抽搐一进一小说| 伊人久久大香线蕉亚洲五| 国产成人啪精品午夜网站| 女生性感内裤真人,穿戴方法视频| 亚洲成人久久性| 欧美极品一区二区三区四区| 99久久无色码亚洲精品果冻| 一区二区三区激情视频| 最新在线观看一区二区三区| 午夜福利免费观看在线| 身体一侧抽搐| 国产av一区在线观看免费| 日日干狠狠操夜夜爽| 亚洲国产日韩欧美精品在线观看 | 免费电影在线观看免费观看| 91大片在线观看| 亚洲熟妇中文字幕五十中出| 成人国产综合亚洲| 国产99久久九九免费精品| 夜夜爽天天搞| 中文字幕人妻丝袜一区二区| 午夜视频精品福利| 精品第一国产精品| 两性午夜刺激爽爽歪歪视频在线观看 | 老司机福利观看| 精品久久久久久久毛片微露脸| 国产午夜精品论理片| 国产97色在线日韩免费| 精品一区二区三区av网在线观看| 亚洲人成电影免费在线| 国产三级在线视频| 国产精品99久久99久久久不卡| 三级毛片av免费| 欧美日韩亚洲综合一区二区三区_| 久久天堂一区二区三区四区| 99国产极品粉嫩在线观看| 在线十欧美十亚洲十日本专区| 国产黄a三级三级三级人| 欧美三级亚洲精品| 校园春色视频在线观看| 2021天堂中文幕一二区在线观| 成人国语在线视频| 国产亚洲精品久久久久5区| 麻豆成人av在线观看| 精品久久久久久久久久免费视频| 麻豆成人午夜福利视频| 亚洲自拍偷在线| 久热爱精品视频在线9| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产精品成人综合色| 久久中文看片网| 91字幕亚洲| 亚洲一区二区三区不卡视频| 老司机福利观看| 日韩高清综合在线| 国产成年人精品一区二区| 一二三四在线观看免费中文在| 国产激情偷乱视频一区二区| 国产精品免费视频内射| 精品国产乱码久久久久久男人| 欧美日本视频| 午夜精品在线福利| 国产蜜桃级精品一区二区三区| 男男h啪啪无遮挡| 日韩欧美在线二视频| 亚洲av熟女| 国产蜜桃级精品一区二区三区| 老熟妇乱子伦视频在线观看| 成人18禁在线播放| 久久这里只有精品中国| 黑人操中国人逼视频| 亚洲精品国产一区二区精华液| 色综合站精品国产| 伊人久久大香线蕉亚洲五| 欧洲精品卡2卡3卡4卡5卡区| 男插女下体视频免费在线播放| 国产欧美日韩精品亚洲av| av中文乱码字幕在线| 黑人巨大精品欧美一区二区mp4| 国产在线精品亚洲第一网站| 久久久久国产一级毛片高清牌| av超薄肉色丝袜交足视频| 国产精品一区二区三区四区久久| 亚洲成人久久性| 三级国产精品欧美在线观看 | 精品国产超薄肉色丝袜足j| 操出白浆在线播放| 男人舔奶头视频| 特级一级黄色大片| 淫妇啪啪啪对白视频| 90打野战视频偷拍视频| 少妇人妻一区二区三区视频| 一级毛片精品| bbb黄色大片| 高清毛片免费观看视频网站| 极品教师在线免费播放| 搡老岳熟女国产| 婷婷精品国产亚洲av| 一级毛片女人18水好多| 久久精品国产99精品国产亚洲性色| 又大又爽又粗| 成在线人永久免费视频| 欧美精品亚洲一区二区| 一区二区三区国产精品乱码| 一边摸一边抽搐一进一小说| 日韩 欧美 亚洲 中文字幕| 黄色a级毛片大全视频| 一进一出抽搐动态| 黄色视频不卡| 欧美乱色亚洲激情| 狂野欧美白嫩少妇大欣赏| 亚洲av成人av| 国产精品久久久久久久电影 | 色综合亚洲欧美另类图片| 中文字幕av在线有码专区| 色av中文字幕| 精品一区二区三区四区五区乱码| 一级作爱视频免费观看| 91麻豆精品激情在线观看国产| 国产爱豆传媒在线观看 | 国产成人啪精品午夜网站| 欧美中文日本在线观看视频| 亚洲人与动物交配视频| 中国美女看黄片| 日韩免费av在线播放| 亚洲欧美日韩东京热| 色综合欧美亚洲国产小说| 非洲黑人性xxxx精品又粗又长| 久久中文字幕人妻熟女| 好看av亚洲va欧美ⅴa在| 成年人黄色毛片网站| 听说在线观看完整版免费高清| 日本五十路高清| 91字幕亚洲| 亚洲精品一卡2卡三卡4卡5卡| 日本熟妇午夜| 天堂影院成人在线观看| 可以在线观看的亚洲视频| 免费搜索国产男女视频| a级毛片a级免费在线| 国产av一区二区精品久久| 身体一侧抽搐| 美女免费视频网站| 色哟哟哟哟哟哟| 免费一级毛片在线播放高清视频| 欧美极品一区二区三区四区| 夜夜夜夜夜久久久久| 午夜福利免费观看在线| 日韩成人在线观看一区二区三区| 天天躁夜夜躁狠狠躁躁| 18禁观看日本| 搞女人的毛片| 精品午夜福利视频在线观看一区| 91国产中文字幕| 老鸭窝网址在线观看| 男插女下体视频免费在线播放| 国模一区二区三区四区视频 | 午夜日韩欧美国产| 少妇裸体淫交视频免费看高清 | 一本大道久久a久久精品| 成人av在线播放网站| 色播亚洲综合网| 中国美女看黄片| 午夜福利视频1000在线观看| 亚洲黑人精品在线| 亚洲五月天丁香| 一本久久中文字幕| 最近视频中文字幕2019在线8| 男人舔女人下体高潮全视频| 日本 av在线| 国产私拍福利视频在线观看| 国产三级中文精品| 精品熟女少妇八av免费久了| 日韩欧美免费精品| 欧美成人性av电影在线观看| 亚洲精品粉嫩美女一区| 久久久久久久久久黄片| 悠悠久久av| 精品少妇一区二区三区视频日本电影| 一区二区三区国产精品乱码| 午夜精品一区二区三区免费看| svipshipincom国产片| 精品国产乱码久久久久久男人| 亚洲18禁久久av| 男女做爰动态图高潮gif福利片| 日本五十路高清| 亚洲精品在线美女| 男人舔奶头视频| 五月伊人婷婷丁香| 又粗又爽又猛毛片免费看| 成人永久免费在线观看视频| 又黄又爽又免费观看的视频| 久久精品亚洲精品国产色婷小说| 波多野结衣巨乳人妻| www.www免费av| 国产午夜精品久久久久久| 制服丝袜大香蕉在线| 久久久久久久午夜电影| 搡老岳熟女国产| 老熟妇乱子伦视频在线观看| 欧美日韩国产亚洲二区| 精品福利观看| 岛国在线观看网站| 91老司机精品| 人妻丰满熟妇av一区二区三区| 国产aⅴ精品一区二区三区波| 在线观看舔阴道视频| 国产高清videossex| 日本 欧美在线| 女生性感内裤真人,穿戴方法视频| 精品电影一区二区在线| 一区二区三区高清视频在线| 精品久久蜜臀av无| 中出人妻视频一区二区| 久久精品国产99精品国产亚洲性色| 欧美在线一区亚洲| 我要搜黄色片| 男女之事视频高清在线观看| 精品久久久久久,| 精品电影一区二区在线| 女同久久另类99精品国产91| 国产精品99久久99久久久不卡| 午夜视频精品福利| av天堂在线播放| 美女大奶头视频| 亚洲,欧美精品.| 亚洲熟妇熟女久久| 一夜夜www| or卡值多少钱| 国产蜜桃级精品一区二区三区| 亚洲成人精品中文字幕电影| 悠悠久久av| 亚洲 国产 在线| 日韩免费av在线播放| 一级黄色大片毛片| 亚洲电影在线观看av| 亚洲一区二区三区色噜噜| 亚洲成人久久性| 婷婷精品国产亚洲av在线| 日韩成人在线观看一区二区三区| 国内揄拍国产精品人妻在线| 国产在线观看jvid| 精品久久久久久,| 波多野结衣高清无吗| 最好的美女福利视频网| 国产乱人伦免费视频| 小说图片视频综合网站| 日本黄大片高清| 国产爱豆传媒在线观看 | 琪琪午夜伦伦电影理论片6080| 亚洲av第一区精品v没综合| 国产精品久久视频播放| 亚洲成a人片在线一区二区| 久久这里只有精品19| 国产97色在线日韩免费| 嫁个100分男人电影在线观看| 一本一本综合久久| 又爽又黄无遮挡网站| 国内精品久久久久精免费| 真人一进一出gif抽搐免费| 黄片大片在线免费观看| 欧美黑人巨大hd| 久久久久久久精品吃奶| 久久久久久久久久黄片| 999精品在线视频| 亚洲精品国产一区二区精华液| 首页视频小说图片口味搜索| 亚洲欧美精品综合久久99| 国产aⅴ精品一区二区三区波| x7x7x7水蜜桃| 久9热在线精品视频| 麻豆国产av国片精品| 一级作爱视频免费观看| 免费看十八禁软件| 久久久国产成人免费| 99久久综合精品五月天人人| 日韩欧美国产在线观看| 精品国产乱子伦一区二区三区| 丁香六月欧美| 欧美成人午夜精品| 国产精品亚洲一级av第二区| 我的老师免费观看完整版| 97人妻精品一区二区三区麻豆| 国产精品一及| 国产黄色小视频在线观看| 免费搜索国产男女视频| 国产av不卡久久| 国产欧美日韩精品亚洲av| 亚洲精华国产精华精| 亚洲人与动物交配视频| av片东京热男人的天堂| www.精华液| 国产激情欧美一区二区| 老司机深夜福利视频在线观看| 久热爱精品视频在线9| 波多野结衣高清无吗| 在线观看一区二区三区| 两个人免费观看高清视频| 婷婷六月久久综合丁香| 亚洲人成77777在线视频| 国产亚洲精品久久久久5区| 母亲3免费完整高清在线观看| 两性夫妻黄色片| aaaaa片日本免费| 精品日产1卡2卡| 亚洲人成网站高清观看| 亚洲av中文字字幕乱码综合| 午夜福利在线观看吧| 黄片大片在线免费观看| 国产激情久久老熟女| 一级毛片女人18水好多| 久久精品亚洲精品国产色婷小说| av福利片在线观看| 18禁观看日本| av欧美777| 欧美av亚洲av综合av国产av| 12—13女人毛片做爰片一| 美女黄网站色视频| 国产成人精品无人区| 欧美黑人巨大hd| av福利片在线| 人妻久久中文字幕网| 国产av麻豆久久久久久久| 成人国产综合亚洲| 亚洲自偷自拍图片 自拍| 精品国产超薄肉色丝袜足j| 无遮挡黄片免费观看| 欧美性猛交黑人性爽| 亚洲精品久久国产高清桃花| 免费在线观看黄色视频的| 91麻豆av在线| 亚洲欧美日韩高清在线视频| 亚洲av五月六月丁香网| 村上凉子中文字幕在线| 变态另类丝袜制服| 亚洲av第一区精品v没综合| 99在线视频只有这里精品首页| www.熟女人妻精品国产| 欧美绝顶高潮抽搐喷水| 国产男靠女视频免费网站| 18禁国产床啪视频网站| 一进一出好大好爽视频| 亚洲av日韩精品久久久久久密| 亚洲av第一区精品v没综合| 18美女黄网站色大片免费观看| 最新在线观看一区二区三区| 99精品在免费线老司机午夜| 精品不卡国产一区二区三区| 亚洲全国av大片| 熟女电影av网| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美在线一区二区| 丰满人妻一区二区三区视频av | 国产精品一区二区精品视频观看| 亚洲全国av大片| 黄频高清免费视频| 免费观看精品视频网站| 国产精品久久久久久亚洲av鲁大| 香蕉丝袜av| 国产精品 欧美亚洲| 国产精品综合久久久久久久免费| 久久人妻福利社区极品人妻图片| 美女 人体艺术 gogo| 欧美不卡视频在线免费观看 | 中亚洲国语对白在线视频| 禁无遮挡网站| 亚洲国产精品久久男人天堂| 免费一级毛片在线播放高清视频| 亚洲精品在线观看二区| 在线视频色国产色| 欧美性猛交黑人性爽| 亚洲成av人片免费观看| 好男人在线观看高清免费视频| 日韩欧美在线二视频| 国产精品一区二区三区四区免费观看 | 亚洲av成人不卡在线观看播放网| 两个人免费观看高清视频| 免费高清视频大片| 亚洲aⅴ乱码一区二区在线播放 | 亚洲av熟女| 精品无人区乱码1区二区| 国产人伦9x9x在线观看| 老司机深夜福利视频在线观看| 日韩欧美 国产精品| 亚洲狠狠婷婷综合久久图片| 欧美中文日本在线观看视频| 欧美日韩精品网址| 久久99热这里只有精品18| а√天堂www在线а√下载| 久久久精品国产亚洲av高清涩受| 久久热在线av| 亚洲av五月六月丁香网| 国产野战对白在线观看| 麻豆国产97在线/欧美 | 午夜福利免费观看在线| 欧美成狂野欧美在线观看| 精品午夜福利视频在线观看一区| 91在线观看av| 欧美日本亚洲视频在线播放| 亚洲av电影在线进入| 国产成人aa在线观看| 人妻久久中文字幕网| 真人做人爱边吃奶动态| 日韩国内少妇激情av| 久久久久亚洲av毛片大全| 99在线人妻在线中文字幕| 久久午夜亚洲精品久久| 亚洲av成人精品一区久久| 国产黄色小视频在线观看| 日韩三级视频一区二区三区| 久久久久九九精品影院| 99久久国产精品久久久| 日韩欧美在线二视频| 婷婷亚洲欧美| 宅男免费午夜| 日本 欧美在线| 欧美国产日韩亚洲一区| 一边摸一边做爽爽视频免费| 草草在线视频免费看| 精品久久久久久久末码| 欧美极品一区二区三区四区| 少妇熟女aⅴ在线视频| 亚洲欧美日韩无卡精品| 国产高清有码在线观看视频 | 十八禁人妻一区二区| 桃色一区二区三区在线观看| 亚洲国产欧美网| av福利片在线观看| 国产免费av片在线观看野外av| 99久久无色码亚洲精品果冻| x7x7x7水蜜桃| 精品国内亚洲2022精品成人| 国产精品一区二区精品视频观看| 巨乳人妻的诱惑在线观看| 欧美成人性av电影在线观看| 成年免费大片在线观看| 久久人妻av系列| 激情在线观看视频在线高清| 美女大奶头视频| 日韩欧美免费精品| 国产亚洲欧美98| 香蕉久久夜色| 日韩精品青青久久久久久| 久久久久久国产a免费观看| 亚洲色图 男人天堂 中文字幕| 久久午夜综合久久蜜桃| av福利片在线观看| 一区福利在线观看| 男女之事视频高清在线观看| 成人国语在线视频| 成人18禁在线播放| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久久久亚洲av鲁大| 欧美一区二区国产精品久久精品 | 色精品久久人妻99蜜桃| 亚洲中文日韩欧美视频| 久久亚洲精品不卡| 国产成人aa在线观看| 麻豆一二三区av精品| 伊人久久大香线蕉亚洲五| or卡值多少钱| 美女 人体艺术 gogo| 真人做人爱边吃奶动态| 午夜福利在线在线| 两个人视频免费观看高清| 国产激情久久老熟女| 欧美日韩精品网址| 亚洲美女黄片视频| 中文字幕熟女人妻在线| 国产成人aa在线观看| 日韩欧美在线乱码| 色噜噜av男人的天堂激情| a在线观看视频网站| 国产激情欧美一区二区| 男人舔奶头视频| 国产一区二区三区在线臀色熟女| 热99re8久久精品国产| 亚洲自偷自拍图片 自拍| 婷婷六月久久综合丁香| 欧美人与性动交α欧美精品济南到| 欧美高清成人免费视频www| 日韩欧美国产一区二区入口| 久久亚洲真实| 欧美另类亚洲清纯唯美| 日韩国内少妇激情av| 欧美av亚洲av综合av国产av| 男女之事视频高清在线观看| 久久人妻av系列| 大型黄色视频在线免费观看| 黄频高清免费视频| 国产av一区在线观看免费| 国产伦在线观看视频一区| 女人高潮潮喷娇喘18禁视频| 一个人观看的视频www高清免费观看 | 国产蜜桃级精品一区二区三区| 麻豆成人午夜福利视频| 特级一级黄色大片| 国产私拍福利视频在线观看| 香蕉国产在线看| svipshipincom国产片| 一二三四在线观看免费中文在| 欧美黑人精品巨大| 免费看美女性在线毛片视频| 亚洲精品国产一区二区精华液| 亚洲av片天天在线观看| 精品久久久久久,| 黄频高清免费视频| 正在播放国产对白刺激| 黄色 视频免费看| 精品国产乱码久久久久久男人| 色av中文字幕| 午夜激情av网站| 亚洲精品美女久久久久99蜜臀| 成人国语在线视频| 精品不卡国产一区二区三区| 国产精品自产拍在线观看55亚洲| 校园春色视频在线观看| 午夜两性在线视频| 长腿黑丝高跟| 99久久精品国产亚洲精品| 白带黄色成豆腐渣| 岛国视频午夜一区免费看| 亚洲人与动物交配视频| 两性午夜刺激爽爽歪歪视频在线观看 | 给我免费播放毛片高清在线观看| 最新美女视频免费是黄的| 97碰自拍视频| 欧美黑人巨大hd| 天堂av国产一区二区熟女人妻 | 国产男靠女视频免费网站| 一级作爱视频免费观看| 香蕉国产在线看| 日韩欧美在线二视频| 国产精品香港三级国产av潘金莲| 欧美中文日本在线观看视频| 亚洲av成人不卡在线观看播放网| 很黄的视频免费| 国产精品影院久久| 国产真实乱freesex| 亚洲中文字幕一区二区三区有码在线看 | 老熟妇乱子伦视频在线观看| 一个人观看的视频www高清免费观看 | 国内毛片毛片毛片毛片毛片| 国产午夜精品论理片| 最近视频中文字幕2019在线8| 欧美高清成人免费视频www| 亚洲中文av在线| 色综合婷婷激情| 国产精品精品国产色婷婷| 黑人操中国人逼视频| 久久中文字幕一级| 免费在线观看视频国产中文字幕亚洲| 精品免费久久久久久久清纯| 老司机深夜福利视频在线观看| 国产成人精品无人区| 激情在线观看视频在线高清| 欧美在线一区亚洲| www日本在线高清视频| 制服丝袜大香蕉在线| 亚洲国产欧美网| 99国产综合亚洲精品| 少妇熟女aⅴ在线视频| 日韩av在线大香蕉| 国产精品免费视频内射| 国产精华一区二区三区| 999精品在线视频| 国产成人aa在线观看| 国产三级在线视频| 少妇的丰满在线观看| 欧美在线黄色| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 免费看日本二区| 搡老岳熟女国产|