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

    赤道太平洋西風(fēng)爆發(fā)事件的時(shí)空特征

    2022-06-22 03:04:12王法明
    海洋科學(xué) 2022年5期
    關(guān)鍵詞:緯向厄爾尼諾西風(fēng)

    莊 鎮(zhèn) , 鄭 建 , , 王法明 ,

    (1. 中國科學(xué)院海洋研究所, 山東 青島 266071; 2. 中國科學(xué)院大學(xué), 北京 100049; 3. 中國科學(xué)院海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266071; 4. 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室 海洋動(dòng)力過程與氣候功能實(shí)驗(yàn)室, 山東 青島 266237)

    西風(fēng)事件是發(fā)生在熱帶的高頻西風(fēng)異常擾動(dòng)。持續(xù)的西風(fēng)異??梢援a(chǎn)生下沉的海洋開爾文波, 從而觸發(fā)和維持厄爾尼諾事件[1]。如果西風(fēng)事件減弱甚至消失, 厄爾尼諾事件就很難形成。例如, 2014年厄爾尼諾事件沒有成長起來, 很大程度上歸因于當(dāng)年弱的西風(fēng)異常[2-3]。由于厄爾尼諾會(huì)通過大氣遙相關(guān)影響全球氣候與天氣異常, 因此研究西風(fēng)事件有著重要科學(xué)意義。

    目前, 關(guān)于西風(fēng)事件并沒有公認(rèn)的嚴(yán)格標(biāo)準(zhǔn),已有的研究都是基于各自經(jīng)驗(yàn)定義西風(fēng)事件。從天氣學(xué)的角度看, 臺風(fēng)、赤道外寒潮的入侵和大氣季節(jié)內(nèi)振蕩都可以促使西風(fēng)事件發(fā)生[4-8]; 而從氣候?qū)W的角度看, 西風(fēng)事件的發(fā)生頻次和強(qiáng)度則受到El Ni?o-Southern Oscillation(ENSO)信號的調(diào)制[8-21]。例如, Eisenman等[9]發(fā)現(xiàn), 厄爾尼諾事件期間, 隨著暖池東邊界向東擴(kuò)展, 西風(fēng)異常區(qū)域也隨之向東移動(dòng)。作為觸發(fā)ENSO的隨機(jī)天氣擾動(dòng), 厘清西風(fēng)事件如何受到氣候年際變化信號的調(diào)控是其中研究的關(guān)鍵。針對這一科學(xué)問題, 本文將基于海洋大氣再分析資料, 定義西風(fēng)事件指數(shù), 診斷熱帶太平洋西風(fēng)異常的時(shí)空特征, 并探討ENSO信號對西風(fēng)事件的影響。

    1 數(shù)據(jù)和方法

    1.1 研究數(shù)據(jù)

    西風(fēng)事件屬于高頻的天氣事件, 進(jìn)行分析時(shí)需采用逐日風(fēng)場資料。為保證統(tǒng)計(jì)結(jié)果的可靠性, 采用兩套風(fēng)場資料數(shù)據(jù)相互驗(yàn)證, 分別是European Centre for Medium-Range Weather Forecas Reanalysis-Interim(以下簡稱ERA-Interim)和National Centers for Environmental Prediction Reanalysis II(以下簡稱NCEP2)。ERA-Interim數(shù)據(jù)分辨率是2.5°×2.5°, NCEP2數(shù)據(jù)分辨率是1.875°×1.875°。風(fēng)場數(shù)據(jù)的時(shí)間范圍是1979年1月1日到2015年12月31日。海表面溫度數(shù)據(jù)采用National Oceanic and Atmospheric Administration(NOAA) Optimum Interpolation (OI) 逐月數(shù)據(jù), 時(shí)間范圍是1982年1月到2015年12月。大氣和海洋數(shù)據(jù)均先減去氣候平均態(tài)得到異常場后再做分析。需要注意的是, 此時(shí)的異常風(fēng)場既包含了大氣自有的高頻天氣變化, 也包含了低頻大尺度海氣耦合信號。這里, 我們借鑒其他研究的做法, 采用帶通濾波提取2~90 d的高頻信號, 在該風(fēng)場異常數(shù)據(jù)中識別西風(fēng)事件。

    1.2 西風(fēng)事件的定義

    西風(fēng)事件由一系列的西風(fēng)異常組成, 它需要達(dá)到一定的空間尺度, 才能促進(jìn)海洋開爾文波的產(chǎn)生。目前, 對于西風(fēng)事件的定義存在較大分歧[6-9]。我們將從個(gè)例分析入手, 導(dǎo)出本文所用的西風(fēng)事件定義。風(fēng)場數(shù)據(jù)的預(yù)處理方法是在2.5°S與2.5°N之間取經(jīng)向平均。該范圍對應(yīng)的是赤道海洋開爾文波的衰減半徑, 可保證影響赤道海洋開爾文波的西風(fēng)異常都涵蓋在此研究區(qū)域內(nèi)。圖1a表示2002年緯向風(fēng)異常隨時(shí)間演變, 其中黑色框?qū)?yīng)圖1b典型西風(fēng)事件。圖1a顯示, 西風(fēng)異常往往集中出現(xiàn),但其持續(xù)時(shí)間和緯向跨度卻各有不同。對于持續(xù)天數(shù)多、作用范圍廣的一系列西風(fēng)異常, 需要給出具體標(biāo)準(zhǔn)定義西風(fēng)事件。以2002年5月西風(fēng)事件為例, 其特征可歸納為緯向?qū)挾却?、持續(xù)時(shí)間長、空間移動(dòng)不顯著、并且具有一定的平均強(qiáng)度。因此, 西風(fēng)事件的具體定義依賴于時(shí)間、空間特征尺度的具體選取。

    圖1 2002年緯向風(fēng)場異常(a), 5月的西風(fēng)事件個(gè)例(b)Fig. 1 Zonal wind anomaly in the year 2002 (a) and westerly wind bursts in May 2002 (b)

    圖2顯示西風(fēng)事件個(gè)數(shù)與選取標(biāo)準(zhǔn)的關(guān)系。緯向?qū)挾葮?biāo)準(zhǔn)過小會(huì)將部分西風(fēng)事件連在一起, 難以區(qū)分; 如果寬度過大會(huì)忽略部分西風(fēng)事件, 從而導(dǎo)致西風(fēng)事件個(gè)數(shù)減少。圖2a表明選取不同的持續(xù)時(shí)間, 緯向?qū)挾葮?biāo)準(zhǔn)為30°時(shí)可以區(qū)分出大部分的西風(fēng)事件。從物理意義上講, 赤道大氣開爾文波的空間特征尺度為30°左右, 因此選擇30°作為西風(fēng)事件緯向?qū)挾鹊臉?biāo)準(zhǔn)具有一定合理性。圖2b表明, 在不同的緯向?qū)挾葮?biāo)準(zhǔn)下, 隨著持續(xù)天數(shù)標(biāo)準(zhǔn)的增加,西風(fēng)事件的個(gè)數(shù)都會(huì)呈現(xiàn)一定的下降趨勢。鑒于大氣開爾文波的特征時(shí)間尺度為一周左右, 選擇5 d作為西風(fēng)事件持續(xù)天數(shù)的標(biāo)準(zhǔn)。另外, 為避免同一時(shí)間段而空間相隔甚遠(yuǎn)的西風(fēng)事件發(fā)生混淆, 需要確定西風(fēng)事件中心緯向位移標(biāo)準(zhǔn)(即第1 d和第2 d西風(fēng)事件中心位置的移動(dòng)距離, 圖1b中兩個(gè)圓圈的間隔)。圖2c表明, 如果西風(fēng)事件中心位移標(biāo)準(zhǔn)設(shè)置得較小, 會(huì)導(dǎo)致西風(fēng)事件個(gè)數(shù)較少, 當(dāng)位移標(biāo)準(zhǔn)逐漸增大時(shí), 西風(fēng)事件個(gè)數(shù)逐漸穩(wěn)定并趨近正常情況。從數(shù)據(jù)分布來看選擇60°作為中心位移的閾值較合適。

    圖2 西風(fēng)事件個(gè)數(shù)與緯向?qū)挾葮?biāo)準(zhǔn)的關(guān)系(a); 持續(xù)天數(shù)標(biāo)準(zhǔn)的關(guān)系(b); 緯向移動(dòng)距離的關(guān)系(c); 隨振幅的概率分布(d)Fig. 2 Relationship between the number of westerly wind bursts and zonal wind criteria (a), duration criteria (b), zonal shift (c),and strength (d)

    赤道海洋開爾文波東傳速度約為2 m/s, 因此可選取3 m/s作為西風(fēng)事件強(qiáng)度的最低閾值。圖2d也表明大部分西風(fēng)事件振幅都大于3 m/s。綜上所述,本文西風(fēng)事件的標(biāo)準(zhǔn)是西風(fēng)異常的緯向?qū)挾却笥?0°, 持續(xù)時(shí)間大于5 d, 中心移動(dòng)距離不超過60°,振幅大于3 m/s。

    2 西風(fēng)事件的時(shí)空特征

    根據(jù)上節(jié)的西風(fēng)事件判定標(biāo)準(zhǔn), 本節(jié)將對赤道太平洋西風(fēng)事件的時(shí)空特征進(jìn)行診斷, 重點(diǎn)從緯向分布特征、季節(jié)變化和年際變化3個(gè)角度開展分析。下文主要采用ERA-Interim數(shù)據(jù)結(jié)果進(jìn)行討論(NCEP2數(shù)據(jù)結(jié)果類似)。

    2.1 西風(fēng)事件的緯向分布

    西風(fēng)事件是一系列西風(fēng)異常, 而西風(fēng)異常方差的緯向分布則代表了天氣噪聲擾動(dòng)的空間差異。圖3表明, 西風(fēng)異常方差集中在140°E附近, 與東風(fēng)異常方差分布相似, 但是西風(fēng)異常的方差要大于東風(fēng)異常。在西風(fēng)事件識別的基礎(chǔ)上, 以事件初期的中心位置(見圖1b圓點(diǎn))來表征西風(fēng)事件的緯向分布。圖4顯示西風(fēng)事件主要集中在140°E和160°W兩個(gè)區(qū)域附近。這與Seiki等[6]西風(fēng)事件分布在140°E附近和170°E~160°W之間的論斷是一致的。由于西風(fēng)異常的能量在140°E附近最大, 因此西風(fēng)事件也會(huì)在這個(gè)區(qū)域集中發(fā)生。而在160°W附近, 西風(fēng)事件的集中存在可能與ENSO信號有關(guān)(見2.3節(jié))。綜合前面的分析, 赤道太平洋西風(fēng)事件主要分布在太平洋中西部, 作用中心分別位于140°E和160°W附近。

    圖3 赤道太平洋西風(fēng)異常和東風(fēng)異常的緯向分布Fig. 3 Zonal distribution of easterly and westerly anomalies in the equatorial Pacific

    圖4 西風(fēng)事件的緯向分布特征Fig. 4 Zonal distributions of westerly wind bursts

    2.2 西風(fēng)事件的時(shí)間特征

    西風(fēng)事件在時(shí)間演化方向上除了隨機(jī)性, 也存在一定的規(guī)律性。例如, 當(dāng)熱帶大氣季節(jié)內(nèi)振蕩(Madden-Julian Oscillation)為濕位相時(shí), 往往伴隨著西風(fēng)事件的產(chǎn)生[10]; 而當(dāng)ENSO為暖相位時(shí), 西風(fēng)事件的數(shù)量增多[8]。下文將從季節(jié)變化和年際變化兩方面分析西風(fēng)事件的時(shí)間特征。

    2.2.1 西風(fēng)事件的季節(jié)特征

    基于前文對西風(fēng)事件的識別, 利用西風(fēng)事件的累積天數(shù)代表西風(fēng)事件的頻率, 西風(fēng)事件中西風(fēng)異常的平均值來代表西風(fēng)事件的平均強(qiáng)度, 將數(shù)據(jù)統(tǒng)一映射到1月—12月的時(shí)間坐標(biāo)上, 以此來表征多年平均下西風(fēng)事件的季節(jié)變化。圖5給出西風(fēng)事件在1979年1月1日到2015年12月31日期間累積的西風(fēng)事件天數(shù)和西風(fēng)事件平均強(qiáng)度。數(shù)據(jù)顯示, 冬季和春季西風(fēng)事件爆發(fā)的頻次較高, 平均強(qiáng)度較大,主要集中在中西太平洋, 與Seiki等[6]研究結(jié)果一致。相比于冬季和春季, 夏季和秋季西風(fēng)事件的頻率較低, 而且其平均強(qiáng)度偏弱。其中, 春季和夏季的西風(fēng)事件對厄爾尼諾有觸發(fā)作用, 可以產(chǎn)生下沉開爾文波, 將海溫異常暖信號向東傳遞; 而秋季和冬季的西風(fēng)事件則對應(yīng)著ENSO信號盛期大尺度背景風(fēng)場對天氣事件的調(diào)制, 使得大氣渦動(dòng)能從低頻向高頻串級[13]。從西風(fēng)事件的緯向分布來看, 其季節(jié)性變化不明顯, 主要集中在140°E~160°W。由上可知, 西風(fēng)事件在冬季和春季較強(qiáng), 夏季和秋季較弱。

    圖5 西風(fēng)事件累積時(shí)間(a)和平均強(qiáng)度的時(shí)間-經(jīng)度剖面圖(b)Fig. 5 Accummulation days (a) and the mean amplitude of westerly wind burst in time-longitude plots (b)

    2.2.2 西風(fēng)事件的年際變化

    與前面的季節(jié)變化一樣, 西風(fēng)事件的年際變化也從頻率和強(qiáng)度兩個(gè)方面展開探討。不同的是, 將每月西風(fēng)事件的累積天數(shù)定義為逐月西風(fēng)事件天數(shù)(N)來表征西風(fēng)事件的頻率, 每月西風(fēng)事件的平均振幅定義為西風(fēng)事件平均振幅(A)來表征西風(fēng)事件的強(qiáng)度,生成各年各月的兩條時(shí)間序列。為了更好的突出年際變化信號, 上述兩個(gè)指標(biāo)均采用6個(gè)月滑動(dòng)平均處理。其中, 西風(fēng)事件月平均振幅的計(jì)算是基于逐日資料, 若 maxU′ 代表西風(fēng)事件時(shí)每天最大西風(fēng)異常值,每一個(gè)月具體的計(jì)算公式如下,

    求和計(jì)算的時(shí)間范圍是從月初到月底。

    圖6a顯示在較強(qiáng)的厄爾尼諾年(1982/1983、1986/1987、1997/1998、2015/2016年)前后, 逐月西風(fēng)事件天數(shù)都有先逐漸增加后減少的趨勢; 而在非厄爾尼諾年, 逐月西風(fēng)事件天數(shù)普遍偏少。這說明西風(fēng)事件的發(fā)生頻次與ENSO有一定的關(guān)聯(lián)。圖7功率譜存在通過95%置信度檢驗(yàn)的譜峰, 說明西風(fēng)事件發(fā)生頻次有準(zhǔn)2 a周期的變化。從西風(fēng)事件的強(qiáng)度來看(圖6b), 在厄爾尼諾年, 平均振幅較大, 最高可以達(dá)到3 m/s; 而在非厄爾尼諾年, 逐月西風(fēng)事件平均振幅則維持在較弱的水平。圖7b功率譜分析結(jié)果也顯示西風(fēng)事件的強(qiáng)度存在類似ENSO的寬譜峰年際變化。總體來說, 西風(fēng)事件的頻次和強(qiáng)度都存在一定的年際變化, 而強(qiáng)度與ENSO信號的周期更吻合。

    圖6 逐月西風(fēng)事件時(shí)間(a)和逐月西風(fēng)事件平均振幅(b)Fig. 6 Interannual variability in monthly westerly wind burst in day numbers (a) and mean amplitude (b)

    圖7 西風(fēng)事件時(shí)間功率譜(a)和西風(fēng)事件平均振幅功率譜(b)Fig. 7 Power spectrum of westerly wind bursts according to day numbers (a) and mean amplitude (b)

    由前述分析可知, 西風(fēng)事件的概率、持續(xù)時(shí)間和緯向?qū)挾鹊忍卣鲄?shù)都隱含了一部分年際變化的信息[9-10]。因此, 厘清海表面溫度異常對西風(fēng)事件的調(diào)制關(guān)系, 對于厄爾尼諾事件的模擬和預(yù)報(bào)都具有重要意義[15-16]。下文將討論西風(fēng)事件的時(shí)空特征如何具體受ENSO信號的影響。

    2.3 西風(fēng)事件時(shí)空特征受ENSO信號的影響

    以逐月緯向風(fēng)異常的方差來表示西風(fēng)事件的高頻大氣變化, 研究西風(fēng)事件時(shí)空特征受ENSO信號的影響時(shí), ERA-Interim和NCEP2的風(fēng)場數(shù)據(jù)給出相似的結(jié)果, 因此本節(jié)僅以ERA-Interim風(fēng)場資料為例進(jìn)行討論。

    以緯向風(fēng)異常逐月方差為左場、逐月海表面溫度異常為右場進(jìn)行奇異值分解, 第一模態(tài)累積解釋協(xié)方差56.6%, 說明海面風(fēng)場和溫度異常場存在顯著的空間相關(guān)。圖8b顯示影響高頻緯向風(fēng)場的海表面溫度異常信號源于ENSO, 而圖8a則表明與ENSO信號相關(guān)的緯向風(fēng)高頻變化主要集中在中太平洋地區(qū)。這種西風(fēng)事件與ENSO信號的相干關(guān)系在不同的研究里都得到了證實(shí)[14], 說明它是熱帶海氣耦合系統(tǒng)的一個(gè)固有特征。圖8c時(shí)間系數(shù)表明, 隨著海表面溫度異常的升高, 緯向風(fēng)異常的高頻變化在逐漸增強(qiáng)。關(guān)于ENSO信號影響大氣高頻變化的物理機(jī)制, Seiki等[7]認(rèn)為ENSO信號影響西風(fēng)事件主要通過改變背景風(fēng)場的水平和豎直剪切從而實(shí)現(xiàn)渦動(dòng)能從低頻向高頻的傳遞。對前文中西風(fēng)事件的位置進(jìn)行區(qū)分, 可以看出厄爾尼諾期間的西風(fēng)事件主要分布在太平洋中部; 而非厄爾尼諾期間, 西風(fēng)事件在西太平洋地區(qū)更集中。綜上所述, 西風(fēng)事件受到ENSO信號的影響, 更多集中在中太平洋地區(qū)。

    圖8 奇異值分解第一模態(tài)的緯向風(fēng)(a), 海面溫度(b), 時(shí)間系數(shù)(c)Fig. 8 First singular value decomposition mode in tropical Pacific zonal winds (a), monthly SST (b), and time series (c)

    為進(jìn)一步研究西風(fēng)事件的年際變化與ENSO信號的關(guān)系, 圖9給出Ni?o3.4指數(shù)與逐月西風(fēng)事件累積天數(shù)和平均振幅的散點(diǎn)圖。其中Ni?o3.4指數(shù)表征ENSO信號, 而西風(fēng)事件年際變化則根據(jù)前文逐月西風(fēng)事件天數(shù)和西風(fēng)事件平均振幅來判別。圖9表明, 隨著Ni?o3.4指數(shù)的增加, 逐月西風(fēng)事件天數(shù)也逐漸增加, 西風(fēng)事件頻率與Ni?o3.4指數(shù)的線性相關(guān)系數(shù)是0.31。圖9a中Ni?o3.4指數(shù)集中分布在非對稱的冷暖區(qū)間–2~3 ℃。在–1~1 ℃范圍內(nèi),逐月西風(fēng)事件天數(shù)與Ni?o3.4指數(shù)的相關(guān)性較弱,離散程度更高。而圖9b表明逐月西風(fēng)事件的平均振幅與Ni?o3.4指數(shù)的線性相關(guān)性更好, 相關(guān)系數(shù)達(dá)0.56。平均振幅與海表面溫度的線性相關(guān)系數(shù)大,說明西風(fēng)事件的強(qiáng)度更容易受到ENSO信號的調(diào)制作用。

    圖9 Ni?o3.4指數(shù)與西風(fēng)事件時(shí)間(a), 西風(fēng)事件平均振幅關(guān)系的散點(diǎn)圖(b)(振幅單位: m/s; Ni?o3.4指數(shù)單位: ℃)Fig. 9 Scatter plots of Ni?o3.4 index and westerly wind burst in day numbers (a) and mean amplitude (b)

    初步猜想ENSO信號影響西風(fēng)事件需要達(dá)到一定的閾值: 當(dāng)Ni?o3.4指數(shù)超過這個(gè)閾值, 大尺度海氣耦合作用明顯, 大氣的低頻能量可以持續(xù)向高頻能量轉(zhuǎn)移, 促使西風(fēng)事件的爆發(fā); 而當(dāng)?shù)陀谶@個(gè)閾值, 低頻能量向高頻的轉(zhuǎn)移不能持續(xù), 使得西風(fēng)事件隨機(jī)發(fā)生, 但是其強(qiáng)度還是會(huì)受到能量傳遞的控制。具體來說, 當(dāng)Ni?o3.4指數(shù)大于1 ℃時(shí), 西風(fēng)事件更易受到ENSO信號的影響, 而當(dāng)Ni?o3.4指數(shù)處于–1~1 ℃時(shí), 西風(fēng)事件的頻率受ENSO信號的影響較小, 呈現(xiàn)很強(qiáng)的隨機(jī)性, 但是西風(fēng)事件強(qiáng)度受到一定的ENSO信號調(diào)控。

    3 結(jié)論

    本文利用ERA-Interim和NCEP2的逐日風(fēng)場資料以及OI逐月海表面溫度資料診斷了西風(fēng)事件爆發(fā)的時(shí)空特征, 并探討了ENSO信號對其發(fā)生頻次和強(qiáng)度的影響, 得出以下初步結(jié)論:

    1) 赤道太平洋西風(fēng)事件集中在140°E和160°W附近。冬季和春季西風(fēng)事件頻率和平均強(qiáng)度較大, 夏季和秋季較小, 緯向?qū)挾鹊募竟?jié)變化不顯著。

    2) 西風(fēng)事件存在明顯的年際變化, 譬如厄爾尼諾期間, 中東太平洋海表面溫度顯著異常升高, 西風(fēng)事件也相應(yīng)的集中在此區(qū)域, 尤以中太平洋為主。當(dāng)Ni?o3.4指數(shù)大于1 ℃時(shí), 西風(fēng)事件的頻率和強(qiáng)度與ENSO信號的相干關(guān)系密切; 當(dāng)Ni?o3.4指數(shù)小于1 ℃時(shí), 西風(fēng)事件頻率的隨機(jī)性增強(qiáng), 但其強(qiáng)度仍能受ENSO信號調(diào)制。

    需要說明的是, 本文西風(fēng)事件時(shí)空特征受ENSO信號影響是統(tǒng)計(jì)結(jié)果, 其中利用帶通濾波分離大尺度海氣耦合信號和高頻天氣噪聲的做法具有一定的局限性。在此基礎(chǔ)上, 如果借助模式實(shí)驗(yàn), 能更準(zhǔn)確地區(qū)分ENSO信號對西風(fēng)事件的影響, 揭示其物理機(jī)制。

    猜你喜歡
    緯向厄爾尼諾西風(fēng)
    古道西風(fēng)(水彩)
    紗線強(qiáng)力對純棉平紋面料強(qiáng)力的影響
    利用掩星溫度數(shù)據(jù)推算大氣月平均緯向風(fēng)場
    秋漸濃,分付西風(fēng)此夜涼
    文苑(2019年20期)2019-11-16 08:52:46
    厄爾尼諾現(xiàn)象橫行中國
    浣溪沙 誰念西風(fēng)獨(dú)自涼
    飛魔幻A(2019年3期)2019-06-12 09:22:35
    溫度對絲綢面料粘襯熱縮率的影響
    絲綢(2018年10期)2018-10-15 09:54:16
    柞蠶絲面料在粘襯過程中的熱縮率變化分析
    厄爾尼諾現(xiàn)象及其影響
    三步圖解厄爾尼諾與拉尼娜
    麻豆国产97在线/欧美| 国产国拍精品亚洲av在线观看| 亚洲美女黄色视频免费看| 亚洲欧美精品专区久久| 天天躁夜夜躁狠狠久久av| 亚洲熟女精品中文字幕| 亚洲怡红院男人天堂| 成人国产av品久久久| 在线观看国产h片| 国产精品嫩草影院av在线观看| 国产在线一区二区三区精| 色5月婷婷丁香| 男人舔奶头视频| 日韩中文字幕视频在线看片 | av福利片在线观看| 日韩制服骚丝袜av| 高清不卡的av网站| 欧美成人精品欧美一级黄| 久久精品久久久久久久性| 日本vs欧美在线观看视频 | 成人高潮视频无遮挡免费网站| 黄色配什么色好看| .国产精品久久| 三级国产精品欧美在线观看| 在线看a的网站| 国产亚洲5aaaaa淫片| 男女下面进入的视频免费午夜| 一本色道久久久久久精品综合| 夫妻午夜视频| 黄色怎么调成土黄色| 亚洲aⅴ乱码一区二区在线播放| 国语对白做爰xxxⅹ性视频网站| 欧美精品亚洲一区二区| 又爽又黄a免费视频| 国产毛片在线视频| 在线观看一区二区三区| 一级毛片我不卡| 欧美成人精品欧美一级黄| 欧美激情国产日韩精品一区| 3wmmmm亚洲av在线观看| 好男人视频免费观看在线| 国产无遮挡羞羞视频在线观看| 极品少妇高潮喷水抽搐| 欧美高清成人免费视频www| 97超碰精品成人国产| 十八禁网站网址无遮挡 | 亚洲欧美精品自产自拍| 久久99精品国语久久久| 日韩伦理黄色片| 亚洲伊人久久精品综合| 一本色道久久久久久精品综合| 777米奇影视久久| 伦理电影免费视频| 国产免费视频播放在线视频| 伊人久久国产一区二区| 午夜福利网站1000一区二区三区| 一级毛片 在线播放| a级毛片免费高清观看在线播放| 少妇人妻 视频| 国产爱豆传媒在线观看| 亚洲美女搞黄在线观看| 国产亚洲一区二区精品| 亚洲不卡免费看| 夜夜爽夜夜爽视频| 伦精品一区二区三区| 欧美区成人在线视频| 菩萨蛮人人尽说江南好唐韦庄| 国产黄片视频在线免费观看| 免费黄网站久久成人精品| 国产伦理片在线播放av一区| 3wmmmm亚洲av在线观看| 免费观看的影片在线观看| 国产精品麻豆人妻色哟哟久久| 午夜福利在线观看免费完整高清在| 免费不卡的大黄色大毛片视频在线观看| 亚洲综合精品二区| 在线观看人妻少妇| 永久免费av网站大全| 在线观看三级黄色| 狂野欧美白嫩少妇大欣赏| 人妻 亚洲 视频| 男女啪啪激烈高潮av片| 少妇裸体淫交视频免费看高清| 亚洲人与动物交配视频| 国产精品人妻久久久久久| 大话2 男鬼变身卡| 亚洲精品色激情综合| 国产男人的电影天堂91| 18+在线观看网站| 亚洲精华国产精华液的使用体验| 国产在线视频一区二区| 国产男女超爽视频在线观看| 亚洲精品视频女| 成人黄色视频免费在线看| 婷婷色综合大香蕉| 中国三级夫妇交换| 搡老乐熟女国产| 最近手机中文字幕大全| 妹子高潮喷水视频| 成年女人在线观看亚洲视频| 汤姆久久久久久久影院中文字幕| 在线观看免费高清a一片| 中文在线观看免费www的网站| 亚洲电影在线观看av| 精品午夜福利在线看| 中文乱码字字幕精品一区二区三区| 99久久精品一区二区三区| 亚洲精品日韩在线中文字幕| 欧美zozozo另类| 久久久久久久久久久丰满| 全区人妻精品视频| 国产精品一区二区在线不卡| 直男gayav资源| 高清午夜精品一区二区三区| 91久久精品国产一区二区三区| 天堂俺去俺来也www色官网| 在线精品无人区一区二区三 | 国产乱来视频区| 久久久久久久久久成人| 一本色道久久久久久精品综合| 国产片特级美女逼逼视频| 久久久精品免费免费高清| 亚洲精品日韩在线中文字幕| 国产成人a∨麻豆精品| 老师上课跳d突然被开到最大视频| 日本黄大片高清| 黑人猛操日本美女一级片| a级毛色黄片| 免费大片黄手机在线观看| 国产伦精品一区二区三区视频9| 日韩欧美精品免费久久| 久久久欧美国产精品| 亚洲欧美日韩无卡精品| 国产精品久久久久成人av| 一级毛片aaaaaa免费看小| 蜜臀久久99精品久久宅男| 青春草亚洲视频在线观看| 五月天丁香电影| 久久人人爽av亚洲精品天堂 | 国产人妻一区二区三区在| 国产久久久一区二区三区| 成年美女黄网站色视频大全免费 | 亚洲国产精品专区欧美| 欧美xxxx性猛交bbbb| 久久av网站| 日本色播在线视频| 亚洲国产色片| 国产精品一区二区三区四区免费观看| 2018国产大陆天天弄谢| 18禁动态无遮挡网站| 久久久久久久久久人人人人人人| 777米奇影视久久| 麻豆精品久久久久久蜜桃| 九草在线视频观看| 国产精品欧美亚洲77777| 干丝袜人妻中文字幕| 波野结衣二区三区在线| 亚洲国产最新在线播放| 午夜免费观看性视频| 亚洲欧美精品专区久久| 99久久精品国产国产毛片| av视频免费观看在线观看| 国产 一区精品| 亚洲高清免费不卡视频| 人妻 亚洲 视频| 欧美日韩精品成人综合77777| 婷婷色综合www| 赤兔流量卡办理| 久久精品国产自在天天线| 高清欧美精品videossex| kizo精华| 国产欧美另类精品又又久久亚洲欧美| 蜜桃亚洲精品一区二区三区| 欧美精品一区二区大全| 18禁在线无遮挡免费观看视频| 成年女人在线观看亚洲视频| 菩萨蛮人人尽说江南好唐韦庄| 下体分泌物呈黄色| 波野结衣二区三区在线| 免费av中文字幕在线| 国产一级毛片在线| 寂寞人妻少妇视频99o| 久久影院123| 亚洲av在线观看美女高潮| 国产高潮美女av| 久久久久久久精品精品| 秋霞伦理黄片| 一区二区av电影网| a级毛色黄片| 欧美最新免费一区二区三区| 一级二级三级毛片免费看| 直男gayav资源| 91午夜精品亚洲一区二区三区| 啦啦啦视频在线资源免费观看| 精品久久久噜噜| 免费人妻精品一区二区三区视频| 欧美成人一区二区免费高清观看| 国产精品人妻久久久久久| 亚洲欧美成人综合另类久久久| 国产淫片久久久久久久久| 在线精品无人区一区二区三 | 尾随美女入室| 啦啦啦在线观看免费高清www| 免费看日本二区| 免费观看av网站的网址| 久久青草综合色| 亚洲第一av免费看| 黄片wwwwww| 国产精品久久久久久精品电影小说 | 精品人妻偷拍中文字幕| 18禁裸乳无遮挡免费网站照片| 亚洲精华国产精华液的使用体验| 欧美精品国产亚洲| 久久精品国产亚洲av天美| 久久热精品热| 欧美日韩国产mv在线观看视频 | 卡戴珊不雅视频在线播放| h视频一区二区三区| 国产一区二区在线观看日韩| a级毛色黄片| 精品国产露脸久久av麻豆| 久久ye,这里只有精品| 亚洲精品一区蜜桃| 一区在线观看完整版| 毛片女人毛片| 成人亚洲精品一区在线观看 | 又黄又爽又刺激的免费视频.| 蜜臀久久99精品久久宅男| 午夜福利影视在线免费观看| 国产成人freesex在线| 午夜激情久久久久久久| 国产乱人偷精品视频| 日日啪夜夜撸| a 毛片基地| 国产黄片美女视频| 天美传媒精品一区二区| 80岁老熟妇乱子伦牲交| 日韩中文字幕视频在线看片 | 寂寞人妻少妇视频99o| 免费观看性生交大片5| 大又大粗又爽又黄少妇毛片口| 日韩伦理黄色片| 51国产日韩欧美| 网址你懂的国产日韩在线| 日韩强制内射视频| 一级毛片aaaaaa免费看小| 国产精品秋霞免费鲁丝片| 建设人人有责人人尽责人人享有的 | 一区二区av电影网| 香蕉精品网在线| 久久久亚洲精品成人影院| 97热精品久久久久久| 午夜视频国产福利| 在线观看三级黄色| 丝袜喷水一区| 久久99热这里只有精品18| 国产视频内射| 国产女主播在线喷水免费视频网站| 国产成人a∨麻豆精品| 欧美老熟妇乱子伦牲交| 日本猛色少妇xxxxx猛交久久| 丝瓜视频免费看黄片| 久久99精品国语久久久| 久久久久国产精品人妻一区二区| 在线看a的网站| 男男h啪啪无遮挡| 国产成人午夜福利电影在线观看| 有码 亚洲区| 涩涩av久久男人的天堂| 精品国产三级普通话版| 成人二区视频| 国产精品成人在线| 久久精品夜色国产| 精品人妻偷拍中文字幕| 下体分泌物呈黄色| 欧美日韩一区二区视频在线观看视频在线| 国产av国产精品国产| 少妇的逼水好多| 最近最新中文字幕免费大全7| 国产v大片淫在线免费观看| 久久精品人妻少妇| 欧美精品亚洲一区二区| 中文天堂在线官网| 在线 av 中文字幕| 成年美女黄网站色视频大全免费 | 亚洲国产精品一区三区| 如何舔出高潮| 老司机影院毛片| 草草在线视频免费看| 国产av国产精品国产| 成年av动漫网址| 国产黄频视频在线观看| 狂野欧美白嫩少妇大欣赏| 欧美xxⅹ黑人| 亚洲av免费高清在线观看| 久久鲁丝午夜福利片| 激情五月婷婷亚洲| 赤兔流量卡办理| 日韩三级伦理在线观看| 国产高清三级在线| 99久久精品热视频| 国产欧美日韩一区二区三区在线 | 国产白丝娇喘喷水9色精品| 深爱激情五月婷婷| 少妇的逼水好多| 国产欧美另类精品又又久久亚洲欧美| av在线app专区| 日韩av在线免费看完整版不卡| 欧美日韩国产mv在线观看视频 | 中文字幕精品免费在线观看视频 | 日韩人妻高清精品专区| 下体分泌物呈黄色| 成人毛片a级毛片在线播放| 能在线免费看毛片的网站| 国产在线免费精品| av国产精品久久久久影院| 尾随美女入室| 超碰av人人做人人爽久久| a级一级毛片免费在线观看| 精品午夜福利在线看| 精品亚洲成a人片在线观看 | 国产又色又爽无遮挡免| 激情五月婷婷亚洲| 国产欧美日韩精品一区二区| 精品久久国产蜜桃| 成人一区二区视频在线观看| 精品久久国产蜜桃| 午夜视频国产福利| 伦理电影免费视频| 91精品一卡2卡3卡4卡| 欧美另类一区| 午夜免费男女啪啪视频观看| 久久ye,这里只有精品| 午夜免费观看性视频| 国产成人精品婷婷| 国产欧美日韩精品一区二区| 国国产精品蜜臀av免费| 国产在线男女| 亚洲精品国产成人久久av| 大又大粗又爽又黄少妇毛片口| 中国三级夫妇交换| 国产成人精品一,二区| 丰满人妻一区二区三区视频av| 国产在线一区二区三区精| 免费播放大片免费观看视频在线观看| 久久久久久人妻| 26uuu在线亚洲综合色| 亚洲精品456在线播放app| 一区二区三区精品91| 少妇 在线观看| av网站免费在线观看视频| 一本一本综合久久| 亚洲av欧美aⅴ国产| 久久久久网色| 久久久亚洲精品成人影院| 简卡轻食公司| 亚洲精品日韩在线中文字幕| 国产精品一二三区在线看| 亚洲一区二区三区欧美精品| 麻豆精品久久久久久蜜桃| 日日啪夜夜撸| 男人狂女人下面高潮的视频| 观看av在线不卡| 亚洲成人一二三区av| 亚洲欧美日韩无卡精品| 中文字幕亚洲精品专区| 搡女人真爽免费视频火全软件| 久久久久久人妻| 日韩精品有码人妻一区| 精品亚洲乱码少妇综合久久| 草草在线视频免费看| 日韩 亚洲 欧美在线| .国产精品久久| 一区二区三区免费毛片| 一个人看视频在线观看www免费| 97超视频在线观看视频| 一个人免费看片子| 大片电影免费在线观看免费| 人妻 亚洲 视频| 丰满乱子伦码专区| 日韩成人伦理影院| 免费观看性生交大片5| 国产欧美亚洲国产| 大片免费播放器 马上看| 免费av不卡在线播放| videossex国产| 大陆偷拍与自拍| 国产免费一区二区三区四区乱码| 亚洲国产高清在线一区二区三| 夜夜骑夜夜射夜夜干| 九九在线视频观看精品| 亚洲欧美清纯卡通| 永久免费av网站大全| 国产精品一区www在线观看| 久久久亚洲精品成人影院| 97精品久久久久久久久久精品| 亚洲一区二区三区欧美精品| 干丝袜人妻中文字幕| 亚洲成人中文字幕在线播放| 精品99又大又爽又粗少妇毛片| 26uuu在线亚洲综合色| 亚洲不卡免费看| 久久久久久伊人网av| 国产高清三级在线| 国产亚洲午夜精品一区二区久久| 小蜜桃在线观看免费完整版高清| 国产精品欧美亚洲77777| 秋霞伦理黄片| 超碰av人人做人人爽久久| 免费看光身美女| 国产免费视频播放在线视频| 欧美日韩综合久久久久久| 国语对白做爰xxxⅹ性视频网站| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美中文字幕日韩二区| 午夜福利在线在线| 18禁动态无遮挡网站| 在线观看一区二区三区| 观看免费一级毛片| 日本猛色少妇xxxxx猛交久久| 国产伦精品一区二区三区四那| 男女国产视频网站| 久久精品国产亚洲av涩爱| 久久久久久久精品精品| 国产高清不卡午夜福利| 天天躁日日操中文字幕| 夫妻性生交免费视频一级片| 热99国产精品久久久久久7| 22中文网久久字幕| 国产无遮挡羞羞视频在线观看| 少妇被粗大猛烈的视频| 中文字幕人妻熟人妻熟丝袜美| 国产欧美另类精品又又久久亚洲欧美| 亚洲人成网站在线观看播放| 精品少妇久久久久久888优播| 亚洲色图av天堂| 国产在视频线精品| 久久精品人妻少妇| 欧美 日韩 精品 国产| 亚洲av成人精品一二三区| 最近最新中文字幕大全电影3| 五月玫瑰六月丁香| 亚洲精品一二三| 亚洲色图综合在线观看| 一级毛片黄色毛片免费观看视频| 精品人妻熟女av久视频| 麻豆成人午夜福利视频| 午夜日本视频在线| 在线亚洲精品国产二区图片欧美 | 亚洲美女黄色视频免费看| 国产91av在线免费观看| 大码成人一级视频| 久久亚洲国产成人精品v| 亚洲色图综合在线观看| 亚洲成人av在线免费| 亚洲精品久久午夜乱码| 人妻一区二区av| 日韩强制内射视频| 国产午夜精品一二区理论片| 18禁裸乳无遮挡动漫免费视频| 国产精品蜜桃在线观看| 汤姆久久久久久久影院中文字幕| 色婷婷久久久亚洲欧美| 国产av码专区亚洲av| 在线观看免费高清a一片| 免费看光身美女| 免费大片18禁| kizo精华| 久久久精品94久久精品| 欧美日韩国产mv在线观看视频 | 国产高清有码在线观看视频| 亚洲欧美成人精品一区二区| 多毛熟女@视频| 夜夜爽夜夜爽视频| 成人国产av品久久久| 日韩欧美精品免费久久| 看免费成人av毛片| 搡老乐熟女国产| 男女无遮挡免费网站观看| 亚洲怡红院男人天堂| 国产一级毛片在线| 交换朋友夫妻互换小说| 久久午夜福利片| 亚洲国产日韩一区二区| 国产高清有码在线观看视频| 男女下面进入的视频免费午夜| 99热全是精品| av网站免费在线观看视频| 色婷婷久久久亚洲欧美| 天堂中文最新版在线下载| 亚洲精品国产av成人精品| 波野结衣二区三区在线| 日韩成人av中文字幕在线观看| 日韩在线高清观看一区二区三区| 国产精品免费大片| 精品久久久噜噜| 爱豆传媒免费全集在线观看| 婷婷色av中文字幕| 欧美亚洲 丝袜 人妻 在线| 亚洲国产色片| 欧美xxⅹ黑人| 高清欧美精品videossex| 国产av精品麻豆| 99久久综合免费| 五月开心婷婷网| 日日摸夜夜添夜夜爱| 欧美激情国产日韩精品一区| 舔av片在线| 久久亚洲国产成人精品v| 18禁裸乳无遮挡动漫免费视频| 观看免费一级毛片| 国产欧美另类精品又又久久亚洲欧美| 最近最新中文字幕大全电影3| 99久久中文字幕三级久久日本| 久久久a久久爽久久v久久| 丰满乱子伦码专区| 亚洲色图av天堂| 性色avwww在线观看| 男女边摸边吃奶| 熟妇人妻不卡中文字幕| 亚洲欧洲日产国产| 五月开心婷婷网| 国产在视频线精品| 制服丝袜香蕉在线| 下体分泌物呈黄色| 日韩一本色道免费dvd| 水蜜桃什么品种好| 亚洲精品久久久久久婷婷小说| 三级国产精品片| 3wmmmm亚洲av在线观看| 国产淫语在线视频| 国精品久久久久久国模美| 国产精品一区二区性色av| 嘟嘟电影网在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品一二三| 午夜老司机福利剧场| 婷婷色麻豆天堂久久| 伊人久久国产一区二区| 久久久久人妻精品一区果冻| 色网站视频免费| 人妻系列 视频| 国产成人免费观看mmmm| 国产毛片在线视频| 2021少妇久久久久久久久久久| 只有这里有精品99| 少妇精品久久久久久久| 美女xxoo啪啪120秒动态图| 人人妻人人澡人人爽人人夜夜| 精品人妻熟女av久视频| 一区二区av电影网| 欧美97在线视频| 欧美国产精品一级二级三级 | 在线看a的网站| 狂野欧美白嫩少妇大欣赏| 街头女战士在线观看网站| 久久久久久久久久久丰满| 国产高清国产精品国产三级 | 99热6这里只有精品| 国产成人91sexporn| 中国美白少妇内射xxxbb| 免费观看的影片在线观看| 亚洲熟女精品中文字幕| 精品视频人人做人人爽| av卡一久久| 2022亚洲国产成人精品| 免费黄频网站在线观看国产| 观看美女的网站| 国产成人精品久久久久久| 中文字幕亚洲精品专区| 免费av中文字幕在线| 久久久久久久久久人人人人人人| 国产视频内射| xxx大片免费视频| 黄色配什么色好看| 久久久久国产精品人妻一区二区| 高清日韩中文字幕在线| 久久久久久久久久成人| 亚洲国产色片| 精品久久久久久电影网| 日日摸夜夜添夜夜添av毛片| 久久久色成人| 成人国产麻豆网| 在线播放无遮挡| 日本欧美国产在线视频| 精品国产乱码久久久久久小说| av视频免费观看在线观看| 性色avwww在线观看| 久久国产乱子免费精品| 亚洲丝袜综合中文字幕| 欧美xxxx黑人xx丫x性爽| 黄片无遮挡物在线观看| 精品视频人人做人人爽| 亚洲婷婷狠狠爱综合网| 亚洲欧美精品专区久久| 久久精品熟女亚洲av麻豆精品| 亚洲婷婷狠狠爱综合网| 黄片无遮挡物在线观看| 精品视频人人做人人爽| 国模一区二区三区四区视频| 深夜a级毛片| 人妻 亚洲 视频| 一本久久精品| 国产一区二区在线观看日韩| 少妇的逼水好多| 久久国内精品自在自线图片| 亚洲欧美精品专区久久| 亚洲丝袜综合中文字幕| 精品久久久精品久久久| 哪个播放器可以免费观看大片| 80岁老熟妇乱子伦牲交| 日本一二三区视频观看| 久久精品国产自在天天线| 在线播放无遮挡| 男人爽女人下面视频在线观看|