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

    深度學習在印度洋偶極子預報中的應用研究

    2022-06-01 07:15:50劉俊唐佑民宋迅殊孫志林
    大氣科學 2022年3期
    關鍵詞:北半球印度洋卷積

    劉俊 唐佑民 宋迅殊 孫志林

    1 浙江大學海洋學院, 杭州310058

    2 河海大學海洋學院自然資源部海洋災害預報重點實驗室, 南京210021

    3 自然資源部第二海洋研究所衛(wèi)星海洋環(huán)境動力學國家重點實驗室, 杭州310012

    4 南方海洋科學與工程廣東省實驗室(珠海), 珠海519082

    1 引言

    在熱帶印度洋存在兩個明顯的年際變率模態(tài):印度洋海盆模態(tài)(Indian Ocean Basin-wide Mode,IOBM)和印度洋偶極子模態(tài)(Indian Ocean Dipole Mode,IODM)。IOBM 表現為整個熱帶印度洋區(qū)域一致增暖的模態(tài),該模態(tài)一般被認為是熱帶印度洋對熱帶太平洋的厄爾尼諾—南方濤動(El Nin?o-Southern Oscillation, ENSO)的響應(Klein et al.,1999; Venzke et al., 2000)。IOD 模態(tài)表現為東西部海表溫度異常(Sea Surface Temperature anomaly,SSTA)符號相反的印度洋偶極子模態(tài)(Saji et al.,1999)。印度洋偶極子的兩極分別位于:熱帶西印度 洋(10°S~10°N,50°~70°E;簡 稱 西 極 子,WIO)與熱帶東南印度洋(10°S~0°,90°~110°E;簡稱東極子,EIO)。兩極子之間SST 距平的區(qū)域平均之差定義為印度洋偶極子指數(Indian Ocean Dipole Mode Index,DMI)。

    IODM 是印度洋氣候異常年際變率的主要模態(tài)(Saji et al., 1999; Ashok et al., 2001),也是全球氣候季節(jié)到年際變化的主要影響因子之一(Saji and Yamagata, 2003)。IOD 一般起始于5~7 月,在9~11 月達到峰值,12 月到來年1、2 月迅速消亡,呈現出明顯的季節(jié)鎖相性(McKenna et al., 2020)。IOD 事件的整個過程受到熱帶印度洋背景場的調控。以正IOD 事件為例,通常在IOD 發(fā)展的初期,即5~7 月,在熱帶印度洋夏季風的作用下,蘇門答臘島沿岸的東南風加強,邊界流引起的Ekman 輸運增強,離岸流造成上升流也增強,使得局地的溫躍層變淺。當溫躍層抬升到一定的深度后,上升流的作用使得蘇門答臘—爪哇島沿岸的SST 變冷。東印度洋的SST 負異常將進一步加強熱帶印度洋東西向的SST 梯度,SST 的梯度反過來導致大氣對流向西流動,從而進一步增強蘇門答臘—爪哇島沿岸和赤道印度洋的東風異常,此東風異常有利于蘇門答臘沿岸的SST 變冷,這個過程也稱為Bjerknes 正反饋。而在來年的冬季,由于季風轉換,背景風場轉向,熱帶東南印度洋的東南風轉為西北風,不再利于上升流的產生,從而關閉利于產生IOD事件的時間窗口,使得IOD 事件迅速消亡(Saji et al., 1999; Saji and Yamagata, 2003; Liu et al., 2017)。

    IOD 會直接或間接影響世界許多地區(qū)的氣候,例如東非、印度尼西亞、澳大利亞、中國、歐洲和日本。IOD 的正相位時期,印度洋東南部地區(qū)的冷SST 異常會影響澳大利亞的氣候,澳大利亞周圍的異常冷水會導致降雨的減少(Ashok et al.,2003);IOD 還可以通過調節(jié)太平洋遙相關對日本的氣候產生影響,并且借由大氣遙相關影響歐洲的夏季(Behera et al., 2013);IOD 還能通過調節(jié)亞洲冬季風影響我國多個地區(qū)的氣溫和降雨(Li and Mu, 2001)。IOD 對許多國家產生了巨大的社會經濟影響,因此,提前精準預測IOD 尤為重要。目前,世界上各大預報中心多使用大氣海洋耦合氣候模型來預測IOD(Luo et al., 2007; Weller and Cai,2013; Kirtman et al., 2014)。這些動力學模型對全球氣候模擬與IOD 預報起到了巨大的推動作用。目前氣候動力學模型對IOD 的實際預測時效約為一個季度左右,而理論分析表明,IOD 的潛在預測時效在兩個季度以上(Zhao and Hendon, 2009; Shi et al., 2012; Becker et al., 2014; Liu et al., 2017; Wang et al., 2017; Wu and Tang, 2019)。因此,IOD 的實際預報技巧仍有較大的提升空間。

    隨著大數據時代的來臨,深度學習由于其高效的數據處理能力(LeCun et al., 2015),在氣候科學中的應用逐漸增多(Sahai et al., 2000; Lee et al.,2018; Rasp et al., 2018; Ham et al., 2019; Ratnam et al., 2020)。深度學習一般由多個處理層來組成計算模型,在經過若干處理層的轉換之后,模型可以學到非常精細的組合模式和空間特征(LeCun et al., 2015; Krizhevsky et al., 2017)。 Ham et al.(2019)利用深度學習中的卷積神經網絡(CNN)模型對ENSO 進行了預報。他們的結果顯示利用深度學習獲得的ENSO 預報技巧要高于現有最先進的耦合氣候動力模式,證明了深度學習在預報氣候事件中的潛力。由于IOD 的預報存在冬春季和夏季預報障礙(Feng et al., 2014; Liu et al.,2018),當起始態(tài)為北半球冬春季時,預報后續(xù)的IOD 發(fā)展最為困難(Wu and Tang, 2019; 雷蕾等,2020)。因此,本文將使用深度學習中的卷積神經網絡處理再分析資料中的SSTA,以春季為起始來預報IOD 事件以及東、西極子的溫度變化。同時,我們應用人工神經網絡(ANN)中的多層感知機(MLP)進行相同的預報試驗,作為經典神經網絡的預測結果。通過比較兩種模型的預報結果,分析CNN 模型的在預報印度洋海溫上的優(yōu)勢。

    本文的主要內容如下:第二節(jié)為本文所用的數據和模型;第三節(jié)為利用CNN 和MLP 模型對DMI、EIO 指數和WIO 指數的預測結果;第四節(jié)為總結與討論。

    2 數據與模型

    2.1 數據

    本文訓練模型選取的觀測資料為美國國家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)的再分析資料Extended Reconstructed Sea Surface Temperature Version 5(ERSSTv5)(Huang et al., 2017),其中,訓練集時間范圍為1854~1989 年;驗證集時間范圍為1990~2019 年(共30 年)。SSTA 所選范圍為整個印度洋區(qū)域(70°S~25°N,25°~130°E)。我們將ERSSTv5 中連續(xù)三個月(1~3 月(JFM)、2~4 月(FMA)、3~5 月(MAM))的SSTA作為CNN 的輸入數據。模型輸出為提前時間(lead month)對應的DMI、EIOI 、WIOI。輸入和輸出的數據均是相對于氣候態(tài)的異常值,并且輸入與輸出的數據都進行了去趨勢化處理。

    2.2 CNN 模型與MLP 模型

    如圖1 所示,本文使用的CNN 模型架構包含兩個卷積池化層,一個全聯(lián)接層和一個輸出層。卷積層的作用是提取輸入SSTA 的空間特征等關鍵信息,池化層對這些信息進行放大或平均。經過兩次卷積池化處理后,全聯(lián)接層將特征圖展開并映射到合理的解空間。每個卷積層中有10~12 個過濾器,用于提取SSTA 的空間信息。過濾器在3×3 網格中提取SSTA 數據特征。平均池化層在2×2 網格中提取特征。全連接層具有8~10 個神經元。因為需要回歸來生成DMI,最終輸出層是一維的。根據卷積層中過濾器和全聯(lián)接層中神經元數目的不同,共有27 種CNN 模型的架構,分別記做C10C10D8,C10C10D9, C10C10D10, C10C11D8, ······,C12C12D10 等(C 和D 分別表示卷積層和全連接層,后面的數字代表過濾器或神經元數量)。數據批次為每次更新CNN 模型讀取的數據量大小。訓練中嘗試了將數據批次設置為12、24、48 和全批次(whole batch)。結果表明,CNN 模型對數據批次并不敏感,因此,數據批次設置為24。

    圖1 CNN 模型框架圖Fig. 1 Architecture of the CNN (convolutional neural network) model

    訓練輸入數據為ERSSTv5 連續(xù)三個月的SSTA,時間范圍為1854~1990 年。原始ERSSTv5網格點由原先的180×89 被重新調整為72×36,目的是避免CNN 模型參數過多導致的過擬合現象。除了限制輸入數據大小外,本文還通過減少神經元數量來盡量避免過擬合現象。如前文所述,CNN模型的神經元數量選取地較少,這本質上也是在限制CNN 模型的參數。CNN 模型的激活函數為Relu 函數,損失函數為常用的平均絕對誤差函數MAE(Mean Absolute Erorr)。

    學習率(Learning Rate,LR)是CNN 模型中相當重要的參數,它代表了模型尋優(yōu)過程的優(yōu)劣。不同CNN 模型架構傾向于不同的LR。因此,對于這27 種CNN 模型架構中的每一種,采用從0.05到0.0001 的LR 范圍進行訓練測試,選取使訓練出的DMI 與觀測值相關系數最高的LR 作為對應CNN模型架構的參數。隨后,在這個模型上,輸入驗證集數據來生成預測的DMI。將27 個不同架構的CNN 模型輸出的平均值用作最終的預測DMI。

    文中針對每個輸入時刻和對應的輸出時刻建立各自的CNN 模型組合(表1)。由于使用了三種預報初始時刻,且分別預報后續(xù)7 個月的DMI,則一共有21 個CNN 模型組合??紤]到不同的27種CNN 模型架構,則總的CNN 模型數量為567(21×27)。同理,預測EIO 和WIO 的CNN 模型的數量也均為567 個。

    表1 CNN 模型、MLP 模型的輸入(SSTA)與輸出(DMI)Table 1 Input (Sea Surface Temperature Anomaly, SSTA) and output (Indian Ocean Dipole Mode Index, DMI) of CNN model and MLP (multi-layer perceptron) model

    MLP 模型是一種前向結構的人工神經網絡(ANN)模型,由輸入層、隱藏層、輸出層組成。MLP 模型所有的神經層均為全聯(lián)接結構,且可以包含多個隱藏層。本文使用的多層感知機MLP 只包含一個隱藏的全聯(lián)接層。與本文所用的CNN 模型相比,MLP 模型缺少了卷積層和池化層,因此提取空間特征的能力相對較弱。為了與CNN 模型的預測結果進行對比,MLP 模型輸入輸出的資料和設置與CNN 一致。

    3 結果

    3.1 IOD 預報技巧

    et al., 2017; Wu and Tang, 2019),利用CNN 模型從JFM、FMA、MAM 三個初始狀態(tài)預報出的DMI 有效預報時效均達到7 個月,而利用MLP 模型所獲得的DMI 有效預報時效則不超過3 個月。同時,利用CNN 模型從三個初始時刻預報出DMI的RMSE 技巧均優(yōu)于MLP 模型。從相關系數和RMSE 這兩種實際預報技巧評價標準來看,CNN模型對DMI 的實際預報技巧遠勝于MLP 模型。這可能是由于CNN 能對海表溫度異常的空間結構進行解析,其卷積與池化能提取SSTA 中的空間特征,并充分考慮SSTA 中的非線性因素,從而能夠將輸入的SSTA 映射到合理的解空間,實現對數據的高效處理。MLP 模型由于不擅長處理結構型數據,對SSTA 空間結構特征的提取能力較弱,因此,MLP 預報的結果較差。現有最先進的動力耦合模式以春季為初始狀態(tài)對DMI 的有效預報時效不到一個季度(Shi et al., 2012),這遠低于CNN 模型的預報技巧,說明CNN 模型在預報DMI 上較現有的動力耦合模式有明顯的優(yōu)勢。RMSE 在預報目標月11 月以后迅速減少,主要是因為IOD 在秋季成熟后迅速衰退所致。RMSE 的大小跟預報量的振幅有直接的關系。

    圖3 為利用CNN 模式預報出的北半球秋季(9 月和10 月)平均的DMI 和對應觀測值隨時間的變化。所有指數均做了標準化處理。我們以一個標準差作為判定IOD 事件的閾值。當秋季平均的DMI 減去多年均值后超過一個標準差時,則認為發(fā)生了正IOD 事件;反之,當小于一個負標準差時,則認為發(fā)生了負IOD 事件。根據此標準,1990~2019 年間共發(fā)生五次正IOD 事件(1994 年、1997 年、2006 年、2015 年、2019 年)和 七 次 負IOD 事件(1990 年、1992 年、1996 年、1998 年、2005 年、2010 年、2016 年)。如果標準化的預報DMI 超過1,則認為預報出了IOD 事件。由于MLP 預報北半球秋季(9 月和10 月)平均DMI 的效果較差,其相關系數不超過0.5(圖2)。因此我們僅重點分析CNN 的結果。

    圖2 1990~2019 年以JFM、FMA、MAM 為起始態(tài),分別使用CNN 模型(實線)與MLP 模型(虛線)預報的DMI 與觀測值的(a)相關系數和(b)均方根誤差(RMSE)Fig. 2 (a) Correlation coefficients and (b) RMSE (root mean square errors) between the forecasted and observed DMI during 1990-2019 using CNN(solid lines) and MLP (dashed lines) models, respectively, with JFM (January-March), FMA (February-April), and MAM (March-May) as the initial conditions

    圖3 1990~2019 年觀測和CNN 模式預報的北半球秋季(9 月、10 月)平均的DMI(標準化的)Fig. 3 Normalized DMI observed and forecasted by CNN model averaged in boreal autumn (September and October) during 1990-2019

    以MAM 為起始態(tài),即提前預報時間為4 個月時,CNN 預報秋季平均DMI 與對應觀測值之間的相關系數為0.72。CNN 模型正確預報出了1997 年和2019 年的強IOD 事件,但低估了1994 年強IOD 事件的強度。對于強負IOD 事件,例如1990 年、1996 年和2016 年,CNN 模型預報的強度與觀測較為符合,但CNN 模型低估了2010 年強負IOD 時間的強度。

    隨著提前時間的增加,預報出的DMI 與觀測的相關系數逐漸減小。以FMA 和JFM 為初始時刻的相關系數預報技巧分別降低至0.67 和0.61。預報出的DMI 在21 世紀初出現較多的誤報。但對于1994 年和1997 年的強正IOD 事件,CNN 模型能夠提前5 個月(以FMA 為初始態(tài))和6 個月(以JFM 為初始態(tài))正確預報出IOD 的強度。而目前動力耦合模式以春季為初始態(tài)對這兩個強事件進行的預報均低估了IOD 事件的強度(Luo et al.,2005; Zhao and Hendon, 2009)。這說明CNN 在預報IOD 事件的強度上比現有的動力耦合模式存在優(yōu)勢。

    3.2 EIO 預報技巧

    印度洋東極子(EIO)在IOD 的發(fā)生發(fā)展中起到重要作用,是IOD 預報的重要區(qū)域(Luo et al.,2005)。圖4 顯示了利用CNN 模型和MLP 模型預報出的DMI 與觀測值的相關系數和RMSE 隨目標月份的變化。MLP 模型以北半球春季為起始態(tài)對EIO 的有效預報能夠提前約3~4 個月做出,而CNN 模型則能夠提前7 個月做出有效預報。同時,CNN 模型預報出EIOI 的RMSE 也明顯低于MLP模式,說明CNN 模型對EIOI 的預報技巧優(yōu)于MLP 模型。目前耦合動力模式以北半球春季為起始態(tài)對EIO 預報的有效預報時效約為4~5 個月(Luo et al., 2007; Zhao and Hendon, 2009)。這一預報技巧優(yōu)于MLP 模型但低于CNN 模型,說明CNN 模型能夠顯著改進目前對EIO 的預報技巧。

    圖4 1990~2019 年以JFM、FMA、MAM 為起始態(tài),分別使用CNN 模型(實線)與MLP 模型(虛線)預報的EIO 指數與觀測值的(a)相關系數和(b)RMSEFig. 4 (a) Correlation coefficients and (b) RMSE between the forecasted and observed EIOI (East Pole Index for Indian Ocean) during 1990-2019 using CNN (solid lines) and MLP (dashed lines) models, respectively, with JFM, FMA, and MAM as the initial conditions

    圖5 是以北半球春季為初始態(tài),CNN 模型預報出的北半球秋季(9 月和10 月)平均的標準化EIO 指數隨時間的變化??梢钥闯?,無論哪個初始態(tài),CNN 模型能夠較好地預報出過去30 年主要的EIO 冷暖事件,例如:1994 年、1997 年的冷事件和1998 年、2010 年的暖事件。但對于2019 年的強冷事件,無論以哪個初始狀態(tài),CNN 模型均低估了其強度。

    圖5 1990~2019 年觀測和CNN 模式預報的北半球秋季(9、10 月)平均EIO 指數(標準化的)Fig. 5 Normalized EIOI index observed and forecasted by CNN mode averaged in boreal autumn (September and October) during 1990-2019

    3.3 WIO 預報技巧

    印度洋西極子(WIO)受到太平洋ENSO 事件的影響,其預報技巧通常比EIO 和DMI 高(Luo et al., 2005, 2007; Zhao and Hendon, 2009; Shi et al., 2012)。圖6 顯示了以北半球春季為初始態(tài),利用CNN 模型和MLP 模型預報出的WIO 指數與觀測的相關系數和RMSE 隨目標月份的變化。MLP 模型對WIO 指數的有效預報時效僅為一個季度不到,這遠低于目前耦合動力模式的有效預報時效(約為6~7 個月,Luo et al., 2007; Zhao and Hendon, 2009)。而CNN 模型對WIO 的有效預報也僅能提前6~7 個月做出,與目前耦合動力模式的預報技巧相近。與EIO 的預報技巧相比,CNN模型對WIO 的預報技巧更低,這與耦合動力模式的預報技巧不一致。這可能是由于我們在構造CNN 模型時僅利用印度洋作為輸入數據,缺少太平洋的ENSO 信號作為先兆因子,從而導致受到ENSO 影響的WIO 的預報技巧降低。

    圖6 1990~2019 年以JFM、FMA、MAM 為起始態(tài),分別使用CNN 模型(實線)與MLP 模型(虛線)預報的WIO 指數與觀測值的(a)相關系數和(b)RMSEFig. 6 (a) Correlation coefficients and (b) RMSE between the forecasted and observed WIOI (West Pole Index for Indian Ocean) during 1990-2019 using CNN and MLP models, respectively, with JFM, FMA, and MAM as the initial conditions

    圖7 給出了以北半球春季為初始態(tài),CNN 模型預報出的北半球秋季(9 月和10 月)平均的標準化WIO 指數隨時間的變化。當提前時間較短時,以MAM 作為初始狀態(tài),CNN 模型基本預報出了主要的WIO 事件,例如,1996 年的冷事件以及2015 年和2019 年的暖事件。隨著提前時間的增長,CNN 模型對WIO 的預報技巧逐漸降低,并且在21 世紀初出現了較多的誤報。對于強WIO 事件(如1996 年、2015 年和2019 年),CNN 模型均能提前5 個月以上(以FMA 為初始條件)預報出。

    圖7 1990~2019 年觀測和CNN 模式預報的北半球秋季(9、10 月)平均WIO 指數(標準化的)Fig. 7 Normalized WIOI observed and forecasted by CNN mode averaged in boreal autumn (September and October) during 1990-2019

    4 結論與討論

    本文使用CNN 模型和MLP 模型處理再分析資料ERSSTv5,并以北半球春季的三個初始態(tài)為起始,對IOD 模態(tài)以及印度洋東西極子的海溫變化進行了預報。結果表明:

    (1)CNN 對DMI、EIOI 和WIOI 的有效預測時效均超過了6 個月。與現在耦合動力模式相比,CNN 模型能夠顯著提升DMI 和EIOI 的預報技巧,但對WIOI 的預報技巧提升有限。

    (2)當預報提前時間為7 個月時,CNN 模型能夠比較準確地預報1994 年、1997 年與2019 年的IOD 事件。

    (3)相比于傳統(tǒng)的神經網絡模型MLP,CNN模型對DMI、EIOI 和WIOI 的預報均有更高的預報技巧。

    IOD 事件在空間結構上的主要特征是東西印度洋的溫度梯度。由于CNN 模型擁有卷積層和池化層,能夠更好地解析圖像數據的空間結構(Zeiler and Fergus, 2014),抓住IOD 事件東西溫度梯度的特征,從而能高效地求解SSTA 與DMI 之間映射關系,因此,CNN 模型對IOD 事件的預報技巧較優(yōu)于傳統(tǒng)的神經網絡MLP 模型。

    CNN 模型對WIOI 的預報不夠理想表明了在訓練模型時,不僅要考慮局地印度洋的信號,而且還要考慮太平洋信號的遙相關作用。今后我們將就這一點進行進一步研究。

    同時,深度學習中超參數的最優(yōu)解會隨具體數據集變化而不同,需根據具體的數據狀況進行逐步調試,應避免過擬合與欠擬合現象發(fā)生。神經元數量、學習率、優(yōu)化器、網絡層數等是模型自身超參數,其調節(jié)一般都是試錯的過程。在算力充足的條件下,可對每個超參數設定一個數值或類型范圍,對它們的組合全部進行訓練。若時間有限,在固定其他超參數為默認值的條件下,可優(yōu)先調節(jié)神經元數量與學習率,同樣能夠取得較為理想的結果。

    總而言之,我們的研究結果表明,CNN 模型在印度洋海溫預報中有著較好的表現,深度學習方法在氣候預報中存在較大的潛力和應用價值。

    猜你喜歡
    北半球印度洋卷積
    北半球最強“星空攝影師”開工啦
    軍事文摘(2023年24期)2023-12-19 06:50:06
    基于3D-Winograd的快速卷積算法設計及FPGA實現
    清涼一夏
    與南亞高壓相聯(lián)的歐亞大陸-印度洋經向環(huán)流
    南北半球天象
    軍事文摘(2019年18期)2019-09-25 08:09:22
    從濾波器理解卷積
    電子制作(2019年11期)2019-07-04 00:34:38
    基于傅里葉域卷積表示的目標跟蹤算法
    去印度洋
    滇池(2018年1期)2018-01-17 21:34:30
    印度洋上接“嫦娥”回家:遠望3號船精確測控探月三期試驗返回器順利再入返回側記
    太空探索(2014年12期)2014-07-12 15:17:10
    一種基于卷積神經網絡的性別識別方法
    電視技術(2014年19期)2014-03-11 15:38:20
    日韩欧美 国产精品| 黄片wwwwww| 又粗又硬又长又爽又黄的视频| 成年女人永久免费观看视频| 亚洲欧美日韩高清专用| 尤物成人国产欧美一区二区三区| 波多野结衣巨乳人妻| 插逼视频在线观看| 国产真实伦视频高清在线观看| 爱豆传媒免费全集在线观看| 男插女下体视频免费在线播放| 热99re8久久精品国产| 日本一本二区三区精品| 免费观看在线日韩| av又黄又爽大尺度在线免费看 | 精品久久久噜噜| 黄色一级大片看看| 神马国产精品三级电影在线观看| АⅤ资源中文在线天堂| 欧美色视频一区免费| 最近最新中文字幕大全电影3| 高清视频免费观看一区二区 | 日韩大片免费观看网站 | 国产成人午夜福利电影在线观看| 男女啪啪激烈高潮av片| 变态另类丝袜制服| 51国产日韩欧美| 国产探花在线观看一区二区| 一级黄色大片毛片| 亚洲欧美成人综合另类久久久 | 人人妻人人澡人人爽人人夜夜 | 色综合亚洲欧美另类图片| 观看美女的网站| 桃色一区二区三区在线观看| 亚洲在久久综合| 只有这里有精品99| 狠狠狠狠99中文字幕| 成人鲁丝片一二三区免费| 国产精品人妻久久久久久| 精品人妻偷拍中文字幕| 久久久精品欧美日韩精品| 午夜激情欧美在线| 国产伦一二天堂av在线观看| 国产不卡一卡二| 99久久无色码亚洲精品果冻| 久99久视频精品免费| 国产单亲对白刺激| 久久亚洲国产成人精品v| 免费无遮挡裸体视频| 亚洲精品成人久久久久久| 18禁在线播放成人免费| 女的被弄到高潮叫床怎么办| 国语自产精品视频在线第100页| 亚洲乱码一区二区免费版| 噜噜噜噜噜久久久久久91| 午夜激情福利司机影院| 尾随美女入室| 中文字幕av成人在线电影| 国产亚洲午夜精品一区二区久久 | 一边亲一边摸免费视频| 久99久视频精品免费| 秋霞伦理黄片| 免费看a级黄色片| 中文字幕熟女人妻在线| 久久久久性生活片| 国产午夜精品论理片| 亚洲成人精品中文字幕电影| 少妇猛男粗大的猛烈进出视频 | 日本欧美国产在线视频| 秋霞伦理黄片| 黄片无遮挡物在线观看| 青春草亚洲视频在线观看| 狂野欧美白嫩少妇大欣赏| 小蜜桃在线观看免费完整版高清| 男的添女的下面高潮视频| 亚洲国产欧美人成| 免费黄色在线免费观看| 国产探花极品一区二区| 看非洲黑人一级黄片| 日日摸夜夜添夜夜爱| 亚洲久久久久久中文字幕| 一个人看的www免费观看视频| av免费观看日本| 成人漫画全彩无遮挡| 女人被狂操c到高潮| 久久久久久久久久成人| 97人妻精品一区二区三区麻豆| 中国国产av一级| 国产美女午夜福利| 亚洲成av人片在线播放无| 国产av在哪里看| 激情 狠狠 欧美| 国产69精品久久久久777片| 久久久久久久久久成人| 在线观看一区二区三区| 午夜精品在线福利| 国产伦理片在线播放av一区| 日韩中字成人| 亚洲最大成人手机在线| 久久久欧美国产精品| 自拍偷自拍亚洲精品老妇| 观看美女的网站| 精品99又大又爽又粗少妇毛片| 国产又色又爽无遮挡免| 午夜a级毛片| videossex国产| 国产高清三级在线| 亚洲av中文av极速乱| av在线天堂中文字幕| 麻豆乱淫一区二区| 国语自产精品视频在线第100页| 色噜噜av男人的天堂激情| 亚洲精品影视一区二区三区av| 国产69精品久久久久777片| 国产免费男女视频| 我的女老师完整版在线观看| 热99re8久久精品国产| 熟妇人妻久久中文字幕3abv| 久久久久久伊人网av| 午夜老司机福利剧场| 国产精品永久免费网站| 1000部很黄的大片| 69人妻影院| 精品久久久久久久久久久久久| 亚洲av二区三区四区| 亚洲av男天堂| 秋霞在线观看毛片| 网址你懂的国产日韩在线| 国产私拍福利视频在线观看| 国产精品国产三级国产专区5o | 美女内射精品一级片tv| 亚洲av中文字字幕乱码综合| www日本黄色视频网| 18禁在线无遮挡免费观看视频| 天堂影院成人在线观看| 午夜日本视频在线| 嫩草影院新地址| 别揉我奶头 嗯啊视频| 久久精品91蜜桃| 精品久久久久久久久av| 欧美成人午夜免费资源| 久久久久久九九精品二区国产| 久久久久久久午夜电影| 午夜激情福利司机影院| 色综合色国产| АⅤ资源中文在线天堂| 99热这里只有精品一区| 国产三级中文精品| 99国产精品一区二区蜜桃av| 国产 一区 欧美 日韩| 久久久亚洲精品成人影院| 天美传媒精品一区二区| 成人一区二区视频在线观看| 日韩欧美精品v在线| 少妇猛男粗大的猛烈进出视频 | 麻豆精品久久久久久蜜桃| 两性午夜刺激爽爽歪歪视频在线观看| 少妇丰满av| 男插女下体视频免费在线播放| 中文字幕久久专区| av播播在线观看一区| 久久韩国三级中文字幕| 亚洲精华国产精华液的使用体验| 3wmmmm亚洲av在线观看| 午夜爱爱视频在线播放| 网址你懂的国产日韩在线| 亚洲国产欧洲综合997久久,| 建设人人有责人人尽责人人享有的 | 亚洲第一区二区三区不卡| 少妇的逼好多水| 国产单亲对白刺激| ponron亚洲| 视频中文字幕在线观看| 夜夜爽夜夜爽视频| 久久久久精品久久久久真实原创| 美女cb高潮喷水在线观看| 国产探花在线观看一区二区| 久久精品久久久久久噜噜老黄 | 久久久久久九九精品二区国产| 亚洲熟妇中文字幕五十中出| a级毛色黄片| 看十八女毛片水多多多| 淫秽高清视频在线观看| 高清毛片免费看| 久久久久久九九精品二区国产| 一级毛片久久久久久久久女| 又爽又黄无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 婷婷六月久久综合丁香| 三级毛片av免费| 看十八女毛片水多多多| 成人漫画全彩无遮挡| 级片在线观看| 国产 一区精品| 久久午夜福利片| 中文字幕av在线有码专区| 男人舔女人下体高潮全视频| 国产精品国产三级国产专区5o | 女人十人毛片免费观看3o分钟| 极品教师在线视频| 国产一级毛片七仙女欲春2| 美女xxoo啪啪120秒动态图| 中文字幕人妻熟人妻熟丝袜美| 国产精品综合久久久久久久免费| 我要搜黄色片| 日韩一区二区视频免费看| 3wmmmm亚洲av在线观看| 国产乱来视频区| 国产精品一区二区三区四区免费观看| 国产国拍精品亚洲av在线观看| 久久精品夜色国产| 最近中文字幕2019免费版| 嫩草影院精品99| 亚洲久久久久久中文字幕| 22中文网久久字幕| 成人午夜高清在线视频| 91精品一卡2卡3卡4卡| 久久久成人免费电影| 国产美女午夜福利| 亚洲中文字幕日韩| 嫩草影院新地址| 国产淫语在线视频| 丰满人妻一区二区三区视频av| 日韩 亚洲 欧美在线| 国产精品1区2区在线观看.| 国产av不卡久久| 特大巨黑吊av在线直播| av在线观看视频网站免费| 一区二区三区免费毛片| 日日摸夜夜添夜夜添av毛片| 天堂√8在线中文| 午夜福利成人在线免费观看| 亚洲欧美精品专区久久| 国产人妻一区二区三区在| 免费无遮挡裸体视频| 九九热线精品视视频播放| 久久精品国产鲁丝片午夜精品| 精品不卡国产一区二区三区| 69人妻影院| 免费人成在线观看视频色| 亚洲欧洲国产日韩| 久久精品国产99精品国产亚洲性色| 国产av码专区亚洲av| 久久精品综合一区二区三区| 亚洲天堂国产精品一区在线| 最新中文字幕久久久久| 成人鲁丝片一二三区免费| 麻豆久久精品国产亚洲av| 亚洲av熟女| 麻豆一二三区av精品| 国产精品久久久久久精品电影| 老司机影院成人| 大香蕉久久网| 日本爱情动作片www.在线观看| 最近最新中文字幕大全电影3| 汤姆久久久久久久影院中文字幕 | 久久久色成人| 三级男女做爰猛烈吃奶摸视频| 九九热线精品视视频播放| 国产亚洲av片在线观看秒播厂 | 久久久久久久久大av| 韩国高清视频一区二区三区| 久久久久久九九精品二区国产| 99国产精品一区二区蜜桃av| 特级一级黄色大片| 久久久久久国产a免费观看| 国产大屁股一区二区在线视频| 国产三级在线视频| 午夜激情福利司机影院| 国产爱豆传媒在线观看| 午夜福利视频1000在线观看| 最近中文字幕2019免费版| 最近中文字幕高清免费大全6| 亚洲在线自拍视频| 高清日韩中文字幕在线| 成年av动漫网址| 国产精品福利在线免费观看| 99热6这里只有精品| 日本免费在线观看一区| 国产精品久久视频播放| 精品酒店卫生间| 中文字幕熟女人妻在线| 看黄色毛片网站| 亚洲精品色激情综合| 免费黄网站久久成人精品| 日产精品乱码卡一卡2卡三| 日本欧美国产在线视频| 久久精品91蜜桃| 性插视频无遮挡在线免费观看| 亚洲人成网站在线观看播放| 少妇的逼好多水| 国产亚洲av片在线观看秒播厂 | 久久久久久久久久黄片| 高清在线视频一区二区三区 | 能在线免费观看的黄片| 九九在线视频观看精品| 久久午夜福利片| 久久鲁丝午夜福利片| 国产精品爽爽va在线观看网站| 精品人妻一区二区三区麻豆| 日韩制服骚丝袜av| 亚洲精品乱码久久久v下载方式| 亚洲精品日韩在线中文字幕| 国产精品美女特级片免费视频播放器| 亚洲国产高清在线一区二区三| 综合色丁香网| av在线亚洲专区| 日本猛色少妇xxxxx猛交久久| 能在线免费观看的黄片| 1000部很黄的大片| 国产淫片久久久久久久久| 熟女电影av网| 能在线免费看毛片的网站| 99久久九九国产精品国产免费| 男女国产视频网站| 成人午夜精彩视频在线观看| 久久精品综合一区二区三区| 波多野结衣巨乳人妻| 亚洲av成人av| 亚洲在线观看片| 亚洲国产精品专区欧美| 日韩欧美国产在线观看| 亚洲av.av天堂| 午夜福利成人在线免费观看| 久久精品久久久久久噜噜老黄 | 国产精品久久久久久精品电影| 桃色一区二区三区在线观看| 亚洲精品,欧美精品| 午夜激情欧美在线| 午夜福利在线在线| 麻豆乱淫一区二区| 亚洲乱码一区二区免费版| 久久人人爽人人片av| 国内少妇人妻偷人精品xxx网站| 天天躁日日操中文字幕| 亚洲最大成人中文| 国产黄片视频在线免费观看| 不卡视频在线观看欧美| 69人妻影院| 特级一级黄色大片| 两个人视频免费观看高清| 精品少妇黑人巨大在线播放 | 网址你懂的国产日韩在线| 一区二区三区四区激情视频| 大香蕉久久网| 精品久久久久久久末码| 亚洲真实伦在线观看| 亚洲欧美一区二区三区国产| 国产三级中文精品| 建设人人有责人人尽责人人享有的 | 美女黄网站色视频| 观看免费一级毛片| 亚洲av免费在线观看| videos熟女内射| 久久久久久伊人网av| 国产av一区在线观看免费| 男女下面进入的视频免费午夜| 日本黄色片子视频| 日韩精品青青久久久久久| 国产精品久久久久久精品电影小说 | 久久精品久久精品一区二区三区| 亚洲人成网站在线观看播放| 成人漫画全彩无遮挡| 日日摸夜夜添夜夜添av毛片| 校园人妻丝袜中文字幕| 国产精品.久久久| 色播亚洲综合网| 精品酒店卫生间| 亚洲国产精品成人综合色| 天天一区二区日本电影三级| 91精品伊人久久大香线蕉| 99久久九九国产精品国产免费| 免费不卡的大黄色大毛片视频在线观看 | 欧美另类亚洲清纯唯美| 久久久久久九九精品二区国产| 偷拍熟女少妇极品色| 午夜久久久久精精品| 亚洲欧美精品专区久久| 大香蕉97超碰在线| 日韩成人伦理影院| 国产精品一区二区三区四区久久| 午夜精品在线福利| 国产精品.久久久| 亚洲精品国产成人久久av| 精品久久久久久成人av| 国产精品一区二区性色av| 男插女下体视频免费在线播放| 夫妻性生交免费视频一级片| 国产亚洲精品久久久com| 日本黄大片高清| 成人av在线播放网站| 国产乱人视频| 啦啦啦啦在线视频资源| www日本黄色视频网| 中文字幕制服av| 久久精品夜夜夜夜夜久久蜜豆| 极品教师在线视频| 亚洲伊人久久精品综合 | 午夜福利网站1000一区二区三区| 好男人在线观看高清免费视频| 中文亚洲av片在线观看爽| 亚洲五月天丁香| 国产av码专区亚洲av| 亚洲最大成人手机在线| 久久久亚洲精品成人影院| 韩国av在线不卡| 日本与韩国留学比较| 一级av片app| 黄片wwwwww| 直男gayav资源| 丰满人妻一区二区三区视频av| 老女人水多毛片| 国产亚洲91精品色在线| 免费av观看视频| 午夜亚洲福利在线播放| 网址你懂的国产日韩在线| 18禁裸乳无遮挡免费网站照片| 99视频精品全部免费 在线| 99久国产av精品国产电影| 91久久精品电影网| 久久久久性生活片| 男人舔女人下体高潮全视频| 乱码一卡2卡4卡精品| 中文精品一卡2卡3卡4更新| 成人一区二区视频在线观看| 久久综合国产亚洲精品| 女人久久www免费人成看片 | 国产女主播在线喷水免费视频网站 | 免费黄色在线免费观看| 哪个播放器可以免费观看大片| 国产视频内射| 免费播放大片免费观看视频在线观看 | 国内精品一区二区在线观看| 一级毛片我不卡| 成人毛片60女人毛片免费| 非洲黑人性xxxx精品又粗又长| 日本免费一区二区三区高清不卡| 久久久久九九精品影院| 国产午夜福利久久久久久| 青春草视频在线免费观看| 麻豆国产97在线/欧美| 国产又色又爽无遮挡免| 亚洲av熟女| 九九在线视频观看精品| 国产黄色视频一区二区在线观看 | 精品久久久久久成人av| 爱豆传媒免费全集在线观看| 免费一级毛片在线播放高清视频| 淫秽高清视频在线观看| 亚洲国产精品合色在线| 99久久精品国产国产毛片| 国产乱人视频| 亚洲av不卡在线观看| 菩萨蛮人人尽说江南好唐韦庄 | av国产免费在线观看| 国产av一区在线观看免费| 国产高清视频在线观看网站| 亚洲国产日韩欧美精品在线观看| 日日干狠狠操夜夜爽| 一级黄色大片毛片| 水蜜桃什么品种好| av免费观看日本| videos熟女内射| 亚洲欧洲日产国产| 人妻制服诱惑在线中文字幕| 99久久精品国产国产毛片| 91午夜精品亚洲一区二区三区| 在线免费观看不下载黄p国产| 亚洲精品,欧美精品| 九九爱精品视频在线观看| 日韩制服骚丝袜av| a级一级毛片免费在线观看| 搡老妇女老女人老熟妇| 成人亚洲欧美一区二区av| av女优亚洲男人天堂| 亚洲欧美精品专区久久| 精品久久久久久久久久久久久| 五月玫瑰六月丁香| 国产伦一二天堂av在线观看| 欧美成人一区二区免费高清观看| 人妻夜夜爽99麻豆av| 成年免费大片在线观看| 久久这里只有精品中国| 秋霞伦理黄片| 99国产精品一区二区蜜桃av| av在线蜜桃| 日本黄色片子视频| av免费在线看不卡| 亚洲成av人片在线播放无| 毛片一级片免费看久久久久| 麻豆久久精品国产亚洲av| 美女高潮的动态| 亚洲综合色惰| 欧美高清成人免费视频www| 亚洲欧美日韩卡通动漫| 人妻夜夜爽99麻豆av| 99热这里只有是精品在线观看| 亚洲乱码一区二区免费版| 亚洲精品一区蜜桃| 国产精品日韩av在线免费观看| 国产成人一区二区在线| 一区二区三区乱码不卡18| 国产乱人视频| 亚洲不卡免费看| 天堂√8在线中文| 亚洲,欧美,日韩| 成人无遮挡网站| 在线观看av片永久免费下载| 国产91av在线免费观看| 精品久久国产蜜桃| 白带黄色成豆腐渣| 精品欧美国产一区二区三| 日日啪夜夜撸| 国产在线一区二区三区精 | 老师上课跳d突然被开到最大视频| 天堂√8在线中文| 一个人看视频在线观看www免费| 女人十人毛片免费观看3o分钟| 美女脱内裤让男人舔精品视频| 亚洲国产精品成人久久小说| 精品久久久久久久久亚洲| 一级爰片在线观看| 精品免费久久久久久久清纯| 亚洲中文字幕一区二区三区有码在线看| 免费黄色在线免费观看| 亚洲伊人久久精品综合 | 在线免费观看的www视频| 波多野结衣巨乳人妻| 性插视频无遮挡在线免费观看| 麻豆久久精品国产亚洲av| 一级黄色大片毛片| 人妻系列 视频| 99久久中文字幕三级久久日本| 精品久久久久久久人妻蜜臀av| 青青草视频在线视频观看| 成年女人看的毛片在线观看| 国产精品伦人一区二区| 久久人人爽人人爽人人片va| 欧美激情久久久久久爽电影| 亚洲欧美中文字幕日韩二区| 国内揄拍国产精品人妻在线| a级一级毛片免费在线观看| av卡一久久| 中文字幕av成人在线电影| av线在线观看网站| 国产精品爽爽va在线观看网站| 少妇高潮的动态图| 久久精品夜色国产| 久久久久久久久大av| 国产伦理片在线播放av一区| 美女内射精品一级片tv| 久久久亚洲精品成人影院| 色尼玛亚洲综合影院| 中文字幕熟女人妻在线| 日本黄色视频三级网站网址| 国产精品伦人一区二区| 亚洲内射少妇av| 白带黄色成豆腐渣| 久久人人爽人人片av| 99久久中文字幕三级久久日本| 国产又色又爽无遮挡免| 国产在视频线精品| 免费观看a级毛片全部| 真实男女啪啪啪动态图| 久久99热这里只频精品6学生 | 亚洲国产精品成人综合色| 国内少妇人妻偷人精品xxx网站| 如何舔出高潮| 久久精品综合一区二区三区| 小说图片视频综合网站| 超碰97精品在线观看| 99在线人妻在线中文字幕| 看免费成人av毛片| 三级毛片av免费| 免费电影在线观看免费观看| 亚洲,欧美,日韩| 久久久国产成人免费| 久久精品久久久久久噜噜老黄 | 久久6这里有精品| 中文天堂在线官网| 内地一区二区视频在线| 综合色av麻豆| 日本-黄色视频高清免费观看| 三级国产精品片| 亚州av有码| 九草在线视频观看| 久久人人爽人人片av| 亚洲精品色激情综合| www.av在线官网国产| 亚洲成人久久爱视频| 国产精品一及| 看免费成人av毛片| 非洲黑人性xxxx精品又粗又长| 看非洲黑人一级黄片| 国产日韩欧美在线精品| 久久韩国三级中文字幕| 国产亚洲5aaaaa淫片| 成人综合一区亚洲| 精品久久久久久久末码| av播播在线观看一区| 成人欧美大片| 久久久久九九精品影院| 国产一区二区在线观看日韩| 18+在线观看网站| 亚洲av免费高清在线观看| 搡女人真爽免费视频火全软件| 看片在线看免费视频| av卡一久久| 日韩欧美精品v在线| 国产激情偷乱视频一区二区| 国产精品久久久久久av不卡| 亚洲欧美精品专区久久| 99久久精品一区二区三区| 男插女下体视频免费在线播放| 深夜a级毛片|