• <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)
    一级黄色大片毛片| 电影成人av| 天天添夜夜摸| 亚洲免费av在线视频| 精品少妇久久久久久888优播| 欧美激情久久久久久爽电影 | 法律面前人人平等表现在哪些方面| 国产在视频线精品| 一进一出抽搐gif免费好疼 | 在线观看午夜福利视频| 成人影院久久| 亚洲五月色婷婷综合| 久久精品熟女亚洲av麻豆精品| 久久中文字幕人妻熟女| 99热网站在线观看| 色综合婷婷激情| 丝袜在线中文字幕| 国产精品九九99| 欧美日韩视频精品一区| 国产在线精品亚洲第一网站| 在线国产一区二区在线| 国产精品 欧美亚洲| 丰满的人妻完整版| 欧美精品高潮呻吟av久久| av天堂在线播放| 老熟妇乱子伦视频在线观看| av电影中文网址| 国产成人精品久久二区二区免费| 久久久久久久久久久久大奶| 19禁男女啪啪无遮挡网站| 亚洲性夜色夜夜综合| 成人亚洲精品一区在线观看| 香蕉丝袜av| 高清毛片免费观看视频网站 | 一区二区三区精品91| 午夜精品久久久久久毛片777| 久久久国产成人免费| 国产精品久久电影中文字幕 | 97人妻天天添夜夜摸| 国产精品免费视频内射| 99国产极品粉嫩在线观看| 久久精品国产亚洲av香蕉五月 | 国产人伦9x9x在线观看| 丁香六月欧美| 天天躁日日躁夜夜躁夜夜| 亚洲精品国产色婷婷电影| 人人妻人人澡人人看| 正在播放国产对白刺激| 国产精品国产高清国产av | 人人妻人人爽人人添夜夜欢视频| 黑丝袜美女国产一区| 亚洲av日韩在线播放| 99香蕉大伊视频| 老司机午夜福利在线观看视频| 国精品久久久久久国模美| 欧美国产精品va在线观看不卡| 12—13女人毛片做爰片一| 久久精品熟女亚洲av麻豆精品| 精品国产一区二区三区久久久樱花| 最新的欧美精品一区二区| 老司机午夜福利在线观看视频| 精品卡一卡二卡四卡免费| 日韩中文字幕欧美一区二区| 99久久人妻综合| 亚洲精品一二三| 俄罗斯特黄特色一大片| 亚洲aⅴ乱码一区二区在线播放 | 国产一卡二卡三卡精品| 国产成人精品在线电影| 日韩人妻精品一区2区三区| 国产区一区二久久| 最近最新中文字幕大全电影3 | 国产成人精品久久二区二区免费| 成年版毛片免费区| 成人国产一区最新在线观看| 欧美午夜高清在线| 美女 人体艺术 gogo| 国产精品美女特级片免费视频播放器 | 亚洲精品自拍成人| 国产成人欧美在线观看 | 美女高潮到喷水免费观看| 18在线观看网站| 婷婷丁香在线五月| 欧美乱码精品一区二区三区| 久久国产精品影院| 成人特级黄色片久久久久久久| 中文字幕制服av| 亚洲精品久久成人aⅴ小说| 午夜福利在线免费观看网站| a在线观看视频网站| 99热国产这里只有精品6| 亚洲精品乱久久久久久| 亚洲在线自拍视频| 老司机亚洲免费影院| 欧美日韩亚洲高清精品| 国产又色又爽无遮挡免费看| 91成年电影在线观看| 一级,二级,三级黄色视频| 亚洲免费av在线视频| 亚洲色图综合在线观看| 国产乱人伦免费视频| 女人高潮潮喷娇喘18禁视频| 99精品欧美一区二区三区四区| 亚洲黑人精品在线| 少妇的丰满在线观看| 精品国产乱子伦一区二区三区| 亚洲av片天天在线观看| 国产精品 国内视频| 欧美亚洲 丝袜 人妻 在线| 99久久精品国产亚洲精品| 成人18禁在线播放| 亚洲精品一卡2卡三卡4卡5卡| 精品人妻1区二区| 91精品三级在线观看| 久久天躁狠狠躁夜夜2o2o| 国产成人啪精品午夜网站| 757午夜福利合集在线观看| 国产精品99久久99久久久不卡| 亚洲欧美色中文字幕在线| 国产淫语在线视频| 国内毛片毛片毛片毛片毛片| 精品一区二区三区四区五区乱码| 正在播放国产对白刺激| 又黄又爽又免费观看的视频| 国产一区有黄有色的免费视频| 午夜亚洲福利在线播放| 精品熟女少妇八av免费久了| 9色porny在线观看| 侵犯人妻中文字幕一二三四区| 欧美黑人精品巨大| 一级,二级,三级黄色视频| 成人三级做爰电影| 欧美乱码精品一区二区三区| av线在线观看网站| 中文字幕人妻丝袜一区二区| 国产欧美日韩精品亚洲av| 51午夜福利影视在线观看| 日本wwww免费看| 久久热在线av| 久久精品国产综合久久久| 亚洲中文字幕日韩| 在线视频色国产色| 深夜精品福利| avwww免费| 国产三级黄色录像| 色94色欧美一区二区| 母亲3免费完整高清在线观看| 欧美日韩一级在线毛片| 在线国产一区二区在线| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻丝袜一区二区| 国产成人影院久久av| 777米奇影视久久| 免费在线观看视频国产中文字幕亚洲| 999精品在线视频| 欧美午夜高清在线| a在线观看视频网站| 国产亚洲欧美精品永久| 香蕉久久夜色| 亚洲美女黄片视频| 99精国产麻豆久久婷婷| 国产精品 国内视频| 亚洲国产中文字幕在线视频| 中文字幕最新亚洲高清| 国产精品久久电影中文字幕 | 不卡av一区二区三区| 一进一出抽搐gif免费好疼 | 丝袜在线中文字幕| 久久天堂一区二区三区四区| 在线观看免费日韩欧美大片| 丰满人妻熟妇乱又伦精品不卡| 亚洲专区字幕在线| 欧美+亚洲+日韩+国产| 亚洲成人国产一区在线观看| 亚洲人成77777在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产日韩一区二区三区精品不卡| 建设人人有责人人尽责人人享有的| 成年动漫av网址| 1024香蕉在线观看| 99国产精品免费福利视频| 精品人妻熟女毛片av久久网站| 女人久久www免费人成看片| 人人妻人人澡人人爽人人夜夜| 日韩视频一区二区在线观看| 青草久久国产| 国产男女内射视频| 免费观看精品视频网站| 9191精品国产免费久久| 欧美 日韩 精品 国产| 免费在线观看完整版高清| 91成年电影在线观看| 1024视频免费在线观看| 国产单亲对白刺激| 欧美日韩国产mv在线观看视频| 九色亚洲精品在线播放| 999久久久国产精品视频| 精品人妻1区二区| 91麻豆精品激情在线观看国产 | 十八禁网站免费在线| 成在线人永久免费视频| 男女下面插进去视频免费观看| 久久久国产成人免费| 欧美乱妇无乱码| 国产黄色免费在线视频| av超薄肉色丝袜交足视频| 精品福利永久在线观看| 久久午夜亚洲精品久久| 波多野结衣一区麻豆| 国产精品秋霞免费鲁丝片| 91麻豆av在线| 国产免费男女视频| 亚洲国产精品sss在线观看 | 欧美日本中文国产一区发布| 一个人免费在线观看的高清视频| 久久久久久久精品吃奶| 人成视频在线观看免费观看| 黄片大片在线免费观看| tocl精华| 99热国产这里只有精品6| 久久久久久亚洲精品国产蜜桃av| 欧美精品啪啪一区二区三区| 日韩成人在线观看一区二区三区| 他把我摸到了高潮在线观看| 国产免费现黄频在线看| 亚洲av第一区精品v没综合| 视频在线观看一区二区三区| 国产熟女午夜一区二区三区| 色婷婷av一区二区三区视频| 久久久国产一区二区| 亚洲av片天天在线观看| a级毛片黄视频| 天堂俺去俺来也www色官网| 婷婷成人精品国产| 丝瓜视频免费看黄片| 日韩制服丝袜自拍偷拍| 午夜福利一区二区在线看| 国产黄色免费在线视频| 久久国产精品人妻蜜桃| 久久精品国产综合久久久| 亚洲专区国产一区二区| 丝瓜视频免费看黄片| 精品国产乱码久久久久久男人| 十八禁人妻一区二区| 日韩有码中文字幕| 91老司机精品| 一二三四在线观看免费中文在| 少妇被粗大的猛进出69影院| 大陆偷拍与自拍| 国产精品久久视频播放| 在线观看免费午夜福利视频| 夜夜夜夜夜久久久久| 丰满人妻熟妇乱又伦精品不卡| 国产激情欧美一区二区| 不卡av一区二区三区| 亚洲中文字幕日韩| 国产男靠女视频免费网站| 一进一出好大好爽视频| 午夜福利在线免费观看网站| 老鸭窝网址在线观看| 日韩有码中文字幕| 搡老熟女国产l中国老女人| 91久久精品国产一区二区成人 | 亚洲人成电影免费在线| 狂野欧美白嫩少妇大欣赏| 夜夜爽天天搞| 日韩欧美免费精品| 观看美女的网站| 天天躁日日操中文字幕| 国产精品嫩草影院av在线观看 | 日日摸夜夜添夜夜添小说| 亚洲一区二区三区不卡视频| 手机成人av网站| 午夜福利在线观看免费完整高清在 | 国产老妇女一区| 嫩草影视91久久| 久久精品国产亚洲av香蕉五月| 亚洲国产欧洲综合997久久,| 一个人免费在线观看电影| 嫁个100分男人电影在线观看| 久久久久性生活片| av黄色大香蕉| 国产伦在线观看视频一区| 99视频精品全部免费 在线| 国产精品影院久久| 免费av毛片视频| 国产精品日韩av在线免费观看| 色播亚洲综合网| 国产成人av教育| 国产av一区在线观看免费| 一区二区三区国产精品乱码| 白带黄色成豆腐渣| 大型黄色视频在线免费观看| 波多野结衣高清无吗| av中文乱码字幕在线| 久久精品国产清高在天天线| 黄色视频,在线免费观看| 俺也久久电影网| 又爽又黄无遮挡网站| 人妻丰满熟妇av一区二区三区| www日本在线高清视频| 我的老师免费观看完整版| 国产一区二区在线观看日韩 | a在线观看视频网站| e午夜精品久久久久久久| 99热这里只有是精品50| 欧美绝顶高潮抽搐喷水| 99久久久亚洲精品蜜臀av| 男女下面进入的视频免费午夜| 国产精品三级大全| 在线观看午夜福利视频| 免费看日本二区| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区免费观看 | 国产久久久一区二区三区| 在线观看免费视频日本深夜| 变态另类丝袜制服| 国产麻豆成人av免费视频| 欧美另类亚洲清纯唯美| 午夜视频国产福利| 女同久久另类99精品国产91| 亚洲乱码一区二区免费版| 亚洲熟妇熟女久久| 欧美成人一区二区免费高清观看| 99久久精品国产亚洲精品| 看免费av毛片| 国产老妇女一区| 色噜噜av男人的天堂激情| 久久久精品欧美日韩精品| 男人和女人高潮做爰伦理| 蜜桃久久精品国产亚洲av| 在线播放无遮挡| 免费在线观看影片大全网站| 国产免费av片在线观看野外av| 极品教师在线免费播放| 欧美日韩一级在线毛片| 中文字幕人妻熟人妻熟丝袜美 | 嫩草影院精品99| 麻豆成人午夜福利视频| 国产aⅴ精品一区二区三区波| 日韩精品中文字幕看吧| 免费看a级黄色片| 琪琪午夜伦伦电影理论片6080| 亚洲成av人片免费观看| 亚洲精品乱码久久久v下载方式 | 色综合站精品国产| 亚洲最大成人手机在线| 国产精品永久免费网站| 国产精品国产高清国产av| 国产野战对白在线观看| 久久久久久大精品| 麻豆国产av国片精品| 国产av一区在线观看免费| 婷婷亚洲欧美| 少妇人妻一区二区三区视频| 午夜福利免费观看在线| 午夜福利视频1000在线观看| 成人18禁在线播放| 精品福利观看| 欧美在线黄色| 免费观看人在逋| 欧美色欧美亚洲另类二区| 午夜免费男女啪啪视频观看 | 国产午夜福利久久久久久| 亚洲激情在线av| 午夜福利在线观看吧| 欧美中文综合在线视频| 国产精品国产高清国产av| 搡老岳熟女国产| 18禁美女被吸乳视频| 亚洲一区高清亚洲精品| 亚洲一区二区三区不卡视频| 欧美日本亚洲视频在线播放| 舔av片在线| 欧美色欧美亚洲另类二区| 午夜免费男女啪啪视频观看 | 日韩精品青青久久久久久| 亚洲内射少妇av| 在线观看日韩欧美| 丁香六月欧美| 欧美一区二区精品小视频在线| 亚洲熟妇熟女久久| 欧美日韩福利视频一区二区| 美女被艹到高潮喷水动态| 国产高潮美女av| 有码 亚洲区| 在线看三级毛片| 成人av一区二区三区在线看| 国产高清有码在线观看视频| 淫妇啪啪啪对白视频| 久久久久免费精品人妻一区二区| 99久久久亚洲精品蜜臀av| 成年女人永久免费观看视频| 在线观看日韩欧美| 小蜜桃在线观看免费完整版高清| 大型黄色视频在线免费观看| 色综合站精品国产| 两人在一起打扑克的视频| 国产亚洲精品一区二区www| 亚洲国产精品sss在线观看| 亚洲av第一区精品v没综合| 欧美日韩亚洲国产一区二区在线观看| 欧美乱色亚洲激情| 亚洲成人中文字幕在线播放| 亚洲七黄色美女视频| 亚洲第一电影网av| 亚洲国产精品sss在线观看| 99热这里只有精品一区| 午夜免费男女啪啪视频观看 | 成人国产一区最新在线观看| 午夜福利免费观看在线| 色视频www国产| 日韩国内少妇激情av| 午夜a级毛片| 午夜福利在线观看吧| 99热这里只有是精品50| 午夜福利高清视频| 舔av片在线| 老司机福利观看| 一级黄色大片毛片| 天堂动漫精品| 国产精品美女特级片免费视频播放器| 亚洲午夜理论影院| 啦啦啦免费观看视频1| 精品久久久久久久末码| 国产成人av激情在线播放| 国产高清视频在线播放一区| 久久久久久大精品| 中出人妻视频一区二区| 免费无遮挡裸体视频| 色吧在线观看| 午夜精品久久久久久毛片777| 91在线观看av| 一级a爱片免费观看的视频| 老鸭窝网址在线观看| 国产 一区 欧美 日韩| ponron亚洲| 香蕉丝袜av| 黄色丝袜av网址大全| 欧美日韩黄片免| 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 麻豆成人av在线观看| 99在线人妻在线中文字幕| 免费在线观看亚洲国产| 久久久久久国产a免费观看| 午夜激情福利司机影院| 国产精品乱码一区二三区的特点| 久久婷婷人人爽人人干人人爱| 青草久久国产| 日韩欧美在线乱码| 一a级毛片在线观看| 给我免费播放毛片高清在线观看| 免费一级毛片在线播放高清视频| 欧美av亚洲av综合av国产av| 一边摸一边抽搐一进一小说| 成人一区二区视频在线观看| 国产野战对白在线观看| 国模一区二区三区四区视频| 欧美最黄视频在线播放免费| 成人三级黄色视频| 十八禁人妻一区二区| 小蜜桃在线观看免费完整版高清| 色噜噜av男人的天堂激情| 日韩亚洲欧美综合| 久久精品国产亚洲av涩爱 | 非洲黑人性xxxx精品又粗又长| 欧美不卡视频在线免费观看| 天堂√8在线中文| 亚洲不卡免费看| 亚洲成人久久性| 日本a在线网址| 精品久久久久久久久久免费视频| 久久久久亚洲av毛片大全| 精品久久久久久久末码| 一区二区三区国产精品乱码| 国产主播在线观看一区二区| 久99久视频精品免费| 日韩精品中文字幕看吧| 床上黄色一级片| 亚洲国产中文字幕在线视频| 国产爱豆传媒在线观看| 老鸭窝网址在线观看| 少妇的逼水好多| 一进一出好大好爽视频| 热99re8久久精品国产| 一级黄色大片毛片| 18禁在线播放成人免费| 亚洲人成网站高清观看| 少妇人妻一区二区三区视频| 亚洲av免费高清在线观看| 国产精品乱码一区二三区的特点| 一二三四社区在线视频社区8| 狂野欧美激情性xxxx| 亚洲熟妇熟女久久| 欧美成人免费av一区二区三区| av在线蜜桃| 成人午夜高清在线视频| av专区在线播放| 男女做爰动态图高潮gif福利片| 很黄的视频免费| 精品国产亚洲在线| 日韩av在线大香蕉| 18+在线观看网站| 亚洲激情在线av| 法律面前人人平等表现在哪些方面| 男女下面进入的视频免费午夜| 99久国产av精品| 少妇丰满av| 亚洲第一欧美日韩一区二区三区| 国产三级在线视频| 成年人黄色毛片网站| 欧美精品啪啪一区二区三区| 久久香蕉国产精品| 在线播放国产精品三级| 亚洲精品一区av在线观看| 欧美乱妇无乱码| 成人国产综合亚洲| 国产精品自产拍在线观看55亚洲| 亚洲人成网站高清观看| 在线观看一区二区三区| 国产日本99.免费观看| 国产99白浆流出| 国产精品久久久久久人妻精品电影| 少妇的逼水好多| 岛国视频午夜一区免费看| 亚洲精品456在线播放app | 久99久视频精品免费| 国产老妇女一区| 成人午夜高清在线视频| 亚洲成av人片在线播放无| 波多野结衣高清作品| 国产精品久久电影中文字幕| 午夜影院日韩av| 可以在线观看的亚洲视频| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久人妻精品电影| 中文资源天堂在线| 女生性感内裤真人,穿戴方法视频| 一本精品99久久精品77| 国产一区二区在线观看日韩 | 国产激情偷乱视频一区二区| 亚洲,欧美精品.| 国产成人啪精品午夜网站| 无遮挡黄片免费观看| 国产中年淑女户外野战色| 色哟哟哟哟哟哟| 国模一区二区三区四区视频| 少妇熟女aⅴ在线视频| 国内毛片毛片毛片毛片毛片| 国产av麻豆久久久久久久| 久久久久久久久久黄片| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品成人久久久久久| 午夜激情欧美在线| 一区二区三区高清视频在线| 亚洲av免费高清在线观看| 九九在线视频观看精品| 欧美黑人巨大hd| 国产精品亚洲美女久久久| 欧美性猛交黑人性爽| 啦啦啦韩国在线观看视频| 亚洲激情在线av| 国内久久婷婷六月综合欲色啪| 一a级毛片在线观看| 国产v大片淫在线免费观看| 亚洲乱码一区二区免费版| 免费在线观看亚洲国产| 久久精品人妻少妇| 一本一本综合久久| 国产伦一二天堂av在线观看| 97超视频在线观看视频| 亚洲不卡免费看| 搡女人真爽免费视频火全软件 | 色在线成人网| 国模一区二区三区四区视频| 搡老熟女国产l中国老女人| 一个人免费在线观看的高清视频| 国产精品野战在线观看| 国产黄a三级三级三级人| 国产精品一区二区免费欧美| 18禁国产床啪视频网站| 尤物成人国产欧美一区二区三区| 国产精品av视频在线免费观看| 亚洲人成网站在线播| av欧美777| 欧洲精品卡2卡3卡4卡5卡区| 男女下面进入的视频免费午夜| 成人亚洲精品av一区二区| 国产一区二区在线观看日韩 | 91九色精品人成在线观看| 变态另类丝袜制服| 国产探花在线观看一区二区| 天堂√8在线中文| 免费人成在线观看视频色| 变态另类丝袜制服| 国产精品野战在线观看| 国产一区二区亚洲精品在线观看| 美女被艹到高潮喷水动态| 国产一区二区在线av高清观看| 日日夜夜操网爽| 成年女人永久免费观看视频| 精品久久久久久久人妻蜜臀av| 少妇的逼好多水| 看片在线看免费视频| 亚洲av电影在线进入| 国产亚洲精品综合一区在线观看| a级毛片a级免费在线| 一区二区三区国产精品乱码| 一个人免费在线观看的高清视频| netflix在线观看网站| 亚洲av一区综合| 手机成人av网站| 成年人黄色毛片网站| 97人妻精品一区二区三区麻豆|