吳玉潔,霍治國(guó),姜會(huì)飛,李建勇,李建平,吳 迪
(1.吉林省氣象科學(xué)研究所/長(zhǎng)白山氣象與氣候變化吉林省重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)春 130062;2.中國(guó)氣象科學(xué)研究院,北京 100081;3.南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,南京 210044;4.中國(guó)農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,北京 100193;5.河南省平輿縣氣象局,駐馬店 463400)
江淮地區(qū)位于我國(guó)東部中緯度地區(qū),受季風(fēng)影響,易出現(xiàn)降水極端天氣事件,是澇漬災(zāi)害重災(zāi)區(qū)。玉米作為區(qū)域主要糧食作物,其生長(zhǎng)季同汛期時(shí)間重疊,故受澇漬災(zāi)害嚴(yán)重。1978年7月24—25日,碭山縣全縣積水農(nóng)田達(dá)2.53萬(wàn)hm2,造成玉米倒伏;1979年7月10—30日,鹽城市積水農(nóng)田達(dá)10萬(wàn)hm2,伴隨臺(tái)風(fēng)影響,陸集鄉(xiāng)600hm2玉米基本倒伏,斷稈率達(dá)30%;1991年,雎寧縣全縣2.7萬(wàn)hm2夏玉米受澇漬危害面積在1萬(wàn)hm2以上,占玉米種植面積的39.3%[1-3]。因此,研究揭示江淮地區(qū)玉米澇漬災(zāi)害時(shí)空變化特征對(duì)科學(xué)開展區(qū)域防洪減災(zāi)具有重要意義。
目前,國(guó)內(nèi)外用于評(píng)估澇漬等級(jí)的指標(biāo)模型形式多樣,農(nóng)業(yè)氣象領(lǐng)域中,始終以氣候指標(biāo)、農(nóng)業(yè)氣象指標(biāo)為研究核心。氣候指標(biāo)如降水距平百分率、降水標(biāo)準(zhǔn)差、SPI指數(shù)、K指數(shù)和Z指數(shù)等較為常用[4,5],基于降水概率、連續(xù)日降水量、暴雨洪災(zāi)次數(shù)、面雨量等指標(biāo)構(gòu)建及劃分方法也較為普遍[6-9]。此類指標(biāo)可作為區(qū)域澇漬災(zāi)害的氣候判別標(biāo)準(zhǔn),但因未考慮土壤和作物因素,用于作物澇漬災(zāi)害的評(píng)估誤差較大。農(nóng)業(yè)氣象指標(biāo)則將氣象要素同作物減產(chǎn)及受災(zāi)情況相結(jié)合,以作物受害的氣象要素影響機(jī)制為建模依據(jù),常用指標(biāo)如作物水分盈虧指數(shù)(相對(duì)濕潤(rùn)度指數(shù))、陰濕系數(shù)、干濕指數(shù)等[10-13],或是以上述指標(biāo)為原型改進(jìn)的指標(biāo)[14,15]。此外,基于過程雨量的作物澇漬災(zāi)害指標(biāo)研究近些年也較為普遍[16,17]。
玉米是一種需水量大但不耐澇的作物,且各階段受澇的敏感性不同。氣象因素中除強(qiáng)降水外,連陰雨天氣同樣會(huì)促使玉米發(fā)生澇漬災(zāi)害。江淮地區(qū)玉米種植季降水顯著,且播種前期常伴有梅雨天氣易導(dǎo)致種子萌發(fā)期間發(fā)生霉變,造成澇漬災(zāi)害的發(fā)生。基于玉米澇漬受害過程及區(qū)域降水特征,文章以作物水分盈余指數(shù)為原型,引入玉米播種前底墑,利用層次分析法[18]確定各階段受澇對(duì)減產(chǎn)率的影響權(quán)重,構(gòu)建適用于江淮地區(qū)玉米澇漬等級(jí)評(píng)估的全生育期綜合指數(shù)。此方法在玉米澇漬指標(biāo)構(gòu)建研究中鮮有報(bào)道,且通過對(duì)江淮地區(qū)春、夏玉米澇漬等級(jí)閾值的精細(xì)劃分,實(shí)現(xiàn)對(duì)玉米澇漬災(zāi)害的精細(xì)化、精準(zhǔn)化研究,可為各級(jí)政府及生產(chǎn)部門指揮玉米生產(chǎn)和防洪抗災(zāi)工作提供更有效的科學(xué)依據(jù)和指導(dǎo)意見。
以蘇皖長(zhǎng)江以北地區(qū)為研究區(qū)域,選取區(qū)域數(shù)據(jù)資料較完整的18個(gè)農(nóng)業(yè)氣象站(圖1)及1981—2010年對(duì)應(yīng)地面氣象觀測(cè)站的數(shù)據(jù)資料進(jìn)行分析。其中,所用農(nóng)業(yè)氣象資料包括:玉米單產(chǎn)、生育期數(shù)據(jù)及災(zāi)情記錄,來源于中國(guó)氣象局氣象檔案館存檔的《作物生長(zhǎng)發(fā)育狀況記錄年報(bào)表(農(nóng)氣表1)》,部分缺失的玉米單產(chǎn)數(shù)據(jù)由“中國(guó)種植業(yè)信息網(wǎng)—縣級(jí)農(nóng)作物數(shù)據(jù)庫(kù)”補(bǔ)充,極少數(shù)缺失的生育期數(shù)據(jù)利用已有站點(diǎn)數(shù)據(jù)插值獲得。所用逐日氣象要素?cái)?shù)據(jù)包括:降水量、平均氣溫、最高氣溫、最低氣溫,來源于國(guó)家級(jí)地面氣象站基本氣象要素日值數(shù)據(jù)集(V3.0),對(duì)極少數(shù)氣象要素的缺測(cè)值根據(jù)歷年同期平均值進(jìn)行補(bǔ)充。
圖1 農(nóng)業(yè)氣象站地理方位分布
減產(chǎn)率是氣象產(chǎn)量相對(duì)化處理后的負(fù)值結(jié)果[19],公式為:
(1)
Y為玉米氣象產(chǎn)量波動(dòng)的百分率,取其中由典型澇漬年引起的負(fù)值序列作為最終減產(chǎn)率序列(%);ya為作物實(shí)際單產(chǎn)(kg/hm2);yt為趨勢(shì)單產(chǎn)(kg/hm2);t為年份(年)。根據(jù)已有研究常用劃分標(biāo)準(zhǔn),認(rèn)為Y≤-5%為減產(chǎn)年,-10% 以作物水分盈余指數(shù)(CWSI)為原型進(jìn)行改進(jìn)[13],在降水量的計(jì)算項(xiàng)中,考慮播種前過度降水對(duì)玉米種子萌發(fā)的不良影響?;谟衩撞煌A段澇漬敏感性差異,計(jì)算玉米各生育時(shí)期的水分盈余指數(shù)(WSIMj)并配以權(quán)重系數(shù),得到玉米綜合水分盈余指數(shù)(CWSIM)。公式為: (2) (3) 式(2)~(3)中,WSIMj為玉米在j階段的水分盈余指數(shù),j取1、2、3、4,分別表示播種—出苗期、出苗—拔節(jié)期、拔節(jié)—吐絲期、吐絲—成熟期;Rn為播種前n天總降水量(mm),即底墑水(Rn僅在第1階段中計(jì)算);Rj為玉米在j階段的總降水量(mm);ETCj為玉米在j階段實(shí)際蒸散總量(mm);CWSIM為玉米綜合水分盈余指數(shù);wj為WSIMj的影響權(quán)重。 農(nóng)田實(shí)際蒸散量常用計(jì)算公式為: ETc=Kc×ET0 (4) (5) 式(4)~(5)中,ETc為實(shí)際蒸散量(mm);ET0為潛在蒸散量[20](mm);Kc為作物系數(shù),取初始生長(zhǎng)期Kc=0.3,生育中期Kc=1.2,收獲期Kc=0.6[21];Tmean、Tmax、Tmin分別為日平均、最高、最低溫度(°C);Ra為大氣頂層輻射[22](MJ/(m2.d));λ為水汽化潛熱,通常取2 450(MJ/kg)。 WSIMj的權(quán)重系數(shù)由層次分析法確定[18],其中,構(gòu)造判斷矩陣中的比較值可由現(xiàn)有玉米淹水試驗(yàn)結(jié)果[23]確定,因苗期受澇時(shí)人工補(bǔ)苗對(duì)結(jié)果的干擾較大,因此該階段的權(quán)重設(shè)定為最小。 利用SPSS軟件,對(duì)CWSIM和其對(duì)應(yīng)減產(chǎn)率序列Y做線性回歸分析,對(duì)通過顯著性檢驗(yàn)的回歸方程,根據(jù)Y的閾值劃分澇漬等級(jí)的標(biāo)準(zhǔn),可以確定對(duì)應(yīng)CWSIM值劃分澇漬等級(jí)的標(biāo)準(zhǔn)。 表1 啟東站春玉米典型澇漬年Y同WSIM相關(guān)性 分別以播種前不同時(shí)間尺度的降水為底墑,并對(duì)典型澇漬年份玉米減產(chǎn)率(Y)同其全生育期/分生育階段水分盈余指數(shù)(WSIM,表示未區(qū)分玉米各生育階段澇漬敏感性的指數(shù))做相關(guān)性分析,得到如下結(jié)果。 表1中,啟東站春玉米Y同WSIM1~WSIM7均呈顯著相關(guān),相關(guān)性排序?yàn)椋篧SIM5>WSIM4>WSIM6>WSIM7>WSIM3>WSIM2>WSIM1。WSIM1最小,說明考慮播種前底墑的全生育期指標(biāo)評(píng)價(jià)效果較不考慮底墑的更好;WSIM8~WSIM11均未通過顯著性檢驗(yàn),說明以玉米某個(gè)關(guān)鍵生育期計(jì)算的指標(biāo),用于評(píng)價(jià)玉米因受澇造成的減產(chǎn)具有其局限性。 表2中,沭陽(yáng)站夏玉米Y同WSIM1~WSIM11相關(guān)性排序有:WSIM9>WSIM5>WSIM3>WSIM1>WSIM4>WSIM2>WSIM7>WSIM6>WSIM11>WSIM8>WSIM10,除WSIM9外,其余均未通過相關(guān)性檢驗(yàn),說明沭陽(yáng)站夏玉米因澇漬害引起的減產(chǎn)主要是受出苗—拔節(jié)期的澇漬害影響,而以常規(guī)的不區(qū)分階段澇漬敏感性計(jì)算的全生育期指數(shù)無(wú)法說明此例情況。 WSIM8同對(duì)應(yīng)Y成正相關(guān)關(guān)系,反觀WSIM8的數(shù)值普遍高于WSIM9,而此階段澇漬敏感性較出苗—拔節(jié)期更大,說明此階段指數(shù)同Y的相關(guān)關(guān)系是異常的。通過對(duì)災(zāi)情資料及實(shí)際農(nóng)民種植行為的分析,這種異常是由苗期強(qiáng)降水過后人為補(bǔ)苗,干擾了該期間指數(shù)同Y的相關(guān)關(guān)系造成的。所以,在構(gòu)建指標(biāo)時(shí)應(yīng)考慮苗期人為因素對(duì)評(píng)估結(jié)果造成的影響。 表2 沭陽(yáng)站夏玉米典型澇漬年Y同WSIM相關(guān)性 根據(jù)上述結(jié)論,對(duì)全生育期玉米澇漬指標(biāo)進(jìn)行分生育階段計(jì)算,確定播種—出苗、出苗—拔節(jié)、拔節(jié)—吐絲、吐絲—成熟期4個(gè)階段的權(quán)重系數(shù)依次為0.06、0.57、0.26、0.11,并結(jié)合不同時(shí)間尺度底墑,計(jì)算玉米綜合水分盈余指數(shù)(CWSIM,表示考慮玉米各生育階段澇漬敏感性的綜合指數(shù)),將其同Y做相關(guān)性分析,得到結(jié)果如下。 表3中,啟東站春玉米所有指標(biāo)均通過顯著性檢驗(yàn),且以CWSIM4即播種前15d降水為底墑?dòng)?jì)算的指標(biāo)值同Y的相關(guān)性最大,相關(guān)系數(shù)為-0.996,通過了極顯著檢驗(yàn)。 表3 啟東站春玉米典型澇漬年Y同CWSIM相關(guān)性 表4中,沭陽(yáng)站夏玉米在改進(jìn)指標(biāo)后,原未通過顯著性檢驗(yàn)的指標(biāo)均通過了檢驗(yàn),CWSIM1、CWSIM5、CWSIM7同Y的相關(guān)性大于僅考慮出苗—拔節(jié)期的情況,且以CWSIM7即播種前30d降水為底墑?dòng)?jì)算的指標(biāo)值同Y的相關(guān)性最大,相關(guān)系數(shù)為-0.937,通過了極顯著檢驗(yàn)。 表4 沭陽(yáng)站夏玉米典型澇漬年Y同CWSIM相關(guān)性 綜合上述結(jié)論有:玉米區(qū)分階段澇漬敏感性構(gòu)建的綜合指標(biāo)更具代表性,對(duì)春、夏玉米均適用,且以考慮底墑的全生育期指標(biāo)優(yōu)于僅考慮玉米關(guān)鍵期和不考慮底墑的全生育期指標(biāo)。因此,采用改進(jìn)的玉米綜合水分盈余指數(shù)(CWSIM)作為評(píng)估江淮地區(qū)玉米澇漬等級(jí)的指標(biāo),且春玉米以播種前15d降水作為底墑,夏玉米以播種前30d降水作為底墑。 將啟東站春玉米、沭陽(yáng)站夏玉米的CWSIM同其對(duì)應(yīng)Y做線性回歸分析,分別得到對(duì)應(yīng)關(guān)系式為:Y=-23.34CWSIM+13.69,R2=0.99;Y=-14.30CWSIM+10.64,R2=0.93。則CWSIM劃分澇漬等級(jí)的閾值結(jié)果如表5。 表5 CWSIM澇漬等級(jí)劃分閾值 表5中,對(duì)比CWSIM閾值,夏玉米對(duì)澇漬的耐受程度高于春玉米。造成該結(jié)果的原因除與玉米品種有關(guān)外,還同底墑水引入的時(shí)間尺度,季節(jié)溫度和日照差異導(dǎo)致玉米需水量的不同等因素有關(guān)。因此,對(duì)同種作物不同季品種的澇漬等級(jí)評(píng)估指標(biāo)的分級(jí)應(yīng)有所區(qū)分。 以上述單站結(jié)果為參考標(biāo)準(zhǔn),對(duì)江淮地區(qū)其他已有記錄的站點(diǎn)進(jìn)行驗(yàn)證(表6),11個(gè)驗(yàn)證樣本均達(dá)標(biāo),說明以單站指標(biāo)為參照標(biāo)準(zhǔn)對(duì)于江淮地區(qū)玉米受澇情況評(píng)價(jià)基本適用。值得注意的是,典型澇漬年,當(dāng)玉米減產(chǎn)率在15%以上達(dá)到重澇級(jí)別時(shí),各站點(diǎn)的指數(shù)同Y的相關(guān)性已不再是簡(jiǎn)單的線性關(guān)系,因?yàn)闃O度降水過后的人為管理對(duì)最終玉米減產(chǎn)的影響是難以預(yù)估的,且當(dāng)降水遠(yuǎn)大于玉米需水時(shí)澇漬災(zāi)害同玉米減產(chǎn)的相關(guān)關(guān)系也已不僅僅是線性相關(guān),但這并不影響澇漬指數(shù)對(duì)玉米澇漬級(jí)別的判斷結(jié)果。 表6 江淮地區(qū)CWSIM驗(yàn)證 利用求算出的已知站點(diǎn)指數(shù),對(duì)區(qū)域進(jìn)行空間插值,得到其春、夏玉米各發(fā)育階段(WSIM)和全生育期(CWSIM)的水分盈余指數(shù)(由于苗前人為因素干擾較大,因此未對(duì)播種—出苗階段的指數(shù)進(jìn)行分析),圖2中3種顏色由淺至深分別表示低值區(qū)、中值區(qū)和高值區(qū),主要結(jié)果如下。 圖2中,春玉米出苗—拔節(jié)、拔節(jié)—吐絲期兩個(gè)階段的WSIM呈由北向南的緯向增加分布,且安徽省西南臨江一帶指數(shù)最大;吐絲—成熟期,安徽省西南臨江一帶、江蘇省淮安市東部、揚(yáng)州市及其周邊區(qū)域?yàn)楦咧祬^(qū),安徽省北部及中部的中心區(qū)域、江蘇省東南邊緣為低值區(qū),其余為中值區(qū)。 (a.出苗—拔節(jié)期;b.拔節(jié)—吐絲期;c.吐絲—成熟期(1981—2010年年均值))圖2 江淮地區(qū)春玉米各階段WSI空間分布 圖3中,夏玉米出苗—拔節(jié)期的WSIM呈由北向南的緯向增加分布,安徽省臨江一帶為高值區(qū);拔節(jié)—吐絲期,江蘇省東北部為高值區(qū),安徽省中部及臨江一帶為低值區(qū),其余為中值區(qū);吐絲—成熟期指數(shù)呈由西至東的經(jīng)向增加分布,江蘇省東部臨海一帶為高值區(qū)。夏玉米各階段WSIM的分布結(jié)果,同前人研究中江淮地區(qū)夏玉米受澇情況分布的結(jié)果大致相同[24]。 (a.出苗—拔節(jié)期;b.拔節(jié)—吐絲期;c.吐絲—成熟期(1981—2010年年均值))圖3 江淮地區(qū)夏玉米各階段WSI空間分布 圖4中,對(duì)1981—2010年30年春玉米CWSIM每10年求1次平均值,分別得到a、b、c所示結(jié)果,d為1981—2010年30年的平均值。從空間看(4d),春玉米CWSIM呈由北向南的緯向增加分布,分布規(guī)律對(duì)a、b、c同樣適用,安徽省西南部沿江一帶為高值區(qū);從時(shí)間看,整體指數(shù)在不同時(shí)期的大小排序?yàn)?991—2000年>1981—1990年>2001—2010年,說明春玉米在30年受澇情況是先加重后減輕。 (a.1981—1990年年均值;b.1991—2000年年均值;c.2001—2010年年均值;d.1981—2010年年均值)圖4 江淮地區(qū)春玉米CWSIM時(shí)空分布 同理,得到夏玉米CWSIM時(shí)空分布圖5。從空間看(圖5.d),夏玉米CWSIM由北向南呈近輻射狀條帶分布(分布規(guī)律對(duì)a、b、c不盡相同),安徽省北部偏東一帶為低值區(qū),安徽省沿江及江蘇省東南邊緣為高值區(qū),其余為中值區(qū)。從時(shí)間看,整體指數(shù)在不同時(shí)期的大小排序?yàn)?991—2000年>2001—2010年>1981—1990年。圖5.a,安徽省沿江一帶為高值區(qū),安徽省北部、江蘇省西北部和沿江中段地帶為低值區(qū),其余為中值區(qū);圖5.b,分布規(guī)律整體同圖5.a相近,圖5.b較圖5.a,安徽省沿江帶的高值區(qū)面積明顯增大,低值區(qū)面積減小且向安徽省和江蘇省的北部中心區(qū)域轉(zhuǎn)移,江蘇省東南部指數(shù)增加;圖5.c,指數(shù)分布規(guī)律與圖5.a、圖5.b均不同,江淮區(qū)域被中值區(qū)整體覆蓋,說明安徽省沿江一帶澇情減輕,但江淮北部地區(qū)加重。 綜合春、夏玉米的情況,安徽省沿江一帶是玉米澇漬害發(fā)生的重災(zāi)區(qū),且1991—2000年較其前后10年澇情最重。 (a.1981—1990年年均值;b.1991—2000年年均值;c.2001—2010年年均值;d.1981—2010年年均值)圖5 江淮地區(qū)夏玉米CWSIM時(shí)空分布 (1)春玉米的輕、中、重澇CWSIM閾值依次為:0.80≤CWSIM<1.01,1.01≤CWSIM<1.23、CWSIM≥1.23。夏玉米依次為:1.09≤CWSIM<1.44,1.44≤CWSIM<1.79,CWSIM≥1.79。 (2)各階段WSIM分布規(guī)律,春玉米為:出苗—拔節(jié)期、拔節(jié)—吐絲期呈由北向南的緯向增加分布;吐絲—成熟期高值區(qū)位于安徽省西南臨江一帶、江蘇省淮安市東部、揚(yáng)州市及其周邊區(qū)域,低值區(qū)位于安徽省北部及中部的中心區(qū)域、江蘇省東南邊緣,其余為中值區(qū)。夏玉米為:出苗—拔節(jié)期呈由北向南的緯向增加分布;拔節(jié)—吐絲期高值區(qū)位于江蘇省東北部,低值區(qū)位于安徽省中部及臨江一帶,其余為中值區(qū);吐絲—成熟期指數(shù)呈由西向東的經(jīng)向增加分布。 (3)春玉米近30年CWSIM時(shí)空分布規(guī)律為:空間上,由北向南呈緯向增加分布。時(shí)間上,整體指數(shù)在不同時(shí)期的排序?yàn)?991—2000年>1981—1990年>2001—2010年;夏玉米的規(guī)律為:空間上,由北向南呈近輻射狀條帶分布。時(shí)間上,整體指數(shù)在不同時(shí)期的排序?yàn)?991—2000年>2001—2010年>1981—1990年。綜合說明,安徽省沿江一帶是玉米澇漬害發(fā)生的重災(zāi)區(qū),且1991—2000年較其前后10年澇情最重。 (1)受資料數(shù)量與質(zhì)量的限制,區(qū)域指標(biāo)研究常以單站結(jié)果作為區(qū)域結(jié)果。但實(shí)地考察證明,即便是鄰近區(qū)域,指標(biāo)的差異也可能非常明顯。這與土壤特性、地勢(shì)、地下水埋深、局地小氣候、人為影響等眾多因素有關(guān)。實(shí)現(xiàn)區(qū)域化研究仍需依靠充足而準(zhǔn)確的數(shù)據(jù)資源。進(jìn)一步完善農(nóng)業(yè)氣象站布局,增設(shè)作物觀察點(diǎn)和觀測(cè)指標(biāo)以提供充足的觀測(cè)數(shù)據(jù),對(duì)于農(nóng)業(yè)氣象的精準(zhǔn)化、精細(xì)化研究是必不可少的。而伴隨著遙感技術(shù)的日益成熟,衛(wèi)星遙感監(jiān)測(cè)及反演技術(shù)也已成為協(xié)助研究的關(guān)鍵方法,可有效提高觀測(cè)效率。 (2)統(tǒng)計(jì)生育期的農(nóng)業(yè)氣象站數(shù)量的局限給區(qū)域分析帶來困難,且單個(gè)站點(diǎn)的統(tǒng)計(jì)數(shù)據(jù)難以實(shí)現(xiàn)連續(xù)長(zhǎng)時(shí)間記錄。構(gòu)建長(zhǎng)時(shí)序列時(shí)需要人為處理補(bǔ)全,以多年均值替代或以單站、鄰近站數(shù)據(jù)替代等方法較為常用,使結(jié)果存在較大誤差,而運(yùn)用積溫理論推測(cè)作物發(fā)育期可有效提高準(zhǔn)確率,減小誤差結(jié)果[25]。 (3)該研究以CWSIM同Y的相關(guān)關(guān)系進(jìn)行分析,可對(duì)全生育期的受澇程度進(jìn)行定量判斷。對(duì)各生育階段WSIM僅對(duì)比分析了指數(shù)的分布情況,一定程度上可反映出玉米不同階段的受澇程度,但若做出定量判斷,還需結(jié)合各階段具體受澇記錄進(jìn)行分析。后續(xù)研究可從各階段玉米受澇次數(shù)、受澇等級(jí)的統(tǒng)計(jì)結(jié)果與階段指數(shù)的相關(guān)關(guān)系方面入手,進(jìn)一步實(shí)現(xiàn)對(duì)玉米澇漬等級(jí)評(píng)估的精細(xì)化研究。1.3 玉米澇漬等級(jí)評(píng)估指標(biāo)構(gòu)建
1.4 農(nóng)田實(shí)際蒸散量計(jì)算
1.5 權(quán)重系數(shù)計(jì)算
1.6 玉米澇漬等級(jí)評(píng)估指標(biāo)閾值劃分
2 結(jié)果與分析
2.1 指標(biāo)構(gòu)建方法的選定
2.2 指標(biāo)閾值劃分及驗(yàn)證
2.3 江淮地區(qū)玉米水分盈余指數(shù)時(shí)空分布特征
3 結(jié)論與討論
3.1 主要結(jié)論
3.2 問題討論
中國(guó)農(nóng)業(yè)資源與區(qū)劃2019年11期