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

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

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

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

    ?

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

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

    (南京農(nóng)業(yè)大學(xué)應(yīng)用生態(tài)研究所,南京210095)

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

    氣候變化;夏玉米-冬小麥;種植模式;產(chǎn)量

    0 引言

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

    1 材料與方法

    1.1 研究區(qū)概況

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

    1.2 資料來源

    數(shù)據(jù)資料主要包括氣象、物候和產(chǎn)量數(shù)據(jù)3個方面。

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

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

    圖2 氣象站點(diǎn)分布圖

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

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

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

    1.3 分析方法

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

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

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

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

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

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

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

    2 結(jié)果

    2.1 夏玉米-冬小麥生長季內(nèi)氣候變化

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

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

    以1992—1996年和2009—2013年首尾兩個時間段表征的黃淮海1992—2013年間不同省(市、區(qū)域)的氣候因子變化情況(表1)。從溫度因子來看,各省、市、區(qū)域除冬小麥生長季最高溫度出現(xiàn)小幅下降外,其他溫度因子均表現(xiàn)為不同程度的增加,溫度上升0.1—0.9℃不等,尤其以最低溫度的上升最為明顯,其中山東省夏玉米生長季最低溫度上升達(dá)0.94℃,溫度升高較為明顯。整個研究區(qū)來看,夏玉米生長季最低溫度、最高溫度和平均溫度均表現(xiàn)為上升的趨勢,分別增加了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ū)冬小麥播種時間均呈現(xiàn)不同程度的推遲,其中山東省推遲達(dá)6 d,而冬小麥?zhǔn)斋@期因不同地區(qū)播種時間推遲的程度不同體現(xiàn)出提前或推遲趨勢,但就冬小麥整個生育進(jìn)程來看,黃淮海各地區(qū)均表現(xiàn)為生育期縮短的趨勢,生育時長平均縮短3 d左右。夏玉米生育期推移情況呈現(xiàn)區(qū)域差異,東南部的江蘇及安徽北部和山東省夏玉米播種時間推遲,而位于黃淮海西部和北部的河南河北則呈現(xiàn)小幅度的提前。夏玉米收獲時間除河南省以外,均呈現(xiàn)推遲趨勢。整個地區(qū)總生育時長未發(fā)生顯著變化,但區(qū)域間的差異較為明顯,這可能與降水的區(qū)域差異有關(guān)。夏玉米以及冬小麥生育時間的變化,顯著改變了夏玉米-冬小麥種植模式的茬口時間,冬小麥生育時長的顯著縮短,使茬口時間出現(xiàn)不同程度的延長,導(dǎo)致了土地和熱量資源的浪費(fèi)。

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

    平均值±標(biāo)準(zhǔn)誤average±standard error

    2.3 夏玉米-冬小麥周年單產(chǎn)變化

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

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

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

    WS:小麥播種期Wheat sowing;WH:小麥?zhǔn)斋@期Wheat harvest;MS:玉米播種期Maize sowing;MH:玉米收獲期Maize harvest

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

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

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

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

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

    3 討論

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

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

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

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

    4 結(jié)論

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

    References

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

    [2] 繆啟龍, 丁園圓, 王勇, 段春鋒. 氣候變暖對中國熱量資源分布的影響分析. 自然資源學(xué)報, 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] 王斌, 顧蘊(yùn)倩, 劉雪, 羅衛(wèi)紅, 戴劍鋒, 張巍, 亓春杰. 中國冬小麥種植區(qū)光熱資源及其配比的時空演變特征分析. 中國農(nóng)業(yè)科學(xué), 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] 潘根興, 高民, 胡國華, 魏欽平, 楊曉光, 張文忠, 周廣勝, 鄒建文. 氣候變化對中國農(nóng)業(yè)生產(chǎn)的影響. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 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] 李正國, 楊鵬, 唐華俊, 吳文斌, 陳仲新, 周清波, 鄒金秋, 張莉. 氣候變化背景下東北三省主要作物典型物候期變化趨勢分析. 中國農(nóng)業(yè)科學(xué), 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ā)育和產(chǎn)量的影響. 大豆科學(xué), 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)農(nóng)業(yè)學(xué)報, 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ū)水稻生育期的影響. 中國農(nóng)業(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īng). 中國農(nóng)業(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] 熊偉, 楊婕, 林而達(dá), 許吟隆. 未來不同氣候變化情景下我國玉米產(chǎn)量的初步預(yù)測. 地球科學(xué)進(jìn)展, 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] 居輝, 熊偉, 許吟隆, 林而達(dá). 氣候變化對我國小麥產(chǎn)量的影響. 作物學(xué)報, 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)農(nóng)業(yè)學(xué)報, 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] 龐新生. 缺失數(shù)據(jù)處理方法的比較. 統(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)學(xué)報, 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] 侯瓊, 郭瑞清, 楊麗桃. 內(nèi)蒙古氣候變化及其對主要農(nóng)作物的影響. 中國農(nóng)業(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] 王力, 李鳳霞, 徐維新, 肖宏斌, 蘇文將, 周萬福. 氣候變化對柴達(dá)木灌區(qū)小麥生育期的影響. 中國農(nóng)業(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)模擬模型. 江蘇農(nóng)業(yè)學(xué)報, 2000, 16(1): 15-21.

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

    [28] 房世波, 譚凱炎, 任三學(xué), 張新時. 氣候變暖對冬小麥生長和產(chǎn)量影響的大田實驗研究. 中國科學(xué)(地球科學(xué)), 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ū)氣候變暖對夏玉米生育進(jìn)程及產(chǎn)量的影響. 中國農(nóng)業(yè)科學(xué), 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)農(nóng)業(yè)學(xué)報, 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] 胡實, 莫興國, 林忠輝. 氣候變化對黃淮海平原冬小麥產(chǎn)量和耗水的影響及品種適應(yīng)性評估. 應(yīng)用生態(tài)學(xué)報, 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ā)育和產(chǎn)量的影響. 生態(tài)學(xué)報, 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] 張吉旺, 董樹亭, 王空軍, 胡昌浩, 劉鵬. 大田增溫對夏玉米產(chǎn)量和品質(zhì)的影響. 應(yīng)用生態(tài)學(xué)報, 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] 王琪, 馬樹慶, 郭建平, 張鐵林, 于海, 徐麗萍. 溫度對玉米生長和產(chǎn)量的影響. 生態(tài)學(xué)雜志, 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] 朱新玉, 張竟竟, 趙文亮, 史本林. 河南省商丘市氣候變化對玉米產(chǎn)量的影響. 湖北農(nóng)業(yè)科學(xué), 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年氣候變化響應(yīng)的敏感性研究. 中國生態(tài)農(nóng)業(yè)學(xué)報, 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] 郝志新, 鄭景云, 陶向新. 氣候增暖背景下的冬小麥種植北界研究——以遼寧省為例. 地理科學(xué)進(jìn)展, 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] 高志強(qiáng), 苗果園, 鄧志鋒. 全球氣候變化與冬麥北移研究. 中國農(nóng)業(yè)科技導(dǎo)報, 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] 楊曉光, 劉志娟, 陳阜. 全球氣候變暖對中國種植制度可能影響Ⅰ.氣候變暖對中國種植制度北界和糧食產(chǎn)量可能影響的分析. 中國農(nóng)業(yè)科學(xué), 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] 霍治國, 李茂松, 王麗, 溫泉沛, 肖晶晶, 黃大鵬, 王春艷. 氣候變暖對中國農(nóng)作物病蟲害的影響. 中國農(nóng)業(yè)科學(xué), 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)

    (責(zé)任編輯 楊鑫浩)

    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

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

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

    猜你喜歡
    產(chǎn)量影響
    2022年11月份我國鋅產(chǎn)量同比增長2.9% 鉛產(chǎn)量同比增長5.6%
    是什么影響了滑動摩擦力的大小
    提高玉米產(chǎn)量 膜下滴灌有效
    哪些顧慮影響擔(dān)當(dāng)?
    世界致密油產(chǎn)量發(fā)展趨勢
    海水稻產(chǎn)量測評平均產(chǎn)量逐年遞增
    2018年我國主要水果產(chǎn)量按?。▍^(qū)、市)分布
    2018上半年我國PVC產(chǎn)量數(shù)據(jù)
    聚氯乙烯(2018年9期)2018-02-18 01:11:34
    沒錯,痛經(jīng)有時也會影響懷孕
    媽媽寶寶(2017年3期)2017-02-21 01:22:28
    擴(kuò)鏈劑聯(lián)用對PETG擴(kuò)鏈反應(yīng)與流變性能的影響
    中國塑料(2016年3期)2016-06-15 20:30:00
    亚洲精品aⅴ在线观看| 新久久久久国产一级毛片| 久久精品久久精品一区二区三区| 99热国产这里只有精品6| 亚洲久久久国产精品| 夜夜爽夜夜爽视频| 一二三四在线观看免费中文在 | 大话2 男鬼变身卡| 考比视频在线观看| a 毛片基地| 亚洲国产av影院在线观看| 一级毛片黄色毛片免费观看视频| 欧美3d第一页| 精品久久久精品久久久| 日本av手机在线免费观看| 七月丁香在线播放| 亚洲精品,欧美精品| 午夜免费男女啪啪视频观看| 18禁观看日本| 日本av手机在线免费观看| 精品一区二区免费观看| 肉色欧美久久久久久久蜜桃| 欧美日本中文国产一区发布| 国产精品久久久av美女十八| 国产日韩欧美视频二区| 毛片一级片免费看久久久久| 日韩三级伦理在线观看| 伊人亚洲综合成人网| 日韩,欧美,国产一区二区三区| 国产一区有黄有色的免费视频| 飞空精品影院首页| 亚洲成人手机| 伦理电影免费视频| 一本色道久久久久久精品综合| 免费人妻精品一区二区三区视频| 国产欧美亚洲国产| 午夜影院在线不卡| 韩国av在线不卡| 亚洲情色 制服丝袜| 亚洲av.av天堂| 成人黄色视频免费在线看| 在线天堂最新版资源| 久久久久久久亚洲中文字幕| 欧美日韩av久久| 美女xxoo啪啪120秒动态图| 18禁裸乳无遮挡动漫免费视频| 日韩不卡一区二区三区视频在线| 午夜福利在线观看免费完整高清在| 免费人成在线观看视频色| 久久久久久久国产电影| 国产精品99久久99久久久不卡 | 国产亚洲精品第一综合不卡 | 香蕉丝袜av| 夜夜爽夜夜爽视频| 精品人妻一区二区三区麻豆| 日韩一区二区三区影片| 亚洲,一卡二卡三卡| 一级毛片电影观看| 国产一区二区在线观看日韩| 亚洲精品乱码久久久久久按摩| 人人妻人人爽人人添夜夜欢视频| 色5月婷婷丁香| 精品亚洲乱码少妇综合久久| 熟女人妻精品中文字幕| av国产精品久久久久影院| 在线观看三级黄色| 精品一区二区三区四区五区乱码 | 巨乳人妻的诱惑在线观看| 成年人午夜在线观看视频| av国产久精品久网站免费入址| 亚洲色图综合在线观看| 在线 av 中文字幕| 丝袜在线中文字幕| 久久久久久久国产电影| 丰满饥渴人妻一区二区三| 中文乱码字字幕精品一区二区三区| 亚洲精品av麻豆狂野| 一级毛片电影观看| 国产成人91sexporn| 制服丝袜香蕉在线| 国产成人精品福利久久| av视频免费观看在线观看| 亚洲图色成人| 婷婷成人精品国产| 国语对白做爰xxxⅹ性视频网站| 久久综合国产亚洲精品| 欧美性感艳星| 国产片内射在线| 欧美另类一区| 午夜免费鲁丝| 一区二区三区四区激情视频| 久久精品国产综合久久久 | 看免费成人av毛片| 午夜福利在线观看免费完整高清在| 欧美日韩综合久久久久久| 成人综合一区亚洲| 97在线人人人人妻| 亚洲国产成人一精品久久久| 美女国产高潮福利片在线看| 亚洲精品国产av成人精品| 免费大片18禁| 亚洲情色 制服丝袜| 欧美国产精品一级二级三级| 高清av免费在线| 亚洲精品乱久久久久久| 亚洲欧洲国产日韩| 街头女战士在线观看网站| 国产女主播在线喷水免费视频网站| 国产有黄有色有爽视频| 97超碰精品成人国产| 亚洲丝袜综合中文字幕| 香蕉精品网在线| 在线免费观看不下载黄p国产| 久久韩国三级中文字幕| √禁漫天堂资源中文www| 亚洲国产av影院在线观看| 精品国产乱码久久久久久小说| 亚洲精品日韩在线中文字幕| 日韩一本色道免费dvd| 国产有黄有色有爽视频| av有码第一页| 两个人免费观看高清视频| 久久国产精品男人的天堂亚洲 | 国产熟女午夜一区二区三区| 波多野结衣一区麻豆| 最新的欧美精品一区二区| 亚洲激情五月婷婷啪啪| 一级,二级,三级黄色视频| 99精国产麻豆久久婷婷| 国产毛片在线视频| 日韩在线高清观看一区二区三区| 国产精品不卡视频一区二区| 亚洲国产看品久久| 在线观看人妻少妇| 男人添女人高潮全过程视频| 日韩电影二区| 精品视频人人做人人爽| 日韩欧美一区视频在线观看| av视频免费观看在线观看| 国产精品不卡视频一区二区| 午夜福利乱码中文字幕| 亚洲丝袜综合中文字幕| 午夜视频国产福利| 久久亚洲国产成人精品v| 午夜免费男女啪啪视频观看| 两个人免费观看高清视频| 精品人妻在线不人妻| 久久毛片免费看一区二区三区| 色5月婷婷丁香| 成人亚洲欧美一区二区av| 大陆偷拍与自拍| 波野结衣二区三区在线| 国产成人精品久久久久久| 国产精品秋霞免费鲁丝片| 欧美日韩精品成人综合77777| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品国产亚洲| 中文精品一卡2卡3卡4更新| 国产精品三级大全| 男女国产视频网站| 久久青草综合色| 国产成人免费无遮挡视频| 一区二区三区四区激情视频| av国产久精品久网站免费入址| 又黄又粗又硬又大视频| 国产成人免费无遮挡视频| 啦啦啦视频在线资源免费观看| 亚洲欧美日韩卡通动漫| 91在线精品国自产拍蜜月| 高清不卡的av网站| 91精品国产国语对白视频| 少妇猛男粗大的猛烈进出视频| 精品第一国产精品| 91精品伊人久久大香线蕉| 精品少妇内射三级| 一本久久精品| 国产精品久久久久成人av| 婷婷色综合www| 日本欧美国产在线视频| 欧美精品亚洲一区二区| 亚洲五月色婷婷综合| 亚洲人与动物交配视频| 国产精品国产三级国产专区5o| 欧美成人午夜免费资源| 精品人妻一区二区三区麻豆| 国产又爽黄色视频| 国产色婷婷99| 人妻 亚洲 视频| 久久久欧美国产精品| 高清视频免费观看一区二区| 十八禁网站网址无遮挡| 看免费成人av毛片| 肉色欧美久久久久久久蜜桃| 亚洲精品aⅴ在线观看| 国产精品无大码| 亚洲国产av影院在线观看| 久久久久久伊人网av| 大片电影免费在线观看免费| 成人漫画全彩无遮挡| 久久久久久久精品精品| 最近手机中文字幕大全| 久久久久国产精品人妻一区二区| 在线免费观看不下载黄p国产| 国产爽快片一区二区三区| 亚洲国产毛片av蜜桃av| 欧美激情国产日韩精品一区| 在线天堂最新版资源| 亚洲精品乱码久久久久久按摩| 免费观看a级毛片全部| 亚洲国产av影院在线观看| 夫妻午夜视频| 视频中文字幕在线观看| 亚洲国产色片| 欧美成人午夜免费资源| 国产高清三级在线| 尾随美女入室| 大陆偷拍与自拍| 国产女主播在线喷水免费视频网站| 香蕉国产在线看| 亚洲精品乱久久久久久| 成人漫画全彩无遮挡| 亚洲高清免费不卡视频| 又黄又粗又硬又大视频| 一边摸一边做爽爽视频免费| 亚洲少妇的诱惑av| 亚洲精华国产精华液的使用体验| 90打野战视频偷拍视频| 纵有疾风起免费观看全集完整版| 制服诱惑二区| 26uuu在线亚洲综合色| 一本色道久久久久久精品综合| 人人澡人人妻人| 久久国内精品自在自线图片| 黄网站色视频无遮挡免费观看| 国产乱来视频区| 精品亚洲成a人片在线观看| 97超碰精品成人国产| 成人亚洲精品一区在线观看| 久久人妻熟女aⅴ| 哪个播放器可以免费观看大片| 色哟哟·www| 99久久综合免费| 最近2019中文字幕mv第一页| 男女午夜视频在线观看 | 日韩不卡一区二区三区视频在线| 久久久精品94久久精品| 91在线精品国自产拍蜜月| 夜夜骑夜夜射夜夜干| 妹子高潮喷水视频| 精品人妻一区二区三区麻豆| 亚洲国产看品久久| 不卡视频在线观看欧美| 99国产综合亚洲精品| 91国产中文字幕| 我的女老师完整版在线观看| 制服人妻中文乱码| 亚洲精品av麻豆狂野| 99热6这里只有精品| 国产色婷婷99| 国产精品.久久久| 1024视频免费在线观看| 各种免费的搞黄视频| 日韩伦理黄色片| 国产精品免费大片| 久久午夜福利片| 精品99又大又爽又粗少妇毛片| 美女中出高潮动态图| 亚洲经典国产精华液单| 18在线观看网站| 黑人高潮一二区| 高清不卡的av网站| 少妇人妻精品综合一区二区| 黑人巨大精品欧美一区二区蜜桃 | 久久久久视频综合| 美女内射精品一级片tv| 99热这里只有是精品在线观看| 极品少妇高潮喷水抽搐| 国产欧美亚洲国产| 亚洲一区二区三区欧美精品| 七月丁香在线播放| av在线播放精品| 国产黄色视频一区二区在线观看| 亚洲av在线观看美女高潮| 日本黄色日本黄色录像| 国产精品成人在线| 黄片播放在线免费| 曰老女人黄片| 成年动漫av网址| 九九在线视频观看精品| 国产极品粉嫩免费观看在线| 成人影院久久| 久久久久国产精品人妻一区二区| 国产日韩欧美在线精品| 国产精品久久久久久久电影| 中文字幕人妻熟女乱码| 亚洲国产精品一区三区| 免费黄色在线免费观看| 中国美白少妇内射xxxbb| 汤姆久久久久久久影院中文字幕| 精品午夜福利在线看| 中文字幕精品免费在线观看视频 | 久久99精品国语久久久| 激情五月婷婷亚洲| 在线免费观看不下载黄p国产| 午夜免费男女啪啪视频观看| 夫妻性生交免费视频一级片| 亚洲美女黄色视频免费看| 制服人妻中文乱码| 国产精品熟女久久久久浪| 涩涩av久久男人的天堂| 成年人免费黄色播放视频| 国产熟女午夜一区二区三区| 91在线精品国自产拍蜜月| 黄色 视频免费看| 久久亚洲国产成人精品v| 男女啪啪激烈高潮av片| 久久久久视频综合| 美国免费a级毛片| 新久久久久国产一级毛片| 欧美激情极品国产一区二区三区 | 五月天丁香电影| 国产免费又黄又爽又色| 国产成人精品久久久久久| 亚洲国产欧美在线一区| 亚洲精品国产av成人精品| 精品人妻一区二区三区麻豆| 精品国产露脸久久av麻豆| 亚洲国产av新网站| 最近中文字幕高清免费大全6| 韩国高清视频一区二区三区| 久久狼人影院| 99视频精品全部免费 在线| 五月伊人婷婷丁香| 9色porny在线观看| 51国产日韩欧美| 欧美变态另类bdsm刘玥| 国产黄色免费在线视频| 精品一区在线观看国产| 精品一区二区三区视频在线| 乱人伦中国视频| 黄色配什么色好看| 国产熟女午夜一区二区三区| 久久ye,这里只有精品| 亚洲精品第二区| 国产成人午夜福利电影在线观看| 婷婷色麻豆天堂久久| 中国国产av一级| 久久精品久久久久久久性| 亚洲欧美成人精品一区二区| 日本与韩国留学比较| 久久久精品94久久精品| 麻豆精品久久久久久蜜桃| av视频免费观看在线观看| 天天操日日干夜夜撸| 一级片'在线观看视频| 日韩 亚洲 欧美在线| 久久精品aⅴ一区二区三区四区 | 久久久久久人妻| 精品国产露脸久久av麻豆| 亚洲av.av天堂| 校园人妻丝袜中文字幕| 国产精品久久久久久精品古装| av视频免费观看在线观看| 狠狠婷婷综合久久久久久88av| 免费不卡的大黄色大毛片视频在线观看| 国产在视频线精品| av黄色大香蕉| 大陆偷拍与自拍| 色吧在线观看| 丰满少妇做爰视频| 日本91视频免费播放| 免费高清在线观看日韩| 少妇人妻久久综合中文| 在线观看免费高清a一片| 久久久久久人妻| 色网站视频免费| 十分钟在线观看高清视频www| 曰老女人黄片| 久久久精品免费免费高清| 成人毛片60女人毛片免费| 男女高潮啪啪啪动态图| 建设人人有责人人尽责人人享有的| 人人妻人人爽人人添夜夜欢视频| 最黄视频免费看| 99久国产av精品国产电影| 一级爰片在线观看| 亚洲精品乱码久久久久久按摩| 黑人欧美特级aaaaaa片| 久久99热6这里只有精品| 久久亚洲国产成人精品v| 最黄视频免费看| 一级毛片我不卡| 青春草视频在线免费观看| 日韩一本色道免费dvd| 久久人人爽av亚洲精品天堂| 久久久久久久国产电影| 香蕉丝袜av| 国产成人91sexporn| 性色avwww在线观看| 秋霞伦理黄片| 亚洲av中文av极速乱| 黑人欧美特级aaaaaa片| 中文字幕制服av| 国产高清不卡午夜福利| 女人精品久久久久毛片| 丝袜脚勾引网站| 免费在线观看完整版高清| 亚洲国产av影院在线观看| 久久国内精品自在自线图片| 久久狼人影院| 街头女战士在线观看网站| 国产1区2区3区精品| 亚洲国产毛片av蜜桃av| 伊人亚洲综合成人网| 黄色怎么调成土黄色| 欧美成人午夜免费资源| 中文字幕人妻丝袜制服| 欧美国产精品一级二级三级| 啦啦啦啦在线视频资源| 夫妻性生交免费视频一级片| 国产精品免费大片| 婷婷成人精品国产| 久久久久国产网址| 亚洲成国产人片在线观看| 九色亚洲精品在线播放| 亚洲,欧美精品.| av片东京热男人的天堂| 国产精品99久久99久久久不卡 | 亚洲一级一片aⅴ在线观看| 下体分泌物呈黄色| 久久久久久久久久成人| 亚洲伊人久久精品综合| 51国产日韩欧美| 欧美最新免费一区二区三区| 2018国产大陆天天弄谢| 天堂俺去俺来也www色官网| 黑丝袜美女国产一区| 又黄又粗又硬又大视频| 中文字幕亚洲精品专区| 亚洲经典国产精华液单| 精品久久国产蜜桃| 如何舔出高潮| 大香蕉久久网| 欧美另类一区| 51国产日韩欧美| 岛国毛片在线播放| 大片电影免费在线观看免费| 51国产日韩欧美| 99精国产麻豆久久婷婷| 在线观看国产h片| 乱码一卡2卡4卡精品| 免费高清在线观看视频在线观看| 狂野欧美激情性xxxx在线观看| 高清欧美精品videossex| kizo精华| 色吧在线观看| 免费黄色在线免费观看| 99视频精品全部免费 在线| 久久人人爽av亚洲精品天堂| 美女中出高潮动态图| 久久久久精品久久久久真实原创| 亚洲国产精品成人久久小说| 欧美激情极品国产一区二区三区 | 丰满饥渴人妻一区二区三| 一级毛片黄色毛片免费观看视频| 欧美性感艳星| 丝袜喷水一区| 午夜av观看不卡| 熟女人妻精品中文字幕| 少妇高潮的动态图| 中文字幕最新亚洲高清| 乱码一卡2卡4卡精品| 欧美精品一区二区免费开放| 免费日韩欧美在线观看| 亚洲精品美女久久av网站| 精品午夜福利在线看| 女人被躁到高潮嗷嗷叫费观| 22中文网久久字幕| 一区二区三区乱码不卡18| 亚洲,欧美精品.| 日日撸夜夜添| 欧美 日韩 精品 国产| 男女下面插进去视频免费观看 | 亚洲精华国产精华液的使用体验| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人一二三区av| 国产永久视频网站| 久久精品国产综合久久久 | 国产精品一区二区在线不卡| 午夜影院在线不卡| 日韩不卡一区二区三区视频在线| 欧美xxxx性猛交bbbb| 五月开心婷婷网| 成人18禁高潮啪啪吃奶动态图| 国产又爽黄色视频| 精品一区二区三区四区五区乱码 | 午夜免费鲁丝| 少妇被粗大的猛进出69影院 | 久久韩国三级中文字幕| 精品一品国产午夜福利视频| 七月丁香在线播放| 999精品在线视频| 亚洲一码二码三码区别大吗| 男男h啪啪无遮挡| 欧美激情国产日韩精品一区| 国产亚洲午夜精品一区二区久久| 亚洲伊人色综图| av福利片在线| 日韩一区二区三区影片| 亚洲欧美一区二区三区黑人 | 亚洲欧洲精品一区二区精品久久久 | 国产亚洲欧美精品永久| 三上悠亚av全集在线观看| 国产黄色免费在线视频| 久久国内精品自在自线图片| 亚洲经典国产精华液单| 天堂8中文在线网| 日韩欧美一区视频在线观看| 欧美精品国产亚洲| 午夜福利视频在线观看免费| 国产麻豆69| videos熟女内射| 自拍欧美九色日韩亚洲蝌蚪91| tube8黄色片| 国产欧美另类精品又又久久亚洲欧美| 午夜精品国产一区二区电影| 大香蕉久久成人网| 一边摸一边做爽爽视频免费| 欧美精品人与动牲交sv欧美| 国产精品熟女久久久久浪| 日韩,欧美,国产一区二区三区| 欧美性感艳星| 国产又色又爽无遮挡免| 大片电影免费在线观看免费| 国产日韩欧美亚洲二区| 国产国语露脸激情在线看| 婷婷成人精品国产| 亚洲精品久久久久久婷婷小说| 亚洲人与动物交配视频| 欧美精品一区二区大全| 亚洲欧美清纯卡通| 亚洲成色77777| 国产1区2区3区精品| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 亚洲欧美成人精品一区二区| 精品久久国产蜜桃| 亚洲三级黄色毛片| 亚洲精品乱久久久久久| 国产乱人偷精品视频| 国产一区二区在线观看日韩| 超碰97精品在线观看| 亚洲精品日韩在线中文字幕| 9热在线视频观看99| 美女脱内裤让男人舔精品视频| 亚洲高清免费不卡视频| 少妇高潮的动态图| 色视频在线一区二区三区| 中国美白少妇内射xxxbb| 亚洲精品国产色婷婷电影| 巨乳人妻的诱惑在线观看| 亚洲精品色激情综合| tube8黄色片| 国产一区亚洲一区在线观看| 免费黄频网站在线观看国产| 99热网站在线观看| 极品少妇高潮喷水抽搐| 久久久久久伊人网av| 97在线人人人人妻| 精品久久久久久电影网| 亚洲国产精品国产精品| 一区二区av电影网| 伊人久久国产一区二区| av片东京热男人的天堂| 亚洲精品一区蜜桃| www.色视频.com| av黄色大香蕉| 日日撸夜夜添| 午夜福利视频在线观看免费| 亚洲av男天堂| 久久久a久久爽久久v久久| 成人国语在线视频| 久久久国产欧美日韩av| 十分钟在线观看高清视频www| 天天影视国产精品| 一边亲一边摸免费视频| 午夜福利乱码中文字幕| √禁漫天堂资源中文www| 18+在线观看网站| 欧美性感艳星| 国产在线视频一区二区| 男男h啪啪无遮挡| 老司机影院毛片| 最近手机中文字幕大全| 免费观看av网站的网址| 热99国产精品久久久久久7| 永久免费av网站大全| 国产精品99久久99久久久不卡 | √禁漫天堂资源中文www| 少妇的逼水好多| 久久热在线av| 久久久久国产网址| 美女内射精品一级片tv| 美女脱内裤让男人舔精品视频| 婷婷色综合www| 亚洲一码二码三码区别大吗| 国产成人aa在线观看| 亚洲精品视频女| 免费看av在线观看网站| 国产一区二区在线观看日韩| av线在线观看网站| 一级毛片 在线播放| 午夜影院在线不卡| 亚洲欧洲日产国产| 伦精品一区二区三区|