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

    海河平原夏玉米主要生育期發(fā)生高溫干旱并發(fā)事件的氣候學分析*

    2021-06-29 06:02:40韓佳昊王麗榮楊再強
    中國農業(yè)氣象 2021年6期
    關鍵詞:全生育期夏玉米生育期

    韓佳昊,張 琪**,王麗榮,楊再強

    (1.南京信息工程大學應用氣象學院,南京 210044;2.河北省氣象災害防御中心,石家莊 050021;3.南京信息工程大學氣象災害預報預警與評估協(xié)同創(chuàng)新中心/江蘇省農業(yè)氣象重點實驗室,南京 210044)

    高溫干旱并發(fā)事件的危害程度要明顯高于單一災害事件[1-2]。Suzuki等[3]指出,美國1980?2012年單一氣候災害損失在10億美元以上,而高溫干旱復合災害造成的損失更為嚴重;21世紀初高溫干旱復合災害導致美國樹木大面積死亡[4];2003年歐洲高溫干旱復合災害給當地社會、農業(yè)生產等造成巨大影響[5];2006年,中國川渝地區(qū)受高溫干旱復合災害影響,直接經濟損失高達216.4億元[6]。玉米作為主要的糧食作物之一,易受到氣溫和降水變化的影響,華北平原玉米種植面積占全國總播種面積的30%,卻產出了全國玉米總產量的50%以上[7],同時華北平原也是中國對氣候變化響應敏感的地帶,受氣候變暖的影響,高溫干旱復合脅迫給夏玉米生產帶來極大的不穩(wěn)定性,已成為威脅糧食安全的重要因素之一[8]。因此,研究華北地區(qū)氣候變化背景下高溫干旱并發(fā)事件時空演變規(guī)律,特別是夏玉米主要生育期的并發(fā)事件發(fā)生規(guī)律具有重要的意義。

    關于華北地區(qū)高溫或干旱單一氣象事件的時空演變規(guī)律的研究結果已較為豐富。研究顯示,全球變暖背景下,華北地區(qū)增溫明顯,是中國增溫最快的地區(qū)之一,作物生育期積溫有所增加但同時極端氣溫事件顯著增多,特別是夏玉米花期高溫熱害發(fā)生頻率明顯增加[9-13]。關于干旱,研究顯示華北地區(qū)自20世紀80年代降水明顯減少,干旱風險增加;夏玉米各生長階段中,抽雄?成熟階段干旱發(fā)生范圍最廣[14-16]。雖然當前國內有較豐富的夏玉米各生育期高溫、干旱單一氣象事件時空演變規(guī)律的研究,但鮮有高溫干旱并發(fā)事件的研究報道。國外學者對高溫干旱并發(fā)事件有一定的研究,如Mazdiyasni等[17]采用累積分布函數(Cumulative Distribution Function,CDFs)研究發(fā)現(xiàn),1990?2010年美國高溫干旱并發(fā)事件發(fā)生風險明顯高于1960?1980年,特別是在美國南部地區(qū);Sharma等[18]也同樣采用累積分布函數發(fā)現(xiàn),印度1981?2010年并發(fā)事件的發(fā)生風險顯著高于1951?1980年,尤其在印度南部地區(qū)。但已有研究較少涉及作物主要生育期內的高溫干旱并發(fā)事件發(fā)生規(guī)律的研究。

    針對當前高溫干旱并發(fā)事件,特別是夏玉米主要生育期并發(fā)事件的時空規(guī)律研究并不充分,本研究選取華北地區(qū)的河海平原為案例研究區(qū),利用研究區(qū)27個氣象站點1960?2019年逐日最高氣溫、降水數據,夏玉米生育期數據,研究夏玉米全生育期及各主要生育期內不同閾值水平的高溫干旱并發(fā)事件時空分布規(guī)律,以期對農業(yè)布局規(guī)劃提供一定理論依據。

    1 材料與方法

    1.1 研究區(qū)概況

    圖1 研究區(qū)及站點分布Fig.1 Study area and meteorological stations

    1.2 數據來源

    利用1960?2019年海河平原27個氣象站逐日最高氣溫和降水量數據識別高溫干旱并發(fā)事件,站點分布如圖1所示。氣象數據源于中國地面氣候資料日值數據集(V3.0),該數據經過嚴格的質量控制和檢查。由于不同站點、不同年份夏玉米生育期不同,主要參考已有文獻中統(tǒng)計的各站點多年平均夏玉米播種日期、拔節(jié)日期、開花和成熟日期(表1)[8]。調查結果顯示,海河平原夏玉米播種期集中在6月1?18日,拔節(jié)期在7月7?24日,開花期為8月1?14日,成熟期為9月17日?10月4日。夏玉米全生育期、播種?拔節(jié)期、拔節(jié)?開花期和開花?成熟期的平均長度分別為104、35、22和48d。

    表1 各站點夏玉米主要生育期多年平均起始時間(月-日)Table 1 Average starting date of main growth periods of summer maize at each station (mm-dd)

    1.3 方法

    1.3.1 高溫干旱并發(fā)事件的識別

    采用標準化降水指數(Standardized Precipitation Index, SPI)識別夏玉米生育期是否發(fā)生干旱事件,將各站點生育期長度作為時間尺度,統(tǒng)計逐年生育期內總降水量,計算SPI,SPI<?0.5為發(fā)生干旱事件[20]。溫度越高、持續(xù)時間越長的高溫事件其影響越嚴重,因此,本研究在確定高溫事件時從溫度高低和持續(xù)時間兩個角度定義,采用百分位法確定溫度高低,研究85百分位和90百分位兩個等級,持續(xù)時間分為3d和5d兩個等級,共4種閾值水平的高溫事件,分別為低高溫閾值,即日最高氣溫超過85百分位且連續(xù)3d及以上確定為一次高溫事件,記為T85%&3d;中高溫閾值,即日最高氣溫超過85百分位且連續(xù)5d及以上確定為一次高溫事件,記為T85%&5d,或日最高氣溫超過90百分位且連續(xù)3d及以上確定為一次高溫事件,記為T90%&3d;高高溫閾值,即日最高氣溫超過90百分位且連續(xù)5d及以上確定為一次高溫事件,記為T90%&5d。只要某生育期出現(xiàn)一次高溫事件就確定為該生育期發(fā)生高溫。若生育期同時發(fā)生了干旱和高溫事件則定義該生育期發(fā)生了高溫干旱并發(fā)事件。

    采用上述方法確定研究區(qū)各站點近60a夏玉米全生育期及各主要生育期(播種?拔節(jié)、拔節(jié)?開花和開花?成熟)在4種高溫閾值水平下是否發(fā)生了并發(fā)事件。將每年發(fā)生高溫干旱并發(fā)事件站點數與總站點數(27個)之比定義為“高溫干旱并發(fā)事件的發(fā)生范圍”,將站點近60a發(fā)生高溫干旱并發(fā)事件的年數與總年數之比定義為“高溫干旱并發(fā)事件的發(fā)生頻率”。

    1.3.2 Cramér-von Mises(CvM)突變檢驗法

    Cramér-von Mises(CvM)突變檢驗法主要用于兩個獨立同分布變量時間序列的差異性檢驗,也可以應用于單變量不同時間序列對比及突變點分析,具體方法是通過比較兩個樣本序列的累積分布函數平均差異水平來計算檢驗統(tǒng)計量[21-23]。CvM突變檢驗方法對時間序列中存在的突變更加敏感,能夠檢驗出其它方法如MK突變檢驗所不能識別出的突變點[17]。具體過程為

    推薦理由:《小老鼠又上燈臺嘍》是“九神鹿繪本館”的一本新作。文字作者袁曉峰在傳統(tǒng)兒歌的基礎上,圍繞“閱讀”這個主題,結合多年閱讀推廣經驗及兒童閱讀特點進行了全新創(chuàng)作,讀起來朗朗上口,好玩風趣。在文字的基礎上,繪圖作者趙曉音用色彩鮮明、稚氣脫俗的表現(xiàn)手法,將老鼠、小貓、大貓、小狗、孩子生動形象地呈現(xiàn)出來,讓繪本充滿了節(jié)奏感和電影的動感。

    式中,I為指標函數,xτ為時間為τ的對應的值[24]。

    在Cramér-von Mises突變檢測方法中,對于所有τ值,要保證兩個子樣本有相等的方差和平均值,使用Anderson等[25]提出的標準化方法實現(xiàn),即

    由式(1)逐步計算CvM統(tǒng)計量Wτ,n,在所有檢驗統(tǒng)計量中挑選最大的Wτ,n所對應的時間τ即為差異最大的點,表明其前后時段內高溫干旱并發(fā)事件序列差異最為明顯,將其確定為突變點所在年份。計算過程中將統(tǒng)計量Wτ,n與檢驗統(tǒng)計表相對比[25],確定突變點的顯著性水平。

    1.3.3 累積分布函數

    采用累積分布函數(CDFs)比較時間序列突變前后數據分布情況,進一步利用雙樣本KS檢驗(Kolmogorov–Smirnov test)確定分布差異的顯著性水平,KS檢驗的統(tǒng)計量P<0.05說明兩個分布存在顯著差異[8,17]。

    2 結果與分析

    2.1 高溫干旱并發(fā)事件發(fā)生范圍的變化趨勢

    由圖2可見,在整個研究期內(1960?2019年),海河平原27個氣象站中,夏玉米全生育期各個閾值水平的高溫干旱并發(fā)事件的平均發(fā)生范圍為0.21(即有21%的站點發(fā)生了高溫干旱并發(fā)事件),其中,低高溫閾值下(T85%&3d)并發(fā)干旱事件的發(fā)生范圍最大,多年平均值為0.29;中高溫閾值下(T85%&5d和T90%&3d)并發(fā)干旱事件的發(fā)生范圍居中,平均值為0.22;高高溫閾值下(T90%&5d)并發(fā)干旱事件的發(fā)生范圍最小,平均值為0.12。3個夏玉米主要生育期中,開花?成熟期并發(fā)事件的發(fā)生范圍最廣,各閾值水平的多年平均值為0.15;播種?拔節(jié)期居中,平均為0.10;拔節(jié)?開花期最小,均值為0.08??梢?,在海河平原夏玉米生長過程中,開花?成熟期容易發(fā)生大范圍的高溫干旱并發(fā)事件。從高溫干旱并發(fā)事件發(fā)生范圍數據的年際變化趨勢看,全生育期以及各主要生育期線性趨勢均未通過0.05水平的顯著性檢驗,說明并發(fā)事件的發(fā)生范圍無明顯長期變化趨勢。

    圖2 1960?2019年夏玉米各生育期不同高溫閾值水平下并發(fā)干旱事件發(fā)生范圍(站點比)的年際變化Fig.2 Interannual variability of occurred scope of concurrent drought events (station ratio) under different extreme heat threshold levels in various growth periods of summer maize from 1960 to 2019

    2.2 高溫干旱并發(fā)事件發(fā)生范圍時間序列突變特征

    2.2.1 突變點識別

    時間序列存在突變點是造成其長期趨勢不顯著的原因之一。采用CvM方法識別海河平原夏玉米生育期高溫干旱并發(fā)事件發(fā)生范圍年際變化是否存在突變點(圖3)。從圖3可知,全生育期的低高溫閾值下(T85%&3d)并發(fā)干旱事件的發(fā)生范圍于1979年附近發(fā)生突變,并達到顯著水平(圖3a),其余3種閾值水平的高溫干旱并發(fā)事件的發(fā)生范圍均在1997年附近存在顯著突變點(圖3b?圖3d);在播種?拔節(jié)期,各閾值水平的高溫干旱并發(fā)事件的發(fā)生范圍均在1997年附近發(fā)生顯著突變(圖3e?圖3h);在拔節(jié)?開花期,僅T85%&3d和T90%&3d這兩種高溫閾值水平下的并發(fā)干旱事件的發(fā)生范圍存在顯著突變點,在1997年附近(圖3i、圖3k);在開花?成熟期,各閾值的高溫干旱并發(fā)事件的發(fā)生范圍均在1991年附近存在顯著突變點(圖3m?圖3p)??傮w來說,海河平原近60a來夏玉米各生育期內高溫干旱并發(fā)事件的發(fā)生范圍在20世紀90年代存在顯著突變,可以用來解釋其長期變化趨勢不顯著。

    圖3 基于CvM的各生育期不同高溫閾值水平下并發(fā)干旱事件發(fā)生范圍時間序列的突變檢驗Fig.3 Mutation test of the time series of the occurred scope of concurrent drought events under different extreme heat threshold levels in the various growth periods based on Cramér-von Mises (CvM) method

    2.2.2 突變檢驗

    可見,雖然各種情景下的高溫干旱并發(fā)事件發(fā)生范圍時間序列的突變點略有不同,但基本都在20世紀90年代,因此,采用CDFs展示各種情景下1960?1989(突變前)和2000?2019(突變后)這兩個時段內并發(fā)事件發(fā)生范圍數據的累積概率分布情況,以檢驗突變前后是否有顯著差異,結果如圖4。由圖可見,各個生育期、各閾值水平的高溫干旱并發(fā)事件發(fā)生范圍的CDFs曲線均為突變后的曲線在突變前的右側,說明突變后并發(fā)事件發(fā)生范圍均高于突變前。全生育期突變后曲線偏右幅度最大,說明發(fā)生范圍較突變前增加最為明顯;播種?拔節(jié)、拔節(jié)?開花與開花?成熟這3個生育期突變后曲線右偏幅度較全生育期小,但開花?成熟期偏移幅度略大于其它兩個生育期,說明開花?成熟期并發(fā)事件發(fā)生范圍自突變后增加幅度最大。對于不同閾值水平的高溫干旱并發(fā)事件,高高溫閾值下(T90%&5d)并發(fā)干旱事件的發(fā)生范圍突變后的曲線右偏最明顯,說明并發(fā)事件發(fā)生范圍增加幅度最大。采用KS檢驗方法研究上述突變后兩個時段內并發(fā)事件發(fā)生范圍數據的累積分布函數的差異,結果表明全生育期及高閾值水平(T90%&5d)的差異性都達到顯著水平(P<0.05),即突變前后兩個時段有顯著差異,檢驗了前面突變點確定的合理性。

    圖4 各生育期不同閾值下并發(fā)事件發(fā)生范圍時間序列突變前后CDFs曲線對比圖Fig.4 Comparison of CDFs curves before and after mutation of the occurred scope of concurrent events under different threshold levels at each growth periods

    2.3 高溫干旱并發(fā)事件發(fā)生頻率及其變化特征

    2.3.1 各站點的發(fā)生頻率

    由圖5可見,在整個分析期內,海河平原27個氣象站中,夏玉米全生育期各個閾值水平的高溫干旱并發(fā)事件的平均發(fā)生頻率達到5a一遇以上(≥0.20),其中,低高溫閾值下(T85%&3d)并發(fā)干旱事件的發(fā)生頻率最大,基本所有地區(qū)為4a一遇以上(≥0.25);中高溫閾值下(T85%&5d和T90%&3d)并發(fā)干旱事件的發(fā)生頻率居中,多數地區(qū)發(fā)生頻率在7a一遇至3a一遇(0.14~0.33);高高溫閾值下(T90%&5d)并發(fā)干旱事件的發(fā)生頻率最小,多集中在6a一遇以下(≤0.17)。3個夏玉米主要生育期中,開花?成熟期各站點并發(fā)事件的發(fā)生頻率最大,各閾值水平的并發(fā)事件發(fā)生頻率多集中在10a一遇至3a一遇(0.10~0.33);播種?拔節(jié)期與拔節(jié)?成熟期并發(fā)事件發(fā)生頻率相似,多集中在6a一遇以下(≤0.17)。對于兩種中高溫閾值(T85%&5d和T90%&3d)并發(fā)干旱事件,各生育期T90%&3d的發(fā)生頻率總體高于T85%&5d,說明不太容易發(fā)生持續(xù)事件較長的高溫并發(fā)干旱事件。從空間分布看,各生育期不同閾值水平的高溫干旱并發(fā)事件發(fā)生頻率總體上由西北部向東南部逐漸減少。

    圖5 各生育期不同高溫閾值水平下并發(fā)干旱事件發(fā)生頻率的空間分布Fig.5 Spatial distribution of the frequency of concurrent drought events under different extreme heat threshold levels in the various growth periods

    2.3.2 突變前后發(fā)生頻率的變化

    統(tǒng)計各站點突變后與突變前高溫干旱并發(fā)事件發(fā)生頻率,圖6為突變后發(fā)生頻率與突變前頻率差值的空間分布。由圖可見,在夏玉米全生育期,低高溫閾值下(T85%&3d)并發(fā)干旱事件的發(fā)生頻率有14個站點(站次比:0.52)發(fā)生頻率差值為正,即有52%的站點突變后并發(fā)事件發(fā)生頻率較突變前增加,這14個站點發(fā)生頻率增加量平均為0.10;高高溫閾值下(T90%&5d)并發(fā)干旱事件的發(fā)生頻率突變后較突變前增加的站點有21個(站次比0.78),發(fā)生頻率增加量平均為0.15。3個夏玉米主要生育期中,開花?成熟期并發(fā)事件的發(fā)生頻率在突變后為增加的站點最多,各閾值水平下平均有19個站點(站次比0.69),發(fā)生頻率平均增加0.11;播種?拔節(jié)期居中,平均約有17個站點(站次比0.63)增加,發(fā)生頻率平均增加0.09;拔節(jié)?開花期最少,平均約有14個站點(站次比:0.52)增加,發(fā)生頻率平均增加0.08。各生育期各閾值水平下高溫干旱并發(fā)事件在突變后增多的站點主要分布在研究區(qū)北部和東部地區(qū),而中南部很多站點并發(fā)事件呈現(xiàn)減少,特別是對于拔節(jié)?開花期及全生育期的低閾值水平下的并發(fā)事件呈現(xiàn)減少的站點更多。

    圖6 各生育期不同閾值水平高溫干旱并發(fā)事件突變后與突變前發(fā)生頻率差值的空間分布Fig.6 Spatial distribution of the difference between the frequency before and after mutation under different extreme heat threshold levels in the various growth periods

    3 討論與結論

    3.1 結論

    (1)海河平原近60a各生育期不同閾值水平的高溫干旱并發(fā)事件發(fā)生范圍均無顯著長期年際變化趨勢,在20世紀90年代存在顯著突變,突變后并發(fā)事件的發(fā)生范圍增加明顯。夏玉米各主要生育期中,在開花?成熟期的高溫干旱并發(fā)事件的發(fā)生范圍最大,且這一發(fā)育期并發(fā)事件發(fā)生范圍在突變后增加最為明顯;各閾值水平下,高高溫閾值下(T90%&5d)并發(fā)干旱事件的發(fā)生范圍雖然較小但突變后增加最為明顯。

    (2)近60a研究區(qū)夏玉米各生育期不同閾值水平的高溫干旱并發(fā)事件發(fā)生頻率在空間上總體呈現(xiàn)出西北部較高向東南減少,且并發(fā)事件發(fā)生頻率較高的地區(qū)同時也是在突變后發(fā)生頻率增多的地區(qū)。研究區(qū)中南部很多站點并發(fā)事件呈現(xiàn)減少趨勢,尤以拔節(jié)?開花期及全生育期低閾值水平下的并發(fā)事件呈現(xiàn)減少的站點更多。

    3.2 討論

    已有研究顯示,全球氣候在20世紀90年代出現(xiàn)明顯突變,主要是由太陽能量的輸出時間變短、火山活動增加、重大厄爾尼諾事件等多方面影響所帶來的[26?27];前人在研究美國高溫干旱并發(fā)事件時也發(fā)現(xiàn)在90年代存在明顯突變[17],本研究得出的海河平原地區(qū)夏玉米生育期的高溫干旱并發(fā)事件突變時間與世界其它地區(qū)相似。90年代以來,中國極端高溫事件顯著增多[28],同時Aghakouchak等[29]認為極端高溫事件增加大大提高了高溫干旱并發(fā)事件發(fā)生頻率,與本研究多數站點自突變后高溫干旱并發(fā)事件發(fā)生頻率呈現(xiàn)增多的結果相一致。

    和驊蕓等[12]在對華北地區(qū)夏玉米各生育期中高溫事件的研究中指出,夏玉米花期高溫發(fā)生頻率最高,萬能涵等[16]對華北地區(qū)夏玉米各生育期干旱事件的研究顯示,在夏玉米花期附近干旱事件發(fā)生較為頻繁,兩種事件都在這個時期頻發(fā),勢必會導致并發(fā)事件頻發(fā)。夏玉米花期是對高溫、干旱最敏感的時期,一旦遭遇高溫、干旱天氣,夏玉米產量會受到極大影響[30]。本研究顯示海河平原地區(qū)夏玉米開花?成熟期高溫干旱并發(fā)事件發(fā)生頻率最大、范圍最廣,并且也是增加幅度最大的階段,因此,當地夏玉米生產面臨較大的挑戰(zhàn),亟需通過調整種植制度、提高田間管理水平等方式降低可能的影響。

    猜你喜歡
    全生育期夏玉米生育期
    小麥收割之后 如何種植夏玉米才能高產
    夏玉米高產的關鍵栽培技術措施
    湖南省2016年審定通過的水稻新品種(下)
    中國稻米(2017年2期)2017-04-28 08:00:06
    2016年靖遠縣春玉米全生育期農業(yè)氣象條件分析
    湖南省2015年審定通過的水稻新品種(4)
    中國稻米(2016年2期)2016-06-29 09:53:29
    基于作物生育期的潛在蒸散的時空演變特征及R/S 分析
    西藏科技(2015年5期)2015-09-26 11:55:25
    黃驊市2013年度冬小麥全生育期農業(yè)氣象條件分析
    2013-2014年度二二二團冬小麥各生育期氣象條件分析
    阿拉爾市2012年棉花生育期氣象條件分析
    棉花生育期調控關鍵技術
    欧美日韩一级在线毛片| 亚洲综合色网址| 成人手机av| 中文欧美无线码| 一区在线观看完整版| 国产精品一国产av| 欧美亚洲 丝袜 人妻 在线| 1024视频免费在线观看| 电影成人av| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品av麻豆av| 在现免费观看毛片| 国产色婷婷99| 亚洲av日韩在线播放| 一级毛片 在线播放| 中文字幕高清在线视频| 国产欧美日韩一区二区三区在线| 亚洲,欧美精品.| 交换朋友夫妻互换小说| 自拍欧美九色日韩亚洲蝌蚪91| 9191精品国产免费久久| 99久久人妻综合| 国产欧美亚洲国产| 在线观看免费日韩欧美大片| 免费久久久久久久精品成人欧美视频| 国产精品久久久久久精品电影小说| 新久久久久国产一级毛片| 日韩熟女老妇一区二区性免费视频| 午夜av观看不卡| 成人18禁高潮啪啪吃奶动态图| 水蜜桃什么品种好| 成人午夜精彩视频在线观看| 久久久久久久久久久久大奶| www.精华液| 国产精品熟女久久久久浪| 精品亚洲成国产av| 国产高清不卡午夜福利| 日日撸夜夜添| 精品久久久精品久久久| 久久久国产欧美日韩av| 肉色欧美久久久久久久蜜桃| 成人国语在线视频| 欧美人与性动交α欧美软件| 午夜影院在线不卡| 大话2 男鬼变身卡| 免费久久久久久久精品成人欧美视频| av女优亚洲男人天堂| 美女中出高潮动态图| 亚洲av日韩在线播放| 国产精品人妻久久久影院| 国产在视频线精品| 久久国产精品男人的天堂亚洲| 肉色欧美久久久久久久蜜桃| 啦啦啦在线观看免费高清www| 欧美人与性动交α欧美软件| 美女中出高潮动态图| 纵有疾风起免费观看全集完整版| 色综合欧美亚洲国产小说| 我要看黄色一级片免费的| 午夜福利视频在线观看免费| 制服丝袜香蕉在线| 日韩av免费高清视频| 国产片特级美女逼逼视频| 国产xxxxx性猛交| 97人妻天天添夜夜摸| 国产精品麻豆人妻色哟哟久久| 最近中文字幕2019免费版| 街头女战士在线观看网站| 欧美精品一区二区免费开放| 蜜桃国产av成人99| 亚洲成人免费av在线播放| kizo精华| 国产高清不卡午夜福利| 国产女主播在线喷水免费视频网站| 成人国语在线视频| av国产久精品久网站免费入址| 欧美日韩av久久| 啦啦啦啦在线视频资源| 亚洲欧美色中文字幕在线| 午夜福利,免费看| 中文字幕另类日韩欧美亚洲嫩草| 两性夫妻黄色片| 午夜av观看不卡| 精品久久久久久电影网| 午夜激情久久久久久久| 少妇猛男粗大的猛烈进出视频| 日韩中文字幕欧美一区二区 | 啦啦啦 在线观看视频| 亚洲成人免费av在线播放| 亚洲情色 制服丝袜| 9色porny在线观看| 一本久久精品| 一二三四中文在线观看免费高清| 久久久久久久精品精品| av.在线天堂| 国产有黄有色有爽视频| 亚洲国产欧美一区二区综合| 视频区图区小说| 久久ye,这里只有精品| 亚洲情色 制服丝袜| 美女福利国产在线| 欧美少妇被猛烈插入视频| 成人漫画全彩无遮挡| 精品少妇内射三级| 99九九在线精品视频| 一边摸一边做爽爽视频免费| 日韩一区二区视频免费看| 啦啦啦 在线观看视频| 午夜老司机福利片| 热99国产精品久久久久久7| 日韩 欧美 亚洲 中文字幕| av网站免费在线观看视频| 丰满乱子伦码专区| 国产成人欧美| 国产av精品麻豆| 中文字幕色久视频| 久久久久国产精品人妻一区二区| 国产探花极品一区二区| 女人精品久久久久毛片| 曰老女人黄片| 欧美人与善性xxx| 欧美日韩一级在线毛片| 欧美人与性动交α欧美软件| 热99久久久久精品小说推荐| 日本色播在线视频| 亚洲国产最新在线播放| 国产淫语在线视频| 又大又爽又粗| 久久99热这里只频精品6学生| 狠狠精品人妻久久久久久综合| 午夜福利在线免费观看网站| 亚洲av福利一区| 国产精品免费视频内射| 亚洲国产日韩一区二区| 国产97色在线日韩免费| 国产 精品1| 欧美黄色片欧美黄色片| 久久这里只有精品19| 人妻人人澡人人爽人人| 久久久欧美国产精品| 91国产中文字幕| 91精品伊人久久大香线蕉| 欧美精品高潮呻吟av久久| 男的添女的下面高潮视频| 久久影院123| 黄色 视频免费看| 成人亚洲欧美一区二区av| 色综合欧美亚洲国产小说| 国产乱来视频区| 九草在线视频观看| 可以免费在线观看a视频的电影网站 | 黄色怎么调成土黄色| netflix在线观看网站| 国产成人欧美在线观看 | 欧美在线一区亚洲| 国产探花极品一区二区| 中文欧美无线码| 多毛熟女@视频| 一二三四中文在线观看免费高清| 亚洲欧美中文字幕日韩二区| 黑人巨大精品欧美一区二区蜜桃| 久久毛片免费看一区二区三区| 国产福利在线免费观看视频| 亚洲免费av在线视频| 黄色视频在线播放观看不卡| 国产精品.久久久| 亚洲成国产人片在线观看| 国产精品99久久99久久久不卡 | 国产成人免费无遮挡视频| 国产精品熟女久久久久浪| 2018国产大陆天天弄谢| 免费观看性生交大片5| 亚洲自偷自拍图片 自拍| 久久久久久人妻| 国产野战对白在线观看| 不卡视频在线观看欧美| 纵有疾风起免费观看全集完整版| 日韩制服骚丝袜av| xxxhd国产人妻xxx| 一级黄片播放器| 秋霞伦理黄片| av网站免费在线观看视频| 观看美女的网站| 亚洲少妇的诱惑av| 日韩大码丰满熟妇| 最近手机中文字幕大全| 91精品三级在线观看| 日韩大片免费观看网站| 国产精品一二三区在线看| 免费高清在线观看日韩| 久久久久人妻精品一区果冻| 精品亚洲成国产av| 男人操女人黄网站| 在线 av 中文字幕| 亚洲久久久国产精品| 桃花免费在线播放| 天堂俺去俺来也www色官网| 人人妻人人澡人人爽人人夜夜| 2018国产大陆天天弄谢| 男人舔女人的私密视频| 伊人亚洲综合成人网| 久久97久久精品| 巨乳人妻的诱惑在线观看| 亚洲色图 男人天堂 中文字幕| 国产精品久久久久成人av| 宅男免费午夜| 国产一区亚洲一区在线观看| 三上悠亚av全集在线观看| 制服丝袜香蕉在线| 国产1区2区3区精品| 国产免费又黄又爽又色| 国产极品粉嫩免费观看在线| 久热爱精品视频在线9| 久久久久久久久久久免费av| 搡老乐熟女国产| 伊人久久国产一区二区| 午夜日本视频在线| 婷婷色综合www| 久久亚洲国产成人精品v| 一级毛片 在线播放| 黄网站色视频无遮挡免费观看| 亚洲av电影在线进入| 国产深夜福利视频在线观看| 欧美亚洲 丝袜 人妻 在线| 一区二区三区乱码不卡18| 亚洲精华国产精华液的使用体验| 中文字幕色久视频| 成年av动漫网址| 伊人亚洲综合成人网| 国产1区2区3区精品| 考比视频在线观看| 成年女人毛片免费观看观看9 | 国产精品熟女久久久久浪| 不卡视频在线观看欧美| 日本vs欧美在线观看视频| 亚洲国产最新在线播放| 亚洲精品美女久久久久99蜜臀 | 女性被躁到高潮视频| 黄色视频不卡| 精品第一国产精品| 精品视频人人做人人爽| 一级毛片 在线播放| 亚洲av综合色区一区| 免费女性裸体啪啪无遮挡网站| 人人妻人人添人人爽欧美一区卜| 久热这里只有精品99| 国产成人91sexporn| 麻豆乱淫一区二区| av网站免费在线观看视频| 精品卡一卡二卡四卡免费| 久久精品aⅴ一区二区三区四区| 国产一区二区在线观看av| 麻豆精品久久久久久蜜桃| 在线天堂最新版资源| 精品一区二区三区av网在线观看 | 日韩不卡一区二区三区视频在线| 丰满饥渴人妻一区二区三| 久久人妻熟女aⅴ| 中文字幕av电影在线播放| 成人午夜精彩视频在线观看| 亚洲av福利一区| 亚洲av成人不卡在线观看播放网 | 男女免费视频国产| 好男人视频免费观看在线| 黄片小视频在线播放| av网站在线播放免费| 亚洲av在线观看美女高潮| 1024视频免费在线观看| 成年人午夜在线观看视频| 成人免费观看视频高清| 91老司机精品| a级片在线免费高清观看视频| 日韩中文字幕视频在线看片| 亚洲av在线观看美女高潮| 街头女战士在线观看网站| 欧美日韩av久久| 丝袜美腿诱惑在线| 中文字幕高清在线视频| 在线天堂中文资源库| 国产精品欧美亚洲77777| 国产亚洲精品第一综合不卡| 少妇人妻 视频| 天堂俺去俺来也www色官网| 纵有疾风起免费观看全集完整版| av不卡在线播放| 亚洲中文av在线| 人人澡人人妻人| 一区福利在线观看| av网站在线播放免费| 男女午夜视频在线观看| 精品国产一区二区久久| 女的被弄到高潮叫床怎么办| 成人手机av| 国产黄色视频一区二区在线观看| 国产福利在线免费观看视频| 热99国产精品久久久久久7| 色94色欧美一区二区| 人人妻人人添人人爽欧美一区卜| 午夜福利,免费看| 欧美国产精品一级二级三级| 欧美97在线视频| 极品人妻少妇av视频| 欧美日韩亚洲高清精品| 制服诱惑二区| 一区二区三区激情视频| 视频在线观看一区二区三区| 亚洲欧美一区二区三区国产| 中文天堂在线官网| 在线 av 中文字幕| 久久99一区二区三区| 人妻人人澡人人爽人人| 人体艺术视频欧美日本| netflix在线观看网站| 国产一区二区激情短视频 | 午夜久久久在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲国产av影院在线观看| 久久久久久人人人人人| av免费观看日本| 国产精品国产三级国产专区5o| 一级片免费观看大全| 超碰成人久久| 在线看a的网站| av电影中文网址| 午夜福利在线免费观看网站| 午夜日韩欧美国产| 成人亚洲精品一区在线观看| 毛片一级片免费看久久久久| 女人久久www免费人成看片| 久久99热这里只频精品6学生| 熟女av电影| 国产亚洲欧美精品永久| 欧美成人精品欧美一级黄| 男女免费视频国产| av一本久久久久| 久久久久久免费高清国产稀缺| 亚洲在久久综合| 亚洲伊人久久精品综合| 天天躁日日躁夜夜躁夜夜| 欧美变态另类bdsm刘玥| bbb黄色大片| 五月天丁香电影| 国产在线视频一区二区| 色婷婷av一区二区三区视频| a 毛片基地| 欧美日本中文国产一区发布| 成人国产av品久久久| 久久久久精品性色| 一本一本久久a久久精品综合妖精| 观看av在线不卡| 成年人午夜在线观看视频| 国产一区二区三区综合在线观看| 大码成人一级视频| 老司机靠b影院| 男女下面插进去视频免费观看| 老司机深夜福利视频在线观看 | 99国产精品免费福利视频| 亚洲av电影在线观看一区二区三区| 午夜免费男女啪啪视频观看| 欧美另类一区| 黄片小视频在线播放| 国产人伦9x9x在线观看| 午夜福利,免费看| 久久人人爽av亚洲精品天堂| 在线天堂中文资源库| 国产av国产精品国产| 不卡视频在线观看欧美| 丁香六月欧美| 波野结衣二区三区在线| 亚洲国产最新在线播放| 51午夜福利影视在线观看| 亚洲综合精品二区| 亚洲国产精品一区二区三区在线| 午夜激情av网站| 人成视频在线观看免费观看| 777久久人妻少妇嫩草av网站| 天堂中文最新版在线下载| 美国免费a级毛片| 欧美日韩一区二区视频在线观看视频在线| 高清欧美精品videossex| 你懂的网址亚洲精品在线观看| 亚洲国产欧美网| 女人久久www免费人成看片| 欧美黑人精品巨大| 久久99精品国语久久久| 日韩中文字幕欧美一区二区 | 日韩制服骚丝袜av| 国产乱人偷精品视频| e午夜精品久久久久久久| 高清在线视频一区二区三区| 久久久久视频综合| 在线观看免费视频网站a站| 午夜激情久久久久久久| 涩涩av久久男人的天堂| 国产 一区精品| 国产精品免费视频内射| 丰满饥渴人妻一区二区三| 一区二区三区激情视频| 午夜福利在线免费观看网站| www.熟女人妻精品国产| 国产免费又黄又爽又色| 欧美人与性动交α欧美精品济南到| 亚洲五月色婷婷综合| 国产福利在线免费观看视频| 人妻人人澡人人爽人人| 国产成人精品在线电影| 日本av手机在线免费观看| 欧美在线黄色| 亚洲熟女精品中文字幕| 狂野欧美激情性bbbbbb| 色网站视频免费| 一区二区三区四区激情视频| 美女中出高潮动态图| 免费在线观看黄色视频的| 国产日韩欧美视频二区| 欧美xxⅹ黑人| 毛片一级片免费看久久久久| 十八禁网站网址无遮挡| 亚洲精品国产一区二区精华液| 中文字幕亚洲精品专区| 伊人亚洲综合成人网| 毛片一级片免费看久久久久| kizo精华| 午夜久久久在线观看| 亚洲专区中文字幕在线 | 国产99久久九九免费精品| 国产精品二区激情视频| 久久久久国产一级毛片高清牌| 久久精品国产a三级三级三级| 国产精品 国内视频| 曰老女人黄片| 亚洲成人av在线免费| 日韩成人av中文字幕在线观看| 亚洲国产看品久久| 久久这里只有精品19| 另类精品久久| 亚洲精品视频女| 国产精品熟女久久久久浪| 美女高潮到喷水免费观看| 中国国产av一级| 在线观看免费视频网站a站| 最近最新中文字幕免费大全7| 狠狠精品人妻久久久久久综合| 亚洲国产欧美一区二区综合| 观看美女的网站| 人妻一区二区av| 卡戴珊不雅视频在线播放| 精品国产乱码久久久久久男人| av女优亚洲男人天堂| 超色免费av| 中文字幕制服av| 成人国产av品久久久| 中国三级夫妇交换| 80岁老熟妇乱子伦牲交| av国产精品久久久久影院| 黑丝袜美女国产一区| 99久国产av精品国产电影| 亚洲一区二区三区欧美精品| 欧美 日韩 精品 国产| 欧美日韩一级在线毛片| 国产成人系列免费观看| 中国国产av一级| av有码第一页| 日本欧美国产在线视频| 亚洲精品,欧美精品| 三上悠亚av全集在线观看| 18禁裸乳无遮挡动漫免费视频| avwww免费| 黑人欧美特级aaaaaa片| 亚洲成人av在线免费| 久久久久精品性色| 日本色播在线视频| videos熟女内射| 欧美中文综合在线视频| 99re6热这里在线精品视频| 久久精品国产综合久久久| 人成视频在线观看免费观看| 美女视频免费永久观看网站| 国产97色在线日韩免费| 免费看不卡的av| 亚洲色图综合在线观看| 亚洲欧美一区二区三区黑人| 日韩一区二区三区影片| 亚洲美女视频黄频| 久久精品熟女亚洲av麻豆精品| 久久久久久久久免费视频了| 丝袜美腿诱惑在线| 麻豆乱淫一区二区| 少妇人妻久久综合中文| 日韩一区二区视频免费看| 最近中文字幕高清免费大全6| 最近中文字幕2019免费版| 九色亚洲精品在线播放| 99热国产这里只有精品6| 久久国产精品男人的天堂亚洲| 国产一区二区 视频在线| 欧美最新免费一区二区三区| 国产 精品1| 别揉我奶头~嗯~啊~动态视频 | 久久人人97超碰香蕉20202| 国产男人的电影天堂91| 亚洲人成网站在线观看播放| 欧美精品av麻豆av| 蜜桃在线观看..| 黑人欧美特级aaaaaa片| 不卡av一区二区三区| 新久久久久国产一级毛片| 日本色播在线视频| 悠悠久久av| 成人国产麻豆网| 午夜福利视频在线观看免费| 大码成人一级视频| 国产日韩欧美在线精品| 最新在线观看一区二区三区 | 午夜日韩欧美国产| 99精国产麻豆久久婷婷| 9191精品国产免费久久| 男女之事视频高清在线观看 | 天天躁日日躁夜夜躁夜夜| 国产精品一区二区在线观看99| 亚洲成人手机| 国产精品一区二区在线观看99| av在线播放精品| 亚洲精品一二三| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美成人精品一区二区| 国产一区二区 视频在线| 中文字幕人妻丝袜一区二区 | 中文字幕人妻丝袜一区二区 | 看免费成人av毛片| 丁香六月欧美| 久久女婷五月综合色啪小说| 王馨瑶露胸无遮挡在线观看| 精品国产露脸久久av麻豆| av不卡在线播放| 黄频高清免费视频| 亚洲一区中文字幕在线| 亚洲av欧美aⅴ国产| 日韩 亚洲 欧美在线| 又大又爽又粗| 韩国精品一区二区三区| 久久久久精品国产欧美久久久 | 精品亚洲成国产av| 亚洲情色 制服丝袜| av.在线天堂| 婷婷成人精品国产| 欧美最新免费一区二区三区| 国产精品久久久久久精品电影小说| 少妇被粗大猛烈的视频| 精品人妻在线不人妻| 亚洲成人手机| 亚洲人成电影观看| 国产精品麻豆人妻色哟哟久久| 人人妻人人添人人爽欧美一区卜| 久久99一区二区三区| 在线亚洲精品国产二区图片欧美| 日本猛色少妇xxxxx猛交久久| 九草在线视频观看| 国产黄色视频一区二区在线观看| 精品一区二区免费观看| 黄片小视频在线播放| 欧美人与性动交α欧美软件| 国产熟女午夜一区二区三区| 国产日韩欧美亚洲二区| 两个人看的免费小视频| 一级毛片黄色毛片免费观看视频| 婷婷成人精品国产| 国产激情久久老熟女| 亚洲成色77777| 国产亚洲av高清不卡| 国产精品亚洲av一区麻豆 | 午夜福利影视在线免费观看| 国产一区二区三区av在线| 精品国产一区二区三区四区第35| 国产免费视频播放在线视频| 亚洲av日韩在线播放| 999久久久国产精品视频| 一级毛片电影观看| 国产熟女欧美一区二区| 美女视频免费永久观看网站| 亚洲欧洲日产国产| 欧美最新免费一区二区三区| 在线观看www视频免费| 大码成人一级视频| 看十八女毛片水多多多| 麻豆精品久久久久久蜜桃| 久久久久国产一级毛片高清牌| 国产成人精品福利久久| 国产精品免费视频内射| 十分钟在线观看高清视频www| 99re6热这里在线精品视频| 国产精品久久久久久精品电影小说| 精品一区二区三区四区五区乱码 | 亚洲av电影在线进入| 韩国高清视频一区二区三区| av又黄又爽大尺度在线免费看| 美女高潮到喷水免费观看| 韩国av在线不卡| 亚洲国产欧美在线一区| 国产精品.久久久| 建设人人有责人人尽责人人享有的| 国产男人的电影天堂91| 久久精品人人爽人人爽视色| av女优亚洲男人天堂| 国产 一区精品| 国产野战对白在线观看| 美女扒开内裤让男人捅视频| 国产人伦9x9x在线观看| 日韩免费高清中文字幕av| 欧美日韩成人在线一区二区| 国产野战对白在线观看| 免费av中文字幕在线| 99re6热这里在线精品视频| 99久久综合免费| 天天躁夜夜躁狠狠久久av|