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

    日本海西南海域現(xiàn)場(chǎng)觀測(cè)和衛(wèi)星高度計(jì)獲取的海面高度距平的比較研究

    2016-07-06 10:30:17徐永生尹寶樹中國(guó)科學(xué)院海洋研究所山東青島266071中國(guó)科學(xué)院大學(xué)北京100049中國(guó)科學(xué)院海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室山東青島266071青島海洋科學(xué)與技術(shù)國(guó)家實(shí)驗(yàn)室海洋動(dòng)力過(guò)程與氣候功能實(shí)驗(yàn)室山東青島266237
    海洋科學(xué) 2016年2期
    關(guān)鍵詞:斜壓高度計(jì)距平

    葛 磊, 徐永生, 尹寶樹(1. 中國(guó)科學(xué)院海洋研究所, 山東 青島 266071; 2. 中國(guó)科學(xué)院大學(xué), 北京 100049; 3. 中國(guó)科學(xué)院 海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266071; 4. 青島海洋科學(xué)與技術(shù)國(guó)家實(shí)驗(yàn)室海洋動(dòng)力過(guò)程與氣候功能實(shí)驗(yàn)室,山東 青島 266237)

    ?

    日本海西南海域現(xiàn)場(chǎng)觀測(cè)和衛(wèi)星高度計(jì)獲取的海面高度距平的比較研究

    葛 磊1, 2, 3, 徐永生1, 3, 4, 尹寶樹1, 3, 4
    (1. 中國(guó)科學(xué)院海洋研究所, 山東 青島 266071; 2. 中國(guó)科學(xué)院大學(xué), 北京 100049; 3. 中國(guó)科學(xué)院 海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266071; 4. 青島海洋科學(xué)與技術(shù)國(guó)家實(shí)驗(yàn)室海洋動(dòng)力過(guò)程與氣候功能實(shí)驗(yàn)室,山東 青島 266237)

    壓力傳感逆式回聲儀(pressure-sensor-equipped inverted echo sounders, PIES)可以用來(lái)測(cè)量海底壓力和聲波從海底到海面的傳播時(shí)間。海底壓力和聲波傳播時(shí)間分別被用來(lái)估計(jì)水體質(zhì)量變化(正壓)和比容變化(斜壓)對(duì)海面高度距平的貢獻(xiàn)。對(duì)由PIES在日本海西南海域現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)得到的海面高度距平(PIES SLA)與衛(wèi)星高度計(jì)海面高度距平(Sat SLA)進(jìn)行了比較研究。利用相關(guān)分析法, 對(duì)PIES SLA和沿軌T/P衛(wèi)星、沿軌ERS-2衛(wèi)星測(cè)得的海面高度距平(TP SLA、ERS-2 SLA)進(jìn)行了比較; 對(duì)PIES SLA 和AVISO網(wǎng)格化海面高度距平進(jìn)行了比較, 估計(jì)可能的誤差來(lái)源, 并分析PIES SLA正壓部分和斜壓部分對(duì)SLA的貢獻(xiàn)。比較發(fā)現(xiàn), PIES SLA和Sat SLA的相關(guān)系數(shù)較高, 且均方根誤差較小, 并且對(duì)特定區(qū)域和特定站點(diǎn)產(chǎn)生誤差可能的原因進(jìn)行了進(jìn)一步的探討。通過(guò)研究, 有以下結(jié)論: (1)相對(duì)于灣流和黑潮地區(qū), 這一區(qū)域正壓部分對(duì)海面高度的貢獻(xiàn)相對(duì)較大; (2)如果再考慮斜壓變化對(duì)海面高度的貢獻(xiàn), PIES SLA和Sat SLA相關(guān)系數(shù)會(huì)有所提升; (3)在高能區(qū)PIES SLA和Sat SLA相關(guān)系數(shù)較高, 符合得相對(duì)比較好。總的來(lái)說(shuō), 在日本海地區(qū), PIES SLA和Sat SLA相關(guān)系數(shù)較高, 具有較高的一致性, 能為我國(guó)海洋二號(hào)(HY-2)等衛(wèi)星高度計(jì)的校驗(yàn)提供一種可靠的方式。該研究對(duì)于PIES的研發(fā)和設(shè)計(jì)以及對(duì)于PIES的布放位置的選擇都有一定的借鑒意義。

    海面高度距平; 壓力傳感逆式回聲儀(PIES); 衛(wèi)星高度計(jì); 聲波傳播時(shí)間; 海水底部壓力

    (全球變化與海氣相互作用)及中國(guó)科學(xué)院知識(shí)創(chuàng)新項(xiàng)目(Y22114101Q) [Foundation: 100 Talents Program of the Chinese Academy of Sciences, No.Y32109101L; National Natural Science Foundation of China, No.41376028; Key Topics in Innovation Engineering of the Chinese Academy of Sciences, No. KZCX2-EW-209; State Oceanic Administration, (Global change and the air sea interaction project), Knowledge Innovation Project of the Chinese Academy of Sciences, No.Y22114101Q]

    回聲探測(cè)儀被廣泛地用于測(cè)量海洋深度等海洋要素, 但是其容易受極端海況、惡劣天氣、繩索傾斜等影響, 從而導(dǎo)致較大的測(cè)量誤差。而逆式回聲探測(cè)儀(inverted echo sounder, IES)與一般的聲波探測(cè)儀正好相反, 它將聲波發(fā)射器置于海底, 向上發(fā)射聲波, 經(jīng)海面反射后回到海底, 從而記錄聲波由海底到海面的往返傳播時(shí)間τ。如果再加上一個(gè)壓力傳感器, 用來(lái)記錄海洋底部壓力Pbot, 便構(gòu)成了壓力傳感逆式回聲儀(pressure-sensor-equipped inverted echo sounders, PIES)。由于具有觀測(cè)周期長(zhǎng)、安全性高、適合陣列式大范圍觀測(cè)的優(yōu)點(diǎn), PIES已發(fā)展成為近代海洋觀測(cè)的一個(gè)重要手段, 在全球海洋環(huán)境監(jiān)測(cè)、海洋環(huán)流動(dòng)力學(xué)研究以及軍事海洋學(xué)等多個(gè)應(yīng)用領(lǐng)域, 全世界大洋多個(gè)海區(qū)得到了成功應(yīng)用。

    自從1969年Rossby[1]發(fā)展逆式回聲儀(inverted echo sounders, IES)以來(lái), 人們一直在發(fā)展對(duì)聲波τ的處理技術(shù), 以便從聲波τ中獲取更多的海洋信息, Meinen 和Watts[2]、Sun和 Watts[3]、Watts等[4]得到了重力經(jīng)驗(yàn)?zāi)7椒?gravest empirical mode, GEM), 能夠利用聲波τ, 來(lái)獲得海水溫度和比容異常的剖面或者廓線。

    近年來(lái), 衛(wèi)星高度計(jì)發(fā)展成為了海洋觀測(cè)衛(wèi)星的基本觀測(cè)儀器之一, 其能夠?qū)θ蚝Q筮M(jìn)行大面積、全天候持續(xù)觀測(cè), 對(duì)全球海洋進(jìn)行精確的測(cè)量。隨著衛(wèi)星測(cè)高技術(shù)的發(fā)展, 衛(wèi)星測(cè)量海面高度(SSH)得到了越來(lái)越廣泛的應(yīng)用, 在我們認(rèn)識(shí)全球環(huán)流分布、中尺度渦、海面動(dòng)力過(guò)程、海洋風(fēng)浪的反演等研究中起著重要的作用。

    Ponte[5]、Guinehut等[6]、Jayne等[7]認(rèn)為海面高度變化由以下兩個(gè)部分構(gòu)成: 由于海水質(zhì)量變化引起的正壓部分(barotropic components)和由于海水溫度、密度等變化引起的斜壓部分(baroclinic components), 有些文獻(xiàn)上也分別稱為質(zhì)量負(fù)載變化(massloading)和比容變化(steric)。Baker-Yeboah等[8]認(rèn)為水柱中質(zhì)量變化(正壓部分)對(duì)SSH的貢獻(xiàn)和純粹由比容變化(斜壓部分)對(duì)SSH的貢獻(xiàn)均可以由PIES記錄數(shù)據(jù)來(lái)得到。也就是說(shuō)PIES能夠通過(guò)同時(shí)測(cè)量聲波τ和底部壓力Pbot, 來(lái)獲得對(duì)應(yīng)SSH的斜壓部分和正壓部分。衛(wèi)星高度計(jì)測(cè)量的SSH, 當(dāng)然也由正壓和斜壓兩部分組成。Fukumori等[9]認(rèn)為跟密度變化相關(guān)的斜壓部分被認(rèn)為是主導(dǎo)SSH低頻變化的主要因素, 而正壓部分則對(duì)高頻帶部分貢獻(xiàn)較大。

    衛(wèi)星觀測(cè)和現(xiàn)場(chǎng)觀測(cè)得到的SSH之間的定量比較少, 因?yàn)樵陂_(kāi)放的海洋中, 難以同時(shí)獲得兩種測(cè)量方式得到SSH的正壓部分和斜壓部分。較早的一些研究中, Gilson等[10]、McCarthy等[11]對(duì)XBT數(shù)據(jù)得到的SSH距平與沿著部分衛(wèi)星測(cè)量的SSH距平進(jìn)行了比較, 以達(dá)到評(píng)估衛(wèi)星測(cè)高的目的。

    這次布放在日本海(Japan/East Sea, JES)西南部區(qū)域的持續(xù)觀測(cè)達(dá)2 a之久24個(gè)PIES, 為我們比較逆式回聲儀測(cè)得的海面高度距平(以下簡(jiǎn)稱為PIES SLA)和衛(wèi)星高度計(jì)海面高度距平(以下簡(jiǎn)稱為Sat SLA)提供了一個(gè)很好的機(jī)會(huì)。本文將利用相關(guān)分析法, 利用衛(wèi)星高度數(shù)據(jù)產(chǎn)品包括AVISO的T/P沿軌(along-track)海面高度距平數(shù)據(jù)產(chǎn)品、ERS-2沿軌(along-track)海面高度距平產(chǎn)品以及網(wǎng)格化的海面高度距平產(chǎn)品, 對(duì)PIES SLA和沿軌T/P衛(wèi)星、沿軌ERS-2測(cè)得的海面高度距平(以下簡(jiǎn)稱為TP SLA、ERS-2 SLA)以及AVISO網(wǎng)格化海面高度距平進(jìn)行比較, 估計(jì)PIES SLA和Sat SLA之間可能的誤差來(lái)源, 并分析PIES SLA正壓部分和斜壓部分對(duì)比較結(jié)果的影響。研究發(fā)現(xiàn)在日本海地區(qū), PIES SLA和衛(wèi)星高度計(jì)SLA相關(guān)系數(shù)較高, 具有較高的一致性, 為我國(guó)海洋二號(hào)(HY-2)等衛(wèi)星高度計(jì)提供了一種比較可靠的校驗(yàn)方式。

    1 數(shù)據(jù)和方法

    在1999年6月到2001年7月期間, 在日本海西南海域的郁陵海盆布放了一個(gè)24個(gè)PIES的陣列, 布放的位置如圖1所示。

    圖1 24個(gè)PIES站點(diǎn)的位置示意圖Fig. 1 Positions of the 24 PIES stations

    對(duì)PIES測(cè)量的聲波τ和Pbot進(jìn)行處理, 獲得時(shí)間分辨率為1 h的數(shù)據(jù)。τ的分辨率是0.05 ms, PIES 的Paroscientific壓力傳感器測(cè)Pbot的精確度±0.01%,分辨率是10 Pa, 大概是0.1 cm。這些數(shù)據(jù)經(jīng)過(guò)剔除異常值、去除漂移(dedrift)、去除跳躍(dejump)、深度校正(leveling)等一系列的處理后[12-14], 并且去除各種表面波的影響, 再進(jìn)行120 h (5 d)的低通濾波處理得到新的時(shí)間分辨率為1 h, 并在每天的00時(shí)和12時(shí)(UTM)重新二次采樣。對(duì)PIES測(cè)得的τ和Pbot的數(shù)據(jù)處理的過(guò)程中, 之所以進(jìn)行深度校正是為了方便以后的科學(xué)分析和研究, 深度校正一般都要選擇一個(gè)參考深度, 這個(gè)參考深度根據(jù)觀測(cè)區(qū)域的不同而不同, 比如在灣流區(qū)域和大西洋的某些區(qū)域,這個(gè)深度能都達(dá)到4 500 m甚至5 000 m, 而Park等[15]等在黑潮延伸體區(qū)域, 這個(gè)深度選擇了4 000 m。因?yàn)楸敬斡^測(cè)深度大都不是太深, 并且500 m以下的斜壓部分就很弱了, 因此參考深度定為500 m。也就是說(shuō)τ 和Pbot要通過(guò)一定的轉(zhuǎn)換關(guān)系, 轉(zhuǎn)換成為τ500和P500。

    下面將討論如何由PIES測(cè)量的τ和Pbot數(shù)據(jù)得到SSH距平。PIES測(cè)量數(shù)據(jù)對(duì)SSH的貢獻(xiàn)可以分為兩個(gè)部分, 一部分是由于水體質(zhì)量負(fù)載變化導(dǎo)致的正壓變化, 這部分可以通過(guò)Pbot的變化得到, 一部分是由于密度變化導(dǎo)致的比容變化, 或者稱為斜壓部分, 這部分可以通過(guò)τ得到。

    我們將SSH記為η, SSH的正壓部分記為ηbt, SSH的斜壓部分記為ηbc, 而η, ηbt, ηbc三者的距平分別記為η′, η′bt, η′bc, 底部壓力記為Pbot, 將聲波傳播時(shí)間記為τ, Baker-Yeboah等[8]給出,

    這里, g代表重力加速度, δ是海水溫度(T)、鹽度(S)和壓強(qiáng)(p)的函數(shù),

    ηbc表示由于位勢(shì)高度異常導(dǎo)致的海面變化, 是由比容或者是跟密度相關(guān)的變化(斜壓)導(dǎo)致的海面變化。

    對(duì)于ηbc來(lái)說(shuō), 去除時(shí)間平均值, 便有η斜壓部分的距平

    和η正壓部分的距平。

    而聲波傳播的時(shí)間是

    這里,sρ代表海表層海水的密度。g代表重力加速度, c代表聲波傳播時(shí)間, c是海水溫度鹽度和壓強(qiáng)的函數(shù),

    而又由于公式(1)的關(guān)系, 所以說(shuō)要由τ確定位勢(shì)高度變化, 就要使用GEM方法, 建立起τ和位勢(shì)高度變化之間的聯(lián)系。

    前面提到, GEM方法首先由Meinen 和Watts[4]為了研究北大西洋海流的垂直結(jié)構(gòu)和輸運(yùn)發(fā)展而來(lái),并且成功地應(yīng)用到副極地鋒面(sub antarctic front)[6],南極繞極流和黑潮小彎曲的發(fā)展[15]的研究當(dāng)中。這些區(qū)域都有兩個(gè)重要的特征, 一是在107Pa處或者更深, 有強(qiáng)的斜壓變化, 導(dǎo)致斜壓信號(hào)比較強(qiáng)。二是不同的水文廓線能夠?qū)?yīng)產(chǎn)生不同的τ。而日本海西南部地區(qū)恰恰不具備這兩點(diǎn), 于是針對(duì)這一海域空間季節(jié)信號(hào)變化和淺溫躍層的特點(diǎn), Mitchell等[16]在GEM的基礎(chǔ)上, 提出了剩余GEM(residual GEM)。把GEM捕捉溫度T的變化和比容(斜壓)的變化分別從70%和64%提升到了89%和84%。

    τ經(jīng)過(guò)深度校正后得到τ500對(duì)SSH距平斜壓部分ΔH的貢獻(xiàn)見(jiàn)圖2。

    圖2 SSH距平斜壓部分ΔH與τ500對(duì)應(yīng)關(guān)系圖Fig. 2 Correspondence between ΔH (baroclinic components of the SSH) and τ500

    這次比較中用到了3種AVISO衛(wèi)星高度計(jì)產(chǎn)品:一個(gè)是T/P衛(wèi)星遲時(shí)沿軌海面高度距平產(chǎn)品, 一個(gè)是ERS-2衛(wèi)星遲時(shí)沿軌海面高度距平產(chǎn)品, 這兩個(gè)產(chǎn)品都是由AVISO制作和分發(fā)的(AVISO DT CorSSH and DT SLA Product Handbook, 2012)。另一個(gè)是“參考”(reference)遲時(shí)的網(wǎng)格化產(chǎn)品(reference maps of sea level anomalies, Ref-MSLA), 空間分辨率是(1/4)°× (1/4)°, 時(shí)間分辨率是7 d, 主要是融合了T/P和ERS-2,以及它們對(duì)應(yīng)后繼的Jason-1和Envisat衛(wèi)星數(shù)據(jù)。T/P衛(wèi)星軌道和ERS-2衛(wèi)星軌道與24個(gè)PIES站點(diǎn)位置示意圖如圖3所示。遲時(shí)的網(wǎng)格化衛(wèi)星數(shù)據(jù)的分辨率和24個(gè)PIES站點(diǎn)位置示意圖如圖4所示。

    圖3 T/P衛(wèi)星軌道和ERS-2衛(wèi)星軌道與24個(gè)PIES站點(diǎn)位置示意圖Fig. 3 Schematic of the T/P and ERS-2 orbits, and the positions of the 24 PIES stations

    2 海面高度的比較

    2.1 PIES SLA與TP SLA的比較

    從圖3中可以看出來(lái), P1-6站點(diǎn)正好在T/P軌道掃過(guò)的區(qū)域內(nèi)。T/P衛(wèi)星周期約為10 d, 相鄰資料的空間間隔(沿軌)是5.8 km, 為了創(chuàng)建能夠與PIES SLA可比的T/P衛(wèi)星時(shí)間序列, T/P衛(wèi)星上軌道上距離P1-6站點(diǎn)最近的點(diǎn)將被確定為可比的數(shù)據(jù)點(diǎn)。P1-6站點(diǎn)處總的PIES SLA和TP SLA的對(duì)比圖見(jiàn)圖5, 從圖5中可以看出總的PIES SLA, 能夠比較好地反映TP SLA的大致趨勢(shì), 吻合度比較高, 計(jì)算兩者的相關(guān)性, 能夠達(dá)到0.914, 計(jì)算二者的均方根誤差為8.7 cm。僅僅考慮PIES BT SLA與TP SLA相比較,從圖上看出二者吻合得也很好, 但是計(jì)算兩者的相關(guān)性,稍微有點(diǎn)降低達(dá)到0.887, 但是均方根誤差為7.0 cm。僅僅考慮斜壓部分PIES BC SLA的與TP SLA的比較, 從圖上看出二者吻合的稍差, 計(jì)算兩者的相關(guān)性, 降低到了0.87, 均方根誤差也稍微增大為7.6 cm。

    圖4 PIES站點(diǎn)位置和網(wǎng)格化衛(wèi)星數(shù)據(jù)空間分辨率的相對(duì)位置圖Fig. 4 Positions of the PIES stations and the grid of the satellite SSH data

    Teague等[17]將PIES數(shù)據(jù)得到的SSH距平與T/P高度計(jì)得到的海面距平做了比較, 結(jié)果顯示兩者符合的比較好, 但是均方根誤差達(dá)到了12 cm, 與他的結(jié)果相比, 無(wú)論是從相關(guān)性來(lái)說(shuō)還是從均方根誤差方面來(lái)講, PIES SLA與T/P SLA的比較結(jié)果都有顯著的改進(jìn)。從圖5還可以看出, 相對(duì)于灣流和黑潮地區(qū)(Book[18]), 在這一海域正壓部分對(duì)海面高度的貢獻(xiàn)相對(duì)較大, 如果再考慮斜壓變化對(duì)海面高度的貢獻(xiàn), PIES SLA和TP SLA相關(guān)系數(shù)將會(huì)進(jìn)一步提升。

    2.2 PIES SLA與ERS-2 SLA的比較

    從圖3中可以看出來(lái), P1-4, P2-2, P2-5, P3-1, P4-1, P4-2站點(diǎn)正好在ERS-2軌道掃過(guò)的區(qū)域內(nèi)。ERS-2周期大約為35 d。為了創(chuàng)建能夠與PIES-SSH可比的ERS-2衛(wèi)星時(shí)間序列, ERS-2衛(wèi)星上軌道上距離上述幾個(gè)站點(diǎn)最近的點(diǎn)將被確定為可比的數(shù)據(jù)點(diǎn)。以P2-5站點(diǎn)為例分析, P2-5站點(diǎn)處ERS-2衛(wèi)星SLA和PIES得到的總的SLA地對(duì)比圖見(jiàn)圖6, 從圖6中可以看出總的PIES TOTAL SLA能夠比較好地反映ERS-2 SLA的大致趨勢(shì), 吻合度比較高, 計(jì)算兩者的相關(guān)性, 能夠達(dá)到0.881, 計(jì)算二者的均方根誤差為7.62 cm。僅僅考慮PIES BT SLA與ERS-2 SLA的相比較, 從圖上看出二者吻合程度也很好, 但是計(jì)算兩者的相關(guān)性, 稍微有點(diǎn)降低達(dá)到0.851, 但是均方根誤差較小為6.8 cm。僅僅考慮PIES BC SLA 與ERS-2 SLA的比較, 從圖上看出二者吻合的稍差,計(jì)算兩者的相關(guān)性, 降低到了0.82, 均方根誤差也稍微增大為7.5 cm。從圖6看出, 相對(duì)于灣流和黑潮地區(qū)(Book[18]), 在這一海域正壓部分對(duì)海面高度的貢獻(xiàn)相對(duì)較大, 如果再考慮斜壓變化對(duì)海面高度的貢獻(xiàn), PIES SLA和ERS-2 SLA相關(guān)系數(shù)將會(huì)進(jìn)一步提升。

    圖5 P1-6站點(diǎn)處T/P衛(wèi)星SLA和PIES得到的SLA地對(duì)比圖Fig. 5 Comparison of the TP SLA and PIES SLA at the P1-6 station

    2.3 PIES SLA與Ref-MSLA的比較

    從圖4可以看出, 大部分PIES站點(diǎn)都不在Ref-MSLA網(wǎng)格上面, 這就需要將Ref-MSLA網(wǎng)格數(shù)據(jù)插值到PIES站點(diǎn)上面去。因?yàn)镽ef-MSLA時(shí)間分辨率是7 d, 為了得到能夠比較的時(shí)間序列, 需要對(duì)PIES得到的SLA數(shù)據(jù)進(jìn)行再采樣。對(duì)于所有的24個(gè)站點(diǎn)計(jì)算相關(guān)性(因P3-2處的可比性不強(qiáng), 故未列出), 整體的相關(guān)系數(shù)從0.75到0.94, 其中相關(guān)系數(shù)低于0.8只有近岸幾個(gè)站點(diǎn), 說(shuō)明大部分站點(diǎn)PIES-SLA和Ref-MSLA吻合的還是很好的。全部站點(diǎn)的相關(guān)系數(shù)見(jiàn)圖7。

    從圖7中可以看出, PIES-SLA和Ref-MSLA相關(guān)系數(shù)低于0.8的站點(diǎn)只有P1-1, P2-2, P3-1和P3-3,并且這4個(gè)站點(diǎn)大都處于近岸地區(qū)。相反, 越是處于深海地區(qū)的其他站點(diǎn), PIES-SLA和Ref-MSLA的相關(guān)系數(shù)反而比較高。在近岸地區(qū), PIES-SLA和Ref-MSLA相關(guān)系數(shù)相對(duì)較低, 究其原因可能有兩個(gè): 一個(gè)是日本海的邊緣海特征, 導(dǎo)致近岸地區(qū)相關(guān)海洋活動(dòng)如浪、流等較為劇烈, 影響了衛(wèi)星高度計(jì)的精確度; 二是我們發(fā)現(xiàn)有些近岸地區(qū)站點(diǎn)P3-1正好位于陸架坡折處, 而陸架坡折處(shelf break)的大地水準(zhǔn)面(geoid)比較陡峭, 衛(wèi)星的橫向不確定性導(dǎo)致衛(wèi)星測(cè)量SLA誤差較大。并且對(duì)于同一個(gè)站點(diǎn)P1-6, 就相關(guān)系數(shù)而言, 同單個(gè)T/P衛(wèi)星相比, Ref-MSLA的精確度確實(shí)有所提升。

    圖6 P2-5站點(diǎn)處ERS-2衛(wèi)星SLA和PIES得到的SLA地對(duì)比圖Fig. 6 Comparison of the ERS-2 SLA and PIES SLA at the P2-5 station

    3 結(jié)論與展望

    本文利用相關(guān)分析法, 對(duì)PIES SLA和沿軌T/P衛(wèi)星、沿軌ERS-2測(cè)得的海面高度距平以及AVISO網(wǎng)格化海面高度距平進(jìn)行了比較, 估計(jì)可能的誤差來(lái)源, 并分析了PIES SLA正壓部分和斜壓部分對(duì)比較結(jié)果的影響。比較發(fā)現(xiàn), PIES SLA和T/P SLA和的相關(guān)系數(shù)高達(dá)0.914, 均方根誤差(Root Mean Squared , RMS)為8.7 cm。PIES SLA和ERS-2 SLA的相關(guān)系數(shù)在0.75~0.90之間, 均方根誤差為7.6~ 10.7 cm。PIES SLA和AVISO網(wǎng)格化海面高度距平的相關(guān)系數(shù)在0.75到0.94之間。這些結(jié)果表明: (1)與灣流和黑潮地區(qū)斜壓部分對(duì)海面高度的貢獻(xiàn)相對(duì)較大不同, 這一區(qū)域正壓部分對(duì)海面高度的貢獻(xiàn)相對(duì)較大; (2)如果再考慮斜壓變化對(duì)海面高度的貢獻(xiàn), PIES SLA和衛(wèi)星高度計(jì)SLA相關(guān)系數(shù)會(huì)有所提升。(3)PIES SLA和衛(wèi)星高度計(jì)SLA的相關(guān)性跟海面變化信號(hào)的強(qiáng)弱有關(guān)系, 一般在高能區(qū), 或者主流區(qū), SLA和衛(wèi)星高度計(jì)SLA相關(guān)系數(shù)相對(duì)較高, 相關(guān)性較強(qiáng)。(4)在近岸地區(qū), 由于陸架坡折處大地水準(zhǔn)面(geoid)比較陡峭, 衛(wèi)星的橫向不確定性導(dǎo)致的測(cè)量誤差, 使得PIES SLA和Sat SLA相關(guān)系數(shù)相對(duì)較低??偟膩?lái)說(shuō), 在日本海地區(qū), PIES SLA和Sat SLA相關(guān)系數(shù)較高, 具有較高的一致性, 能為海洋二號(hào)衛(wèi)星(HY-2)高度計(jì)的校驗(yàn)提供一種可靠的方式。

    以往大都用驗(yàn)潮站、海上石油平臺(tái)和GPS浮標(biāo)來(lái)校驗(yàn)衛(wèi)星高度計(jì)數(shù)據(jù), 而驗(yàn)潮站數(shù)據(jù)大都處于近岸,因而校驗(yàn)的效果相對(duì)較差, 海上石油平臺(tái)需要高額的運(yùn)行維護(hù)費(fèi)用, 且存在諸多試驗(yàn)限制, 操作要求高, GPS浮標(biāo)容易受到海浪、風(fēng)的影響, 因此具有觀測(cè)周期長(zhǎng)、安全性高、適合陣列式大范圍觀測(cè)的優(yōu)點(diǎn)的PIES正好為我們提供了一個(gè)較好的校驗(yàn)機(jī)會(huì)。在高能區(qū)和主流區(qū), SLA和衛(wèi)星高度計(jì)SLA相關(guān)系數(shù)相對(duì)較高,相關(guān)性較強(qiáng), 這就提醒我們, 在我們使用PIES校驗(yàn)海洋二號(hào)(HY-2)衛(wèi)星數(shù)據(jù)時(shí), 可以選擇PIES布放在高能區(qū)和斜壓部分對(duì)海面貢獻(xiàn)較大的區(qū)域來(lái)校驗(yàn)。

    圖7 全部PIES站點(diǎn)PIES-SLA和Ref-MSLA的相關(guān)系數(shù)示意圖Fig. 7 Correlation coefficient of the gridded Ref-MSLA and the PIES-SLA

    下一步工作可以考慮在日本海區(qū)域利用PIES數(shù)據(jù)來(lái)進(jìn)行時(shí)間和空間相關(guān)性尺度的分析(temporal and spatial correlation scale)。

    [1] Rossby T. On monitoring depth variations of the main thermocline acoustically[J]. Journal of Geophysical Research Atmospheres, 1969, 74(23): 5542-5546.

    [2] Meinen C S, Watts D R. Vertical structure and transport on a transect across the North Atlantic Current near 42°N: Time series and mean[J]. Journal of Geophysical Research, 2000, 105(C9): 21869-21891.

    [3] Sun Che, and Watts D R. Acircumpolar gravest empirical mode for the Southern Ocean hydrography[J]. Journal of Geophysical Research, 2001, 106(C2): 2833- 2855.

    [4] Watts D R, Qian Xiaoshu, Tracey K L. Mapping abyssal current and pressure fields under the meandering Gulf Stream[J]. Journal of Atmospheric and Oceanic Technology, 2001, 18(6): 1052-1067.

    [5] Ponte R M. A preliminary model study of the largescale seasonal cycle in bottom pressure over the global ocean[J]. Journal of Geophysical Research, 1999, 104(C1): 1289-1300.

    [6] Guinehut S, Traon P-YL, Larnicol G. What can we learn from global altimetry/hydrography comparisons?[J]. Geophysical Research Letters, 2006, 33(10): 229-237.

    [7] Jayne S R, Wahr J M, Bryan F O. Observing ocean heat content using satellite gravity and altimetry[J]. Journal of Geophysical Research Atmospheres, 2003, 108(C2): 295-314.

    [8] Baker-Yeboah S, Watts D R, Byrne D A, et al. Measurements of sea surface height variability in the Easten South Atlantu from pressure sensor-equipped inverted echo sounders: baroclinic and barotropù[J]. Journal of Atmospheric & Oceanic Technology, 2009, 26(12): 2593.

    [9] Fukumori I, Raghunath R, Fu L-L. The nature of global largescale sea-level variability in relation to atmospheric forcing: a modeling study[J]. Journal of Geophysical Research, 1998, 103: 5493-5512.

    [10] Gilson J, Roemmich D, Cornuelle B. Relationship of TOPEX/Poseidon altimetric height to steric height and circulation in theNorth Pacific[J]. Journal of Geophysical Research, 1998, 103(C12): 27947-27965.

    [11] McCarthy M C, Talley L D, Roemmich D. Seasonal to interannual variability from expendable bathythermographand TOPEX/Poseidon altimeter data in the South Pacific subtropicalgyre[J]. Journal of Geophysical Research, 2000, 105(C8): 19535-19550

    [12] Xu Yongsheng, Watts D R, Wimbush M. Coupled patterns between fields of dynamic height and bottom pressure in the Japan/East Sea[J]. Ocean Science Journal, 2009, 44(1): 35-42.

    [13] Xu Yongsheng, Watts D R, Wimbush M, et al. Fundamentalmode basin oscillations in the Japan/East Sea[J]. Geophysical Research Letters, 2007, 34(4): 545-559.

    [14] Xu Yongsheng, Watts D R, Park J H. De-aliasing of large-scale high-frequency barotropic signals from satellite altimetry in the Japan/East Sea[J]. Journal of Atmospheric and Oceanic Technology, 2008, 25(9): 1703-1709.

    [15] Park J H, Watts D R, Donohue K A, et al. Comparisons of sea surface height variability observed by pressurerecording inverted echo sounders and satellite altimetry in the Kuroshio Extension[J]. Journal of Oceanography, 2012, 68(3): 401-416.

    [16] Mitchell D A, Wimbush M, Watts D R, et al. The residual GEM technique and its application to the southwestern Japan/East sea[J]. Journal of Atmospheric and Oceanic Technology, 2004, 21(12): 1895-1909.

    [17] Teague W J, Hallock Z R, Jacobs G A, et al. Kuroshio sea surface height fluctuations observed simultaneously with inverted echo sounders and TOPEX/Poseidon[J]. Journal of Geophysical Research, 1995, 100(C12): 24987–24994.

    [18] Book J W, Wimbush M, Imawaki S, et al. Kuroshio temporal and spatial variations south of Japan determined from inverted echo sounder measurements[J]. Journal of Geophysical Research-Oceans, 2002, 107(C9): 4-1–4-12.

    (本文編輯: 劉珊珊)

    Comparison of sea surface height anomalies derived by pressure-sensor-equipped inverted echo sounders and satellite altimetry in the Southwest Japan/East Sea

    GE Lei1, 2, 3, XU Yong-sheng1, 3, 4, YIN Bao-shu1, 3, 4
    (1. Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. Key Laboratory of Ocean Circulation and Waves, Chinese Academy of Sciences, Qingdao 266071, China; 4. Laboratory of Ocean and Climate Dynamies, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266237, China)

    Mar., 28, 2013

    sea surface height; pressure-sensor-equipped inverted echo sounder (PIES); satellite altimetry; acoustic echo time; ocean bottom pressure

    Pressure-sensor-equipped inverted echo sounders (PIES) can be used to measure sea bottom pressure and the acoustic-wave propagation time from the seabed to the surface. Bottom pressure and the acoustic propagation time are used to estimate the contribution of the sea surface height anomaly (SLA) to water quality changes (positive pressure) and hematocrit changes (baroclinic), respectively. In this paper, we compare SLAs derived by PIES with SLAs derived from satellite altimetry data in the Southwest Japan/East Sea. Using correlation analysis, we compared the PIES SLAs with SLAs measured by the along-track T/P satellite (TP SLA) and the ERS-2 satellite (ERS-2 SLA). We then compared the PIES SLA and AVISO gridded SLAs to estimate possible sources of error. We also analyzed the contributions of the PIES SLA positive pressure and baroclinic aspects of the SLA. The comparison results show that the correlation coefficient of the PIES SLA and Sat SLA is relatively high, and the root mean square error is relatively small. Further, we studied data from specific regions to identify possible causes of site-specific errors. Based on our results, we draw the following conclusions: (1) with respect to the Gulf Stream and the Kuroshio region, the contribution of positive sea surface height is relatively large; (2) if we consider the contribution of baroclinic changes in sea surface height, the PIES SLA and Sat SLA correlation coefficient will improve; (3) in the high-energy zone, the PIES SLA and Sat SLA correlation coefficient is relatively high and relatively better. Overall, in the Japan/East Sea, the PIES SLA and satellite altimeter SLA correlation coefficient is relatively high with a high level of consistency, providing a reliable way for checking the operation of Ocean II satellite (HY-2) altimeters. The research results have significance for the development and design of PIES, and for the selection of the placement of PIES.

    P733.1

    A

    1000-3096(2016)02-0128-08

    10.11759/hykx20130328001

    2013-03-28;

    2013-04-04

    中國(guó)科學(xué)院“百人計(jì)劃”(Y32109101L); 國(guó)家自然科學(xué)基金項(xiàng)目(41376028); 中國(guó)科學(xué)院方向項(xiàng)目(KZCX2-EW-209); 國(guó)家海洋局

    葛磊(1987-), 男, 河南商丘人, 碩士研究生, 電話: 15806599166, E-mail: gexinlei@126.com; 徐永生(1970-), 通信作者, 男, 山東青島人, 博士, 研究員, 從事衛(wèi)星海洋學(xué), 海洋遙感, 海洋動(dòng)力學(xué)等方面的研究工作, 電話: 86-13656488817, E-mail: yongsheng.xu@qdio.ac.cn ; 尹寶樹(1963-), 男, 山東青島人, 博士, 研究員, 主要從事海洋動(dòng)力過(guò)程及數(shù)值模擬和海洋災(zāi)害研究, 電話: 86-532-82898502, E-mail: bsyin@qdio.ac.cn

    猜你喜歡
    斜壓高度計(jì)距平
    颶風(fēng)Edouard(2014)暖心結(jié)構(gòu)的多資料對(duì)比分析
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計(jì)有效波高研究
    MIMU/GNSS/ODO/高度計(jì)/航姿儀組合導(dǎo)航微系統(tǒng)硬件設(shè)計(jì)
    航天控制(2020年5期)2020-03-29 02:10:36
    同化衛(wèi)星高度計(jì)觀測(cè)對(duì)CAS-ESM-C上層海洋溫度模擬的改進(jìn)
    斜壓渦度的變化與臺(tái)風(fēng)暴雨的關(guān)系研究
    近40年阿里地區(qū)云量和氣溫的年際變化
    西藏科技(2018年9期)2018-10-17 05:51:30
    合成孔徑雷達(dá)高度計(jì)與傳統(tǒng)高度計(jì)精度比對(duì)分析與機(jī)載試驗(yàn)驗(yàn)證
    甘肅省降水和冰雹天氣氣候分析
    Rossby波的線性穩(wěn)定性
    科技資訊(2016年23期)2016-05-30 20:25:49
    2005—2009年、2011年和2013年南海東北部120°E斷面秋季體積輸運(yùn)的年際變化*
    海洋與湖沼(2016年1期)2016-01-15 03:50:46
    国产亚洲精品av在线| 国产精品免费视频内射| 欧美一区二区精品小视频在线| 搡老岳熟女国产| 亚洲av成人av| 中国美女看黄片| av片东京热男人的天堂| 日韩免费av在线播放| 欧美亚洲日本最大视频资源| 国产精品免费一区二区三区在线| 国产精品 国内视频| 久久久久久亚洲精品国产蜜桃av| 波多野结衣巨乳人妻| 在线观看66精品国产| 精品久久久久久,| 成人亚洲精品一区在线观看| 一级作爱视频免费观看| 国产激情偷乱视频一区二区| 好男人电影高清在线观看| 国产一卡二卡三卡精品| 国产激情欧美一区二区| 国语自产精品视频在线第100页| 欧美一区二区精品小视频在线| 欧美黄色片欧美黄色片| 免费观看人在逋| 亚洲五月婷婷丁香| 一进一出抽搐动态| 国产视频内射| 狂野欧美激情性xxxx| 国产片内射在线| 色播在线永久视频| 国产午夜福利久久久久久| 国产三级在线视频| 老汉色av国产亚洲站长工具| 一级片免费观看大全| 三级毛片av免费| 99久久国产精品久久久| 日本三级黄在线观看| 真人做人爱边吃奶动态| 黄色成人免费大全| 男人舔女人下体高潮全视频| 1024香蕉在线观看| 免费看美女性在线毛片视频| av视频在线观看入口| 欧美色视频一区免费| 欧美成人免费av一区二区三区| 国产成人av激情在线播放| 伦理电影免费视频| 在线国产一区二区在线| 悠悠久久av| 国产不卡一卡二| 亚洲精品久久成人aⅴ小说| 久久国产精品男人的天堂亚洲| 亚洲专区国产一区二区| 成人欧美大片| 国产aⅴ精品一区二区三区波| 正在播放国产对白刺激| 在线播放国产精品三级| 成人18禁在线播放| 日韩欧美国产在线观看| 欧美大码av| 丰满人妻熟妇乱又伦精品不卡| www.精华液| 法律面前人人平等表现在哪些方面| 麻豆av在线久日| 亚洲男人天堂网一区| 又黄又粗又硬又大视频| 天堂影院成人在线观看| 1024手机看黄色片| 亚洲五月天丁香| 俄罗斯特黄特色一大片| 美女高潮到喷水免费观看| 国产欧美日韩一区二区三| 亚洲国产毛片av蜜桃av| 伊人久久大香线蕉亚洲五| 亚洲 欧美 日韩 在线 免费| 女性被躁到高潮视频| 人人妻人人看人人澡| 香蕉久久夜色| 2021天堂中文幕一二区在线观 | 日韩欧美国产一区二区入口| 亚洲久久久国产精品| 日本a在线网址| 麻豆成人午夜福利视频| 正在播放国产对白刺激| 欧美 亚洲 国产 日韩一| 欧美亚洲日本最大视频资源| 香蕉av资源在线| 国产精品自产拍在线观看55亚洲| 亚洲色图 男人天堂 中文字幕| 国产又黄又爽又无遮挡在线| 中文亚洲av片在线观看爽| 成人精品一区二区免费| 亚洲成人免费电影在线观看| 国产精品久久久久久亚洲av鲁大| 中文字幕另类日韩欧美亚洲嫩草| 成人亚洲精品av一区二区| 757午夜福利合集在线观看| 国产在线观看jvid| 欧美成人性av电影在线观看| 久久伊人香网站| 侵犯人妻中文字幕一二三四区| 久久精品aⅴ一区二区三区四区| 国产成人精品无人区| 色av中文字幕| 亚洲 欧美 日韩 在线 免费| 欧美日韩中文字幕国产精品一区二区三区| 欧美色欧美亚洲另类二区| 欧美日本视频| 村上凉子中文字幕在线| 亚洲精品久久国产高清桃花| 一级毛片女人18水好多| 亚洲真实伦在线观看| 黄色视频,在线免费观看| 最新在线观看一区二区三区| 欧美精品亚洲一区二区| 国产私拍福利视频在线观看| 亚洲成人久久爱视频| 在线国产一区二区在线| 一a级毛片在线观看| 亚洲国产精品成人综合色| 日本三级黄在线观看| 国产视频一区二区在线看| 国产精品免费一区二区三区在线| 免费在线观看黄色视频的| 麻豆一二三区av精品| 欧美又色又爽又黄视频| 欧美在线一区亚洲| www.熟女人妻精品国产| 国产爱豆传媒在线观看 | 午夜福利在线观看吧| 黑人欧美特级aaaaaa片| 波多野结衣高清无吗| 嫩草影院精品99| 狂野欧美激情性xxxx| 亚洲色图av天堂| 国产av一区二区精品久久| 亚洲五月天丁香| 久久人妻av系列| 少妇粗大呻吟视频| 日韩av在线大香蕉| 人妻丰满熟妇av一区二区三区| 18禁观看日本| 久久久久久久精品吃奶| 国产在线精品亚洲第一网站| 精品久久久久久,| 精品久久久久久久末码| 久久久久久久久中文| 美女午夜性视频免费| 丝袜人妻中文字幕| 亚洲av成人一区二区三| 国产一卡二卡三卡精品| 夜夜夜夜夜久久久久| 成年版毛片免费区| 色综合站精品国产| 天堂影院成人在线观看| 久久久久亚洲av毛片大全| 亚洲精品中文字幕一二三四区| 国产片内射在线| 午夜福利高清视频| 午夜亚洲福利在线播放| 国产人伦9x9x在线观看| 精品电影一区二区在线| 成人午夜高清在线视频 | a级毛片a级免费在线| 久久久久国产一级毛片高清牌| 又紧又爽又黄一区二区| 亚洲中文av在线| 欧美中文日本在线观看视频| 国产99白浆流出| 亚洲无线在线观看| 观看免费一级毛片| 亚洲国产精品sss在线观看| 中国美女看黄片| 一级毛片高清免费大全| 国产蜜桃级精品一区二区三区| 18禁裸乳无遮挡免费网站照片 | 精品国产一区二区三区四区第35| 1024视频免费在线观看| 欧美乱色亚洲激情| 黄色视频不卡| 国产人伦9x9x在线观看| 丁香欧美五月| 好男人在线观看高清免费视频 | 亚洲午夜理论影院| 人人妻,人人澡人人爽秒播| 色播亚洲综合网| 国产精品一区二区免费欧美| 国产亚洲精品第一综合不卡| 少妇裸体淫交视频免费看高清 | 老司机午夜十八禁免费视频| 国产免费男女视频| 他把我摸到了高潮在线观看| 真人一进一出gif抽搐免费| 日日爽夜夜爽网站| 成人三级黄色视频| 欧美乱妇无乱码| 99精品在免费线老司机午夜| 97碰自拍视频| 亚洲av熟女| 国产激情久久老熟女| 亚洲国产高清在线一区二区三 | 欧美成狂野欧美在线观看| 成人18禁高潮啪啪吃奶动态图| 精品国产超薄肉色丝袜足j| 亚洲成av片中文字幕在线观看| 欧美乱妇无乱码| 国产亚洲精品一区二区www| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全免费视频| √禁漫天堂资源中文www| 精品不卡国产一区二区三区| 琪琪午夜伦伦电影理论片6080| www国产在线视频色| 欧美不卡视频在线免费观看 | 免费高清在线观看日韩| 美女扒开内裤让男人捅视频| 又紧又爽又黄一区二区| 国产亚洲欧美在线一区二区| a级毛片在线看网站| 俺也久久电影网| 成人国产一区最新在线观看| 成人三级黄色视频| 亚洲精品av麻豆狂野| 岛国在线观看网站| www.自偷自拍.com| 国产精品久久久av美女十八| 国产真人三级小视频在线观看| 日韩欧美一区二区三区在线观看| 天堂√8在线中文| 欧美zozozo另类| 在线观看免费午夜福利视频| 夜夜夜夜夜久久久久| 制服人妻中文乱码| 亚洲国产日韩欧美精品在线观看 | 久久青草综合色| 亚洲人成网站在线播放欧美日韩| 精华霜和精华液先用哪个| 99精品欧美一区二区三区四区| 首页视频小说图片口味搜索| 午夜两性在线视频| 叶爱在线成人免费视频播放| 黑人欧美特级aaaaaa片| 一级a爱片免费观看的视频| 久久久国产精品麻豆| 成人三级做爰电影| 国产成人欧美| 中国美女看黄片| 黄色毛片三级朝国网站| 亚洲成人精品中文字幕电影| 操出白浆在线播放| 国产1区2区3区精品| 久久人人精品亚洲av| 一区二区三区精品91| 看黄色毛片网站| 日韩精品免费视频一区二区三区| 亚洲成av人片免费观看| 亚洲人成网站在线播放欧美日韩| 90打野战视频偷拍视频| 国产成人啪精品午夜网站| 丝袜在线中文字幕| 亚洲成av人片免费观看| 午夜日韩欧美国产| 香蕉av资源在线| 日本一区二区免费在线视频| 亚洲专区国产一区二区| 欧美日韩乱码在线| 制服人妻中文乱码| 嫩草影视91久久| 亚洲精品中文字幕一二三四区| 日本在线视频免费播放| 久久久国产欧美日韩av| 又大又爽又粗| 在线观看舔阴道视频| 国产又爽黄色视频| 可以免费在线观看a视频的电影网站| 婷婷精品国产亚洲av| 成人免费观看视频高清| 精品久久久久久成人av| 精品国产一区二区三区四区第35| aaaaa片日本免费| 欧美性长视频在线观看| 免费女性裸体啪啪无遮挡网站| 在线看三级毛片| 久久青草综合色| 成人国产一区最新在线观看| 黄色视频,在线免费观看| 一本精品99久久精品77| 91成人精品电影| 午夜亚洲福利在线播放| 亚洲专区字幕在线| 国产亚洲欧美在线一区二区| 男女那种视频在线观看| 欧美日本亚洲视频在线播放| bbb黄色大片| 美女国产高潮福利片在线看| 亚洲第一青青草原| 亚洲国产中文字幕在线视频| 午夜免费成人在线视频| 非洲黑人性xxxx精品又粗又长| 亚洲国产精品合色在线| 久久久国产精品麻豆| 精品国产乱子伦一区二区三区| 久久国产乱子伦精品免费另类| 日韩有码中文字幕| 国产精品免费视频内射| 99国产综合亚洲精品| 久久这里只有精品19| 变态另类成人亚洲欧美熟女| 亚洲九九香蕉| 精品国产国语对白av| 一进一出好大好爽视频| 少妇的丰满在线观看| 此物有八面人人有两片| 麻豆av在线久日| 国产亚洲精品一区二区www| 琪琪午夜伦伦电影理论片6080| 国产激情久久老熟女| 精品电影一区二区在线| 成人午夜高清在线视频 | 欧美激情极品国产一区二区三区| 中出人妻视频一区二区| 日韩成人在线观看一区二区三区| 亚洲色图 男人天堂 中文字幕| 在线十欧美十亚洲十日本专区| 亚洲成人久久性| 黄网站色视频无遮挡免费观看| 桃色一区二区三区在线观看| 两个人看的免费小视频| 老司机靠b影院| 欧美日韩乱码在线| 老司机福利观看| 国产一区二区三区视频了| 国产蜜桃级精品一区二区三区| 18禁美女被吸乳视频| 麻豆成人午夜福利视频| 制服诱惑二区| 亚洲av片天天在线观看| 婷婷精品国产亚洲av在线| 亚洲aⅴ乱码一区二区在线播放 | 变态另类丝袜制服| 欧美性猛交╳xxx乱大交人| 俺也久久电影网| 一卡2卡三卡四卡精品乱码亚洲| 国产在线精品亚洲第一网站| 日韩欧美国产一区二区入口| 极品教师在线免费播放| 精品午夜福利视频在线观看一区| 精品少妇一区二区三区视频日本电影| 岛国视频午夜一区免费看| 一本久久中文字幕| 免费在线观看日本一区| 18禁观看日本| 欧美日韩中文字幕国产精品一区二区三区| 久久精品成人免费网站| 成人一区二区视频在线观看| 国产精品综合久久久久久久免费| 色综合站精品国产| 好男人在线观看高清免费视频 | 在线观看舔阴道视频| 在线视频色国产色| tocl精华| 黄片播放在线免费| 亚洲全国av大片| 俄罗斯特黄特色一大片| 中文在线观看免费www的网站 | 国产亚洲精品久久久久5区| 啦啦啦韩国在线观看视频| 99re在线观看精品视频| 午夜久久久在线观看| 一区二区三区高清视频在线| а√天堂www在线а√下载| 在线观看免费日韩欧美大片| 亚洲av日韩精品久久久久久密| 99精品久久久久人妻精品| 久久人妻av系列| 午夜影院日韩av| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久| 国产精品国产高清国产av| 他把我摸到了高潮在线观看| 18美女黄网站色大片免费观看| 天堂影院成人在线观看| 亚洲一区中文字幕在线| 亚洲国产欧洲综合997久久, | 成人精品一区二区免费| 最好的美女福利视频网| av视频在线观看入口| 成人欧美大片| 免费在线观看影片大全网站| 琪琪午夜伦伦电影理论片6080| 国产麻豆成人av免费视频| 一二三四在线观看免费中文在| 99精品欧美一区二区三区四区| 日本一区二区免费在线视频| 亚洲一区二区三区不卡视频| 日韩欧美一区二区三区在线观看| 18美女黄网站色大片免费观看| 久久热在线av| 男人舔奶头视频| 午夜亚洲福利在线播放| 精品久久久久久久久久免费视频| 日本精品一区二区三区蜜桃| 国产精品永久免费网站| 搡老熟女国产l中国老女人| 人妻久久中文字幕网| 国产精品国产高清国产av| 欧美av亚洲av综合av国产av| 99国产综合亚洲精品| av片东京热男人的天堂| 欧美中文日本在线观看视频| 手机成人av网站| www.999成人在线观看| e午夜精品久久久久久久| 黄网站色视频无遮挡免费观看| 女生性感内裤真人,穿戴方法视频| 久久久国产成人精品二区| 亚洲一区二区三区色噜噜| 成人亚洲精品一区在线观看| 午夜成年电影在线免费观看| 精品人妻1区二区| 欧美黑人巨大hd| 看片在线看免费视频| 亚洲激情在线av| 精品熟女少妇八av免费久了| 国产成+人综合+亚洲专区| 一区二区日韩欧美中文字幕| 午夜福利在线观看吧| 欧美三级亚洲精品| 在线观看免费日韩欧美大片| 日本免费一区二区三区高清不卡| 一级作爱视频免费观看| 国产男靠女视频免费网站| 热99re8久久精品国产| 亚洲成人精品中文字幕电影| 亚洲一区高清亚洲精品| 久久婷婷人人爽人人干人人爱| 亚洲av电影不卡..在线观看| 国产男靠女视频免费网站| 婷婷丁香在线五月| 成人特级黄色片久久久久久久| 成人手机av| 国产精品1区2区在线观看.| 男人的好看免费观看在线视频 | 国产黄片美女视频| 中亚洲国语对白在线视频| 亚洲免费av在线视频| www.精华液| 亚洲一区中文字幕在线| 日本 av在线| 母亲3免费完整高清在线观看| 国产亚洲精品av在线| 99re在线观看精品视频| 好男人在线观看高清免费视频 | 久久久国产精品麻豆| 国产精品电影一区二区三区| 1024香蕉在线观看| 精品高清国产在线一区| 久久久水蜜桃国产精品网| netflix在线观看网站| 午夜福利免费观看在线| 曰老女人黄片| 国产黄色小视频在线观看| av免费在线观看网站| 欧美黑人巨大hd| 亚洲精品国产区一区二| 国产亚洲av嫩草精品影院| 韩国精品一区二区三区| 老司机午夜十八禁免费视频| 男女下面进入的视频免费午夜 | 老熟妇仑乱视频hdxx| 青草久久国产| 免费无遮挡裸体视频| cao死你这个sao货| 欧美成人午夜精品| 精品欧美国产一区二区三| 中文字幕人妻丝袜一区二区| 在线国产一区二区在线| 日韩av在线大香蕉| 亚洲一区二区三区色噜噜| x7x7x7水蜜桃| 男人舔奶头视频| 色av中文字幕| 午夜a级毛片| 高潮久久久久久久久久久不卡| 黄色女人牲交| 久久 成人 亚洲| 黄色视频不卡| 欧美国产精品va在线观看不卡| 黄频高清免费视频| 亚洲精品在线美女| 1024手机看黄色片| 亚洲精品国产区一区二| 精品熟女少妇八av免费久了| 国产在线观看jvid| 亚洲男人的天堂狠狠| 岛国视频午夜一区免费看| 91老司机精品| 久久中文字幕人妻熟女| 国产亚洲精品一区二区www| 午夜福利在线观看吧| 白带黄色成豆腐渣| 欧美人与性动交α欧美精品济南到| av中文乱码字幕在线| 成熟少妇高潮喷水视频| 波多野结衣高清作品| 日日摸夜夜添夜夜添小说| 性色av乱码一区二区三区2| 亚洲最大成人中文| 制服诱惑二区| 欧美成人午夜精品| 国产日本99.免费观看| 欧美黑人精品巨大| 亚洲人成网站高清观看| 91av网站免费观看| 国产精品乱码一区二三区的特点| 青草久久国产| 久久久久久九九精品二区国产 | 91大片在线观看| 国产三级在线视频| 成人18禁在线播放| 精品乱码久久久久久99久播| 韩国av一区二区三区四区| 中文字幕另类日韩欧美亚洲嫩草| 欧美三级亚洲精品| 可以在线观看的亚洲视频| avwww免费| 日本三级黄在线观看| 久久久国产精品麻豆| 香蕉久久夜色| 亚洲五月婷婷丁香| 999精品在线视频| 国产极品粉嫩免费观看在线| 中文字幕人成人乱码亚洲影| 国产av一区在线观看免费| 视频区欧美日本亚洲| 国产亚洲精品一区二区www| 亚洲性夜色夜夜综合| 侵犯人妻中文字幕一二三四区| 成人手机av| 国产成人欧美| 亚洲男人的天堂狠狠| 18禁国产床啪视频网站| 好看av亚洲va欧美ⅴa在| 51午夜福利影视在线观看| 免费无遮挡裸体视频| 人人妻人人澡人人看| 99国产极品粉嫩在线观看| aaaaa片日本免费| 亚洲成人国产一区在线观看| 99精品在免费线老司机午夜| 日本a在线网址| 亚洲欧美精品综合一区二区三区| 青草久久国产| 国产一区在线观看成人免费| 中文字幕最新亚洲高清| 国产精品av久久久久免费| 女性被躁到高潮视频| 亚洲五月天丁香| av中文乱码字幕在线| 老司机福利观看| 久久精品人妻少妇| 国产高清videossex| 国产精品一区二区三区四区久久 | 午夜福利在线在线| 成人18禁在线播放| 岛国视频午夜一区免费看| 国产成人av激情在线播放| 国产成人精品无人区| 精品少妇一区二区三区视频日本电影| 欧美黄色淫秽网站| 日韩精品免费视频一区二区三区| 国产精品一区二区免费欧美| 狠狠狠狠99中文字幕| 变态另类成人亚洲欧美熟女| 成人手机av| 国产又色又爽无遮挡免费看| 制服人妻中文乱码| 日日干狠狠操夜夜爽| 黄色成人免费大全| 亚洲色图 男人天堂 中文字幕| 国产精品香港三级国产av潘金莲| 国产aⅴ精品一区二区三区波| 长腿黑丝高跟| 亚洲色图av天堂| 青草久久国产| av欧美777| 精品久久久久久久久久免费视频| 亚洲第一av免费看| 国产精品一区二区免费欧美| 女生性感内裤真人,穿戴方法视频| 成人欧美大片| 18禁裸乳无遮挡免费网站照片 | 19禁男女啪啪无遮挡网站| 中文资源天堂在线| 久久午夜亚洲精品久久| 757午夜福利合集在线观看| 又大又爽又粗| 成年版毛片免费区| 成年免费大片在线观看| 97碰自拍视频| 麻豆国产av国片精品| 日韩欧美一区二区三区在线观看| 一级毛片女人18水好多| 国产精品久久视频播放| 国产99白浆流出| 欧美乱妇无乱码| 国产男靠女视频免费网站| 国产精品亚洲av一区麻豆| 国产欧美日韩一区二区三| 亚洲av片天天在线观看| 免费av毛片视频| 欧美精品啪啪一区二区三区| 俺也久久电影网| 国产精品香港三级国产av潘金莲|