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

    氣候變化對黃淮海地區(qū)夏玉米-冬小麥種植模式和產量的影響

    2017-10-13 20:34:05孫新素龍致煒宋廣鵬陳長青
    中國農業(yè)科學 2017年13期
    關鍵詞:產量影響

    孫新素,龍致煒,宋廣鵬,陳長青

    ?

    氣候變化對黃淮海地區(qū)夏玉米-冬小麥種植模式和產量的影響

    孫新素,龍致煒,宋廣鵬,陳長青

    (南京農業(yè)大學應用生態(tài)研究所,南京210095)

    【目的】全球氣候正以變暖為主要特征發(fā)生顯著變化,探究氣候變化對黃淮海地區(qū)夏玉米-冬小麥種植制度的影響,為制定合理的應對措施提供理論依據?!痉椒ā客ㄟ^氣象站點觀測值的加權平均和一元線性回歸分析黃淮海各省市地區(qū)1992—2013年來的氣候變化特征。利用農業(yè)氣象站點多年長期觀察的夏玉米-冬小麥物候數據,通過加權求平均,分析氣候變暖背景下夏玉米-冬小麥的生育期和茬口推移情況。采用一元線性回歸分析1992—2013年來黃淮海地區(qū)夏玉米-冬小麥周年產量變化。同時利用非線性回歸分析法和面板數據敏感性分析法分析氣候變化對黃淮海地區(qū)夏玉米-冬小麥周年產量的影響。【結果】1992—2013年來,黃淮海地區(qū)溫度整體呈現波動上升趨勢,降水總量變化趨勢不明顯,但區(qū)域差異顯著。在氣候變化的背景下,黃淮海地區(qū)夏玉米-冬小麥種植模式發(fā)生明顯改變:冬小麥播種時間推遲,生育期存在縮短趨勢,不同地區(qū)縮短2—5 d不等;夏玉米播種時間南部推遲而北部提前,收獲時間總體呈現推遲趨勢,整個黃淮海地區(qū)生育時長未發(fā)生明顯變化。茬口時間因夏玉米-冬小麥生育期的推移呈現不同程度延長,造成了氣候和土地資源的浪費。1992—2013年間黃淮海地區(qū)夏玉米-冬小麥單產呈顯著上升趨勢,多數省份達到顯著水平。非線性敏感性分析表明,最低溫度、最高溫度和平均溫度對夏玉米-冬小麥產量的影響基本表現為同時增產或同時減產的一致性。冬小麥產量受最低溫度的影響最為顯著,東南部的江蘇省和山東省減產明顯,而北部河北省和西部河南省表現為增產。溫度升高除對河南省夏玉米有增產作用外,其他省份夏玉米產量均出現不同程度的降低,這可能與溫度升高的幅度不同和降水的區(qū)域性差異有關。降水量對夏玉米-冬小麥產量影響存在地區(qū)差異??傮w上氣候變暖對周年單產影響表現為北部增產,而南部減產,因而選擇適宜早播且生育期長的夏玉米品種對保障周年產量具有重要意義?!窘Y論】氣候變暖背景下,黃淮海地區(qū)冬小麥播種時間推遲,生育期縮短,夏玉米生育期北部延長而南部縮短,生育期的推移導致茬口時間延長,造成了氣候資源和土地資源的浪費。1992—2013年間夏玉米-冬小麥周年產量顯著提高。溫度升高和降水增加對產量的影響存在區(qū)域差異,整個區(qū)域平均來看升溫使夏玉米減產,冬小麥增產;降水增加有利于黃淮海北部地區(qū)夏玉米的產量形成,對南部地區(qū)夏玉米產量則存在不利影響,而對黃淮海大部分地區(qū)冬小麥的產量形成不利。

    氣候變化;夏玉米-冬小麥;種植模式;產量

    0 引言

    【研究意義】全球氣候正以變暖為主要特征發(fā)生顯著變化。IPCC第五次評估報告指出,過去130年(1883—2012年)全球升溫0.85℃,預計21世紀末,全球地表溫度將比21世紀初升高0.3—4.8℃[1]。氣候變化改變了熱量資源的時空格局分布[2],對農業(yè)生產影響顯著,不僅影響到農產品數量,還關系到國家糧食生產和安全供應。黃淮海平原是中國重要的農業(yè)生產基地,夏玉米-冬小麥是該區(qū)的主要復種方式。近10年來,黃淮海地區(qū)冬小麥播種面積占全國冬小麥播種面積的38%,冬小麥產量占全國的50%左右,而夏玉米的播種面積占全國夏玉米播種面積的28%,產量約占全國的30%。隨著氣候變化日益加劇,傳統(tǒng)的夏玉米-冬小麥種植模式受到嚴峻挑戰(zhàn)。研究夏玉米-冬小麥種植模式在氣候變化背景下生育期和茬口推移以及周年產量變化,探究氣候因子與種植模式變化的內在聯系與機制,緩解乃至消除氣候變化對該區(qū)夏玉米-冬小麥生產的消極影響,對保障黃淮海地區(qū)農業(yè)的可持續(xù)發(fā)展具有重要意義?!厩叭搜芯窟M展】全球氣候變化的不斷加劇,使關于氣候變化的研究成為學術界的熱點[3]。前人關于氣候變化對農業(yè)生產影響的研究已相當廣泛,王斌等[4]研究表明氣候變暖使冬小麥全生育期平均縮短2.2 d·(10a)-1,北方冬麥區(qū)生育期縮短最明顯。其他研究也表明華北、淮北平原冬小麥返青期、開花期和成熟期表現為提前的趨勢[5-8]。薛昌穎[9]指出氣候變暖使水稻潛在生育期延長,孟林[10]在對華北平原夏玉米生育期的研究時發(fā)現華北平原夏玉米全生育期天數對氣候變暖呈極顯著增加的響應趨勢。在作物產量方面,熊偉等[11]認為氣候變化將導致中國玉米主產區(qū)單產降低,總產下降。居輝等[12]也對中國不同地區(qū)冬小麥產量對氣候變化的響應進行了研究,結果以減產為主,不同地區(qū)存在明顯的區(qū)域差異?!颈狙芯壳腥朦c】以往的研究多集中在某一作物上[13-14],而對作物種植模式系統(tǒng)的研究卻相對較少,無法準確判斷氣候變化對作物系統(tǒng)的整體影響[15]。而僅僅基于歷史氣象數據的模擬模型分析,很難全面揭示氣候變化與作物系統(tǒng)的實際關系[16]。因此,本文利用近年來黃淮海地區(qū)氣象、物候和產量數據,采用回歸分析和敏感性分析,研究氣候變化對夏玉米-冬小麥兩熟制種植模式系統(tǒng)和周年產量的影響?!緮M解決的關鍵問題】利用近年來黃淮海地區(qū)的氣象數據以及夏玉米-冬小麥物候和產量數據,探究氣候變化對夏玉米-冬小麥作物系統(tǒng)的影響,以期在氣候變化背景下為黃淮海夏玉米-冬小麥種植模式的管理調整、合理布局提供基礎性的理論依據。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于東經110°21′—122°42′、北緯29°24′—42°36′,覆蓋中國黃淮海平原的絕大部分地區(qū),主要包括安徽省和江蘇省的北部、河南省、河北省、山東省以及京津6個省市地區(qū),區(qū)域總面積達到5×107hm2,其中農田面積約為2.6×106hm2,是中國農業(yè)生產核心區(qū)之一(圖1)。區(qū)域以暖溫帶季風性氣候為主,表現出明顯的雨熱同季特點,春冬干燥寒冷,旱情較重,夏季高溫多雨,常伴有洪澇的發(fā)生?!?℃積溫為4 100—5 400℃,≥10℃積溫為3 700—4 700℃,年降水量為500—900 mm,降水年際變化較大。年平均溫度和年降水量均隨著緯度的增加表現為減少的趨勢。

    1.2 資料來源

    數據資料主要包括氣象、物候和產量數據3個方面。

    圖1 研究區(qū)域在中國的位置

    1.2.1 氣象數據 氣象數據主要包括逐日溫度(最高溫度、平均溫度和最低溫度)和逐日降水量,來源于中國氣象局國家氣象信息中心。本研究共選擇有完整數據記錄的82個黃淮海地區(qū)氣象觀測站,氣象站在研究區(qū)的分布情況為安徽省北部6個,江蘇省北部13個,河北省20個,河南省18個,山東省19個,京津地區(qū)6個(圖2)。所有氣象數據用于計算夏玉米生育進程(夏玉米播種期-收獲期)和冬小麥生育進程(冬小麥播種期-收獲期)不同年份各氣象指標值。

    圖2 氣象站點分布圖

    1.2.2 物候數據 物候數據同樣整理自國家氣象信息中心,本文選擇了研究區(qū)的91個農業(yè)氣象站,記錄了1992—2013年夏玉米和冬小麥的具體生育期數值。

    1.2.3 產量數據 產量數據主要包括1992—2013年黃淮海各?。ㄊ?,地區(qū))的夏玉米和冬小麥播種面積、單產和總產量,整理摘錄自中國農業(yè)信息中心。

    本文對所有的數據進行初步的篩選和質量控制,對極個別缺失的數據,例如某一氣象觀測站點個別日期的氣象數據,采用世界氣象組織規(guī)定的氣候標準時段(1980—2010年)的平均值進行填充,并與該缺失值前后5年觀測值的平均值進行比對。對于物候數據某一站點某一生育時期數據的缺失,采用該區(qū)域內其他站點的平均值進行估算,對該缺失數據進行替換,剔除缺失數據的影響[17]。缺失數據量所占比例極小,氣象數據不超過數據總量的0.1‰,物候觀測數據約占數據量的1%,產量數據不存在缺失。因津京地區(qū)農業(yè)生產受城市化進程沖擊較大,記錄數據不能準確反映農業(yè)生產實際情況,故不對其進行過多研究。

    1.3 分析方法

    本研究主要運用了數理統(tǒng)計的相關分析方法。

    1.3.1 氣候變化 氣象數據處理是對82個氣象觀測站點,根據其附近農業(yè)氣象站點物候觀測多年的平均值(播種期和成熟期),分別計算出該站點夏玉米和冬小麥生長季內各年的氣象指標平均值;根據氣象站點所在不同省份分類,計算出各省夏玉米和冬小麥生長季內各氣象指標的年際平均數值,再根據各省的面積比加權計算出黃淮海區(qū)域夏玉米和冬小麥生長季的年際氣象數值,通過簡單一元線性模型分別對1992—2013年夏玉米生長季和冬小麥生長季溫度和降水量變化趨勢進行擬合分析。

    1.3.2 周年單產變化 產量數據依據記錄的夏玉米和冬小麥年單產進行年際間的一元線性回歸分析,分析夏玉米-冬小麥種植模式的周年產量變化趨勢[18]。

    1.3.3 生育進程 各省年際間夏玉米和冬小麥物候數據平均值是通過計算不同年份研究區(qū)域內不同省份各觀測站點的不同生育期的平均值獲得,而整個研究區(qū)域物候數據是按各省面積比加權計算出的平均值。氣候變化是一個長期而緩慢的過程,為了突顯這種變化,將1992—2013年這22年分為1992—1996年和2009—2013年首尾兩個5年時間段進行比較分析,表征1992—2013年的氣候變化情況。通過處理得到的黃淮海地區(qū)不同省份的年際氣象數據和物候數據,分別計算1992—1996年和2009—2013年首尾兩個5年時間段的5年平均值,同樣將黃淮海地區(qū)氣象數據和物候數據分不同省進行同樣處理,以分析研究區(qū)內不同地區(qū)氣候變化對夏玉米-冬小麥的生育期及茬口推移影響情況。

    1.3.4 氣候變化對周年產量的影響 作物產量變化除受氣候變化影響外,政策導向和科技的進步等非氣候因素的影響同樣不容忽視。為了更準確的分析氣候變化和作物產量之間的具體關系,必須剔除影響因素中的非氣候因子[19]。研究認為年份的線性趨勢和二次趨勢可以看成剔除非氣候因素對產量的貢獻[14,20],Lobell的多元非線性回歸模型[21]考慮到了溫度、降水、技術進步、政策等因子,以年份的二次回歸,分析出技術進步和管理措施改進帶來的產量增加趨勢。該模型是非線性的,不同于線性模型只利用偏平均值的方式,而是同時利用年際間的歷史氣象變量值(溫度和降水的平均值)來鑒定模型系數。因此在農業(yè)氣候變化的背景下,這被看作是分析氣候變化對作物產量影響的有效舉措[22-23]。本研究首先利用Lobell的多元非線性回歸模型對夏玉米-冬小麥產量進行去趨勢處理,再進一步對面板數據進行分析,確定氣候變化對周年產量的影響。即:

    Log(Y)=c+d×+d×2+×X+

    式中,為研究區(qū)的不同地區(qū);=1992,…,2013,為年份時間變量;為固定效應,為時間趨勢系數,為自變量參數,X為自變量,包括平均氣溫、最高氣溫、最低氣溫和降水等氣候因子,Y為作物產量,ε為誤差項。

    2 結果

    2.1 夏玉米-冬小麥生長季內氣候變化

    1992—2013年黃淮海平原夏玉米-冬小麥生長季內氣溫整體呈波動上升趨勢。從圖3-A、B可以看出,除冬小麥生長季最高溫度出現小幅下降趨勢外,其他溫度因子均表現為不同程度的增加。夏玉米生長季內,最低氣溫、最高氣溫和平均氣溫分別上漲0.44、0.10、0.20℃·(10a)-1,其中最低氣溫的上升幅度達到顯著水平。冬小麥生長季最低氣溫和平均氣溫分別上漲0.46和0.11℃·(10a)-1,最低氣溫上升幅度達顯著水平。溫度升高趨勢表現為最低溫度上升最明顯,其次為平均溫度,而最高溫度并未出現明顯變化。

    降水量表現為小幅度的上升趨勢(圖3-C)。玉米生長季降水量以27.87 mm·(10a)-1的速度增加,小麥生長季降水量增幅較小,為4.61 mm·(10a)-1。整個夏玉米-冬小麥內,降水量以32.48 mm·(10a)-1的速度緩慢增加,但年際間的波動性很大。

    以1992—1996年和2009—2013年首尾兩個時間段表征的黃淮海1992—2013年間不同?。ㄊ?、區(qū)域)的氣候因子變化情況(表1)。從溫度因子來看,各省、市、區(qū)域除冬小麥生長季最高溫度出現小幅下降外,其他溫度因子均表現為不同程度的增加,溫度上升0.1—0.9℃不等,尤其以最低溫度的上升最為明顯,其中山東省夏玉米生長季最低溫度上升達0.94℃,溫度升高較為明顯。整個研究區(qū)來看,夏玉米生長季最低溫度、最高溫度和平均溫度均表現為上升的趨勢,分別增加了0.86℃、0.52℃、0.61℃,冬小麥生長季最低溫度和平均溫度分別增加了0.55℃、0.15℃,最高溫度則降低了0.22℃。整個黃淮海地區(qū)降水總量并未發(fā)生顯著的變化,但不同區(qū)域降水量的空間差異較為明顯,年際間降水量的不穩(wěn)定性較為突出。

    *表示顯著(P<0.05) *represents significant at P<0.05

    2.2 夏玉米-冬小麥種植模式生育時間及茬口推移變化

    1992—2013年來,黃淮海地區(qū)近5年(2009—2013年)與初始5年(1992—1996年)相比,夏玉米-冬小麥種植模式生育時間推移情況如圖4所示。氣候變暖的背景下,黃淮海各地區(qū)冬小麥播種時間均呈現不同程度的推遲,其中山東省推遲達6 d,而冬小麥收獲期因不同地區(qū)播種時間推遲的程度不同體現出提前或推遲趨勢,但就冬小麥整個生育進程來看,黃淮海各地區(qū)均表現為生育期縮短的趨勢,生育時長平均縮短3 d左右。夏玉米生育期推移情況呈現區(qū)域差異,東南部的江蘇及安徽北部和山東省夏玉米播種時間推遲,而位于黃淮海西部和北部的河南河北則呈現小幅度的提前。夏玉米收獲時間除河南省以外,均呈現推遲趨勢。整個地區(qū)總生育時長未發(fā)生顯著變化,但區(qū)域間的差異較為明顯,這可能與降水的區(qū)域差異有關。夏玉米以及冬小麥生育時間的變化,顯著改變了夏玉米-冬小麥種植模式的茬口時間,冬小麥生育時長的顯著縮短,使茬口時間出現不同程度的延長,導致了土地和熱量資源的浪費。

    表1 黃淮海不同地區(qū)夏玉米-冬小麥生長季氣象因子變化

    平均值±標準誤average±standard error

    2.3 夏玉米-冬小麥周年單產變化

    1992—2013年來,黃淮海地區(qū)夏玉米-冬小麥種植模式單產除江蘇和安徽北部的夏玉米外,均呈現年際間顯著增加的趨勢(表2)。無論夏玉米還是冬小麥,都表現為較大幅度的增產,其中冬小麥較夏玉米增產幅度更大,但增產表現不如夏玉米穩(wěn)定。不同省份地區(qū)來看,除江蘇和安徽夏玉米外,其他地區(qū)無論夏玉米還是冬小麥,增產均達到極顯著水平。其中夏玉米單位產量增加速度以山東省最為明顯,達到77.94 kg·hm-2·a-1,比黃淮海地區(qū)平均值高出25.39 kg·hm-2·a-1;冬小麥單產和夏玉米-冬小麥周年單產增速均以河南省最高,增幅分別為122.25和104.36 kg·hm-2·a-1,分別比黃淮海地區(qū)平均值高出36.28和29.01 kg·hm-2·a-1,表現出很明顯的增產趨勢,但變異度也最高。

    2.4 氣候變暖對夏玉米-冬小麥種植模式周年單產的影響

    農作物產量不僅受到氣候變化調節(jié),品種、科學技術進步等也對其有不可忽視的影響。依據面板數據分析,對黃淮海各省(市、地區(qū))的周年產量對溫度和降水的敏感性進行分析。溫度每升高1℃、降水每增加100 mm,夏玉米-冬小麥周年單產變化結果如圖5所示。

    WS:小麥播種期Wheat sowing;WH:小麥收獲期Wheat harvest;MS:玉米播種期Maize sowing;MH:玉米收獲期Maize harvest

    表2 1992—2013年黃淮海地區(qū)夏玉米-冬小麥單產變化

    **表示在0.01水平上差異顯著;NS 表示無顯著性差異 **, significant difference at the 0.01 level; NS, no significance

    溫度升高對整個黃淮海地區(qū)夏玉米單產的影響并不明顯,但不同地區(qū)間的差異較為顯著。河南省對溫度升高的響應表現為正,最低溫度升高1℃產量上升達7.27%,而江蘇及安徽北部、山東和河北均表現為減產的負效應,尤以安徽省北部最為明顯,最低溫度升高1℃減產8.04%。冬小麥產量對溫度升高的響應以最低溫度最為明顯。不同地區(qū)來看,河南、河北和安徽北部表現為正效應,其中安徽北部增產幅度最大,最低溫度上升1℃增產3.95%,山東和江蘇北部表現為負效應,減產幅度基本相同。整個黃淮海地區(qū)小麥對溫度的整體響應表現為正效應,溫度上升1℃產量增加1%左右,3個溫度因子的影響程度基本相同。溫度升高對黃淮海地區(qū)夏玉米-冬小麥周年單產有小幅度的促進作用,這主要是由于河南和安徽較大程度的增產效應抵消了山東和江蘇的負面影響。

    圖5 黃淮海地區(qū)夏玉米-冬小麥種植模式周年單產對溫度(+1℃)和降水(+100 mm)敏感性的面板數據分析

    降水增加對夏玉米產量影響的區(qū)域性差異較大,對夏季較為干旱的山東和河北,降水的增產效果明顯,每增加100 mm降水,河北省夏玉米單產增加6.47%;而對處于南部的河南省以及安徽和江蘇北部,降水增加則會導致不同程度的減產。整個黃淮海地區(qū)夏玉米對降水的響應表現為正效應,每增加100 mm降水,產量增加1.96%。而冬小麥單產對降水增加的反饋在黃淮海不同地區(qū)表現基本相同,除河北表現為小幅增產外,其他地區(qū)均表現為減產,其中安徽及江蘇北部最為明顯,減產幅度分別達到8.03%和4.73%。整個黃淮海地區(qū)夏玉米-冬小麥種植模式周年單產對降水的響應整體表現為每增加100 mm降水減產0.52%,且不同地區(qū)的差異較大,位于黃淮海南部的河南、江蘇和安徽影響為負,北部的河北和山東影響為正,這和南北降水差異有關。

    3 討論

    以氣候變暖為主要特征的全球氣候變化已得到各界學者的公認。黃淮海平原作為中國氣候變化敏感區(qū)之一,農業(yè)氣候資源在氣候變化背景下發(fā)生顯著變化[24]。本研究發(fā)現,1992—2013年黃淮海平原冬小麥和夏玉米生長季最高溫度和平均溫度變化不明顯,但最低溫度呈現顯著上升趨勢。不同地區(qū)溫度上升幅度存在區(qū)域差異,以山東省溫度升高最為明顯,西部和北部的河南和河北溫度升高幅度相對較低,但均表現為最低氣溫上升幅度顯著高于最高溫度和平均溫度。這與全球氣候溫度上升的區(qū)域差異性和最低溫度較最高溫度增加明顯的特征相一致。

    氣候變暖對作物生育進程的影響因作物類型的不同而存在顯著差異[25]。以往的研究中,冬小麥對氣候變化的響應表現為生育期的明顯縮短[4,26],這與本研究的結果相一致,不同地區(qū)冬小麥全生育期縮短2—5 d不等,這與不同地區(qū)背景溫度以及溫度增加幅度的不同有關。而夏玉米生育期在本研究中整體推移變化并不明顯,且存在明顯的區(qū)域差異。這一方面是由于黃淮海地區(qū)夏季溫度基數大,小幅度的溫度增加效果不明顯,且夏玉米生長的積溫需求比冬小麥高得多,溫度升高對夏玉米生育期的影響相對冬小麥較?。涣硪环矫嬗衩鬃鳛槎倘照兆魑飳庵芷谳^為敏感[27],日長效應也需要考慮,同時夏季降水在時間和空間上的差異也增加了黃淮海地區(qū)夏玉米生育期推移研究的復雜性。氣候變化對冬小麥和夏玉米生育期的影響,尤其是冬小麥全生育期的縮短,使夏玉米-冬小麥輪作系統(tǒng)茬口時間出現一定時間的延長,造成了氣候資源和土地資源的浪費,而溫度升高有利于喜溫和晚熟品種的種植[28],因此選擇適合早播、生育期長的夏玉米品種,調整夏玉米-冬小麥的播期和種植結構,對黃淮海地區(qū)夏玉米-冬小麥生產具有重要意義。本研究只研究了夏玉米-冬小麥全生育期的變化推移情況,對于夏玉米-冬小麥的營養(yǎng)生長和生殖生長的具體生育階段沒有進行過多研究,而以往的研究表明無論是夏玉米還是冬小麥,營養(yǎng)生長和生殖生長不同階段對氣候變化的響應存在顯著差異[29-30],因此之后的研究方向可集中在黃淮海地區(qū)夏玉米-冬小麥系統(tǒng)不同生育階段上。

    氣候變暖對黃淮海平原夏玉米-冬小麥種植模式產量的影響既是機遇又是挑戰(zhàn)。溫度升高使夏玉米-冬小麥生長季積溫升高,有利于作物生產[31],但溫度升高提前滿足了作物積溫需求,導致物候提前和生育期縮短[32]。因此在這種綜合影響下,增溫對冬小麥產量的影響目前尚無定論[33]。在本研究中,除安徽北部和河南省小麥產量對溫度升高表現為正效應外,山東、河北、江蘇北部均表現為減產的負效應。從地理位置來看,安徽北部和河南省較其他3個地區(qū)更加靠近內陸,地理位置的差異導致了小麥產量對氣候變暖響應的差異。許多研究表明氣候變暖對夏玉米產量的影響表現為下降[34-35],本研究中除河南省外,其他地區(qū)玉米產量對溫度升高的響應均表現為負效應,這與以往的研究結果相類似。而河南省夏玉米季溫度升高的幅度在所有黃淮海地區(qū)中是最低的,因此較低水平的增溫是否對夏玉米產量具有一定的促進作用有待探究,尤其在水分適宜的情況下,有研究表明溫度升高對玉米產量有促進作用[36],暖濕氣候有利于玉米的生產[37]。氣候變暖使降水的區(qū)域差異更加明顯,而水分條件是夏玉米-冬小麥生長發(fā)育和產量形成的基礎[38],降水在時間和空間上的改變對研究區(qū)夏玉米-冬小麥種植模式是一個嚴峻考驗。本研究中夏玉米-冬小麥種植模式周年單產對降水的敏感性分析也表明了氣候變暖帶來的降水區(qū)域性影響差異,降水增加對黃淮海北部較為干旱的河北和山東有促進作用,對較為濕潤的南部地區(qū)則相反。此外,有研究表明,氣候變暖使中國多熟制的北界發(fā)生不同程度的北移,雨養(yǎng)區(qū)冬小麥-夏玉米作物系統(tǒng)穩(wěn)產北界因該區(qū)降水的減少向東南方向移動[39-41]。而本研究結果表明黃淮海地區(qū)在氣候變暖的背景下降水量呈現小幅的增加趨勢,這對黃淮海地區(qū)夏玉米-冬小麥種植模式的生產十分有利。

    氣候變暖對作物系統(tǒng)的影響不僅表現為平均態(tài)影響,非常態(tài)影響更為突出,如氣候波動加劇,極端氣候天氣增加,災害性事件頻率和強度增加等[42]。本研究僅就溫度和降水的平均變化作了探討,在一定程度上可定量反映氣候變化對夏玉米-冬小麥種植模式的影響,但對于氣候因子在時間分布、以及極端氣候事件等方面并沒有涉及,還存在一定的不確定性。從全球變化和區(qū)域糧食生產的角度來看,作物系統(tǒng)的適應性調整是一個長期而緩慢的過程,只有深刻認識氣候變化特征及作物的響應機理,進而調整作物栽培和管理模式,才能不斷提高農業(yè)生產適應氣候變暖和應對氣候災害的能力,保障區(qū)域和國家糧食安全。

    4 結論

    近年來氣候變化對黃淮海地區(qū)夏玉米-冬小麥種植模式產生了較大影響。氣候變暖使黃淮海地區(qū)夏玉米-冬小麥生長季內溫度升高,尤以最低溫度上升最為明顯,同時改變了降水在時間和空間上的分布。氣候變暖使冬小麥播種時間推遲,生育期縮短,夏玉米生育時長無明顯變化,但區(qū)域差異顯著,生育進程的改變延長了夏玉米-冬小麥種植模式的茬口時間,造成了氣候資源和土地資源的浪費。1992—2013年黃淮海平原夏玉米-冬小麥周年單產顯著提高,溫度升高和降水增加對產量的影響存在區(qū)域差異,但整體來看升溫使夏玉米減產,冬小麥增產;降水增加對黃淮海北部夏玉米產量有促進作用,對南部則有抑制作用,對研究區(qū)內冬小麥基本表現為產量的抑制。選擇適合早播且生育期長的夏玉米品種可有效保障夏玉米-冬小麥種植模式周年產量。

    References

    [1] IPCC. Climate change 2013: The physical science basis//. Cambridge: Cambridge University Press, 2013.

    [2] 繆啟龍, 丁園圓, 王勇, 段春鋒. 氣候變暖對中國熱量資源分布的影響分析. 自然資源學報, 2009(5): 934-944.

    MIAO Q L, DING Y Y, WANG Y, DUAN C F. Impact of climate warming on the distribution of China's thermal resources., 2009(5): 934-944. (in Chinese)

    [3] Alexandrov V A, Hoogenboom G. The impact of climate variability and change on crop yield in Bulgaria., 2000, 104(4): 315-327.

    [4] 王斌, 顧蘊倩, 劉雪, 羅衛(wèi)紅, 戴劍鋒, 張巍, 亓春杰. 中國冬小麥種植區(qū)光熱資源及其配比的時空演變特征分析. 中國農業(yè)科學, 2012, 45(2): 228-238.

    WANG B, GU Y Q, LIU X, LUO W H, DAI J F, ZHANG W, QI C J. Analysis of the temporal and spatial changes of photo-thermal resources in winter wheat growing regions in China., 2012, 45(2): 228-238. (in Chinese)

    [5] 潘根興, 高民, 胡國華, 魏欽平, 楊曉光, 張文忠, 周廣勝, 鄒建文. 氣候變化對中國農業(yè)生產的影響. 農業(yè)環(huán)境科學學報, 2011, 30(9): 1698-1706.

    PAN G X, GAO M, HU G H, WEI Q P, YANG X G, ZHANG W Z, ZHOU G S, ZOU J W. Impacts of climate change on agricultural production of China., 2011, 30(9): 1698-1706. (in Chinese)

    [6] 李正國, 楊鵬, 唐華俊, 吳文斌, 陳仲新, 周清波, 鄒金秋, 張莉. 氣候變化背景下東北三省主要作物典型物候期變化趨勢分析. 中國農業(yè)科學, 2011, 44(20): 4180-4189.

    LI Z G, YANG P, TANG H J, WU W B, CHEN Z X, ZHOU Q B, ZOU J Q, ZHANG L. Trend analysis of typical phenophases of major crops under climate change in the three provinces of Northeast China., 2011, 44(20): 4180-4189. (in Chinese)

    [7] 姜麗霞, 李帥, 李秀芬, 張立群, 杜春英. 黑龍江省近三十年氣候變化對大豆發(fā)育和產量的影響. 大豆科學, 2011, 30(6): 921-926.

    JIANG L X, LI S, LI X F, ZHANG L Q, DU C Y. Impacts of climate change on development and yield of soybean over past 30 years in Heilongjiang province.2011, 30(6): 921-926. (in Chinese)

    [8] 肖登攀, 陶福祿. 過去30年氣候變化對華北平原冬小麥物候的影響研究. 中國生態(tài)農業(yè)學報, 2012, 20(11): 1539-1545.

    XIAO D P, TAO F L. Impact of climate change in 1981-2009 on winter wheat phenology in the North China Plain., 2012, 20(11): 1539-1545. (in Chinese)

    [9] 薛昌穎, 劉榮花, 吳騫. 氣候變暖對信陽地區(qū)水稻生育期的影響. 中國農業(yè)氣象, 2010, 31(3): 353-357.

    XUE C Y, LIU R H, WU Q. Effect of climate warming on rice growing stages in Xinyang., 2010, 31(3): 353-357. (in Chinese)

    [10] 孟林, 劉新建, 鄔定榮, 王春乙. 華北平原夏玉米主要生育期對氣候變化的響應. 中國農業(yè)氣象, 2015, 36(4): 375-382.

    MENG L, LIU X J, WU D R, WANG C Y. Responses of summer maize main phenology to climate change in the North China Plain., 2015, 36(4): 375-382. (in Chinese)

    [11] 熊偉, 楊婕, 林而達, 許吟隆. 未來不同氣候變化情景下我國玉米產量的初步預測. 地球科學進展, 2008, 23(10): 1092-1101.

    XIONG W, YANG J, LIN E D, XU Y L. The projection of maize yield in China under climate change scenarios., 2008, 23(10): 1092-1101. (in Chinese)

    [12] 居輝, 熊偉, 許吟隆, 林而達. 氣候變化對我國小麥產量的影響. 作物學報, 2005, 31(10): 1340-1343.

    JU H, XIONG W, XU Y L, LIN E D. Impacts of climate change on wheat yield in China., 2005, 31(10): 1340-1343. (in Chinese)

    [13] Fulco L, Senthold A. Climate change impacts on wheat production in a Mediterranean environment in Western Australia., 2006, 90(1/3): 159-179.

    [14] Urban D, Roberts M J, Schlenker W, Lobell D B. Projected temperature changes indicate significant increase in interannual variability of U.S. maize yields., 2012, 112(2): 525-533.

    [15] 齊永青, 孫宏勇, 沈彥俊. 太行山山前平原近50年氣候變暖特征及其對冬小麥-夏玉米作物系統(tǒng)的影響. 中國生態(tài)農業(yè)學報, 2011, 19(5): 1048-1053.

    Qi Y Q, SUN H Y, SHEN Y J. Characteristics and effects of climate warming on winter wheat/summer maize cropping system in recent 50 years in the piedmont of Mount Taihang., 2011, 19(5): 1048-1053. (in Chinese)

    [16] Lobell D B, Asner G P. Climate and management contributions to recent trends in U.S. agricultural yield., 2009, 299(5609): 1032.

    [17] 龐新生. 缺失數據處理方法的比較. 統(tǒng)計與決策, 2010(24): 152-155.

    PANG X S. The comparison of different missing data processing method., 2010(24): 152-155. (in Chinese)

    [18] Katz R W. Assessing the impact of climatic change on food production.1977, 1(1): 85-96.

    [19] Hu Q, Buyanovsky G. Climate effects on corn yield in Missouri., 2003, 42(11): 1626-1635.

    [20] Schlenker W, Lobell D B. Robust negative impacts of climate change on African agriculture., 2010, 5(1): 1-8.

    [21] Lobell D B, Schlenker W, Costa-Roberts J. Climate trends and global crop production since 1980., 2011, 333(6042): 616-620.

    [22] Shi W, Tao F, Zhang Z. A review on statistical models for identifying climate contributions to crop yields., 2013, 23(3): 567-576.

    [23] Liu B, Asseng S, Müller C, Ewert F, Elliott J, Lobell D B, Martre P, Ruane A C, Wallach D, Jones J W, Rosenzweig C, Aggarwal P K, Alderman P D, Anothai J, Basso B, Biernath C, Cammarano D, Challinor A, Deryng D, Sanctis G D, Doltra J, Fereres E, Folberth C, Garcia-Vila M, Gayler S, Hoogenboom G, Hunt L A, Izaurralde R C, Jabloun M, Jones C D, Kersebaum K C, Kimball B A, Koehler A, Kumar S N, Nendel C, O’Leary G, Olesen J E, Ottman M J, Palosuo T, Prasad P V, Priesack E, Pugh T A M, Reynolds M, Rezaei E E, R?tter R P, Schmid E, Semenov M A, Shcherbak I, Stehfest E, St?ckle C O, Stratonovitch P, Streck T, Supit I, Tao F L, Thorburn P, Waha K, Wall G W, Wang E L, White J W, Wolf J, Zhao Z G, Zhu Y. Similar estimates of temperature impacts on global wheat yield by three independent methods.2016, 6(12): 1130-1136.

    [24] 馬潔華, 劉園, 楊曉光, 王文峰, 薛昌穎, 張曉煜. 全球氣候變化背景下華北平原氣候資源變化趨勢. 生態(tài)學報, 2010, 30(14): 3818-3827.

    MA J H, LIU Y, YANG X H, WANG W F, XUE C Y, ZHANG X Y. Characteristics of climate resources under global climate change in the North China Plain., 2010, 30(14): 3818-3827. (in Chinese)

    [25] 侯瓊, 郭瑞清, 楊麗桃. 內蒙古氣候變化及其對主要農作物的影響. 中國農業(yè)氣象, 2009, 30(4): 560-564.

    HOU Q, GUO R Q, YANG L T. Climate change and its impact on main crops in Inner Mongolia., 2009, 30(4): 560-564. (in Chinese)

    [26] 王力, 李鳳霞, 徐維新, 肖宏斌, 蘇文將, 周萬福. 氣候變化對柴達木灌區(qū)小麥生育期的影響. 中國農業(yè)氣象, 2010, 31(s1): 81-83.

    WANG L, LI F X, XU W X, XIAO H B, SU W J, ZHOU W F. Impact of climatic change on growth of wheat in the Irrigation District of Qaidam Oasis., 2010, 31(S1): 81-83. (in Chinese)

    [27] 鄭國清, 高亮之. 玉米發(fā)育期動態(tài)模擬模型. 江蘇農業(yè)學報, 2000, 16(1): 15-21.

    ZHENG G Q, GAO L Z. Simulation model of maize phenology., 2000, 16(1): 15-21. (in Chinese)

    [28] 房世波, 譚凱炎, 任三學, 張新時. 氣候變暖對冬小麥生長和產量影響的大田實驗研究. 中國科學(地球科學), 2012, 42(7): 1069-1075.

    Fang S B, Tan K Y, Ren S X, Zhang X S. Fields experiments in North China show no decrease in winter wheat yields with night temperature increased by 2.0-2.5℃.(), 2012, 42(7): 1069-1075. (in Chinese)

    [29] 李德. 近30年淮北平原冬小麥物候期演變特征. 氣象科技, 2009, 37(5): 607-612.

    LI D. Variation characteristics of winter wheat phenophases in Huaibei Plain in recent 30 years., 2009, 37(5): 607-612. (in Chinese)

    [30] 陸偉婷, 于歡, 曹勝男, 陳長青. 近20年黃淮海地區(qū)氣候變暖對夏玉米生育進程及產量的影響. 中國農業(yè)科學, 2015, 48(16): 3132-3145.

    LU W T, YU H, CAO S N, CHEN C Q. Effects of climate warming on growth process and yield of summer maize in Huang-Huai-Hai Plain in last 20 years., 2015, 48(16): 3132-3145. (in Chinese)

    [31] 王占彪, 王猛, 尹小剛, 張海林, 褚慶全, 文新亞, 陳阜. 氣候變化背景下華北平原夏玉米各生育期水熱時空變化特征. 中國生態(tài)農業(yè)學報, 2015, 23(4): 473-481.

    WANG Z B, WANG M, YIN X G, ZHANG H L, CHU Q Q, WEN X Y, CHEN F. Spatiotemporal characteristics of heat and rainfall changes in summer maize season under climate change in the North China Plain., 2015, 23(4): 473-481. (in Chinese)

    [32] Fang S B, Cammarano D, Zhou G S, TAN K Y, REN S X. Effects of increased day and night temperature with supplemental infrared heating on winter wheat growth in North China., 2015, 64: 67-77.

    [33] 胡實, 莫興國, 林忠輝. 氣候變化對黃淮海平原冬小麥產量和耗水的影響及品種適應性評估. 應用生態(tài)學報, 2015, 26(4): 1153-1161.

    HU S, MO X G, LIN Z H. Evaluating the response of yield and evapotranspiration of winter wheat and the adaptation by adjusting crop variety to climate change in Huang-Huai-Hai Plain.2015, 26(4): 1153-1161. (in Chinese)

    [34] 馬樹慶, 王琪, 羅新蘭. 基于分期播種的氣候變化對東北地區(qū)玉米()生長發(fā)育和產量的影響. 生態(tài)學報, 2008, 28(5): 2131-2139.

    MA S Q, WANG Q, LUO X L. Effect of climate change on maize()growth and yield based on stage sowing.2008, 28(5): 2131-2139. (in Chinese)

    [35] 張吉旺, 董樹亭, 王空軍, 胡昌浩, 劉鵬. 大田增溫對夏玉米產量和品質的影響. 應用生態(tài)學報, 2007, 18(1): 52-56.

    ZHANG J W, DONG S T, WANG K J, HU C H, LIU P. Effects of high field temperature on summer maize grain yield and quality.2007, 18(1): 52-56. (in Chinese)

    [36] 王琪, 馬樹慶, 郭建平, 張鐵林, 于海, 徐麗萍. 溫度對玉米生長和產量的影響. 生態(tài)學雜志, 2009, 28(2): 255-260.

    WANG Q, MA S Q, GUO J P, ZHANG T L, YU H, XU L P. Effects of air temperature on maize growth and its yield.2009, 28(2): 255-260. (in Chinese)

    [37] 朱新玉, 張竟竟, 趙文亮, 史本林. 河南省商丘市氣候變化對玉米產量的影響. 湖北農業(yè)科學, 2012, 51(11): 2198-2200.

    ZHU X Y, ZHANG J J, ZHAO W L, SHI B L. Impacts of climate change on maize yield in Shangqiu, Henan, China.2012, 51(11): 2198-2200. (in Chinese)

    [38] 肖登攀, 陶福祿, 沈彥俊, 劉劍鋒, 王仁德. 華北平原冬小麥對過去30年氣候變化響應的敏感性研究. 中國生態(tài)農業(yè)學報, 2014, 22(4): 430-438.

    XIAO D P, TAO F L, SHEN Y J, LIU J F, WANG R D. Sensitivity of response of winter wheat to climate change in the North China Plain in the last three decades., 2014, 22(4): 430-438. (in Chinese)

    [39] 郝志新, 鄭景云, 陶向新. 氣候增暖背景下的冬小麥種植北界研究——以遼寧省為例. 地理科學進展, 2001, 20(3): 253-260.

    HAO Z X, ZHENG J Y, TAO X X. A Study on northern boundary of winter wheat during climate warming: A case study in Liaoning province.2001, 20(3): 253-260. (in Chinese)

    [40] 高志強, 苗果園, 鄧志鋒. 全球氣候變化與冬麥北移研究. 中國農業(yè)科技導報, 2004, 6(1): 9-13.

    GAO Z Q, MIAO G Y, DENG Z F. Global climate change and winter wheat planting in northern region., 2004, 6(1): 9-13. (in Chinese)

    [41] 楊曉光, 劉志娟, 陳阜. 全球氣候變暖對中國種植制度可能影響Ⅰ.氣候變暖對中國種植制度北界和糧食產量可能影響的分析. 中國農業(yè)科學, 2010, 43(2): 329-336.

    YANG X G, LIU Z J, CHEN F. The possible effects of global warming on cropping systems in ChinaⅠ. The possible effects of climate warming on northern limits of cropping systems and crop yields in China., 2010, 43(2): 329-336. (in Chinese)

    [42] 霍治國, 李茂松, 王麗, 溫泉沛, 肖晶晶, 黃大鵬, 王春艷. 氣候變暖對中國農作物病蟲害的影響. 中國農業(yè)科學, 2012, 45(10): 1926-1934.

    HUO Z G, LI M S, WANG L, WEN Q P, XIAO J J, HUANG D P, WANG C Y. Impacts of climate warming on crop diseases and pests in China., 2012, 45(10): 1926-1934. (in Chinese)

    (責任編輯 楊鑫浩)

    Effects of Climate Change on Cropping Pattern and Yield of Summer Maize-Winter Wheat in Huang-Huai-Hai Plain

    SUN XinSu, LONG ZhiWei, SONG GuangPeng, CHEN ChangQing

    (Institute of Applied Ecology, Nanjing Agricultural University, Nanjing 210095)

    【Objective】Global climate is changing significantly characterized by warming. Exploring the influence of climate warming on winter wheat-summer maize cropping system in Huang-Huai-Hai Plain (HHH) to provide a theoretical foundation for formulating reasonable counter measures. 【Method】The meteorological data, phenological data and yield data of HHH during 1992-2013 were used to analyze the effect of climate change on summer maize-winter wheat system by the mathematical statistics analysis method. The meteorological information including the temperature and precipitation data were obtained from China Meteorological Administration. Climate change of HHH area during 1992-2013 were analyzed by monadic linear regression. Through the weighted average of phenological data, the growth period and stubbles changes of summer maize-winter wheat were analyzed. The summer maize-winter wheat yields during 1992-2013 were analyzed by monadic linear regression. What’s more, nonlinear regression analysis and panel data sensitivity analysis were used to analyze the effect of climate change on annual yield of summer maize-winter wheat.【Result】The temperature of HHH showed a rising trend while precipitation generally remained unchange in the 1992-2013, but the regional difference was obvious. Climate change significantly changed the summer maize and winter wheat planting system: The whole growth stages of winter wheat observably shortened by 2 to 5 days in different regions with sowing time delayed. Sowing time of maize showed a difference of southern delay but northern advance while mature time delayed in general. What’s more, there was no significant change of summer maize growth period in the whole area. Stubbles time of summer maize-winter wheat was postponed because of the change of growth period, which led to the waste of climatic resources and land resources. In recent years, the anniversary yields of summer maize and winter wheat showed a significant rising trend in HHH. Climate sensitivity analysis showed that the impact of low temperature, high temperature and average temperature on winter wheat yield was alike, the yield in Jiangsu and Shandong yield decreased while that in Hebei and Henan increased, which the minimum temperature showed a significant effect. In the whole HHH plain, summer maize yield was reduced in different degrees when temperature rising except Henan province because of the difference of temperature rise range and precipitation. There were regional difference of annual yield affected by climate changes due to the regional difference of precipitation. Effect of climate change on annual yield was different that yield in north was increased but decreased in south, so choosing appropriate summer maize varieties is meaningful to guarantee the annual yield. 【Conclusion】Global warming delayed the winter wheat sowing time and shortened winter wheat growing period. It also made summer maize growth period changed that is extended in the north but that is shortened in the south. The prolonged stubbles owing to the shortened growth period led to the waste of climatic resources and land resources. Actual anniversary yield showed a significantly increasing trend in recent years. Effect of increasing temperature and precipitation had regional differences. Generally speaking, temperature rise made maize yield reduced and wheat yield increased, while precipitation made maize yield in northern HHH increased and maize yield reduced in the southern HHH as well as wheat production decreased in different degrees.

    climate change; summer maize-winter wheat; cropping pattern; yield

    2016-10-13;接受日期:2017-05-10

    國家重點基礎研究發(fā)展計劃(2010CB951501)、“十二五”國家科技支撐計劃(2011BAD16B14)

    孫新素,E-mail:15851800979@163.com。通信作者陳長青,E-mail:c0828@126.com

    猜你喜歡
    產量影響
    2022年11月份我國鋅產量同比增長2.9% 鉛產量同比增長5.6%
    是什么影響了滑動摩擦力的大小
    提高玉米產量 膜下滴灌有效
    哪些顧慮影響擔當?
    當代陜西(2021年2期)2021-03-29 07:41:24
    世界致密油產量發(fā)展趨勢
    海水稻產量測評平均產量逐年遞增
    2018年我國主要水果產量按?。▍^(qū)、市)分布
    2018上半年我國PVC產量數據
    聚氯乙烯(2018年9期)2018-02-18 01:11:34
    沒錯,痛經有時也會影響懷孕
    媽媽寶寶(2017年3期)2017-02-21 01:22:28
    擴鏈劑聯用對PETG擴鏈反應與流變性能的影響
    中國塑料(2016年3期)2016-06-15 20:30:00
    国产大屁股一区二区在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 在线免费十八禁| 免费观看人在逋| 最新中文字幕久久久久| 国产老妇女一区| 久久久久久久久中文| 99热这里只有精品一区| 一级黄片播放器| 青春草国产在线视频| 午夜精品在线福利| 国产亚洲精品久久久com| 六月丁香七月| 99久久无色码亚洲精品果冻| 女人十人毛片免费观看3o分钟| 国国产精品蜜臀av免费| 国产精品人妻久久久影院| 亚洲av成人精品一二三区| 午夜福利高清视频| 国产老妇女一区| 日本免费在线观看一区| 中文亚洲av片在线观看爽| 国产毛片a区久久久久| 国产激情偷乱视频一区二区| 国产精品久久电影中文字幕| 国产69精品久久久久777片| 久久99热6这里只有精品| 午夜精品在线福利| 99国产精品一区二区蜜桃av| 国产精品女同一区二区软件| 变态另类丝袜制服| 国产精品国产三级国产av玫瑰| 日本色播在线视频| 自拍偷自拍亚洲精品老妇| 欧美又色又爽又黄视频| 边亲边吃奶的免费视频| 高清视频免费观看一区二区 | 青春草国产在线视频| 国产精品一区www在线观看| 久久久久九九精品影院| 国产精品久久久久久av不卡| 又爽又黄a免费视频| 51国产日韩欧美| 一区二区三区免费毛片| 午夜亚洲福利在线播放| 日韩强制内射视频| 国产高清视频在线观看网站| 久久亚洲精品不卡| 人人妻人人澡欧美一区二区| 亚洲精品aⅴ在线观看| 国产老妇女一区| 国产老妇女一区| 亚洲国产高清在线一区二区三| 女人久久www免费人成看片 | 色吧在线观看| 成人国产麻豆网| 国产午夜精品久久久久久一区二区三区| 亚洲精品aⅴ在线观看| 亚洲成人av在线免费| 91狼人影院| 欧美变态另类bdsm刘玥| 国产又黄又爽又无遮挡在线| 欧美成人午夜免费资源| 久久精品久久久久久久性| 国产人妻一区二区三区在| 91精品伊人久久大香线蕉| 成人午夜高清在线视频| 国产精品精品国产色婷婷| 久久久成人免费电影| 婷婷色麻豆天堂久久 | 观看美女的网站| 日韩大片免费观看网站 | 我要看日韩黄色一级片| 国产av不卡久久| 我要看日韩黄色一级片| 天美传媒精品一区二区| 日本爱情动作片www.在线观看| 嫩草影院精品99| 国产精品爽爽va在线观看网站| 国产成人福利小说| 国产成年人精品一区二区| 亚洲怡红院男人天堂| 国产精品av视频在线免费观看| 亚洲国产精品sss在线观看| 91久久精品国产一区二区成人| 嫩草影院新地址| 69人妻影院| 搡老妇女老女人老熟妇| 91精品伊人久久大香线蕉| 精品久久久久久久末码| 最近手机中文字幕大全| 国产在视频线精品| 亚洲欧美清纯卡通| 一本一本综合久久| 桃色一区二区三区在线观看| 久久草成人影院| 观看免费一级毛片| 亚洲色图av天堂| 波野结衣二区三区在线| 亚洲电影在线观看av| 国产爱豆传媒在线观看| 亚洲精品自拍成人| 中文欧美无线码| eeuss影院久久| 欧美bdsm另类| 国产一区二区在线观看日韩| 精品人妻视频免费看| 久久99精品国语久久久| 汤姆久久久久久久影院中文字幕 | 最近最新中文字幕大全电影3| 最近最新中文字幕大全电影3| 高清在线视频一区二区三区 | 日韩在线高清观看一区二区三区| 亚洲电影在线观看av| 国产高清不卡午夜福利| 一个人看视频在线观看www免费| 久久久久久久亚洲中文字幕| 国产色婷婷99| 秋霞伦理黄片| 99久久精品一区二区三区| av在线亚洲专区| 欧美精品一区二区大全| 午夜日本视频在线| 深夜a级毛片| 国产精品一及| 亚洲va在线va天堂va国产| 国产真实乱freesex| 又爽又黄无遮挡网站| 十八禁国产超污无遮挡网站| 国产伦在线观看视频一区| 亚洲欧美精品自产自拍| 亚洲精品自拍成人| 亚洲欧洲国产日韩| 日韩精品有码人妻一区| 久久久久久久久大av| 精品一区二区免费观看| 国产精品国产三级国产av玫瑰| 国产成人91sexporn| 九九久久精品国产亚洲av麻豆| 又粗又硬又长又爽又黄的视频| 69人妻影院| 欧美xxxx黑人xx丫x性爽| 欧美色视频一区免费| 日韩 亚洲 欧美在线| 欧美成人午夜免费资源| 青春草视频在线免费观看| 真实男女啪啪啪动态图| 精品久久久噜噜| 亚洲丝袜综合中文字幕| 赤兔流量卡办理| 色5月婷婷丁香| 久久久久久久久久成人| 国产成人freesex在线| 国产老妇女一区| 精品熟女少妇av免费看| 日韩欧美 国产精品| 亚洲精品亚洲一区二区| 亚洲精品日韩av片在线观看| 人人妻人人澡欧美一区二区| 亚洲自拍偷在线| www日本黄色视频网| 色噜噜av男人的天堂激情| 成人国产麻豆网| 国产精品国产三级国产av玫瑰| 亚洲自偷自拍三级| 欧美一区二区国产精品久久精品| 天堂√8在线中文| 国产高清国产精品国产三级 | 亚洲真实伦在线观看| 精品久久久久久电影网 | 蜜桃久久精品国产亚洲av| 99热这里只有是精品50| 91狼人影院| av天堂中文字幕网| 亚洲中文字幕一区二区三区有码在线看| 深夜a级毛片| 日韩三级伦理在线观看| 淫秽高清视频在线观看| 能在线免费看毛片的网站| 亚洲欧美一区二区三区国产| 亚洲不卡免费看| 一级毛片aaaaaa免费看小| 亚洲最大成人中文| 在线观看66精品国产| 麻豆久久精品国产亚洲av| 26uuu在线亚洲综合色| 国产亚洲91精品色在线| 嫩草影院精品99| 午夜福利在线观看吧| 乱系列少妇在线播放| 久久草成人影院| 国产高清三级在线| 尾随美女入室| 欧美+日韩+精品| 久久精品91蜜桃| 深夜a级毛片| 国产 一区精品| eeuss影院久久| 亚洲国产欧美在线一区| 中国国产av一级| 亚洲av熟女| 日韩 亚洲 欧美在线| av在线老鸭窝| 中文资源天堂在线| 国产白丝娇喘喷水9色精品| 22中文网久久字幕| 永久免费av网站大全| 成人国产麻豆网| 一区二区三区四区激情视频| 亚洲精品456在线播放app| 日韩欧美在线乱码| 日韩精品青青久久久久久| 亚洲av免费在线观看| 男人和女人高潮做爰伦理| 欧美极品一区二区三区四区| 一个人看视频在线观看www免费| 中文字幕av成人在线电影| 免费电影在线观看免费观看| 午夜福利在线观看吧| 精品99又大又爽又粗少妇毛片| 天天一区二区日本电影三级| 亚洲精品日韩av片在线观看| 一夜夜www| 亚洲国产成人一精品久久久| 日韩,欧美,国产一区二区三区 | 91久久精品国产一区二区成人| 国产精品一区www在线观看| 国产成人一区二区在线| 久久久久久九九精品二区国产| 长腿黑丝高跟| 在线免费十八禁| 国产精品野战在线观看| 亚洲中文字幕日韩| 一边亲一边摸免费视频| 春色校园在线视频观看| 欧美一级a爱片免费观看看| 一区二区三区高清视频在线| 国产亚洲午夜精品一区二区久久 | 国产精品,欧美在线| 亚洲精品影视一区二区三区av| 中文字幕人妻熟人妻熟丝袜美| 一个人看的www免费观看视频| 九九久久精品国产亚洲av麻豆| 日日干狠狠操夜夜爽| 男女视频在线观看网站免费| 日韩亚洲欧美综合| 99久国产av精品| 成人无遮挡网站| 免费在线观看成人毛片| 亚洲成av人片在线播放无| 精品人妻视频免费看| 99久久精品国产国产毛片| 联通29元200g的流量卡| 国产三级在线视频| 国产精品久久久久久av不卡| 十八禁国产超污无遮挡网站| 神马国产精品三级电影在线观看| 视频中文字幕在线观看| 亚洲国产成人一精品久久久| 夫妻性生交免费视频一级片| 少妇的逼好多水| 国产亚洲一区二区精品| 亚洲人成网站高清观看| 国产一级毛片七仙女欲春2| 国产亚洲最大av| 亚洲国产成人一精品久久久| 九九在线视频观看精品| 高清av免费在线| 少妇熟女欧美另类| 免费看光身美女| 男人舔女人下体高潮全视频| 国产精品一区www在线观看| 精品酒店卫生间| 亚洲熟妇中文字幕五十中出| 久久精品久久久久久久性| 国语自产精品视频在线第100页| 九色成人免费人妻av| 女的被弄到高潮叫床怎么办| 在线a可以看的网站| 日韩强制内射视频| 直男gayav资源| 91久久精品电影网| 免费大片18禁| 天堂网av新在线| 亚洲欧美精品专区久久| 国产黄片视频在线免费观看| av黄色大香蕉| 亚洲精品久久久久久婷婷小说 | 18禁在线无遮挡免费观看视频| 天堂网av新在线| av免费观看日本| 美女cb高潮喷水在线观看| 国产精品人妻久久久久久| 少妇丰满av| 亚洲国产精品成人久久小说| 九草在线视频观看| 国产亚洲91精品色在线| 亚洲av电影不卡..在线观看| 久久久久久九九精品二区国产| 女人十人毛片免费观看3o分钟| 免费人成在线观看视频色| 日本午夜av视频| 国产真实乱freesex| 三级经典国产精品| 综合色av麻豆| 国产大屁股一区二区在线视频| 蜜臀久久99精品久久宅男| 亚洲国产最新在线播放| 看黄色毛片网站| 国产私拍福利视频在线观看| 非洲黑人性xxxx精品又粗又长| 在线观看66精品国产| 舔av片在线| 亚洲一区高清亚洲精品| 寂寞人妻少妇视频99o| 一级二级三级毛片免费看| 老女人水多毛片| 亚洲真实伦在线观看| 久久午夜福利片| 国产精品野战在线观看| www日本黄色视频网| 国产老妇伦熟女老妇高清| 国产精品乱码一区二三区的特点| 长腿黑丝高跟| 国产av一区在线观看免费| 日韩一本色道免费dvd| 一级毛片久久久久久久久女| 国产av一区在线观看免费| 男人舔奶头视频| 亚洲国产最新在线播放| 国产午夜精品论理片| 高清日韩中文字幕在线| 麻豆乱淫一区二区| 国产精品.久久久| 一边亲一边摸免费视频| 国产v大片淫在线免费观看| 人体艺术视频欧美日本| 综合色av麻豆| 亚洲人成网站在线播| 在线天堂最新版资源| 国产高清国产精品国产三级 | 老师上课跳d突然被开到最大视频| 波多野结衣高清无吗| 欧美zozozo另类| 久久久久久久亚洲中文字幕| 日日啪夜夜撸| 一区二区三区乱码不卡18| 18+在线观看网站| 久久久成人免费电影| 免费黄色在线免费观看| 国产熟女欧美一区二区| 国产高清国产精品国产三级 | 亚洲色图av天堂| 国产亚洲91精品色在线| 99热这里只有精品一区| 22中文网久久字幕| 一区二区三区四区激情视频| 黄色欧美视频在线观看| 五月伊人婷婷丁香| 99久久精品国产国产毛片| 欧美激情国产日韩精品一区| 日韩三级伦理在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久久成人免费电影| 国产在线一区二区三区精 | 成人鲁丝片一二三区免费| 国产三级在线视频| 如何舔出高潮| 日本wwww免费看| 蜜桃久久精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有精品一区| 亚洲av.av天堂| 国产av一区在线观看免费| 精品一区二区免费观看| 亚洲精品乱码久久久v下载方式| 26uuu在线亚洲综合色| 午夜精品国产一区二区电影 | 高清在线视频一区二区三区 | 女的被弄到高潮叫床怎么办| 亚洲色图av天堂| 长腿黑丝高跟| 九九爱精品视频在线观看| 成人毛片60女人毛片免费| av视频在线观看入口| 天堂网av新在线| 人人妻人人澡欧美一区二区| 国产精华一区二区三区| 欧美激情国产日韩精品一区| 国产探花在线观看一区二区| 国产女主播在线喷水免费视频网站 | 亚洲av男天堂| 干丝袜人妻中文字幕| 精品一区二区三区人妻视频| 国产免费视频播放在线视频 | 亚洲三级黄色毛片| 亚洲av男天堂| 国产久久久一区二区三区| 欧美bdsm另类| 久久久国产成人精品二区| 久久热精品热| 午夜免费激情av| 精品国产露脸久久av麻豆 | 毛片一级片免费看久久久久| 国产成人freesex在线| 国国产精品蜜臀av免费| 直男gayav资源| 亚洲精品成人久久久久久| 韩国av在线不卡| 国产精品久久久久久久电影| 成人午夜高清在线视频| 麻豆成人av视频| 一级二级三级毛片免费看| 波多野结衣巨乳人妻| 亚洲一级一片aⅴ在线观看| 欧美成人精品欧美一级黄| 亚洲av中文av极速乱| 亚洲国产色片| 中文亚洲av片在线观看爽| 亚洲婷婷狠狠爱综合网| 三级经典国产精品| 欧美激情在线99| 免费看光身美女| 嫩草影院精品99| 看黄色毛片网站| 97在线视频观看| 啦啦啦韩国在线观看视频| 看非洲黑人一级黄片| 国产69精品久久久久777片| 22中文网久久字幕| 99久国产av精品国产电影| 日本熟妇午夜| 天堂√8在线中文| 久久久久久久亚洲中文字幕| 成人三级黄色视频| 久久精品久久久久久久性| 亚洲精品自拍成人| 精品人妻熟女av久视频| 国产午夜福利久久久久久| 日本爱情动作片www.在线观看| 成人午夜精彩视频在线观看| 日韩欧美三级三区| 亚洲国产精品久久男人天堂| 亚洲精品乱码久久久久久按摩| 国产精品熟女久久久久浪| 小说图片视频综合网站| 亚洲av电影在线观看一区二区三区 | 精品久久久久久久人妻蜜臀av| 亚洲精品成人久久久久久| 99久国产av精品国产电影| 青春草视频在线免费观看| 久久久久久九九精品二区国产| 国语对白做爰xxxⅹ性视频网站| 亚洲国产精品合色在线| 视频中文字幕在线观看| 欧美日韩国产亚洲二区| 久久人人爽人人爽人人片va| 亚洲美女搞黄在线观看| 最新中文字幕久久久久| 欧美成人免费av一区二区三区| 亚洲精品日韩av片在线观看| 国产一级毛片在线| 寂寞人妻少妇视频99o| 亚洲aⅴ乱码一区二区在线播放| 国产探花在线观看一区二区| 日韩av在线免费看完整版不卡| 欧美xxxx性猛交bbbb| ponron亚洲| 久99久视频精品免费| 国产成人freesex在线| av视频在线观看入口| 色综合色国产| 18禁在线无遮挡免费观看视频| 国语自产精品视频在线第100页| a级毛片免费高清观看在线播放| 亚洲在线观看片| 亚洲自偷自拍三级| 日韩视频在线欧美| 菩萨蛮人人尽说江南好唐韦庄 | 小说图片视频综合网站| 国模一区二区三区四区视频| 免费av不卡在线播放| 最近手机中文字幕大全| 亚洲电影在线观看av| 亚洲av免费高清在线观看| 国产免费福利视频在线观看| 日韩欧美精品免费久久| 中文亚洲av片在线观看爽| 成人无遮挡网站| 全区人妻精品视频| 国产精品熟女久久久久浪| 久久精品国产亚洲av涩爱| 搡老妇女老女人老熟妇| 国产精品乱码一区二三区的特点| 亚洲av电影在线观看一区二区三区 | 中文字幕av成人在线电影| 春色校园在线视频观看| 国产熟女欧美一区二区| 春色校园在线视频观看| 全区人妻精品视频| 国产免费视频播放在线视频 | 日本-黄色视频高清免费观看| 午夜免费激情av| av福利片在线观看| 日韩av在线大香蕉| 日本熟妇午夜| 中文字幕熟女人妻在线| 国产 一区 欧美 日韩| 日本午夜av视频| 欧美3d第一页| 国产三级中文精品| 午夜福利在线观看吧| 亚洲国产精品合色在线| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久午夜电影| 免费观看精品视频网站| 91午夜精品亚洲一区二区三区| 国产真实伦视频高清在线观看| 少妇人妻精品综合一区二区| 日本黄色视频三级网站网址| 美女cb高潮喷水在线观看| 97超碰精品成人国产| h日本视频在线播放| av在线天堂中文字幕| 在线播放无遮挡| 国产精品一及| 天堂影院成人在线观看| 国产一区二区在线观看日韩| 1024手机看黄色片| 一级黄色大片毛片| 18禁在线播放成人免费| 国产成人午夜福利电影在线观看| 九色成人免费人妻av| 一级毛片aaaaaa免费看小| 亚洲精品国产成人久久av| 久久久久免费精品人妻一区二区| 综合色丁香网| 看免费成人av毛片| 亚洲婷婷狠狠爱综合网| h日本视频在线播放| 国语对白做爰xxxⅹ性视频网站| 日本-黄色视频高清免费观看| 麻豆成人午夜福利视频| 日产精品乱码卡一卡2卡三| 女人被狂操c到高潮| 精华霜和精华液先用哪个| 亚洲久久久久久中文字幕| 久久精品久久久久久噜噜老黄 | av.在线天堂| 欧美色视频一区免费| 日本五十路高清| 亚洲人成网站在线播| 国产男人的电影天堂91| 99久久精品一区二区三区| 国产精品女同一区二区软件| 黄色一级大片看看| 嫩草影院入口| 精品熟女少妇av免费看| 在线免费观看不下载黄p国产| 精品人妻视频免费看| 高清在线视频一区二区三区 | 久久久久网色| 欧美又色又爽又黄视频| 久久精品国产亚洲网站| 91aial.com中文字幕在线观看| 老司机福利观看| 亚洲国产精品专区欧美| 午夜精品一区二区三区免费看| 国产一区二区亚洲精品在线观看| 亚洲精华国产精华液的使用体验| 久久精品综合一区二区三区| 麻豆成人av视频| 深爱激情五月婷婷| 不卡视频在线观看欧美| 好男人视频免费观看在线| 亚洲精品aⅴ在线观看| 啦啦啦啦在线视频资源| 免费av毛片视频| 国产精品野战在线观看| 我要看日韩黄色一级片| 色5月婷婷丁香| 一个人观看的视频www高清免费观看| 波野结衣二区三区在线| 欧美成人a在线观看| 国产一级毛片七仙女欲春2| 麻豆乱淫一区二区| 精品国产三级普通话版| 亚洲av成人精品一区久久| 嫩草影院精品99| www.av在线官网国产| 亚洲四区av| 国产三级在线视频| 99久久成人亚洲精品观看| 变态另类丝袜制服| 日日摸夜夜添夜夜添av毛片| 搡老妇女老女人老熟妇| 九九在线视频观看精品| 少妇人妻精品综合一区二区| 97人妻精品一区二区三区麻豆| 国产精品蜜桃在线观看| 我的女老师完整版在线观看| 国产免费又黄又爽又色| 国产视频首页在线观看| 黄色一级大片看看| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 久久久国产成人精品二区| 久久久成人免费电影| 联通29元200g的流量卡| 免费看美女性在线毛片视频| 亚洲av二区三区四区| 好男人在线观看高清免费视频| 深夜a级毛片| 日韩视频在线欧美| 搡老妇女老女人老熟妇|