黑龍江省松嫩平原南部表層土壤肥力質(zhì)量評價研究

徐曉嘉,雷國平,張 慧,王元輝

（東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,哈爾濱 150030）

土壤肥力質(zhì)量是土壤提供植物養(yǎng)分和生產(chǎn)生物物質(zhì)的能力,是保障糧食生產(chǎn)的根本[1]。土壤肥力質(zhì)量主要取決于土壤的自然組成部分,同時也受到來自利用管理等人類行為的影響。作為一個復(fù)雜的功能實體,土壤肥力質(zhì)量不能夠直接測定,而是主要取決于其地理位置、土壤類型以及土壤本身各種元素間的相互作用,因此,對影響土壤肥力質(zhì)量的各種元素進行空間變異特征分析后,通過定量的土壤肥力質(zhì)量指標(biāo)來計算得到土壤肥力質(zhì)量等級,能夠更加直觀地揭示土壤養(yǎng)分含量及肥力質(zhì)量分布狀況。近年來采用數(shù)學(xué)的方法對土壤肥力質(zhì)量進行數(shù)量化評價,取得了一定的效果,Fuzzy綜合評判法、指數(shù)和法、聚類分析法、灰色關(guān)聯(lián)分析法等方法的研究與應(yīng)用,使土壤肥力質(zhì)量的研究取得了顯著進展[2]。同時,土壤空間變異性研究越來越受到人們的重視,特別是在土壤養(yǎng)分的管理上,以精準(zhǔn)農(nóng)業(yè)為中心的養(yǎng)分管理正成為土壤養(yǎng)分管理研究的重點[3]。

松嫩平原位于東北平原的中部,是我國最好的一熟制作物種植區(qū)和國家最為重要的商品農(nóng)業(yè)基地之一,該地區(qū)的糧食生產(chǎn)對于保障國家糧食安全具有舉足輕重的地位[4]。本文選取黑龍江省松嫩平原南部地區(qū)表層土壤作為研究對象應(yīng)用地理信息系統(tǒng)軟件對研究區(qū)內(nèi)土壤肥力元素全量進行空間變異特征分析,利用隸屬度函數(shù)模型和加法模型對研究區(qū)表層土壤進行肥力質(zhì)量綜合評價,以期能夠為研究區(qū)土地合理施肥、實現(xiàn)精準(zhǔn)農(nóng)業(yè)以及進行區(qū)域內(nèi)縣級土地質(zhì)量評價提供借鑒。

1 研究區(qū)概況

研究區(qū)位于黑龍江省松嫩平原南部,地理坐標(biāo)為東經(jīng) 122°25′-127°45′,北緯 45°35′-48°00′,總面積為8.15萬km2。該區(qū)地處中緯度亞洲大陸東岸,屬于中溫帶。冬季在極地大陸氣團控制下,氣候嚴(yán)寒干燥;夏季受副熱帶海洋氣團的影響,降水集中,氣候溫暖濕潤,具有明顯的大陸性季風(fēng)氣候特征,因而一年中的太陽高度角變化及各季節(jié)間的太陽輻射量的變化均較大,致使氣候年變化大,四季分明,冬季漫長而嚴(yán)寒,夏季短促而炎熱。研究區(qū)內(nèi)共有土壤類型9種,分別是暗棕壤、白漿土、黑土、黑鈣土、草甸土、沼澤土、堿土、風(fēng)砂土、水稻土,其中以黑土、黑鈣土和草甸土為主,這三類土壤約占松嫩平原區(qū)面積的90%以上。

2 土壤樣品采集及數(shù)據(jù)分析

2.1 土壤樣品采集

土壤樣品采集的工作比例尺為1∶250 000,采樣點以1個點/km2的密度按照網(wǎng)格化布置,按0-20 cm的采樣深度采集表層土壤樣本20 929件。采樣點分布如圖1所示。

圖1 樣點分布圖

土壤樣品的編號、采樣深度、采樣方法、定點與標(biāo)繪、采樣記錄、重復(fù)樣采集、樣品重量、樣品加工、樣品分析及質(zhì)量監(jiān)控等各項內(nèi)容均按照中國地質(zhì)調(diào)查局《多目標(biāo)區(qū)域地球化學(xué)調(diào)查規(guī)范(1∶250 000)》DD2005-01[5]中土壤表層樣品規(guī)范的要求進行。應(yīng)用X射線熒光光譜法等方法對18種土壤肥力元素全量進行分析。

2.2 土壤元素空間變異特征分析

土壤作為歷史自然體,受氣候、生物、母質(zhì)、地形、成土?xí)r間等成土因素的影響,具有復(fù)雜性和時空變異性,其隨空間位置發(fā)生的變化,被稱為土壤的空間變異[6]。

應(yīng)用Kriging插值方法,分別為各元素全量選擇擬合狀況最好的變異函數(shù)模型對其進行空間統(tǒng)計分析。模型選擇合適與否的判斷標(biāo)準(zhǔn)為:預(yù)測誤差的均值越接近于0,標(biāo)準(zhǔn)均方根預(yù)測誤差越接近于1,其他值越小時,其模型擬合狀況越好[7]。

表1為應(yīng)用球形、高斯和指數(shù)等模型,對研究區(qū)內(nèi)表層土壤元素全量進行模型擬合所得到的相應(yīng)參數(shù)。其中塊金值C0通常反映由實驗誤差和小于取樣尺度的施肥、作物種植、管理水平等隨機性因素引起的變異,結(jié)構(gòu)方差C是由土壤母質(zhì)、地形、氣候等區(qū)域性因素引起的變異,變程a反映了變量空間自相關(guān)的尺度,塊基比C0/C0+C則反映變量的空間自相關(guān)程度,塊基比越小,說明元素具有越顯著的空間自相關(guān)性。即利用該元素的原始數(shù)據(jù)和半方差函數(shù)的結(jié)構(gòu)性,對未采樣點區(qū)域化變量的取值進行線性無偏估值的效果越理想。

由表1可看出松嫩平原南部表層土壤絕大部分元素全量的變異受人為因素干擾情況較少,而主要是受其土壤母質(zhì)、地形、氣候等自然因素的影響。

3 土壤肥力質(zhì)量評價指標(biāo)體系構(gòu)建

研究區(qū)表層土壤肥力指標(biāo)的篩選、指標(biāo)權(quán)重的賦值、隸屬度函數(shù)的選擇計算以及肥力等別的劃分,均參照中國地質(zhì)調(diào)查局《土地質(zhì)量地球化學(xué)評價技術(shù)要求(試行)》DD2007-XX[8]的標(biāo)準(zhǔn),按照植物生長對土壤元素的不同需求,將所有元素分為必需大量、必需中量、必需微量和有益元素4個指標(biāo)層。同時,考慮到研究區(qū)內(nèi)表層土壤pH值偏高,土壤鹽堿化嚴(yán)重,在評價指標(biāo)的篩選過程中,將pH值對元素有效量的干擾列為篩選評價指標(biāo)的另一要因,對各元素全量進行綜合分析。

3.1 評價指標(biāo)的篩選

按照第二次土壤普查養(yǎng)分分級標(biāo)準(zhǔn)將N,P,K,Corg4種元素全量劃分為六級;采用全國土壤背景值(A層)順序統(tǒng)計量的25%,50%,75%,90%作為分級標(biāo)準(zhǔn),將 Ca,Mg,S,Fe,Mn,Zn,Cu,B,Mo,Cl,Si,Co,Ni,Na共14種元素全量劃分為五級。按照以上分級統(tǒng)計判斷各元素含量豐缺度,結(jié)合pH值對肥力元素有效性的影響,綜合考慮各元素空間自相關(guān)性以及元素間相關(guān)性選取評價指標(biāo)。

表1 黑龍江省松嫩平原南部表層土壤指標(biāo)空間變異特征統(tǒng)計

必需大量元素中P元素全量為嚴(yán)重缺乏,Corg次之;P元素全量具有強烈的空間自相關(guān)性;P元素全量在6≤pH≤7時其有效量最高,最易被植物所吸收,而當(dāng)pH值低于或高于此閾值范圍時,P元素有效量開始不斷降低,不易被植物吸收。研究區(qū)內(nèi)僅有19.61%的P元素處于有效量最佳,即從植物營養(yǎng)供給的角度分析,P元素處于嚴(yán)重缺乏狀態(tài);最后參考各必需大量元素間的相關(guān)關(guān)系,最終確定必需大量元素中P和Corg為參評指標(biāo)。

必需中量元素中S元素全量相對缺乏;Ca,Mg元素全量具有強烈的空間自相關(guān)性,S元素全量具有中等空間自相關(guān)性;僅有不到50%的Ca,Mg元素全量處在6.5≤pH≤8.5的有效量最佳狀態(tài),而S元素有效量則受pH值干擾較小;最后考慮到S與Ca之間相關(guān)性十分低,S與Mg之間相關(guān)性較高,因此,確定S和Ca為參評指標(biāo)。

必需微量元素中Zn,Cu,B,Mo等元素全量缺乏區(qū)面積占到了50%以上,為較嚴(yán)重缺乏;Fe,Mn,Zn,Cu,Cl等元素全量均具有強烈的空間自相關(guān)性,全量B,Mo兩元素則具有中等空間自相關(guān)性;全量的Fe,Mn,Zn,Cu四種元素,只有 19.61%的樣本處于6≤pH≤7的有效量最佳閾值范圍內(nèi),在極強酸的土壤中,其可溶性變大,而使植物受到毒害,在石灰性和堿性土壤中植物又往往發(fā)生缺乏現(xiàn)象,全量的B元素只有不到25%的樣本處于5＜pH≤7的有效量最佳閾值范圍,而全量的Mo元素受pH值影響則較小;處于有效量最大狀態(tài);最后考慮到以上各必需微量元素間相相關(guān)關(guān)系,最終確定參評指標(biāo)為Mo,B,Fe,Cl。

有益元素中只有Co和Ni相對較為缺乏;各元素均具有強烈的空間自相關(guān)性;參考各有益元素間相關(guān)關(guān)系,可知Co,Ni相關(guān)性極高;Si和Na,Co為一般相關(guān),Si與Ni相關(guān)性稍高;Na與Ni,Co的相關(guān)性均較高。最終確定Co,Na為參評指標(biāo)。

通過以上綜合分析,最終選取對研究區(qū)土壤肥力質(zhì)量影響至關(guān)重要的元素作為評價指標(biāo),分別是必需大量元素中P與Corg;必需中量元素中的S與Ca;必需微量元素中的Mo,B,Fe和Cl;有益元素中的Co與Na,共10種元素。

3.2 評價指標(biāo)權(quán)重的確定

結(jié)合研究區(qū)表層土壤元素的豐缺度分布特征及空間變異性,采用經(jīng)驗法和層次分析法對各指標(biāo)進行權(quán)重賦值。將選取的評價指標(biāo)構(gòu)建遞階層次結(jié)構(gòu),參考相關(guān)研究中的指標(biāo)權(quán)重值,由下及上對各判斷矩陣中的評價指標(biāo)進行影響程度的兩兩比較,由判斷矩陣計算得到所比較元素的單層次權(quán)重,并由準(zhǔn)則層權(quán)重和指標(biāo)層權(quán)重計算得到評價指標(biāo)的最終組合權(quán)重值,結(jié)果如表2。

表2 土地地球化學(xué)評價不同層次指標(biāo)權(quán)重

經(jīng)檢驗,中間層要素判斷矩陣的一致性比例CI=0.0123＜0.1;必需大量元素判斷矩陣的一致性比例CI=0.0000＜0.1;必需中量元素判斷矩陣的一致性比例CI=0.0000＜0.1;必需微量元素判斷矩陣的一致性比例CI=0.0226＜0.1;有益元素判斷矩陣的一致性比例CI=0.0000＜0.1。各判斷矩陣均符合一致性檢驗。

3.3 隸屬度函數(shù)計算

首先進行評價指標(biāo)的標(biāo)準(zhǔn)化處理,為使評價的因子具有可比性,需要對評價因子進行無量綱的標(biāo)準(zhǔn)化處理[9]。根據(jù)參評元素與土壤肥力質(zhì)量成正相關(guān)關(guān)系,選擇戒上型隸屬度模型對各評價指標(biāo)值進行標(biāo)準(zhǔn)化,隸屬度函數(shù)曲線如圖2所示。

圖2 松嫩平原南部表層土壤肥力指標(biāo)的隸屬度函數(shù)曲線

相應(yīng)的隸屬度函數(shù)為

式中:U——評估指標(biāo)的上限值;L——評估指標(biāo)的下限值;x——評估指標(biāo)的測定值。

當(dāng)規(guī)定測定值高于U時f(x)=1.0,表明該評價指標(biāo)值能完全滿足植物生長所需肥力;當(dāng)規(guī)定測定值低于L時f(x)=0.1,表明該評價指標(biāo)對植物生長已表現(xiàn)不出肥力。

對各元素全量進行異常剔除后,應(yīng)用累計頻率曲線法按照20%,40%,60%,80%的累積頻率值進行指標(biāo)分級,分別確定20%和80%的累積頻率值為界限值L和U。

3.4 土壤肥力質(zhì)量綜合指數(shù)確定

應(yīng)用隸屬度函數(shù)模型計算結(jié)果和各元素權(quán)重閾值,采用加法模型,對各評價指標(biāo)的實測值進行權(quán)重

式中:P——表層土壤肥力質(zhì)量地球化學(xué)評價的綜合指數(shù);fi——第i個評價指標(biāo)的隸屬函數(shù)值;Ci——第i個評價指標(biāo)的權(quán)重。

為實現(xiàn)土壤肥力質(zhì)量評價結(jié)果的全國可比,參照中國地質(zhì)調(diào)查局《土地質(zhì)量地球化學(xué)評價技術(shù)要求(試行)》DD2007-XX[8]中對土壤肥力質(zhì)量等的劃分標(biāo)準(zhǔn),以0.3和0.7為等別分界值,將表層土壤肥力質(zhì)量分為三等。并與黑龍江省行政區(qū)劃圖疊加,生成1∶200萬的表層土壤肥力質(zhì)量評價等別分布圖,如圖3所示。和隸屬度計算,獲得最終的表層土壤肥力質(zhì)量地球化學(xué)評價綜合指數(shù)。

圖3 松嫩平原南部表層土壤肥力質(zhì)量評價等別分布圖

其中,肥力質(zhì)量為豐富等的土壤占研究區(qū)總面積的19.19%,適量等占64.56%,缺乏等占16.25%。從表層土壤肥力質(zhì)量的分布中可以看出,表層土壤肥力豐富等主要集中在研究區(qū)中部,其他區(qū)域也有零星分布,肥力適量等主要圍繞肥力豐富等分布,集中于研究區(qū)東部區(qū)域內(nèi),肥力缺乏等則較零散地分布于研究區(qū)西側(cè),其中以肇源縣西側(cè)和杜蒙縣分布較為集中。

4 結(jié)論

應(yīng)用地理信息系統(tǒng)軟件的統(tǒng)計分析和空間Kriging插值功能,對研究區(qū)表層土壤元素地球化學(xué)特征進行了分析,了解各元素的含量豐缺度和空間變異性特征,研究表明,P,Corg,Zn,B,Mo等元素為嚴(yán)重缺乏,Si元素含量豐富,其他元素含量較為適中;研究區(qū)內(nèi)各元素均具有自相關(guān)性,其中Mo元素塊基比最大為40.5%,仍具有中等自相關(guān)性,說明研究區(qū)內(nèi)各養(yǎng)分元素的變異受自然因素的影響較大。

將pH值對各元素有效性的影響,一同作為指標(biāo)篩選的判斷依據(jù),避免了只考慮元素全量豐缺度造成的指標(biāo)選擇偏差。同時應(yīng)用地理信息系統(tǒng)與數(shù)學(xué)模型相結(jié)合的方法,對表層土壤肥力質(zhì)量進行評價,使得定量化的評價過程更加直觀,評價結(jié)果更加客觀和具有參考價值。

將松嫩平原南部表層土壤肥力質(zhì)量劃分為3個等別,肥力豐富等占研究區(qū)面積的19.19%,肥力適量等占64.56%,肥力缺乏等占16.25%。

[1]劉占鋒,傅伯杰,劉國華,等.土壤質(zhì)量與土壤質(zhì)量指標(biāo)及其評價[J].生態(tài)學(xué)報,2006,26(3):901-913.

[2]顏雄,張楊珠,劉晶.土壤肥力質(zhì)量評價的研究進展[J].湖南農(nóng)業(yè)科學(xué),2008(5):82-85.

[3]趙穎麗,吳克寧,黃勤,等.鄉(xiāng)級土壤有效態(tài)元素空間變異特征及質(zhì)量評價[J].農(nóng)業(yè)工程學(xué)報,2008,24(1):48-52.

[4]王宗明,張柏,宋開山,等.松嫩平原土地利用變化對區(qū)域生態(tài)系統(tǒng)服務(wù)價值的影響研究[J].中國人口?資源與環(huán)境,2008,18(1):149-154.

[5]中國地質(zhì)調(diào)查局.《多目標(biāo)區(qū)域地球化學(xué)調(diào)查規(guī)范(1∶250000)》(DD2005-01)[S].北京:中國地質(zhì)調(diào)查局.2005:3-19.

[6]李亮亮,依艷麗,凌國鑫,等.地統(tǒng)計學(xué)在土壤空間變異研究中的應(yīng)用[J].土壤通報,2005,36(2):265-268.

[7]陳彥,呂新.基于GIS和地統(tǒng)計學(xué)的土壤養(yǎng)分空間變異特征研究[J].中國農(nóng)學(xué)通報,2005,21(7):289-405.

[8]中國地質(zhì)調(diào)查局.《土地質(zhì)量地球化學(xué)評價技術(shù)要求(試行)》(DD2007-XX)[S].北京:中國地質(zhì)調(diào)查局,2007:8-10.

[9]賈士靖,楊廣林.創(chuàng)業(yè)農(nóng)場農(nóng)田土壤肥力評價[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究,2003,19(2):142-144.