• <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
    99久久无色码亚洲精品果冻| 999久久久精品免费观看国产| 大型黄色视频在线免费观看| 日本熟妇午夜| 亚洲国产高清在线一区二区三| 精华霜和精华液先用哪个| 黄色视频,在线免费观看| 久久久久九九精品影院| 成人三级做爰电影| 手机成人av网站| 亚洲美女黄片视频| 制服丝袜大香蕉在线| 久久这里只有精品19| 真实男女啪啪啪动态图| 欧美日韩综合久久久久久 | tocl精华| 中文在线观看免费www的网站| 亚洲最大成人中文| 一a级毛片在线观看| 性色av乱码一区二区三区2| 18禁黄网站禁片午夜丰满| 熟女少妇亚洲综合色aaa.| 熟妇人妻久久中文字幕3abv| 天天躁狠狠躁夜夜躁狠狠躁| 两个人视频免费观看高清| 97超级碰碰碰精品色视频在线观看| 真人做人爱边吃奶动态| 18禁观看日本| 长腿黑丝高跟| 国产成人一区二区三区免费视频网站| 日韩欧美在线乱码| 日韩欧美精品v在线| 叶爱在线成人免费视频播放| xxxwww97欧美| 亚洲欧美日韩东京热| 久久人人精品亚洲av| 亚洲无线在线观看| 欧美3d第一页| 亚洲专区字幕在线| 精品久久蜜臀av无| 网址你懂的国产日韩在线| 国产一区二区在线观看日韩 | 免费看十八禁软件| 日本一本二区三区精品| 成人av一区二区三区在线看| 亚洲精品粉嫩美女一区| 精品无人区乱码1区二区| 麻豆一二三区av精品| 国产精品综合久久久久久久免费| 噜噜噜噜噜久久久久久91| 久久这里只有精品19| 色综合婷婷激情| 大型黄色视频在线免费观看| 少妇裸体淫交视频免费看高清| 色视频www国产| 亚洲国产日韩欧美精品在线观看 | 黄色视频,在线免费观看| 免费在线观看亚洲国产| 欧美精品啪啪一区二区三区| 日韩欧美一区二区三区在线观看| 国产亚洲av嫩草精品影院| 中文在线观看免费www的网站| 日本成人三级电影网站| 十八禁网站免费在线| 真人一进一出gif抽搐免费| 亚洲欧美日韩高清在线视频| 99久久久亚洲精品蜜臀av| 哪里可以看免费的av片| 99久久国产精品久久久| 天天一区二区日本电影三级| 在线永久观看黄色视频| 久久久色成人| 亚洲午夜理论影院| 在线观看一区二区三区| 又大又爽又粗| 午夜精品一区二区三区免费看| 亚洲一区二区三区色噜噜| 日韩欧美一区二区三区在线观看| 麻豆国产av国片精品| 日日夜夜操网爽| 18禁观看日本| 日日干狠狠操夜夜爽| 一卡2卡三卡四卡精品乱码亚洲| 又黄又爽又免费观看的视频| 黄片小视频在线播放| 毛片女人毛片| 在线国产一区二区在线| 亚洲自拍偷在线| 日日干狠狠操夜夜爽| 国产精品久久电影中文字幕| 亚洲中文字幕一区二区三区有码在线看 | 90打野战视频偷拍视频| 特大巨黑吊av在线直播| 成年版毛片免费区| 免费av毛片视频| 亚洲精品粉嫩美女一区| 久久久国产精品麻豆| 中文字幕人妻丝袜一区二区| 黄频高清免费视频| 欧美成人免费av一区二区三区| 老司机在亚洲福利影院| 老汉色av国产亚洲站长工具| 国产精品亚洲av一区麻豆| 成人高潮视频无遮挡免费网站| 国产不卡一卡二| 18美女黄网站色大片免费观看| 观看免费一级毛片| 18禁观看日本| 99久久无色码亚洲精品果冻| 怎么达到女性高潮| 亚洲成人精品中文字幕电影| 老鸭窝网址在线观看| 黑人欧美特级aaaaaa片| 日本五十路高清| 一个人免费在线观看电影 | 天堂动漫精品| 欧美黄色淫秽网站| 国产精品亚洲一级av第二区| 成年女人永久免费观看视频| 午夜福利成人在线免费观看| 我的老师免费观看完整版| 琪琪午夜伦伦电影理论片6080| 欧美色欧美亚洲另类二区| 亚洲最大成人中文| 老汉色av国产亚洲站长工具| 日本 欧美在线| 亚洲精品456在线播放app | 欧美日本视频| 亚洲男人的天堂狠狠| 一个人免费在线观看的高清视频| 夜夜夜夜夜久久久久| 欧美日韩瑟瑟在线播放| АⅤ资源中文在线天堂| 婷婷丁香在线五月| 久久精品91蜜桃| 亚洲精品国产精品久久久不卡| 日本撒尿小便嘘嘘汇集6| 在线十欧美十亚洲十日本专区| 热99在线观看视频| 亚洲成人久久爱视频| 久久精品亚洲精品国产色婷小说| 国内揄拍国产精品人妻在线| 黄色片一级片一级黄色片| 又黄又粗又硬又大视频| 亚洲18禁久久av| 在线免费观看不下载黄p国产 | 欧美黑人欧美精品刺激| 桃色一区二区三区在线观看| 黄色日韩在线| 九九热线精品视视频播放| 精品久久蜜臀av无| 国产一级毛片七仙女欲春2| 成人特级av手机在线观看| 久久精品国产99精品国产亚洲性色| av天堂中文字幕网| 成在线人永久免费视频| 久久伊人香网站| 三级国产精品欧美在线观看 | 欧美日韩精品网址| 亚洲国产欧洲综合997久久,| 欧美日韩亚洲国产一区二区在线观看| 成人av一区二区三区在线看| 老鸭窝网址在线观看| 脱女人内裤的视频| 无遮挡黄片免费观看| 亚洲精品久久国产高清桃花| 亚洲一区二区三区不卡视频| 中文字幕人妻丝袜一区二区| 美女扒开内裤让男人捅视频| 此物有八面人人有两片| 俺也久久电影网| 久久久久久久久中文| 99re在线观看精品视频| 黄色片一级片一级黄色片| 国内精品一区二区在线观看| 久久精品综合一区二区三区| 欧美中文综合在线视频| av天堂在线播放| 夜夜爽天天搞| 欧美日韩国产亚洲二区| 久久天堂一区二区三区四区| 性欧美人与动物交配| 中文亚洲av片在线观看爽| 噜噜噜噜噜久久久久久91| 欧洲精品卡2卡3卡4卡5卡区| 国产av在哪里看| 欧美成人一区二区免费高清观看 | 日日干狠狠操夜夜爽| 久久久久久久久免费视频了| 日本与韩国留学比较| 18美女黄网站色大片免费观看| 国产伦精品一区二区三区四那| 精华霜和精华液先用哪个| 毛片女人毛片| 老司机深夜福利视频在线观看| 91麻豆精品激情在线观看国产| 国产淫片久久久久久久久 | 午夜福利成人在线免费观看| 久久久久亚洲av毛片大全| 嫩草影院入口| 国产精品99久久久久久久久| 亚洲精品国产精品久久久不卡| 人妻久久中文字幕网| 国产高潮美女av| 男女之事视频高清在线观看| 99久久精品一区二区三区| www日本黄色视频网| 国产成人影院久久av| 欧美日本视频| 国产午夜精品久久久久久| 精品日产1卡2卡| 国产一区二区三区视频了| 在线观看美女被高潮喷水网站 | 国产欧美日韩一区二区精品| 亚洲欧美精品综合一区二区三区| 男女之事视频高清在线观看| 午夜日韩欧美国产| 亚洲av成人av| 丰满人妻熟妇乱又伦精品不卡| 欧美+亚洲+日韩+国产| 午夜福利在线观看吧| netflix在线观看网站| 久久天堂一区二区三区四区| 久久草成人影院| 三级国产精品欧美在线观看 | 三级国产精品欧美在线观看 | 成人18禁在线播放| 99国产综合亚洲精品| 精品一区二区三区视频在线 | 非洲黑人性xxxx精品又粗又长| 丁香六月欧美| 麻豆久久精品国产亚洲av| 中文字幕人妻丝袜一区二区| 男人和女人高潮做爰伦理| 99精品在免费线老司机午夜| 99视频精品全部免费 在线 | 欧美大码av| 特级一级黄色大片| 三级男女做爰猛烈吃奶摸视频| 日韩欧美一区二区三区在线观看| 桃色一区二区三区在线观看| 午夜福利在线观看免费完整高清在 | 久久欧美精品欧美久久欧美| 男人和女人高潮做爰伦理| 亚洲国产精品合色在线| 日韩欧美免费精品| 小蜜桃在线观看免费完整版高清| 国产精品久久久人人做人人爽| 热99re8久久精品国产| 久久国产精品影院| 99久久无色码亚洲精品果冻| 亚洲最大成人中文| 日韩欧美免费精品| 国产精品野战在线观看| 天堂动漫精品| 午夜福利欧美成人| 欧美色欧美亚洲另类二区| 亚洲精品在线美女| 久久久久国内视频| 网址你懂的国产日韩在线| 听说在线观看完整版免费高清| svipshipincom国产片| 午夜免费成人在线视频| 操出白浆在线播放| 熟女电影av网| 老鸭窝网址在线观看| 国产午夜精品论理片| 91麻豆精品激情在线观看国产| 久久久国产成人精品二区| 一区二区三区高清视频在线| 免费观看的影片在线观看| 成人国产综合亚洲| 日日干狠狠操夜夜爽| 一区二区三区国产精品乱码| 变态另类丝袜制服| 夜夜躁狠狠躁天天躁| 男人舔女人下体高潮全视频| 精品久久蜜臀av无| 午夜福利18| 一级作爱视频免费观看| 国产99白浆流出| 国产精品 国内视频| 熟妇人妻久久中文字幕3abv| 俄罗斯特黄特色一大片| 91在线观看av| 露出奶头的视频| 午夜精品一区二区三区免费看| 国产精品av视频在线免费观看| 亚洲成a人片在线一区二区| 一区二区三区高清视频在线| 最近在线观看免费完整版| 精品不卡国产一区二区三区| 中亚洲国语对白在线视频| 老司机深夜福利视频在线观看| 久久久久久久午夜电影| 国产伦精品一区二区三区四那| 88av欧美| 真实男女啪啪啪动态图| 黄片小视频在线播放| 法律面前人人平等表现在哪些方面| 亚洲乱码一区二区免费版| 在线观看免费午夜福利视频| 免费在线观看影片大全网站| 久久久色成人| 人妻夜夜爽99麻豆av| 免费看十八禁软件| 在线观看免费视频日本深夜| av视频在线观看入口| 国产亚洲精品av在线| 19禁男女啪啪无遮挡网站| 最新美女视频免费是黄的| 亚洲美女视频黄频| 国产精品永久免费网站| 久久婷婷人人爽人人干人人爱| 欧美一级毛片孕妇| 巨乳人妻的诱惑在线观看| 免费在线观看成人毛片| 男人舔奶头视频| 精品一区二区三区视频在线 | 黄片小视频在线播放| 国产久久久一区二区三区| 一级a爱片免费观看的视频| 99久久精品国产亚洲精品| 午夜福利高清视频| 操出白浆在线播放| 国产三级中文精品| 美女 人体艺术 gogo| 人妻夜夜爽99麻豆av| 给我免费播放毛片高清在线观看| 日本黄色视频三级网站网址| 这个男人来自地球电影免费观看| 国产毛片a区久久久久| 日本免费a在线| 久久精品国产综合久久久| 国产一区二区激情短视频| 夜夜躁狠狠躁天天躁| 在线免费观看的www视频| 色尼玛亚洲综合影院| 少妇丰满av| 99久国产av精品| 欧美av亚洲av综合av国产av| 精品午夜福利视频在线观看一区| 18禁黄网站禁片免费观看直播| 色视频www国产| 99国产极品粉嫩在线观看| 国产精品美女特级片免费视频播放器 | 亚洲七黄色美女视频| 欧美乱色亚洲激情| 亚洲av第一区精品v没综合| 成人永久免费在线观看视频| 国产黄色小视频在线观看| 好看av亚洲va欧美ⅴa在| 成在线人永久免费视频| 国产免费男女视频| 岛国在线免费视频观看| 国产黄片美女视频| 91久久精品国产一区二区成人 | 动漫黄色视频在线观看| 香蕉国产在线看| 999久久久精品免费观看国产| 日本一二三区视频观看| 亚洲一区二区三区色噜噜| 两性夫妻黄色片| 成人特级黄色片久久久久久久| www.www免费av| 国产视频一区二区在线看| 成人无遮挡网站| 小说图片视频综合网站| 久久九九热精品免费| 一区二区三区激情视频| 精品一区二区三区视频在线 | 视频区欧美日本亚洲| 国产1区2区3区精品| 国产一区二区在线av高清观看| 国产探花在线观看一区二区| 亚洲精品国产精品久久久不卡| 国产欧美日韩精品一区二区| 久久久久久大精品| 亚洲专区字幕在线| 国产精品久久久久久人妻精品电影| 亚洲色图av天堂| 母亲3免费完整高清在线观看| 亚洲午夜精品一区,二区,三区| 亚洲av熟女| 手机成人av网站| 嫩草影院入口| 动漫黄色视频在线观看| 熟女人妻精品中文字幕| 免费看a级黄色片| 欧美丝袜亚洲另类 | 天天躁日日操中文字幕| 国产综合懂色| 欧美又色又爽又黄视频| 国产成人影院久久av| 午夜两性在线视频| 男人和女人高潮做爰伦理| 人妻夜夜爽99麻豆av| 在线观看日韩欧美| 最好的美女福利视频网| 欧美成人免费av一区二区三区| 色噜噜av男人的天堂激情| 国产成人啪精品午夜网站| 久久婷婷人人爽人人干人人爱| 高清在线国产一区| 十八禁网站免费在线| 日日摸夜夜添夜夜添小说| 欧美日韩中文字幕国产精品一区二区三区| 久久久久免费精品人妻一区二区| 少妇裸体淫交视频免费看高清| 成人高潮视频无遮挡免费网站| 免费人成视频x8x8入口观看| 欧美日韩综合久久久久久 | 国产精品1区2区在线观看.| 18美女黄网站色大片免费观看| 熟女电影av网| 欧美日韩一级在线毛片| 最好的美女福利视频网| 国产精品久久电影中文字幕| 亚洲精品色激情综合| 成人一区二区视频在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲国产日韩欧美精品在线观看 | 熟女少妇亚洲综合色aaa.| 在线永久观看黄色视频| 中文字幕久久专区| 听说在线观看完整版免费高清| 两性午夜刺激爽爽歪歪视频在线观看| 中亚洲国语对白在线视频| 精品国产亚洲在线| 狂野欧美白嫩少妇大欣赏| 亚洲真实伦在线观看| 亚洲精品456在线播放app | 真人做人爱边吃奶动态| 国产黄色小视频在线观看| 免费无遮挡裸体视频| 18禁裸乳无遮挡免费网站照片| 999久久久精品免费观看国产| 日本a在线网址| 亚洲av五月六月丁香网| 日本一二三区视频观看| 嫩草影视91久久| 又爽又黄无遮挡网站| 久久精品夜夜夜夜夜久久蜜豆| 天堂影院成人在线观看| 欧洲精品卡2卡3卡4卡5卡区| 男女做爰动态图高潮gif福利片| av黄色大香蕉| 国产一区二区激情短视频| 韩国av一区二区三区四区| av视频在线观看入口| 黑人巨大精品欧美一区二区mp4| 久久午夜亚洲精品久久| 神马国产精品三级电影在线观看| 色av中文字幕| 亚洲av美国av| 国产亚洲精品av在线| 69av精品久久久久久| 99视频精品全部免费 在线 | avwww免费| 两个人看的免费小视频| 久久久国产精品麻豆| 亚洲va日本ⅴa欧美va伊人久久| 亚洲五月天丁香| av在线蜜桃| 老司机午夜十八禁免费视频| 一a级毛片在线观看| 最近最新中文字幕大全电影3| 国产精品一区二区三区四区久久| 在线观看舔阴道视频| 国产成人精品无人区| 九九在线视频观看精品| 熟女人妻精品中文字幕| www日本黄色视频网| 国产aⅴ精品一区二区三区波| 亚洲中文av在线| 日本熟妇午夜| 两性夫妻黄色片| 午夜免费成人在线视频| 日本在线视频免费播放| 女生性感内裤真人,穿戴方法视频| 成人永久免费在线观看视频| www.www免费av| 男人舔女人的私密视频| 国产v大片淫在线免费观看| 色吧在线观看| 午夜日韩欧美国产| 最近最新免费中文字幕在线| 老熟妇仑乱视频hdxx| 久久久国产成人精品二区| 欧美午夜高清在线| 亚洲中文日韩欧美视频| 男女视频在线观看网站免费| 欧美一区二区精品小视频在线| 国产单亲对白刺激| 欧美黄色片欧美黄色片| netflix在线观看网站| 欧美黑人巨大hd| 色哟哟哟哟哟哟| 日本撒尿小便嘘嘘汇集6| 国产精品av视频在线免费观看| 精品国内亚洲2022精品成人| 51午夜福利影视在线观看| 午夜两性在线视频| xxxwww97欧美| 1024香蕉在线观看| 日韩欧美在线乱码| 18美女黄网站色大片免费观看| 大型黄色视频在线免费观看| 法律面前人人平等表现在哪些方面| 天堂√8在线中文| 特级一级黄色大片| 国产亚洲精品久久久com| 欧美又色又爽又黄视频| 国产av一区在线观看免费| 视频区欧美日本亚洲| 国产乱人伦免费视频| 欧美日韩福利视频一区二区| 国产一区二区在线观看日韩 | 免费在线观看视频国产中文字幕亚洲| 国产亚洲av嫩草精品影院| 久久人人精品亚洲av| 国产又黄又爽又无遮挡在线| 久久中文字幕人妻熟女| 亚洲熟女毛片儿| 在线免费观看不下载黄p国产 | 99久久精品一区二区三区| 亚洲av熟女| 成人av一区二区三区在线看| 欧美极品一区二区三区四区| 亚洲熟妇熟女久久| 欧美不卡视频在线免费观看| 一级作爱视频免费观看| 久久精品夜夜夜夜夜久久蜜豆| 日日摸夜夜添夜夜添小说| 91久久精品国产一区二区成人 | 99re在线观看精品视频| 欧美最黄视频在线播放免费| 欧美zozozo另类| 12—13女人毛片做爰片一| 成熟少妇高潮喷水视频| 成年女人永久免费观看视频| 亚洲精品一卡2卡三卡4卡5卡| bbb黄色大片| 久久精品aⅴ一区二区三区四区| 波多野结衣巨乳人妻| 国产精品 国内视频| 亚洲av成人不卡在线观看播放网| 国产高清videossex| 黄色视频,在线免费观看| 午夜两性在线视频| 性欧美人与动物交配| 国产精品免费一区二区三区在线| 欧美日本视频| 午夜免费激情av| 免费在线观看影片大全网站| 精品久久久久久成人av| 欧美日韩中文字幕国产精品一区二区三区| 国产乱人伦免费视频| 午夜激情福利司机影院| 欧美午夜高清在线| 91九色精品人成在线观看| 午夜免费激情av| 成人三级黄色视频| 夜夜躁狠狠躁天天躁| 欧美zozozo另类| 悠悠久久av| 99热这里只有精品一区 | 久久人妻av系列| 免费大片18禁| 黄色 视频免费看| 成人三级做爰电影| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av五月六月丁香网| 非洲黑人性xxxx精品又粗又长| 午夜免费观看网址| 国产精品99久久久久久久久| 午夜精品久久久久久毛片777| 男人舔女人下体高潮全视频| 人人妻人人看人人澡| 国模一区二区三区四区视频 | 熟女少妇亚洲综合色aaa.| 欧美日韩中文字幕国产精品一区二区三区| 黑人巨大精品欧美一区二区mp4| 日本黄色片子视频| 国产精品一区二区三区四区久久| 成人国产综合亚洲| 很黄的视频免费| 国产精品久久久久久精品电影| netflix在线观看网站| 美女cb高潮喷水在线观看 | 精品一区二区三区av网在线观看| 99国产精品一区二区蜜桃av| 嫩草影视91久久| 两个人看的免费小视频| 精品免费久久久久久久清纯| 日韩人妻高清精品专区| 国产又黄又爽又无遮挡在线| 亚洲欧美一区二区三区黑人| 每晚都被弄得嗷嗷叫到高潮| 成人18禁在线播放| 一个人免费在线观看电影 | 欧美中文日本在线观看视频| 热99在线观看视频| 亚洲在线自拍视频| 欧美一区二区国产精品久久精品| 日本三级黄在线观看| 1024香蕉在线观看| 日日夜夜操网爽| 极品教师在线免费播放| 动漫黄色视频在线观看| www.精华液| 亚洲av电影在线进入| 国产精品久久久久久久电影 | 久久精品亚洲精品国产色婷小说| 操出白浆在线播放|