利用黑土層厚度評價東北黑土區(qū)土壤生產(chǎn)力
——以鶴北小流域?yàn)槔?/h1>
2014-09-13 10:06:56段興武
水土保持通報 2014年1期
李 驁， 段興武
(1．華中農(nóng)業(yè)大學(xué) 430070， 湖北 武漢； 2．云南大學(xué) 亞洲國際河流研究中心， 云南 昆明 650091)
東北黑土區(qū)是中國重要的糧食生產(chǎn)基地。據(jù)農(nóng)業(yè)部農(nóng)業(yè)信息網(wǎng)公布的統(tǒng)計資料，2011年東北三省的糧食總產(chǎn)量為1.07×1011kg，約占全國總產(chǎn)的18.86%，其中玉米產(chǎn)量占全國的33.06%，大豆占46.78%，人均糧食占有量是全國平均水平的2.4倍，在保障國家糧食安全問題上起著舉足輕重的作用，被譽(yù)為中國糧食生產(chǎn)的“穩(wěn)壓器”[1-2]。
黑土養(yǎng)分含量豐富、結(jié)構(gòu)良好、土壤肥力極高[3]。自從20世紀(jì)初大面積開墾以來，黑土土壤理化性質(zhì)發(fā)生了顯著變化[4-8]，土壤生產(chǎn)力顯著降低[7-9]，引起黑土土壤生產(chǎn)力變化的原因比較復(fù)雜，可主要分為施肥和侵蝕兩種。為提高作物產(chǎn)量，自20世紀(jì)70年代末80年代初開始，中國開始大量施用化工肥料[10]，化工肥料的施用能顯著提高作物產(chǎn)量[11]，同時也能改變土壤的理化性質(zhì)[12-13]，從而影響到作物品質(zhì)的提高和土壤肥力的可持續(xù)利用[14]。另一方面，在高強(qiáng)度的人為耕作下，黑土區(qū)土壤侵蝕嚴(yán)重[1-4]。據(jù)水利部松遼水利委員會2002年的統(tǒng)計報告，該區(qū)水土流失面積占全區(qū)總面積的27%，發(fā)育大型切溝2.5×104條余，吞噬耕地面積約3.9 hm2[5]。每年僅黑龍江省流失的土壤就變化在2.0×108～3.0×108m3，流失的養(yǎng)分折合標(biāo)準(zhǔn)化肥5.0×109～6.0×109kg，相當(dāng)于全省化肥施用總量[6]，嚴(yán)重的侵蝕造成土壤養(yǎng)分損失、土壤質(zhì)地粗化、持水能力降低，直接導(dǎo)致土壤生產(chǎn)力退化[13]。在這樣的前提下，研究簡便易行的土壤生產(chǎn)力評價指標(biāo)可為黑土資源的合理利用和保護(hù)提供理論依據(jù)和方法參考。
土壤生產(chǎn)力評價方法可分為定性和定量兩種[15]。定性評價是通過對田間土壤的實(shí)地調(diào)查和診斷，對土壤宏觀性狀進(jìn)行定性描述，多帶有經(jīng)驗(yàn)性，很難量化[16]。定量評價主要有3種：一是單指標(biāo)評價；二是多指標(biāo)表土綜合評價，即選取表層土壤多個理化性質(zhì)指標(biāo)評價表土生產(chǎn)力；三是多指標(biāo)土體綜合評價，即將整個土體看作是作物生長的環(huán)境，選擇各個土層的土壤性質(zhì)指標(biāo)，將不同深度土層賦予不同的權(quán)重，綜合、立體地評價土壤生產(chǎn)力，比較有代表性的是土壤生產(chǎn)力指數(shù)模型(productivity index, PI)[17]，該模型在國內(nèi)外獲得了廣泛應(yīng)用[18-20]。段興武等[15]在該模型基礎(chǔ)上，根據(jù)東北黑土土壤理化性質(zhì)和作物產(chǎn)量調(diào)查結(jié)果，構(gòu)建了適用于東北黑土區(qū)的修正的生產(chǎn)力指數(shù)模型(modified productivity Index, MPI)并得到多次驗(yàn)證[21-22]。但該模型仍然需要采集剖面土壤樣本、測定多項(xiàng)土壤理化指標(biāo)，不利于野外快速和大范圍土壤生產(chǎn)力調(diào)查與評價。表層黑土養(yǎng)分含量豐富、土壤結(jié)構(gòu)良好，表層黑土厚度(簡稱黑土層厚度)便于野外調(diào)查，對黑土土壤生產(chǎn)力貢獻(xiàn)較大[3]，被認(rèn)為是衡量黑土退化程度的重要指標(biāo)[7-8,23]。但時至今日，鮮有探討黑土層厚度與土壤生產(chǎn)力的研究案例。
本文以位于東北典型黑土區(qū)北端的鶴北小流域?yàn)榘咐齾^(qū)，通過調(diào)查和分析流域內(nèi)不同等級土地的土壤厚度及土壤理化性質(zhì)，構(gòu)建了黑土層厚度與土壤生產(chǎn)力間的關(guān)系，并據(jù)此評價典型地塊的土壤生產(chǎn)力空間特征及其影響因素。
1 資料與方法
1.1 研究區(qū)概況
研究區(qū)位于黑龍江省農(nóng)墾總局九三分局鶴山農(nóng)場境內(nèi)的鶴北小流域，地理坐標(biāo)為48°59′0″—49°3′0″N，125°16′0″—125°121′0″E，屬于寒溫帶大陸性季風(fēng)氣候，冬季漫長，夏季短促，秋季涼爽初霜早，春季風(fēng)大，年平均氣溫-1.5～0.4 ℃，歷年最高氣溫37.9 ℃，最低氣溫-46.5 ℃，年降水量在400～500 mm，集中于6—9月。該區(qū)位于小興安嶺南麓，是小興安嶺向嫩江平原過渡地帶，屬漫崗丘陵區(qū)，坡長一般800～1 500 m，坡度一般1°～4°。自然植被屬于溫帶半濕潤草甸草原，土壤屬中國發(fā)生學(xué)分類中的黑土，屬中國土壤系統(tǒng)分類中的均腐土[24]。
鶴北小流域的面積為28 km2，大部分已開墾為農(nóng)地，且開墾歷史都在60 a以上。開墾之初，通過種植人工林帶將各個地塊隔開，且無論是起壟或行播均沿林帶方向?qū)嵭袡C(jī)械化作業(yè)。流域內(nèi)共包括41個地塊，種植的農(nóng)作物主要有大豆、小麥和玉米，尤以大豆為主。農(nóng)業(yè)機(jī)械化程度高，整地、播種、噴藥、施肥、除草和收割等全部實(shí)行機(jī)械化作業(yè)。
1.2 樣本采集
本研究根據(jù)地形狀況，選擇流域內(nèi)20個典型地塊調(diào)查土壤厚度和剖面土壤理化性質(zhì)。
(1)黑土厚度調(diào)查。從調(diào)查地塊的坡頂開始每10 m高差布設(shè)一個厚度樣點(diǎn)，共計布設(shè)厚度樣點(diǎn)122個。利用半圓鉆(長1 m，可加延長桿，直徑4 cm)采集土樣，將樣本以10 cm為間隔按照厚度順序放入鋁盒。為防止打鉆過程中將土壤壓實(shí)，每次打鉆厚度控制在30 cm以內(nèi)，直至采到母質(zhì)層。全剖面樣采完后，參照土壤發(fā)生分層標(biāo)準(zhǔn)劃分土層，并利用比色卡判定黑土厚度[25]?；厥覂?nèi)風(fēng)干后，再判定一次。
(2)剖面土樣采集。根據(jù)厚度分析結(jié)果，選擇其中18個地塊參照標(biāo)準(zhǔn)方法采集剖面土樣[25]，具體方法如下：在地塊中部地形變化不大且黑土層厚度接近地塊平均厚度的部位挖取剖面，以20 cm間隔分層采樣，直至母質(zhì)下20 cm為止，18個地塊共獲取分層土樣104個。每層采集3個標(biāo)準(zhǔn)環(huán)刀原狀土樣和1個混合土樣。
(3)土壤理化性質(zhì)測定。環(huán)刀樣本用于測定土壤容重?；旌贤寥缼Щ厥覂?nèi)風(fēng)干后參照相關(guān)標(biāo)準(zhǔn)[25]測定主要理化指標(biāo)，包括：有機(jī)質(zhì)含量(OM，g/100g)，機(jī)械組成(砂粒(2～0.02 mm)、粉粒(0.02～0.002 mm)和黏粒(<0.002 mm)，土壤凋萎濕度PWP(%)，田間持水量FC(%)、有效含水量AWC(%)以及pH值。OM含量用重鉻酸鉀氧化—外加熱法測定；機(jī)械組成采用吸管法測定；凋萎濕度和田間持水量用壓力膜法測定(壓力值分別設(shè)為1.5 MPa以及33 kPa)，有效持水量取二者之差；pH值用電位法測定，水土比采用2.5∶1。
1.3 土壤生產(chǎn)力指數(shù)計算
土壤生產(chǎn)力采用段興武等[15]修訂的東北黑土區(qū)土壤生產(chǎn)力評價模型MPI計算。該模型根據(jù)東北黑土的特性，剔除了原模型的容重因子，增加了土壤黏粒含量和有機(jī)質(zhì)含量因子，具體如下：
式中：MPI——土壤生產(chǎn)力指數(shù),范圍在0～1，數(shù)值越大，表明土壤的生產(chǎn)力越高；Ai——第i層土壤有效水含量的適宜性指數(shù)；Di——第i層土壤pH值的適宜性指數(shù)；Oi——第i層土壤有機(jī)質(zhì)含量的適宜性指數(shù)； CLi——第i層土壤黏粒含量的適宜性指數(shù)； WFi——第i層土壤的權(quán)重，反映不同土層理化性質(zhì)對土壤生產(chǎn)力的影響差異。各模型因子的適宜性指數(shù)計算參照段興武等[15]的研究結(jié)果。
為了評價模型的精度，收集了鶴山農(nóng)場調(diào)查地塊1990—2008年的大豆產(chǎn)量統(tǒng)計資料。其中2007年為重旱災(zāi)年份，產(chǎn)量明顯比往年偏低。為了剔除氣候?yàn)?zāi)害的影響，在結(jié)果分析中沒有采用2007年的產(chǎn)量數(shù)據(jù)。
2 結(jié)果與分析
2.1 土壤生產(chǎn)力指數(shù)與黑土厚度的關(guān)系
利用MPI模型計算了18個地塊的土壤生產(chǎn)力指數(shù)，各地塊的土壤剖面平均理化性質(zhì)及其生產(chǎn)力指數(shù)如表1所示。案例區(qū)土壤以壤土為主，有機(jī)質(zhì)含量平均2.04%，18個剖面的黑土層厚度變化在22～90 cm，平均41.56 cm，從整個黑土區(qū)上看，較開墾初期顯著降低[7-8，26-28]。土壤生產(chǎn)力指數(shù)變化在0.32～0.82，平均0.64，不同地塊間差異顯著。
表1 典型地塊土壤剖面平均理化性質(zhì)指標(biāo)和土壤生產(chǎn)力指數(shù)MPI
MPI與黑土厚度有顯著正相關(guān)，相關(guān)系數(shù)0.525(p<0.05)。MPI指數(shù)隨黑土厚度增加而增加，但在40—50 cm處有一個拐點(diǎn)：當(dāng)黑土厚度小于40—50 cm時，MPI隨黑土厚度增加速度較快；當(dāng)黑土厚度大于40—50 cm時，MPI隨黑土厚度增加速度有所減緩，二者間的關(guān)系可用對數(shù)函數(shù)表達(dá)(圖1)。這種變化可能與土壤剖面構(gòu)成有關(guān)，40—50 cm處土壤性質(zhì)有較為明顯的變化，下部土壤生產(chǎn)力較低，導(dǎo)致增加速度較小，因此，40—50 cm黑土層厚度可能是黑土保護(hù)的一個閾值厚度。
圖1 典型地塊黑土厚度與土壤生產(chǎn)力指數(shù)MPI間的關(guān)系
2.2 利用黑土厚度評價地塊土壤生產(chǎn)力
利用建立的MPI與黑土厚度之間的關(guān)系，計算了20個地塊122個厚度樣點(diǎn)的MPI及其平均值(表2)。為了驗(yàn)證黑土厚度評價土壤生產(chǎn)力的可行性，統(tǒng)計了各地塊MPI與1990—2008年平均大豆單產(chǎn)及當(dāng)?shù)赝恋氐燃壷g的關(guān)系。結(jié)果表明，地塊平均MPI與大豆單產(chǎn)為顯著正相關(guān)，與當(dāng)?shù)赝恋氐燃墳樨?fù)相關(guān)(圖2)，MPI隨大豆產(chǎn)量增加而增加，隨土地等級增高(1等地土壤質(zhì)量最高，5等地土壤質(zhì)量最低)而降低，表明利用黑土厚度評價黑土土壤生產(chǎn)力可行。
表2 典型地塊平均黑土厚度及MPI變化
圖2 土壤生產(chǎn)力指數(shù)MPI與大豆產(chǎn)量及土地等級間的關(guān)系
為了分析地塊土壤生產(chǎn)力的空間差異，在ArcGIS中利用克里金插值方法得到各個地塊的黑土厚度及MPI的空間分布(圖3)。122個樣點(diǎn)黑土厚度變化在15～125 cm，平均厚度為36.84 cm，其中85%以上樣點(diǎn)的黑土厚度小于50 cm。根據(jù)本文研究結(jié)果，隨著侵蝕程度的加劇，這些樣點(diǎn)的土壤生產(chǎn)力將迅速降低，未來土壤侵蝕導(dǎo)致土壤生產(chǎn)力下降不容忽視。為了維持生產(chǎn)力水平，必須采取行之有效的水土保持措施，如等高耕作、免耕等。20個地塊MPI變化在0.48～0.84，平均0.62，與東北典型黑土區(qū)的研究結(jié)果相比[15]，本案例區(qū)生產(chǎn)力水平變幅較小，均值略高，屬中等生產(chǎn)力水平。
黑土厚度與MPI空間分布一致：總體從坡頂?shù)狡碌?，土壤厚度增大，土壤生產(chǎn)力水平增加。高生產(chǎn)力土壤主要分布在底部地塊，如6-2-1和6-2-2號地塊；坡頂?shù)貕K土壤生產(chǎn)力最低，如6-5和9-5號地塊；位于坡中的大部分地塊生產(chǎn)力水平居中。從地塊內(nèi)部土壤生產(chǎn)力差異看，等高或接近等高耕作地塊的土壤生產(chǎn)力水平差異較小，如9-7，9-7-1，9-7-2，6-2-1和6-2-2號地塊等，而以順坡、接近順坡或與等高線有明顯交叉耕作的地塊，土壤生產(chǎn)力差異較大，如6-3-2，6-3-4，6-4-2和6-4-5號地塊等。
圖3 鶴北流域典型地塊黑土厚度及土壤生產(chǎn)力指數(shù)MPI空間分布
已有研究表明，盡管黑土區(qū)水土流失面積較大，但侵蝕的土壤物質(zhì)大部分沒有進(jìn)入河道，而是在低洼部位堆積下來，形成侵蝕和堆積共存的土壤物質(zhì)再分布景觀特征[29]。閻百興等[9]的研究結(jié)果也表明，漫崗區(qū)黑土從坡腳到坡頂侵蝕程度逐漸增加，由此造成的土壤質(zhì)量也逐漸降低。這些與本研究結(jié)果一致。流域內(nèi)農(nóng)地坡度大部分小于3°，除等高耕作外沒有其他水土保持措施。研究發(fā)現(xiàn)，等高耕作地塊土壤生產(chǎn)力較高，且空間差異不大，順坡或與等高線交叉耕作地塊土壤生產(chǎn)力較低且空間差異較大，說明耕作方式對土壤侵蝕影響的重要性。順坡耕作導(dǎo)致坡上黑土加速侵蝕，使結(jié)構(gòu)良好、土壤肥力水平高的表層黑土沉積于坡度較緩的坡腳地區(qū)[29-31]。作為三北防護(hù)林體系重要組成部分的東北地區(qū)，在20世紀(jì)70—80年代建設(shè)了大面積的防護(hù)林帶[32]。這些防護(hù)林逐漸成為黑土區(qū)農(nóng)地的天然地塊邊界，且決定了地塊的耕作方向(一般耕作方向都是平行于地塊長邊)。案例區(qū)當(dāng)前林帶的順坡布設(shè)，導(dǎo)致小流域內(nèi)地塊的耕作以順坡為主，這無疑加劇了土壤侵蝕速率和生產(chǎn)力退化進(jìn)程，因此，將來在考慮林帶更新時，應(yīng)注意林帶的等高布設(shè)，以便于等高耕作實(shí)施，從而減緩侵蝕速率，保護(hù)黑土生產(chǎn)力。
3 結(jié) 論
(1)黑土厚度與土壤生產(chǎn)力呈對數(shù)關(guān)系，隨著黑土厚度增加，土壤生產(chǎn)力水平逐漸增加，但在40—50 cm處有較明顯的變化，小于該厚度時，生產(chǎn)力隨厚度增加顯著；大于該厚度時，生產(chǎn)力隨厚度增加速度減緩，從而推測40—50 cm黑土厚度可能是黑土生產(chǎn)力保護(hù)的厚度閾值，應(yīng)進(jìn)一步深入研究。
(2)利用黑土厚度評價土壤生產(chǎn)力簡易可行，結(jié)果得到了大豆單產(chǎn)及當(dāng)?shù)赝恋刭|(zhì)量等級的驗(yàn)證。案例區(qū)土壤生產(chǎn)力總體從坡頂?shù)狡履_有明顯增加趨勢，等高或接近等高耕作地塊的土壤生產(chǎn)力較高，且地塊內(nèi)生產(chǎn)力空間差異小，順坡、接近順坡或與等高線交叉耕作的地塊，土壤生產(chǎn)力空間差異較大，耕作方式是導(dǎo)致該區(qū)黑土生產(chǎn)力水平退化的重要原因。建議將來林帶更新時，應(yīng)將林帶位置及方向等的選擇與機(jī)械化水土保持耕作措施的實(shí)施相結(jié)合，以減緩黑土侵蝕速率，保護(hù)黑土生產(chǎn)力。
[參考文獻(xiàn)]
[1]劉寶元,閻百興,沈波,等.東北黑土區(qū)農(nóng)地水土流失現(xiàn)狀與綜合治理對策[J].中國水土保持科學(xué),2008,6(1):1-8.
[2]劉興土,閻百興.東北黑土區(qū)水土流失與糧食安全[J].中國水土保持,2009(1):17-19.
[3]全國土壤普查辦公室.中國土壤[M].北京:中國農(nóng)業(yè)出版社,1998:318-320.
[4]劉振舉,蔡強(qiáng)國,李國強(qiáng).東北黑土區(qū)水土流失脅迫下的生態(tài)安全評價指標(biāo)體系的建立[J].水土保持通報,2009，29(6):64-69.
[5]李發(fā)鵬,李景玉,徐宗學(xué).東北黑土區(qū)土壤退化及水土流失研究現(xiàn)狀[J].水土保持研究,2006,31(3):50-54.
[6]劉丙友.典型黑土區(qū)土壤退化及可持續(xù)利用問題探討[J].中國水土保持,2003(12):28-29.
[7]丁瑞興,劉樹桐.黑土開墾后肥力演變的研究[J].土壤學(xué)報,1980,17(1):20-30.
[8]辛剛,顏麗,汪景寬,等.不同開墾年限黑土有機(jī)質(zhì)變化的研究[J].土壤通報,2002,33(5):332-335.
[9]閻百興,湯潔.黑土侵蝕速率及其對土壤質(zhì)量的影響[J].地理研究,2005,24(4):499-506.
[10]唐近春.中國土壤肥料工作的成就與任務(wù)[J].土壤學(xué)報,1994,31(4):341-347.
[11]李忠芳,徐明崗,張會民,等.長期施肥下中國主要糧食作物產(chǎn)量的變化[J].中國農(nóng)業(yè)科學(xué),2009,42(7):2407-2414.
[12]林葆,林繼雄,李家康.長期施肥的作物產(chǎn)量和土壤肥力變化[J].植物營養(yǎng)與肥料學(xué)報,1994(1):6-18.
[13]劉恩科,趙秉強(qiáng),胡昌浩,等.長期施氮、磷、鉀化肥對玉米產(chǎn)量及其土壤肥力的影響[J].植物營養(yǎng)與肥料學(xué)報,2007,13(5):789-941.
[14]徐明崗,梁國慶,張夫道.中國土壤肥力演變[M].北京：中國農(nóng)業(yè)科學(xué)技術(shù)出版社,2006.
[15]段興武,謝云,馮艷杰,等.東北黑土區(qū)土壤生產(chǎn)力評價方法研究[J].中國農(nóng)業(yè)科學(xué),2009,42(5):1656-1664.
[16]段興武,謝云,張玉平,等.PI 模型在東北松嫩黑土區(qū)土壤生產(chǎn)力評價中的應(yīng)用[J].中國農(nóng)學(xué)通報,2010,26(8):179-188.
[17]Pierce F J, Larson W E, Dowdy R H, et al. Productivity of soil: Assessing of long-term changes due to erosion[J] . J. of Soil and Water Cons., 1983,38(1):39-44.
[18]Yang J, Hammer R D, Thompson A L, et al. Predicting soybean yields in dry and wet year using a soil productivity index[J]. Plant and Soil, 2003,250(2):175-182.
[19]Lobo D, Lozano Z, Delgado F.Water erosion risk assessment and impact on productivity of a Venezuelan soil[J].Catena, 2005,64(3):297-306.
[20]Garcia-paredes J D, Olson K R, Lang J M. Predicting corn and soybean productivity for Illinois soils[J]. Agricultural Systems, 2000,64(4):151-170.
[21]謝云,段興武,劉寶元,等.東北黑土區(qū)主要黑土土種的容許土壤流失量[J].地理學(xué)報,2011,66(7):940-952.
[22]Duan Xingwu, Xie Yun, Ou Tinghai,et al.Effects of soil erosion on long-term soil productivity in the black soil region of northeastern China[J]. Catena, 2011,87(2):268-275.
[23]張之一.黑土開墾后黑土層厚度的變化[J].黑龍江八一農(nóng)墾大學(xué)學(xué)報,2010,22(5):1-3.
[24]中國科學(xué)院南京土壤研究所系統(tǒng)分類課題組.中國土壤系統(tǒng)分類:修訂方案[M].北京:中國農(nóng)業(yè)科技出版社,1995:1-218.
[25]劉光菘.土壤理化分析與剖面描述[M].北京:中國標(biāo)準(zhǔn)出版社,1996.
[26]汪景寬,王鐵宇,張旭東,等.黑土土壤質(zhì)量演變初探(I): 不同開墾年限黑土主要質(zhì)量指標(biāo)演變規(guī)律[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報,2002,33(1):43-47.
[27]黃健,張惠琳,傅文玉,等.東北黑土區(qū)土壤肥力變化特征的分析[J].土壤通報,2005,36(5):659-663.
[28]汪景寬,李雙異,張旭東,等.20年來東北典型黑土地區(qū)土壤肥力質(zhì)量變化[J].中國生態(tài)農(nóng)業(yè)學(xué)報,2007,15(1):19-24.
[29]張曉平,梁愛珍,申艷，等.東北黑土水土流失特點(diǎn)[J].地理科學(xué),2006,26(6):687-692.
[30]方華軍,楊學(xué)明,張曉平,等.137Cs示蹤技術(shù)研究坡耕地黑土侵蝕和沉積特征[J].生態(tài)學(xué)報,2005,25(6):1376-1382.
[31]梁愛珍,張曉平,楊學(xué)明,等.耕作方式對耕層黑土有機(jī)碳庫儲量的短期影響[J].中國農(nóng)業(yè)科學(xué),2006,39(6):1287-1293.
[32]張光輝,于艷,劉寶元.東北黑土區(qū)防護(hù)林帶對土壤水分影響的研究[J].水土保持學(xué)報,2006,20(6):54-57.
猜你喜歡
自然資源部：加強(qiáng)黑土耕地保護(hù)中國化肥信息(2022年8期)2022-12-05 20:01:27 基于Sentinel-2遙感影像的黑土區(qū)土壤有效磷反演中國農(nóng)業(yè)信息(2022年2期)2022-07-15 09:07:16 把“風(fēng)景”轉(zhuǎn)化成“生產(chǎn)力紅土地(2016年3期)2017-01-15 13:45:17 耕作深度對紫色土坡地旋耕機(jī)耕作侵蝕的影響水土保持研究(2016年4期)2016-10-21 02:28:59 玉米保護(hù)性耕作的技術(shù)要領(lǐng)新農(nóng)業(yè)(2016年15期)2016-08-16 03:40:48 草地耕作技術(shù)在澳大利亞的應(yīng)用土壤與作物(2015年3期)2015-12-08 00:46:58 人口紅利與提高生產(chǎn)力海峽姐妹(2015年9期)2015-02-27 15:13:07 典型黑土區(qū)不同尺度觀測場地融雪徑流水土保持通報(2014年5期)2014-06-09 08:27:26 多菌靈在酸化黑土中遷移性能的抑制水土保持通報(2014年5期)2014-06-09 08:27:10 西洞庭湖區(qū)免耕耕作模式及其配套技術(shù)作物研究(2014年6期)2014-03-01 03:39:12
李 驁， 段興武
(1．華中農(nóng)業(yè)大學(xué) 430070， 湖北 武漢； 2．云南大學(xué) 亞洲國際河流研究中心， 云南 昆明 650091)
東北黑土區(qū)是中國重要的糧食生產(chǎn)基地。據(jù)農(nóng)業(yè)部農(nóng)業(yè)信息網(wǎng)公布的統(tǒng)計資料，2011年東北三省的糧食總產(chǎn)量為1.07×1011kg，約占全國總產(chǎn)的18.86%，其中玉米產(chǎn)量占全國的33.06%，大豆占46.78%，人均糧食占有量是全國平均水平的2.4倍，在保障國家糧食安全問題上起著舉足輕重的作用，被譽(yù)為中國糧食生產(chǎn)的“穩(wěn)壓器”[1-2]。
黑土養(yǎng)分含量豐富、結(jié)構(gòu)良好、土壤肥力極高[3]。自從20世紀(jì)初大面積開墾以來，黑土土壤理化性質(zhì)發(fā)生了顯著變化[4-8]，土壤生產(chǎn)力顯著降低[7-9]，引起黑土土壤生產(chǎn)力變化的原因比較復(fù)雜，可主要分為施肥和侵蝕兩種。為提高作物產(chǎn)量，自20世紀(jì)70年代末80年代初開始，中國開始大量施用化工肥料[10]，化工肥料的施用能顯著提高作物產(chǎn)量[11]，同時也能改變土壤的理化性質(zhì)[12-13]，從而影響到作物品質(zhì)的提高和土壤肥力的可持續(xù)利用[14]。另一方面，在高強(qiáng)度的人為耕作下，黑土區(qū)土壤侵蝕嚴(yán)重[1-4]。據(jù)水利部松遼水利委員會2002年的統(tǒng)計報告，該區(qū)水土流失面積占全區(qū)總面積的27%，發(fā)育大型切溝2.5×104條余，吞噬耕地面積約3.9 hm2[5]。每年僅黑龍江省流失的土壤就變化在2.0×108～3.0×108m3，流失的養(yǎng)分折合標(biāo)準(zhǔn)化肥5.0×109～6.0×109kg，相當(dāng)于全省化肥施用總量[6]，嚴(yán)重的侵蝕造成土壤養(yǎng)分損失、土壤質(zhì)地粗化、持水能力降低，直接導(dǎo)致土壤生產(chǎn)力退化[13]。在這樣的前提下，研究簡便易行的土壤生產(chǎn)力評價指標(biāo)可為黑土資源的合理利用和保護(hù)提供理論依據(jù)和方法參考。
土壤生產(chǎn)力評價方法可分為定性和定量兩種[15]。定性評價是通過對田間土壤的實(shí)地調(diào)查和診斷，對土壤宏觀性狀進(jìn)行定性描述，多帶有經(jīng)驗(yàn)性，很難量化[16]。定量評價主要有3種：一是單指標(biāo)評價；二是多指標(biāo)表土綜合評價，即選取表層土壤多個理化性質(zhì)指標(biāo)評價表土生產(chǎn)力；三是多指標(biāo)土體綜合評價，即將整個土體看作是作物生長的環(huán)境，選擇各個土層的土壤性質(zhì)指標(biāo)，將不同深度土層賦予不同的權(quán)重，綜合、立體地評價土壤生產(chǎn)力，比較有代表性的是土壤生產(chǎn)力指數(shù)模型(productivity index, PI)[17]，該模型在國內(nèi)外獲得了廣泛應(yīng)用[18-20]。段興武等[15]在該模型基礎(chǔ)上，根據(jù)東北黑土土壤理化性質(zhì)和作物產(chǎn)量調(diào)查結(jié)果，構(gòu)建了適用于東北黑土區(qū)的修正的生產(chǎn)力指數(shù)模型(modified productivity Index, MPI)并得到多次驗(yàn)證[21-22]。但該模型仍然需要采集剖面土壤樣本、測定多項(xiàng)土壤理化指標(biāo)，不利于野外快速和大范圍土壤生產(chǎn)力調(diào)查與評價。表層黑土養(yǎng)分含量豐富、土壤結(jié)構(gòu)良好，表層黑土厚度(簡稱黑土層厚度)便于野外調(diào)查，對黑土土壤生產(chǎn)力貢獻(xiàn)較大[3]，被認(rèn)為是衡量黑土退化程度的重要指標(biāo)[7-8,23]。但時至今日，鮮有探討黑土層厚度與土壤生產(chǎn)力的研究案例。
本文以位于東北典型黑土區(qū)北端的鶴北小流域?yàn)榘咐齾^(qū)，通過調(diào)查和分析流域內(nèi)不同等級土地的土壤厚度及土壤理化性質(zhì)，構(gòu)建了黑土層厚度與土壤生產(chǎn)力間的關(guān)系，并據(jù)此評價典型地塊的土壤生產(chǎn)力空間特征及其影響因素。
1 資料與方法
1.1 研究區(qū)概況
研究區(qū)位于黑龍江省農(nóng)墾總局九三分局鶴山農(nóng)場境內(nèi)的鶴北小流域，地理坐標(biāo)為48°59′0″—49°3′0″N，125°16′0″—125°121′0″E，屬于寒溫帶大陸性季風(fēng)氣候，冬季漫長，夏季短促，秋季涼爽初霜早，春季風(fēng)大，年平均氣溫-1.5～0.4 ℃，歷年最高氣溫37.9 ℃，最低氣溫-46.5 ℃，年降水量在400～500 mm，集中于6—9月。該區(qū)位于小興安嶺南麓，是小興安嶺向嫩江平原過渡地帶，屬漫崗丘陵區(qū)，坡長一般800～1 500 m，坡度一般1°～4°。自然植被屬于溫帶半濕潤草甸草原，土壤屬中國發(fā)生學(xué)分類中的黑土，屬中國土壤系統(tǒng)分類中的均腐土[24]。
鶴北小流域的面積為28 km2，大部分已開墾為農(nóng)地，且開墾歷史都在60 a以上。開墾之初，通過種植人工林帶將各個地塊隔開，且無論是起壟或行播均沿林帶方向?qū)嵭袡C(jī)械化作業(yè)。流域內(nèi)共包括41個地塊，種植的農(nóng)作物主要有大豆、小麥和玉米，尤以大豆為主。農(nóng)業(yè)機(jī)械化程度高，整地、播種、噴藥、施肥、除草和收割等全部實(shí)行機(jī)械化作業(yè)。
1.2 樣本采集
本研究根據(jù)地形狀況，選擇流域內(nèi)20個典型地塊調(diào)查土壤厚度和剖面土壤理化性質(zhì)。
(1)黑土厚度調(diào)查。從調(diào)查地塊的坡頂開始每10 m高差布設(shè)一個厚度樣點(diǎn)，共計布設(shè)厚度樣點(diǎn)122個。利用半圓鉆(長1 m，可加延長桿，直徑4 cm)采集土樣，將樣本以10 cm為間隔按照厚度順序放入鋁盒。為防止打鉆過程中將土壤壓實(shí)，每次打鉆厚度控制在30 cm以內(nèi)，直至采到母質(zhì)層。全剖面樣采完后，參照土壤發(fā)生分層標(biāo)準(zhǔn)劃分土層，并利用比色卡判定黑土厚度[25]?；厥覂?nèi)風(fēng)干后，再判定一次。
(2)剖面土樣采集。根據(jù)厚度分析結(jié)果，選擇其中18個地塊參照標(biāo)準(zhǔn)方法采集剖面土樣[25]，具體方法如下：在地塊中部地形變化不大且黑土層厚度接近地塊平均厚度的部位挖取剖面，以20 cm間隔分層采樣，直至母質(zhì)下20 cm為止，18個地塊共獲取分層土樣104個。每層采集3個標(biāo)準(zhǔn)環(huán)刀原狀土樣和1個混合土樣。
(3)土壤理化性質(zhì)測定。環(huán)刀樣本用于測定土壤容重?；旌贤寥缼Щ厥覂?nèi)風(fēng)干后參照相關(guān)標(biāo)準(zhǔn)[25]測定主要理化指標(biāo)，包括：有機(jī)質(zhì)含量(OM，g/100g)，機(jī)械組成(砂粒(2～0.02 mm)、粉粒(0.02～0.002 mm)和黏粒(<0.002 mm)，土壤凋萎濕度PWP(%)，田間持水量FC(%)、有效含水量AWC(%)以及pH值。OM含量用重鉻酸鉀氧化—外加熱法測定；機(jī)械組成采用吸管法測定；凋萎濕度和田間持水量用壓力膜法測定(壓力值分別設(shè)為1.5 MPa以及33 kPa)，有效持水量取二者之差；pH值用電位法測定，水土比采用2.5∶1。
1.3 土壤生產(chǎn)力指數(shù)計算
土壤生產(chǎn)力采用段興武等[15]修訂的東北黑土區(qū)土壤生產(chǎn)力評價模型MPI計算。該模型根據(jù)東北黑土的特性，剔除了原模型的容重因子，增加了土壤黏粒含量和有機(jī)質(zhì)含量因子，具體如下：
式中：MPI——土壤生產(chǎn)力指數(shù),范圍在0～1，數(shù)值越大，表明土壤的生產(chǎn)力越高；Ai——第i層土壤有效水含量的適宜性指數(shù)；Di——第i層土壤pH值的適宜性指數(shù)；Oi——第i層土壤有機(jī)質(zhì)含量的適宜性指數(shù)； CLi——第i層土壤黏粒含量的適宜性指數(shù)； WFi——第i層土壤的權(quán)重，反映不同土層理化性質(zhì)對土壤生產(chǎn)力的影響差異。各模型因子的適宜性指數(shù)計算參照段興武等[15]的研究結(jié)果。
為了評價模型的精度，收集了鶴山農(nóng)場調(diào)查地塊1990—2008年的大豆產(chǎn)量統(tǒng)計資料。其中2007年為重旱災(zāi)年份，產(chǎn)量明顯比往年偏低。為了剔除氣候?yàn)?zāi)害的影響，在結(jié)果分析中沒有采用2007年的產(chǎn)量數(shù)據(jù)。
2 結(jié)果與分析
2.1 土壤生產(chǎn)力指數(shù)與黑土厚度的關(guān)系
利用MPI模型計算了18個地塊的土壤生產(chǎn)力指數(shù)，各地塊的土壤剖面平均理化性質(zhì)及其生產(chǎn)力指數(shù)如表1所示。案例區(qū)土壤以壤土為主，有機(jī)質(zhì)含量平均2.04%，18個剖面的黑土層厚度變化在22～90 cm，平均41.56 cm，從整個黑土區(qū)上看，較開墾初期顯著降低[7-8，26-28]。土壤生產(chǎn)力指數(shù)變化在0.32～0.82，平均0.64，不同地塊間差異顯著。
表1 典型地塊土壤剖面平均理化性質(zhì)指標(biāo)和土壤生產(chǎn)力指數(shù)MPI
MPI與黑土厚度有顯著正相關(guān)，相關(guān)系數(shù)0.525(p<0.05)。MPI指數(shù)隨黑土厚度增加而增加，但在40—50 cm處有一個拐點(diǎn)：當(dāng)黑土厚度小于40—50 cm時，MPI隨黑土厚度增加速度較快；當(dāng)黑土厚度大于40—50 cm時，MPI隨黑土厚度增加速度有所減緩，二者間的關(guān)系可用對數(shù)函數(shù)表達(dá)(圖1)。這種變化可能與土壤剖面構(gòu)成有關(guān)，40—50 cm處土壤性質(zhì)有較為明顯的變化，下部土壤生產(chǎn)力較低，導(dǎo)致增加速度較小，因此，40—50 cm黑土層厚度可能是黑土保護(hù)的一個閾值厚度。
圖1 典型地塊黑土厚度與土壤生產(chǎn)力指數(shù)MPI間的關(guān)系
2.2 利用黑土厚度評價地塊土壤生產(chǎn)力
利用建立的MPI與黑土厚度之間的關(guān)系，計算了20個地塊122個厚度樣點(diǎn)的MPI及其平均值(表2)。為了驗(yàn)證黑土厚度評價土壤生產(chǎn)力的可行性，統(tǒng)計了各地塊MPI與1990—2008年平均大豆單產(chǎn)及當(dāng)?shù)赝恋氐燃壷g的關(guān)系。結(jié)果表明，地塊平均MPI與大豆單產(chǎn)為顯著正相關(guān)，與當(dāng)?shù)赝恋氐燃墳樨?fù)相關(guān)(圖2)，MPI隨大豆產(chǎn)量增加而增加，隨土地等級增高(1等地土壤質(zhì)量最高，5等地土壤質(zhì)量最低)而降低，表明利用黑土厚度評價黑土土壤生產(chǎn)力可行。
表2 典型地塊平均黑土厚度及MPI變化
圖2 土壤生產(chǎn)力指數(shù)MPI與大豆產(chǎn)量及土地等級間的關(guān)系
為了分析地塊土壤生產(chǎn)力的空間差異，在ArcGIS中利用克里金插值方法得到各個地塊的黑土厚度及MPI的空間分布(圖3)。122個樣點(diǎn)黑土厚度變化在15～125 cm，平均厚度為36.84 cm，其中85%以上樣點(diǎn)的黑土厚度小于50 cm。根據(jù)本文研究結(jié)果，隨著侵蝕程度的加劇，這些樣點(diǎn)的土壤生產(chǎn)力將迅速降低，未來土壤侵蝕導(dǎo)致土壤生產(chǎn)力下降不容忽視。為了維持生產(chǎn)力水平，必須采取行之有效的水土保持措施，如等高耕作、免耕等。20個地塊MPI變化在0.48～0.84，平均0.62，與東北典型黑土區(qū)的研究結(jié)果相比[15]，本案例區(qū)生產(chǎn)力水平變幅較小，均值略高，屬中等生產(chǎn)力水平。
黑土厚度與MPI空間分布一致：總體從坡頂?shù)狡碌?，土壤厚度增大，土壤生產(chǎn)力水平增加。高生產(chǎn)力土壤主要分布在底部地塊，如6-2-1和6-2-2號地塊；坡頂?shù)貕K土壤生產(chǎn)力最低，如6-5和9-5號地塊；位于坡中的大部分地塊生產(chǎn)力水平居中。從地塊內(nèi)部土壤生產(chǎn)力差異看，等高或接近等高耕作地塊的土壤生產(chǎn)力水平差異較小，如9-7，9-7-1，9-7-2，6-2-1和6-2-2號地塊等，而以順坡、接近順坡或與等高線有明顯交叉耕作的地塊，土壤生產(chǎn)力差異較大，如6-3-2，6-3-4，6-4-2和6-4-5號地塊等。
圖3 鶴北流域典型地塊黑土厚度及土壤生產(chǎn)力指數(shù)MPI空間分布
已有研究表明，盡管黑土區(qū)水土流失面積較大，但侵蝕的土壤物質(zhì)大部分沒有進(jìn)入河道，而是在低洼部位堆積下來，形成侵蝕和堆積共存的土壤物質(zhì)再分布景觀特征[29]。閻百興等[9]的研究結(jié)果也表明，漫崗區(qū)黑土從坡腳到坡頂侵蝕程度逐漸增加，由此造成的土壤質(zhì)量也逐漸降低。這些與本研究結(jié)果一致。流域內(nèi)農(nóng)地坡度大部分小于3°，除等高耕作外沒有其他水土保持措施。研究發(fā)現(xiàn)，等高耕作地塊土壤生產(chǎn)力較高，且空間差異不大，順坡或與等高線交叉耕作地塊土壤生產(chǎn)力較低且空間差異較大，說明耕作方式對土壤侵蝕影響的重要性。順坡耕作導(dǎo)致坡上黑土加速侵蝕，使結(jié)構(gòu)良好、土壤肥力水平高的表層黑土沉積于坡度較緩的坡腳地區(qū)[29-31]。作為三北防護(hù)林體系重要組成部分的東北地區(qū)，在20世紀(jì)70—80年代建設(shè)了大面積的防護(hù)林帶[32]。這些防護(hù)林逐漸成為黑土區(qū)農(nóng)地的天然地塊邊界，且決定了地塊的耕作方向(一般耕作方向都是平行于地塊長邊)。案例區(qū)當(dāng)前林帶的順坡布設(shè)，導(dǎo)致小流域內(nèi)地塊的耕作以順坡為主，這無疑加劇了土壤侵蝕速率和生產(chǎn)力退化進(jìn)程，因此，將來在考慮林帶更新時，應(yīng)注意林帶的等高布設(shè)，以便于等高耕作實(shí)施，從而減緩侵蝕速率，保護(hù)黑土生產(chǎn)力。
3 結(jié) 論
(1)黑土厚度與土壤生產(chǎn)力呈對數(shù)關(guān)系，隨著黑土厚度增加，土壤生產(chǎn)力水平逐漸增加，但在40—50 cm處有較明顯的變化，小于該厚度時，生產(chǎn)力隨厚度增加顯著；大于該厚度時，生產(chǎn)力隨厚度增加速度減緩，從而推測40—50 cm黑土厚度可能是黑土生產(chǎn)力保護(hù)的厚度閾值，應(yīng)進(jìn)一步深入研究。
(2)利用黑土厚度評價土壤生產(chǎn)力簡易可行，結(jié)果得到了大豆單產(chǎn)及當(dāng)?shù)赝恋刭|(zhì)量等級的驗(yàn)證。案例區(qū)土壤生產(chǎn)力總體從坡頂?shù)狡履_有明顯增加趨勢，等高或接近等高耕作地塊的土壤生產(chǎn)力較高，且地塊內(nèi)生產(chǎn)力空間差異小，順坡、接近順坡或與等高線交叉耕作的地塊，土壤生產(chǎn)力空間差異較大，耕作方式是導(dǎo)致該區(qū)黑土生產(chǎn)力水平退化的重要原因。建議將來林帶更新時，應(yīng)將林帶位置及方向等的選擇與機(jī)械化水土保持耕作措施的實(shí)施相結(jié)合，以減緩黑土侵蝕速率，保護(hù)黑土生產(chǎn)力。
[參考文獻(xiàn)]
[1]劉寶元,閻百興,沈波,等.東北黑土區(qū)農(nóng)地水土流失現(xiàn)狀與綜合治理對策[J].中國水土保持科學(xué),2008,6(1):1-8.
[2]劉興土,閻百興.東北黑土區(qū)水土流失與糧食安全[J].中國水土保持,2009(1):17-19.
[3]全國土壤普查辦公室.中國土壤[M].北京:中國農(nóng)業(yè)出版社,1998:318-320.
[4]劉振舉,蔡強(qiáng)國,李國強(qiáng).東北黑土區(qū)水土流失脅迫下的生態(tài)安全評價指標(biāo)體系的建立[J].水土保持通報,2009，29(6):64-69.
[5]李發(fā)鵬,李景玉,徐宗學(xué).東北黑土區(qū)土壤退化及水土流失研究現(xiàn)狀[J].水土保持研究,2006,31(3):50-54.
[6]劉丙友.典型黑土區(qū)土壤退化及可持續(xù)利用問題探討[J].中國水土保持,2003(12):28-29.
[7]丁瑞興,劉樹桐.黑土開墾后肥力演變的研究[J].土壤學(xué)報,1980,17(1):20-30.
[8]辛剛,顏麗,汪景寬,等.不同開墾年限黑土有機(jī)質(zhì)變化的研究[J].土壤通報,2002,33(5):332-335.
[9]閻百興,湯潔.黑土侵蝕速率及其對土壤質(zhì)量的影響[J].地理研究,2005,24(4):499-506.
[10]唐近春.中國土壤肥料工作的成就與任務(wù)[J].土壤學(xué)報,1994,31(4):341-347.
[11]李忠芳,徐明崗,張會民,等.長期施肥下中國主要糧食作物產(chǎn)量的變化[J].中國農(nóng)業(yè)科學(xué),2009,42(7):2407-2414.
[12]林葆,林繼雄,李家康.長期施肥的作物產(chǎn)量和土壤肥力變化[J].植物營養(yǎng)與肥料學(xué)報,1994(1):6-18.
[13]劉恩科,趙秉強(qiáng),胡昌浩,等.長期施氮、磷、鉀化肥對玉米產(chǎn)量及其土壤肥力的影響[J].植物營養(yǎng)與肥料學(xué)報,2007,13(5):789-941.
[14]徐明崗,梁國慶,張夫道.中國土壤肥力演變[M].北京：中國農(nóng)業(yè)科學(xué)技術(shù)出版社,2006.
[15]段興武,謝云,馮艷杰,等.東北黑土區(qū)土壤生產(chǎn)力評價方法研究[J].中國農(nóng)業(yè)科學(xué),2009,42(5):1656-1664.
[16]段興武,謝云,張玉平,等.PI 模型在東北松嫩黑土區(qū)土壤生產(chǎn)力評價中的應(yīng)用[J].中國農(nóng)學(xué)通報,2010,26(8):179-188.
[17]Pierce F J, Larson W E, Dowdy R H, et al. Productivity of soil: Assessing of long-term changes due to erosion[J] . J. of Soil and Water Cons., 1983,38(1):39-44.
[18]Yang J, Hammer R D, Thompson A L, et al. Predicting soybean yields in dry and wet year using a soil productivity index[J]. Plant and Soil, 2003,250(2):175-182.
[19]Lobo D, Lozano Z, Delgado F.Water erosion risk assessment and impact on productivity of a Venezuelan soil[J].Catena, 2005,64(3):297-306.
[20]Garcia-paredes J D, Olson K R, Lang J M. Predicting corn and soybean productivity for Illinois soils[J]. Agricultural Systems, 2000,64(4):151-170.
[21]謝云,段興武,劉寶元,等.東北黑土區(qū)主要黑土土種的容許土壤流失量[J].地理學(xué)報,2011,66(7):940-952.
[22]Duan Xingwu, Xie Yun, Ou Tinghai,et al.Effects of soil erosion on long-term soil productivity in the black soil region of northeastern China[J]. Catena, 2011,87(2):268-275.
[23]張之一.黑土開墾后黑土層厚度的變化[J].黑龍江八一農(nóng)墾大學(xué)學(xué)報,2010,22(5):1-3.
[24]中國科學(xué)院南京土壤研究所系統(tǒng)分類課題組.中國土壤系統(tǒng)分類:修訂方案[M].北京:中國農(nóng)業(yè)科技出版社,1995:1-218.
[25]劉光菘.土壤理化分析與剖面描述[M].北京:中國標(biāo)準(zhǔn)出版社,1996.
[26]汪景寬,王鐵宇,張旭東,等.黑土土壤質(zhì)量演變初探(I): 不同開墾年限黑土主要質(zhì)量指標(biāo)演變規(guī)律[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報,2002,33(1):43-47.
[27]黃健,張惠琳,傅文玉,等.東北黑土區(qū)土壤肥力變化特征的分析[J].土壤通報,2005,36(5):659-663.
[28]汪景寬,李雙異,張旭東,等.20年來東北典型黑土地區(qū)土壤肥力質(zhì)量變化[J].中國生態(tài)農(nóng)業(yè)學(xué)報,2007,15(1):19-24.
[29]張曉平,梁愛珍,申艷，等.東北黑土水土流失特點(diǎn)[J].地理科學(xué),2006,26(6):687-692.
[30]方華軍,楊學(xué)明,張曉平,等.137Cs示蹤技術(shù)研究坡耕地黑土侵蝕和沉積特征[J].生態(tài)學(xué)報,2005,25(6):1376-1382.
[31]梁愛珍,張曉平,楊學(xué)明,等.耕作方式對耕層黑土有機(jī)碳庫儲量的短期影響[J].中國農(nóng)業(yè)科學(xué),2006,39(6):1287-1293.
[32]張光輝,于艷,劉寶元.東北黑土區(qū)防護(hù)林帶對土壤水分影響的研究[J].水土保持學(xué)報,2006,20(6):54-57.
猜你喜歡
中國化肥信息(2022年8期)2022-12-05 20:01:27
中國農(nóng)業(yè)信息(2022年2期)2022-07-15 09:07:16
紅土地(2016年3期)2017-01-15 13:45:17
水土保持研究(2016年4期)2016-10-21 02:28:59
新農(nóng)業(yè)(2016年15期)2016-08-16 03:40:48
土壤與作物(2015年3期)2015-12-08 00:46:58
海峽姐妹(2015年9期)2015-02-27 15:13:07
水土保持通報(2014年5期)2014-06-09 08:27:26
水土保持通報(2014年5期)2014-06-09 08:27:10
作物研究(2014年6期)2014-03-01 03:39:12
