• <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日韩在线播放| 国产1区2区3区精品| 9热在线视频观看99| 国产在线一区二区三区精| 一本一本久久a久久精品综合妖精| 男的添女的下面高潮视频| 七月丁香在线播放| 久久精品久久精品一区二区三区| 中文字幕制服av| 高清视频免费观看一区二区| 91国产中文字幕| 2018国产大陆天天弄谢| 高清欧美精品videossex| 欧美 日韩 精品 国产| 熟女av电影| 日韩人妻精品一区2区三区| 一区二区av电影网| 9191精品国产免费久久| 丝瓜视频免费看黄片| 人人妻人人爽人人添夜夜欢视频| 十八禁人妻一区二区| √禁漫天堂资源中文www| 亚洲四区av| 91成人精品电影| 亚洲av日韩精品久久久久久密 | 国产麻豆69| 久久99热这里只频精品6学生| 国产亚洲一区二区精品| 亚洲视频免费观看视频| 久久天躁狠狠躁夜夜2o2o | 下体分泌物呈黄色| 日韩一区二区视频免费看| 亚洲七黄色美女视频| 亚洲精品乱久久久久久| 欧美日韩精品网址| 免费在线观看完整版高清| 伊人亚洲综合成人网| 亚洲综合色网址| 欧美人与性动交α欧美软件| 欧美黑人精品巨大| 欧美亚洲 丝袜 人妻 在线| 免费av中文字幕在线| 亚洲一区中文字幕在线| 国产又色又爽无遮挡免| 亚洲精品,欧美精品| 精品一区在线观看国产| 免费黄色在线免费观看| 午夜福利,免费看| 又大又爽又粗| 日韩成人av中文字幕在线观看| 建设人人有责人人尽责人人享有的| 亚洲国产av新网站| 亚洲国产最新在线播放| 搡老岳熟女国产| 一级毛片黄色毛片免费观看视频| 亚洲精品国产一区二区精华液| 亚洲一区中文字幕在线| 国产av一区二区精品久久| 中文天堂在线官网| 欧美中文综合在线视频| 国产精品一二三区在线看| 久久天堂一区二区三区四区| 国产极品天堂在线| 久久精品久久精品一区二区三区| 涩涩av久久男人的天堂| 亚洲精品aⅴ在线观看| 中文精品一卡2卡3卡4更新| 又黄又粗又硬又大视频| 男人添女人高潮全过程视频| 亚洲精品久久久久久婷婷小说| 国产精品一区二区在线不卡| 韩国高清视频一区二区三区| 欧美最新免费一区二区三区| 久久久久视频综合| 国产亚洲av高清不卡| 精品人妻在线不人妻| 午夜日本视频在线| 91老司机精品| 国产乱人偷精品视频| 在线天堂最新版资源| 中文字幕色久视频| 亚洲欧美日韩另类电影网站| 国产探花极品一区二区| 99热网站在线观看| 一二三四在线观看免费中文在| 国产又色又爽无遮挡免| 精品国产乱码久久久久久小说| 91精品国产国语对白视频| 999久久久国产精品视频| 伦理电影大哥的女人| 人人妻人人澡人人看| av卡一久久| 久久精品国产a三级三级三级| 免费高清在线观看日韩| 可以免费在线观看a视频的电影网站 | 国产在线免费精品| 久久久久国产一级毛片高清牌| 最近手机中文字幕大全| 考比视频在线观看| 91aial.com中文字幕在线观看| av国产精品久久久久影院| 亚洲精品中文字幕在线视频| 老司机影院毛片| 女人精品久久久久毛片| 国产黄色视频一区二区在线观看| 亚洲在久久综合| 国产精品无大码| 日本爱情动作片www.在线观看| 在线亚洲精品国产二区图片欧美| 亚洲第一区二区三区不卡| 自拍欧美九色日韩亚洲蝌蚪91| 一本大道久久a久久精品| 亚洲第一av免费看| 只有这里有精品99| 中文字幕高清在线视频| 国产成人系列免费观看| 一区二区三区精品91| 精品第一国产精品| 黄色 视频免费看| 国产精品.久久久| 精品久久久久久电影网| 久热爱精品视频在线9| 国产精品二区激情视频| 亚洲国产av影院在线观看| 国产精品久久久久成人av| h视频一区二区三区| 亚洲av中文av极速乱| 精品第一国产精品| 日本91视频免费播放| 你懂的网址亚洲精品在线观看| 欧美日韩亚洲国产一区二区在线观看 | 丰满少妇做爰视频| 少妇被粗大的猛进出69影院| 亚洲熟女精品中文字幕| 亚洲四区av| 国产精品国产av在线观看| 亚洲美女黄色视频免费看| 男男h啪啪无遮挡| 日本欧美视频一区| 国产精品一国产av| 亚洲国产最新在线播放| 黄片小视频在线播放| 极品人妻少妇av视频| 久久韩国三级中文字幕| 国产精品久久久人人做人人爽| 亚洲人成电影观看| 一本久久精品| 亚洲人成网站在线观看播放| 国产成人免费观看mmmm| 麻豆乱淫一区二区| 人人妻人人澡人人爽人人夜夜| 日日撸夜夜添| 青春草国产在线视频| 看十八女毛片水多多多| 亚洲欧洲日产国产| 亚洲伊人色综图| 午夜福利网站1000一区二区三区| 黄色怎么调成土黄色| 亚洲av成人不卡在线观看播放网 | 老司机影院毛片| 国产亚洲欧美精品永久| 国产欧美日韩一区二区三区在线| 国产福利在线免费观看视频| 亚洲精品av麻豆狂野| 无遮挡黄片免费观看| 亚洲精品视频女| 看免费av毛片| 秋霞伦理黄片| 99热网站在线观看| 亚洲三区欧美一区| 免费女性裸体啪啪无遮挡网站| 欧美另类一区| 十八禁高潮呻吟视频| 日韩人妻精品一区2区三区| 久久97久久精品| 日本欧美视频一区| 一区在线观看完整版| 一级毛片我不卡| 两个人免费观看高清视频| 日本av免费视频播放| 精品久久久精品久久久| 亚洲久久久国产精品| 少妇人妻 视频| 欧美日韩av久久| 国产探花极品一区二区| 午夜久久久在线观看| 蜜桃国产av成人99| 交换朋友夫妻互换小说| 男人操女人黄网站| 在现免费观看毛片| 国产免费现黄频在线看| 久久ye,这里只有精品| 午夜精品国产一区二区电影| 色精品久久人妻99蜜桃| 99九九在线精品视频| 欧美黑人精品巨大| 中国三级夫妇交换| 飞空精品影院首页| 最新的欧美精品一区二区| 99久久99久久久精品蜜桃| 亚洲精品久久成人aⅴ小说| 91aial.com中文字幕在线观看| 大片免费播放器 马上看| 欧美在线黄色| 精品一区二区免费观看| 另类精品久久| 男的添女的下面高潮视频| 高清不卡的av网站| 久久久欧美国产精品| 日本wwww免费看| 97精品久久久久久久久久精品| 亚洲成色77777| av一本久久久久| 免费黄网站久久成人精品| 亚洲国产欧美一区二区综合| 丰满少妇做爰视频| 九草在线视频观看| 1024视频免费在线观看| 亚洲婷婷狠狠爱综合网| 日韩,欧美,国产一区二区三区| 又粗又硬又长又爽又黄的视频| 如何舔出高潮| 亚洲一区二区三区欧美精品| 久久久久网色| 国产成人精品在线电影| 欧美日韩视频高清一区二区三区二| 精品亚洲成a人片在线观看| 十八禁人妻一区二区| 亚洲精品国产av成人精品| 国产成人一区二区在线| 欧美中文综合在线视频| 国产精品麻豆人妻色哟哟久久| 大香蕉久久成人网| 欧美国产精品va在线观看不卡| 成年人免费黄色播放视频| 亚洲国产欧美在线一区| 97在线人人人人妻| av福利片在线| 建设人人有责人人尽责人人享有的| 欧美日韩一区二区视频在线观看视频在线| 亚洲,一卡二卡三卡| 日韩精品免费视频一区二区三区| 人人妻人人澡人人看| 一级黄片播放器| 丝袜美腿诱惑在线| 国产av国产精品国产| 久热这里只有精品99| 日本wwww免费看| 激情视频va一区二区三区| 精品视频人人做人人爽| 日韩大片免费观看网站| 国产无遮挡羞羞视频在线观看| 美女脱内裤让男人舔精品视频| 黄色怎么调成土黄色| 精品一品国产午夜福利视频| 中国三级夫妇交换| 国产乱人偷精品视频| 国产一区二区三区av在线| 亚洲欧美中文字幕日韩二区| 一级毛片黄色毛片免费观看视频| 王馨瑶露胸无遮挡在线观看| 看免费成人av毛片| www.熟女人妻精品国产| 午夜福利乱码中文字幕| 又粗又硬又长又爽又黄的视频| 久久人人97超碰香蕉20202| 亚洲精品自拍成人| 亚洲精品日本国产第一区| 久久精品人人爽人人爽视色| 国产99久久九九免费精品| 久久久精品免费免费高清| 日韩大片免费观看网站| 日日啪夜夜爽| 国产国语露脸激情在线看| 亚洲精品第二区| 国产有黄有色有爽视频| 高清不卡的av网站| 亚洲在久久综合| 免费av中文字幕在线| 国产成人系列免费观看| 亚洲国产av新网站| 久久人人97超碰香蕉20202| 亚洲精品,欧美精品| 中文乱码字字幕精品一区二区三区| 老司机影院毛片| xxxhd国产人妻xxx| 免费在线观看完整版高清| 亚洲精品国产一区二区精华液| 97在线人人人人妻| 99热国产这里只有精品6| 99热全是精品| 日韩av在线免费看完整版不卡| 亚洲欧美清纯卡通| 老司机影院毛片| 国产99久久九九免费精品| 1024香蕉在线观看| 国产精品 欧美亚洲| 天天影视国产精品| 亚洲一级一片aⅴ在线观看| 熟妇人妻不卡中文字幕| 久久久久精品久久久久真实原创| 久久精品久久久久久久性| 国产成人精品福利久久| 日韩成人av中文字幕在线观看| 成人午夜精彩视频在线观看| 免费久久久久久久精品成人欧美视频| 搡老岳熟女国产| 菩萨蛮人人尽说江南好唐韦庄| 晚上一个人看的免费电影| 这个男人来自地球电影免费观看 | 女人高潮潮喷娇喘18禁视频| 国产精品秋霞免费鲁丝片| 国产成人欧美| 欧美人与性动交α欧美精品济南到| 国产日韩欧美视频二区| 国产精品国产三级专区第一集| 欧美日韩精品网址| 国产精品 国内视频| 亚洲欧美一区二区三区久久| 国产av精品麻豆| av卡一久久| 19禁男女啪啪无遮挡网站| 91老司机精品| 国产探花极品一区二区| 看免费av毛片| 黄色一级大片看看| 少妇猛男粗大的猛烈进出视频| 国产探花极品一区二区| av在线播放精品| 国产精品二区激情视频| 亚洲精品久久午夜乱码| 91精品三级在线观看| 国产人伦9x9x在线观看| 久久天堂一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 考比视频在线观看| 中文字幕最新亚洲高清| 国产伦理片在线播放av一区| 亚洲情色 制服丝袜| 少妇被粗大的猛进出69影院| 中文字幕最新亚洲高清| 欧美黄色片欧美黄色片| 国产成人av激情在线播放| 久久av网站| 久久影院123| 亚洲视频免费观看视频| 欧美久久黑人一区二区| 一本一本久久a久久精品综合妖精| 国产高清国产精品国产三级| 久久av网站| 亚洲国产精品一区三区| 99re6热这里在线精品视频| 如何舔出高潮| 国产精品久久久久久人妻精品电影 | 秋霞伦理黄片| 一级片'在线观看视频| 久久国产精品男人的天堂亚洲| 七月丁香在线播放| 一级爰片在线观看| 免费久久久久久久精品成人欧美视频| 国产淫语在线视频| av一本久久久久| 国产精品国产av在线观看| 一级毛片电影观看| 久久亚洲国产成人精品v| 青春草亚洲视频在线观看| 亚洲人成网站在线观看播放| 亚洲av电影在线观看一区二区三区| 两个人免费观看高清视频| 国产精品欧美亚洲77777| 精品人妻在线不人妻| 欧美另类一区| 午夜福利,免费看| 黑人欧美特级aaaaaa片| 亚洲在久久综合| 看免费av毛片| 亚洲成色77777| 精品卡一卡二卡四卡免费| 日韩熟女老妇一区二区性免费视频| 一区二区三区四区激情视频| 亚洲中文av在线| 久久 成人 亚洲| 久久久久久久久久久免费av| 一本久久精品| 亚洲成色77777| 9191精品国产免费久久| 精品国产乱码久久久久久男人| 无限看片的www在线观看| 亚洲欧美精品自产自拍| 水蜜桃什么品种好| 久久性视频一级片| 欧美成人精品欧美一级黄| 成人黄色视频免费在线看| 悠悠久久av| 日本一区二区免费在线视频| 丁香六月天网| 日本猛色少妇xxxxx猛交久久| 成人午夜精彩视频在线观看| 久久综合国产亚洲精品| 日韩,欧美,国产一区二区三区| 一二三四中文在线观看免费高清| 人人妻人人澡人人看| 国产在线一区二区三区精| 十分钟在线观看高清视频www| 国产日韩欧美亚洲二区| 大片电影免费在线观看免费| 国产午夜精品一二区理论片| 日韩成人av中文字幕在线观看| 19禁男女啪啪无遮挡网站| 美女福利国产在线| 免费女性裸体啪啪无遮挡网站| 老司机亚洲免费影院| 国产成人精品无人区| 久久久久网色| 欧美亚洲 丝袜 人妻 在线| 亚洲成人一二三区av| 欧美人与性动交α欧美精品济南到| 亚洲精品国产av蜜桃| 亚洲av成人精品一二三区| av线在线观看网站| 色综合欧美亚洲国产小说| 久热这里只有精品99| 中国国产av一级| 不卡视频在线观看欧美| 99九九在线精品视频| 精品国产露脸久久av麻豆| 老司机深夜福利视频在线观看 | 99精国产麻豆久久婷婷| 久久久久网色| 又大又爽又粗| 纯流量卡能插随身wifi吗| 色婷婷av一区二区三区视频| 欧美97在线视频| 多毛熟女@视频| 午夜福利,免费看| 日日摸夜夜添夜夜爱| 最近中文字幕2019免费版| 最近手机中文字幕大全| 一级,二级,三级黄色视频| 妹子高潮喷水视频| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久av美女十八| 国产精品三级大全| 亚洲精品美女久久久久99蜜臀 | 精品人妻一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 两个人免费观看高清视频| 久久综合国产亚洲精品| 久久女婷五月综合色啪小说| av女优亚洲男人天堂| 丝袜人妻中文字幕| 在线亚洲精品国产二区图片欧美| 亚洲免费av在线视频| 美女午夜性视频免费| www.熟女人妻精品国产| 欧美激情极品国产一区二区三区| tube8黄色片| 国产男女内射视频| 欧美人与善性xxx| 男男h啪啪无遮挡| 51午夜福利影视在线观看| 最近的中文字幕免费完整| 一区在线观看完整版| 高清在线视频一区二区三区| √禁漫天堂资源中文www| 欧美黑人欧美精品刺激| 亚洲欧美一区二区三区久久| 亚洲国产欧美一区二区综合| 国产精品蜜桃在线观看| 成人国产av品久久久| 大码成人一级视频| 热99久久久久精品小说推荐| 久久国产亚洲av麻豆专区| 亚洲第一区二区三区不卡| 精品少妇一区二区三区视频日本电影 | 精品一区二区三区av网在线观看 | 国产99久久九九免费精品| 国产精品av久久久久免费| 一边摸一边抽搐一进一出视频| 最新的欧美精品一区二区| 午夜精品国产一区二区电影| 亚洲av福利一区| 国产伦理片在线播放av一区| 亚洲欧美一区二区三区久久| 色精品久久人妻99蜜桃| 国产爽快片一区二区三区| 在线观看三级黄色| 亚洲久久久国产精品| 欧美乱码精品一区二区三区| 亚洲精品自拍成人| 免费黄网站久久成人精品| 色综合欧美亚洲国产小说| 国产成人一区二区在线| 美国免费a级毛片| 亚洲一区二区三区欧美精品| 国产爽快片一区二区三区| 亚洲综合色网址| 欧美人与善性xxx| 一级毛片我不卡| 1024视频免费在线观看| 嫩草影视91久久| 午夜久久久在线观看| 亚洲国产日韩一区二区| 老司机影院成人| 香蕉国产在线看| 涩涩av久久男人的天堂| 男女国产视频网站| 波野结衣二区三区在线| 97人妻天天添夜夜摸| av福利片在线| 国产淫语在线视频| 国产av国产精品国产| 亚洲av日韩在线播放| 亚洲av综合色区一区| 女人精品久久久久毛片| 熟女av电影| 日韩欧美精品免费久久| 十八禁人妻一区二区| 国产乱人偷精品视频| av一本久久久久| 90打野战视频偷拍视频| 久久久久国产精品人妻一区二区| 一区二区av电影网| av视频免费观看在线观看| 久久韩国三级中文字幕| 日韩一本色道免费dvd| 免费黄网站久久成人精品| 国产激情久久老熟女| 岛国毛片在线播放| a 毛片基地| 日韩精品有码人妻一区| 秋霞伦理黄片| 精品少妇内射三级| 天天躁夜夜躁狠狠久久av| 国产欧美日韩一区二区三区在线| 免费久久久久久久精品成人欧美视频| 国产一区二区激情短视频 | 好男人视频免费观看在线| 日韩av在线免费看完整版不卡| 韩国精品一区二区三区| 欧美精品av麻豆av| 精品国产一区二区三区久久久樱花| 亚洲五月色婷婷综合| 国产成人91sexporn| 久久天堂一区二区三区四区| 日韩欧美精品免费久久| 国产精品成人在线| 搡老岳熟女国产| 国产99久久九九免费精品| 午夜福利网站1000一区二区三区| 久久久久久久国产电影| 黑人猛操日本美女一级片| 性高湖久久久久久久久免费观看| 一级a爱视频在线免费观看| 国产99久久九九免费精品| 亚洲一码二码三码区别大吗| 欧美成人精品欧美一级黄| 高清av免费在线| 国产精品.久久久| 久久久久国产精品人妻一区二区| 国产成人91sexporn| 久久av网站| 欧美国产精品va在线观看不卡| 婷婷色综合大香蕉| 日韩精品免费视频一区二区三区| 美女福利国产在线| 这个男人来自地球电影免费观看 | 国产精品一区二区在线不卡| 一级毛片黄色毛片免费观看视频| 美女视频免费永久观看网站| 中文精品一卡2卡3卡4更新| 黄色一级大片看看| www日本在线高清视频| 日韩伦理黄色片| 午夜福利一区二区在线看| 午夜91福利影院| 中文欧美无线码| 久久狼人影院| 男男h啪啪无遮挡| 免费高清在线观看视频在线观看| 欧美精品高潮呻吟av久久| 久久97久久精品| 麻豆av在线久日| 丝袜喷水一区| 日日摸夜夜添夜夜爱| 亚洲一区二区三区欧美精品| 免费黄频网站在线观看国产| av在线播放精品| 麻豆av在线久日| 免费黄频网站在线观看国产| 国产av国产精品国产| 最近最新中文字幕大全免费视频 | 亚洲人成77777在线视频| 久久精品久久久久久久性| 国产亚洲欧美精品永久| 午夜福利在线免费观看网站| 精品国产超薄肉色丝袜足j| 熟女少妇亚洲综合色aaa.| 夫妻性生交免费视频一级片| 妹子高潮喷水视频| 亚洲国产中文字幕在线视频| 精品国产一区二区三区久久久樱花| 亚洲精品中文字幕在线视频| 国产亚洲午夜精品一区二区久久| 99国产综合亚洲精品| 飞空精品影院首页| 国产片内射在线| 成年美女黄网站色视频大全免费| 天天躁夜夜躁狠狠躁躁| 成人漫画全彩无遮挡| 男女下面插进去视频免费观看| 亚洲av在线观看美女高潮| 午夜免费鲁丝| 日韩av在线免费看完整版不卡| 亚洲国产毛片av蜜桃av|