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

    索馬里東岸溫度鋒季節(jié)分布特征分析

    2017-12-27 02:16:34劉鵬張永剛曹震卿巨澤光
    海洋預(yù)報(bào) 2017年6期
    關(guān)鍵詞:上升流鋒面溫度梯度

    劉鵬,張永剛,曹震卿,巨澤光

    (1.海軍大連艦艇學(xué)院,遼寧大連116018;2.北海艦隊(duì)航空兵航管氣象處,山東青島266000)

    索馬里東岸溫度鋒季節(jié)分布特征分析

    劉鵬1,張永剛1,曹震卿1,巨澤光2

    (1.海軍大連艦艇學(xué)院,遼寧大連116018;2.北海艦隊(duì)航空兵航管氣象處,山東青島266000)

    基于50 a的SODA的溫度數(shù)據(jù),利用絕對(duì)梯度分析法求出鋒面的概率分布,對(duì)索馬里東岸的溫度鋒面的分布特征進(jìn)行了分析。此外利用WOA13的溫鹽數(shù)據(jù),對(duì)鋒區(qū)出現(xiàn)概率較大海區(qū)的聲速剖面進(jìn)行了分析。結(jié)果表明:溫度鋒面在淺水處的概率分布與深水處存在異同,鋒面主要集中分布在6—10月、6°~10°N沿岸附近,通過(guò)研究SODA海流數(shù)據(jù),發(fā)現(xiàn)鋒面分布概率較大的海域有著明顯的上升流。利用EOF分析,發(fā)現(xiàn)溫度鋒面的年際變化并不強(qiáng),主要是年內(nèi)的季節(jié)變化較為明顯。聲速剖面在有無(wú)鋒面存在的位置處存在明顯差異,上升流帶來(lái)的底層冷水導(dǎo)致聲速要明顯小于其他區(qū)域,鋒面同一位置處,由于夏季鋒強(qiáng)最大,鋒區(qū)的聲速也明顯小于其他季節(jié)。

    SODA;WOA13;溫度鋒;聲速

    1 引言

    海洋鋒一般是指特性明顯不同的兩種或幾種水體之間的狹窄過(guò)渡帶。狹義而言,有人將其定義為水團(tuán)之間的邊界線。廣義地說(shuō),可泛指任一種海洋環(huán)境參數(shù)的躍變帶,因而出現(xiàn)了諸如水溫鋒、鹽度鋒、密度鋒、聲速鋒、水色鋒、透明度鋒,以及海水化學(xué)、生物等等要素的海洋鋒的稱(chēng)謂[1]。

    索馬里半島位于非洲東北部,其東部沿海地區(qū)存在著隨季節(jié)變化的季風(fēng)及海流,因此該海域的水文環(huán)境異常復(fù)雜。隨著中國(guó)在中東地區(qū)的商業(yè)貿(mào)易以及護(hù)航任務(wù)的日趨頻繁,其戰(zhàn)略位置也愈發(fā)重要。目前國(guó)內(nèi)外對(duì)于這片海域的研究大多都集中在對(duì)這一海域的海溫年際變化、流系的三維結(jié)構(gòu)、海流對(duì)季風(fēng)的響應(yīng)等方面,如曹宗元[2]利用SODA數(shù)據(jù)對(duì)索馬里流系的三維結(jié)構(gòu)及其季節(jié)變化和年際變化進(jìn)行了分析,趙芳[3]利用Argo海溫資料和SODA資料對(duì)索馬里東岸的海溫時(shí)空變化進(jìn)行了研究,Sschott[4]對(duì)索馬里上升流與季風(fēng)的相應(yīng)做了研究分析。然而,對(duì)于索馬里東海岸的溫度鋒面研究較少,也很少研究鋒面對(duì)水聲傳播的影響。本文利用SODA數(shù)據(jù),分析了該海域50 a的溫度鋒面逐月平均分布特征,并利用WOA13數(shù)據(jù)分析了鋒區(qū)聲速剖面的分布。

    圖1 索馬里地圖

    2 數(shù)據(jù)及方法介紹

    2.1 數(shù)據(jù)介紹

    本文用來(lái)分析鋒面的數(shù)據(jù)是來(lái)自美國(guó)馬里蘭大學(xué)(http://dsrs.atmos.umd.edu/DATA/)的簡(jiǎn)單海洋同化數(shù)據(jù)(Simple Ocean Data Assimilation,SODA)的月平均數(shù)據(jù)資料。SODA海洋數(shù)據(jù)集包含的變量有溫度、鹽度、海流速度、海表風(fēng)應(yīng)力、海洋上層0~500 m熱含量、海洋上層0~125 m熱含量、海平面高度等。SODA數(shù)據(jù)的時(shí)間分布為1871年1月—2010年12月,水平分辨率為0.5°×0.5°,垂向上覆蓋了水深5.01~5 375 m不等的40層。本文所用的是1960年1月—2009年12月的溫度數(shù)據(jù),以及2000年1月—2004年12月的海流數(shù)據(jù)。

    用來(lái)分析聲速斷面的是來(lái)自美國(guó)國(guó)家海洋大氣局(National Oceanic and Atmospheric Administration,NOAA)的國(guó)家海洋數(shù)據(jù)中心海洋氣象實(shí)驗(yàn)室的海洋氣候態(tài)數(shù)據(jù)集產(chǎn)品WOA13(World Ocean Atlas13)數(shù)據(jù)(http://www.nodc.noaa.gov/OC5/indprod.html)。該數(shù)據(jù)包括了全球從1955—2012年的溫度、鹽度、溶解氧、磷酸鹽、硅酸鹽等海洋要素?cái)?shù)據(jù),時(shí)間上分為年平均、季節(jié)平均、月平均數(shù)據(jù),空間分辨率有5°、1°、1/4°等3種,在深度上,利用內(nèi)插值的方法,從表層到最大深度5 500 m分為了102層,其中0~100 m每5 m一層,100~500 m每25 m一層,500~2 000 m每50 m一層,2 000~5 500 m每100 m一層。本文選用的是1955—2012年、空間分辨率1/4°的季節(jié)平均的溫度和鹽度數(shù)據(jù)。

    2.2 方法介紹

    鋒面的提取,國(guó)內(nèi)外有許多人做過(guò)研究,主要是溫度梯度法,如劉建斌等[5-6]曾用絕對(duì)梯度分析法提取了印尼西岸的穩(wěn)鹽鋒面和阿爾沃蘭海鋒面,國(guó)外的Moore等[7]、Park等[8]、Wang等[9]利用溫度梯度法分析了不同海區(qū)的海洋鋒面并獲得較好的效果,本文同樣采取該方法提取了溫度鋒面,其中絕對(duì)梯度的定義為:

    式中:T為所要算的要素,本文為海水溫度。

    趙寶宏等[10]和朱鳳芹等[11]利用鋒面有效概率分布來(lái)分析鋒面,都能夠較好地判斷出鋒面隨時(shí)間的分布規(guī)律和主要分布區(qū)域。根據(jù)索馬里東海岸溫度梯度的平均分布,并考慮資料的差異性,本文選取0.015℃/km作為溫度鋒面的判斷標(biāo)準(zhǔn)。即某點(diǎn)的梯度大于判斷標(biāo)準(zhǔn),則把該點(diǎn)看做出現(xiàn)點(diǎn)。按照式(2)計(jì)算各點(diǎn)在50 a中的各月出現(xiàn)鋒的概率P,并選取概率大于0.5的區(qū)域作為鋒出現(xiàn)的區(qū)域。

    式中:N表示50 a中該點(diǎn)有效觀測(cè)數(shù)據(jù)次數(shù),n表示50 a中該點(diǎn)出現(xiàn)鋒的次數(shù)。

    2.3 聲速計(jì)算

    關(guān)于聲速計(jì)算公式目前已經(jīng)建立了幾十種,這些公式大多包含從幾項(xiàng)到幾十項(xiàng)不等的參數(shù),各公式對(duì)于海水溫度、鹽度和壓力的適用范圍也有所不同。本文基于WOA13溫鹽數(shù)據(jù),使用Mackenzie聲速經(jīng)驗(yàn)公式[12]計(jì)算了該海域的聲速剖面。

    3 分析與討論

    3.1 鋒面分布特征分析

    根據(jù)前文所介紹方法,作出在水深500 m以淺的逐月溫度梯度分布圖,圖2分別為1960年1—12月在水深96.92 m、171.4 m的溫度梯度分布圖。由圖可知:溫度鋒面在水深96.92 m處主要存在于6—10月,11月之后鋒面逐漸減弱,到次年1月鋒面消失,6月開(kāi)始鋒面又開(kāi)始逐漸增強(qiáng)。溫度鋒面存在的主要位置在索馬里沿岸地區(qū)6°~10°N、50°~53°E之間。其中6、7、8月鋒面分布范圍較廣,但是鋒面強(qiáng)度較??;9、10月鋒面集中在8°~10°N、51°~52°E沿岸海域,但是鋒面強(qiáng)度較大。在水深171.4 m處,溫度鋒面主要存在于7—12月,且出現(xiàn)一個(gè)環(huán)狀的鋒面,這與該區(qū)域存在的一個(gè)反氣旋型渦旋有關(guān)[2]。另外觀察其他更深水深層的溫度梯度可以發(fā)現(xiàn)從水深268.46 m開(kāi)始,溫度鋒面消失(圖略)。

    圖2 1960年1—12月不同水深溫度梯度分布

    分別作出水深15.07、57.98、96.92 m處的鋒面概率分布特征圖(見(jiàn)圖3—5)。圖中只給出了鋒面在50 a中出現(xiàn)概率大于50%,即50 a中能出現(xiàn)25次以上的區(qū)域,這些區(qū)域能夠很好地代表溫度鋒面存在的主要區(qū)域。

    圖3 1—12月水深15.07 m處出現(xiàn)概率大于50%的溫度鋒

    水深15.07 m處的溫度鋒50 a的概率分布如圖3所示,由圖可知:1—3月幾乎沒(méi)有鋒面存在,4月份開(kāi)始,從索馬里東北角12°~6°N的沿岸海域,逐漸出現(xiàn)穩(wěn)定的溫度鋒面,7月份鋒面范圍達(dá)到最大,且概率分布值大多達(dá)到90%以上,7月以后,鋒面逐漸減弱,到11月鋒面分布范圍逐漸趨近于0。總體來(lái)說(shuō),鋒面出現(xiàn)概率較大的區(qū)域在6—9月,8°~11°N沿岸海域。

    水深57.98 m處的溫度鋒50 a的概率分布如圖4所示,由圖可知:對(duì)比水深15.07 m處的鋒面概率分布圖,可以發(fā)現(xiàn)水深57.98 m處溫度鋒面的分布范圍明顯要更大,可見(jiàn)冷暖水團(tuán)在此處比淺層混合更劇烈。1—2月溫度鋒基本不存在,3月開(kāi)始,在10°N沿岸出現(xiàn)小范圍的鋒面,且概率分布值只在50%左右,4月以后鋒面分布逐漸向南部擴(kuò)散,6月份的時(shí)候在2°N沿岸附近也開(kāi)始出現(xiàn)了少量溫度鋒面,7月鋒面分布范圍達(dá)到最廣且在8°N沿岸附近概率分布值接近于100%,8月以后鋒面分布逐漸減小且在11°N、55°E出現(xiàn)一個(gè)概率較小的鋒面,9月2°N處鋒面消失,但10°N以北鋒面向東部延伸到55°E,10月索馬里東岸的鋒面主要集中在8°~12°N沿岸,11—12月鋒面逐漸減弱直至消失??傮w來(lái)說(shuō),鋒面出現(xiàn)概率較大的區(qū)域在6—9月,6°~10°N沿岸海域。

    圖4 1—12月水深57.98 m處出現(xiàn)概率大于50%的溫度鋒

    水深96.92 m處的溫度鋒50 a的概率分布如圖5所示,由圖可知:與淺水處不一樣的是,水深96.92 m處在1—3月也存在溫度鋒面,1月鋒面主要包括9°~12°N和3°~5°N兩處,2月鋒面范圍比1月小,只在10°~12°N和4°N附近存在兩個(gè)范圍很小的鋒面,3月南部的鋒面消失,只在8°~11°N處存在鋒面,4月鋒面范圍逐漸開(kāi)始擴(kuò)散,且分布概率值也開(kāi)始增大,4、5月鋒面主要存在于7°~12°N,6月鋒面向南移動(dòng)1°W緯度的距離,7月鋒面向東南擴(kuò)散且范圍增大,存在2°~10°N和2°N沿岸兩處鋒面,8月鋒面向東北擴(kuò)散,存在6°~12°N和2°N沿岸兩處鋒面,北部的鋒面向東擴(kuò)散到56°E附近,9月鋒面開(kāi)始收縮,集中在6°~12°N這一海域,10—12月鋒面分布逐漸減弱,集中在8°~12°N沿岸海域??傮w來(lái)說(shuō),鋒面出現(xiàn)概率較大的區(qū)域在6—10月,6°~10°N沿岸海域。

    圖5 1—12月水深96.92 m處出現(xiàn)概率大于50%的溫度鋒

    圖6 51.2°E、9.2°N處2000—2004年垂向流隨時(shí)間的變化(流向向上為正值,向下為負(fù)值)

    圖7 水深96.92 m處夏季(7—10月)溫度梯度的EOF

    圖8 四季聲速剖面圖(藍(lán)、紅、綠線分別代表(9°N、52°E)、(9°N、53°E)、(9°N、54°E)3處的聲速剖面)

    圖9 9°N、52°E四季聲速剖面圖

    總的說(shuō)來(lái),鋒面的概率分布圖與溫度梯度圖是比較吻合的,這也說(shuō)明了利用鋒面概率分布可以有效的預(yù)測(cè)出鋒面可能出現(xiàn)的位置和時(shí)間。此外根據(jù)各水深層溫度鋒面的分布特點(diǎn),可以發(fā)現(xiàn)鋒面主要集中分布在6—9月、7°~9°N沿岸海域,根據(jù)以往的資料顯示,該海域存在較強(qiáng)的上升流,如圖6所示為51.2°E、9.2°N處2000—2004年垂向流隨時(shí)間的變化,由圖可見(jiàn):上升流主要存在于6—9月、水深700 m以淺。在夏季的時(shí)候,上升流帶來(lái)的底層冷水與海表的高溫海水的混合,從而造成了兩個(gè)不同水團(tuán)的過(guò)渡帶,即形成了溫度鋒面。由溫度梯度圖可知鋒面主要存在區(qū)域在水深300 m以淺,可見(jiàn)并不是上升流所存在區(qū)域全都存在鋒面,而是上升流所帶來(lái)的冷水與表層暖水的過(guò)渡區(qū)域存在鋒面。

    3.2 EOF分析結(jié)果

    對(duì)鋒面較明顯的96.92 m水深夏季(7—9月)的50 a(1960—2009年)溫度梯度距平場(chǎng)做EOF分解(見(jiàn)圖7)。由圖可見(jiàn):前3個(gè)模態(tài)貢獻(xiàn)率之和有60.4%,第一模態(tài)在9°~10°N的沿岸海區(qū)的很小的一個(gè)區(qū)域內(nèi)存在一個(gè)正位相中心,時(shí)間序列的分布變化較為劇烈,其中時(shí)間系數(shù)為正值時(shí),代表正位相中心所在區(qū)域溫度梯度增大,時(shí)間系數(shù)為負(fù)值時(shí),正位相中心所在區(qū)域溫度梯度減小,總體上看第一模態(tài)的周期性并不是很明顯,但是可以看出在1976—1986年和1994—2005年這兩個(gè)時(shí)間段,溫度梯度總體呈現(xiàn)增大趨勢(shì)。第二模態(tài)在8°~9°N的沿岸海區(qū)存在一個(gè)正位相中心,時(shí)間序列的變化大致為一年正一年負(fù)的循環(huán)變化,變化幅度不是很大。第三模態(tài)在4°~8°N沿岸海區(qū)存在一個(gè)正位相中心,時(shí)間序列在1991—2009年呈上升趨勢(shì),可見(jiàn)該海區(qū)的溫度梯度在有逐漸增大的趨勢(shì)??傮w分析這3個(gè)模態(tài),可以知道各個(gè)模態(tài)分別對(duì)應(yīng)不同海區(qū)的溫度鋒面的時(shí)空變化特征,且總的來(lái)看索馬里東岸海區(qū)的溫度梯度變化的年際變化并不是很明顯。

    3.3 聲速剖面分析

    基于WOA13的溫鹽數(shù)據(jù),選擇鋒面上強(qiáng)度較大的點(diǎn)9°N、51°E,分析其所在位置聲速剖面,并與東部另外兩個(gè)鋒強(qiáng)較弱和無(wú)鋒面的兩個(gè)位置對(duì)比,由圖8可知:在鋒強(qiáng)較小的冬春秋季,這3處位置聲速剖面分布差別不大,主要差異在于表層聲速的差異,近岸處由于水溫較低,聲速較小。在夏季,上升流導(dǎo)致的溫度鋒面處聲速明顯異于無(wú)鋒面處,其聲速在鋒面所在水深0~200 m處要低于其他兩處,且一直延伸到水深1000 m處,聲速仍然最小。為了更好對(duì)比該點(diǎn)聲速剖面隨時(shí)間變化的差異,做出如圖9所示的9°N、52°E四季的聲速剖面,可以發(fā)現(xiàn)夏季的聲速在水深0~200 m處是明顯低于其他季節(jié)的,秋冬兩季表層聲速差異較小,夏季表層聲速與春季表層聲速相差了8 m/s,水深250 m附近四季聲速大小達(dá)到相同值,水深250 m以深,聲速差異較小。

    4 結(jié)論

    通過(guò)利用50 a的SODA溫度數(shù)據(jù)和WOA13溫鹽數(shù)據(jù),對(duì)索馬里東岸水深5~100 m處的溫度鋒的分布特征與鋒區(qū)聲速剖面進(jìn)行了分析,得出了以下結(jié)論:

    (1)溫度鋒面所在的區(qū)域,水深淺處的概率分布與較深處存在異同,且鋒面主要集中分布在6—10月,其他月份的鋒面的分布范圍及概率分布均較??;

    (2)鋒面主要集中分布在6°~10°N沿岸附近,且在6—10月鋒面分布概率較大的海域,其上升流也較為強(qiáng)烈,上升流帶來(lái)的底層冷水團(tuán)與海表的高溫水團(tuán)的混合,形成了溫度鋒面,可見(jiàn)上升流對(duì)該海域溫度鋒面有著極其重要的影響;

    (3)根據(jù)溫度梯度的EOF分析,空間場(chǎng)的前3個(gè)模態(tài)分別對(duì)應(yīng)溫度鋒面從北到南的3個(gè)不同海區(qū),時(shí)間場(chǎng)上,年際變化不是很明顯,第二模態(tài)的時(shí)間序列周期性較強(qiáng),呈現(xiàn)一年正一年負(fù)的循環(huán)。這也說(shuō)明此片海區(qū)的溫度鋒面的年際變化并不強(qiáng),主要的變化為年內(nèi)的季節(jié)變化;

    (4)聲速剖面在有無(wú)鋒面位置處存在明顯差異,上升流帶來(lái)的底層冷水導(dǎo)致聲速要明顯小于其他區(qū)域。鋒面同一位置處,在鋒強(qiáng)最大的夏季,鋒區(qū)的聲速也明顯小于其他季節(jié)。

    [1]馮士筰,李鳳歧,李少菁.海洋科學(xué)導(dǎo)論[M].北京:高等教育出版社,1999:470-479.

    [2]曹宗元.索馬里流系的三維結(jié)構(gòu)及其季節(jié)與年際變化[D].青島:中國(guó)海洋大學(xué),2013.

    [3]趙芳.索馬里流系區(qū)域海溫時(shí)空變化研究[D].青島:中國(guó)海洋大學(xué),2013.

    [4]Schott F.Monsoon Response of the Somali Current and Associated Upwelling[J].Progress in Oceanography,1983,12(3):357-381.

    [5]劉建斌,張永剛.印度尼西亞西岸溫鹽鋒時(shí)空變化特征[J].海洋預(yù)報(bào),2015,32(5):14-23.

    [6]劉建斌,張永剛.阿爾沃蘭海海洋鋒時(shí)空氣候變化特征[J].海洋預(yù)報(bào),2016,33(1):37-44.

    [7]Moore J K,Abbott M R,Richman J G.Variability in the Location of the Antarctic Polar Front(90°~20°W)from Satellite Sea Surface Temperature Data[J].Journal of Geophysics Research,1997,102(C13):27825-27833.

    [8]Park K A,Chung J Y,Kim K.Sea Surface Temperature Fronts in the East(Japan)Sea and Temporal Variations[J].Geophysical Research Letters,2004,31(7):L07304.

    [9]Wang D X,Liu Y,Qi Y Q,et al.Seasonal Variability of Thermal Fronts in the northern South China Sea from Satellite Data[J].Geophysical Research Letters,2001,28(20):3963-3966.

    [10]趙寶宏,劉宇迪,趙加華,等.南海海洋鋒季節(jié)分布特征初探[C]//2011年第二十八屆中國(guó)氣象學(xué)會(huì)年會(huì)論文集.廈門(mén):中國(guó)氣象學(xué)會(huì),2011.

    [11]朱鳳芹,謝玲玲,成印河.南海溫度鋒的分布特征及季節(jié)變化[J].海洋與湖沼,2014,45(4):695-702.

    [12]Mackenzie K V.Nine-Term Equation for Sound Speed in the Oceans[J].The Journal of the Acoustical Society of America,1981,70(3):807-812.

    Analysis on seasonal distribution of temperature front in Somalia East Coast

    LIU Peng1,ZHANG Yong-gang1,CAO Zhen-qing1,JU Ze-guang2
    (1.Dalian Naval Academy,Dalian 116018 China;2.North ocean fleet air traffic control meteorological office,Qingdao 266000 China)

    Based on the 50 years temperature data from SODA,the probability distribution of the temperature fronts was obtained by the absolute gradient analysis method.The characteristics of the temperature fronts at the east coast of Somalia were analyzed.Based on the temperature and salinity data of WOA13,sound velocity profiles were analyzed.The results showed that the probability distribution of the temperature fronts in the shallow water was a little different from that in the deep water,and the fronts were mainly distributed in the vicinity of 6—10 months and 6°~10°N.There was a significant increase in the area with large probability of frontal distribution flow.The sound velocity profile in the front area was obviously different from that of other areas.The cold water from the bottom of the upwelling current caused the sound velocity to be significantly smaller than the other areas.At the same position,the sound velocity of the front area is obviously smaller than that of other seasons.

    SODA;WOA13;Temperature fronts;Sound velocity

    P731.11

    A

    1003-0239(2017)06-0048-09

    10.11737/j.issn.1003-0239.2017.06.006

    2017-06-21;

    2017-07-17。

    劉鵬(1993-),男,碩士在讀,主要從事世界大洋海洋鋒面研究。E-mail:515540650@qq.com

    猜你喜歡
    上升流鋒面溫度梯度
    夏季舟山上升流長(zhǎng)期遙感觀測(cè)與分析
    溫度對(duì)中小跨徑混凝土梁橋應(yīng)力及變形的影響研究
    熱聲耦合燃燒振蕩中火焰鋒面識(shí)別分析
    煤氣與熱力(2022年2期)2022-03-09 06:29:20
    2019年夏季長(zhǎng)江口及鄰近海域鋒面控制下葉綠素a的分布特征及其環(huán)境影響因素分析
    遼東半島頂端海域上升流長(zhǎng)期變化特征及影響因素*
    海洋與湖沼(2020年1期)2020-02-06 06:01:38
    嚴(yán)寒地區(qū)混凝土箱梁實(shí)測(cè)溫度梯度分析
    山西建筑(2019年21期)2019-12-02 02:35:08
    基于核心素養(yǎng)的高中地理“問(wèn)題式教學(xué)”——以“鋒面氣旋”為例
    溫度梯度場(chǎng)對(duì)聲表面波器件影響研究
    電子制作(2018年23期)2018-12-26 01:01:20
    瓊東上升流的年際變化及長(zhǎng)期變化趨勢(shì)*
    海洋與湖沼(2016年1期)2016-01-15 03:50:49
    高速鐵路CRTSⅢ型板式無(wú)砟軌道溫度梯度試驗(yàn)研究
    鐵道建筑(2015年3期)2015-12-26 09:41:44
    老女人水多毛片| 亚洲综合色惰| 少妇丰满av| 2021少妇久久久久久久久久久| 国语对白做爰xxxⅹ性视频网站| 免费黄色在线免费观看| 久久女婷五月综合色啪小说| 日韩欧美 国产精品| 久久国内精品自在自线图片| 少妇高潮的动态图| 大片电影免费在线观看免费| 精品亚洲乱码少妇综合久久| 在线观看免费日韩欧美大片 | 欧美日本中文国产一区发布| 91久久精品国产一区二区三区| 国产亚洲av片在线观看秒播厂| 一级毛片 在线播放| 乱人伦中国视频| 精品久久久久久电影网| 亚洲精品一区蜜桃| av免费观看日本| 亚洲人与动物交配视频| 亚洲精品乱码久久久久久按摩| a 毛片基地| 午夜日本视频在线| 欧美xxxx性猛交bbbb| 这个男人来自地球电影免费观看 | 在线观看免费视频网站a站| 国产一区亚洲一区在线观看| 新久久久久国产一级毛片| 亚洲自偷自拍三级| 黄色欧美视频在线观看| 亚洲av.av天堂| 日韩电影二区| 久久久久人妻精品一区果冻| 伦精品一区二区三区| a级毛片免费高清观看在线播放| 久久久久人妻精品一区果冻| 午夜91福利影院| 人妻少妇偷人精品九色| 又粗又硬又长又爽又黄的视频| 九九爱精品视频在线观看| 国产免费一级a男人的天堂| 久久久久精品性色| 亚洲av国产av综合av卡| 婷婷色av中文字幕| 日本色播在线视频| 纯流量卡能插随身wifi吗| 美女国产视频在线观看| 涩涩av久久男人的天堂| 免费久久久久久久精品成人欧美视频 | 成人国产av品久久久| 最黄视频免费看| 精品亚洲成国产av| 成人特级av手机在线观看| 精品人妻熟女毛片av久久网站| 亚洲成人手机| 永久免费av网站大全| 成人毛片60女人毛片免费| 亚洲精品久久久久久婷婷小说| 久久国产精品大桥未久av | 爱豆传媒免费全集在线观看| 国产永久视频网站| 国产一区二区三区av在线| 99久久人妻综合| 有码 亚洲区| 成人免费观看视频高清| 亚洲欧美一区二区三区国产| 日本与韩国留学比较| 97在线人人人人妻| 男女国产视频网站| 久久久久国产网址| 好男人视频免费观看在线| 亚洲国产av新网站| 精品久久久久久电影网| av播播在线观看一区| 国产av码专区亚洲av| 久热这里只有精品99| 99久久人妻综合| 国产深夜福利视频在线观看| 丰满迷人的少妇在线观看| 日韩熟女老妇一区二区性免费视频| 97精品久久久久久久久久精品| 少妇人妻 视频| 日韩在线高清观看一区二区三区| 国产成人免费无遮挡视频| 亚洲精品乱码久久久久久按摩| 国产欧美日韩一区二区三区在线 | 亚洲精品,欧美精品| 久久久久久人妻| 国产亚洲最大av| 国产白丝娇喘喷水9色精品| 国产成人aa在线观看| 亚洲美女视频黄频| 人人妻人人看人人澡| 一级毛片我不卡| 日韩av不卡免费在线播放| 如何舔出高潮| 一级a做视频免费观看| 国产精品一区二区在线观看99| 亚洲精品第二区| 日韩中文字幕视频在线看片| 国产探花极品一区二区| 熟女人妻精品中文字幕| 日日撸夜夜添| 国产精品不卡视频一区二区| 大片免费播放器 马上看| 日本与韩国留学比较| 免费人成在线观看视频色| 91精品伊人久久大香线蕉| 极品人妻少妇av视频| 亚洲精品一二三| 欧美日韩国产mv在线观看视频| 国产真实伦视频高清在线观看| 日本91视频免费播放| 一区二区三区四区激情视频| 寂寞人妻少妇视频99o| 另类精品久久| 免费av中文字幕在线| 日本-黄色视频高清免费观看| 国产无遮挡羞羞视频在线观看| 日本猛色少妇xxxxx猛交久久| 黄色一级大片看看| 三级国产精品片| 一本—道久久a久久精品蜜桃钙片| 国产欧美亚洲国产| 亚洲成人一二三区av| 欧美另类一区| 国产免费福利视频在线观看| 亚洲国产欧美日韩在线播放 | 乱人伦中国视频| 国产成人免费无遮挡视频| 如日韩欧美国产精品一区二区三区 | 肉色欧美久久久久久久蜜桃| 蜜桃久久精品国产亚洲av| 青青草视频在线视频观看| 国产成人精品久久久久久| 欧美日韩亚洲高清精品| 久久影院123| 一级,二级,三级黄色视频| 大又大粗又爽又黄少妇毛片口| 六月丁香七月| 亚洲美女黄色视频免费看| 久久久欧美国产精品| 久久国内精品自在自线图片| 久久精品国产自在天天线| 免费看光身美女| 久久精品国产亚洲网站| 自拍欧美九色日韩亚洲蝌蚪91 | 成人18禁高潮啪啪吃奶动态图 | 免费看日本二区| 欧美性感艳星| 亚洲精品久久久久久婷婷小说| 80岁老熟妇乱子伦牲交| 亚洲人与动物交配视频| 夫妻午夜视频| 日韩av不卡免费在线播放| 男女免费视频国产| 精品久久久久久久久亚洲| 男女边摸边吃奶| 日韩欧美 国产精品| 国产在线男女| 久久久久久久久久久丰满| 超碰97精品在线观看| 高清不卡的av网站| 亚洲怡红院男人天堂| 国产白丝娇喘喷水9色精品| 97在线人人人人妻| av不卡在线播放| 日韩av免费高清视频| 亚洲精品乱码久久久久久按摩| 精品一品国产午夜福利视频| 国产精品99久久99久久久不卡 | av天堂中文字幕网| 成人国产av品久久久| 成人毛片a级毛片在线播放| 五月玫瑰六月丁香| 久久国内精品自在自线图片| 久久久久久久亚洲中文字幕| 岛国毛片在线播放| 亚洲欧美精品专区久久| 一边亲一边摸免费视频| 男人和女人高潮做爰伦理| 超碰97精品在线观看| 亚洲精品,欧美精品| 午夜福利网站1000一区二区三区| 少妇人妻一区二区三区视频| 亚洲国产精品一区三区| 久久99热6这里只有精品| 中文资源天堂在线| 中文字幕人妻丝袜制服| xxx大片免费视频| 人妻系列 视频| 色5月婷婷丁香| 精品久久久久久久久亚洲| 中文字幕av电影在线播放| 精品亚洲乱码少妇综合久久| 久久久久久久久久人人人人人人| 蜜臀久久99精品久久宅男| 亚洲内射少妇av| 老司机亚洲免费影院| 亚洲国产毛片av蜜桃av| av.在线天堂| 亚洲av男天堂| 少妇猛男粗大的猛烈进出视频| 女的被弄到高潮叫床怎么办| 成人午夜精彩视频在线观看| 少妇的逼水好多| 纵有疾风起免费观看全集完整版| 中文欧美无线码| 国内少妇人妻偷人精品xxx网站| 女人精品久久久久毛片| av在线播放精品| 一本久久精品| 高清毛片免费看| 亚洲在久久综合| 只有这里有精品99| 夫妻午夜视频| 麻豆精品久久久久久蜜桃| 免费人成在线观看视频色| 两个人的视频大全免费| 日韩强制内射视频| 亚洲欧美中文字幕日韩二区| 午夜免费男女啪啪视频观看| 自拍偷自拍亚洲精品老妇| 伦理电影免费视频| 99热这里只有是精品在线观看| 日韩一区二区视频免费看| 日本黄大片高清| 久久人人爽人人爽人人片va| 日韩 亚洲 欧美在线| 99re6热这里在线精品视频| 观看av在线不卡| 18禁在线无遮挡免费观看视频| 热re99久久精品国产66热6| 国产精品一区二区三区四区免费观看| 国内揄拍国产精品人妻在线| 婷婷色综合www| 免费人成在线观看视频色| 欧美日韩视频精品一区| 亚洲婷婷狠狠爱综合网| 国产欧美日韩一区二区三区在线 | 伊人久久国产一区二区| 精品久久久噜噜| 黄色一级大片看看| 亚洲国产av新网站| 欧美日韩国产mv在线观看视频| 久久久久久久亚洲中文字幕| 在线天堂最新版资源| 一级片'在线观看视频| 黄色视频在线播放观看不卡| 精品酒店卫生间| 大片电影免费在线观看免费| 久久久a久久爽久久v久久| 亚洲真实伦在线观看| 亚洲精品国产av成人精品| 欧美高清成人免费视频www| av免费在线看不卡| 嘟嘟电影网在线观看| 亚洲电影在线观看av| 国精品久久久久久国模美| 国产午夜精品一二区理论片| 黑人高潮一二区| 色视频www国产| 国产成人精品无人区| 卡戴珊不雅视频在线播放| 日韩电影二区| 精品午夜福利在线看| 婷婷色av中文字幕| 欧美精品一区二区大全| 中文字幕av电影在线播放| 国产又色又爽无遮挡免| 啦啦啦在线观看免费高清www| 狠狠精品人妻久久久久久综合| 亚洲av中文av极速乱| 精品卡一卡二卡四卡免费| 欧美人与善性xxx| 丰满人妻一区二区三区视频av| 久久久久久久大尺度免费视频| 亚洲一区二区三区欧美精品| 亚洲国产精品成人久久小说| 少妇高潮的动态图| av国产精品久久久久影院| 久久狼人影院| 啦啦啦在线观看免费高清www| 99久久人妻综合| 免费人妻精品一区二区三区视频| 最近2019中文字幕mv第一页| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦中文免费视频观看日本| 午夜免费观看性视频| 一本色道久久久久久精品综合| 精品99又大又爽又粗少妇毛片| 国内少妇人妻偷人精品xxx网站| 午夜福利影视在线免费观看| 老司机影院毛片| 久久久久人妻精品一区果冻| 国产成人午夜福利电影在线观看| 精品午夜福利在线看| 亚洲国产精品国产精品| 久久精品国产自在天天线| 麻豆精品久久久久久蜜桃| 国产免费福利视频在线观看| 午夜老司机福利剧场| 亚洲国产日韩一区二区| 美女xxoo啪啪120秒动态图| 亚洲三级黄色毛片| 午夜激情久久久久久久| 国产成人精品一,二区| 欧美bdsm另类| 中文字幕免费在线视频6| 国产午夜精品久久久久久一区二区三区| 日本欧美国产在线视频| 夫妻性生交免费视频一级片| 日日摸夜夜添夜夜爱| 97超碰精品成人国产| 在线观看三级黄色| 国产一区有黄有色的免费视频| 少妇人妻 视频| 香蕉精品网在线| 亚洲欧美精品专区久久| 秋霞在线观看毛片| 91久久精品电影网| 亚洲欧美中文字幕日韩二区| 亚洲欧美一区二区三区黑人 | 久久免费观看电影| 51国产日韩欧美| 国产精品一区二区性色av| 中文字幕精品免费在线观看视频 | 日韩中字成人| 欧美日韩亚洲高清精品| 欧美bdsm另类| 亚洲不卡免费看| 欧美日韩精品成人综合77777| 亚洲精品成人av观看孕妇| 国产成人精品无人区| 亚洲四区av| 日韩熟女老妇一区二区性免费视频| 亚洲精品久久午夜乱码| 亚洲欧美日韩东京热| 天堂8中文在线网| 尾随美女入室| 青春草亚洲视频在线观看| 中文字幕久久专区| 久久人人爽人人爽人人片va| 日本欧美视频一区| 777米奇影视久久| 国产一区有黄有色的免费视频| 七月丁香在线播放| 又爽又黄a免费视频| 国产伦精品一区二区三区视频9| 18禁动态无遮挡网站| 国产精品成人在线| 美女视频免费永久观看网站| 69精品国产乱码久久久| 免费黄色在线免费观看| 交换朋友夫妻互换小说| 亚洲人与动物交配视频| 国产精品国产三级专区第一集| 国产片特级美女逼逼视频| 插逼视频在线观看| 一级毛片我不卡| 看十八女毛片水多多多| 国产免费一区二区三区四区乱码| 全区人妻精品视频| 热re99久久国产66热| 少妇高潮的动态图| 国模一区二区三区四区视频| 99精国产麻豆久久婷婷| 国内精品宾馆在线| 大话2 男鬼变身卡| 五月天丁香电影| 人妻夜夜爽99麻豆av| 成年人午夜在线观看视频| 我的老师免费观看完整版| 99久久精品国产国产毛片| 久久女婷五月综合色啪小说| 亚洲精品国产av成人精品| 国产高清不卡午夜福利| 精品卡一卡二卡四卡免费| 久久国产精品大桥未久av | 久久久亚洲精品成人影院| av不卡在线播放| 久久99一区二区三区| 老司机亚洲免费影院| 纵有疾风起免费观看全集完整版| 日韩 亚洲 欧美在线| 亚洲精品乱码久久久久久按摩| 乱系列少妇在线播放| 人妻夜夜爽99麻豆av| 丝袜脚勾引网站| 国内少妇人妻偷人精品xxx网站| 亚洲欧美一区二区三区国产| 亚洲精品日韩在线中文字幕| 狂野欧美白嫩少妇大欣赏| 日韩一本色道免费dvd| 水蜜桃什么品种好| 亚洲精品aⅴ在线观看| 春色校园在线视频观看| 亚洲欧美一区二区三区黑人 | 哪个播放器可以免费观看大片| 丰满饥渴人妻一区二区三| 国产一区二区三区av在线| 妹子高潮喷水视频| 视频区图区小说| av线在线观看网站| 久久久久久伊人网av| 99久久人妻综合| 国产精品秋霞免费鲁丝片| av黄色大香蕉| av又黄又爽大尺度在线免费看| 美女中出高潮动态图| 午夜久久久在线观看| 少妇的逼水好多| 美女xxoo啪啪120秒动态图| a级片在线免费高清观看视频| 亚洲情色 制服丝袜| 看免费成人av毛片| 精品少妇内射三级| 精品熟女少妇av免费看| 在现免费观看毛片| 午夜日本视频在线| 欧美丝袜亚洲另类| 最黄视频免费看| 久久久久久久久久久丰满| 亚洲国产精品专区欧美| 婷婷色av中文字幕| 国内少妇人妻偷人精品xxx网站| 老司机影院成人| 汤姆久久久久久久影院中文字幕| 美女cb高潮喷水在线观看| 日韩电影二区| 婷婷色综合大香蕉| 少妇熟女欧美另类| 伦理电影大哥的女人| 嫩草影院新地址| 免费人成在线观看视频色| 亚洲欧美成人综合另类久久久| 国产精品国产三级国产专区5o| 三级国产精品欧美在线观看| 91久久精品国产一区二区成人| 国产伦理片在线播放av一区| 大码成人一级视频| 国产乱人偷精品视频| 这个男人来自地球电影免费观看 | 欧美精品国产亚洲| 在线观看人妻少妇| 国产伦精品一区二区三区视频9| 如何舔出高潮| 黄片无遮挡物在线观看| 免费av不卡在线播放| 18禁在线无遮挡免费观看视频| 黑丝袜美女国产一区| 女性被躁到高潮视频| 亚洲一区二区三区欧美精品| 午夜精品国产一区二区电影| 久久久久视频综合| 秋霞伦理黄片| 噜噜噜噜噜久久久久久91| 人妻制服诱惑在线中文字幕| 精品久久久精品久久久| 国产永久视频网站| 国产精品蜜桃在线观看| 国产成人91sexporn| 欧美精品高潮呻吟av久久| 色哟哟·www| 亚洲久久久国产精品| 一区二区三区免费毛片| 精品人妻偷拍中文字幕| 中文字幕制服av| 日本vs欧美在线观看视频 | 日韩成人av中文字幕在线观看| 日本黄色日本黄色录像| 高清欧美精品videossex| 亚洲精品自拍成人| 亚洲av综合色区一区| 日本wwww免费看| 一个人免费看片子| 久久影院123| 亚洲成色77777| 美女cb高潮喷水在线观看| 男人舔奶头视频| 亚洲国产欧美日韩在线播放 | 国产精品不卡视频一区二区| 一区二区三区四区激情视频| 欧美日韩综合久久久久久| 2021少妇久久久久久久久久久| 22中文网久久字幕| 丰满饥渴人妻一区二区三| 哪个播放器可以免费观看大片| 国产深夜福利视频在线观看| 国产精品三级大全| 国产欧美亚洲国产| 久久久久网色| 免费人妻精品一区二区三区视频| 十八禁高潮呻吟视频 | 国产在线免费精品| 一本一本综合久久| av黄色大香蕉| 噜噜噜噜噜久久久久久91| 成人亚洲精品一区在线观看| 水蜜桃什么品种好| 国产真实伦视频高清在线观看| 蜜桃久久精品国产亚洲av| 国产亚洲欧美精品永久| 精品少妇黑人巨大在线播放| 久久国内精品自在自线图片| 日日啪夜夜撸| 99热这里只有精品一区| 国内揄拍国产精品人妻在线| 精品亚洲乱码少妇综合久久| 人人妻人人澡人人看| 免费黄网站久久成人精品| 如日韩欧美国产精品一区二区三区 | 国产 一区精品| 国产精品久久久久久精品电影小说| 国产熟女欧美一区二区| 亚洲欧美日韩卡通动漫| 欧美一级a爱片免费观看看| 国产亚洲午夜精品一区二区久久| 在线播放无遮挡| 久久97久久精品| 九草在线视频观看| 高清不卡的av网站| 国产在线免费精品| av国产精品久久久久影院| 女人久久www免费人成看片| 免费观看无遮挡的男女| 黄色日韩在线| 青春草亚洲视频在线观看| 夜夜爽夜夜爽视频| 丰满人妻一区二区三区视频av| tube8黄色片| 国产亚洲欧美精品永久| 国产欧美日韩精品一区二区| 欧美少妇被猛烈插入视频| 日韩在线高清观看一区二区三区| 男女边吃奶边做爰视频| 日韩精品免费视频一区二区三区 | 成年av动漫网址| 日本午夜av视频| 国产精品久久久久久精品古装| 美女主播在线视频| 国产爽快片一区二区三区| 亚洲怡红院男人天堂| 国产深夜福利视频在线观看| 99热网站在线观看| 老熟女久久久| 国产精品国产三级国产av玫瑰| 久久久久久久亚洲中文字幕| 欧美性感艳星| 永久网站在线| 亚洲国产成人一精品久久久| 啦啦啦在线观看免费高清www| 成年人免费黄色播放视频 | 插逼视频在线观看| 久久人人爽人人片av| 最新中文字幕久久久久| 午夜免费观看性视频| 六月丁香七月| 男女边摸边吃奶| 中文字幕人妻熟人妻熟丝袜美| 欧美性感艳星| 三级国产精品欧美在线观看| 在线亚洲精品国产二区图片欧美 | 韩国av在线不卡| 亚洲怡红院男人天堂| 国产伦精品一区二区三区四那| 18禁在线无遮挡免费观看视频| a级毛片在线看网站| 国产乱人偷精品视频| 伦理电影免费视频| 一级毛片黄色毛片免费观看视频| 一级av片app| 日本猛色少妇xxxxx猛交久久| 黄色一级大片看看| 久久久久国产精品人妻一区二区| 七月丁香在线播放| 日韩成人伦理影院| 一级毛片我不卡| 亚洲精品日韩在线中文字幕| 亚洲精品一区蜜桃| 一级毛片电影观看| 全区人妻精品视频| 特大巨黑吊av在线直播| 久久久精品94久久精品| 波野结衣二区三区在线| 久久精品久久久久久久性| 国产精品国产三级国产av玫瑰| 国产伦在线观看视频一区| 亚洲精品日本国产第一区| 欧美bdsm另类| 97超碰精品成人国产| 一区二区三区乱码不卡18| 日产精品乱码卡一卡2卡三| 亚洲精品一区蜜桃| 久久午夜综合久久蜜桃| 寂寞人妻少妇视频99o| 丝袜在线中文字幕| 色哟哟·www| 国产有黄有色有爽视频| 亚洲精品中文字幕在线视频 | 又黄又爽又刺激的免费视频.| 亚洲色图综合在线观看| 日韩人妻高清精品专区| 免费黄频网站在线观看国产| 高清不卡的av网站| 国产亚洲5aaaaa淫片| 精品亚洲成a人片在线观看| av国产精品久久久久影院| 亚洲av中文av极速乱| 乱码一卡2卡4卡精品| 高清不卡的av网站| 免费高清在线观看视频在线观看| 一级爰片在线观看|