內(nèi)蒙古東北部大豆氣候適宜度等級(jí)及種植區(qū)劃研究*

王彥平, 陰秀霞, 張 昉, 張煦明, 烏長(zhǎng)順

?

內(nèi)蒙古東北部大豆氣候適宜度等級(jí)及種植區(qū)劃研究*

王彥平1, 陰秀霞1, 張 昉2, 張煦明1, 烏長(zhǎng)順3

(1. 內(nèi)蒙古呼倫貝爾市氣象局 海拉爾 021008; 2. 內(nèi)蒙古呼倫貝爾扎蘭屯市氣象局 扎蘭屯 162650; 3. 內(nèi)蒙古呼倫貝爾市新巴爾虎右旗氣象局 阿拉坦額莫勒 021300)

為高效利用內(nèi)蒙古東北部地區(qū)水熱資源, 發(fā)展冷涼地區(qū)大豆種植業(yè), 利用呼倫貝爾市、興安盟及周邊黑龍江地區(qū)近30年的氣象觀測(cè)資料和大豆發(fā)育期資料, 采用氣候適宜度分析方法, 計(jì)算了該區(qū)域近30年逐年大豆氣候適宜度, 建立了氣候適宜度評(píng)價(jià)指標(biāo), 進(jìn)行了基于GIS的呼倫貝爾市大豆適宜種植區(qū)劃。研究結(jié)果表明, 溫度適宜度和水熱綜合適宜度≥0.6為適宜, 0.4~0.6為較適宜, ≤0.4為不適宜; 降水適宜度≥0.7為適宜, 0.6~0.7為較適宜, ≤0.6為不適宜。大豆適宜種植區(qū)域?yàn)榇笈d安嶺農(nóng)區(qū)東南端; 較適宜種植區(qū)域?yàn)榇笈d安嶺東麓, 適宜種植區(qū)的外沿; 不適宜種植區(qū)位于呼倫貝爾市北部林區(qū)及西部牧區(qū)。在林區(qū)南部的農(nóng)林交錯(cuò)帶中, 該區(qū)1 900~2 100 ℃積溫區(qū)域?yàn)橛衩撞贿m宜種植地區(qū), 但可以大力發(fā)展大豆種植業(yè), 在調(diào)整種植業(yè)結(jié)構(gòu)時(shí)可考慮在這一區(qū)域增加大豆種植面積。

內(nèi)蒙古東北部; 氣候適宜度; 評(píng)價(jià)指標(biāo); 大豆種植; 區(qū)劃方法

內(nèi)蒙古東北部是全國(guó)重要的優(yōu)質(zhì)大豆[(L.) Merrill]生產(chǎn)基地, 也是東北大豆產(chǎn)區(qū)的重要組成部分, 其種植面積僅次于玉米(L.)。2016—2017年, 內(nèi)蒙古自治區(qū)根據(jù)市場(chǎng)需求變化, 對(duì)種植業(yè)進(jìn)行了科學(xué)調(diào)整, 破解種植業(yè)的玉米作物“一糧獨(dú)大”問(wèn)題, 逐步增加?xùn)|北部地區(qū)大豆的種植面積。

氣候適宜度方法近幾年被廣泛用于農(nóng)作物和氣候條件之間適應(yīng)程度的研究, 主要應(yīng)用領(lǐng)域?yàn)殚_展氣候變化背景下農(nóng)業(yè)氣候資源的適宜性評(píng)估、農(nóng)業(yè)應(yīng)對(duì)氣候變化及合理開發(fā)和高效利用農(nóng)業(yè)氣候資源等[1]。

在氣候變化對(duì)農(nóng)作物適宜度的評(píng)價(jià)方面, 陳炳東[2]總結(jié)了氣候資源對(duì)農(nóng)作物適宜度影響的一般研究思路、研究?jī)?nèi)容和研究進(jìn)展, 指出進(jìn)行氣候變化對(duì)于農(nóng)作物適宜度影響的研究中應(yīng)當(dāng)首先進(jìn)行農(nóng)業(yè)氣候資源適宜度模式的建立, 然后根據(jù)具體研究方法結(jié)合農(nóng)作物情況進(jìn)行相關(guān)研究計(jì)算。隨著農(nóng)業(yè)氣候適宜度的提出, 國(guó)內(nèi)不少學(xué)者對(duì)不同尺度區(qū)域的農(nóng)業(yè)生態(tài)氣候適度進(jìn)行了系統(tǒng)的分析探討。劉實(shí)等[3]研究了近50年?yáng)|北地區(qū)熱量資源變化特征, 得出氣溫和積溫均持續(xù)增加利于農(nóng)作物生長(zhǎng)和產(chǎn)量提高; 冶明珠等[4-5]研究了氣候變化背景下東北地區(qū)熱量資源及玉米溫度適宜度變化特征, 得出東北地區(qū)熱量資源顯著增加, 玉米溫度適宜度也逐漸升高。在國(guó)外, 氣候變化對(duì)農(nóng)作物影響的聚集點(diǎn)都是對(duì)氣候變化潛在的或者已經(jīng)呈現(xiàn)出來(lái)的對(duì)農(nóng)作物產(chǎn)量、水平衡、農(nóng)業(yè)系統(tǒng)組成成分影響等問(wèn)題展開的分析。Matthews等[6]對(duì)比了1989年在菲律賓調(diào)查的氣候變化和水稻產(chǎn)量關(guān)系的研究, 利用模型模擬量化氣候變化和水稻產(chǎn)量的關(guān)系, 最后擴(kuò)展到中尺度甚至大尺度研究; Mestre-Sanchís等[7]研究表明, 氣候變化會(huì)導(dǎo)致某些區(qū)域作物生境惡化, 隨著深入研究進(jìn)而發(fā)現(xiàn)氣候變化不僅影響著自然環(huán)境及人類社會(huì), 同時(shí)還影響了農(nóng)作物生長(zhǎng)的適宜程度。

2000年以來(lái), 國(guó)內(nèi)學(xué)者在氣候適宜度領(lǐng)域的研究主要集中在4方面: 1)氣候要素的適宜度曲線、適宜態(tài)與隸屬函數(shù)的建立。王麗霞等[8]對(duì)山西大同市農(nóng)業(yè)生態(tài)氣候適宜度進(jìn)行評(píng)價(jià), 擬合出符合區(qū)域特點(diǎn)的典型氣候要素適宜度隸屬函數(shù), 分區(qū)計(jì)算了農(nóng)業(yè)生態(tài)資源指數(shù)、效能指數(shù)及利用系數(shù)。2)農(nóng)業(yè)氣候適宜度模型的建立。魏瑞江等[9]運(yùn)用模糊數(shù)學(xué)理論, 結(jié)合前人研究成果, 建立了河北省冬小麥(L.)溫度、降水、日照隸屬度模型和氣候適宜度模型; 俞芬等[10]對(duì)淮河流域水稻(L.)的氣候適宜及其變化趨勢(shì)進(jìn)行分析, 得出淮河流域水稻溫度、降水、日照和氣候適宜度的區(qū)域差異, 將淮河流域分為最適宜區(qū)、適宜性區(qū)、次適宜區(qū); 劉國(guó)成等[11]建立了基于模糊數(shù)學(xué)的農(nóng)業(yè)氣候適宜度劃分方法, 結(jié)合專家系統(tǒng)和地理信息系統(tǒng)技術(shù)對(duì)吉林省農(nóng)業(yè)氣候生態(tài)適宜度以鄉(xiāng)鎮(zhèn)為最小單位進(jìn)行劃分; 黃淑娥等[12]通過(guò)對(duì)江西省雙季水稻生長(zhǎng)季溫、光、水的需求及當(dāng)?shù)貧夂驐l件進(jìn)行分析, 建立了雙季水稻各生育期光、溫、水及氣候適宜度評(píng)價(jià)模型。3)不同尺度農(nóng)業(yè)生態(tài)氣候資源的分區(qū)與評(píng)價(jià)研究。任玉玉等[13]對(duì)河南省棉花(sp.)氣候適宜度變化趨勢(shì)進(jìn)行分析, 得出氣候適宜度總體呈下降的變化趨勢(shì), 各地的變化趨勢(shì)依據(jù)熱量帶和地形地貌的不同有明顯的地域差異, 并結(jié)合各地的氣候資料對(duì)各類型的氣候適宜度變化原因進(jìn)行了初步分析。4)運(yùn)用農(nóng)業(yè)生態(tài)氣候適宜度理論指導(dǎo)農(nóng)業(yè)生產(chǎn)實(shí)際方面的研究。馬興祥等[14]將生態(tài)氣候適宜度應(yīng)用到甘肅省春小麥種植區(qū)劃中, 用生態(tài)氣候適應(yīng)性對(duì)春小麥的影響程度和依賴度及產(chǎn)量, 確定生態(tài)氣候區(qū)劃綜合指標(biāo)體系將春小麥劃分出5級(jí)適生種植區(qū)劃, 并提出利用和開發(fā)生態(tài)氣候資源有效途徑; 魏瑞江等[15]研究了基于氣候適宜度的玉米產(chǎn)量動(dòng)態(tài)預(yù)報(bào)方法, 建立了夏玉米氣候適宜度模型和不同時(shí)段產(chǎn)量預(yù)報(bào)模型; 千懷遂等[16]以河南省冬小麥為例研究了全球氣候變化對(duì)作物氣候適宜性的影響, 分析發(fā)現(xiàn), 受近年來(lái)全球變暖的影響, 冬小麥拔節(jié)期溫度與降水適宜度明顯提高。

目前對(duì)大豆適應(yīng)性研究的報(bào)道較少。卜坤等[17]和Fischer等[18]依據(jù)FAO/IIASA發(fā)布的GAEZ研究框架, 對(duì)三江平原大豆種植進(jìn)行土壤適宜性評(píng)價(jià), 引入了適宜度指數(shù)(SI), 對(duì)三江平原大豆種植的土地適宜度進(jìn)行了空間上的比較。東北春大豆種植區(qū)是我國(guó)大豆第一大產(chǎn)區(qū), 也是我國(guó)大豆商品基地, 包括黑龍江、吉林、遼寧3省和內(nèi)蒙古東北部地區(qū), 因此對(duì)東北地區(qū)大豆氣候適宜性的研究對(duì)我國(guó)大豆生產(chǎn)具有十分重要的意義。楊顯峰等[19]提出春大豆氣候適宜性數(shù)學(xué)表達(dá)方式, 根據(jù)春大豆的生育特點(diǎn)及其對(duì)環(huán)境條件的要求確定了適宜春大豆生長(zhǎng)發(fā)育的適宜降水區(qū)間、適宜氣溫區(qū)間和適宜日照區(qū)間, 為運(yùn)用模糊數(shù)學(xué)隸屬函數(shù)的方法建立東北春大豆與溫度、日照、降水等氣候要素之間的適宜性關(guān)系奠定了基礎(chǔ)。沈國(guó)強(qiáng)[20]運(yùn)用MaxEnt模型, 給出了影響東北地區(qū)大豆氣候適宜性的氣候因子綜合貢獻(xiàn)率, 得到基于所有氣候因子的東北地區(qū)大豆氣候適宜性模型, 進(jìn)而得出大豆氣候適宜性分布的年代際變化。運(yùn)用隸屬函數(shù)建立氣候大豆適宜度評(píng)價(jià)模型, 并以適宜度指標(biāo)作為種植區(qū)劃因子, 結(jié)合地理信息系統(tǒng)進(jìn)行種植區(qū)劃, 是總結(jié)以往研究成果進(jìn)行種植區(qū)劃研究農(nóng)作物氣候適宜性分布的發(fā)展趨勢(shì)。

本文以內(nèi)蒙古東北部大豆種植基地呼倫貝爾市作為研究對(duì)象, 其氣候條件和種植面積均可代表內(nèi)蒙古東北部大豆種植現(xiàn)狀。從當(dāng)?shù)剞r(nóng)業(yè)部門了解到, 大豆氣候資源利用目前存在兩大問(wèn)題: 一是遵循傳統(tǒng)種植習(xí)慣, 大豆種植區(qū)域和玉米種植區(qū)域沒有明確劃分, 利用≥10 ℃積溫劃分大豆種植地區(qū), 遠(yuǎn)遠(yuǎn)不能滿足生產(chǎn)要求, 缺少專門針對(duì)大豆的種植區(qū)劃; 二是目前大豆種植業(yè)中熱量和水分資源利用效率較低, 不能充分利用呼倫貝爾大豆種植的豐富水熱資源。本文擬對(duì)內(nèi)蒙古東北部大豆種植區(qū)的呼倫貝爾市大豆氣候適宜度進(jìn)行分析, 建立大豆氣候適宜度劃分指標(biāo), 并在此基礎(chǔ)上進(jìn)行基于GIS的大豆種植區(qū)劃研究, 形成較為科學(xué)的區(qū)劃結(jié)論。通過(guò)本區(qū)劃可以改變多年以玉米種植區(qū)域劃分大豆種植區(qū)域的模糊種植區(qū)劃, 高效利用呼倫貝爾市的熱量和水分資源, 為發(fā)展冷涼地區(qū)大豆種植業(yè), 指導(dǎo)大豆生產(chǎn)提供一定理論依據(jù)。

1 研究區(qū)域概況及資料來(lái)源

1.1 研究區(qū)自然特征

內(nèi)蒙古東北部的呼倫貝爾市地處呼倫貝爾高原, 介于115°31′~126°04′E, 47°05′~53°20′N, 境內(nèi)地形復(fù)雜多樣, 地勢(shì)起伏多變, 包括大興安嶺山地、呼倫貝爾高原、河谷平原低地3個(gè)較大的地形單元, 大興安嶺呈北東北—南西南走向, 貫穿本區(qū)中部, 構(gòu)成呼倫貝爾市地形的整體。以大興安嶺為界, 全市氣溫東西差異較大, 大興安嶺山地為寒冷作物帶, 西部草原為溫涼作物帶, 嶺東為溫和作物帶, 且氣候資源具有明顯的水平地帶性(緯向)和垂直地帶性, 其中熱量資源體現(xiàn)更為明顯。

1.2 資料來(lái)源

1.2.1 地理信息資料

地理信息數(shù)據(jù)來(lái)源于國(guó)家測(cè)繪局1∶25萬(wàn)地理信息數(shù)據(jù), 包括呼倫貝爾市鄉(xiāng)鎮(zhèn)行政邊界、道路、數(shù)字高程?hào)鸥駭?shù)據(jù)等地理信息資料; 2007年版各市/縣、各鎮(zhèn)的行政區(qū)劃圖件以及各氣象觀測(cè)站點(diǎn)的經(jīng)緯度坐標(biāo)。

1.2.2 氣象數(shù)據(jù)及大豆發(fā)育期資料

為了建立小網(wǎng)格推算模型, 本文氣象數(shù)據(jù)以呼倫貝爾市(16個(gè)臺(tái)站)為主, 黑龍江省(富裕、甘南、加格達(dá)奇、龍江、訥河、齊齊哈爾、依安)和興安盟(胡爾勒、阿爾山、索倫、音德爾)等11個(gè)臺(tái)站為輔(主要參與建立模型), 利用連續(xù)30年歷史氣象數(shù)據(jù), 包括1987—2016年27個(gè)氣象觀測(cè)站點(diǎn)的日平均氣溫、日降水量數(shù)據(jù); 呼倫貝爾市和興安盟大豆各發(fā)育期資料以扎蘭屯農(nóng)業(yè)氣象監(jiān)測(cè)站資料為準(zhǔn), 黑龍江各站以龍江農(nóng)業(yè)氣象監(jiān)測(cè)站資料為準(zhǔn)。

2 分析方法

2.1 地理信息資料的處理

應(yīng)用1∶25萬(wàn)地理信息數(shù)據(jù), 在地理信息系統(tǒng)平臺(tái)下, 進(jìn)行各層地理信息數(shù)據(jù)的導(dǎo)入、拼接、拓?fù)溥\(yùn)算和投影、配準(zhǔn)等圖層處理, 最終生成呼倫貝爾市的海拔高度、經(jīng)度和緯度柵格圖層。

2.2 精細(xì)化種植區(qū)劃方法

區(qū)劃采用氣候適宜度分析方法, 計(jì)算27個(gè)氣象監(jiān)測(cè)站30年氣候適宜度平均值, 建立氣候適宜度與地理信息的小網(wǎng)格推算模型, 充分利用GIS強(qiáng)大的空間分析功能, 結(jié)合農(nóng)業(yè)氣候相似理論和地域分布規(guī)律, 實(shí)現(xiàn)由點(diǎn)到面的空間合并, 最終生成基于氣候適宜度的大豆種植區(qū)劃圖, 并與傳統(tǒng)的≥10 ℃積溫區(qū)劃圖進(jìn)行對(duì)比分析, 進(jìn)行分區(qū)評(píng)述。

2.2.1 氣候適宜度評(píng)價(jià)模型的建立

2.2.1.1 溫度適宜度函數(shù)

以往研究表明, Beta函數(shù)能較好地反映作物生長(zhǎng)與溫度的關(guān)系, 并且具有普適性, 其值為0~1[21]。本研究采用該方法進(jìn)行溫度適宜度函數(shù)的建立, 計(jì)算方法見公式(1)和(2)。

(2)

式中:()為某一階段溫度適宜度;為大豆某發(fā)育階段日平均氣溫;l、h和o分別為大豆各發(fā)育期所需的下限溫度、上限溫度和適宜溫度, 參考東北地區(qū)春大豆指標(biāo)體系[19], 結(jié)合內(nèi)蒙古東部呼倫貝爾市和興安盟生產(chǎn)實(shí)踐, 確定大豆各生育階段三基點(diǎn)溫度(表1)。

表1 大豆各生育階段三基點(diǎn)溫度指標(biāo)

2.2.1.2 降水適宜度函數(shù)

降水是旱作地區(qū)大豆需水的主要來(lái)源, 是制約產(chǎn)量的主要?dú)庀笠蜃?。降水適宜度表征了種植地區(qū)大豆各生長(zhǎng)發(fā)育階段內(nèi)的降水量對(duì)其生長(zhǎng)發(fā)育和產(chǎn)量形成的適宜程度。研究認(rèn)為, 用作物正常生長(zhǎng)需水量作為作物生長(zhǎng)的適宜水量標(biāo)準(zhǔn)是可行的。

采用以往研究成果, 以0.7倍為旱與正常的分界線, 作物某一生育階段降水量/需水量<0.7時(shí), 出現(xiàn)輕旱; ≥0.7時(shí)為降水適宜大豆生長(zhǎng)發(fā)育[21]。因此, 降水適宜度模型見式(3)。

式中:()為大豆生育期旬降水適宜度,為某階段的降水量, ETc為大豆相應(yīng)階段的理論需水量(表2)。

2.2.1.3 綜合氣候適宜度模型

由于呼倫貝爾市大豆生長(zhǎng)季日照充足, 日照時(shí)數(shù)不是其生長(zhǎng)發(fā)育和產(chǎn)量形成的主導(dǎo)因子, 因此進(jìn)行綜合適宜度模型建立時(shí)只考慮溫度和降水的影響,所建模型為水熱綜合適宜度模型。采用層次分析方法(AHP)確定大豆各發(fā)育階段溫度和降水的權(quán)重值(表2)。AHP是將與決策總是有關(guān)的元素分解成目標(biāo)、準(zhǔn)則、方案等層次, 在此基礎(chǔ)上進(jìn)行定性和定量分析的決策方法[22-24]。根據(jù)文獻(xiàn), 建立了大豆氣溫和降水綜合氣候適宜度模型, 見公式(4)。

表2 大豆各生育期理論需水量

式中:m為某一生育期(具體為該生育期普期至下一生育期普期之間的時(shí)間間隔)氣候適宜度,m()為該生育期溫度適宜度,bt為該生育期溫度適宜度的權(quán)重系數(shù),m()為該生育期降水適宜度,br為該階段降水適宜度的權(quán)重系數(shù)(表3)。

表3 呼倫貝爾地區(qū)大豆各生育期氣候適宜度權(quán)重系數(shù)

b為該發(fā)育期在整個(gè)生育期的權(quán)重系數(shù);bt為該發(fā)育期溫度的權(quán)重系數(shù);br為該發(fā)育期降水的權(quán)重系數(shù)。b: weight coefficient of the growth phaserelative to the whole growth period;bt: weight coefficient of temperature of the growth phase;br: weight coefficient of precipitation in of the growth phase.

2.3 氣候適宜度推算模型建立及區(qū)劃圖制作

2.3.1 呼倫貝爾市各地氣候適宜度分析結(jié)果

計(jì)算呼倫貝爾市16個(gè)臺(tái)站溫度適宜度、降水適宜度及水熱綜合適宜度30年平均值, 統(tǒng)計(jì)造成各地大豆生長(zhǎng)發(fā)育不適宜年份并進(jìn)行原因分析(表4)。

2.3.2 呼倫貝爾市大豆種植氣候適宜度評(píng)價(jià)指標(biāo)

根據(jù)表4分析結(jié)果, 結(jié)合呼倫貝爾地區(qū)大豆種植生產(chǎn)實(shí)際及相應(yīng)單產(chǎn)數(shù)據(jù)進(jìn)行綜合分析、評(píng)判, 確定了呼倫貝爾地區(qū)大豆種植氣候適宜度劃分等級(jí)標(biāo)準(zhǔn), 見公式(5)~(7)。

表4 呼倫貝爾市各地大豆氣候適宜度平均值及適宜度分析(1987—2016)

溫度適宜度評(píng)價(jià)模型:

降水適宜度評(píng)價(jià)指標(biāo):

水熱綜合適宜度評(píng)價(jià)指標(biāo)

當(dāng)溫度適宜度或降水適宜度評(píng)價(jià)指標(biāo)為不適宜時(shí), 水熱適宜度即為不適宜。

2.3.3 呼倫貝爾市大豆種植氣候適宜度推算模型

以往研究表明, 氣象要素與經(jīng)度、緯度和海拔的數(shù)量關(guān)系, 使用三元線性回歸模型較三元二次多項(xiàng)式回歸模型穩(wěn)定性強(qiáng), 因此, 本研究使用多元線性模型。利用27個(gè)站點(diǎn)地理信息和溫度適宜度、降水適宜度及水熱綜合適宜度建立多元線性模型(表5)。

表5 呼倫貝爾市大豆氣候適宜度與地理信息的三元線性模型

方程中為經(jīng)度,為緯度,為海拔高度。is longitude,is latitudes, andis altitude in the equations.

2.3.4 呼倫貝爾大豆種植精細(xì)區(qū)劃圖制作

在GIS系統(tǒng)下, 根據(jù)生成的經(jīng)度、緯度和高度柵格圖層, 應(yīng)用表3的推算方程及適宜度評(píng)價(jià)模型, 繪制出溫度適宜度、降水適宜度及水熱綜合適宜度區(qū)劃圖。

3 結(jié)果與分析

3.1 呼倫貝爾大豆種植的溫度適宜度變化及區(qū)劃

由圖1a可知, 從溫度適宜度分析, 呼倫貝爾可滿足大豆生長(zhǎng)發(fā)育需要的區(qū)域明顯分為兩大部分。第1部分為呼倫貝爾市東部的農(nóng)區(qū), 該區(qū)域大豆生長(zhǎng)季(5—9月)平均氣溫在14.7~17.8 ℃, 基本滿足大豆生長(zhǎng)發(fā)育需要。其中扎蘭屯市東南部邊緣、阿榮旗東南部、莫旗南部溫度適宜度在0.6以上, 處于大豆生長(zhǎng)發(fā)育的溫度適宜等級(jí); 扎蘭屯市中部、阿榮旗中部、莫旗大部和鄂倫春旗東南部地區(qū)溫度適宜度為0.4~0.6, 處于較適宜等級(jí)。第2部分為大興安嶺以西的牧區(qū),該區(qū)域5—9月平均氣溫在15.6~17.3℃, 可以滿足大豆生長(zhǎng)期對(duì)溫度的要求。其中新右旗中部和南部地區(qū)溫度適宜度為適宜等級(jí), 新右旗北部、新左旗、陳旗大部、鄂溫克旗及海拉爾區(qū)部分地區(qū)處于較適宜等級(jí)。

不能滿足大豆生長(zhǎng)發(fā)育對(duì)溫度需求的區(qū)域主要位于大興安嶺山地的廣大林區(qū)和嶺西牧區(qū)東部地區(qū), 該區(qū)域5—9月平均氣溫為12.2~15.0 ℃, 大部地區(qū)無(wú)霜期不足100 d, 種植大豆不能成熟。

3.2 降水適宜度區(qū)劃結(jié)果分析

由圖1b可知, 從降水適宜度分析, 呼倫貝爾大豆生長(zhǎng)季降水適宜度由東向西逐漸減小, 地帶性明顯。第1部分為大興安嶺東麓的東北部地區(qū), 包括鄂倫春旗東北部和莫旗大部地區(qū), 該區(qū)域年降水量420~490 mm, 降水適宜度在0.7以上, 處于適宜等級(jí), 出現(xiàn)干旱的幾率較小, 水分條件完全可以滿足大豆生長(zhǎng)發(fā)育需要。第2部分為大興安嶺東麓東南部的農(nóng)區(qū)和大興安嶺山地的林區(qū), 包括扎蘭屯市、阿榮旗、根河市、牙克石市以及額爾古納市東部和鄂溫克旗東部, 該區(qū)域年降水量400 mm左右, 降水適宜度為0.6~0.7, 處于較適宜等級(jí), 水分條件基本可以滿足大豆生長(zhǎng)發(fā)育需要, 但出現(xiàn)干旱的幾率較大。第3部分為大興安嶺以西的牧區(qū), 該區(qū)域年降水量不足300 mm, 降水適宜度在0.6以下, 處于不適宜等級(jí), 不能滿足大豆生長(zhǎng)發(fā)育對(duì)水分的需求。

3.3 水熱綜合適宜度區(qū)劃結(jié)果分析

由圖2可見, 從水熱綜合適宜度分析, 呼倫貝爾市大豆適宜種植區(qū)域?yàn)榇笈d安嶺東麓農(nóng)區(qū)東南端, 與溫度適宜度區(qū)域基本重合, 包括扎蘭屯市東南端、阿榮旗東部和莫旗南部, 水熱綜合適宜度在0.6以上, 該地區(qū)土壤以黑鈣土、黑土、暗色草甸土、暗棕壤為主, 地勢(shì)以平原為主, 土壤肥沃, 水資源豐富。該區(qū)的主要?dú)夂蛱攸c(diǎn)是光熱資源最豐富, 雨熱同季, 無(wú)霜期較長(zhǎng), 年降水量490~520 mm, ≥10 ℃活動(dòng)積溫2 300~2 500 ℃×d, ≥2 ℃無(wú)霜期120~130 d?？蓾M足大多數(shù)早、中熟大豆品種的生長(zhǎng)發(fā)育需要, 多數(shù)年份可獲得較好的收成, 且所產(chǎn)大豆產(chǎn)量和品質(zhì)都比較高。

圖1 呼倫貝爾市大豆種植溫度適宜度(a)和降水適宜度(b)區(qū)劃

圖2 呼倫貝爾市大豆種植水熱綜合適宜度區(qū)劃

大豆較適宜種植區(qū)域地處大興安嶺東麓, 適宜種植區(qū)的外沿, 水熱綜合適宜度在0.4~0.6, 該地區(qū)土壤以暗棕壤為主, 土壤較為肥沃, 水資源豐富, 主要受熱量條件限制。本區(qū)氣候特點(diǎn)是熱量偏少, 雨熱同季, 年降水量490~560 mm, ≥10 ℃活動(dòng)積溫1 900~2 300 ℃×d, 只能滿足部分早熟品種大豆的生長(zhǎng)發(fā)育需求, ≥2 ℃無(wú)霜期較短, 為110~120 d, 常會(huì)受到低溫和霜凍的侵襲, 因此該區(qū)域內(nèi)大豆不宜大面積種植, 推廣時(shí)需要密切關(guān)注當(dāng)年的氣候預(yù)測(cè), 合理選擇大豆品種。

大豆不適宜種植區(qū)位于呼倫貝爾市北部林區(qū)及西部牧區(qū), 水熱綜合適宜度在0.4以下。北部林區(qū)主要是由于氣溫偏低, 該地區(qū)≥2 ℃無(wú)霜期小于115 d, ≥10 ℃活動(dòng)積溫小于1 900 ℃×d, 熱量條件不能滿足大豆生長(zhǎng)的基本要求, 多數(shù)年份無(wú)法獲得收成。該地區(qū)以山地為主, 地勢(shì)崎嶇, 森林覆蓋度高, 不利于大豆栽培。呼倫貝爾牧區(qū)的南部, 新左旗、新右旗的熱量條件雖然可以達(dá)到較適宜程度, 但年降水量較少, 不能滿足大豆生育期內(nèi)水分要求, 在無(wú)灌溉條件下, 不適宜大豆生長(zhǎng), 且此地區(qū)以牧業(yè)生產(chǎn)為主, 不宜推廣大豆種植。牧區(qū)北部的熱量條件不能滿足大豆生長(zhǎng)的基本條件, 不宜種植大豆。

4 討論與結(jié)論

呼倫貝爾地區(qū)大豆種植沒有明確的區(qū)劃方法和區(qū)劃圖, 多年來(lái)以玉米種植積溫帶界定大豆種植區(qū)域, 認(rèn)為適宜種玉米的區(qū)域也適宜大豆種植。目前, 大豆種植區(qū)劃采用傳統(tǒng)的≥10 ℃積溫區(qū)劃和分級(jí)判別區(qū)劃方法, 前者沒有考慮水熱資源的綜合效應(yīng), 后者沒有體現(xiàn)水熱條件對(duì)各生育期階段的不同效應(yīng), 本研究利用內(nèi)蒙古和黑龍江大豆主產(chǎn)區(qū)的氣候適宜度確定的評(píng)價(jià)指標(biāo), 對(duì)內(nèi)蒙古呼倫貝爾市大豆產(chǎn)區(qū)進(jìn)行了基于適宜度方法的種植區(qū)劃。

4.1 大豆氣候適宜度區(qū)劃結(jié)果與≥10 ℃積溫區(qū)劃結(jié)果對(duì)比分析

大豆適宜度種植區(qū)劃圖(圖3a)與80%保證率下≥10 ℃活動(dòng)積溫圖[25](圖3b)對(duì)比分析發(fā)現(xiàn), 呼倫貝爾市大豆適宜種植區(qū)域處于2 300~2 500 ℃積溫帶, 具體分布于扎蘭屯市、阿榮旗南部大部地區(qū); ≥2 500 ℃積溫帶, 分布于扎蘭屯市、阿榮旗及莫旗最東南端邊緣, 這兩個(gè)積溫帶為春播玉米的中熟區(qū)和中早熟區(qū), 為玉米的最適宜和適宜種植區(qū), 這一區(qū)域也為大豆適宜種植區(qū)域, 與傳統(tǒng)種植理念十分吻合。

大豆較適宜種植區(qū)域較為廣泛, 處于2 100~ 2 300 ℃積溫帶, 包括扎蘭屯市、阿榮旗北部鄉(xiāng)鎮(zhèn)、莫旗大部鄉(xiāng)鎮(zhèn)及鄂倫春旗東南部鄉(xiāng)鎮(zhèn), 這些地區(qū)也是玉米較適宜種植區(qū)域; 1 900~2 100 ℃積溫帶, 包括扎蘭屯市、阿榮旗北部邊界鄉(xiāng)鎮(zhèn)及鄂倫春旗中部鄉(xiāng)鎮(zhèn), 由于熱量資源欠缺和早霜的危害, 該區(qū)域?yàn)橛衩撞贿m宜種植地區(qū), 但可以大力發(fā)展大豆種植業(yè), 在不大面積縮減玉米種植面積的同時(shí)以增加大豆種植面積。

圖3 呼倫貝爾市大豆氣候適宜度區(qū)劃(a)與≥10 ℃積溫區(qū)劃(b)對(duì)比圖

4.2 大豆氣候適宜度區(qū)劃結(jié)果與分級(jí)判別區(qū)劃結(jié)果對(duì)比分析

參考玉米分級(jí)判別區(qū)劃方法[25], 選用日平均溫度≥10 ℃活動(dòng)積溫、年降水量、日最低氣溫≥2 ℃無(wú)霜期作為大豆區(qū)劃指標(biāo)進(jìn)行區(qū)劃與本研究區(qū)劃結(jié)果對(duì)比分析。由圖4可知, 與分級(jí)判別區(qū)劃結(jié)果相比, 氣候適宜度判別區(qū)劃結(jié)果中, 大豆適宜種植區(qū)域邊緣向東南縮減, 區(qū)域面積有所減小, 較適宜種植區(qū)域邊緣向北延伸, 區(qū)域面積明顯增加。出現(xiàn)差別的原因是分級(jí)判別雖然也綜合考慮了水熱的綜合影響效應(yīng), 但是以整個(gè)發(fā)育期的氣象條件作為適宜區(qū)劃因子, 沒有考慮作物各個(gè)發(fā)育期對(duì)氣象條件需求的動(dòng)態(tài)變化, 而適宜度區(qū)劃則是在綜合考慮作物各個(gè)發(fā)育期水熱適宜程度的基礎(chǔ)上, 確定各發(fā)育期適宜指標(biāo)及權(quán)重系數(shù), 再確定整個(gè)發(fā)育階段的適宜指標(biāo), 在此基礎(chǔ)上進(jìn)行區(qū)劃, 其結(jié)果更能體現(xiàn)大豆整個(gè)生長(zhǎng)期對(duì)氣候條件的適宜程度, 對(duì)擴(kuò)大種植面積有更科學(xué)的理論依據(jù)和實(shí)用價(jià)值。

圖4 呼倫貝爾市大豆氣候適宜度區(qū)劃(a)與分級(jí)判別區(qū)劃(b)對(duì)比圖

4.3 氣候適宜度區(qū)劃方法的利弊分析

氣候適宜度區(qū)劃方法不僅綜合考慮了熱量和水分條件對(duì)大豆生長(zhǎng)發(fā)育的共同影響, 而且考慮了大豆各個(gè)生育階段及全生育期對(duì)水熱的要求, 所得區(qū)劃結(jié)論比≥10 ℃活動(dòng)積溫區(qū)劃方法和分級(jí)判別區(qū)劃方法更能真實(shí)體現(xiàn)大豆生長(zhǎng)全過(guò)程對(duì)氣象條件的具體要求。由于氣候適宜度存在諸多計(jì)算過(guò)程, 需要大豆生長(zhǎng)的溫度和降水指標(biāo)作為計(jì)算的基礎(chǔ), 因此存在指標(biāo)確定是否合理的弊端, 溫度三基點(diǎn)指標(biāo)及理論需水量指標(biāo)仍需在生產(chǎn)實(shí)際中進(jìn)行驗(yàn)證, 區(qū)劃結(jié)論還需不斷完善。

4.4 結(jié)論

呼倫貝爾市溫度適宜度和水熱綜合適宜度≥0.6為適宜, 0.4~0.6為較適宜, ≤0.4為不適宜; 降水適宜度≥0.7為適宜, 0.6~0.7為較適宜, ≤0.6為不適宜。大豆適宜種植區(qū)域?yàn)榇笈d安嶺東麓農(nóng)區(qū)東南端, 水熱綜合適宜度在0.6以上; 較適宜種植區(qū)域?yàn)榇笈d安嶺東麓, 適宜種植區(qū)的外沿, 水熱綜合適宜度為0.4~0.6; 不適宜種植區(qū)位于呼倫貝爾市北部林區(qū)及西部牧區(qū), 水熱綜合適宜度在0.4以下, 且1 900~2 100 ℃積溫區(qū)域不適宜種植玉米地區(qū), 可以大力發(fā)展大豆種植業(yè)。

[1] 宮麗娟, 李宇光, 王萍, 等. 黑龍江省玉米氣候適宜度變化分析[J]. 吉林農(nóng)業(yè)科學(xué), 2012, 37(5): 75–80 GONG L J, LI Y G, WANG P, et al. Variation of climate suitability of maize in Heilongjiang Province[J]. Journal of Jilin Agricultural Sciences, 2012, 37(5): 75–80

[2] 陳炳東. 氣候變化對(duì)農(nóng)作物適宜度影響研究進(jìn)展[J]. 北京農(nóng)業(yè), 2012, (18): 153 CHEN B D. Impacts of climate change on crop suitability influence research progress[J]. Beijing Agriculture, 2012, (18): 153

[3] 劉實(shí), 王勇, 繆啟龍, 等. 近50年?yáng)|北地區(qū)熱量資源變化特征[J]. 應(yīng)用氣象學(xué)報(bào), 2010, 21(3): 266–278 LIU S, WANG Y, MIAO Q L, et al. Variation characteristics of thermal resources in Northeast China in recent 50 years[J]. Journal of Applied Meteorological Science, 2010, 21(3): 266–278

[4] 冶明珠, 郭建平, 袁彬, 等. 氣候變化背景下東北地區(qū)熱量資源及玉米溫度適宜度[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(10): 2786–2794 YE M Z, GUO J P, YUAN B, et al. Thermal resources and maize temperature suitability in Northeast China under climate change[J]. Chinese Journal of Applied Ecology, 2012, 23(10): 2786–2794

[5] 冶明珠, 郭建平, 蔣躍林, 等. 氣候變化對(duì)農(nóng)作物氣候適宜度影響研究進(jìn)展[J]. 安徽農(nóng)業(yè)科學(xué), 2011, 39(15): 9104–9105 YE M Z, GUO J P, JIANG Y L, et al. Research advances in impact of climate changes on crop climatic suitability[J]. Journal of Anhui Agricultural Sciences, 2011, 39(15): 9104–9105

[6] MATTHEWS R, WASSMANN R. Modelling the impacts of climate change and methane emission reductions on rice production: A review[J]. European Journal of Agronomy, 2003, 19(4): 573–598

[7] MESTRE-SANCHíS F, FEIJóO-BELLOM L. Climate change and its marginalizing effect on agriculture[J]. Ecological Economics, 2009, 68(3): 896–904

[8] 王麗霞, 任志遠(yuǎn). 山西省大同市農(nóng)業(yè)生態(tài)氣候適宜度評(píng)價(jià)[J].地理研究, 2007, 26(1): 53–59 WANG L X, REN Z Y. Appraisal and analysis of agriculture eco-climatic amenity in Datong City of Shanxi Province[J]. Geographical Research, 2007, 26(1): 53–59

[9] 魏瑞江, 張文宗, 康西言, 等. 河北省冬小麥氣候適宜度動(dòng)態(tài)模型的建立及應(yīng)用[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2007, 25(6): 5–9 WEI R J, ZHANG W Z, KANG X Y, et al. Application and establishment of climatic suitability dynamic model of winter wheat in Hebei Province[J]. Agricultural Research in the Arid Areas, 2007, 25(6): 5–9

[10] 俞芬, 千懷遂, 段海來(lái). 淮河流域水稻的氣候適宜度及其變化趨勢(shì)分析[J]. 地理科學(xué), 2008, 28(4): 537–542 YU F, QIAN H S, DUAN H L. Climate suitability of rice and its changes in Huaihe watershed[J]. Scientia Geographica Sinica, 2008, 28(4): 537–542

[11] 劉國(guó)成, 揚(yáng)長(zhǎng)保, 劉萬(wàn)崧, 等. 基于模糊數(shù)學(xué)的農(nóng)業(yè)氣候適宜度劃分研究及應(yīng)用[J]. 吉林農(nóng)業(yè)大學(xué)學(xué)報(bào), 2007, 29(4): 460–463 LIU G C, YANG C B, LIU W S, et al. Study and application of cropper climate suitability partition based on fuzzy mathematics[J]. Journal of Jilin Agricultural University, 2007, 29(4): 460–463

[12] 黃淑娥, 田俊, 吳慧峻. 江西省雙季水稻生長(zhǎng)季氣候適宜度評(píng)價(jià)分析[J]. 中國(guó)農(nóng)業(yè)氣象, 2012, 33(4): 527–533 HUANG S E, TIAN J, WU H J. Evaluation on climatic suitability during growth season of double rice in Jiangxi Province[J]. Chinese Journal of Agrometeorology, 2012, 33(4): 527–533

[13] 任玉玉, 千懷遂. 河南省棉花氣候適宜度變化趨勢(shì)分析[J]. 應(yīng)用氣象學(xué)報(bào), 2006, 17(1): 87–93 REN Y Y, QIAN H S. Climatic suitability of cotton and its changes in Henan Province[J]. Journal of Applied Meteorological Science, 2006, 17(1): 87–93

[14] 馬興祥, 鄧振鏞, 李棟梁, 等. 甘肅省春小麥生態(tài)氣候適宜度在適生種植區(qū)劃中的應(yīng)用[J]. 應(yīng)用氣象學(xué)報(bào), 2005, 16(6): 820–827 MA X X, DENG Z Y, LI D L, et al. Study on eco-climate applicability of spring wheat for condign planting division in Gansu Province[J]. Journal of Applied Meteorological Science, 2005, 16(6): 820–827

[15] 魏瑞江, 宋迎波, 王鑫. 基于氣候適宜度的玉米產(chǎn)量動(dòng)態(tài)預(yù)報(bào)方法[J]. 應(yīng)用氣象學(xué)報(bào), 2009, 20(5): 622–627 WEI R J, SONG Y B, WANG X. Method for dynamic forecast of corn yield based on climatic suitability[J]. Journal of Applied Meteorological Science, 2009, 20(5): 622–627

[16] 千懷遂, 焦士興, 趙峰. 全球氣候變化對(duì)作物氣候適宜性的影響——以河南省冬小麥為例[J]. 農(nóng)業(yè)現(xiàn)代化研究, 2004, 25(2): 106–110 QIAN H S, JIAO S X, ZHAO F. Influence of global climate change on climatic suitability of crop: A case study of winter wheat in Henan Province[J]. Research of Agricultural Modernization, 2004, 25(2): 106–110

[17] 卜坤, 王治良, 張樹文, 等. 三江平原大豆種植的土地適宜性評(píng)價(jià)[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2017, 25(3): 419–428 BU K, WANG Z L, ZHANG S W, et al. Evaluation of agricultural land suitability for soybean cultivation in the Sanjiang Plain, Northeast China[J]. Chinese Journal of Eco-Agriculture, 2017, 25(3): 419–428

[18] FISCHER G, SHAH M, VAN VELTHUIZEN H, et al. Agro-ecological zones assessments[M]//VERHEYE W H. Land Use, Land Cover and Soil Sciences. Oxford, UK: EOLSS Publishers, 2006

[19] 楊顯峰, 楊德光. 東北春大豆氣候適宜性指標(biāo)體系的建立研究初步[J]. 種子世界, 2009, (11): 36–38 YANG X F, YANG D G. Preliminary studies on establishing climatic adaptability index system for spring soybean in Northeast China[J]. Seed World, 2009, (11): 36–38

[20] 沈國(guó)強(qiáng). 適應(yīng)氣候變化的農(nóng)作物分布格局研究——以東北地區(qū)為例[D]. 北京: 中國(guó)科學(xué)院大學(xué)(中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所), 2017: 52–57 SHEN G Q. Research on crops pattern for adapting to climate change — A case study in Northeast China[D]. Beijing: University of Chinese Academy of Science (Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences), 2017: 52–57

[21] 李昊宇, 王建林, 鄭昌玲, 等. 氣候適宜度在華北冬小麥發(fā)育期預(yù)報(bào)中的應(yīng)用[J]. 氣象, 2012, 38(12): 1554–1559 LI H Y, WANG J L, ZHENG C L, et al. The development period prediction of winter wheat based on climatic suitability in North China[J]. Meteorological Monthly, 2012, 38(12): 1554–1559

[22] 王鶯, 王勁松, 姚玉璧, 等. 甘肅省河?xùn)|地區(qū)氣象干旱災(zāi)害風(fēng)險(xiǎn)評(píng)估與區(qū)劃[J]. 中國(guó)沙漠, 2014, 34(4): 1115–1124 WANG Y, WANG J S, YAO Y B, et al. Assessment and regionalization on meteorological drought disaster risk in the Hedong Area of Gansu Province, China[J]. Journal of Desert Research, 2014, 34(4): 1115–1124

[23] 張麗娟, 陳紅, 張金峰, 等. 黑龍江省雷暴災(zāi)害損失風(fēng)險(xiǎn)評(píng)估與區(qū)劃[J]. 自然災(zāi)害學(xué)報(bào), 2011, 20(6): 117–123 ZHANG L J, CHEN H, ZHANG J F, et al. Risk assessment and zoning of thunderstorm loss in Heilongjiang Province[J]. Journal of Natural Disasters, 2011, 20(6): 117–123

[24] 李文亮, 張麗娟, 張冬有. 黑龍江省低溫冷害風(fēng)險(xiǎn)評(píng)估與區(qū)劃研究[J]. 干旱區(qū)資源與環(huán)境, 2009, 23(10): 149–153 LI W L, ZHANG L J, ZHANG D Y. Study on risk assessment and zoning of low temperature and cold damage in Heilongjiang Province[J]. Journal of Arid Land Resources and Environment, 2009, 23(10): 149–153

[25] 王希平, 趙慧穎. 內(nèi)蒙古呼倫貝爾市林牧農(nóng)業(yè)氣候資源與區(qū)劃[M]. 北京: 氣象出版社, 2006 WANG X P, ZHAO H Y. Climatic Resources and Regionalization of Forestry and Pastoral Farming in Hulun Buir[M]. Beijing: China Meteorological Press, 2006

Climate suitability grading and planting zoning of soybean in Northeast Inner Mongolia*

WANG Yanping1, YIN Xiuxia1, ZHANG Fang2, ZHANG Xuming1, WU Changshun3

(1. Hulun Buir Meteorological Office, Hailar 021008, China; 2. Zhalantun Buir Meteorological Office, Zhalantun 162650, China; 3. Xinba’erhuyouqi Buir Meteorological Office, Alatan’emole 021300, China)

In northeast of Inner Mongolia, soybean production is faced with two problems. The first is that soybean planting area is zoned with ≥10 ℃accumulative temperature without non clear division between soybean and maize region. The second is the use efficiencies of heat and water in the area are low in the area. In order to efficiently make use of water and heat resources in soybean cultivation and promote soybean production in cool area, we used meteorological observation data and soybean crop growing phrases data in years from 1987 to 2016 in Hulun Buir, Xing’an League in Northeast Inner Mongolia, and surrounding areas in Heilongjiang Province to analyze the climate suitability of soybean in Hulun Buir. In the study, temperature suitability, precipitation suitability of soybean were calculated, respectively, with the “three basis points” of temperature and water demand in each growth phase of soybean, while the climate suitability was calculated with the comprehensive hydrothermal suitability model. Furthermore, planting zoning of soybean cultivation in Hulun Buir was conducted according the comprehensive hydrothermal suitability level in GIS environment. The results showed that for the temperature suitability and comprehensive hydrothermal suitability, ≥0.6, 0.4-0.6 and ≤ 0.4 were classified as most suitable, more suitable and unsuitable. For precipitation suitability, ≥ 0.7 was most suitable, 0.6-0.7 more suitable and ≤ 0.6 unsuitable. The most suitable areas for planting soybean were in the southeast of Daxing’anling farming region, more suitable areas at the edge of the planting region in eastern Daxing’anling, then unsuitable areas were in the north woodland and western pastoral of Hulun Buir. Agroforestry ecotone in the south of forest region where accumulated temperature exceeded 10 ℃ was 1 900-2 100 ℃, which was not suitable for cultivating corn, but could be developed vigorously for cultivating soybean. This area was considered for increased soybean acreage during structural crop adjustment. The results provided reference for soybean production in the cool area in Northeast China.

Northeastern Inner Mongolia; Climate suitability; Evaluation index; Soybean cultivation; Regionalization method

, WANG Yanping, E-mail: fei_ping@sohu.com

Dec. 4, 2017;

Mar. 19, 2018

P49

A

1671-3990(2018)07-0948-10

10.13930/j.cnki.cjea.171128

2017-12-04

2018-03-19

* This study was supported by the Special Climate Change Project of the China Meteorological Administration (CCSF201606), China Meteorological Administration’s Core Business Development Project (CMAHX20160205) and China Meteorological Administration Lanzhou Institute of Drought Meteorological Science Research Foundation (IAM201706).

* 中國(guó)氣象局氣候變化專項(xiàng)(CCSF201606)、中國(guó)氣象局核心業(yè)務(wù)發(fā)展專項(xiàng)(CMAHX20160205)和中國(guó)氣象局蘭州干旱氣象研究所干旱氣象科學(xué)研究基金(IAM201706)資助

王彥平, 主要研究方向?yàn)檗r(nóng)牧業(yè)氣象服務(wù)與研究。E-mail: fei_ping@sohu.com

王彥平, 陰秀霞, 張昉, 張煦明, 烏長(zhǎng)順. 內(nèi)蒙古東北部大豆氣候適宜度等級(jí)及種植區(qū)劃研究[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2018, 26(7): 948-957

WANG Y P, YIN X X, ZHANG F, ZHANG X M, WU C S. Climate suitability grading and planting zoning of soybean in Northeast Inner Mongolia[J]. Chinese Journal of Eco-Agriculture, 2018, 26(7): 948-957