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

    寧夏中部干旱帶玉米生育期熱量條件時空特征分析

    2016-02-13 06:15:06陳璐劉靜王連喜李琪許晨純
    關鍵詞:海原積溫寧夏

    陳璐,劉靜,王連喜,李琪,許晨純

    (1.江蘇省農(nóng)業(yè)氣象重點實驗室,江蘇南京210044;2.南京信息工程大學應用氣象學院,江蘇南京210044;3.寧夏氣象防災減災重點實驗室,寧夏銀川750002)

    寧夏中部干旱帶玉米生育期熱量條件時空特征分析

    陳璐1,2,劉靜3,王連喜1,2,李琪1,2,許晨純1,2

    (1.江蘇省農(nóng)業(yè)氣象重點實驗室,江蘇南京210044;2.南京信息工程大學應用氣象學院,江蘇南京210044;3.寧夏氣象防災減災重點實驗室,寧夏銀川750002)

    利用寧夏中部干旱帶鹽池、同心及海原三地1988—2013年逐日氣象資料建立玉米全生育期及各主要生育階段溫度適宜度模型,并對寧夏中部干旱帶玉米生育期溫度適宜度及≥10℃積溫進行時空變化分析。結(jié)果表明:寧夏中部干旱帶玉米全生育期溫度適宜度以及≥10℃積溫均呈波動上升趨勢,三地的溫度適宜度從大到小依次為同心、鹽池、海原,數(shù)值分別為:0.966,0.911和0.726。溫度適宜度在播種~出苗,出苗~拔節(jié),拔節(jié)~抽雄,開花吐絲~成熟階段均呈上升趨勢,其中拔節(jié)~抽雄階段上升速率最快,速率為0.03·10a-1。鹽池玉米的熱量關鍵期為抽雄~開花吐絲階段,同心為播種~出苗階段,而海原為拔節(jié)~抽雄階段。

    玉米生育期;溫度適宜度;積溫;熱量時空變化;寧夏中部干旱帶

    近百年來,以全球氣候變暖為主要特征的氣候變化對我國農(nóng)業(yè)生態(tài)環(huán)境、作物生長及產(chǎn)量造成了重大影響。尤其是我國西北地區(qū)生態(tài)環(huán)境脆弱,農(nóng)業(yè)生產(chǎn)對氣候變化尤為敏感[1]。

    玉米是我國三大糧食作物之一,已有的研究表明:氣候變暖導致玉米種植北界北移,種植面積擴大,但也加快了玉米生育進程,縮短生育期而導致玉米減產(chǎn)[2-4]。對寧夏中部干旱帶的研究表明:氣溫對玉米生長影響較大,總體來講,生育前期的高溫與后期低溫都將導致玉米減產(chǎn);降水量是影響本地區(qū)春玉米生產(chǎn)潛力的制約因子[5-6]。氣候環(huán)境適宜度直接制約著玉米的生長發(fā)育,影響當?shù)赜衩桩a(chǎn)量。目前,氣候環(huán)境適宜度分析已經(jīng)逐步從定性化發(fā)展為定量化[7-10],氣候適宜度通過模糊數(shù)學中隸屬度函數(shù)的方法把氣候因子(光照、溫度及降水)數(shù)量變化轉(zhuǎn)化為對作物生長發(fā)育、產(chǎn)量形成的適宜程度,國內(nèi)不少專家對不同區(qū)域尺度玉米氣候適宜度進行了研究分析,但研究區(qū)域多集中在華中、華北及東北三省,對西北尤其是寧夏自治區(qū)的研究較少[11-15]。

    寧夏中部干旱帶位于寧夏中部地區(qū),包括鹽池縣、同心縣、海原縣等地,海拔1 300~2 100m,是我國北方地區(qū)與天然草地牧區(qū)接壤的過渡帶[16]。玉米生長喜溫喜光,在生長關鍵期需水量大。寧夏中部干旱帶日照充足,降水少,在當?shù)貙嶋H生產(chǎn)中,玉米需水關鍵期通常通過補充灌溉來保證玉米生長發(fā)育的水分需求,因此,在該地區(qū)的大田生產(chǎn)中,日照與降水對玉米生長影響不大。本文以當?shù)赜衩椎臒崃織l件作為研究對象,分析1988—2013年寧夏中部干旱帶玉米溫度適宜度及≥10℃積溫的時空分布特征,以對當?shù)赜衩追N植適應氣候變化決策建議提供依據(jù)。

    1 材料與方法

    1.1 材料來源

    本文氣象資料來自同心縣、鹽池縣和海原縣三個氣象站觀測的1988—2013年的逐日觀測資料。產(chǎn)量數(shù)據(jù)包括同心縣、鹽池縣及海原縣1988—2013年玉米單產(chǎn),由寧夏自治區(qū)統(tǒng)計局提供提供。

    取這三個站點溫度適宜度平均值作為該地區(qū)溫度適宜度的時間序列值,各個縣的單產(chǎn)逐年平均值作為該地區(qū)單產(chǎn)的時間序列值。

    1.2 研究方法

    1.2.1 分離氣象產(chǎn)量玉米作物產(chǎn)量組成包括趨勢產(chǎn)量Yt和氣象產(chǎn)量Yw,所以

    利用鹽池縣、同心縣及海原縣統(tǒng)計產(chǎn)量,運用11 a直線滑動平均法提取趨勢產(chǎn)量Yt,再運用公式(2)分離氣象產(chǎn)量。該方法為滑動平均法和線性回歸法的結(jié)合,以11 a為步長,在每個11 a將產(chǎn)量看作線性函數(shù),隨著時間滑動,直線不斷改變位置,后延滑動,依次建立各11 a內(nèi)的線性回歸方程,并得到模擬值,然后每一年的線性回歸模擬值的平均值就是當年的趨勢產(chǎn)量。用該方法既不用主觀判定產(chǎn)量歷史演變的曲線類型,也不損失樣本量[18]。分離得到的氣象產(chǎn)量將代表玉米產(chǎn)量中受氣象條件影響的那部分產(chǎn)量,即表示對氣象條件變化的響應情況。

    1.2.2 溫度適宜度計算方法根據(jù)相關文獻研究[19-20],結(jié)合寧夏地區(qū)的實際情況,各個生育階段玉米的溫度適宜度函數(shù)可以表達為:

    式中,S(T)為平均氣溫為T時溫度適宜度;T0為玉米發(fā)育的最適氣溫(℃),不同發(fā)育階段最適氣溫不同;T1為玉米發(fā)育的最低氣溫(℃),低于這一氣溫,發(fā)育速率為0;T2為玉米發(fā)育的最高氣溫(℃),高于這一氣溫,發(fā)育速率為0。根據(jù)文獻對不同生育階段選取玉米三基點溫度[14-15](表1)。

    表1 玉米不同生育階段T0,T1及T2/℃Table 1 T0,T1and T2values inmaize growth period

    玉米在各個生育階段對熱量條件的要求不同,同時各個生育階段熱量條件對玉米生長發(fā)育及產(chǎn)量的影響也不同。所以,為客觀反映不同生育階段熱量條件對玉米生長發(fā)育的影響強度,需要為玉米各個生育階段的溫度適宜度設定權(quán)重。本文采用加權(quán)平均法,確定各個生育階段溫度適宜度的權(quán)重[20],結(jié)果見表2。

    式中,bi為第i個生育期溫度適宜度的權(quán)重系數(shù);ai為第i個生育期平均氣溫與玉米氣象產(chǎn)量的相關系數(shù);Sm(T)為全生育期玉米的溫度適宜度;Si(T)為第i個生育期的溫度適宜度。

    2 結(jié)果與分析

    根據(jù)多年觀測資料,寧夏中部干旱帶玉米種植在4—9月份,其中播種~出苗期大致在4月中旬—5月上旬,出苗~拔節(jié)期大約在5月中旬—6月中旬,拔節(jié)~抽雄階段在6月下旬—7月上旬,抽雄~開花吐絲階段時間較短,在7月中旬,開花吐絲~成熟階段在7月中旬—9月上旬。

    播種~拔節(jié)階段為玉米生長的苗期,該階段根系發(fā)育快,地上部生長較慢,對環(huán)境要求不高,但怕冷不怕熱,怕澇不怕旱,輕度干旱有益于根系的發(fā)育和下扎。拔節(jié)~抽雄階段為穗期,該階段是玉米生長發(fā)育最旺盛的時期,對水分要求迫切;抽雄前10 d至抽雄后20 d,是水分臨界期、穗分化及開花期對水分要求敏感。抽雄~成熟階段為花粒期,此時葉片、株高已定,營養(yǎng)生長基本停止,該階段對水分需求量減少。

    寧夏中部干旱帶1988—2003年玉米氣象產(chǎn)量呈波動持平狀態(tài),在1988—1997年這10年中呈增加趨勢,2000年氣象產(chǎn)量為28年中最低,2005—2013年中部干旱帶玉米氣象產(chǎn)量呈平穩(wěn)下降趨勢,但趨勢不顯著。鹽池、同心和海原玉米氣象產(chǎn)量變化情況與中部干旱帶整體相似,同心玉米氣象產(chǎn)量變化幅度較小,鹽池于1995年氣象產(chǎn)量表現(xiàn)最高;海原玉米氣象產(chǎn)量變化幅度較為劇烈,并且于1998年和2000年分別出現(xiàn)氣象產(chǎn)量最高和最低值。

    三個地區(qū)氣象產(chǎn)量變化情況大體一致,但最值出現(xiàn)年份不同,因為三個地區(qū)所在的地理位置和當?shù)貧庀髼l件存在差異(見圖1)。

    2.1 玉米生育期平均氣溫時空分布特征

    寧夏中部干旱帶玉米生育期從4月中旬開始,到9月上旬結(jié)束,用鹽池、同心及海原三地氣溫的平均值代表寧夏中部干旱帶的平均氣溫,1988—2013年全生育期平均氣溫呈周期性波動,并以0.35℃· 10a-1(R=0.487,P<0.01)的速率升高。從空間分布上看,生育期內(nèi)同心的平均氣溫最高而海原最低。三地均呈周期性波動增溫趨勢。其中,鹽池以0.22℃·10a-1速率增溫,同心以0.40℃·10a-1(R=0.519,P<0.01)速率增溫,海原以0.44℃·10a-1(R=0.563,P<0.01)速率增溫,可以看出寧夏中部干旱帶在玉米生育期中部氣溫最高,南部氣溫低于北部,但南部增溫速率最高(見圖2)。

    平均氣溫的時間變化在不同生育期趨勢不同,在抽雄~開花階段平均氣溫呈減小趨勢,但趨勢不明顯;在播種~出苗階段,中部干旱帶平均氣溫以0.81℃·10a-1(R=0.458,P<0.05)的速率增溫,同心增溫最快,海原增溫速率較慢。出苗~拔節(jié)階段中部干旱帶平均氣溫以0.48℃·10a-1(R=0.449,P<0.05)的速率增溫,增溫趨勢明顯;其中同心增溫速率最快,其次是海原,但鹽池增溫不明顯。在拔節(jié)~抽雄階段,海原增溫明顯,增溫速率為0.61℃· 10a-1(R=0.390,P<0.05)??傮w來說,播種~出苗階段增溫速率最快,抽雄~開花吐絲階段平均氣溫呈減小趨勢(表3)。

    表2 各個生育期溫度適宜度權(quán)重系數(shù)(bi)Table 2 Weight coefficient in different growth period(bi)

    圖1 中部干旱帶玉米氣象產(chǎn)量時間變化Fig.1 Variations of meteorological yield of maize in middle arid region

    圖2 寧夏中部干旱帶玉米全生育期平均氣溫時間變化Fig.2 Variations of average temperature in middle arid region of whole growth period

    表3 近26年玉米各生育期平均氣溫增溫速率/(℃·10a-1)Table 3 Increase rate of average temperature in different growth period in the last 26 years

    2.2 玉米生育期≥10℃積溫時空變化特征

    ≥10℃積溫(即活動積溫)反映了生育期內(nèi)玉米對熱量的要求。從圖2可以看出,1988—2013年近26年寧夏中部干旱帶不同地區(qū)玉米全生育期內(nèi)≥10℃積溫時間變化趨勢大致相同,中部干旱帶≥10℃積溫呈波動上升趨勢,以69.3℃·10a-1(R=0.515,P<0.01)的速率增加。海原為三地積溫增溫速率最高的區(qū)域,其增溫速率為88.5℃·10a-1(R=0.542,P<0.01);其次為同心,增溫速率為75.2℃·10a-1(R=0.549,P<0.01);鹽池增溫速率為44.0℃·10a-1,增溫趨勢不明顯(圖3)。

    利用MK-突變檢驗方法對寧夏中部干旱帶玉米全生育期≥10℃積溫的時間序列進行分析,結(jié)果表明:寧夏中部干旱帶玉米全生育期≥10℃積溫在1994年左右發(fā)生突變,積溫由增加不明顯突變?yōu)樵黾用黠@。鹽池全生育期積溫于1991年由增加不明顯突變?yōu)樵黾用黠@,但僅在2000—2002年增溫達到顯著性水平;積溫于2007年由增大趨勢突變?yōu)闇p少趨勢。同心及海原兩地玉米全生育期≥10℃積溫突變情況相似,在1994年突變?yōu)樵黾用黠@,并且從2006年開始積溫增溫趨勢達到顯著性水平(圖4)。

    圖3 寧夏中部干旱帶玉米全生育期≥10℃積溫時間變化Fig.3 Variations of active accumulated temperature in whole growth period inmiddle arid region

    表4的統(tǒng)計結(jié)果表明:近26年來玉米不同生育階段≥10℃積溫變化趨勢表現(xiàn)各不相同。播種~出苗及出苗~拔節(jié)階段≥10℃積溫波動增加,中部干旱帶在播種~出苗階段增溫速率較大,為48.1℃· 10a-1(R=0.542,P<0.01),其中海原增溫速率最大,為53.2℃·10a-1(R=0.563,P<0.01),其次為同心,增溫速率為47.8℃·10a-1(R=0.527,P<0.01),最后為鹽池,增溫速率為46.3℃·10a-1(R=0.483,P<0.05)。出苗~拔節(jié)階段中部干旱帶增溫速率為19.4℃·10a-1(R=0.404,P<0.05),其中同心增溫速率最大,其次為海原;雖然鹽池在該階段也表現(xiàn)出增溫趨勢,但趨勢不明顯。中部干旱帶在拔節(jié)~抽雄階段積溫呈波動增加趨勢,但除海原以12.3℃·10a-1(R=0.386,P<0.05)明顯增溫外,其它地區(qū)增溫趨勢不明顯。抽雄~開花吐絲及開花吐絲~成熟階段中部干旱帶大部分地區(qū)積溫呈波動下降趨勢,趨勢不明顯。

    2.3 全生育期溫度適宜度時空分布

    中部干旱帶自1988—2013年,玉米全生育期溫度適宜,適宜度為0.869,且呈波動上升趨勢,上升速率為0.012·10a-1,波動較為平緩。從距平上看,中部干旱帶玉米全生育期溫度適宜度距平基本在±0.1內(nèi)波動,1988—1999年正負距平交替出現(xiàn),進入21世紀后,正距平增多,以正距平為主。

    空間分布上,鹽池、同心及海原三地年際變化情況基本一致,其中同心全生育期溫度適宜度最高,平均為0.966,呈波動上升趨勢,上升速率為0.005· 10a-1。距平波動范圍較小,1988—1999年負距平較多,進入21世紀后正距平增多,轉(zhuǎn)變?yōu)橐哉嗥綖橹?。鹽池溫度適宜度次之,多年平均為0.911,呈波動下降趨勢,下降速率為0.005·10a-1。年際波動情況與同心相似,但波動幅度較同心大。海原全生育期溫度適宜度較低,平均為0.726,多年呈波動上升趨勢,上升速率為三地最快,為0.037·10a-1。海原溫度適宜度波動最劇烈,其距平在±0.2間波動,且出現(xiàn)較強距平現(xiàn)象頻繁,進入21世紀后正距平增多,由以負距平為主轉(zhuǎn)變?yōu)橐哉嗥綖橹鳎▓D5)。

    圖4 玉米全生育期≥10℃積溫MK-突變檢驗Fig.4 MK-mutation test of active accumulated temperature during whole growth period

    表4 近26年玉米各生育期≥10℃積溫增溫速率/(℃·10a-1)Table 4 Increase rate of active accumulated temperature in different growth period in the last26 years

    圖5 玉米全生育期溫度適宜度變化Fig.5 Variations of temperature suitability degrees during thewhole growth period

    2.4 不同生育階段溫度適宜度時空分布

    2.4.1 播種~出苗階段溫度適宜度時空分布寧夏中部干旱帶玉米播種~出苗階段溫度適宜度較高,平均為0.965,波動較平緩,以0.004·10a-1的速率上升。距平在-0.15~0.05之間波動,且正距平較多。由于1991年該生育階段中部干旱帶平均氣溫為26年最低值,其溫度適宜度明顯偏低,負距平為26年來最大負距平值。

    鹽池、同心及海原三地在該生育階段溫度適宜度均在0.96以上,其中鹽池溫度適宜度最高,為0.971,波動平緩,以0.002·10a-1的速率上升。除較冷的1991年外,距平均在-0.1~0.05間小幅波動,21世紀后正距平逐漸增多。同心次之,平均為0.963,以0.013·10a-1速率減少;從2004年開始,負距平增多,且有增大的趨勢。海原為0.962,波動較劇烈,以0.021·10a-1的速率上升;在最冷的1991年,負距平較強,1992—2000年距均為正值,21世紀后負距平增多(圖6)。

    2.4.2 出苗~拔節(jié)階段溫度適宜度時空分布近26年,寧夏中部干旱帶玉米出苗~拔節(jié)階段溫度適宜度較好,平均適宜度為0.963,以0.014·10a-1(R=0.543,P<0.01)速率平緩上升。距平在-0.05~0.04間波動,21世紀以前以負距平為主,較強負距平出現(xiàn)在1988年,21世紀之后變?yōu)橐哉嗥綖橹?,距平值有上升趨勢?/p>

    鹽池與同心該生育階段溫度適宜度較高。同心適宜度波動平緩,為最高0.992,且以0.008·10a-1(R=0.609,P<0.01)速率呈顯著上升。鹽池較好,平均為0.988,波動幅度最小,以0.004·10a-1速率上升。海原最低,為0.910,波動較劇烈,以0.031· 10a-1(R=0.549,P<0.01)速率顯著上升。三地距平變化趨勢與中部干旱帶整體變化一致(圖7)。

    2.4.3 拔節(jié)~抽雄階段溫度適宜度時空分布寧夏中部干旱帶玉米拔節(jié)~抽雄階段溫度適宜度為0.880,以0.03·10a-1速率劇烈波動上升。1994年適宜度最低,2008年適宜度最高,距平變化與出苗~拔節(jié)階段變化相同。

    該階段同心溫度適宜度最高,為0.961,近26年以0.012·10a-1速率平緩波動上升,2004—2010年連續(xù)正距平,1994年氣溫偏低,故溫度適宜度最低,出現(xiàn)較強負距平。鹽池次之,為0.957,以0.008·10a-1速率平緩波動上升,其中2003年最低,出現(xiàn)較強負距平。海原最低,為0.720,以0.072·10a-1速率快速波動上升趨勢,且波動幅度較大;溫度適宜度波動情況與整體波動情況保持一致(圖8)。

    2.4.4 抽雄~開花吐絲階段溫度適宜度時空分布

    中部干旱帶溫度適宜度在玉米抽雄~開花吐絲階段較低,適宜度為0.780,近26年來呈波動下降趨勢,波動幅度較大,但0.007·10a-1的下降速率緩慢。距平波動幅度較大,變化趨勢與出苗~拔節(jié)階段相似,1992—1995年與2008年出現(xiàn)較強負距平,其中1995年溫度適宜度最低;2000—2003年及2005—2006年出現(xiàn)較強正距平。

    該階段同心溫度適宜度最高,為0.861,以0.008 ·10a-1速率緩慢波動下降。鹽池次之,為0.840,以0.014·10a-1速率波動下降。海原最低,為0.600,整體呈波動持平狀態(tài),波動激烈。三地適宜度波動情況與中部干旱帶整體波動情況一致(圖9)。

    2.4.5 開花吐絲~成熟階段溫度適宜度時空分布

    中部干旱帶玉米抽雄~開花吐絲階段溫度適宜度最低,為0.897,以0.012·10a-1速率波動上升,波動幅度較小。其距平變化與出苗~拔節(jié),抽雄~開花吐絲階段變化趨勢相同;在1993—1995年及1999年出現(xiàn)較強負距平。

    圖6 播種~出苗階段溫度適宜度變化Fig.6 Variations of temperature suitability degrees in the sowing to emergence period

    圖7 出苗~拔節(jié)階段溫度適宜度變化Fig.7 Variations of temperature suitability degrees in the emergence to jointing period

    圖8 拔節(jié)~抽雄階段溫度適宜度變化Fig.8 Variations of temperature suitability degrees of jointing to tasseling period

    該生育階段同心溫度適宜度最高,平均0.978,以0.007·10a-1速率平緩波動上升。次高為鹽池,平均0.953,以0.003·10a-1速率平緩波動上升。海原較低,為0.762,以0.012·10a-1速率劇烈波動上升。三地溫度適宜度時間波動情況與中部干旱帶整體波動情況一致(圖10)。

    2.5 熱量指標與玉米產(chǎn)量的關系

    玉米生長不僅受氣象因子影響,還受土壤類型、作物品種、病蟲害、農(nóng)業(yè)管理等其它因素影響。本文以氣象產(chǎn)量作為對象,討論寧夏中部干旱帶玉米生育期溫度適宜度及≥10℃積溫對玉米氣象產(chǎn)量的影響。

    表5為不同生育期溫度適宜度與氣象產(chǎn)量之間的相關系數(shù),結(jié)果表明,在鹽池地區(qū)玉米抽雄~開花吐絲階段為玉米生長的關鍵時期,溫度的高低直接影響當?shù)赜衩咨L發(fā)育,通過計算可知:該階段平均氣溫每升高1℃,玉米氣象產(chǎn)量增加10.165 kg· hm-2。鹽池在該生育階段玉米產(chǎn)量與溫度適宜度呈顯著正相關,溫度適宜度每上升0.01,氣象產(chǎn)量將增加1.731 kg·hm-2(R=0.530,P<0.01);并且全生育期溫度適宜度與氣象產(chǎn)量呈顯著正相關。同心播種~出苗階段為玉米熱量的關鍵時期,溫度適宜度每增加0.01,氣象產(chǎn)量將增加4.389 kg·hm-2(R=0.442,P<0.05)。海原玉米熱量條件關鍵期為拔節(jié)~抽雄階段,溫度適宜度每上升0.01,氣象產(chǎn)量將增加1.698 kg·hm-2(R=0.326,P<0.1)。

    圖9 抽雄~開花吐絲階段溫度適宜度變化Fig.9 Variations of temperature suitability degrees of tasseling to flowering silking period

    圖10 開花吐絲~成熟階段溫度適宜度變化Fig.10 Variations of temperature suitability degrees of flowering silking tomature period

    表5 玉米溫度適宜度與氣象產(chǎn)量的相關系數(shù)Table 5 Correlative coefficient between temperature suitability and meteorological yield ofmaize

    溫度對玉米的影響表現(xiàn)在兩方面:一是直接影響玉米的生長發(fā)育及產(chǎn)量;二是影響玉米的生育進程及發(fā)育時間,從而影響產(chǎn)量。玉米拔節(jié)~抽雄階段溫度過高,會縮短幼穗各分化期,減少分化小穗和小花數(shù),果穗變小。該階段溫度適宜度對海原玉米產(chǎn)量起明顯正作用,這段時間海原氣溫呈由低溫向最適溫度T0顯著增加趨勢,所以該階段海原溫度適宜度上升,有利于提高玉米產(chǎn)量。抽雄到開花吐絲階段要求平均氣溫25℃~28℃,溫度適宜,授粉良好;如果溫度達到32℃~35℃,而相對濕度低于30%,則花粉將迅速失去生命力,花絲也將枯萎導致受精不良,造成禿尖。但寧夏中部干旱帶在該階段氣溫偏低,且低于最適溫度,故造成玉米禿尖率增大的主要原因是干旱。該階段鹽池溫度適宜度對玉米產(chǎn)量起顯著正作用,而鹽池氣溫呈降低趨勢,故該階段鹽池溫度適宜度呈下降趨勢,玉米產(chǎn)量降低。籽粒形成和灌漿期,要求氣溫20℃~24℃為宜,若氣溫過高或過低,將影響酶的活動,玉米將延期成熟,造成籽粒不飽滿,產(chǎn)量下降;如再遇上干旱,易受高溫逼熟,使千粒重下降[7]。

    計算不同生育階段≥10℃積溫與氣象產(chǎn)量之間的相關系數(shù)(表6),得到結(jié)果與溫度適宜度與氣象產(chǎn)量之間相關關系類似。即,鹽池在抽雄~開花吐絲階段,積溫與玉米氣象產(chǎn)量呈顯著正相關,積溫每升高1℃,氣象產(chǎn)量將增加1.017 kg·hm-2(R=0.453,P<0.05);海原在拔節(jié)~抽雄階段積溫每升高1℃,氣象產(chǎn)量將增加1.154 kg·hm-2(R=0.372,P<0.1)。

    海原玉米氣象產(chǎn)量對熱量條件的相關性并不高,原因可能為:海原大部分地區(qū)覆膜種植玉米,由于覆膜影響,在沒有考慮覆膜的增溫保墑作用對玉米生長造成的影響情況下,導致出現(xiàn)如此結(jié)果。

    表6 玉米不同生育階段≥10℃積溫與氣象產(chǎn)量的相關系數(shù)Table 6 Correlative coefficient between active accumulated temperature and meteorological yield of maize

    表6為玉米不同生育階段≥10℃積溫與氣象產(chǎn)量的相關系數(shù)計算結(jié)果,通過與表5的對比可以看出:雖然溫度適宜度與≥10℃積溫均為描述熱量條件的指標,但≥10℃積溫只考慮了玉米的生物學下限溫度,而忽略了上限溫度;并且對于不同生育階段及全生育期,考慮的下限溫度都相同,從而導致了在與氣象產(chǎn)量的相關性上出現(xiàn)差異。而根據(jù)已有研究,玉米不同生育階段下限溫度都不相同,所以與溫度適宜度比較,≥10℃活動積溫的涵蓋范圍太大。從與氣象產(chǎn)量的關系上看,溫度適宜度能夠更好地描述熱量條件與氣象產(chǎn)量之間的關系。

    3 結(jié)論

    1)寧夏中部干旱帶玉米全生育期近26年以0.35℃·10a-1的速率明顯增溫,其中增溫速率海原>同心>鹽池。不同生育階段平均氣溫變化趨勢不同,播種~出苗階段增溫明顯且速率最快,其次是出苗~拔節(jié)階段。抽雄~開花吐絲階段平均氣溫呈減少趨勢,其它生育階段有增溫趨勢,但趨勢均不明顯。

    2)寧夏中部干旱帶玉米全生育期≥10℃積溫近26年以69.3℃·10a-1明顯上升,增溫速率同心>海原,鹽池增溫不明顯。不同生育階段≥10℃積溫變化趨勢不同,播種~出苗階段增溫速率大于出苗~拔節(jié)階段,而抽雄~開花吐絲及開花吐絲~成熟階段≥10℃積溫呈下降趨勢。

    3)寧夏中部干旱帶玉米全生育期溫度適宜度多年呈波動上升趨勢,上升速率拔節(jié)~抽雄>出苗~拔節(jié)>開花吐絲~成熟>播種~出苗,在抽雄~開花吐絲階段溫度適宜度呈下降趨勢。從適宜度上,播種~出苗>出苗~拔節(jié)>開花吐絲~成熟>拔節(jié)~抽雄>抽雄~開花吐絲。

    4)鹽池玉米熱量關鍵期為抽雄~開花吐絲階段,同心為播種~出苗階段,海原為拔節(jié)~抽雄階段。與≥10℃積溫相比,溫度適宜度考慮了玉米生物學上限溫度,并對每個生育階段的上、下限溫度作了細化,能夠更好地描述熱量條件與氣象產(chǎn)量之間的關系。

    4 討論

    寧夏中部干旱帶溫度適宜度總體呈上升趨勢,且播種~出苗,出苗~拔節(jié)階段溫度適宜度較高,影響玉米的關鍵期為播種~出苗階段,抽雄~開花吐絲階段,這與已有對西北地區(qū)環(huán)境適宜度對旱作玉米的影響研究得到的溫度適宜度關鍵期相符[15]。從全生育期熱量條件的變化趨勢上可以看出,寧夏中部干旱帶熱量條件在增加,更加有利于玉米生長,但從各個生育階段來看,抽雄~成熟階段熱量條件有下降趨勢,這也與目前分析氣候變化特征的主要結(jié)論即夏季增溫不明顯相符[21],故在該區(qū)域的玉米播期調(diào)整、品種選擇及玉米生產(chǎn)管理中應充分考慮到該區(qū)域玉米生育期內(nèi)熱量條件的變化與不同,以合理利用氣候資源,適應氣候變化。而溫度適宜度較低的海原,其適宜度上升速率較快,在不考慮降水的前提下,未來氣候變化將使海原更適合種植玉米,并提高其單產(chǎn)。由于玉米生育期內(nèi)熱量條件的增加,寧夏中部干旱帶玉米種植范圍將擴大,而根據(jù)姚小英等人[15]對西北地區(qū)的研究,結(jié)合西北地區(qū)暖干化趨勢,預計玉米生育期水分適宜度將下降,所需灌溉量將增大,所以玉米種植風險將有所增加,由熱量增加帶來的有利條件可能將無法抵消缺水所造成的影響。

    本文只討論了寧夏中部干旱帶的熱量條件,而并未涉及降水及日照。在寧夏中部干旱帶,光照充沛,但干旱少雨,故玉米種植需依賴補充灌溉,補充灌溉可能出現(xiàn)灌溉不及時或未達到理想灌溉量而影響玉米產(chǎn)量的結(jié)果,同時沒有涉及溫度、水、日照和水汽等變量的綜合作用對玉米的影響以及不同玉米品種對氣候條件的不同要求。同時,由于對中部干旱帶玉米種植觀測資料的不足,本文對玉米生育期的劃分主要是根據(jù)引黃灌區(qū)玉米觀測資料以旬為單位進行了生育期的劃分,使得生育期表達不夠精確。此外,覆膜作為重要的農(nóng)業(yè)技術(shù)在玉米生產(chǎn)中得到應用,今后應在量化單一氣象要素和覆膜技術(shù)的增溫效應對產(chǎn)量的影響方面以及更為適用的適宜度計算公式方面加強研究,以進一步提高適宜度的表達準確度。

    [1]王馥棠,趙宗慈,王石立,等.氣候變化對農(nóng)業(yè)生態(tài)的影響[M].北京:氣象出版社,2003:96-126.

    [2]張強,鄧振鏞,趙映東,等.全球氣候變化對我國西北地區(qū)農(nóng)業(yè)的影響[J].生態(tài)學報,2008,(3):1210-1219.

    [3]郭慶海.中國玉米主產(chǎn)區(qū)的演變與發(fā)展[J].玉米科學,2010,18(1):139-145.

    [5]Watson R T,Albritton D L,Barker T.Climate Change 2001—Synthesis Report.United States of America:Cambridge University Press,2001:46.

    [6]賴榮生,余海龍,黃菊瑩.寧夏中部干旱帶氣候變化及其對春玉米氣候生產(chǎn)潛力的影響[J].中國農(nóng)業(yè)大學學報,2014,19(3):108-114.

    [7]馬樹慶.吉林省農(nóng)業(yè)氣候研究[M].北京:氣象出版社,1994:33.

    [8]白永平.西北地區(qū)(甘寧青)農(nóng)業(yè)生態(tài)氣候資源量化與評價[J].自然資源學報,2000,15(3):218-224.

    [9]趙峰,千懷遂,焦士興.農(nóng)作物氣候適宜度模型研究——以河南省冬小麥為例[J].資源科學,2003,25(6):77-82.

    [10]趙峰,千懷遂.全球變暖影響下農(nóng)作物氣候適宜性研究進展[J].中國生態(tài)農(nóng)業(yè)學報,2004,12(2):134-137.

    [11]郝君明,趙軍,劉紅兵,等.基于FIS與GIS方法的山西省玉米氣候適宜性評價[J].安徽農(nóng)業(yè)科學,2008,36(23):10129-10130.

    [12]劉偉昌,陳懷亮,余衛(wèi)東,等.基于氣候適宜度指數(shù)的冬小麥動態(tài)產(chǎn)量預報技術(shù)研究[J].氣象與環(huán)境科學,2008,(31):21-24.

    [13]賈建英,郭建平.東北地區(qū)近46年玉米氣候資源變化研究[J].中國農(nóng)業(yè)氣象,2009,30(3):302-307.

    [14]宮麗娟,王晨軼,朱海霞,等.玉米氣候適宜度變化分析[J].玉米科學,2013,21(5):140-146.

    [15]姚小英,蒲金涌,姚茹莘,等.氣候暖干化背景下甘肅旱作區(qū)玉米氣候適宜性變化[J].地理學報,2011,66(1):59-67.

    [16]何永平,徐秀梅.寧夏中部干旱帶氣候變化對玉米產(chǎn)量的影響—以鹽池縣、同心縣為例[J].寧夏農(nóng)林科技,2012,53(4):1-3,17.

    [17]魏瑞江,宋迎波,王鑫.基于氣候適宜度的玉米產(chǎn)量動態(tài)預報方法[J].應用氣象學報,2009,20(5):622-627.

    [18]范磊,鄭國清,趙全志,等.河南省水稻氣象災害因子分析及主要災害發(fā)生區(qū)劃研究[J].中國農(nóng)業(yè)資源與區(qū)劃,2013,34(6);162-168.

    [19]千懷遂,焦士興,趙峰.全球氣候變化對作物氣候適宜性的影響—以河南省冬小麥為例[J].農(nóng)業(yè)現(xiàn)代化研究,2004,25(2):106-110.

    [20]魏瑞江,張文宗,康西言,等.河北省冬小麥氣候適宜度動態(tài)模型的建立及應用[J].干旱地區(qū)農(nóng)業(yè)研究,2007,25(6):5-15.

    [21]陳曉光,蘇占勝,鄭廣芬,等.寧夏氣候變化的事實分析[J].干旱區(qū)資源與環(huán)境,2005,19(6):43-47.

    Variation of heat condition of maize in the middle arid region of Ningxia

    CHEN Lu1,2,LIU Jing3,WANG Lian-xi1,2LIQi1,2,XU Chen-chun1,2
    (1.Jiangsu Key Laboratory of Agricultural Meteorology,Nanjing University of Information Science&Technology,Nanjing,Jiangsu 210044,China;2.Ningxia Key Laboratory for Meteorological Disaster prevention and Reduction,Yinchuan Ningxia 750002,China)

    Maize growth period data and daily climate data of Yanchi,Tongxin and Haiyuan counties in the middle arid region of Ningxia were used to build the maize temperature suitability degrees.And spatio-temporal analysis of these degrees and active accumulated temperature≥10℃.The result indicated that:the maize temperature suitability degrees and active accumulated temperature≥10℃were in ascent,and the suitability degrees in these three counties from high to low listed as Tongxin county,Yanchi county and Haiyuan county,and the suitability degrees are 0.966,0.911 and 0.726 respectively.The temperature suitability degrees in sowing to emergence,emergence to jointing,jointing to tasseling and silking to maturity periods represented an increase,while the degrees in jointing to tasseling period increased the fast,and the degrees are 0.03·10a-1.Active accumulated temperature≥10℃in sowing to emergence and emergence to jointing periods represented an increase.

    maize growth period;temperature suitability degrees;active accumulated temperature;temporal and spatial variation;themiddle arid region of Ningxia

    S513;S161.2

    A

    1000-7601(2016)06-0257-09

    10.7606/j.issn.1000-7601.2016.06.39

    2016-01-24

    中國氣象局氣候變化專項“氣候變化對寧夏旱作農(nóng)業(yè)的影響與適應對策”(CCSF201315);國家科技支撐計劃(2012BAH29B03)

    陳璐(1991—),女,江蘇常州人,研究生,研究方向為農(nóng)業(yè)氣象。E-mail:844313487@qq.com。

    王連喜(1959—),男,河南省鄭州人,教授。E-mail:wlx4533@sina.com。

    猜你喜歡
    海原積溫寧夏
    大寫的“人”字(外一篇)
    湟源縣30年0℃活動積溫變化特征分析
    踏訪“寰球大震”遺址
    寧夏
    寧夏畫報(2019年5期)2019-09-19 14:42:26
    石河子地區(qū)近50年積溫變化特征分析
    寧夏煤電博弈
    能源(2017年5期)2017-07-06 09:25:50
    淺談“花兒”的保護
    淺談引起海原水土流失的人為因素
    寧夏
    寧夏畫報(2016年5期)2016-06-28 05:48:08
    寧夏:首辦立法論辯會
    浙江人大(2014年4期)2014-03-20 16:20:15
    少妇 在线观看| 久久精品国产亚洲网站| 亚洲一级一片aⅴ在线观看| 亚洲国产精品999| 建设人人有责人人尽责人人享有的| 黑人猛操日本美女一级片| 自拍偷自拍亚洲精品老妇| 国产在线视频一区二区| 最近2019中文字幕mv第一页| 国产亚洲午夜精品一区二区久久| 一级,二级,三级黄色视频| 欧美+日韩+精品| 内地一区二区视频在线| 春色校园在线视频观看| 九色成人免费人妻av| 国产永久视频网站| 色5月婷婷丁香| 九九爱精品视频在线观看| 一本大道久久a久久精品| 国产亚洲最大av| 曰老女人黄片| 亚洲精品国产成人久久av| 国产一区二区在线观看日韩| 国产有黄有色有爽视频| 一级爰片在线观看| 国产精品嫩草影院av在线观看| 欧美精品国产亚洲| 深夜a级毛片| 亚洲精品视频女| 天天操日日干夜夜撸| 久久狼人影院| 亚洲人成网站在线播| 亚洲自偷自拍三级| 大话2 男鬼变身卡| 午夜福利,免费看| 免费人成在线观看视频色| 五月天丁香电影| 亚洲欧洲国产日韩| 五月开心婷婷网| 久久久欧美国产精品| 亚洲人成网站在线播| 在现免费观看毛片| 高清欧美精品videossex| 最黄视频免费看| 欧美精品人与动牲交sv欧美| 精品亚洲成a人片在线观看| 99九九在线精品视频 | 国产 精品1| 免费久久久久久久精品成人欧美视频 | 久久人人爽av亚洲精品天堂| 日韩伦理黄色片| 亚洲第一区二区三区不卡| 久久久欧美国产精品| 晚上一个人看的免费电影| 一边亲一边摸免费视频| av播播在线观看一区| 中文字幕人妻丝袜制服| 岛国毛片在线播放| 久久亚洲国产成人精品v| 免费少妇av软件| 久热久热在线精品观看| 婷婷色综合www| 一区二区三区四区激情视频| 少妇人妻 视频| 亚洲精品,欧美精品| 男女边摸边吃奶| 亚洲第一av免费看| 国产中年淑女户外野战色| 久久6这里有精品| 久久久精品94久久精品| 婷婷色综合大香蕉| 国产极品天堂在线| 日韩人妻高清精品专区| 久久婷婷青草| 精品亚洲乱码少妇综合久久| 简卡轻食公司| 国产成人a∨麻豆精品| 国产亚洲午夜精品一区二区久久| 欧美国产精品一级二级三级 | 国产国拍精品亚洲av在线观看| 哪个播放器可以免费观看大片| 一级毛片黄色毛片免费观看视频| a级片在线免费高清观看视频| 天堂俺去俺来也www色官网| 狂野欧美激情性xxxx在线观看| 亚洲av中文av极速乱| 国产成人精品福利久久| 国产在线视频一区二区| 波野结衣二区三区在线| 少妇被粗大的猛进出69影院 | 日日啪夜夜爽| 久久人人爽人人片av| 久久精品夜色国产| 亚洲自偷自拍三级| 欧美三级亚洲精品| 亚洲精品一二三| 久久国产精品男人的天堂亚洲 | .国产精品久久| 内地一区二区视频在线| 在线观看人妻少妇| 午夜福利视频精品| 在线观看三级黄色| 黄片无遮挡物在线观看| 欧美精品一区二区大全| 国产精品麻豆人妻色哟哟久久| 少妇熟女欧美另类| 欧美 亚洲 国产 日韩一| 久久精品久久精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 五月天丁香电影| 少妇猛男粗大的猛烈进出视频| 噜噜噜噜噜久久久久久91| 日韩中文字幕视频在线看片| 日韩伦理黄色片| 亚洲精品乱码久久久久久按摩| 亚洲成人av在线免费| 国产精品99久久久久久久久| 精品视频人人做人人爽| 国产深夜福利视频在线观看| 性色av一级| 男人爽女人下面视频在线观看| 一区二区三区四区激情视频| 中文精品一卡2卡3卡4更新| 97在线视频观看| 国产一区亚洲一区在线观看| 国产亚洲最大av| 一级毛片电影观看| 日韩三级伦理在线观看| 久久免费观看电影| 夫妻性生交免费视频一级片| 五月开心婷婷网| 成年女人在线观看亚洲视频| 特大巨黑吊av在线直播| 免费在线观看成人毛片| 少妇人妻精品综合一区二区| 乱系列少妇在线播放| 91精品伊人久久大香线蕉| 热re99久久国产66热| 黄色毛片三级朝国网站 | 秋霞在线观看毛片| 2021少妇久久久久久久久久久| 丰满人妻一区二区三区视频av| 欧美日韩视频精品一区| 日本av免费视频播放| 国产黄频视频在线观看| h日本视频在线播放| 在线看a的网站| 久久精品国产亚洲av涩爱| 精品视频人人做人人爽| 中文字幕久久专区| 亚洲av成人精品一区久久| 高清在线视频一区二区三区| 自线自在国产av| 欧美成人午夜免费资源| 亚洲美女搞黄在线观看| 三上悠亚av全集在线观看 | 纵有疾风起免费观看全集完整版| 天堂中文最新版在线下载| 一级毛片黄色毛片免费观看视频| 欧美区成人在线视频| 看十八女毛片水多多多| 成人亚洲精品一区在线观看| 亚洲成色77777| 夫妻性生交免费视频一级片| 黄色日韩在线| 大陆偷拍与自拍| 亚洲成色77777| 久久久久久久久久人人人人人人| 久久久久久久久久久免费av| 国产精品一区二区三区四区免费观看| 国产欧美亚洲国产| 久久久久久久久久久免费av| 国产熟女欧美一区二区| 欧美日韩精品成人综合77777| 国产一区二区三区av在线| 一级片'在线观看视频| 99九九线精品视频在线观看视频| 日日爽夜夜爽网站| 国产一级毛片在线| 涩涩av久久男人的天堂| 中文资源天堂在线| 久久久亚洲精品成人影院| 亚洲精品国产成人久久av| 国内精品宾馆在线| 久久久国产欧美日韩av| 十八禁网站网址无遮挡 | 国产极品粉嫩免费观看在线 | 日韩视频在线欧美| a级片在线免费高清观看视频| 新久久久久国产一级毛片| 欧美日韩亚洲高清精品| 国产av一区二区精品久久| 亚洲精品一二三| 成年人免费黄色播放视频 | 边亲边吃奶的免费视频| 婷婷色麻豆天堂久久| 亚洲av成人精品一二三区| 有码 亚洲区| 久久韩国三级中文字幕| 欧美精品亚洲一区二区| 国产精品伦人一区二区| 美女视频免费永久观看网站| 欧美成人午夜免费资源| 99久久人妻综合| 亚洲欧美日韩另类电影网站| 性色avwww在线观看| 青青草视频在线视频观看| 亚洲av福利一区| 国产男人的电影天堂91| 中国国产av一级| 亚洲av福利一区| 多毛熟女@视频| 99久久综合免费| 曰老女人黄片| 99热这里只有精品一区| 成人美女网站在线观看视频| 久久久a久久爽久久v久久| 赤兔流量卡办理| 五月开心婷婷网| 成人毛片a级毛片在线播放| 不卡视频在线观看欧美| 久久国产精品男人的天堂亚洲 | 蜜臀久久99精品久久宅男| 亚洲精品自拍成人| 国产黄色免费在线视频| 亚洲av国产av综合av卡| 日韩伦理黄色片| 少妇人妻久久综合中文| 亚洲内射少妇av| 欧美日韩国产mv在线观看视频| 亚洲成人手机| 亚洲va在线va天堂va国产| 曰老女人黄片| 国产精品久久久久久精品古装| 九九爱精品视频在线观看| 久久国产精品男人的天堂亚洲 | 中国三级夫妇交换| 亚洲内射少妇av| 日本av手机在线免费观看| 好男人视频免费观看在线| 亚洲精品日本国产第一区| 欧美精品亚洲一区二区| 特大巨黑吊av在线直播| 欧美精品高潮呻吟av久久| 日本黄色日本黄色录像| 国产免费一级a男人的天堂| 国产精品福利在线免费观看| 精品国产乱码久久久久久小说| 桃花免费在线播放| 噜噜噜噜噜久久久久久91| 国产成人a∨麻豆精品| 麻豆成人午夜福利视频| 国产精品久久久久久久久免| 国产乱人偷精品视频| 午夜免费观看性视频| 国产在线视频一区二区| 免费观看av网站的网址| 国产黄色免费在线视频| 日本wwww免费看| 黑人猛操日本美女一级片| 在线观看免费高清a一片| 2021少妇久久久久久久久久久| 特大巨黑吊av在线直播| 国产色婷婷99| 成年人免费黄色播放视频 | 久久久午夜欧美精品| 一级毛片 在线播放| 久久 成人 亚洲| 国产成人aa在线观看| 一个人免费看片子| 亚洲国产毛片av蜜桃av| 国产欧美日韩精品一区二区| 国产成人精品福利久久| 天天操日日干夜夜撸| 欧美bdsm另类| 多毛熟女@视频| 99久久综合免费| av女优亚洲男人天堂| 女的被弄到高潮叫床怎么办| 欧美精品高潮呻吟av久久| 久久久久国产网址| 午夜av观看不卡| 成人国产麻豆网| 成人漫画全彩无遮挡| 中文乱码字字幕精品一区二区三区| 秋霞在线观看毛片| 极品少妇高潮喷水抽搐| 精品视频人人做人人爽| 免费大片黄手机在线观看| 人妻夜夜爽99麻豆av| 自拍偷自拍亚洲精品老妇| 国产男人的电影天堂91| 国产免费一区二区三区四区乱码| 男人爽女人下面视频在线观看| 妹子高潮喷水视频| 99热这里只有是精品50| 欧美日本中文国产一区发布| 日本黄大片高清| 国产免费一级a男人的天堂| 精品人妻偷拍中文字幕| 国产一区二区在线观看av| 欧美性感艳星| 男人和女人高潮做爰伦理| 嘟嘟电影网在线观看| 国产在线男女| 国产精品久久久久久久电影| 黄片无遮挡物在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 久久国产精品男人的天堂亚洲 | 偷拍熟女少妇极品色| 国产片特级美女逼逼视频| 日本91视频免费播放| 免费人成在线观看视频色| 日韩熟女老妇一区二区性免费视频| 91午夜精品亚洲一区二区三区| 久久人人爽人人爽人人片va| 日韩亚洲欧美综合| 一级爰片在线观看| 日韩欧美 国产精品| 国产亚洲欧美精品永久| 久久久久网色| 久久久久久人妻| 黑人巨大精品欧美一区二区蜜桃 | 如日韩欧美国产精品一区二区三区 | 日韩av在线免费看完整版不卡| 欧美bdsm另类| 成人18禁高潮啪啪吃奶动态图 | 久久国产乱子免费精品| 偷拍熟女少妇极品色| 少妇 在线观看| 三级国产精品片| 欧美日韩亚洲高清精品| 国产免费又黄又爽又色| 欧美日韩一区二区视频在线观看视频在线| 国产一区二区三区av在线| 亚洲内射少妇av| 色婷婷久久久亚洲欧美| 黑人巨大精品欧美一区二区蜜桃 | 日本wwww免费看| 免费大片18禁| 免费在线观看成人毛片| 国产欧美日韩精品一区二区| 日本-黄色视频高清免费观看| 精品国产一区二区久久| 在线观看美女被高潮喷水网站| 免费观看的影片在线观看| 亚洲图色成人| 黑人巨大精品欧美一区二区蜜桃 | 国产黄片视频在线免费观看| 国语对白做爰xxxⅹ性视频网站| 熟妇人妻不卡中文字幕| 久久久久国产网址| 菩萨蛮人人尽说江南好唐韦庄| 丰满少妇做爰视频| 99热这里只有是精品50| 国产av码专区亚洲av| 精品国产乱码久久久久久小说| 国产精品免费大片| 一本大道久久a久久精品| 日韩视频在线欧美| 亚洲av欧美aⅴ国产| 久久久国产一区二区| 最近手机中文字幕大全| 亚洲精品国产av蜜桃| 亚洲国产日韩一区二区| 国产真实伦视频高清在线观看| 青春草亚洲视频在线观看| 亚洲精品中文字幕在线视频 | 夫妻性生交免费视频一级片| 人人妻人人澡人人看| 国模一区二区三区四区视频| 日本黄色片子视频| h视频一区二区三区| 色94色欧美一区二区| 日韩中字成人| 青青草视频在线视频观看| 熟女人妻精品中文字幕| 日韩电影二区| 婷婷色麻豆天堂久久| 青青草视频在线视频观看| 99热这里只有是精品在线观看| 午夜视频国产福利| 亚洲一区二区三区欧美精品| 国产亚洲午夜精品一区二区久久| 丝袜脚勾引网站| 人妻 亚洲 视频| 亚洲精品自拍成人| 成人影院久久| 麻豆精品久久久久久蜜桃| 免费看不卡的av| 男女边摸边吃奶| 欧美另类一区| 一级黄片播放器| 人人妻人人爽人人添夜夜欢视频 | 啦啦啦在线观看免费高清www| 高清欧美精品videossex| 免费少妇av软件| 这个男人来自地球电影免费观看 | 少妇的逼水好多| av在线观看视频网站免费| 日韩成人av中文字幕在线观看| 黑人巨大精品欧美一区二区蜜桃 | √禁漫天堂资源中文www| 老司机亚洲免费影院| 久久久精品免费免费高清| 久久久久久伊人网av| 麻豆精品久久久久久蜜桃| 女性生殖器流出的白浆| 久久精品国产自在天天线| 精品人妻熟女av久视频| 99视频精品全部免费 在线| 一级毛片久久久久久久久女| 男女无遮挡免费网站观看| 少妇人妻久久综合中文| 亚洲精品中文字幕在线视频 | 精品一区在线观看国产| 欧美精品高潮呻吟av久久| 26uuu在线亚洲综合色| 简卡轻食公司| 夜夜看夜夜爽夜夜摸| 最近中文字幕高清免费大全6| 成人美女网站在线观看视频| 一区二区三区精品91| 久久国产亚洲av麻豆专区| 91精品一卡2卡3卡4卡| xxx大片免费视频| 精品亚洲乱码少妇综合久久| 亚洲天堂av无毛| 少妇熟女欧美另类| 我的老师免费观看完整版| 大片免费播放器 马上看| 伊人久久精品亚洲午夜| 亚洲精品国产成人久久av| 成人综合一区亚洲| 亚洲av.av天堂| 免费在线观看成人毛片| 日韩欧美 国产精品| 国产精品嫩草影院av在线观看| 国产成人精品福利久久| 免费大片18禁| 久久狼人影院| 五月伊人婷婷丁香| 大香蕉97超碰在线| 在线播放无遮挡| 99久久中文字幕三级久久日本| 在线观看www视频免费| 欧美区成人在线视频| 9色porny在线观看| 成年av动漫网址| 欧美三级亚洲精品| 偷拍熟女少妇极品色| 国产亚洲av片在线观看秒播厂| 亚洲av国产av综合av卡| 国产免费一区二区三区四区乱码| 亚洲欧美成人精品一区二区| 久久久久精品久久久久真实原创| 国产淫片久久久久久久久| 内射极品少妇av片p| 国产男人的电影天堂91| 极品教师在线视频| 老司机影院毛片| 国产成人精品婷婷| 美女视频免费永久观看网站| 国产熟女欧美一区二区| 精品久久国产蜜桃| 伊人久久国产一区二区| 亚洲av成人精品一二三区| 日本av手机在线免费观看| 国产爽快片一区二区三区| 中文欧美无线码| 少妇人妻精品综合一区二区| 亚洲av欧美aⅴ国产| 嫩草影院入口| 欧美97在线视频| 69精品国产乱码久久久| 国内精品宾馆在线| 高清视频免费观看一区二区| 一级毛片 在线播放| 国产视频内射| 一级毛片黄色毛片免费观看视频| 国产亚洲最大av| 在线亚洲精品国产二区图片欧美 | 亚洲在久久综合| 在线观看三级黄色| 97在线人人人人妻| 国产av码专区亚洲av| 丰满迷人的少妇在线观看| 久久久久久人妻| 久久99精品国语久久久| 日韩免费高清中文字幕av| 欧美+日韩+精品| 欧美国产精品一级二级三级 | 看非洲黑人一级黄片| 男人爽女人下面视频在线观看| 亚洲成色77777| 最近中文字幕高清免费大全6| kizo精华| 欧美少妇被猛烈插入视频| 一个人免费看片子| 日韩不卡一区二区三区视频在线| 午夜激情福利司机影院| 在线观看av片永久免费下载| 亚洲精品成人av观看孕妇| 亚洲真实伦在线观看| 一级黄片播放器| 日本黄大片高清| 好男人视频免费观看在线| 欧美激情国产日韩精品一区| 一级片'在线观看视频| 日日啪夜夜爽| av免费观看日本| 精品人妻熟女av久视频| 国产在线免费精品| 国产高清国产精品国产三级| 国产免费又黄又爽又色| 极品人妻少妇av视频| 晚上一个人看的免费电影| 亚洲欧美精品专区久久| 国产中年淑女户外野战色| 精品熟女少妇av免费看| 亚洲,一卡二卡三卡| 国产毛片在线视频| 97在线视频观看| 国产日韩欧美在线精品| 国产伦精品一区二区三区视频9| 亚洲欧美日韩另类电影网站| 亚洲电影在线观看av| 91精品国产国语对白视频| 色5月婷婷丁香| 伦理电影免费视频| 极品少妇高潮喷水抽搐| 国产一区二区在线观看日韩| 亚洲欧美日韩另类电影网站| 永久网站在线| 日韩欧美精品免费久久| 大香蕉久久网| 黄色怎么调成土黄色| 欧美高清成人免费视频www| 国产欧美另类精品又又久久亚洲欧美| 亚洲av不卡在线观看| 国产极品天堂在线| 麻豆成人午夜福利视频| 99视频精品全部免费 在线| 日日爽夜夜爽网站| 亚洲精品一二三| 精品久久久久久久久av| av黄色大香蕉| 久久久久精品久久久久真实原创| 2018国产大陆天天弄谢| 亚洲精品国产av蜜桃| 国产精品国产三级专区第一集| 亚洲av免费高清在线观看| 国产高清不卡午夜福利| 亚洲色图综合在线观看| 亚洲无线观看免费| 成人特级av手机在线观看| 久久国产精品男人的天堂亚洲 | 老司机亚洲免费影院| 成人亚洲欧美一区二区av| 亚洲欧美日韩卡通动漫| 熟女av电影| 国产精品欧美亚洲77777| 久久国产精品男人的天堂亚洲 | 日本黄色片子视频| 午夜日本视频在线| 男人爽女人下面视频在线观看| 一级毛片黄色毛片免费观看视频| 美女福利国产在线| 免费大片黄手机在线观看| 国产综合精华液| 久久毛片免费看一区二区三区| 亚洲精品乱久久久久久| 又爽又黄a免费视频| 草草在线视频免费看| 狠狠精品人妻久久久久久综合| 深夜a级毛片| 秋霞在线观看毛片| 日韩av在线免费看完整版不卡| 国产淫片久久久久久久久| 成年美女黄网站色视频大全免费 | 精品少妇内射三级| 岛国毛片在线播放| 99九九线精品视频在线观看视频| 亚洲精华国产精华液的使用体验| 欧美性感艳星| 高清视频免费观看一区二区| 人人澡人人妻人| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 国产成人freesex在线| 日日摸夜夜添夜夜爱| 国产白丝娇喘喷水9色精品| 亚洲色图综合在线观看| 精品酒店卫生间| 国产成人免费无遮挡视频| 大陆偷拍与自拍| 中文乱码字字幕精品一区二区三区| 制服丝袜香蕉在线| 五月天丁香电影| av网站免费在线观看视频| 亚洲欧美一区二区三区国产| a级毛片在线看网站| 亚洲精品成人av观看孕妇| 汤姆久久久久久久影院中文字幕| 国产精品国产三级国产专区5o| 国产精品偷伦视频观看了| 超碰97精品在线观看| 国产精品无大码| 久久久久久久大尺度免费视频| 69精品国产乱码久久久| 午夜精品国产一区二区电影| 极品人妻少妇av视频| 又黄又爽又刺激的免费视频.| 久久人人爽人人爽人人片va| 青春草视频在线免费观看|