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

    滁州市冬小麥生育期作物需水量和降水變化特征

    2017-09-12 13:22:52金華星胡姍姍龔年祖張鑫童
    水土保持研究 2017年2期
    關(guān)鍵詞:滁州市需水量冬小麥

    金華星, 胡姍姍, 龔年祖, 張鑫童

    (滁州市氣象局, 安徽 滁州239000)

    滁州市冬小麥生育期作物需水量和降水變化特征

    金華星, 胡姍姍, 龔年祖, 張鑫童

    (滁州市氣象局, 安徽 滁州239000)

    基于1961—2013年滁州市氣象數(shù)據(jù),運用集合經(jīng)驗?zāi)B(tài)分解和交叉小波分析,探究了冬小麥生育期需水量和灌溉需水量的變化特征。結(jié)果表明:1961—2013年滁州市冬小麥生育期需水量和灌溉需水量均呈增加趨勢(增勢分別為:2.80,1.48 mm/a);年代變化上,冬小麥生育期需水量逐漸增大,灌溉需水量大致呈現(xiàn)先減少后增加的變化特征。冬小麥生育期作物需水量和灌溉需水量均呈現(xiàn)明顯的多尺度變化特征,表現(xiàn)為多尺度的周期變化特征和趨勢變化特征。作物需水量存在2.65,6.63,13.25,26.50 a的周期,灌溉需水量存在3.11,5.89,10.60,26.5 a的周期。研究時段,作物需水量和灌溉需水量與潛在蒸散量呈同相位分布,與降水量呈顯著的反相位。

    滁州市; 作物需水量; 灌溉蓄水量; 冬小麥生育期; 多尺度變化

    作物需水量是作物全生育期消耗于自身蒸騰和棵間土壤蒸發(fā)的水量總和,它是灌溉決策和水資源規(guī)劃的重要基礎(chǔ)[1]。因此,國內(nèi)外學(xué)者就作物需水量變化特征及其成因開展了一些研究,主要基于野外定位觀測和氣象數(shù)據(jù),結(jié)合作物系數(shù)探究不同作物的需水量、灌溉需求量和灌溉需求指數(shù)的變化等。我國學(xué)者就此開展了深入研究。陳秋帆等[2]基于野外滴灌和觀測試驗探究了云南春作馬鈴薯作物系數(shù)及需水規(guī)律研究。宋妮[3]、楊曉琳[4]和吳燕鋒[5]等分別基于聯(lián)合國糧農(nóng)組織(FAO)推薦的參考作物蒸散計算方法和作物系并結(jié)合氣象數(shù)據(jù)探究了河南省、黃淮海農(nóng)作區(qū)和石河子地區(qū)冬小麥作物需水量時空變化及其氣候要素的影響。符娜等[6]、李搖勇[7]等分析了氣候變化下我國西南地區(qū)和長江中下游地區(qū)水稻不同生育期灌溉需水量的變異性及灌溉需求指數(shù)分布。武勇利等[8]基于改進的SEBS模型對山西省進行作物需水量研究。牛紀(jì)蘋等[9]利用大氣環(huán)流模型Had CM3的輸出和SDSM統(tǒng)計降尺度模型,探究了兩種情境下未來石羊河流域凈灌溉需水量和耗水量。這些學(xué)者們大多采用傳統(tǒng)的趨勢分析法、相關(guān)分析法進行需水量的變化特征及其成因分析,為作物需水量研究提供了大量的理論基礎(chǔ)。但是,氣候要素具有非線性特征,氣候因素驅(qū)動下的作物需水量也可能具有非線性和多尺度特征,而目前從多尺度角度探究需水量的演變特征的研究較少;其次,基于交叉小波分析探究作物需水量與氣象因子相位關(guān)系的研究也鮮有報道。

    滁州市位于安徽省東部,是全國重要的商品糧基地之一。流域內(nèi)河流多且短,河道坡度較大,徑流歷時短,對徑流量滯留極少。因此,農(nóng)田主要靠塘、庫、壩蓄水灌溉[10]。其次,降水的年際變化較大,旱災(zāi)頻繁發(fā)生,防旱減災(zāi)的任務(wù)仍很艱巨[11]。前人就滁州市降水量、氣溫、日照時數(shù)等變化特征開展了研究,而對滁州市冬小麥不同生育期灌溉需水研究較少,且氣候變化下農(nóng)作物需水量和灌溉需水量對區(qū)域農(nóng)業(yè)發(fā)展有重要意義。鑒于此,本研究基于滁州市1961—2013年的氣象數(shù)據(jù)和前人有關(guān)冬小麥作物系數(shù)的研究,探究冬小麥生育期需水量和灌溉需水量的變化特征及其與降水量和潛在蒸散量的關(guān)系,從而為氣候變化下冬小麥需水量和灌溉需水量以及灌溉排水決策提供理論基礎(chǔ)。

    1 數(shù)據(jù)來源與研究方法

    1.1 數(shù)據(jù)來源

    本研究利用滁州市氣象站(32.30N,118.30E),海拔:275 m,逐日氣溫、降水、日照時數(shù)、風(fēng)速、相對濕度等資料實測資料??紤]到資料的完整性,選用1961—2013年的氣象數(shù)據(jù)。

    1.2 研究方法

    1.2.1 冬小麥需水量與灌溉需水量 冬小麥需水量(ETc)計算公式如下[12]:

    ETc=ET0×Kc

    (1)

    式中:ET0為參考作物需水量(mm),采用FAO56推薦的Penman-Monteith 公式[12]計算;Kc為作物系數(shù)。滁州市冬小麥生育期為10月表示次年6月,因此本研究主要探究冬小麥生育期作物需水量的變化特征。參照FAO-56給出的標(biāo)準(zhǔn)條件下作物系數(shù)[13]及已有的相關(guān)研究[14],最終確定冬小麥的作物系數(shù)為1.35。

    由于土壤水和地下水對大氣降水具有巨大調(diào)蓄作用,作物需水量和有效降水量的差值定義為灌溉需水量(IR)。冬小麥灌溉需水量計算公式如下:

    IR=ETc—Pe

    (2)

    式中:IR表示作物灌溉需水量(mm);Pe表示作物生育期有效降雨量(mm)。

    1.2.2 集合經(jīng)驗?zāi)B(tài)分解 本研究采用對冬小麥生育期作物需水量和灌溉需水量序列進行多尺度分析。集合經(jīng)驗?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition,EEMD)是針對EMD方法的不足,提出的一種噪聲輔助數(shù)據(jù)分析方法。該方法目前在水文學(xué)和氣象學(xué)領(lǐng)域得到廣泛應(yīng)用[15]。吳燕鋒等[16]提出利用白噪聲在時間尺度均勻分布的統(tǒng)計特性,在原始信號中加入一定比例的白噪聲,對原始信號和白噪聲組成的總體進行EMD分解,經(jīng)過多次計算,最后取平均使加入的白噪聲相互抵消,加噪聲的次數(shù)越多,平均后的結(jié)果所含的噪聲越小,結(jié)果越接近真實值。因此,EEMD不僅保留了原信號的信息,還可以克服模態(tài)混淆問題。其具體計算步驟為:首先將原始信號加上強度很低的白噪聲,然后將附加噪聲的信號進行EMD分解,得到一組IMF分量。重復(fù)執(zhí)行以上兩部過程,每次添加不同的白噪聲序列,共重復(fù)N次;為了能夠消除人工添加噪聲的影響,將N次EMD分解得到的IMF分量相加,取其均值作為最終的IMF分量:

    (1)

    式中:cj(t)表示對原始信號進行EEMD分解后所得到的第j個IMF分量。本文在EEMD過程中,每次添加白噪聲的振幅為合成序列標(biāo)準(zhǔn)差的20%,集合平均次數(shù)為200次。采用鏡像拓延的方法對數(shù)據(jù)進行處理,以減小邊界效應(yīng)的影響[17]。

    1.2.3 交叉小波分析 利用小波變換探討降水量與氣溫的多時間尺度的變化特征和相關(guān)性。小波變換是把稱為基本小波的函數(shù)做位移后,再在不同的尺度域下對待分析信號做內(nèi)積[18]。交叉小波譜是將小波變換與交叉譜分析相結(jié)合產(chǎn)生的一種信號分析技術(shù),它能夠揭示出兩個時間序列在多時間尺度上的相互關(guān)系。與傳統(tǒng)相關(guān)函數(shù)的只能從時域考察兩個時間序列的相關(guān)關(guān)系,交叉譜能夠在時頻兩域都具有表征氣候信號局部化特征的作用[19]。

    2 結(jié)果與分析

    2.1 冬小麥生育期作物需水量和灌溉需水量年代際變化特征

    圖1為滁州市1961—2013年冬小麥生育期需水量和灌溉需水量年代際變化曲線,可以看出,冬小麥生育期需水量和灌溉需水量均呈增加趨勢。其中,需水量增勢較輕弱(2.80 mm/a),灌溉需水量略有增勢(1.48 mm/a)。過去53 a,滁州市冬小麥生育期需水量均值為869.12 mm。年際變化上,冬小麥生育期需水量在1964年達到最小值,為748.44 mm;2005年達到最大值,為1 015.21 mm;年代變化上,冬小麥生育期需水量自1960s至2011—2014年逐漸增大,在2011—2014年得到年代最大值,為947.19 mm。冬小麥生育期灌溉需水量均值為267.54 mm,變幅較大,其值為11.24~548.40 mm,最大值和最小值分別發(fā)生在2006年和1994年;年代變化上,灌溉需水量波動較大,大致呈現(xiàn)先減少后增加的變化特征,在1970s達到年代最小,為195.92 mm,在2011—2014年達到最大,為382.83 mm。

    圖11961-2013年冬小麥生育期需水量和灌溉需水量年代際變化特征

    2.2 冬小麥生育期作物需水量和灌溉需水量多尺度變化特征

    圖2為過去53 a年作物需水量距平序列EEMD分析結(jié)果,分解為4個IMF和1個殘差(RES)序列。IMF1和IMF2分量表征距平序列的年際變化特征,IMF3和IMF4序列表征距平序列的年代變化特征,RES表征距平序列的變化趨勢特征。表1為冬小麥生育期作物需水量和灌溉需水量的IMF分量的周期特征??梢钥闯觯^去53 a冬小麥生育期作物需水量和灌溉需水量均呈現(xiàn)明顯的多尺度變化特征,表現(xiàn)為多尺度的周期變化特征和趨勢變化特征(RES)。作物需水量年際波動較大,存在2.65 a和6.63 a的周期特征;年代變化尺度上存在13.25 a和26.50 a的周其特征;趨勢項表明冬小麥作物需水量呈明顯的增加趨勢。冬小麥灌溉需水量在年際和年代尺度變化上波動較大,年際尺度上存在3.11 a和5.89 a的周期,年代尺度上存在10.60 a和26.5 a的周期;趨勢項可以看出,冬小麥灌溉需水量呈現(xiàn)先增加、后趨于穩(wěn)定、最后大幅增加的變化特征。

    表1 冬小麥生育期作物需水量和灌溉需水量的IMF分量的周期特征

    2.3 冬小麥生育期作物需水量和灌溉需水量與降水和潛在蒸散量的關(guān)系

    作物需水量和灌溉需水量與氣象因子演變密切相關(guān)[5,7]。由于降水量直接關(guān)乎區(qū)域干濕的演變和水分來源,多種氣象因子通過潛在蒸散量的變化影響冬小麥需水量和土壤供水量,因此,本研究主要探究作物生育期作物需水量和灌溉需水量與降水和潛在蒸散量(ET0)的關(guān)系。圖3是1961—2013年冬小麥生育期降水量和潛在蒸散量變化曲線,可以看出,過去53 a,冬小麥生育期潛在蒸散量和降水量均呈增加趨勢,其中,潛在蒸散量增勢明顯(氣候傾向率為2.07 mm/a),降水量增勢較弱(氣候傾向率為0.22 mm/a)。

    為進一步探究冬小麥生育期作物需水量和灌溉需水量與降水和潛在蒸散量的關(guān)系,基于交叉小波分析分別繪制了冬小麥生育期降水和潛在蒸散量與作物需水量和灌溉需水量的小波交叉能量時頻分布圖(圖4,5)。圖中黑色粗線的區(qū)域表示通過95%紅噪聲檢驗的區(qū)域,箭頭向左和向右分別表示反向位和同相位。作物需水量與潛在蒸散量呈同相位分布,即隨著潛在蒸散量呈增加趨勢,小麥需水量也呈增加趨勢;但兩者的關(guān)系未達到顯著水平(圖4A)。作物需水量與降水量呈顯著的反相位,顯著時頻分布主要表現(xiàn)在:1~2 a周期尺度的1980s中后期,2~4 a周期尺度的2000s,5.5~12.5 a的1960s—1970s前半期,6~12 a的1997—2014年(圖4B)。冬小麥灌溉需水量與潛在蒸散量呈顯著的同相位,即是隨著潛在蒸散量的增加,冬小麥灌溉需水量也呈增加趨勢(圖5A);兩者顯著的時頻分布主要表現(xiàn)在:2~3.5 a周期的1980—2014年、4~10 a周期的2000—2014年、5~13 a周期的1961—1980年和16 a周期的1987—2010年。冬小麥灌溉需水量與降水量呈顯著的反相位,即是降水增加引起灌溉需水量的減少;兩者顯著的時頻分布主要表現(xiàn)在3.8~8 a周期的1970—2010年。

    圖2 1961-2013年冬小麥生育期作物需水量和灌溉需水量的多尺度演變特征

    圖3 1961-2013年冬小麥生育期降水和潛在蒸散量變化特征

    圖4 冬小麥生育期作物需水量與ET0和降水量交叉小波分析

    圖5 冬小麥生育期灌溉需水量與ET0和降水量交叉小波分析

    3 結(jié) 論

    (1) 1961—2013年滁州市冬小麥生育期需水量和灌溉需水量均呈增加趨勢,增勢分別為:2.80 mm/a和1.48 mm/a;年代變化上,冬小麥生育期需水量逐漸增大,灌溉需水量大致呈現(xiàn)先減少后增加的變化特征。

    (2) 冬小麥生育期作物需水量和灌溉需水量均呈現(xiàn)明顯的多尺度變化特征,表現(xiàn)為多尺度的周期變化特征和趨勢變化特征。作物需水量存在2.65,6.63,13.25,26.50 a的周期,灌溉需水量存在3.11,5.89,10.60,26.5 a的周期。

    (3) 研究時段,冬小麥生育期潛在蒸散量和降水量均呈增加趨勢;作物需水量和灌溉需水量與潛在蒸散量呈同相位分布,與降水量呈顯著的反相位。

    作物需水量的和灌溉需水量的變化關(guān)乎作物水分虧缺,對建立合理的灌溉制度有重要的意義,這不僅需要對作物全生育期水分虧缺進行探究,還需要進一步分析不同生育期作物水分需求和灌溉需求。本研究僅僅探究了冬小麥整個生育期的變化特征,未從各生育期分別做深入研究。其次,作物需水量的和灌溉需水量與氣象因子關(guān)系復(fù)雜,兩者關(guān)系有待進一步研究。

    [1] 閆映宇.塔里木灌區(qū)膜下滴灌的棉花需水量及節(jié)水效益[J].水土保持研究,2016,23(1):123-127.

    [2] 陳秋帆,陳勁松,代興梅,等.云南春作馬鈴薯作物系數(shù)及需水規(guī)律研究[J].湖北農(nóng)業(yè)科學(xué),2016,55(3):564-566.

    [3] 宋妮,孫景生,王景雷,等.河南省冬小麥需水量的時空變化及影響因素[J].應(yīng)用生態(tài)學(xué)報,2014,25(6):1693-1700.

    [4] 楊曉琳,宋振偉,王宏,等.黃淮海農(nóng)作區(qū)冬小麥需水量時空變化特征及氣候影響因素分析[J].中國生態(tài)農(nóng)業(yè)學(xué)報,2012,20(3):356-362.

    [5] 吳燕鋒,巴特爾·巴克,羅那那,等.石河子地區(qū)冬小麥生育期需水量變化特征及其氣候成因[J].水土保持通報,2016,36(1):69-74.

    [6] 符娜,劉小剛,張巖石,等.西南地區(qū)水稻灌溉需水量變化規(guī)律[J].生態(tài)學(xué)雜志,2014,33(7):1895-1901.

    [7] 李搖勇,楊曉光,葉搖清,等.1961—2007年長江中下游地區(qū)水稻需水量的變化特征[J].農(nóng)業(yè)工程學(xué)報,2011,27(9):175-183.

    [8] 武永利,田國珍,王云峰.基于改進的SEBS模型的作物需水量研究[J].干旱地區(qū)農(nóng)業(yè)研究,2015,33(1):113-118.

    [9] 牛紀(jì)蘋,粟曉玲,唐澤軍.氣候變化條件下石羊河流域農(nóng)業(yè)灌溉需水量的模擬與預(yù)測[J].干旱地區(qū)農(nóng)業(yè)研究,2016,34(1):206-212.

    [10] 徐天勇.滁州市清流河水環(huán)境容量研究[D].南京:南京大學(xué),2006.

    [11] 徐志龍,曹陽,楊敏.1951—2005年海河流域汛期降水量的時空變化特征分析[J].水文,2009,29(1):85-88.

    [12] Karimov A K,nek J, Hanjra M A, et al. Effects of the shallow water table on water use of winter wheat and ecosystem health:Implications for unlocking the potential of groundwater in the Fergana Valley(Central Asia)[J]. Agricultural Water Management, 2014,131(1):57-69.

    [13] Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration-Guidelines for computing crop water requirements-FAO Irrigation and drainage paper56[M]. Rome:FAO, 1998.

    [14] 何斌.中國農(nóng)業(yè)干旱風(fēng)險定量評價與空間格局分析[D].北京:中國科學(xué)院大學(xué)(中國科學(xué)院地理科學(xué)與資源研究所),2010.

    [15] Yang Z, Yu Z, Xie C, et al. Application of Hilbert-Huang Transform to acoustic emission signal for burn feature extraction in surface grinding process[J]. Measurement, 2014,47:14-21.

    [16] 吳燕鋒,巴特爾,巴克,等.基于EEMD的杜尚別市1950—2013年降水多尺度分析[J].干旱區(qū)資源與環(huán)境,2015,29(6):152-157.

    [17] 張余慶,陳昌春,姚鑫,等.江西省信江流域極端降水時空變化特征[J].水土保持研究,2015,22(4):189-194.

    [18] Valdés-Galicia J F, Velasco V M. Variations of mid-term periodicities in solar activity physical phenomena [J]. Advances in Space Research, 2008,41(2):297-305.

    [19] Wu Y, Bake B, Zhang J, et al. Spatio-temporal patterns of drought in North Xinjiang, China, 1961—2012 based on meteorological drought index[J]. Journal of Arid Land, 2015,7(4):527-543.

    VariationofWaterRequirementandPrecipitationinWheatGrowthStagesfrom1961to2013

    JIN Huaxing, HU Shanshan, GONG Nianzu, ZHANG Xintong

    (MeteorologicalBureauofChuzhouCity,Chuzhou,Anhui239000,China)

    Based on the data from Meteorological Bureau of Chuzhou City and using ensemble empirical mode decomposition and cross wavelet analysis, the crop and irrigation water requirements during the wheat growth stages were analyzed. The results showed that the crop and irrigation water requirements experienced increasing trend with the rate of 2.80 mm/year and 1.48 mm/year from 1961 to 2013, respectively. For decadal change, the crop water requirement presented gradually increasing trend, while the irrigation water requirement exhibited decrease first and then increase. The crop water requirement has the periodicities of 2.65 years, 6.63 years, 13.25 years and 26.50 years, the irrigation water requirement has the periodicities of 3.11 years, 5.89 years, 10.60 years and 26.5 years. The crop and irrigation water requirements have the in-phase relationship with evapotranspiration and significant anti-phase relationship with precipitation.

    Chouzhou City; crop water requirement; irrigation water requirement; wheat growth stages; multiscale variation

    2016-05-04

    :2016-05-11

    安徽省氣象科技發(fā)展基金“滁州市短期溫度降水多模式集成預(yù)報”(KM201404)

    金華星(1973—),男,安徽鳳陽人,工程師,本科,主要研究方向:綜合氣象觀測業(yè)務(wù)、氣象裝備保障和人工影響天氣等工作。E-mail:jhx210@aliyun.com

    S512.1+1

    :A

    :1005-3409(2017)02-0178-06

    猜你喜歡
    滁州市需水量冬小麥
    安徽省滁州市瑯琊山銅礦床成礦模式研究
    南京城市中心區(qū)生態(tài)需水量探析
    習(xí)作點評
    中國篆刻(2019年9期)2019-09-19 00:48:18
    習(xí)作點評
    中國篆刻(2019年1期)2019-01-26 01:37:40
    淺析如何提高滁州市城市防洪標(biāo)準(zhǔn)
    組合預(yù)測模型在區(qū)域需水量預(yù)測中的應(yīng)用
    甘肅冬小麥田
    植物保護(2017年1期)2017-02-13 06:44:34
    帕滿灌區(qū)供需水量平衡分析
    大洋河生態(tài)需水量分析與評價
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    一区二区三区激情视频| 午夜福利乱码中文字幕| 蜜桃在线观看..| 黄片播放在线免费| 免费日韩欧美在线观看| 亚洲欧美精品综合一区二区三区| 妹子高潮喷水视频| 91字幕亚洲| 搡老岳熟女国产| 欧美午夜高清在线| 亚洲中文av在线| 国产有黄有色有爽视频| 亚洲精品日韩在线中文字幕| 十分钟在线观看高清视频www| 日本猛色少妇xxxxx猛交久久| 9色porny在线观看| 在线观看舔阴道视频| 丝袜美足系列| 50天的宝宝边吃奶边哭怎么回事| 亚洲自偷自拍图片 自拍| 黄色视频在线播放观看不卡| 欧美日韩亚洲高清精品| 99国产极品粉嫩在线观看| 午夜免费观看性视频| 乱人伦中国视频| 国产精品麻豆人妻色哟哟久久| 丝袜喷水一区| 亚洲,欧美精品.| 90打野战视频偷拍视频| 国产伦人伦偷精品视频| 国产黄色免费在线视频| 午夜老司机福利片| 欧美中文综合在线视频| 日日爽夜夜爽网站| 亚洲一区中文字幕在线| 国产一级毛片在线| 丝瓜视频免费看黄片| 黄片小视频在线播放| 国产精品自产拍在线观看55亚洲 | 热re99久久国产66热| 午夜福利影视在线免费观看| 一区二区三区乱码不卡18| 欧美乱码精品一区二区三区| 中文字幕制服av| 国产在线观看jvid| 亚洲成av片中文字幕在线观看| 一个人免费在线观看的高清视频 | 香蕉国产在线看| 夫妻午夜视频| 免费日韩欧美在线观看| 久久久久久久国产电影| 久久久精品免费免费高清| 国产欧美亚洲国产| 51午夜福利影视在线观看| 国产精品偷伦视频观看了| a级毛片在线看网站| 精品欧美一区二区三区在线| 黑人操中国人逼视频| 一级黄色大片毛片| 精品人妻1区二区| 在线永久观看黄色视频| 久久精品人人爽人人爽视色| 国产精品偷伦视频观看了| 爱豆传媒免费全集在线观看| 久久青草综合色| 五月天丁香电影| 美女视频免费永久观看网站| 宅男免费午夜| 啪啪无遮挡十八禁网站| 欧美日韩福利视频一区二区| 国产人伦9x9x在线观看| 日韩免费高清中文字幕av| 少妇 在线观看| 久久精品国产a三级三级三级| 天天影视国产精品| 五月开心婷婷网| 电影成人av| bbb黄色大片| 亚洲国产中文字幕在线视频| 不卡av一区二区三区| 久久久久久久大尺度免费视频| 亚洲精品美女久久av网站| 久久久国产成人免费| 亚洲色图综合在线观看| 男男h啪啪无遮挡| 我要看黄色一级片免费的| videos熟女内射| 叶爱在线成人免费视频播放| 两人在一起打扑克的视频| 大片电影免费在线观看免费| 黑人操中国人逼视频| 动漫黄色视频在线观看| 美女中出高潮动态图| 久久女婷五月综合色啪小说| 一区在线观看完整版| 2018国产大陆天天弄谢| 亚洲精品国产av成人精品| 热re99久久精品国产66热6| 亚洲第一青青草原| 国产av精品麻豆| 成人国产一区最新在线观看| 亚洲久久久国产精品| 日本91视频免费播放| 999久久久精品免费观看国产| 亚洲中文日韩欧美视频| 久久久久久久久免费视频了| 久久ye,这里只有精品| 各种免费的搞黄视频| 视频区图区小说| 亚洲精品久久久久久婷婷小说| 十八禁人妻一区二区| 国产精品av久久久久免费| 岛国毛片在线播放| 久久狼人影院| 性高湖久久久久久久久免费观看| 中文字幕人妻熟女乱码| 国产日韩一区二区三区精品不卡| 12—13女人毛片做爰片一| 他把我摸到了高潮在线观看 | 久久狼人影院| 天堂中文最新版在线下载| 三级毛片av免费| 中国国产av一级| 国产精品熟女久久久久浪| 久久人人97超碰香蕉20202| a级片在线免费高清观看视频| 51午夜福利影视在线观看| 国产精品秋霞免费鲁丝片| 国产精品免费大片| 动漫黄色视频在线观看| kizo精华| 日韩人妻精品一区2区三区| 女人精品久久久久毛片| 国产成人免费无遮挡视频| 麻豆乱淫一区二区| 国产又色又爽无遮挡免| 国产成人精品久久二区二区91| 国产成人av教育| 久久综合国产亚洲精品| 精品少妇一区二区三区视频日本电影| 天天躁狠狠躁夜夜躁狠狠躁| 80岁老熟妇乱子伦牲交| 色94色欧美一区二区| 一区在线观看完整版| 欧美大码av| 亚洲国产精品成人久久小说| 男人添女人高潮全过程视频| 51午夜福利影视在线观看| 欧美黄色淫秽网站| 国产成人精品久久二区二区免费| 老熟妇仑乱视频hdxx| 亚洲精品一区蜜桃| 黄片小视频在线播放| 99re6热这里在线精品视频| bbb黄色大片| 熟女少妇亚洲综合色aaa.| 丰满少妇做爰视频| 免费黄频网站在线观看国产| 51午夜福利影视在线观看| 后天国语完整版免费观看| 亚洲国产中文字幕在线视频| 侵犯人妻中文字幕一二三四区| 久久国产精品男人的天堂亚洲| 亚洲久久久国产精品| 久久免费观看电影| 久久久国产一区二区| 亚洲精品av麻豆狂野| 超碰97精品在线观看| 人妻 亚洲 视频| 交换朋友夫妻互换小说| 亚洲欧洲日产国产| 老熟妇乱子伦视频在线观看 | videos熟女内射| 一个人免费在线观看的高清视频 | 欧美乱码精品一区二区三区| 国产深夜福利视频在线观看| 国产精品自产拍在线观看55亚洲 | 在线十欧美十亚洲十日本专区| 午夜免费成人在线视频| 黑丝袜美女国产一区| 国产欧美日韩一区二区精品| 日韩免费高清中文字幕av| 另类亚洲欧美激情| 免费观看a级毛片全部| 欧美久久黑人一区二区| 99国产精品99久久久久| 亚洲国产av新网站| 伦理电影免费视频| 老司机在亚洲福利影院| 女人高潮潮喷娇喘18禁视频| 一本一本久久a久久精品综合妖精| 电影成人av| 人人妻人人澡人人看| 国产成人免费观看mmmm| 2018国产大陆天天弄谢| 一区二区三区乱码不卡18| 久久久久久久国产电影| 少妇被粗大的猛进出69影院| 婷婷色av中文字幕| 天天影视国产精品| 色94色欧美一区二区| 黄片大片在线免费观看| 夫妻午夜视频| 老司机亚洲免费影院| 纯流量卡能插随身wifi吗| 精品人妻熟女毛片av久久网站| 亚洲精品中文字幕一二三四区 | 一级毛片电影观看| 啦啦啦中文免费视频观看日本| 12—13女人毛片做爰片一| 久久影院123| 久久国产精品男人的天堂亚洲| 久久精品国产a三级三级三级| 中文字幕制服av| 一级a爱视频在线免费观看| 一二三四在线观看免费中文在| 丰满少妇做爰视频| 久久精品亚洲熟妇少妇任你| 一本综合久久免费| 亚洲国产欧美网| 国产99久久九九免费精品| 精品高清国产在线一区| 手机成人av网站| 老汉色∧v一级毛片| 精品亚洲乱码少妇综合久久| 日韩制服丝袜自拍偷拍| 99国产综合亚洲精品| 少妇人妻久久综合中文| 99久久99久久久精品蜜桃| 国产主播在线观看一区二区| 90打野战视频偷拍视频| 成人国产av品久久久| 成在线人永久免费视频| 少妇的丰满在线观看| 老司机深夜福利视频在线观看 | 日本a在线网址| 大码成人一级视频| 亚洲专区字幕在线| 午夜福利视频在线观看免费| 最近最新免费中文字幕在线| 人人澡人人妻人| 国产一级毛片在线| 天天添夜夜摸| 各种免费的搞黄视频| 男人爽女人下面视频在线观看| 亚洲精品在线美女| 捣出白浆h1v1| 麻豆国产av国片精品| 成年人午夜在线观看视频| 在线永久观看黄色视频| 美女主播在线视频| 午夜91福利影院| 国产在线观看jvid| 亚洲国产欧美一区二区综合| 久久久久视频综合| 成年av动漫网址| 国产亚洲精品久久久久5区| 日本撒尿小便嘘嘘汇集6| 亚洲自偷自拍图片 自拍| 一区福利在线观看| 国产一区二区三区av在线| 亚洲成国产人片在线观看| 曰老女人黄片| 成在线人永久免费视频| 丁香六月欧美| 中文字幕人妻丝袜一区二区| av一本久久久久| 涩涩av久久男人的天堂| 黄色片一级片一级黄色片| 少妇 在线观看| 精品国产一区二区三区久久久樱花| 久久热在线av| 日韩欧美一区二区三区在线观看 | 中文字幕人妻丝袜一区二区| av天堂久久9| 啦啦啦免费观看视频1| 午夜老司机福利片| 12—13女人毛片做爰片一| 又大又爽又粗| 久久国产精品影院| 久久女婷五月综合色啪小说| 高潮久久久久久久久久久不卡| 亚洲五月色婷婷综合| 美女扒开内裤让男人捅视频| 永久免费av网站大全| 在线精品无人区一区二区三| 午夜免费成人在线视频| 一级黄色大片毛片| cao死你这个sao货| 久久久水蜜桃国产精品网| 国产欧美日韩一区二区三 | 美女福利国产在线| www.熟女人妻精品国产| 日韩大码丰满熟妇| 1024视频免费在线观看| 麻豆av在线久日| 国产免费现黄频在线看| 国产一区二区三区综合在线观看| 黄色视频不卡| 亚洲国产日韩一区二区| 亚洲精品中文字幕在线视频| 老司机在亚洲福利影院| 大片电影免费在线观看免费| 久久久水蜜桃国产精品网| 一边摸一边抽搐一进一出视频| 中文字幕最新亚洲高清| 国产精品秋霞免费鲁丝片| a级毛片在线看网站| cao死你这个sao货| 亚洲国产精品一区二区三区在线| 91成年电影在线观看| 精品亚洲成国产av| 色婷婷久久久亚洲欧美| 国产精品免费视频内射| 国产无遮挡羞羞视频在线观看| 老熟女久久久| 国产免费现黄频在线看| 精品欧美一区二区三区在线| 日本a在线网址| 真人做人爱边吃奶动态| 欧美黄色片欧美黄色片| 午夜老司机福利片| 久久久久网色| 一本一本久久a久久精品综合妖精| 久久精品久久久久久噜噜老黄| 精品久久久精品久久久| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久久久久婷婷小说| 老司机福利观看| 亚洲第一欧美日韩一区二区三区 | av福利片在线| 狠狠狠狠99中文字幕| 亚洲国产精品999| 免费黄频网站在线观看国产| 妹子高潮喷水视频| 欧美 亚洲 国产 日韩一| www.自偷自拍.com| 免费黄频网站在线观看国产| 婷婷色av中文字幕| 日韩一区二区三区影片| 久久久久网色| 丰满饥渴人妻一区二区三| 91精品国产国语对白视频| 国产亚洲av高清不卡| 捣出白浆h1v1| 成人三级做爰电影| 国产成人免费无遮挡视频| 日本五十路高清| 大型av网站在线播放| 十分钟在线观看高清视频www| 亚洲精品第二区| 大香蕉久久网| 亚洲第一欧美日韩一区二区三区 | 曰老女人黄片| 亚洲性夜色夜夜综合| 精品一品国产午夜福利视频| 高清欧美精品videossex| 欧美激情久久久久久爽电影 | 我的亚洲天堂| 纯流量卡能插随身wifi吗| 丝袜美腿诱惑在线| 日本撒尿小便嘘嘘汇集6| 国产免费视频播放在线视频| 亚洲专区中文字幕在线| 午夜福利视频在线观看免费| 日韩熟女老妇一区二区性免费视频| 中文欧美无线码| 久久人人爽av亚洲精品天堂| 国产日韩欧美视频二区| 青春草视频在线免费观看| 成人免费观看视频高清| 国产区一区二久久| 99九九在线精品视频| 亚洲专区中文字幕在线| 日本精品一区二区三区蜜桃| 免费在线观看黄色视频的| 成人18禁高潮啪啪吃奶动态图| 中文字幕制服av| 日韩熟女老妇一区二区性免费视频| 自拍欧美九色日韩亚洲蝌蚪91| a级毛片黄视频| 国产一区二区激情短视频 | 亚洲九九香蕉| 在线十欧美十亚洲十日本专区| 香蕉丝袜av| av电影中文网址| 精品久久久久久久毛片微露脸 | 成年人免费黄色播放视频| 精品一区二区三区av网在线观看 | 美女脱内裤让男人舔精品视频| 丝袜美足系列| 少妇裸体淫交视频免费看高清 | 国产av精品麻豆| 女人爽到高潮嗷嗷叫在线视频| 最近中文字幕2019免费版| 亚洲av片天天在线观看| 国产精品.久久久| 美女福利国产在线| 精品久久久久久久毛片微露脸 | 啦啦啦 在线观看视频| av在线播放精品| 日韩,欧美,国产一区二区三区| 精品国产乱子伦一区二区三区 | 国产成+人综合+亚洲专区| 天天添夜夜摸| 免费日韩欧美在线观看| 色精品久久人妻99蜜桃| 精品国产超薄肉色丝袜足j| 嫁个100分男人电影在线观看| 亚洲精品成人av观看孕妇| 黄片大片在线免费观看| 少妇粗大呻吟视频| 人人妻人人澡人人看| 国产在线一区二区三区精| 国产精品秋霞免费鲁丝片| 99国产极品粉嫩在线观看| 精品国产一区二区三区久久久樱花| 亚洲国产日韩一区二区| 老司机影院成人| 国产高清视频在线播放一区 | a级毛片在线看网站| 国产激情久久老熟女| 超碰97精品在线观看| 90打野战视频偷拍视频| 91麻豆av在线| 九色亚洲精品在线播放| av线在线观看网站| 男女下面插进去视频免费观看| 黑人猛操日本美女一级片| 久久久久国产精品人妻一区二区| 国产成人欧美| 波多野结衣av一区二区av| 少妇精品久久久久久久| 日本黄色日本黄色录像| 国产成人影院久久av| 成人18禁高潮啪啪吃奶动态图| 老汉色∧v一级毛片| 精品人妻熟女毛片av久久网站| 免费观看a级毛片全部| 淫妇啪啪啪对白视频 | 精品一品国产午夜福利视频| 桃红色精品国产亚洲av| 国产99久久九九免费精品| 极品少妇高潮喷水抽搐| 亚洲一码二码三码区别大吗| 中国国产av一级| 国产深夜福利视频在线观看| 狂野欧美激情性xxxx| 纯流量卡能插随身wifi吗| 国产一区二区在线观看av| 欧美av亚洲av综合av国产av| www.熟女人妻精品国产| 中文字幕人妻丝袜一区二区| 丝袜美足系列| 国产成人精品无人区| 亚洲精品一卡2卡三卡4卡5卡 | a级毛片黄视频| 18在线观看网站| 丝袜在线中文字幕| 日韩制服骚丝袜av| 国产又爽黄色视频| 丰满人妻熟妇乱又伦精品不卡| 宅男免费午夜| 男男h啪啪无遮挡| 天堂中文最新版在线下载| 亚洲精品一卡2卡三卡4卡5卡 | 一本综合久久免费| 亚洲av日韩在线播放| 一边摸一边抽搐一进一出视频| 男女边摸边吃奶| 国产成人欧美| 日韩大片免费观看网站| av视频免费观看在线观看| 精品少妇久久久久久888优播| 性高湖久久久久久久久免费观看| 欧美性长视频在线观看| 国产亚洲精品久久久久5区| 水蜜桃什么品种好| 久9热在线精品视频| 欧美 日韩 精品 国产| 国产成人免费观看mmmm| av视频免费观看在线观看| 国产高清videossex| 国产精品久久久人人做人人爽| 久久热在线av| 99国产精品免费福利视频| 精品少妇内射三级| 国产福利在线免费观看视频| 999久久久精品免费观看国产| 亚洲国产av影院在线观看| 精品国产一区二区三区四区第35| av天堂久久9| 亚洲av日韩精品久久久久久密| 久久精品亚洲熟妇少妇任你| 国产黄频视频在线观看| 99国产极品粉嫩在线观看| 亚洲avbb在线观看| 精品少妇黑人巨大在线播放| 国产成+人综合+亚洲专区| 久久天躁狠狠躁夜夜2o2o| 不卡一级毛片| 一本—道久久a久久精品蜜桃钙片| 考比视频在线观看| 91精品国产国语对白视频| 欧美日韩亚洲综合一区二区三区_| 欧美日韩av久久| 久9热在线精品视频| 无遮挡黄片免费观看| svipshipincom国产片| 欧美激情极品国产一区二区三区| 一本大道久久a久久精品| 丝袜人妻中文字幕| 乱人伦中国视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲天堂av无毛| 丝瓜视频免费看黄片| 亚洲avbb在线观看| www.av在线官网国产| 老司机亚洲免费影院| 国产色视频综合| 精品亚洲成a人片在线观看| √禁漫天堂资源中文www| 国产欧美亚洲国产| 亚洲欧美清纯卡通| 丝袜人妻中文字幕| 亚洲少妇的诱惑av| 搡老乐熟女国产| 成人亚洲精品一区在线观看| av免费在线观看网站| 女警被强在线播放| 久热这里只有精品99| 窝窝影院91人妻| 狠狠婷婷综合久久久久久88av| 侵犯人妻中文字幕一二三四区| 国产成人一区二区三区免费视频网站| 青春草视频在线免费观看| 精品国内亚洲2022精品成人 | 国产真人三级小视频在线观看| 人人澡人人妻人| 肉色欧美久久久久久久蜜桃| 亚洲成人免费电影在线观看| 最新在线观看一区二区三区| 免费在线观看完整版高清| 国产亚洲精品一区二区www | 伊人亚洲综合成人网| 黄色 视频免费看| 国产成人欧美在线观看 | 久久国产精品影院| 亚洲av国产av综合av卡| 午夜免费成人在线视频| 精品一区在线观看国产| 午夜福利在线免费观看网站| 国产亚洲av高清不卡| 99热网站在线观看| kizo精华| 国产在视频线精品| 黄片大片在线免费观看| 亚洲少妇的诱惑av| 亚洲av电影在线观看一区二区三区| 久久精品人人爽人人爽视色| 精品少妇一区二区三区视频日本电影| 免费不卡黄色视频| 亚洲第一青青草原| 热99国产精品久久久久久7| 午夜精品久久久久久毛片777| h视频一区二区三区| 免费一级毛片在线播放高清视频 | 日本猛色少妇xxxxx猛交久久| 国产精品九九99| 丝袜脚勾引网站| 亚洲国产av新网站| 精品乱码久久久久久99久播| 在线十欧美十亚洲十日本专区| 精品熟女少妇八av免费久了| av欧美777| 啦啦啦中文免费视频观看日本| 51午夜福利影视在线观看| 精品国产超薄肉色丝袜足j| 自拍欧美九色日韩亚洲蝌蚪91| 美女高潮到喷水免费观看| 我要看黄色一级片免费的| 亚洲欧美一区二区三区黑人| av免费在线观看网站| 国产熟女午夜一区二区三区| 亚洲av成人不卡在线观看播放网 | 黑人巨大精品欧美一区二区蜜桃| 久久这里只有精品19| 啦啦啦免费观看视频1| 成人影院久久| 亚洲国产精品999| 人人妻人人添人人爽欧美一区卜| 国产一区二区三区综合在线观看| 亚洲成人手机| 久久人妻熟女aⅴ| 欧美日韩av久久| 中文欧美无线码| 亚洲伊人色综图| 亚洲国产中文字幕在线视频| 一本一本久久a久久精品综合妖精| 中文字幕av电影在线播放| 大陆偷拍与自拍| 老司机在亚洲福利影院| 欧美成狂野欧美在线观看| 好男人电影高清在线观看| 精品欧美一区二区三区在线| 亚洲成人国产一区在线观看| 国产成人免费无遮挡视频| 欧美xxⅹ黑人| 一区二区三区激情视频| 日本wwww免费看| 欧美黑人精品巨大| 亚洲色图综合在线观看| 亚洲七黄色美女视频| 丰满饥渴人妻一区二区三| 少妇人妻久久综合中文|