• <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
    91九色精品人成在线观看| 麻豆成人av在线观看| 嫩草影院精品99| 99热只有精品国产| 岛国视频午夜一区免费看| 日韩欧美一区二区三区在线观看| 美女大奶头视频| av中文乱码字幕在线| 久久人人精品亚洲av| 99re在线观看精品视频| 99re在线观看精品视频| 国产精品香港三级国产av潘金莲| 黄色女人牲交| 久久国产乱子伦精品免费另类| 欧美黑人精品巨大| 国产精品久久久人人做人人爽| 国产成年人精品一区二区 | 性色av乱码一区二区三区2| 久久久久久大精品| 亚洲av电影在线进入| 亚洲精品国产精品久久久不卡| 99国产综合亚洲精品| 交换朋友夫妻互换小说| 在线观看午夜福利视频| 国产黄色免费在线视频| 嫁个100分男人电影在线观看| 久久午夜亚洲精品久久| 国产欧美日韩综合在线一区二区| 日韩高清综合在线| 在线看a的网站| 1024视频免费在线观看| 99riav亚洲国产免费| 麻豆国产av国片精品| 国产色视频综合| 真人一进一出gif抽搐免费| bbb黄色大片| tocl精华| 日韩三级视频一区二区三区| av中文乱码字幕在线| 岛国在线观看网站| 最新在线观看一区二区三区| 国产免费现黄频在线看| 精品人妻1区二区| 亚洲欧美精品综合一区二区三区| 黄色怎么调成土黄色| 国产精品美女特级片免费视频播放器 | 夜夜看夜夜爽夜夜摸 | 大香蕉久久成人网| 91精品三级在线观看| 午夜精品国产一区二区电影| 日韩大尺度精品在线看网址 | 亚洲精品一二三| 涩涩av久久男人的天堂| 国产深夜福利视频在线观看| 巨乳人妻的诱惑在线观看| 乱人伦中国视频| 欧美不卡视频在线免费观看 | 欧美中文综合在线视频| 亚洲第一欧美日韩一区二区三区| 国产97色在线日韩免费| 欧美另类亚洲清纯唯美| √禁漫天堂资源中文www| 黄色丝袜av网址大全| 丰满人妻熟妇乱又伦精品不卡| 级片在线观看| xxxhd国产人妻xxx| 精品高清国产在线一区| 不卡av一区二区三区| 亚洲欧美精品综合一区二区三区| 国产精品免费视频内射| 国产亚洲av高清不卡| 女人被躁到高潮嗷嗷叫费观| 99久久人妻综合| 中国美女看黄片| 级片在线观看| 91老司机精品| 人妻丰满熟妇av一区二区三区| 在线观看免费高清a一片| 亚洲情色 制服丝袜| 国产区一区二久久| 中文字幕色久视频| 欧美黄色淫秽网站| 日本三级黄在线观看| 日韩大码丰满熟妇| 日韩三级视频一区二区三区| 嫩草影院精品99| 香蕉丝袜av| 男女之事视频高清在线观看| 天堂中文最新版在线下载| 色精品久久人妻99蜜桃| 午夜精品在线福利| 波多野结衣一区麻豆| 99热国产这里只有精品6| 精品久久蜜臀av无| 久久伊人香网站| 91老司机精品| 日日摸夜夜添夜夜添小说| 亚洲专区字幕在线| 又黄又爽又免费观看的视频| 亚洲精品久久午夜乱码| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久久5区| 男女下面进入的视频免费午夜 | 99久久精品国产亚洲精品| 90打野战视频偷拍视频| 男人舔女人下体高潮全视频| 午夜激情av网站| 亚洲欧洲精品一区二区精品久久久| 三级毛片av免费| 欧美激情高清一区二区三区| 一二三四社区在线视频社区8| 国产精品久久久av美女十八| 麻豆久久精品国产亚洲av | 99精品久久久久人妻精品| 精品日产1卡2卡| 日日爽夜夜爽网站| 亚洲av五月六月丁香网| 男人舔女人的私密视频| 无遮挡黄片免费观看| avwww免费| 国产亚洲精品久久久久5区| 久9热在线精品视频| 午夜久久久在线观看| 黄色怎么调成土黄色| 国产免费男女视频| 身体一侧抽搐| 亚洲激情在线av| 亚洲性夜色夜夜综合| 久久午夜亚洲精品久久| 精品日产1卡2卡| 中文字幕av电影在线播放| 亚洲午夜精品一区,二区,三区| 纯流量卡能插随身wifi吗| 欧美黑人欧美精品刺激| 亚洲一码二码三码区别大吗| 欧美日韩视频精品一区| 亚洲成国产人片在线观看| 精品一区二区三区av网在线观看| 久久精品亚洲熟妇少妇任你| 久久婷婷成人综合色麻豆| 男女下面进入的视频免费午夜 | 中出人妻视频一区二区| 99精品在免费线老司机午夜| 成年版毛片免费区| 韩国精品一区二区三区| 两个人免费观看高清视频| 99热只有精品国产| 久久久久久久精品吃奶| 日本wwww免费看| 精品欧美一区二区三区在线| 国产精品一区二区免费欧美| 亚洲欧美一区二区三区久久| 亚洲成人免费电影在线观看| 巨乳人妻的诱惑在线观看| 五月开心婷婷网| 亚洲人成电影观看| 校园春色视频在线观看| 亚洲欧美一区二区三区久久| 黄色 视频免费看| 国产又色又爽无遮挡免费看| 亚洲狠狠婷婷综合久久图片| 波多野结衣一区麻豆| 亚洲精品国产色婷婷电影| 国产精品乱码一区二三区的特点 | 女同久久另类99精品国产91| 成人av一区二区三区在线看| 91麻豆精品激情在线观看国产 | 黑人操中国人逼视频| 老汉色∧v一级毛片| 91麻豆精品激情在线观看国产 | 美女福利国产在线| 热re99久久国产66热| 青草久久国产| 午夜久久久在线观看| 亚洲欧美一区二区三区久久| 亚洲 欧美 日韩 在线 免费| 国产av一区二区精品久久| 俄罗斯特黄特色一大片| 青草久久国产| xxx96com| 亚洲三区欧美一区| 午夜福利欧美成人| 久久中文字幕人妻熟女| 欧美日韩福利视频一区二区| 欧美中文综合在线视频| 久久精品亚洲精品国产色婷小说| 色哟哟哟哟哟哟| 黄色女人牲交| 大型av网站在线播放| 精品一区二区三区四区五区乱码| 国产精品一区二区免费欧美| 欧美精品亚洲一区二区| 亚洲精品一区av在线观看| 亚洲精品国产精品久久久不卡| 中文欧美无线码| 国产精品久久电影中文字幕| 亚洲欧美一区二区三区久久| 久久人人97超碰香蕉20202| 在线免费观看的www视频| 久久精品国产亚洲av高清一级| 久99久视频精品免费| 黑人猛操日本美女一级片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人一区二区三区免费视频网站| 午夜成年电影在线免费观看| 中文字幕高清在线视频| 悠悠久久av| 日本精品一区二区三区蜜桃| 国产黄色免费在线视频| 夫妻午夜视频| 老司机在亚洲福利影院| 国产无遮挡羞羞视频在线观看| 老熟妇仑乱视频hdxx| 午夜福利,免费看| 满18在线观看网站| 亚洲欧美激情综合另类| 制服人妻中文乱码| 一级黄色大片毛片| 国产成人系列免费观看| 男人舔女人下体高潮全视频| 黄色视频不卡| 中文欧美无线码| 国内久久婷婷六月综合欲色啪| 757午夜福利合集在线观看| 一边摸一边抽搐一进一出视频| 国产真人三级小视频在线观看| 午夜福利在线免费观看网站| 别揉我奶头~嗯~啊~动态视频| 国产av一区在线观看免费| 亚洲va日本ⅴa欧美va伊人久久| 嫩草影视91久久| 国产一区二区三区综合在线观看| 国产av又大| 久久伊人香网站| 九色亚洲精品在线播放| 欧美精品一区二区免费开放| 亚洲 欧美一区二区三区| 天天影视国产精品| 午夜日韩欧美国产| 色婷婷av一区二区三区视频| 久久中文字幕人妻熟女| 久久久久久久久中文| 精品熟女少妇八av免费久了| 国产一卡二卡三卡精品| 曰老女人黄片| 亚洲欧美精品综合久久99| 少妇裸体淫交视频免费看高清 | 午夜日韩欧美国产| 热99re8久久精品国产| 天天添夜夜摸| 乱人伦中国视频| 欧美人与性动交α欧美软件| 精品久久蜜臀av无| 欧美日韩中文字幕国产精品一区二区三区 | 久久亚洲真实| 欧美午夜高清在线| 国产不卡一卡二| 午夜视频精品福利| 亚洲一区二区三区色噜噜 | 精品一区二区三区av网在线观看| 国产aⅴ精品一区二区三区波| 极品人妻少妇av视频| 亚洲一区二区三区欧美精品| 一夜夜www| 国产精品二区激情视频| 欧美在线黄色| 亚洲av成人不卡在线观看播放网| 国产伦一二天堂av在线观看| 国产高清videossex| 亚洲成a人片在线一区二区| 在线观看日韩欧美| 1024视频免费在线观看| 嫩草影院精品99| 成在线人永久免费视频| 亚洲熟妇熟女久久| 另类亚洲欧美激情| 色在线成人网| 老熟妇仑乱视频hdxx| 亚洲专区国产一区二区| 岛国在线观看网站| 国产亚洲精品一区二区www| 亚洲va日本ⅴa欧美va伊人久久| 久久 成人 亚洲| 99精品在免费线老司机午夜| 精品久久久久久成人av| 99久久国产精品久久久| 99riav亚洲国产免费| 男女之事视频高清在线观看| 亚洲av电影在线进入| 日韩三级视频一区二区三区| 女生性感内裤真人,穿戴方法视频| 久久久久久久午夜电影 | 精品福利永久在线观看| 乱人伦中国视频| 精品一区二区三区视频在线观看免费 | 精品福利永久在线观看| 日韩视频一区二区在线观看| 女性生殖器流出的白浆| 国产精品免费一区二区三区在线| 久久精品aⅴ一区二区三区四区| 波多野结衣av一区二区av| 激情视频va一区二区三区| 曰老女人黄片| 中文字幕精品免费在线观看视频| 色综合欧美亚洲国产小说| 别揉我奶头~嗯~啊~动态视频| 国产精品秋霞免费鲁丝片| 国产成人av教育| 国产aⅴ精品一区二区三区波| 一级毛片精品| 国产精品98久久久久久宅男小说| 亚洲成人免费av在线播放| 在线免费观看的www视频| 免费女性裸体啪啪无遮挡网站| 国产精品av久久久久免费| 中亚洲国语对白在线视频| 在线观看日韩欧美| 黑人猛操日本美女一级片| 黑人操中国人逼视频| 久久香蕉精品热| 亚洲 国产 在线| 日韩欧美一区视频在线观看| 一个人免费在线观看的高清视频| 国产欧美日韩一区二区精品| 黑人欧美特级aaaaaa片| 香蕉久久夜色| 欧美成人午夜精品| 日本五十路高清| 深夜精品福利| 18美女黄网站色大片免费观看| 熟女少妇亚洲综合色aaa.| 啦啦啦在线免费观看视频4| 亚洲av五月六月丁香网| 在线视频色国产色| 啦啦啦 在线观看视频| 大香蕉久久成人网| 久久热在线av| 女同久久另类99精品国产91| 日韩成人在线观看一区二区三区| 长腿黑丝高跟| 午夜精品在线福利| 精品久久久久久久久久免费视频 | 国产欧美日韩一区二区三区在线| 99国产精品免费福利视频| 这个男人来自地球电影免费观看| 91麻豆av在线| 黄色毛片三级朝国网站| 欧美日韩瑟瑟在线播放| 亚洲av成人不卡在线观看播放网| 在线永久观看黄色视频| 欧美日韩乱码在线| 18禁黄网站禁片午夜丰满| 丰满迷人的少妇在线观看| 亚洲一区二区三区不卡视频| 97人妻天天添夜夜摸| 亚洲专区国产一区二区| 91九色精品人成在线观看| 免费观看精品视频网站| 国产成人av激情在线播放| 国产精品久久电影中文字幕| 国产成人欧美在线观看| 精品熟女少妇八av免费久了| 这个男人来自地球电影免费观看| 午夜福利免费观看在线| 国产一区二区在线av高清观看| 88av欧美| 亚洲视频免费观看视频| 在线永久观看黄色视频| av欧美777| 丁香欧美五月| 涩涩av久久男人的天堂| 欧美日韩亚洲高清精品| 老汉色∧v一级毛片| 曰老女人黄片| 女性被躁到高潮视频| 在线观看免费视频日本深夜| 精品一品国产午夜福利视频| 一区在线观看完整版| 欧美日韩国产mv在线观看视频| 99re在线观看精品视频| 国产精华一区二区三区| 亚洲av成人一区二区三| 在线观看免费视频网站a站| 18禁黄网站禁片午夜丰满| 又黄又粗又硬又大视频| 欧美午夜高清在线| 99久久国产精品久久久| 免费在线观看视频国产中文字幕亚洲| 亚洲视频免费观看视频| 人人妻人人澡人人看| 国产成+人综合+亚洲专区| 少妇被粗大的猛进出69影院| 日韩高清综合在线| 国产精品电影一区二区三区| 日韩 欧美 亚洲 中文字幕| 欧美精品一区二区免费开放| 欧美乱码精品一区二区三区| 亚洲,欧美精品.| 在线av久久热| 久久久久久免费高清国产稀缺| 日韩欧美免费精品| 黄色丝袜av网址大全| 99久久久亚洲精品蜜臀av| 黄频高清免费视频| 黑人猛操日本美女一级片| 国产精品永久免费网站| 国产精品国产高清国产av| 一进一出好大好爽视频| 免费在线观看黄色视频的| 亚洲国产毛片av蜜桃av| 日韩人妻精品一区2区三区| 美女福利国产在线| 两个人免费观看高清视频| 99久久国产精品久久久| 亚洲av五月六月丁香网| 亚洲熟妇熟女久久| 欧洲精品卡2卡3卡4卡5卡区| ponron亚洲| 亚洲黑人精品在线| 国产1区2区3区精品| 另类亚洲欧美激情| 亚洲欧美激情综合另类| 亚洲成人免费电影在线观看| 黑丝袜美女国产一区| 亚洲精品国产色婷婷电影| 亚洲全国av大片| 中亚洲国语对白在线视频| 国产亚洲精品一区二区www| 国产成年人精品一区二区 | 欧美亚洲日本最大视频资源| а√天堂www在线а√下载| 日韩大尺度精品在线看网址 | 黄色女人牲交| e午夜精品久久久久久久| 亚洲熟女毛片儿| 久久久久久久精品吃奶| 亚洲中文日韩欧美视频| 亚洲激情在线av| 精品欧美一区二区三区在线| 黄色 视频免费看| 久久天躁狠狠躁夜夜2o2o| 91在线观看av| 另类亚洲欧美激情| 中出人妻视频一区二区| 亚洲一区二区三区色噜噜 | 国产成人欧美在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久视频播放| av电影中文网址| 亚洲熟妇中文字幕五十中出 | 国产亚洲欧美98| 日本wwww免费看| 91九色精品人成在线观看| 在线视频色国产色| а√天堂www在线а√下载| 久久久国产一区二区| 欧美激情久久久久久爽电影 | 午夜福利一区二区在线看| 精品无人区乱码1区二区| 日韩中文字幕欧美一区二区| 国内毛片毛片毛片毛片毛片| 怎么达到女性高潮| 在线永久观看黄色视频| 成年女人毛片免费观看观看9| 91字幕亚洲| 久久精品亚洲av国产电影网| 99国产精品一区二区三区| 淫秽高清视频在线观看| 日韩大尺度精品在线看网址 | 黄片小视频在线播放| 欧美成人免费av一区二区三区| 成人国产一区最新在线观看| 丝袜人妻中文字幕| 欧美性长视频在线观看| 亚洲熟妇熟女久久| 久9热在线精品视频| 在线播放国产精品三级| 日日夜夜操网爽| 欧美乱色亚洲激情| 97人妻天天添夜夜摸| 97人妻天天添夜夜摸| 欧美性长视频在线观看| 热99re8久久精品国产| 精品久久蜜臀av无| 热re99久久国产66热| 久久久久久久午夜电影 | 久久国产精品男人的天堂亚洲| 久久香蕉激情| 国产亚洲精品综合一区在线观看 | 精品福利永久在线观看| 欧美乱色亚洲激情| av在线天堂中文字幕 | 国产精品影院久久| 男人操女人黄网站| 久久精品91蜜桃| 黄片大片在线免费观看| 少妇 在线观看| av网站免费在线观看视频| 不卡av一区二区三区| 女性生殖器流出的白浆| 人妻丰满熟妇av一区二区三区| 精品国产美女av久久久久小说| 十分钟在线观看高清视频www| 一进一出抽搐gif免费好疼 | 热99国产精品久久久久久7| 女生性感内裤真人,穿戴方法视频| 亚洲情色 制服丝袜| 另类亚洲欧美激情| 午夜久久久在线观看| 男女下面进入的视频免费午夜 | 麻豆久久精品国产亚洲av | 久久精品国产综合久久久| 中文字幕人妻丝袜一区二区| 国产激情久久老熟女| 免费少妇av软件| 成年女人毛片免费观看观看9| 村上凉子中文字幕在线| 日韩一卡2卡3卡4卡2021年| 嫩草影视91久久| 免费在线观看黄色视频的| 99久久久亚洲精品蜜臀av| 日韩有码中文字幕| 日韩免费av在线播放| 久久午夜亚洲精品久久| 欧美日韩av久久| 日韩视频一区二区在线观看| 亚洲色图综合在线观看| 亚洲久久久国产精品| 国产高清videossex| 日韩 欧美 亚洲 中文字幕| 狂野欧美激情性xxxx| 久久精品91无色码中文字幕| 99久久99久久久精品蜜桃| 亚洲精品国产精品久久久不卡| 怎么达到女性高潮| 人人妻人人爽人人添夜夜欢视频| 91老司机精品| 亚洲熟妇中文字幕五十中出 | 亚洲男人天堂网一区| 91大片在线观看| 亚洲av美国av| 高清在线国产一区| 在线永久观看黄色视频| 精品第一国产精品| 91麻豆精品激情在线观看国产 | 欧美av亚洲av综合av国产av| 国产欧美日韩综合在线一区二区| 国产高清视频在线播放一区| 久9热在线精品视频| 12—13女人毛片做爰片一| 别揉我奶头~嗯~啊~动态视频| 每晚都被弄得嗷嗷叫到高潮| 日韩有码中文字幕| 欧美日韩亚洲高清精品| 久久精品成人免费网站| 在线观看66精品国产| 亚洲五月天丁香| 色老头精品视频在线观看| 狠狠狠狠99中文字幕| 久久精品亚洲熟妇少妇任你| 亚洲精品中文字幕一二三四区| 欧美性长视频在线观看| 日日干狠狠操夜夜爽| www日本在线高清视频| 亚洲精品中文字幕一二三四区| 超碰成人久久| 大型av网站在线播放| 在线十欧美十亚洲十日本专区| 交换朋友夫妻互换小说| 少妇 在线观看| 如日韩欧美国产精品一区二区三区| 老司机靠b影院| 18禁黄网站禁片午夜丰满| 91成人精品电影| 亚洲欧美精品综合久久99| 婷婷六月久久综合丁香| 91大片在线观看| 亚洲精品国产精品久久久不卡| 91精品国产国语对白视频| 天天添夜夜摸| 老司机深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 在线十欧美十亚洲十日本专区| cao死你这个sao货| 侵犯人妻中文字幕一二三四区| a在线观看视频网站| 国产一区二区三区视频了| 99热国产这里只有精品6| 成人国产一区最新在线观看| 操美女的视频在线观看| 久久人人精品亚洲av| 校园春色视频在线观看| 黄片小视频在线播放| 可以在线观看毛片的网站| 成在线人永久免费视频| 女性被躁到高潮视频| 男人的好看免费观看在线视频 | 最好的美女福利视频网| 国产区一区二久久| 国产又爽黄色视频| 久久久国产成人精品二区 | 久久久久久亚洲精品国产蜜桃av| av在线播放免费不卡| 黄片播放在线免费| av有码第一页| 日韩国内少妇激情av| 人成视频在线观看免费观看| 精品国产亚洲在线| 99在线视频只有这里精品首页| 亚洲精品在线观看二区| 操美女的视频在线观看| 一a级毛片在线观看| 成人18禁高潮啪啪吃奶动态图| 一二三四社区在线视频社区8| 韩国av一区二区三区四区|