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

    東北三省西部春玉米適應(yīng)氣候變化的高產(chǎn)高效灌溉方案分析

    2020-11-09 05:19:12黃秋婉劉志娟楊曉光白帆劉濤張鎮(zhèn)濤孫爽趙錦
    中國農(nóng)業(yè)科學(xué) 2020年21期
    關(guān)鍵詞:氣候區(qū)東北三省降水

    黃秋婉,劉志娟,楊曉光,白帆,劉濤,張鎮(zhèn)濤,孫爽,趙錦

    東北三省西部春玉米適應(yīng)氣候變化的高產(chǎn)高效灌溉方案分析

    黃秋婉,劉志娟,楊曉光,白帆,劉濤,張鎮(zhèn)濤,孫爽,趙錦

    (中國農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,北京 100193)

    【】東北三省是我國重要的商品糧生產(chǎn)基地之一,同時也是對氣候變化最敏感的地區(qū),因此明確氣候變化背景下東北三省西部干旱區(qū)春玉米的適宜灌溉措施,對于當(dāng)?shù)卮河衩赘弋a(chǎn)穩(wěn)產(chǎn)水資源高效利用有重要意義。依據(jù)春玉米生長季積溫和水分虧缺率將東北三省春玉米潛在種植區(qū)劃分為10個氣候區(qū),以東北三省西部5個水分虧缺率>0的氣候區(qū)為研究區(qū)域,基于1981—2017年的氣象資料、農(nóng)業(yè)氣象觀測站春玉米試驗數(shù)據(jù)和土壤資料,對農(nóng)業(yè)生產(chǎn)系統(tǒng)模型(APSIM-Maize)相關(guān)參數(shù)進(jìn)行調(diào)試并驗證其適用性。設(shè)置不同灌溉情景,利用驗證后的模型模擬各氣候區(qū)不同灌溉情景下的春玉米產(chǎn)量,結(jié)合水分利用效率明確各氣候區(qū)不同年代的適宜灌溉措施及產(chǎn)量提升幅度。(1)近37年(1981—2017年)5個氣候區(qū)有效積溫均呈顯著上升趨勢,降水量呈下降趨勢。從過去37年平均來看,第一和第三氣候區(qū)降水對春玉米產(chǎn)量的限制程度較小,分別為0—27%和0—9%,通過灌溉對產(chǎn)量提升的貢獻(xiàn)較小,但能有效提高產(chǎn)量的穩(wěn)定性(可使第一氣候區(qū)產(chǎn)量變異系數(shù)由0.24降低到0.11,第三氣候區(qū)產(chǎn)量變異系數(shù)由0.14降低到0.12);第五、七和九氣候區(qū)降水對春玉米產(chǎn)量的限制程度較大,分別為27%—69%,15%—35%,31%—51%,灌溉不僅可提升當(dāng)?shù)赜衩桩a(chǎn)量,同時可使3個氣候區(qū)的產(chǎn)量變異系數(shù)由0.54降低到0.15,0.46降低到0.13,0.65降低到0.13。表明在東北三省西部干旱區(qū)通過灌溉能達(dá)到高產(chǎn)穩(wěn)產(chǎn)的目的。(2)第一和第三氣候區(qū)大部分年代春玉米高產(chǎn)高效適宜灌溉量為40 mm,且灌溉時間對春玉米產(chǎn)量和水分利用效率的影響較?。坏谖?、七和九氣候區(qū)大部分年代高產(chǎn)高效適宜灌溉量為60—80 mm,3個氣候區(qū)適宜灌溉時間分別為吐絲到吐絲后20 d、拔節(jié)到拔節(jié)后10 d、拔節(jié)到拔節(jié)后10 d。(3)與雨養(yǎng)條件相比,不同氣候區(qū)適宜灌溉措施條件下增產(chǎn)幅度不同。其中第五、七和九氣候區(qū)增產(chǎn)幅度較大,年代際變化范圍為33%—86%、24%—46%和50%—77%,第一和三氣候區(qū)增產(chǎn)幅度較小,年代際變化范圍為5%—43%和9%—19%。東北三省西部地區(qū)春玉米適宜灌溉量隨緯度的升高呈減少的趨勢,適宜灌溉時間隨緯度的升高呈推遲趨勢,且隨年代的推移,氣候變暖,各氣候區(qū)適宜灌溉時間呈提前趨勢。與雨養(yǎng)條件相比,各氣候區(qū)適宜灌溉措施條件下春玉米可增產(chǎn)0—86%,其中第五、七和九氣候區(qū)增產(chǎn)幅度較第一和第三氣候區(qū)更大。

    春玉米;灌溉措施;APSIM-Maize模型;水分利用效率;干旱;氣候區(qū)

    0 引言

    【研究意義】東北三省是我國重要的商品糧生產(chǎn)基地之一,根據(jù)中國農(nóng)業(yè)統(tǒng)計年鑒數(shù)據(jù),近10年(2008—2017年)東北三省玉米總產(chǎn)量約占東北三省糧食總產(chǎn)量的60.8%,約占全國玉米總產(chǎn)量的29.6%[1]。由此可見,玉米作為東北三省第一大糧食作物[2],其產(chǎn)量的提升對于全國糧食產(chǎn)量提升有重要意義。東北三省由于緯度較高,使其成為對氣候變化影響最為敏感的地區(qū)之一[3]。已有研究表明,近50年來(1960—2010年)東北三省平均、最低及最高氣溫均呈明顯上升趨勢,平均氣溫每10年升高0.38℃[4-5],使得玉米生長季內(nèi)積溫明顯增加;全年及玉米生長季內(nèi)降水量均呈減少趨勢且波動性增大[4,6-8],使該地區(qū)玉米生產(chǎn)的穩(wěn)定性受到極大威脅[9-10],特別是西部水資源匱乏地區(qū)[11]。因此研究氣候變化背景下,東北三省西部地區(qū)春玉米適宜灌溉措施,對于充分合理利用當(dāng)?shù)貧夂蛸Y源,應(yīng)對氣候變化,保障我國玉米高產(chǎn)穩(wěn)產(chǎn)具有重要理論和現(xiàn)實意義。【前人研究進(jìn)展】在不考慮適應(yīng)措施的前提下,玉米生長發(fā)育階段氣溫升高將使玉米生育期縮短,不利于玉米干物質(zhì)積累和產(chǎn)量形成。LIU等[12]利用APSIM模型解析了氣候變化對東北三省春玉米生長發(fā)育及產(chǎn)量形成的影響,研究結(jié)果表明,在不考慮適應(yīng)措施條件下,氣候變暖使東北三省春玉米全生育期每10年縮短2.3—4.8 d,生長季內(nèi)平均氣溫每升高1℃,產(chǎn)量將下降4%—22%。然而,在氣候變化背景下可以通過采取適當(dāng)?shù)倪m應(yīng)措施有效減緩氣候變化對作物的負(fù)面影響,如調(diào)整種植制度、更換品種、優(yōu)化播期、改善水肥管理等。已有研究表明,氣候變暖使得多熟種植北界北移,界線變化區(qū)域內(nèi)調(diào)整種植制度可使作物單產(chǎn)增加25%—106%[13],從而使我國三大糧食作物總產(chǎn)增加2.2%[14]。在東北地區(qū)更換生育期較長的玉米品種可使產(chǎn)量增加13%—38%[12]。調(diào)整播期可使華北地區(qū)夏玉米產(chǎn)量提升2%—10%[15],東北地區(qū)春玉米產(chǎn)量提升4%[12]?!颈狙芯壳腥朦c】優(yōu)化水肥管理措施亦可有效應(yīng)對氣候變化對作物生產(chǎn)產(chǎn)生的負(fù)面影響[16-18]。由于氣候、環(huán)境等的空間差異性和時間變異性,不同區(qū)域適宜的高產(chǎn)管理措施存在較大的差異。我國東北地區(qū)春玉米生長季內(nèi)降水量呈現(xiàn)經(jīng)向分布即由西北向東南逐漸遞增,空間差異較大(變化范圍為298—880 mm),且目前東北地區(qū)春玉米種植以雨養(yǎng)為主,如不采取相應(yīng)措施,干旱造成的玉米減產(chǎn)程度增加且各地減產(chǎn)幅度變化復(fù)雜,特別是對于西部降水偏少區(qū)域,水分是限制該地區(qū)春玉米產(chǎn)量的主要因子[19]。為此針對東北西部地區(qū)不同氣候區(qū)開展應(yīng)對氣候變化的適宜灌溉措施研究對保障該地區(qū)糧食高產(chǎn)穩(wěn)產(chǎn)有重要意義。【擬解決的關(guān)鍵問題】本文將東北三省劃分為10個氣候區(qū),選擇西部地區(qū)5個水分虧缺較為明顯的氣候區(qū)為研究區(qū)域,利用驗證后的農(nóng)業(yè)生產(chǎn)系統(tǒng)模型(agricultural production system simulator,APSIM-Maize),明確不同氣候區(qū)降水對春玉米產(chǎn)量的限制程度,結(jié)合玉米產(chǎn)量和水分利用效率提出不同氣候區(qū)應(yīng)對氣候變化的高產(chǎn)高效適宜灌溉模式,為東北三省西部地區(qū)春玉米高產(chǎn)高效生產(chǎn)提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域

    本文研究區(qū)域為我國東北三省,地理位置介于北緯38°26′—53°24′,東經(jīng)118°30′—135°8′,是我國緯度最高的地區(qū),屬于溫帶濕潤半濕潤大陸性季風(fēng)氣候,大部分農(nóng)區(qū)≥10℃積溫為1 600—3 600℃·d,無霜期140—170 d,降水量500—800 mm,干燥度0.75—1.5,80%的降水集中在5—9月份,與作物生長季吻合度較好。由于熱量資源略顯不足,東北大部分地區(qū)的種植模式為一年一熟。該地區(qū)土質(zhì)以黑土為主,是世界三大黑土帶之一,土地肥沃,有機(jī)質(zhì)含量高,為春玉米的種植提供了良好的條件。

    東北三省部分地區(qū)由于熱量資源的限制,不能種植春玉米,本文將穩(wěn)定通過10℃界限溫度范圍內(nèi)≥10℃的活動積溫達(dá)到2 100℃·d的地區(qū)界定為春玉米的潛在種植區(qū)。有研究表明東北三省西部地區(qū)降水對春玉米生產(chǎn)的限制達(dá)到14%—48%[20],因此筆者基于前人研究,依據(jù)積溫和水分虧缺率將研究區(qū)域劃分為10個氣候區(qū)(climate zone,CZ),每200℃·d劃分一個積溫帶[18],選擇水分虧缺較為明顯(>0)且位于東北三省西部的5個氣候區(qū)(CZ1、CZ3、CZ5、CZ7和CZ9)作為本文的研究區(qū)域。并從5個氣候區(qū)中分別選擇1個代表站點進(jìn)行研究(圖1),代表站點春玉米生長季內(nèi)的氣候資源特征見表1??梢钥闯?,各代表站點1981—2017年有效積溫呈上升趨勢,降水量均呈下降趨勢。

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

    氣象資料來源于中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://cdc.cma.gov.cn),包括研究區(qū)域5個站點1981—2017年的逐日地面氣象觀測數(shù)據(jù),站點分布如圖1所示,氣象要素有平均本站氣壓、平均氣溫、最高氣溫、最低氣溫、降水量、平均相對濕度、日照時數(shù)和平均風(fēng)速;模型所需的太陽輻射采用Penman- Monteith 公式[21]計算。

    作物資料來源于東北三省農(nóng)業(yè)氣象觀測站(1981—2007年)試驗數(shù)據(jù),包括品種、播種日期、播種密度、播種深度、行距、施肥和灌溉措施、生育期(出苗期、開花期和成熟期等)、地上部生物量和產(chǎn)量,用于作物模型品種參數(shù)的校準(zhǔn)。

    土壤資料來自于農(nóng)業(yè)氣象觀測站土壤數(shù)據(jù),包括分層土壤容重(BD)、飽和含水量(SAT)、田間持水量(DUL)、凋萎系數(shù)(LL15)等。

    表1 各氣候區(qū)代表站點春玉米生長季內(nèi)氣候資源特征(1981—2017年)

    **表示通過<0.01的顯著性檢驗 ** indicates significantly different at<0.01

    審圖號:GS(2020)5270號

    1.3 研究方法

    1.3.1 農(nóng)業(yè)生產(chǎn)系統(tǒng)模型 農(nóng)業(yè)生產(chǎn)系統(tǒng)模型(APSIM)是一種可用于模擬農(nóng)業(yè)生產(chǎn)系統(tǒng)中各主要組成部分的機(jī)理模型,它是由澳大利亞的聯(lián)邦科工組織以及昆士蘭州農(nóng)業(yè)生產(chǎn)系統(tǒng)組(APSRU)聯(lián)合開發(fā)的作物模型[22-23],目前已在全世界不同國家的地區(qū)得到廣泛驗證[24-29]。

    APSIM模型主要由生物物理模塊、管理模塊、輸入輸出模塊3部分組成,并由中心引擎來控制引導(dǎo),生物物理模塊用于模擬農(nóng)業(yè)系統(tǒng)中各種生物和物理過程,用戶可以利用管理模塊定義栽培管理措施并控制模擬過程,輸入輸出模塊管理數(shù)據(jù)的調(diào)用、輸入和輸出[30]。

    本文使用APSIM7.6版本,應(yīng)用的主要模塊為氣象模塊(Met)、時間模塊(Clock)、土壤模塊(Soil)、地表有機(jī)質(zhì)模塊(Surface organic matter)、施肥模塊(Fertiliser)、玉米模塊(Maize)、管理措施模塊(Manager folder)、灌溉模塊(Irrigation)和結(jié)果輸出模塊(Outputfile)。

    1.3.2 模型評價方法 本文利用實測的氣象、土壤及作物資料對APSIM-Maize模型進(jìn)行調(diào)參和驗證,通過模型模擬結(jié)果與實測結(jié)果的圖形比較以及各項評價指標(biāo)來評價APSIM-Maize模型在東北三省的適用性。采用以下統(tǒng)計量作為檢驗?zāi)P偷闹笜?biāo):模擬值與實測值之間的決定系數(shù)(R)、均方根誤差()、歸一化均方根誤差()[31]、一致性指標(biāo)(指標(biāo))[32]、平均絕對誤差()。其中,R和指標(biāo)可反映模擬值與實測值之間的一致性,其值愈接近1說明模擬效果愈好;相對于R,指標(biāo)對系統(tǒng)模擬誤差的響應(yīng)更敏感;和值反映了模擬值與實測值之間的絕對誤差和相對誤差,其值愈小,表明誤差值愈小。一般認(rèn)為值小于10%時,模擬值與實測值一致性非常好,10%—20%時較好[33]。

    式中,S為模擬值,O為實測值,為實測平均值,為樣本數(shù)。

    1.3.3 模型灌溉情景設(shè)置 本文利用調(diào)參驗證后的APSIM-Maize模型,模擬分析東北三省西部干旱區(qū)各氣候區(qū)典型站點春玉米產(chǎn)量對灌溉措施(灌溉時間和灌溉量)的響應(yīng),情景設(shè)置如下:

    (1)充分灌溉和雨養(yǎng)條件設(shè)置

    為明確各氣候區(qū)典型站點降水對春玉米產(chǎn)量的限制,需要模擬各氣候區(qū)代表性品種的潛在產(chǎn)量和氣候生產(chǎn)潛力。潛在產(chǎn)量可以代表一個地區(qū)充分灌溉條件下所能達(dá)到的最大產(chǎn)量,氣候生產(chǎn)潛力可以代表一個地區(qū)雨養(yǎng)條件下的產(chǎn)量。本文模擬潛在產(chǎn)量時設(shè)定春玉米生長過程中水肥充足,而模擬氣候生產(chǎn)潛力時設(shè)定養(yǎng)分充足但無灌溉。

    (2)灌溉情景的設(shè)置

    灌溉時間設(shè)定從出苗開始至玉米成熟,每隔10 d設(shè)定一個灌溉時間,為方便理解,結(jié)合玉米關(guān)鍵生育期(出苗、開花和成熟)對灌溉時期進(jìn)行命名,如E10代表灌溉時期為出苗后10 d,F(xiàn)10代表灌溉時期為開花后10 d,具體可參見圖2。每個灌溉時期分別設(shè)置5個灌溉量,即40、50、60、70和80 mm,分別模擬每年不同灌溉時間不同灌溉量的產(chǎn)量,然后按照1980s(1981—1990年)、1990s(1991—2000年)、2000s(2001—2010年)和2011—2017年4個年代分別統(tǒng)計相應(yīng)結(jié)果。各站點選擇已通過調(diào)參驗證的當(dāng)?shù)卦耘嗥贩N,基于農(nóng)業(yè)氣象觀測站春玉米實際栽培管理措施資料,設(shè)定播種深度為8 cm,行距為0.7 m,播種密度為50 000 株/hm2,播期為當(dāng)?shù)仄骄テ凇?/p>

    1.3.4 降水對產(chǎn)量的限制程度的計算 降水對產(chǎn)量的限制程度是指在雨養(yǎng)條件下作物產(chǎn)量相對于潛在產(chǎn)量的損失程度。

    式中,RC為降水對產(chǎn)量的限制程度(%);Y為潛在條件下的產(chǎn)量(kg·hm-2);Y為雨養(yǎng)條件下的產(chǎn)量(kg·hm-2)。

    1.3.5水分利用效率(water use efficiency,WUE) 水分利用效率是指消耗單位體積降水量和灌溉水所制造的干物質(zhì)量,計算公式如下:

    式中,WUE為水分利用效率(kg·hm-2·mm-1);YI為灌溉后的春玉米產(chǎn)量(kg·hm-2);P為生育期降水量(mm);Ir為灌溉量(mm)。

    2 結(jié)果

    2.1 APSIM-Maize模型有效性驗證

    本文代表性站點的作物參數(shù)校驗采用的是每個區(qū)一個品種,年代間品種不更替,即所有年代使用一個代表性品種,以農(nóng)業(yè)氣象觀測資料數(shù)據(jù)為基礎(chǔ),采用“試錯法”對各氣候區(qū)春玉米代表性品種的參數(shù)加以調(diào)試,目的是使模擬值與實測值之差盡可能小,最終確定各站點代表性品種的參數(shù)。代表性玉米品種的選擇基于農(nóng)業(yè)氣象站數(shù)據(jù),每個代表站點選擇一個品種,選取數(shù)據(jù)資料完整的6年,其中3年數(shù)據(jù)用于調(diào)參,3年數(shù)據(jù)用于結(jié)果驗證。生育期主要調(diào)試的參數(shù)為出苗到營養(yǎng)生長期結(jié)束、開花到灌漿和開花到成熟的熱時數(shù)(thermal time),產(chǎn)量主要調(diào)試每穗最大籽粒數(shù)和潛在灌漿速率,根據(jù)模型評價指標(biāo)判斷使用調(diào)試后的品種參數(shù)模擬值與實測值的擬合效果。圖3和圖4分別為各氣候區(qū)代表品種模擬與實測生育期天數(shù)、產(chǎn)量的比較結(jié)果。表2為各氣候區(qū)品種驗證結(jié)果評價。結(jié)果表明,播種—開花和播種—成熟生育階段模擬天數(shù)與實測天數(shù)基本吻合,決定系數(shù)(2)依次為0.72—0.92和0.80—0.89,值的范圍依次為0.89—0.96和0.77— 0.96,歸一化均方根誤差()在3%以下??芍P湍鼙容^準(zhǔn)確地模擬東北地區(qū)不同氣候區(qū)的春玉米生長發(fā)育進(jìn)程。

    模擬產(chǎn)量和實測產(chǎn)量的2值范圍為0.60—0.95,值范圍為0.81—0.88,歸一化均方根誤差()在16%以下,可以看出全區(qū)春玉米實測產(chǎn)量與模擬產(chǎn)量有很好的一致性,模型能比較準(zhǔn)確地模擬東北三省各氣候區(qū)的春玉米產(chǎn)量的形成。

    A:播種到開花,B:播種到成熟 A: from sowing to flowering, B: from sowing to maturity

    圖4 東北三省西部地區(qū)各氣候區(qū)春玉米產(chǎn)量實測值和模擬值驗證結(jié)果

    2.2 降水對研究區(qū)域內(nèi)代表站點春玉米產(chǎn)量的限制程度

    各氣候區(qū)全生育階段總降水量基本能滿足作物對水分的需求(第五氣候區(qū)除外),但是由于降水量年際波動較大且在玉米生長季內(nèi)分布不均,導(dǎo)致不同氣候區(qū)玉米不同生育階段出現(xiàn)不同程度的缺水。5個氣候區(qū)玉米出苗—拔節(jié)期和抽雄—成熟期大部分年份降水量可以滿足玉米需水,而在玉米水分關(guān)鍵期(拔節(jié)—抽雄期)各氣候區(qū)降水與需水的匹配程度不一,其中第一和第三氣候區(qū)降水量與需水量的差距不明顯,而第五、七和九氣候區(qū)降水量不能完全滿足作物需水(圖5)。

    A:播種到拔節(jié),B:拔節(jié)到抽雄,C:抽雄到成熟,D:全生育期

    基于各氣候區(qū)調(diào)參驗證后的APSIM-Maize模型,分別模擬1981—2017年春玉米潛在產(chǎn)量(充分灌溉)和雨養(yǎng)產(chǎn)量(無灌溉),采用公式(5)計算不同氣候區(qū)降水對玉米產(chǎn)量的限制程度,結(jié)果如圖6所示。結(jié)果表明,從過去37年平均來看,第一、三、五、七、九氣候區(qū)春玉米潛在產(chǎn)量分別為9 258、7 057、8 546、8 676和7 940 kg·hm-2,年代際之間差異較小。雨養(yǎng)產(chǎn)量分別為8 205、6 789、5 022、6 223和4 566 kg·hm-2,但年代際之間差異較大,變化范圍為 7 103—9 549、6 000—7 730、2 439—6 656、5 730—6 738和3 881— 5 659 kg·hm-2。其中第五、第七和第九氣候區(qū)雨養(yǎng)產(chǎn)量變異系數(shù)明顯高于第一和三氣候區(qū),這可能是由于第一和三氣候區(qū)在作物水分關(guān)鍵期降水基本能滿足作物需水,而第五、第七和第九氣候區(qū)在作物水分關(guān)鍵期降水并不能滿足作物需水(圖5)。

    5個氣候區(qū)降水對產(chǎn)量的限制程度分別為0— 27%、0—9%、27%—69%、15%—35%和31%—51%。因此可以看出第一和第三氣候區(qū)降水對產(chǎn)量的限制程度較小,通過灌溉對產(chǎn)量提升的貢獻(xiàn)較小,但能有效提高產(chǎn)量的穩(wěn)定性(第一氣候區(qū)產(chǎn)量變異系數(shù)由0.24降低到0.11,第三氣候區(qū)產(chǎn)量變異系數(shù)由0.14降低到0.12),表明在這2個氣候區(qū)采取適當(dāng)?shù)墓喔却胧┠苓_(dá)到穩(wěn)產(chǎn)目的。第五、七和九氣候區(qū)無論在哪個年代,春玉米的潛在產(chǎn)量與雨養(yǎng)產(chǎn)量差異明顯,說明降水對春玉米產(chǎn)量的限制較大,同時灌溉使得3個氣候區(qū)的產(chǎn)量變異系數(shù)由0.54降低到0.15,0.46降低到0.13,0.65降低到0.13。表明在東北三省西部干旱區(qū)通過灌溉能達(dá)到高產(chǎn)穩(wěn)產(chǎn)目的。

    表2 研究區(qū)域5個氣候區(qū)APSIM-Maize模型驗證結(jié)果評價

    圖6 東北三省西部地區(qū)各氣候區(qū)春玉米潛在產(chǎn)量、雨養(yǎng)產(chǎn)量及水分對春玉米產(chǎn)量的限制程度

    2.3 灌溉對春玉米產(chǎn)量和水分利用效率的影響

    利用調(diào)參驗證后的APSIM-Maize模型,分別模擬不同灌溉量(圖7)和灌溉時間(圖8)條件下各氣候區(qū)的春玉米產(chǎn)量和水分利用效率。

    總體而言,隨著灌溉量的增加,各氣候區(qū)春玉米產(chǎn)量呈現(xiàn)增長趨勢,但由于不同氣候區(qū)降水對春玉米產(chǎn)量的限制程度的差異,增加相同灌溉量對玉米產(chǎn)量的提升幅度不同。其中,第一和第三氣候區(qū)產(chǎn)量增加幅度較小,因此玉米的水分利用效率隨灌溉量的增加而減小。第五、七和九氣候區(qū)玉米的水分利用效率總體變化趨勢為先增大后減小。不僅不同氣候區(qū)之間存在差異,同一氣候區(qū)不同年代灌溉量對產(chǎn)量和水分利用效率的影響也存在差異(圖7)。分析不同灌溉時間對春玉米的產(chǎn)量和水分利用效率的影響發(fā)現(xiàn),第一和第三氣候區(qū)灌溉時間對春玉米產(chǎn)量和水分利用效率的影響較小,第五、七和九氣候區(qū)影響較大(圖8)。造成氣候區(qū)間差異的原因主要是由于第一和第三氣候區(qū)在作物水分關(guān)鍵期降水基本能滿足作物需水,而第五、第七和第九氣候區(qū)在作物水分關(guān)鍵期降水并不能滿足作物需水。

    圖7 灌溉量對東北三省西部地區(qū)各氣候區(qū)春玉米產(chǎn)量和水分利用效率的影響

    圖8 灌溉時期對東北三省西部地區(qū)各氣候區(qū)春玉米產(chǎn)量和水分利用效率的影響

    不同氣候區(qū)不同年代高產(chǎn)高效的春玉米適宜灌溉措施如圖7和圖8所示,具體而言,結(jié)合圖5可以看出,在1980s,第一和第三氣候區(qū)各生育階段降水量與生育期內(nèi)需水量匹配較好,適宜的灌溉量為40 mm,且灌溉時期影響較小;第五氣候區(qū)適宜的灌溉量為60 mm,適宜的灌溉時期在吐絲后20 d;第七氣候區(qū)適宜的灌溉量為60 mm,適宜的灌溉時期在拔節(jié)期;第九氣候區(qū)適宜的灌溉量為60 mm,適宜的灌溉時期在拔節(jié)后10 d。在1990s,第一氣候區(qū)適宜的灌溉量為60 mm,灌溉時期影響較??;第三氣候區(qū)適宜的灌溉量為40 mm,灌溉時期的影響較小;第五氣候區(qū)適宜的灌溉量為60 mm,適宜的灌溉時期在吐絲后20 d;第七氣候區(qū)適宜的灌溉量為40 mm,適宜的灌溉時期在拔節(jié)期;第九氣候區(qū)適宜的灌溉量為60 mm,適宜的灌溉時期在拔節(jié)后10 d。在2000s,第一氣候區(qū)適宜的灌溉量為60 mm,灌溉時期影響較?。坏谌龤夂騾^(qū)適宜的灌溉量為40 mm,灌溉時期的影響較小;第五氣候區(qū)適宜的灌溉量為80 mm,適宜的灌溉時期在吐絲期;第七氣候區(qū)適宜的灌溉量為60 mm,適宜的灌溉時期在拔節(jié)后10 d;第九氣候區(qū)適宜的灌溉量為80 mm,適宜的灌溉時期在拔節(jié)期。在2011—2017年,第一氣候區(qū)適宜的灌溉量為40 mm,灌溉時期影響較??;第三氣候區(qū)適宜的灌溉量為40 mm,灌溉時期的影響較??;第五氣候區(qū)適宜的灌溉量為50 mm,適宜的灌溉時期在拔節(jié)后10 d;第七氣候區(qū)適宜的灌溉量為60 mm,適宜的灌溉時期在拔節(jié)期;第九氣候區(qū)適宜的灌溉量為60 mm,適宜的灌溉時期在拔節(jié)期(圖5,7—8)。

    綜上,2000s適宜灌溉量普遍高于其他年代,第一和三氣候區(qū)受灌溉時間影響較小,第五氣候區(qū)適宜灌溉時間在吐絲期前后,第七和九氣候區(qū)適宜灌溉時間在拔節(jié)期前后,可看出隨緯度的增大,適宜灌溉時間呈推遲趨勢,且隨年代的推移,各氣候區(qū)適宜灌溉時間呈提前趨勢。

    2.4 不同灌溉方案較雨養(yǎng)條件的增產(chǎn)幅度

    將適宜灌溉措施和雨養(yǎng)2種條件下對應(yīng)的模擬產(chǎn)量按氣候區(qū)進(jìn)行匯總分析,得出2種條件下各氣候區(qū)各年代的增產(chǎn)幅度(圖9)。

    審圖號:GS(2020)5270號

    與雨養(yǎng)條件相比,各區(qū)適宜灌溉措施增產(chǎn)幅度如下:第一氣候區(qū)1980s、1990s、2000s和2011—2017這4個研究時段增產(chǎn)幅度分別為5%、13%、43%和16%;第三氣候區(qū)分別為9%、19%、15%和9%;第五氣候區(qū)分別為42%、48%、86%和33%;第七氣候區(qū)分別為46%、24%、34%和39%;第九氣候區(qū)分別為77%、50%、76%和56%??梢姴扇∵m宜的灌溉措施較雨養(yǎng)條件下可實現(xiàn)不同程度的增產(chǎn),但不同氣候區(qū)之間增產(chǎn)效果存在差異。第一和第三氣候區(qū)降水條件較好,作物需水關(guān)鍵期受水分脅迫較小,因此適宜灌溉措施條件下增產(chǎn)幅度較小。第五、七和九氣候區(qū)降水條件較差,作物生長季受水分脅迫較嚴(yán)重,因此適宜灌溉措施條件下增產(chǎn)幅度較大。不同年代適宜灌溉措施較雨養(yǎng)條件的增產(chǎn)幅度也存在差異,從圖5中可看出,2000s降水條件較差,降水集中時間與作物水分關(guān)鍵期匹配度較差,因此相較于其他3個時段,2000s適宜灌溉措施條件下增產(chǎn)幅度較大。

    3 討論

    前人采用田間試驗、統(tǒng)計分析及作物模型的方法開展了作物適宜灌溉措施研究[34-38],其中,田間試驗受限于人力物力、氣候條件等因素,無法全面反映灌溉措施在不同區(qū)域不同年代對作物的影響;統(tǒng)計分析主要是依據(jù)前人研究和實際數(shù)據(jù)進(jìn)行數(shù)理統(tǒng)計分析,缺少生理依據(jù);而作物模型方法很好地解決了田間試驗代表性差和統(tǒng)計分析缺乏生理基礎(chǔ)的問題。

    本文使用調(diào)參驗證后的APSIM-Maize模型,模擬了東北三省西部不同氣候區(qū)不同灌溉情境下春玉米的產(chǎn)量,綜合考慮產(chǎn)量和水分利用效率2個因素,提出了各氣候區(qū)的高產(chǎn)高效適宜灌溉措施,可為東北三省西部春玉米應(yīng)對氣候變化提供科學(xué)依據(jù)。

    基于玉米產(chǎn)量和水分利用效率確定灌溉量時,2000s適宜灌溉量普遍高于其他年代,這主要是因為該年代降水較少且并沒有集中在玉米的需水關(guān)鍵期。選擇適宜灌溉時間時,第一和第三氣候區(qū)受灌溉時間影響較小,第五氣候區(qū)適宜灌溉時間在吐絲期前后,第七、九氣候區(qū)適宜灌溉時間在拔節(jié)期前后,可看出隨緯度的增大,適宜灌溉時間呈推遲趨勢,且隨年代的推移,各氣候區(qū)適宜灌溉時間呈提前趨勢,這主要是因為喜溫作物生長季內(nèi)東北三省呈暖干趨勢[12],氣溫升高會使玉米生育期縮短,水分關(guān)鍵期提前,從而加重玉米春旱。

    東北三省西部5個氣候區(qū)水分虧缺率均大于0,表明5個氣候區(qū)均存在不同程度的水分虧缺情況,然而不同氣候區(qū)水分虧缺程度存在較大差異,因此適宜的高產(chǎn)高效灌溉方案存在時空差異性。同時東部地區(qū)水分條件較好,大部分年份玉米關(guān)鍵生育期降水可滿足玉米水分的需求,高產(chǎn)高效的灌溉方案應(yīng)與西部地區(qū)存在較大差異,例如宮亮等[39]在遼寧省瓦房店采用田間試驗方法研究發(fā)現(xiàn),不同灌溉時間對玉米生長發(fā)育及產(chǎn)量有不同影響,研究結(jié)果表明該試驗區(qū)最佳灌溉時間為生育后期,抽雄期和灌漿期各補(bǔ)充灌溉一次為最佳。因此在分析適宜灌溉措施時應(yīng)考慮氣候背景的差異性。本文采用模型模擬作物產(chǎn)量時選用各氣候區(qū)代表品種進(jìn)行研究,由于資料有限,并未能針對不同品種提出適宜的灌溉措施,然而不同品種不同生育期對水分的需求差異較大,未來應(yīng)針對不同品種開展相應(yīng)研究。同時,本研究發(fā)現(xiàn),在東北三省西部適宜灌溉時間隨緯度增大而呈推遲趨勢,該趨勢是否適用于其他地區(qū)及灌溉與緯度之間是否存在相關(guān)性,仍需進(jìn)一步驗證。此外,未來研究中還應(yīng)根據(jù)不同降水年型提出適宜的灌溉措施。

    4 結(jié)論

    本文以東北三省春玉米為研究對象,基于1981— 2017年的氣候資料、1981—2007年春玉米的作物資料和土壤資料對APSIM-Maize模型進(jìn)行調(diào)參和驗證,通過模型模擬結(jié)果與實測結(jié)果的圖形比較以及各項評價指標(biāo)評價了APSIM-Maize模型在東北地區(qū)的適用性,表明該模型可用于東北三省西部春玉米產(chǎn)量的模擬研究。

    通過設(shè)置不同的灌溉情景,基于驗證后的APSIM-Maize模型,模擬不同灌溉方案對春玉米產(chǎn)量的影響,確定了研究區(qū)域內(nèi)不同氣候區(qū)適宜的高產(chǎn)高效灌溉模式。主要結(jié)論如下:第一氣候區(qū)1980s、1990s、2000s和2011—2017年這4個研究時段內(nèi)適宜的灌溉量分別為40 mm、60 mm、60 mm、40 mm;第三氣候區(qū)4個研究時段內(nèi)適宜的灌溉量為40 mm、40 mm、40 mm、40 mm;第五氣候區(qū)4個研究時段內(nèi)適宜的灌溉量為60 mm、60 mm、80 mm、50 mm;第七氣候區(qū)4個研究時段內(nèi)適宜的灌溉量為60 mm、40 mm、60 mm、60 mm;第九氣候區(qū)4個研究時段內(nèi)適宜的灌溉量為60 mm、60 mm、80 mm、60 mm。第一和三氣候區(qū)灌溉時期對玉米產(chǎn)量和水分利用效率影響較小,第五氣候區(qū)4個研究時段適宜的灌溉時期分別為吐絲后20 d、吐絲后20 d、吐絲期、拔節(jié)后10 d;第七氣候區(qū)4個研究時段適宜的灌溉時期分別為拔節(jié)期、拔節(jié)期、拔節(jié)后10 d、拔節(jié)期;第九氣候區(qū)4個研究時段適宜的灌溉時期分別為拔節(jié)后10 d、拔節(jié)后10 d、拔節(jié)期、拔節(jié)期。各氣候區(qū)適宜灌溉措施條件下獲得的產(chǎn)量較無灌溉條件下產(chǎn)量均有不同程度提高,增產(chǎn)幅度為0—86%,其中第五、七和九氣候區(qū)增產(chǎn)幅度較第一和第三氣候區(qū)更大。

    [1] 中華人民共和國國家統(tǒng)計局. 中國統(tǒng)計年鑒. 北京: 中國統(tǒng)計出版社, 2009-2018.

    National Bureau of Statistics of the People’s Republic of China.. Beijing: China Statistics Press, 2009-2018. (in Chinese)

    [2] 陳群, 耿婷, 侯雯嘉, 陳長青. 近20年東北氣候變暖對春玉米生長發(fā)育及產(chǎn)量的影響. 中國農(nóng)業(yè)科學(xué), 2014, 47(10): 1904-1916.

    CHEN Q, GENG T, HOU W J, CHEN C Q. Impacts of climate warming on growth and yield of spring maize in recent 20 years in Northeast China., 2014, 47(10): 1904-1916. (in Chinese)

    [3] 趙秀蘭. 近50年中國東北地區(qū)氣候變化對農(nóng)業(yè)的影響. 東北農(nóng)業(yè)大學(xué)學(xué)報, 2010, 41(9): 144-149.

    ZHAO X L. Influence of climate change on agriculture in Northeast China in recent 50 years.2010, 41(9): 144-149. (in Chinese)

    [4] 劉志娟, 楊曉光, 王文峰, 李克南, 張曉煜. 氣候變化背景下我國東北三省農(nóng)業(yè)氣候資源變化特征. 應(yīng)用生態(tài)學(xué)報. 2009, 20(9): 2199-2206.

    LIU Z J, YANG X G, WANG W F, LI K N, ZHANG X Y. Characteristics of agricultural climate resources in three provinces of Northeast China under global climate change., 2009, 20(9): 2199-2206. (in Chinese)

    [5] 董滿宇, 吳正方. 近50年來東北地區(qū)氣溫變化時空特征分析. 資源科學(xué), 2008, 30(7) :1093-1099.

    DONG M Y, WU Z F. Analysis of temporal and spatial characteristics of temperature change over the last 50 years in Northeast China.2008, 30(7) :1093-1099. (in Chinese)

    [6] 付長超, 劉吉平, 劉志明. 近60年東北地區(qū)氣候變化時空分異規(guī)律的研究. 干旱區(qū)資源與環(huán)境, 2009, 23(12): 60-65.

    FU C C, LIU J P, LIU Z M. Spacial and temporal differentiation rule of the climate change in Northeast China in about 60 years.2009, 23(12): 60-65. (in Chinese)

    [7] 姜曉艷, 劉樹華, 馬明敏, 張菁, 宋軍. 東北地區(qū)近百年降水時間序列變化規(guī)律的小波分析. 地理研究, 2009, 28(2): 354-362.

    JIANG X Y, LIU S H, MA M M, ZHANG J, SONG J. A wavelet analysis of the precipitation time series in Northeast China during the last 100 years.2009, 28(2): 354-362. (in Chinese)

    [8] 孫鳳華, 吳志堅, 楊素英. 東北地區(qū)近50年來極端降水和干燥事件時空演變特征. 生態(tài)學(xué)雜志, 2006, 25(7): 779-784.

    SUN F H, WU Z J, YANG S Y. Temporal and spatial variations of extreme precipitation and dryness events in Northeast China in last 50 years.2006, 25(7): 779-784. (in Chinese)

    [9] 顧莉麗, 郭慶海. 中國糧食主產(chǎn)區(qū)的演變與發(fā)展研究. 農(nóng)業(yè)經(jīng)濟(jì)問題, 2011, 32(8): 4-9, 110.

    GU L L, GUO Q H. Evolution and development of China’s major grain producing areas.2011, 32(8): 4-9, 110. (in Chinese)

    [10] 劉?;? 陳新平, 崔振嶺, 孟慶鋒, 趙明. 三大糧食作物產(chǎn)量潛力與產(chǎn)量差研究進(jìn)展. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2015, 23(5): 525-534.

    LIU B H, CHEN X P, CUI Z L, MENG Q F, ZHAO M. Research advance in yield potential and yield gap of three major cereal crops.2015, 23(5): 525-534. (in Chinese)

    [11] 安剛, 孫力, 廉毅. 東北地區(qū)可利用降水資源的初步分析. 氣候與環(huán)境研究, 2005, 10(1): 132-139.

    AN G, SUN L, LIAN Y. A primary analysis of utilizable precipitation in Northeast China.2005, 10(1): 132-139. (in Chinese)

    [12] LIU Z, HUBBARD K G, LIN X M , YANG X G. Negative effects of climate warming on maize yield are reversed by the changing of sowing date and cultivar selection in Northeast China., 2013, 19(11): 3481-3492.

    [13] 楊曉光, 劉志娟, 陳阜. 全球氣候變暖對中國種植制度可能影響Ⅰ.氣候變暖對中國種植制度北界和糧食產(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)

    [14] YANG X G, CHEN F, LIN X M, LIU Z J, ZHANG H L, ZHAO J, LI K N, YE Q, LI Y, LV S, YANG P, WU W B, LI Z G, LAL R, TANG H J. Potential benefits of climate change for crop productivity in China., 2015, 208: 76-84.

    [15] 張鎮(zhèn)濤, 楊曉光, 高繼卿, 王曉煜, 白帆, 孫爽, 劉志娟, 明博, 謝瑞芝, 王克如, 李少昆. 氣候變化背景下華北平原夏玉米適宜播期分析. 中國農(nóng)業(yè)科學(xué), 2018, 51(17): 3258-3274.

    ZHANG Z T, YANG X G, GAO J Q, WANG X Y, BAI F, SUN S, LIU Z J, MING B, XIE R Z, WANG K R, LI S K. Analysis of suitable sowing date for summer maize in north China plain under climate change.2018, 51(17): 3258-3274. (in Chinese)

    [16] SUN S, YANG X G, LIN X M, GRETCHEN F S, LI K N. Climate-smart management can further improve winter wheat yield in China., 2018, 162: 10-18.

    [17] 王曉煜, 楊曉光, Tao L, 張?zhí)煲? 劉濤, 項洪濤, 孫翌晨, 劉志娟. 東北三省水稻干濕交替灌溉模式適宜性分區(qū). 農(nóng)業(yè)工程學(xué)報, 2018, 34(6): 111-120.

    WANG X Y, YANG X G, TAO L, ZHANG T Y, LIU T, XIANG H T, SUN Y C, LIU Z J. Rice suitability zoning of alternative wetting and drying irrigation mode in three provinces of Northeast China., 2018, 34(6): 111-120. (in Chinese)

    [18] LIU Z J, YANG X G, LIN X M, HUBBARD K D, Lü S, WANG J. Narrowing the agronomic yield gaps of maize by improved soil, cultivar, and agricultural management practices in different climate zones of Northeast China., 2016, 20(12): 1-18.

    [19] WANG X H, PENG L Q, ZHANG X P, YIN G D. Divergence of climate impacts on maize yield in northeast china., 2014, 196 (7522): 51-58.

    [20] LIU Z J, YANG X G, LIN X M, GOWDA P, LV S, WANG J. Climate zones determine where substantial increases of maize yields can be attained in Northeast China., 2018, 149: 473-487.

    [21] ALLEN R G, PEREIRA L S, RAES D, SMITH M.. Rome: Food and Agriculture Organization of the United Nations, 1998: 56.

    [22] ASSENG S, KEATING B A, IRP F, GREGORY P J. Performance of the APSIM-wheat model in western Australia., 1998, 57(2): 163-179.

    [23] ASSENG S, KEULEN H V, STOL W. Performance and application of the APSIM Nwheat model in the Netherlands., 2000, 12(1): 37-54.

    [24] 王琳, 鄭有飛, 于強(qiáng), 王恩利. APSIM模型對華北平原小麥-玉米連作系統(tǒng)的適用性. 應(yīng)用生態(tài)學(xué)報, 2007, 18(11): 2480-2486.

    WANG L, ZHENG Y F, YU Q, WANG E L. Applicability of Agricultural Production Systems Simulator (APSIM) in simulating the production and water use of wheat-maize continuous cropping system in North China Plain.2007, 18(11): 2480-2486. (in Chinese)

    [25] 孫寧, 馮利平. 利用冬小麥作物生長模型對產(chǎn)量氣候風(fēng)險的評估. 農(nóng)業(yè)工程學(xué)報, 2005, 21(2): 106-110.

    SUN N, FENG L P. Assessing the climatic risk to crop yield of winter wheat using crop growth models., 2005, 21(2): 106-110. (in Chinese)

    [26] 李艷, 薛昌穎, 劉園, 楊曉光, 王靖, 王恩利. APSIM模型對冬小麥生長模擬的適應(yīng)性研究. 氣象, 2008, 34(增刊): 271-279.

    LI Y, XUE C Y, LIU Y, YANG X G, WANG J, WANG E L. Adaptability of APSIM model to winter wheat growth simulation.2008, 34(Supplement): 271-279. (in Chinese)

    [27] 劉志娟, 楊曉光, 王靜, 呂碩, 李克南, 荀欣, 王恩利. APSIM玉米模型在東北地區(qū)的適應(yīng)性. 作物學(xué)報, 2012, 38(4): 740-746.

    LIU Z J, YANG X G, WANG J, Lü S, LI K N, XUN X, WANG E L. Adaptability of APSIM Maize model in Northeast China.2012, 38(4): 740-746. (in Chinese)

    [28] 王靜, 楊曉光, 呂碩, 劉志娟, 李克南, 荀欣, 劉園, 王恩利. 黑龍江省春玉米產(chǎn)量潛力及產(chǎn)量差的時空分布特征. 中國農(nóng)業(yè)科學(xué), 2012, 45(10): 1914-1925.

    WANG J, YANG X G, Lü S, LIU Z J, LI K N, XUN X, LIU Y, WANG E L. Spatial-temporal characteristics of potential yields and yield gaps of spring maize in Heilongjiang province.2012, 45(10): 1914-1925. (in Chinese)

    [29] LOBELL D B, HAMMER G L, MCLEAN G. The critical role of extreme heat for maize production in the United States., 2013, 3(5): 497-501.

    [30] KEATING B A, CARBERRY P S, HAMMER G L, PROBERT M E, ROBERTSON M J, HOLZWORTH D, HUTH N I, HARGREAVES J N G, MEINKE H, HOCHMAN Z, MCLEAN G, VERBURG K, SNOW V, DIMES J P, SILBURN M, WANG E. An overview of APSIM, a model designed for farming systems simulation., 2003, 18(3): 267-288.

    [31] WALLACH D, GOFFINET B. Mean squared error of prediction in models for studying ecological and agronomic systems., 1987, 43(3): 561-573.

    [32] WILLMOTT C J. Some comments on the evaluation of model performance., 1982, 63(11): 1309-1313

    [33] 曹衛(wèi)星, 朱艷, 湯亮, 荊奇, 田永超, 姚霞, 劉小軍. 數(shù)字農(nóng)作技術(shù). 北京: 科學(xué)出版社, 2008.

    Cao W X, Zhu Y, Tang L, Jing Q, Tian Y C, Yao X, Liu X J.. Beijing: Science Press, 2008. (in Chinese)

    [34] 曹彩云, 黨紅凱, 鄭春蓮, 郭麗, 馬俊永, 李科江. 不同灌溉模式對小麥產(chǎn)量、耗水及水分利用效率的影響. 華北農(nóng)學(xué)報, 2016, 31(S1): 17-24.

    CAO C Y, DANG H K, ZHENG C L, GUO L, MA J Y, LI K J. Effects of different irrigation regime on yield, water consumption and water use efficiency of winter wheat.2016, 31(S1): 17-24. (in Chinese)

    [35] 金欣欣, 張喜英, 陳素英, 孫宏勇, 王彥梅, 邵立威, 高麗娜. 不同灌溉次數(shù)和灌溉量對冬小麥氮素吸收轉(zhuǎn)移的影響. 華北農(nóng)學(xué)報, 2009, 24(4): 112-118.

    JIN X X, ZHANG X Y, CHEN S Y, SUN H Y, WANG Y M, SHAO L W, GAO L N. Effect of different irrigation frequency and amount on nitrogen uptake, translocation of winter wheat.2009, 24(4): 112-118. (in Chinese)

    [36] 陳林, 王磊, 宋乃平, 張慶霞, 裴雪雁. 灌溉量和灌溉次數(shù)對紫花苜蓿耗水特性和生物量的影響. 水土保持學(xué)報, 2009, 23(4): 91-95.

    CHEN L, WANG L, SONG N P, ZHANG Q X, PEI X Y. Effects of irrigation amount and times on water consumption characteristics and biomass of alfalfa.2009, 23(4): 91-95. (in Chinese)

    [37] 郭金路, 尹光華, 谷健, 劉作新. 基于CROPWAT模型的阜新地區(qū)春玉米灌溉制度的確定. 生態(tài)學(xué)雜志, 2016, 35(12): 3428-3434.

    GUO J L, YIN G H, GU J, LIU Z X. Determination of irrigation scheduling of spring maize in different hydrological years in Fuxin, Liaoning province based on CROPWAT model.2016, 35(12): 3428-3434. (in Chinese)

    [38] 李良濤, 毛振強(qiáng), 牛靈安, 郝晉珉, 宇振榮. 基于作物生長模型的冬小麥灌溉方案研究. 灌溉排水學(xué)報, 2010, 29(2): 22-27.

    LI L T, MAO Z Q, NIU L A, HAO J M, YU Z R. Irrigation scenarios of winter wheat production aided by crop growth simulation model.2010, 29(2): 22-27. (in Chinese)

    [39] 宮亮, 孫文濤, 雋英華, 尹同波, 劉玉軍. 補(bǔ)充灌溉對玉米生理指標(biāo)及水分利用效率的影響. 節(jié)水灌溉, 2017(1): 9-11, 15.

    GONG L, SUN W T, TUAN Y H, YIN T B, LIU Y J. Effects of supplementary irrigation on physiological indices and water use efficiency of maize., 2017(1): 9-11, 15. (in Chinese)

    Analysis of Suitable Irrigation Schemes with High-Production and High-Efficiency for Spring Maize to Adapt to Climate Change in the West of Northeast China

    HUANG QiuWan, LIU ZhiJuan, YANG XiaoGuang, BAI Fan, LIU Tao, ZHANG ZhenTao, SUN Shuang, ZHAO Jin

    (College of Resources and Environmental Sciences, China Agricultural University, Beijing 100193)

    【】 The three provinces of Northeast China are not only an important grain production commodity in the country, but also the most sensitive areas to climate change. Thus, it is critical to clarify the suitable irrigation schemes for spring maize in arid areas of the west of Northeast China, which may be benefit on spring maize yield and its stability with higher water use efficiency under climate change. 【】 Based on the accumulated temperature and water deficit rate () during the growing season in 1981-2017, the potential planting areas of spring maize in Northeast China were divided into 10 climate zones (CZs). Five of them (>0) in the west portion were selected as research areas.With meteorological data, experimental data, and soil data, maize yield potential was assessed by the well-calibrated and validated agricultural production system model (APSIM-Maize) in each CZ under different irrigation scenarios. According to the comprehensive analysis of both yield and water use efficiency, the suitable irrigation measures and the yield increment in different decades in each CZ were identified. 【】 (1) In the past 37 years, the water limitation on spring maize yield in CZ1 and CZ3 were lessthan that in threeother CZs, with a range of 0-27% and 0-9%, respectively. Irrigation contributed little to yield increase, but it could improve the yield stability. The coefficient of yield variation reduced from 0.24 to 0.11 in CZ1, respectively, and reduced from 0.14 to 0.12 in CZ3, respectively. In CZ5, CZ7 and CZ9, more water limitation was found on maize yield, with a range of 27%-69%, 15%-35%, and 31%-51%,respectively. Moreover, irrigation also reduced the coefficients of yield variation by 0.39, 0.33 and 0.52 in the three CZs. The results indicated that irrigation could lead to a high and stable maize yield in the arid areas of Northeast China. (2) The suitable irrigation amount was 40 mm, which could produce high spring maize yield with high water use efficiency in CZ1 and CZ3. However, irrigation time had a little influence on the yield and water use efficiency of spring maize. Meanwhile, the suitable irrigation amounts for high yield and high water use efficiency were 60-80 mm in CZ5, CZ7 and CZ9, and the suitable irrigation times were from silking to 20 days after silking, jointing to 10 days after jointing and jointing to 10 days after jointing. (3)Compared with the rain-fed conditions, the yield increments varied in different CZs under suitable irrigation measures, which ranged from 33% to 86%, 24% to 46% and 50% to 77% in CZ5, CZ7 and CZ9, respectively. Lower yield increments were found in CZ1 and CZ3, with ranges of 5% to 43% and 9% to 19%, respectively. 【】 The suitable irrigation amount for spring maize decreased with the latitude increased, and the suitable irrigation time delayed with the increased latitude.In addition, the suitable irrigation time in each CZ advanced because of warming climate. Compared with the rain-fed conditions, spring maize yield could be increased by 0-86% under suitable irrigation measures in each CZ. In particular, the yield increments in CZ5, CZ7 and CZ9 were greater than those in CZ1 and CZ3.

    spring maize; irrigation measures; APSIM-Maize; water use efficiency; drought; climate zone

    10.3864/j.issn.0578-1752.2020.21.015

    2019-05-16;

    2019-06-25

    國家重點研發(fā)計劃(2016YFD0300101,2017YFD0300301)

    黃秋婉,E-mail:hqw893784946@163.com。通信作者劉志娟,E-mail:zhijuanliu@cau.edu.cn

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

    猜你喜歡
    氣候區(qū)東北三省降水
    2003 年—2022 年東北三省地區(qū)生產(chǎn)總值及一般公共預(yù)算收支情況表
    初中地理《東北三省》教學(xué)案例與反思
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    各氣候區(qū)被動建筑節(jié)能技術(shù)實用性分析
    科技資訊(2017年23期)2017-09-09 13:08:43
    冷涼氣候區(qū)不同形態(tài)氮肥對蘋果根系氮代謝的影響
    濕熱氣候區(qū)建筑防熱研究進(jìn)展
    降水現(xiàn)象儀模擬軟件設(shè)計與實現(xiàn)
    東北三省滑雪產(chǎn)業(yè)的空間拓展路徑
    冰雪運動(2016年5期)2016-04-16 05:55:30
    中國不同氣候區(qū)土壤濕度特征及其氣候響應(yīng)
    久久精品国产综合久久久| 一边摸一边做爽爽视频免费| 免费高清在线观看日韩| 激情视频va一区二区三区| 中文字幕色久视频| 叶爱在线成人免费视频播放| 97在线人人人人妻| 久久 成人 亚洲| 人人妻人人澡人人看| 高清视频免费观看一区二区| 男人舔女人的私密视频| 亚洲av片天天在线观看| 99热国产这里只有精品6| 纵有疾风起免费观看全集完整版| 成人免费观看视频高清| 久久国产精品影院| 久久热在线av| 国产成人一区二区三区免费视频网站| 国产精品自产拍在线观看55亚洲 | 女性生殖器流出的白浆| 精品亚洲成国产av| 精品久久久久久电影网| 亚洲熟妇熟女久久| 亚洲专区国产一区二区| 伦理电影免费视频| 十八禁人妻一区二区| 91麻豆精品激情在线观看国产 | 天堂动漫精品| 亚洲一码二码三码区别大吗| 亚洲成人国产一区在线观看| 国产一卡二卡三卡精品| kizo精华| 啦啦啦 在线观看视频| 国产真人三级小视频在线观看| 三上悠亚av全集在线观看| 少妇 在线观看| 大码成人一级视频| 亚洲欧美精品综合一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 99在线人妻在线中文字幕 | 咕卡用的链子| 在线亚洲精品国产二区图片欧美| 国产精品成人在线| 久久久久久人人人人人| av网站在线播放免费| 视频区欧美日本亚洲| 亚洲av成人一区二区三| 久久这里只有精品19| av欧美777| 国产成人免费无遮挡视频| av有码第一页| 欧美激情久久久久久爽电影 | 国产欧美日韩一区二区精品| 国产一区二区激情短视频| 自线自在国产av| 99国产精品一区二区三区| 欧美中文综合在线视频| 国产成人影院久久av| a在线观看视频网站| 亚洲七黄色美女视频| 国产精品美女特级片免费视频播放器 | 久久久久视频综合| av国产精品久久久久影院| 下体分泌物呈黄色| 久久久久久久大尺度免费视频| 欧美日韩亚洲高清精品| 一二三四在线观看免费中文在| 国产精品1区2区在线观看. | 又黄又粗又硬又大视频| 午夜两性在线视频| 天堂俺去俺来也www色官网| 欧美av亚洲av综合av国产av| 久久久精品免费免费高清| 欧美日本中文国产一区发布| 天堂中文最新版在线下载| 亚洲天堂av无毛| 亚洲综合色网址| 99香蕉大伊视频| 精品国产一区二区久久| 大片电影免费在线观看免费| 我的亚洲天堂| 久久青草综合色| xxxhd国产人妻xxx| 一级毛片女人18水好多| 国产成人av激情在线播放| 亚洲天堂av无毛| 超碰97精品在线观看| 天天躁日日躁夜夜躁夜夜| 一级毛片电影观看| 黄色片一级片一级黄色片| 国产成人精品久久二区二区91| 99精品欧美一区二区三区四区| 免费看a级黄色片| 侵犯人妻中文字幕一二三四区| 丝袜美足系列| 久久久久精品国产欧美久久久| 一进一出抽搐动态| 最近最新免费中文字幕在线| av有码第一页| 国产免费av片在线观看野外av| 午夜精品国产一区二区电影| 手机成人av网站| 精品福利观看| 下体分泌物呈黄色| 另类精品久久| 丁香六月天网| 国产精品1区2区在线观看. | 在线观看舔阴道视频| 露出奶头的视频| 十八禁网站免费在线| 国产又色又爽无遮挡免费看| 国产高清激情床上av| 国产精品久久久人人做人人爽| 亚洲国产毛片av蜜桃av| 丝瓜视频免费看黄片| 欧美 日韩 精品 国产| 另类精品久久| 亚洲精品中文字幕在线视频| 久久久久网色| 精品一区二区三区av网在线观看 | 国产高清视频在线播放一区| 免费一级毛片在线播放高清视频 | 一级,二级,三级黄色视频| 大香蕉久久网| 久久毛片免费看一区二区三区| av天堂在线播放| 妹子高潮喷水视频| 一边摸一边抽搐一进一出视频| 精品少妇一区二区三区视频日本电影| 亚洲精品国产色婷婷电影| 亚洲国产av影院在线观看| 国产成人影院久久av| 国产伦人伦偷精品视频| 19禁男女啪啪无遮挡网站| 少妇被粗大的猛进出69影院| 极品少妇高潮喷水抽搐| 不卡一级毛片| 国产精品98久久久久久宅男小说| netflix在线观看网站| 黄色怎么调成土黄色| 久久久久精品国产欧美久久久| bbb黄色大片| 免费不卡黄色视频| 欧美另类亚洲清纯唯美| 久久久久精品国产欧美久久久| 丰满迷人的少妇在线观看| 天堂8中文在线网| videosex国产| 99国产综合亚洲精品| 亚洲成人国产一区在线观看| 在线观看一区二区三区激情| av线在线观看网站| 女人高潮潮喷娇喘18禁视频| 国产高清视频在线播放一区| 亚洲成人免费av在线播放| 在线观看免费视频网站a站| 国产伦理片在线播放av一区| 一本一本久久a久久精品综合妖精| 一个人免费在线观看的高清视频| 少妇被粗大的猛进出69影院| 两个人免费观看高清视频| av视频免费观看在线观看| 99re在线观看精品视频| 中文字幕制服av| 中文亚洲av片在线观看爽 | 高潮久久久久久久久久久不卡| 大片电影免费在线观看免费| 色婷婷av一区二区三区视频| 精品国产乱码久久久久久小说| 91国产中文字幕| 一级片免费观看大全| 久久中文字幕一级| 十八禁网站免费在线| 他把我摸到了高潮在线观看 | 亚洲人成电影免费在线| 欧美日韩视频精品一区| 麻豆av在线久日| 久久精品国产亚洲av高清一级| 久久中文字幕一级| 国产精品av久久久久免费| 国产亚洲欧美在线一区二区| 国产欧美日韩精品亚洲av| 亚洲,欧美精品.| 91麻豆精品激情在线观看国产 | 日本一区二区免费在线视频| 后天国语完整版免费观看| 欧美乱妇无乱码| 69精品国产乱码久久久| 嫁个100分男人电影在线观看| 黄色片一级片一级黄色片| 国产精品美女特级片免费视频播放器 | bbb黄色大片| 嫁个100分男人电影在线观看| 亚洲成a人片在线一区二区| 一级毛片电影观看| 高清av免费在线| 夫妻午夜视频| 国产亚洲午夜精品一区二区久久| √禁漫天堂资源中文www| 脱女人内裤的视频| 久热爱精品视频在线9| 午夜两性在线视频| 精品少妇一区二区三区视频日本电影| 九色亚洲精品在线播放| 亚洲一区中文字幕在线| 制服人妻中文乱码| 欧美国产精品一级二级三级| 欧美国产精品一级二级三级| 男女高潮啪啪啪动态图| 欧美精品av麻豆av| 黄色丝袜av网址大全| 成人影院久久| 午夜成年电影在线免费观看| 男女高潮啪啪啪动态图| 成人18禁高潮啪啪吃奶动态图| 欧美人与性动交α欧美软件| 成人av一区二区三区在线看| 三上悠亚av全集在线观看| 亚洲国产欧美日韩在线播放| 每晚都被弄得嗷嗷叫到高潮| 亚洲专区中文字幕在线| 日韩欧美国产一区二区入口| 国产日韩欧美亚洲二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产真人三级小视频在线观看| av欧美777| 久久久国产精品麻豆| 日韩欧美免费精品| 91国产中文字幕| 国产单亲对白刺激| 亚洲 欧美一区二区三区| 变态另类成人亚洲欧美熟女 | 动漫黄色视频在线观看| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区三区在线| 久久久国产欧美日韩av| 国产一卡二卡三卡精品| 视频区欧美日本亚洲| 久久久国产精品麻豆| 欧美激情高清一区二区三区| 在线看a的网站| 又大又爽又粗| 99香蕉大伊视频| 操美女的视频在线观看| 亚洲成国产人片在线观看| 国产精品免费大片| 国产日韩一区二区三区精品不卡| xxxhd国产人妻xxx| 王馨瑶露胸无遮挡在线观看| 手机成人av网站| 国产精品免费视频内射| 亚洲精品中文字幕在线视频| 高清av免费在线| 免费女性裸体啪啪无遮挡网站| videosex国产| 男女边摸边吃奶| 一区二区三区精品91| av福利片在线| 一本色道久久久久久精品综合| 夫妻午夜视频| 日日夜夜操网爽| 一本大道久久a久久精品| 亚洲欧美日韩另类电影网站| 亚洲精品在线美女| 国产精品影院久久| 成年动漫av网址| 一二三四社区在线视频社区8| 91成年电影在线观看| 国产成人精品无人区| 欧美亚洲 丝袜 人妻 在线| 国产极品粉嫩免费观看在线| 欧美国产精品一级二级三级| 青青草视频在线视频观看| 成人永久免费在线观看视频 | 99国产综合亚洲精品| 成年人黄色毛片网站| 成人三级做爰电影| 久久久久国产一级毛片高清牌| 久久久久久亚洲精品国产蜜桃av| 看免费av毛片| 欧美日韩黄片免| 亚洲美女黄片视频| av欧美777| 国产精品九九99| 久久99一区二区三区| 国产成人欧美在线观看 | 亚洲熟妇熟女久久| 宅男免费午夜| 亚洲欧美一区二区三区黑人| 97在线人人人人妻| 久久久久久久国产电影| 国产xxxxx性猛交| 日本精品一区二区三区蜜桃| 香蕉丝袜av| 制服诱惑二区| 成人三级做爰电影| 久久青草综合色| 欧美日韩精品网址| 亚洲专区中文字幕在线| 国产淫语在线视频| 美女高潮到喷水免费观看| 交换朋友夫妻互换小说| 如日韩欧美国产精品一区二区三区| 日本vs欧美在线观看视频| 国产av国产精品国产| 亚洲精品在线美女| 国产精品九九99| 大型av网站在线播放| 亚洲伊人久久精品综合| 熟女少妇亚洲综合色aaa.| 欧美性长视频在线观看| 国产日韩欧美亚洲二区| 日韩欧美免费精品| bbb黄色大片| www.熟女人妻精品国产| 高清欧美精品videossex| 午夜精品国产一区二区电影| 日本精品一区二区三区蜜桃| 高清毛片免费观看视频网站 | 18禁裸乳无遮挡动漫免费视频| 菩萨蛮人人尽说江南好唐韦庄| www.精华液| 男男h啪啪无遮挡| 动漫黄色视频在线观看| 天堂动漫精品| 欧美乱妇无乱码| 国产在线精品亚洲第一网站| 五月天丁香电影| 亚洲精品中文字幕在线视频| 久久国产精品男人的天堂亚洲| www.999成人在线观看| 满18在线观看网站| 可以免费在线观看a视频的电影网站| 日韩免费高清中文字幕av| 另类精品久久| 精品视频人人做人人爽| 国产免费现黄频在线看| 黄片大片在线免费观看| aaaaa片日本免费| 国产精品一区二区免费欧美| 中文字幕精品免费在线观看视频| 欧美 日韩 精品 国产| 男女无遮挡免费网站观看| 人人妻人人澡人人爽人人夜夜| 久久精品亚洲av国产电影网| 午夜激情久久久久久久| 亚洲av片天天在线观看| 999久久久精品免费观看国产| 免费黄频网站在线观看国产| av有码第一页| 日本av免费视频播放| 免费av中文字幕在线| 两性午夜刺激爽爽歪歪视频在线观看 | 老鸭窝网址在线观看| 亚洲第一欧美日韩一区二区三区 | 久久久久久免费高清国产稀缺| 国产一卡二卡三卡精品| 亚洲七黄色美女视频| 99久久精品国产亚洲精品| 免费一级毛片在线播放高清视频 | 国产成人av激情在线播放| 不卡av一区二区三区| 人人妻,人人澡人人爽秒播| 精品人妻熟女毛片av久久网站| 一边摸一边抽搐一进一小说 | 后天国语完整版免费观看| 伦理电影免费视频| 国产精品久久久久久精品电影小说| 国产在线一区二区三区精| 黄色视频不卡| 成人手机av| 久久影院123| 手机成人av网站| 日本欧美视频一区| 黑丝袜美女国产一区| 欧美日韩福利视频一区二区| 亚洲欧洲日产国产| 天堂8中文在线网| 久热爱精品视频在线9| 欧美午夜高清在线| 日本a在线网址| 欧美日韩视频精品一区| 国产精品 欧美亚洲| 日日爽夜夜爽网站| 欧美乱码精品一区二区三区| 超碰97精品在线观看| 一进一出好大好爽视频| 一夜夜www| 91麻豆av在线| 亚洲精品一卡2卡三卡4卡5卡| 精品一区二区三卡| 日韩大码丰满熟妇| 亚洲色图综合在线观看| 侵犯人妻中文字幕一二三四区| 精品国产国语对白av| 久久精品熟女亚洲av麻豆精品| 国产黄色免费在线视频| 如日韩欧美国产精品一区二区三区| 亚洲av第一区精品v没综合| 人人妻人人澡人人看| 亚洲色图综合在线观看| 老司机靠b影院| 欧美日本中文国产一区发布| 亚洲熟女精品中文字幕| 精品国产乱子伦一区二区三区| 纯流量卡能插随身wifi吗| 大香蕉久久网| 精品久久久久久电影网| 精品国产国语对白av| 国产精品久久久av美女十八| 人成视频在线观看免费观看| 成年女人毛片免费观看观看9 | 中文字幕人妻熟女乱码| 免费日韩欧美在线观看| 久久亚洲精品不卡| 亚洲国产欧美一区二区综合| 国产高清激情床上av| 男女免费视频国产| 看免费av毛片| 波多野结衣av一区二区av| 中亚洲国语对白在线视频| 黄色a级毛片大全视频| 国产欧美日韩一区二区三区在线| 亚洲成av片中文字幕在线观看| 欧美精品一区二区免费开放| 淫妇啪啪啪对白视频| 久久久精品区二区三区| 老汉色av国产亚洲站长工具| 无人区码免费观看不卡 | 亚洲,欧美精品.| 国产精品亚洲av一区麻豆| 中文亚洲av片在线观看爽 | 狠狠精品人妻久久久久久综合| 99在线人妻在线中文字幕 | 国产福利在线免费观看视频| 国产亚洲欧美精品永久| 下体分泌物呈黄色| 一进一出好大好爽视频| 中文欧美无线码| 啦啦啦视频在线资源免费观看| 99久久人妻综合| 国产一区二区三区在线臀色熟女 | 首页视频小说图片口味搜索| 中文字幕另类日韩欧美亚洲嫩草| 欧美亚洲 丝袜 人妻 在线| 如日韩欧美国产精品一区二区三区| 亚洲伊人久久精品综合| www.精华液| 国产1区2区3区精品| 国产老妇伦熟女老妇高清| 精品国产超薄肉色丝袜足j| 国产高清激情床上av| 久久久久国内视频| 国产一区二区 视频在线| 美女视频免费永久观看网站| 成年人免费黄色播放视频| 香蕉国产在线看| 久久国产亚洲av麻豆专区| 午夜免费成人在线视频| 欧美激情高清一区二区三区| 免费av中文字幕在线| 成人特级黄色片久久久久久久 | 肉色欧美久久久久久久蜜桃| 亚洲精品国产色婷婷电影| 亚洲av成人不卡在线观看播放网| 久久久久久久久久久久大奶| av国产精品久久久久影院| 男女午夜视频在线观看| 国产免费福利视频在线观看| 男女高潮啪啪啪动态图| 欧美日韩黄片免| 少妇裸体淫交视频免费看高清 | 啦啦啦免费观看视频1| 又大又爽又粗| 丰满少妇做爰视频| 亚洲国产看品久久| 亚洲五月婷婷丁香| 啦啦啦视频在线资源免费观看| 精品人妻1区二区| 大型av网站在线播放| 久久热在线av| 亚洲av成人不卡在线观看播放网| 久久人妻福利社区极品人妻图片| 熟女少妇亚洲综合色aaa.| 激情视频va一区二区三区| videosex国产| 欧美另类亚洲清纯唯美| 久久久久久免费高清国产稀缺| 成人亚洲精品一区在线观看| 大片电影免费在线观看免费| 日本五十路高清| 午夜精品久久久久久毛片777| 在线观看免费视频日本深夜| 一级a爱视频在线免费观看| 欧美黑人欧美精品刺激| 成人国产av品久久久| 女警被强在线播放| 老司机午夜福利在线观看视频 | 黄色 视频免费看| 99国产精品一区二区蜜桃av | 黑人猛操日本美女一级片| 午夜精品国产一区二区电影| 交换朋友夫妻互换小说| netflix在线观看网站| 午夜福利,免费看| 久久精品国产99精品国产亚洲性色 | 大型黄色视频在线免费观看| 日本a在线网址| 国产免费视频播放在线视频| 一区福利在线观看| 人人妻人人澡人人看| 一级,二级,三级黄色视频| 一级片'在线观看视频| 亚洲精品中文字幕在线视频| 男女之事视频高清在线观看| 亚洲精品国产一区二区精华液| 999精品在线视频| 狠狠婷婷综合久久久久久88av| 正在播放国产对白刺激| 亚洲全国av大片| 电影成人av| 一区二区三区国产精品乱码| 国产主播在线观看一区二区| 免费在线观看黄色视频的| 国产精品电影一区二区三区 | 亚洲欧美激情在线| 少妇裸体淫交视频免费看高清 | 亚洲性夜色夜夜综合| 国产日韩一区二区三区精品不卡| 午夜免费成人在线视频| 这个男人来自地球电影免费观看| 咕卡用的链子| 91精品三级在线观看| 人妻久久中文字幕网| 久久人人97超碰香蕉20202| 麻豆av在线久日| av电影中文网址| 欧美日本中文国产一区发布| 午夜成年电影在线免费观看| 日韩欧美国产一区二区入口| 桃红色精品国产亚洲av| 国产精品熟女久久久久浪| 国产亚洲欧美在线一区二区| 中文字幕最新亚洲高清| 精品一区二区三卡| 国产97色在线日韩免费| 老汉色∧v一级毛片| 亚洲精品自拍成人| 在线 av 中文字幕| 国产精品av久久久久免费| 美女福利国产在线| 亚洲精品粉嫩美女一区| 亚洲欧洲精品一区二区精品久久久| 天天影视国产精品| 国产97色在线日韩免费| 国产深夜福利视频在线观看| 国产精品成人在线| 黄色视频在线播放观看不卡| 黄色视频,在线免费观看| 精品一区二区三区av网在线观看 | 男女床上黄色一级片免费看| 高清在线国产一区| 久久毛片免费看一区二区三区| www.精华液| av视频免费观看在线观看| 男女高潮啪啪啪动态图| 国产xxxxx性猛交| 精品国产一区二区久久| 黄色 视频免费看| 欧美亚洲 丝袜 人妻 在线| 国产精品国产高清国产av | 一进一出好大好爽视频| 91麻豆av在线| 新久久久久国产一级毛片| 亚洲人成电影免费在线| 激情视频va一区二区三区| 亚洲黑人精品在线| 99国产精品一区二区三区| 久久国产精品大桥未久av| 啦啦啦在线免费观看视频4| 午夜福利,免费看| 在线播放国产精品三级| 啦啦啦免费观看视频1| 久久久久久亚洲精品国产蜜桃av| 天天躁夜夜躁狠狠躁躁| 成人黄色视频免费在线看| 亚洲精品中文字幕在线视频| 日本一区二区免费在线视频| 99国产极品粉嫩在线观看| 精品一区二区三区四区五区乱码| 成人av一区二区三区在线看| 18在线观看网站| 亚洲精品一卡2卡三卡4卡5卡| 国产成人精品无人区| 老司机亚洲免费影院| 欧美日韩亚洲综合一区二区三区_| 美女福利国产在线| 成人免费观看视频高清| 脱女人内裤的视频| 在线天堂中文资源库| 午夜免费鲁丝| 国产不卡av网站在线观看| 免费av中文字幕在线| 在线观看一区二区三区激情| 操美女的视频在线观看| 亚洲av成人不卡在线观看播放网| 在线观看免费午夜福利视频| 超碰97精品在线观看| 欧美精品一区二区免费开放| 丰满饥渴人妻一区二区三| 超色免费av| 亚洲三区欧美一区| 欧美日韩成人在线一区二区| 在线观看www视频免费| 久久九九热精品免费|