亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于積溫產(chǎn)量模型確定山東夏玉米拔節(jié)前后的極端高溫閾值*

        2017-12-18 06:49:06馮一淳張藝璇蔡成良
        中國(guó)農(nóng)業(yè)氣象 2017年12期
        關(guān)鍵詞:有效積溫日數(shù)夏玉米

        張 琪,唐 婕,馮一淳,張藝璇,蔡成良

        ?

        基于積溫產(chǎn)量模型確定山東夏玉米拔節(jié)前后的極端高溫閾值*

        張 琪1,2,唐 婕1,馮一淳1,張藝璇1,蔡成良1

        (1.南京信息工程大學(xué)應(yīng)用氣象學(xué)院,南京 210024;2.江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室,南京 210024)

        參考已有的積溫-產(chǎn)量模型方法,利用山東省10個(gè)氣象站點(diǎn)1983-2012年逐日最高氣溫、最低氣溫?cái)?shù)據(jù)模擬其間逐小時(shí)溫度,分拔節(jié)前和拔節(jié)后兩個(gè)時(shí)段確定山東省夏玉米的極端高溫閾值,并對(duì)研究區(qū)極端高溫事件的時(shí)空分布特征進(jìn)行分析。結(jié)果表明:夏玉米生育期溫度主要分布在17~35℃區(qū)間,拔節(jié)后不同溫度持續(xù)時(shí)長(zhǎng)較大;有效積溫轉(zhuǎn)化為實(shí)際產(chǎn)量的比例系數(shù)有隨時(shí)間增大的趨勢(shì);山東省夏玉米拔節(jié)前、后兩個(gè)階段極端高溫閾值分別為35.2℃和34.5℃;基于此閾值統(tǒng)計(jì)的極端高溫發(fā)生日數(shù)由西部平原地區(qū)向東部沿海半島地區(qū)減少,拔節(jié)后較拔節(jié)前頻發(fā),研究期內(nèi)極端高溫發(fā)生日數(shù)有隨時(shí)間增多的趨勢(shì)。

        極端高溫閾值;積溫產(chǎn)量模型;夏玉米;拔節(jié)

        在以變暖為主要特征的全球氣候變化背景下,各類極端天氣事件尤以極端高溫,呈現(xiàn)增多增強(qiáng)的趨勢(shì),這種趨勢(shì)在未來(lái)一段時(shí)間很可能繼續(xù)延續(xù)[1-3]。而極端高溫除直接影響農(nóng)作物生長(zhǎng)發(fā)育外,持續(xù)高溫還會(huì)誘發(fā)或加劇干旱進(jìn)而影響農(nóng)業(yè)生產(chǎn)。近年來(lái),持續(xù)高溫造成的農(nóng)業(yè)損失越來(lái)越突出,如2016年華北、西北7-8月連續(xù)≥35℃的高溫對(duì)玉米開(kāi)花灌漿造成不利影響,降低了結(jié)實(shí)率,導(dǎo)致減產(chǎn)[4]。相對(duì)于其它農(nóng)業(yè)氣象災(zāi)害,極端高溫是一個(gè)重要但較少被關(guān)注的領(lǐng)域[5],其相關(guān)研究亟待加強(qiáng)。

        閾值確定是研究極端高溫影響的前提和關(guān)鍵,以往研究多采用百分位法確定極端事件閾值[6-8],該方法從數(shù)據(jù)本身結(jié)構(gòu)出發(fā)并未考慮承載體的特征,然而不同的承載體對(duì)災(zāi)害敏感程度不同,高溫閾值也有差異。目前極端高溫對(duì)農(nóng)業(yè)影響的研究中,多參考大田試驗(yàn)或以往經(jīng)驗(yàn)確定閾值[9-11],通常認(rèn)為在營(yíng)養(yǎng)生長(zhǎng)階段一般36℃以上會(huì)對(duì)玉米產(chǎn)量帶來(lái)影響,生殖生長(zhǎng)階段日最高溫度超過(guò)32~35℃玉米生長(zhǎng)會(huì)受到影響;尹小剛等[12]參考Lobell等[13]研究美國(guó)玉米生產(chǎn)時(shí)確定了30℃作為閾值,研究東北地區(qū)春玉米各生育階段極端高溫發(fā)生的時(shí)間變化特征。但由于不同地區(qū)玉米品種、氣候以及下墊面情況的不同,借鑒其它地區(qū)的閾值進(jìn)行研究將給結(jié)果帶來(lái)一定偏差。

        針對(duì)以往極端高溫災(zāi)害研究報(bào)道較少且閾值確定忽略承載體特征等問(wèn)題,本研究以夏玉米主產(chǎn)區(qū)同時(shí)也是高溫災(zāi)害頻發(fā)區(qū)的山東省為例,擬通過(guò)積溫產(chǎn)量模型的方法確定該地區(qū)夏玉米不同生育階段極端高溫災(zāi)害的閾值,并分析極端高溫災(zāi)害的時(shí)空變化特征,以期為極端事件閾值確定開(kāi)拓新思路,促進(jìn)氣候變化背景下極端事件對(duì)農(nóng)業(yè)影響的研究。

        1 材料與方法

        1.1 研究區(qū)概況

        以華北平原的山東省為研究區(qū),該區(qū)為主要夏玉米產(chǎn)區(qū)(圖1),暖溫帶大陸性氣候,降水主要集中在6-9月,雨熱同期,利于夏玉米的生長(zhǎng)。但降水分布不均勻、高溫頻發(fā)導(dǎo)致夏玉米產(chǎn)量波動(dòng)較大。本研究以拔節(jié)為界將玉米整個(gè)生育期分為兩個(gè)階段:拔節(jié)前,玉米主要進(jìn)行營(yíng)養(yǎng)生長(zhǎng);拔節(jié)后,開(kāi)始有生殖生長(zhǎng)活動(dòng)。山東地區(qū)夏玉米拔節(jié)主要發(fā)生在7月下旬,生育期平均長(zhǎng)度98d左右,拔節(jié)前持續(xù)38d,拔節(jié)后60d左右。全球變暖影響下,研究區(qū)夏玉米播種日期有推遲的趨勢(shì),導(dǎo)致拔節(jié)前的天數(shù)有減少趨勢(shì)。

        圖1 研究區(qū)及站點(diǎn)分布

        1.2 資料來(lái)源

        研究區(qū)10個(gè)農(nóng)業(yè)氣象站點(diǎn)1983-2012年逐年玉米發(fā)育期觀測(cè)資料,主要包括播種、拔節(jié)、成熟日期,來(lái)源于國(guó)家氣象信息中心;1983-2012年逐日最高氣溫、最低氣溫?cái)?shù)據(jù)為10個(gè)農(nóng)業(yè)氣象站對(duì)應(yīng)的氣象臺(tái)站觀測(cè)數(shù)據(jù),以該數(shù)據(jù)為基礎(chǔ)模擬其間逐小時(shí)溫度。

        1.3 積溫產(chǎn)量模型

        研究氣溫對(duì)作物產(chǎn)量的影響時(shí),一般認(rèn)為溫度對(duì)產(chǎn)量的影響是一個(gè)非線性的逐漸累積的過(guò)程,產(chǎn)量與整個(gè)生育階段有效積溫成比例關(guān)系。基于這一思路并參考Schlenker等研究氣溫對(duì)美國(guó)玉米產(chǎn)量影響時(shí)建立的回歸方程,用非線性函數(shù)表示實(shí)際產(chǎn)量與生育期溫度之間的關(guān)系[14],即在站點(diǎn)i第t年的夏玉米實(shí)際產(chǎn)量yit為

        Hawks等[15]將8℃和32℃作為玉米有效積溫的上下限,由式(2)計(jì)算。例如10℃,其貢獻(xiàn)了2℃·h的有效積溫,若生育期內(nèi)10℃持續(xù)時(shí)長(zhǎng)為3h,則該溫度在生育期內(nèi)貢獻(xiàn)了6℃的有效積溫。

        由一日內(nèi)最高氣溫和最低氣溫構(gòu)建正弦曲線,模擬氣溫日變化以獲取逐小時(shí)氣溫[16],即

        2 結(jié)果與分析

        2.2 比例系數(shù)k的統(tǒng)計(jì)分析

        k為積溫產(chǎn)量模型中有效積溫轉(zhuǎn)化為實(shí)際產(chǎn)量的比例系數(shù),通過(guò)式(1)得到歷年各站點(diǎn)拔節(jié)前和拔節(jié)后兩個(gè)時(shí)段有效積溫轉(zhuǎn)化為產(chǎn)量的比例系數(shù)k,并計(jì)算歷年山東省10個(gè)站點(diǎn)的平均值,結(jié)果見(jiàn)圖3。由圖中可以看出,夏玉米全生育期、拔節(jié)前和拔節(jié)后的k值均隨著時(shí)間呈增大趨勢(shì),上升趨勢(shì)達(dá)到P<0.01的極顯著水平。整體上,1993年之前k值變化不明顯,之后有較為明顯的升高趨勢(shì),即同樣的有效積溫轉(zhuǎn)化成的實(shí)際產(chǎn)量增多,說(shuō)明這期間農(nóng)業(yè)投入的增加等使夏玉米對(duì)溫度的利用效率提高。

        圖2 1983-2012年山東省夏玉米拔節(jié)前(a)、后(b)不同溫度持續(xù)時(shí)長(zhǎng)

        注:每個(gè)溫度對(duì)應(yīng)的矩形盒內(nèi)部橫線代表所有站點(diǎn)研究期內(nèi)該溫度持續(xù)時(shí)長(zhǎng)的中位數(shù),矩形盒的頂線和底線分別代表上、下四分位,上部標(biāo)記短線為最大值,下部標(biāo)記短線為最小值

        Note: In the box-plot of each temperature, the horizontal line inside the box shows the median duration exposure to this temperature among all stations and years, the top and bottom of the box show the upper quartile and lower quartile. The upper short lines are the maximum values, and the bottom short lines are the minimum values

        圖3 1983-2012年山東省夏玉米產(chǎn)量與有效積溫比例系數(shù)(k)的年際變化

        2.3 極端高溫閾值的確定

        分夏玉米拔節(jié)前和拔節(jié)后兩個(gè)時(shí)段確定。圖4中灰色曲線為所有站點(diǎn)各年份利用式(1)的積分函數(shù)得到的曲線,黑色點(diǎn)線為所有曲線的平均。曲線與橫軸圍成的面積為當(dāng)年的實(shí)際產(chǎn)量,縱坐標(biāo)為各溫度貢獻(xiàn)的產(chǎn)量,生育期內(nèi)所有溫度貢獻(xiàn)的產(chǎn)量和即為當(dāng)年的實(shí)際產(chǎn)量。由圖可以看出,溫度貢獻(xiàn)的產(chǎn)量隨著溫度的升高而升高,達(dá)到產(chǎn)量峰值后,隨著溫度的繼續(xù)升高其帶來(lái)的產(chǎn)量貢獻(xiàn)逐漸下降。拔節(jié)前、拔節(jié)后和全生育期分別在35.2℃、34.5℃、35.3℃時(shí)溫度帶來(lái)的產(chǎn)量貢獻(xiàn)為0。此后隨著溫度的繼續(xù)升高,溫度帶來(lái)的產(chǎn)量貢獻(xiàn)為負(fù)值,可以認(rèn)為氣溫對(duì)作物生長(zhǎng)和產(chǎn)量形成產(chǎn)生負(fù)效應(yīng)。因此,將該產(chǎn)量貢獻(xiàn)值為0時(shí)所對(duì)應(yīng)的溫度作為玉米各發(fā)育階段極端高溫的閾值。采用該方法識(shí)別出拔節(jié)后的極端高溫事件的閾值較拔節(jié)前要低,說(shuō)明夏玉米拔節(jié)后對(duì)高溫更為敏感。

        圖4 玉米不同發(fā)育階段溫度-產(chǎn)量積分函數(shù)圖像

        2.4 極端高溫日數(shù)的分布

        將研究期(1983-2012年)分成1983-1992、1993-2002和2003-2012年3個(gè)時(shí)期,分夏玉米拔節(jié)前和拔節(jié)后兩個(gè)時(shí)段統(tǒng)計(jì)山東省極端高溫發(fā)生日數(shù)的空間分布,結(jié)果見(jiàn)圖5。由圖中可見(jiàn),在各個(gè)時(shí)期、各生長(zhǎng)階段極端高溫的空間分布均呈現(xiàn)出沿海半島地區(qū)年均極端高溫日數(shù)較少,向西部平原地區(qū)逐漸增多,菏澤、濟(jì)陽(yáng)、濟(jì)寧為極端高溫高發(fā)地區(qū),膠州、萊陽(yáng)發(fā)生較少。在時(shí)期2和時(shí)期3,全省夏玉米拔節(jié)后極端高溫發(fā)生日數(shù)普遍高于拔節(jié)前;在時(shí)期1,東部地區(qū)拔節(jié)后極端高溫日數(shù)也較拔節(jié)前多。不同時(shí)期之間比較可以發(fā)現(xiàn),時(shí)期3的極端高溫日數(shù)較多的站點(diǎn)面積最大,時(shí)期1面積最小,多數(shù)站點(diǎn)極端高溫日數(shù)有隨時(shí)間增加的趨勢(shì)。

        圖5 不同時(shí)期山東省夏玉米不同發(fā)育階段年均極端高溫日數(shù)的分布

        3 結(jié)論與討論

        (1)基于積溫產(chǎn)量模型方法確定了山東省夏玉米拔節(jié)前、后兩個(gè)階段的極端高溫閾值分別為35.2℃和34.5℃,拔節(jié)后極端高溫閾值較低,可見(jiàn)該階段對(duì)高溫更為敏感,容易造成損失。與從玉米生理特征角度看,拔節(jié)后的開(kāi)花灌漿期對(duì)高溫最為敏感的研究結(jié)論相一致[17]。Schlenker等[14]基于積溫產(chǎn)量模型方法研究時(shí),將貢獻(xiàn)的產(chǎn)量開(kāi)始下降時(shí)的溫度定義為極端高溫閾值(溫度-產(chǎn)量積分函數(shù)曲線的頂點(diǎn)),但超過(guò)該值的溫度貢獻(xiàn)的產(chǎn)量仍為正數(shù),即溫度對(duì)產(chǎn)量形成仍是有益的。而本研究將貢獻(xiàn)的產(chǎn)量開(kāi)始為零時(shí)的溫度定義為極端高溫閾值,則更為合理。

        (2)山東省夏玉米生長(zhǎng)季極端高溫日數(shù)呈現(xiàn)出西部?jī)?nèi)陸地區(qū)較高,自西向東逐漸降低的趨勢(shì);研究期內(nèi),極端高溫日數(shù)隨時(shí)間呈增多的趨勢(shì);夏玉米拔節(jié)后極端高溫發(fā)生較拔節(jié)前更為頻繁。說(shuō)明山東省中西部?jī)?nèi)陸地區(qū)面臨更嚴(yán)重高溫災(zāi)害威脅,且在氣候變暖背景下有加重的趨勢(shì),拔節(jié)后更為嚴(yán)重。同時(shí),拔節(jié)后也是玉米受極端高溫影響的敏感時(shí)期,極端高溫給玉米生產(chǎn)帶來(lái)更嚴(yán)重的威脅。

        [1]IPCC.Climate change 2007:impacts,adaptation and vulnerability contribution of working group II to the fourth assessment report of the intergovernmental panel on climate change[M]. Cambridge,UK and New York,USA:Cambridge University Press,2007.

        [2]Zhang Y,Xu Y,Dong W,et al.A future climate scenario of regional changes in extreme climate events over China using the PRECIS climate model[J].Geophysical Research Letters, 2006,33:L24702.

        [3]徐金芳,鄧振鏞,陳敏.中國(guó)高溫?zé)崂宋:μ卣鞯难芯烤C述[J].干旱氣象,2009,27(2):163-168.

        Xu J F,Deng Z Y,Chen M.A summary of studying on characteristicsof high temperature and heat wave damage in China[J].Journal of Arid Meteorology,2009,27(2):163-168.(in Chinese)

        [4]鄭昌玲,宋迎波,張艷紅,等.2016年秋收作物生長(zhǎng)季農(nóng)業(yè)氣象條件評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)氣象,2017,38(3):196-200.

        Zheng C L,Song Y B,Zhang Y H,et al.The assessment of agro-meteorological condition for autumn harvesting crops in 2016[J].Chinese Journal of Agrometeorology,2017,38(3):196- 200.(in Chinese)

        [5]Lobell D B,Sibley A,Ortiz-Monasterio J I.Extreme heat effects on wheat senescence in India[J].Nature Climate Change,2012, 2:186-189.

        [6]葛非凡,毛克彪,蔣躍林,等.華東地區(qū)夏季極端高溫特征及其對(duì)植被的影響[J].中國(guó)農(nóng)業(yè)氣象,2017,38(1):42-51.

        Ge F F,Mao K B,Jiang Y L,et al.Extreme high-temperature in summer and its impacts on vegetation in east China[J].Chinese Journal of Agrometeorology,2017,38(1):42-51.(in Chinese)

        [7]李慶祥,黃嘉佑.對(duì)我國(guó)極端高溫事件閾值的探討[J].應(yīng)用氣象學(xué)報(bào),2011,22(2):138-144.

        Li Q X,Huang J Y.Threshold values on extreme high temperature events in China[J].Journal of Applied Meteorological Science, 2011,22(2):138-144.(in Chinese)

        [8]周夢(mèng)子,王會(huì)軍,霍治國(guó).極端高溫天氣對(duì)玉米產(chǎn)量的影響及其與大氣環(huán)流和海溫的關(guān)系[J].氣候與環(huán)境研究,2017,22 (2):134-148.

        Zhou M Z,Wang H J,Huo Z G.The influence of heat stress on maize yield and its association with atmospheric general circulation and sea surface temperature[J].Climatic and Environmental Research,2017,22(2):134-148.(in Chinese)

        [9]馮春梅,郭東,馮愛(ài)國(guó).不同溫度對(duì)玉米產(chǎn)量的影響[J].現(xiàn)代農(nóng)業(yè),2014,(6):48-49.

        Feng C M,Guo D,Feng A G.Impacts of different temperature on corn yield[J].Modern Agriculture,2014,(6):48-49.(in Chinese)

        [10]劉哲,喬紅興,趙祖亮.黃淮海夏播玉米花期高溫?zé)岷臻g分布規(guī)律研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(7):272-279.

        Liu Z,Qiao H X,Zhao Z L.Spatial distribution of high temperature stress at corn flowering stage in Huang-Huai- Hai plain of China[J].Transactions of the Chinese Society for Agricultural Machinery,2015,46(7):272-279.(in Chinese)

        [11]李德,孫義,孫有豐.淮北平原夏玉米花期高溫?zé)岷C合氣候指標(biāo)研究[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2015,23(8):1035-1044.

        Li D,Sun Y,Sun Y F.Use of integrated climatic index to determine high temperature damage to summer maize at florescence in the Huaibei Plain[J].Chinese Journal of Eco-Agriculture,2015,23(8):1035-1044.(in Chinese)

        [12]尹小剛,王猛,孔箐鋅,等.東北地區(qū)高溫對(duì)玉米生產(chǎn)的影響及對(duì)策[J].應(yīng)用生態(tài)學(xué)報(bào),2015,26(1):186-198.

        Yin X G,Wang M,Kong Q X,et al.Impacts of high temperature on maize production and adaptation measures in Northeast China[J].Chinese Journal of Applied Ecology,2015,26(1): 186-198.(in Chinese)

        [13]Lobell D B,Hammer G L,McLean G,et al.The critical role of extreme heat for maize production in the United States[J]. Nature Climate Change,2013,(3):497-501.

        [14]Schlenker W,Roberts M J.Nonlinear temperature effects indicate severe damages to U.S.crop yields under climate change[J]. Proceedings of the National Academy of Sciences,2009,106 (37):15594-15598.

        [15]Hawks J,Ritchie J T.Modeling plant and soil systems[M]. American Society of Agronomy,Madison,1991:5-29.

        [16]姜會(huì)飛,溫德永,李楠.利用正弦分段法模擬氣溫日變化[J].氣象與減災(zāi)研究,2010,33(3):61-65.

        Jiang H F,Wen D Y,Li N.A new simulation method for the diurnal variation of temperature-sub-sine simulation[J]. Meteorology and Disaster Reduction Research,2010,33(3): 61-65.(in Chinese)

        [17]郭慶法,王慶成,汪黎明.中國(guó)玉米栽培學(xué)[M].上海:上??茖W(xué)技術(shù)出版社,2004.

        Guo Q F,Wang Q C,Wang L M.Maize culture science in China[M].Shanghai:Shanghai Science and Technology Press, 2004.(in Chinese)

        Determination of Extreme High Temperature Thresholds before and after SummerCorn Jointing Stage in Shandong Based on Accumulated Temperature-Yield Model

        ZHANG Qi1, 2, TANG Jie1, FENG Yi-chun1, ZHANG Yi-xuan1, CAI Cheng-liang1

        (1. College of Applied Meteorology, Nanjing University of Information Science and Technology, Nanjing 210024, China; 2. Jiangsu Key Laboratory of Agricultural Meteorology, Nanjing 210024)

        By using accumulated temperature-yield model, the hourly temperature was simulated based on daily temperature from 10 meteorological stations in Shandong province from 1983 to 2012. Then the extreme high temperature thresholds both before and after summer corn jointing stage were determined, and the spatial and temporal distribution of extreme high temperature events was analyzed. The results showed that temperature were in 17-35℃ during summer corn growing season, and the duration of each temperature after jointing was larger than that of before jointing. The proportional coefficient of effective accumulated temperature converting into actual yield increased during 1983-2012. The threshold of extreme high temperature before and after summer corn jointing stage was 35.2℃ and 34.5℃ respectively. The regional extreme high temperature days decreased from western to eastern in Shandong province, which was higher after jointing stage than that of before jointing stage. The extreme high temperature days increased gradually during the research period.

        Extreme high temperature threshold; Accumulated temperature-yield model; Summer corn; Jointing stage

        10.3969/j.issn.1000-6362.2017.12.006

        張琪,唐婕,馮一淳,等.基于積溫產(chǎn)量模型確定山東夏玉米拔節(jié)前后的極端高溫閾值[J].中國(guó)農(nóng)業(yè)氣象,2017,38(12):795-800

        2017-06-04

        江蘇省基礎(chǔ)研究計(jì)劃(自然科學(xué)基金)青年基金項(xiàng)目(BK20150898);公益性行業(yè)(氣象)科研專項(xiàng)(GYHY201506001-06);國(guó)家自然科學(xué)基金項(xiàng)目(41501553;41571491);公益性行業(yè)(氣象)科研專項(xiàng)儲(chǔ)備項(xiàng)目(2015001)

        張琪(1986-),女,博士,研究方向?yàn)檗r(nóng)業(yè)氣象災(zāi)害風(fēng)險(xiǎn)評(píng)估。E-mail:zhangq861206@126.com

        猜你喜歡
        有效積溫日數(shù)夏玉米
        漢江上游漢中區(qū)域不同等級(jí)降水日數(shù)的氣候變化特征分析
        綠色科技(2022年16期)2022-09-15 03:04:46
        有效積溫與不同供氮水平夏玉米干物質(zhì)和氮素積累定量化研究
        收斂式有效積溫與馬鈴薯植株性狀變化的相關(guān)性
        天津市濱海新區(qū)塘沽地域雷暴日數(shù)變化規(guī)律及特征分析
        天津科技(2020年2期)2020-03-03 05:09:48
        貴州省中東部水稻有效積溫對(duì)氣候變化的響應(yīng)
        貴州省中東部水稻有效積溫對(duì)氣候變化的響應(yīng)
        小麥?zhǔn)崭钪?如何種植夏玉米才能高產(chǎn)
        夏玉米高產(chǎn)的關(guān)鍵栽培技術(shù)措施
        海南省雷暴日數(shù)年代際變化特征
        ESSENTIAL NORMS OF PRODUCTS OF WEIGHTED COMPOSITION OPERATORS AND DIFFERENTIATION OPERATORS BETWEEN BANACH SPACES OF ANALYTIC FUNCTIONS?
        国产精品亚洲精品日产久久久| 色哟哟网站在线观看| 中文字幕在线亚洲日韩6页手机版| 国产美女精品AⅤ在线老女人| 男人的天堂av你懂得| 久久久久亚洲av成人网人人软件| 久久久久久免费毛片精品| 亚洲精品国产品国语在线app| 久久婷婷国产色一区二区三区| 亚洲综合另类小说色区| 无码中文字幕日韩专区视频| 国产小屁孩cao大人| 午夜国产精品一区二区三区| 99国产精品久久久久久久成人热| 中文字幕日韩精品无码内射| 一本大道久久a久久综合| 老司机在线免费视频亚洲| 国产成人午夜高潮毛片| 久久久www成人免费无遮挡大片 | 中文亚洲欧美日韩无线码| 国产综合自拍| 亚洲一区二区三区厕所偷拍| 狂猛欧美激情性xxxx大豆行情| 精品国产人妻一区二区三区| 亚洲—本道中文字幕久久66| 青青草免费观看视频免费| 性无码一区二区三区在线观看| 老湿机香蕉久久久久久| 国产av天堂亚洲国产av麻豆| 一本色道久久88加勒比| 人妻少妇偷人精品无码| 国产av天堂成人网| 女同中文字幕在线观看| 2019nv天堂香蕉在线观看| 欧美人与动牲交a欧美精品| 国产精品无码久久久久久蜜臀AV| 国产免费操美女逼视频| 亚洲熟妇无码一区二区三区导航| 日韩中文无线码在线视频观看| 成人综合激情自拍视频在线观看| 久青草影院在线观看国产 |