• <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高清涩受| 9色porny在线观看| 久久人人爽人人片av| av在线老鸭窝| 亚洲av日韩在线播放| 免费在线观看黄色视频的| 国产片内射在线| 成年人午夜在线观看视频| 十八禁高潮呻吟视频| 美女高潮到喷水免费观看| 欧美日韩成人在线一区二区| 美女扒开内裤让男人捅视频| a 毛片基地| 国产在线观看jvid| 男女之事视频高清在线观看 | 一级毛片 在线播放| 男女国产视频网站| 一级黄色大片毛片| 人妻 亚洲 视频| 欧美xxⅹ黑人| 国产亚洲欧美精品永久| 免费日韩欧美在线观看| tube8黄色片| 午夜福利免费观看在线| 免费不卡黄色视频| 夜夜骑夜夜射夜夜干| 天天躁夜夜躁狠狠躁躁| 国产成人精品在线电影| 日本一区二区免费在线视频| 王馨瑶露胸无遮挡在线观看| 极品少妇高潮喷水抽搐| 欧美日韩亚洲综合一区二区三区_| 黄色片一级片一级黄色片| 午夜免费鲁丝| 亚洲精品中文字幕在线视频| 国产伦理片在线播放av一区| 男人操女人黄网站| 男男h啪啪无遮挡| 亚洲国产成人一精品久久久| 日本黄色日本黄色录像| 国产精品秋霞免费鲁丝片| 亚洲精品av麻豆狂野| 美女视频免费永久观看网站| 国产在线免费精品| 两人在一起打扑克的视频| 精品国产国语对白av| 老熟女久久久| 黄色视频不卡| 久久这里只有精品19| 一区二区三区四区激情视频| 日韩制服丝袜自拍偷拍| 国产一区二区激情短视频 | 两个人看的免费小视频| 老司机靠b影院| 99re6热这里在线精品视频| 最新的欧美精品一区二区| 男女国产视频网站| 国产精品一二三区在线看| 80岁老熟妇乱子伦牲交| 1024视频免费在线观看| 一边亲一边摸免费视频| 老司机在亚洲福利影院| 久久毛片免费看一区二区三区| 老鸭窝网址在线观看| 精品久久蜜臀av无| 亚洲一区中文字幕在线| 久久天堂一区二区三区四区| 两性夫妻黄色片| 久久国产精品人妻蜜桃| 美女视频免费永久观看网站| 中文欧美无线码| 欧美中文综合在线视频| 精品国产超薄肉色丝袜足j| 国产91精品成人一区二区三区 | 黄频高清免费视频| 午夜激情久久久久久久| 手机成人av网站| 成人午夜精彩视频在线观看| 成人手机av| 国产一区二区 视频在线| 久久久久久免费高清国产稀缺| 免费在线观看完整版高清| 十八禁人妻一区二区| 9色porny在线观看| 亚洲伊人色综图| 欧美日韩综合久久久久久| 99热网站在线观看| 最近手机中文字幕大全| 80岁老熟妇乱子伦牲交| 婷婷色麻豆天堂久久| 欧美大码av| 韩国高清视频一区二区三区| 国产一区二区激情短视频 | 一区二区三区乱码不卡18| 少妇精品久久久久久久| 中国美女看黄片| 中国国产av一级| 丰满饥渴人妻一区二区三| www.999成人在线观看| 久久九九热精品免费| 黄色a级毛片大全视频| 亚洲国产成人一精品久久久| 男女之事视频高清在线观看 | 久久精品国产a三级三级三级| 亚洲成色77777| 十分钟在线观看高清视频www| 女性被躁到高潮视频| 久久久精品免费免费高清| 免费看十八禁软件| 91精品伊人久久大香线蕉| 新久久久久国产一级毛片| 国产伦人伦偷精品视频| 成人18禁高潮啪啪吃奶动态图| 一区二区三区激情视频| 成人亚洲欧美一区二区av| 亚洲成国产人片在线观看| 亚洲一区中文字幕在线| 亚洲成人国产一区在线观看 | 在线观看人妻少妇| 亚洲av日韩在线播放| 国产片内射在线| 国产无遮挡羞羞视频在线观看| 19禁男女啪啪无遮挡网站| 极品人妻少妇av视频| 国产精品熟女久久久久浪| 少妇被粗大的猛进出69影院| 午夜老司机福利片| 老司机在亚洲福利影院| 国产一区二区激情短视频 | 精品少妇内射三级| h视频一区二区三区| 精品福利永久在线观看| 国产精品一区二区免费欧美 | 免费高清在线观看日韩| 免费av中文字幕在线| 在线观看www视频免费| 岛国毛片在线播放| 国产精品免费大片| 日韩中文字幕欧美一区二区 | 在线 av 中文字幕| 成年人免费黄色播放视频| 亚洲国产av影院在线观看| 国产精品熟女久久久久浪| 亚洲视频免费观看视频| 亚洲伊人久久精品综合| 97精品久久久久久久久久精品| 亚洲久久久国产精品| 精品国产一区二区久久| 九色亚洲精品在线播放| 免费高清在线观看日韩| 久久国产精品人妻蜜桃| 欧美精品高潮呻吟av久久| 久久精品成人免费网站| xxxhd国产人妻xxx| 国产在线观看jvid| 欧美日韩一级在线毛片| 手机成人av网站| 99re6热这里在线精品视频| av片东京热男人的天堂| 妹子高潮喷水视频| 亚洲欧美精品自产自拍| 天堂8中文在线网| 999久久久国产精品视频| 欧美激情高清一区二区三区| 我的亚洲天堂| 青春草视频在线免费观看| 亚洲五月婷婷丁香| 老司机深夜福利视频在线观看 | 丝袜人妻中文字幕| 一级毛片我不卡| 人妻 亚洲 视频| 99精品久久久久人妻精品| 一本久久精品| 国产欧美日韩精品亚洲av| 国产成人免费观看mmmm| 日韩 亚洲 欧美在线| www.av在线官网国产| 中文字幕色久视频| 蜜桃国产av成人99| 伊人亚洲综合成人网| 午夜免费鲁丝| 天天躁狠狠躁夜夜躁狠狠躁| 在线观看免费午夜福利视频| 久久久久国产精品人妻一区二区| 午夜福利在线免费观看网站| 成人18禁高潮啪啪吃奶动态图| 亚洲精品一区蜜桃| 两性夫妻黄色片| 午夜福利影视在线免费观看| 久久鲁丝午夜福利片| 18禁裸乳无遮挡动漫免费视频| 亚洲中文字幕日韩| 国产成人精品无人区| 午夜免费鲁丝| 国产精品久久久人人做人人爽| 中文字幕人妻熟女乱码| 欧美精品亚洲一区二区| 国产熟女欧美一区二区| 国产一区二区 视频在线| 97在线人人人人妻| 国产色视频综合| 国产在视频线精品| 国产在线免费精品| 国产激情久久老熟女| 纯流量卡能插随身wifi吗| 亚洲黑人精品在线| 热re99久久国产66热| h视频一区二区三区| 欧美日韩一级在线毛片| 九色亚洲精品在线播放| 亚洲成色77777| 丰满人妻熟妇乱又伦精品不卡| 成年av动漫网址| 国产精品 国内视频| 免费av中文字幕在线| 自线自在国产av| 亚洲精品自拍成人| 五月天丁香电影| 九草在线视频观看| 国产成人啪精品午夜网站| 狠狠精品人妻久久久久久综合| 久久久久久久久久久久大奶| 我要看黄色一级片免费的| 欧美黑人欧美精品刺激| 欧美中文综合在线视频| 男女午夜视频在线观看| 美女主播在线视频| 亚洲熟女毛片儿| 最近中文字幕2019免费版| 久久人人爽av亚洲精品天堂| 精品人妻在线不人妻| 亚洲伊人久久精品综合| 亚洲精品国产av蜜桃| 男女边摸边吃奶| 国语对白做爰xxxⅹ性视频网站| 久久亚洲国产成人精品v| 女人精品久久久久毛片| 男人舔女人的私密视频| 国产精品一国产av| 亚洲专区国产一区二区| 亚洲欧美精品自产自拍| 女人爽到高潮嗷嗷叫在线视频| 国产伦理片在线播放av一区| 久久久精品94久久精品| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产av新网站| 国产高清视频在线播放一区 | 777米奇影视久久| 飞空精品影院首页| 超色免费av| 中文乱码字字幕精品一区二区三区| 九草在线视频观看| 黄色a级毛片大全视频| 男女无遮挡免费网站观看| 国产国语露脸激情在线看| 久久久久国产精品人妻一区二区| 久久99一区二区三区| 黄色毛片三级朝国网站| 91麻豆精品激情在线观看国产 | 中文乱码字字幕精品一区二区三区| 美女高潮到喷水免费观看| 亚洲欧洲精品一区二区精品久久久| 十八禁高潮呻吟视频| 男人操女人黄网站| 国产亚洲一区二区精品| 1024视频免费在线观看| 黄色视频在线播放观看不卡| 精品免费久久久久久久清纯 | 高清黄色对白视频在线免费看| 男男h啪啪无遮挡| 亚洲男人天堂网一区| 免费久久久久久久精品成人欧美视频| 国产一区有黄有色的免费视频| 男人操女人黄网站| 啦啦啦啦在线视频资源| 又大又爽又粗| 国产精品一区二区精品视频观看| 大话2 男鬼变身卡| 欧美人与性动交α欧美软件| 叶爱在线成人免费视频播放| 久久精品国产综合久久久| 亚洲久久久国产精品| 欧美少妇被猛烈插入视频| 无遮挡黄片免费观看| 精品国产乱码久久久久久小说| 国产精品国产三级国产专区5o| 99国产精品一区二区蜜桃av | 色精品久久人妻99蜜桃| 夜夜骑夜夜射夜夜干| 日韩熟女老妇一区二区性免费视频| 久久热在线av| 日本欧美视频一区| 爱豆传媒免费全集在线观看| av在线app专区| 制服诱惑二区| 建设人人有责人人尽责人人享有的| 在线看a的网站| 欧美乱码精品一区二区三区| videosex国产| 啦啦啦视频在线资源免费观看| 热99久久久久精品小说推荐| 国产男女内射视频| 精品视频人人做人人爽| 一级片'在线观看视频| 国产野战对白在线观看| 色婷婷av一区二区三区视频| 97在线人人人人妻| 女人精品久久久久毛片| 国产亚洲av片在线观看秒播厂| 亚洲av成人精品一二三区| 国产无遮挡羞羞视频在线观看| 亚洲激情五月婷婷啪啪| 侵犯人妻中文字幕一二三四区| 成年美女黄网站色视频大全免费| 一本综合久久免费| 中文字幕另类日韩欧美亚洲嫩草| 无限看片的www在线观看| 久久久久精品人妻al黑| 日日夜夜操网爽| 深夜精品福利| 伦理电影免费视频| 国产91精品成人一区二区三区 | 精品久久蜜臀av无| 搡老岳熟女国产| 日韩大片免费观看网站| 少妇精品久久久久久久| 久久久久国产一级毛片高清牌| 久久精品国产a三级三级三级| 欧美日韩一级在线毛片| 18禁观看日本| 日本wwww免费看| 狠狠精品人妻久久久久久综合| 后天国语完整版免费观看| 国产精品一二三区在线看| 亚洲精品自拍成人| 亚洲成人免费电影在线观看 | 成人国产一区最新在线观看 | 日韩 欧美 亚洲 中文字幕| 成人国语在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久蜜臀av无| 日韩电影二区| www.熟女人妻精品国产| 免费看不卡的av| 欧美另类一区| a级片在线免费高清观看视频| 欧美另类一区| a级片在线免费高清观看视频| 亚洲精品国产色婷婷电影| 午夜免费观看性视频| 日本猛色少妇xxxxx猛交久久| 日韩 欧美 亚洲 中文字幕| av又黄又爽大尺度在线免费看| av在线老鸭窝| 高清av免费在线| 男人添女人高潮全过程视频| 国产亚洲午夜精品一区二区久久| 少妇被粗大的猛进出69影院| 久久影院123| 国产午夜精品一二区理论片| av天堂久久9| 欧美激情 高清一区二区三区| 久久影院123| 黄片小视频在线播放| 精品第一国产精品| 精品人妻熟女毛片av久久网站| 国产日韩一区二区三区精品不卡| 欧美人与善性xxx| 人人澡人人妻人| 深夜精品福利| 丁香六月欧美| 最近中文字幕2019免费版| 女人爽到高潮嗷嗷叫在线视频| 人体艺术视频欧美日本| 欧美性长视频在线观看| 美女大奶头黄色视频| a 毛片基地| 日韩人妻精品一区2区三区| 亚洲国产精品成人久久小说| 日韩一本色道免费dvd| 亚洲欧美中文字幕日韩二区| 日韩av不卡免费在线播放| 男人爽女人下面视频在线观看| 肉色欧美久久久久久久蜜桃| www日本在线高清视频| 久久99一区二区三区| 欧美 亚洲 国产 日韩一| 久久久国产精品麻豆| 国产1区2区3区精品| 黄色一级大片看看| 精品少妇一区二区三区视频日本电影| 一区二区三区精品91| 亚洲精品在线美女| 又大又黄又爽视频免费| 久久久国产欧美日韩av| 亚洲综合色网址| 欧美在线黄色| 精品一区在线观看国产| 亚洲成av片中文字幕在线观看| 久久人人爽av亚洲精品天堂| 亚洲中文av在线| 亚洲av片天天在线观看| 亚洲七黄色美女视频| 男人爽女人下面视频在线观看| 一级片免费观看大全| 精品少妇一区二区三区视频日本电影| 后天国语完整版免费观看| 纯流量卡能插随身wifi吗| 精品国产一区二区久久| 在线观看一区二区三区激情| 国产精品欧美亚洲77777| 热99久久久久精品小说推荐| 久久99热这里只频精品6学生| 五月开心婷婷网| 精品国产国语对白av| 狠狠婷婷综合久久久久久88av| www日本在线高清视频| 母亲3免费完整高清在线观看| 少妇 在线观看| 亚洲美女黄色视频免费看| 久久久久久人人人人人| 性少妇av在线| 日韩精品免费视频一区二区三区| 精品第一国产精品| 看十八女毛片水多多多| 热99国产精品久久久久久7| 亚洲av欧美aⅴ国产| 亚洲熟女毛片儿| 日本午夜av视频| 亚洲av在线观看美女高潮| 爱豆传媒免费全集在线观看| 赤兔流量卡办理| 久久人人爽av亚洲精品天堂| 人人妻人人澡人人爽人人夜夜| 亚洲精品国产av成人精品| 久久精品国产亚洲av高清一级| 黄色视频在线播放观看不卡| 香蕉国产在线看| 18禁观看日本| 亚洲国产欧美在线一区| 午夜免费男女啪啪视频观看| 人人澡人人妻人| 看十八女毛片水多多多| 欧美国产精品va在线观看不卡| 人成视频在线观看免费观看| 久久人人爽人人片av| 国产熟女午夜一区二区三区| 91成人精品电影| 久久久久国产一级毛片高清牌| 考比视频在线观看| 欧美日韩视频高清一区二区三区二| 精品一区在线观看国产| 精品少妇黑人巨大在线播放| 欧美乱码精品一区二区三区| 电影成人av| 久久久精品区二区三区| 亚洲伊人久久精品综合| 午夜福利乱码中文字幕| 人人澡人人妻人| 午夜福利免费观看在线| 亚洲国产精品成人久久小说| 午夜av观看不卡| 韩国精品一区二区三区| av又黄又爽大尺度在线免费看| 午夜免费观看性视频| 91老司机精品| 午夜福利在线免费观看网站| 欧美日韩av久久| 欧美亚洲 丝袜 人妻 在线| 女人被躁到高潮嗷嗷叫费观| 日韩大片免费观看网站| 1024香蕉在线观看| 丝袜美足系列| 国产女主播在线喷水免费视频网站| 国产亚洲av高清不卡| 中文字幕人妻熟女乱码| 波多野结衣av一区二区av| 激情五月婷婷亚洲| 九草在线视频观看| 色婷婷久久久亚洲欧美| 如日韩欧美国产精品一区二区三区| 老熟女久久久| 亚洲综合色网址| 亚洲,欧美,日韩| 91九色精品人成在线观看| 男的添女的下面高潮视频| 性色av一级| 色婷婷av一区二区三区视频| 久久毛片免费看一区二区三区| 亚洲,一卡二卡三卡| 成年美女黄网站色视频大全免费| 天堂俺去俺来也www色官网| 午夜老司机福利片| 在线观看免费视频网站a站| 91老司机精品| 肉色欧美久久久久久久蜜桃| 亚洲欧美中文字幕日韩二区| 国产成人影院久久av| 成年av动漫网址| 大片电影免费在线观看免费| 一二三四社区在线视频社区8| 亚洲欧美清纯卡通| 国产精品偷伦视频观看了| 97在线人人人人妻| 成人影院久久| 久久女婷五月综合色啪小说| 中文乱码字字幕精品一区二区三区| 成年人黄色毛片网站| 国产野战对白在线观看| 亚洲av电影在线观看一区二区三区| 波野结衣二区三区在线| 久久久久国产精品人妻一区二区| 精品国产一区二区久久| 考比视频在线观看| 大片免费播放器 马上看| a级毛片在线看网站| 亚洲 国产 在线| 国产精品免费大片| 看十八女毛片水多多多| 国产精品免费视频内射| 亚洲国产av新网站| 国产又色又爽无遮挡免| 亚洲国产最新在线播放| 免费高清在线观看视频在线观看| 国产日韩欧美亚洲二区| 国产日韩一区二区三区精品不卡| 亚洲九九香蕉| 水蜜桃什么品种好| 汤姆久久久久久久影院中文字幕| 热re99久久精品国产66热6| 操出白浆在线播放| 涩涩av久久男人的天堂| 午夜激情久久久久久久| 19禁男女啪啪无遮挡网站| 色婷婷久久久亚洲欧美| 最黄视频免费看| 国产xxxxx性猛交| 校园人妻丝袜中文字幕| 国产伦人伦偷精品视频| 免费看av在线观看网站| 亚洲国产欧美在线一区| 在线天堂中文资源库| 中文乱码字字幕精品一区二区三区| 69精品国产乱码久久久| 国产淫语在线视频| 91精品伊人久久大香线蕉| 久久天躁狠狠躁夜夜2o2o | 2018国产大陆天天弄谢| av有码第一页| 亚洲精品第二区| 国产成人免费观看mmmm| 中文字幕人妻丝袜一区二区| 欧美 亚洲 国产 日韩一| 日本色播在线视频| 一区福利在线观看| 真人做人爱边吃奶动态| 国产黄色免费在线视频| 精品人妻1区二区| 两个人看的免费小视频| 久久精品国产亚洲av高清一级| 欧美 日韩 精品 国产| 午夜久久久在线观看| 人人妻人人澡人人爽人人夜夜| 秋霞在线观看毛片| 18禁裸乳无遮挡动漫免费视频| 91精品三级在线观看| 久久精品亚洲熟妇少妇任你| 在线观看免费午夜福利视频| 赤兔流量卡办理| 一本久久精品| 69精品国产乱码久久久| 两个人免费观看高清视频| 亚洲成人手机| 久久精品亚洲av国产电影网| 亚洲国产毛片av蜜桃av| 国产99久久九九免费精品| 一级片免费观看大全| 日本五十路高清| 天堂8中文在线网| 欧美在线一区亚洲| 好男人电影高清在线观看| 搡老乐熟女国产| 精品免费久久久久久久清纯 | 国产精品免费大片| 老司机影院毛片| 欧美日韩黄片免| 一边亲一边摸免费视频| 国产男女超爽视频在线观看| 热re99久久精品国产66热6| 亚洲精品国产一区二区精华液| 90打野战视频偷拍视频| 精品国产超薄肉色丝袜足j| 国产亚洲欧美精品永久| 国产男女超爽视频在线观看| 超碰97精品在线观看| 国产精品三级大全| 99精国产麻豆久久婷婷| 免费在线观看影片大全网站 | 欧美在线黄色| 亚洲成人手机| 制服诱惑二区| 欧美日韩亚洲国产一区二区在线观看 | 99热网站在线观看| 免费少妇av软件| 激情五月婷婷亚洲| 黄色视频不卡| 9色porny在线观看| 欧美激情极品国产一区二区三区| 久久精品国产综合久久久| 香蕉丝袜av| 下体分泌物呈黄色| 成人国产av品久久久| 久久精品成人免费网站| 亚洲精品中文字幕在线视频|