徐 宏， 劉海娥，王世來，楊子凡

（1.甘肅省文縣農(nóng)業(yè)技術(shù)推廣中心，甘肅 文縣 746400； 2.甘肅省農(nóng)業(yè)廣播電視學(xué)校文縣分校，甘肅文縣 746400； 3.甘肅農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，甘肅 蘭州 730070）

土壤養(yǎng)分是植物營養(yǎng)的主要來源，在農(nóng)業(yè)生產(chǎn)中為了滿足農(nóng)作物對營養(yǎng)的需求，在土壤養(yǎng)分不足時需補充施肥，而經(jīng)濟合理地施肥要以土壤養(yǎng)分含量為依據(jù)［1～8］。土壤養(yǎng)分的空間分布是精準農(nóng)業(yè)管理土壤養(yǎng)分和合理施肥的基礎(chǔ)，自20世紀70年代以來，土壤養(yǎng)分空間分布及其變化定量化研究一直是研究關(guān)注的熱點。

文縣位于甘肅省東南麓，地處西秦嶺山脈南秦嶺山脈帶，東經(jīng)104°16′～105°27′，北緯32°35′～33°20′，面積5 000 km2，屬亞熱帶向暖溫帶過渡區(qū)，垂直氣候差異明顯，形成了亞熱帶、溫帶、寒帶疊次鑲嵌的不同氣候類型區(qū)。全境高山、溝谷交錯，現(xiàn)有耕地39 427.48 hm2，人均0.183 hm2。為了解整個縣域境內(nèi)農(nóng)田內(nèi)耕層土壤養(yǎng)分豐缺狀況，指導(dǎo)科學(xué)施肥和服務(wù)精準農(nóng)業(yè)，我們選擇有機質(zhì)、堿解氮、有效磷、速效鉀作為主要指標，根據(jù)空間插值的基本原理［9～11］，使用Arcgis軟件中的kriging插值法，對文縣耕層土壤養(yǎng)分的空間分布進行了研究，現(xiàn)將結(jié)果報道如下。

1 材料與方法

1.1 資料來源

研究資料來源于2008年文縣測土配方施肥2390個實際測量數(shù)據(jù)以及縣域耕地資源管理信息系統(tǒng)空間數(shù)據(jù)庫。

1.2 方法

根據(jù)2 390個實際土壤測試基礎(chǔ)數(shù)據(jù)，利用Arcgis軟件提供的Histogram和QQPlot數(shù)據(jù)分布檢驗方法，分別對各養(yǎng)分化驗數(shù)據(jù)進行檢驗，對檢驗出的離群值、極大極小值進行校正和剔除，以避免對空間插值結(jié)果的干擾，降低誤差。最終分別選用有效數(shù)據(jù)（有機質(zhì)2 384個、堿解氮2 378個、有效磷2 388個和速效鉀2 364個）進行空間插值，通過Arcgis軟件空間分析功能下的kriging插值法獲得文縣耕層土壤的有機質(zhì)、堿解氮、有效磷及速效鉀等主要土壤養(yǎng)分含量的空間分布圖。并根據(jù)分布圖對各鄉(xiāng)鎮(zhèn)土壤養(yǎng)分的豐缺分布進行分析，同時使用Arcgis軟件的分級功能，按照甘肅省耕層土壤養(yǎng)分分級標準對文縣耕層土壤有機質(zhì)、堿解氮、有效磷和速效鉀空間分布進行分析及分級。

kriging法是在空間插值中直接使用擬合半方差圖［9］，先通過由對一組聯(lián)立方程的求解得到權(quán)重：

式中，γ（hij）是已知點i和j間的半方差；γ（hj）是已知點和未知點之間的半方差；λ是拉格朗日系數(shù)。

式中Z0是已知點的值，Wx是與每個已知點關(guān)聯(lián)的權(quán)重，s是用于估算的已知點的數(shù)目［10］。

2 結(jié)果與分析

2.1 耕層土壤有機質(zhì)空間分布及分級

由圖1可以看出，文縣北部、南部等靠近林場保護區(qū)的鄉(xiāng)鎮(zhèn)耕層土壤有機質(zhì)含量高，中部地區(qū)含量較低，耕層土壤有機質(zhì)分布的地帶性明顯。沿江、沿河低階地比旱坡地、旱地梯田高。北部堡子壩鄉(xiāng)、橋頭鄉(xiāng)部分區(qū)域和南部劉家坪鄉(xiāng)有機質(zhì)含量最高，其中堡子壩鄉(xiāng)和劉家坪鄉(xiāng)中部有機質(zhì)含量為25.00～30.00 g/kg；中部城關(guān)鎮(zhèn)、尚德鎮(zhèn)、西部鐵樓鄉(xiāng)等鄉(xiāng)鎮(zhèn)有機質(zhì)含量低，大約在10.00～15.00 g/kg。根據(jù)甘肅省耕層土壤養(yǎng)分分級標準（表1），全縣有機質(zhì)含量由高水平到低水平依次為25.00～30.00、20.00～25.00、15.00～20.00、10.00～15.00 g/kg，分別對應(yīng)甘肅省有機質(zhì)含量的Ⅱ、Ⅲ、Ⅳ、Ⅴ級，Ⅵ級和Ⅶ級在文縣沒有分布，大于30.00 g/kg的Ⅰ級水平在文縣也沒有分布。綜合來看，文縣耕層土壤有機質(zhì)含量較高，屬于有機質(zhì)較豐區(qū)，但有機質(zhì)空間分布不均，部分區(qū)域有機質(zhì)含量較低。

2.2 耕層土壤堿解氮空間分布及分級

文縣耕層土壤堿解氮的空間分布見圖2。大多數(shù)地區(qū)的堿解氮含量在80～100mg/kg，含量大于120mg/kg的區(qū)域在北部、中部和中東部地區(qū)廣泛分布，范壩鄉(xiāng)、碧口鎮(zhèn)、中廟鄉(xiāng)、玉磊鄉(xiāng)和尚德鎮(zhèn)等鄉(xiāng)鎮(zhèn)的堿解氮含量高；土壤堿解氮含量最低的在40～60mg/kg，主要分布在中西部地區(qū)的鐵樓、石雞壩和橋頭等鄉(xiāng)。對照甘肅省耕層土壤養(yǎng)分分級標準（表1），堿解氮含量由高水平到低水平依次為100～150、50～100、＜50mg/kg，分別應(yīng)對甘肅省堿解氮的Ⅴ、Ⅵ、Ⅶ級?？梢娢目h耕層土壤堿解氮含量整體偏低，屬堿解氮欠缺地區(qū)。

圖1 文縣耕層土壤有機質(zhì)空間分布

圖2 文縣耕層土壤堿解氮空間分布

表1 甘肅省耕層土壤養(yǎng)分分級標準（主要養(yǎng)分）

2.3 耕層土壤有效磷空間分布及分級

由圖3可以看出，文縣有效磷含量普遍較低，基本處于8～11mg/kg水平段。含量最高的主要是地處南部的丹堡鄉(xiāng)、劉家坪鄉(xiāng)，有效磷含量為14～17mg/kg；西北部的石雞壩鄉(xiāng)、保護區(qū)，東部的臨江鄉(xiāng)，中部的橋頭鄉(xiāng)等地有效磷含量也相對較高，為12～14mg/kg；中北部的堡子壩鄉(xiāng)、城關(guān)鎮(zhèn)，東部的尚德鎮(zhèn)、玉磊鄉(xiāng)等地有效磷含量低，為3～8 mg/kg。根據(jù)甘肅省耕層土壤養(yǎng)分分級標準（表1），文縣有效磷含量由高水平到低水平依次為15～20、10～15、5～10、≤5 mg/kg，分別對應(yīng)甘肅省有效磷的Ⅳ、Ⅴ、Ⅵ、Ⅶ級，可見文縣耕層土壤有效磷含量低且變幅大，屬有效磷缺乏地區(qū)。

圖3 文縣耕層土壤有效磷空間分布

2.4 耕層土壤速效鉀空間分布分析及分級

圖4 文縣耕層土壤速效鉀空間分布

文縣耕層土壤速效鉀分布見圖4。大多數(shù)地區(qū)的速效鉀含量在150～200mg/kg，處于甘肅省的中等水平；含量小于50mg/kg的速效鉀極低區(qū)域主要集中在城關(guān)鎮(zhèn)周圍，耕層土壤速效鉀含量250～300mg/kg的較高區(qū)，分布在文縣東北部的舍書鄉(xiāng)，南部的鐵樓鄉(xiāng)、保護區(qū)和劉家坪鄉(xiāng)以及東部的中廟鄉(xiāng)等區(qū)域。對照甘肅省耕層土壤養(yǎng)分分級標準（表1），速效鉀含量由高水平到低水平分為5級，依次為250～300、200～250、150～200、100～150、50～100mg/kg，分別對應(yīng)甘肅省速效鉀的Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ級，Ⅶ級在文縣沒有分布，可見文縣耕層土壤富含速效鉀。

3 小結(jié)與討論

1） 分析研究可知，文縣耕層土壤有機質(zhì)含量水平對應(yīng)甘肅省的Ⅱ、Ⅲ、Ⅳ、Ⅴ級，堿解氮對應(yīng)甘肅省的Ⅴ、Ⅵ、Ⅶ級，有效磷對應(yīng)甘肅省的Ⅳ、Ⅴ、Ⅵ、Ⅶ級，速效鉀對應(yīng)甘肅省的Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ級，總體上，文縣耕層土壤有機質(zhì)、速效鉀含量豐富，堿解氮、有效磷缺乏。

2） 從文縣耕層土壤養(yǎng)分分析結(jié)果看，文縣耕層土壤堿解氮含量偏低?？赡芘c氮肥的變異程度和喜氮作物對氮素的大量吸收有較大關(guān)系，因此在施肥過程需密切關(guān)注氮肥的變化。速效鉀含量雖說較高，但對于喜鉀作物和灌區(qū)不應(yīng)該忽視施用鉀肥。文縣境內(nèi)水資源豐富，對于海拔較高、質(zhì)地黏重、土層淺薄、漏水漏肥的耕地，需要采取少量多次的施肥方法。

［1］ 王 祎，蔡立群，張興嘉，等.清水縣耕層土壤主要養(yǎng)分空間分布與變異研究［J］．甘肅農(nóng)業(yè)大學(xué)學(xué)報，2010，47（5）：121-128.

［2］ MISHRA TK，BANERJEE SK．Spatial variability of soil pHand organicmatter under Shorearobusta in lateritic region［J］．Indian Journal of Forestry，1995，18（2）：144-152.

［3］ 楊玉玲，田長彥，盛建東，等．灌淤土壤有機質(zhì)、全量氮磷鉀空間變異性初探［J］．干旱區(qū)農(nóng)業(yè)研究，2002，20（3）：26－30.

［4］ 張淑娟，何 勇，方 慧．基于GPS和GIS的田間土壤特性空間變異性的研究［J］．農(nóng)業(yè)工程學(xué)報，2003，19（2）：39－44.

［5］ 王 軍，傅伯杰，邱 揚，等．黃土高原小流域土壤養(yǎng)分的間異質(zhì)性［J］．生態(tài)學(xué)報，2002，22（8）：1173－1178.

［6］ 高 峻，黃元仿，李保國．農(nóng)田土壤顆粒組成及其剖面分層的空間變異分析［J］．植物營養(yǎng)與肥料學(xué)報，2003，9（2）：151－157.

［7］ 徐 宏，楊子凡，劉海娥，等.基于縣域耕地資源管理信息系統(tǒng)的文縣紋黨適應(yīng)性布局研究［J］．甘肅農(nóng)業(yè)科技，2012（6）：22-24.

［8］ 張 環(huán)，劉淑英.基于綜合指數(shù)評價模型的清水縣土地適應(yīng)性評價［J］．甘肅農(nóng)業(yè)科技，2012（6）：31-33.

［9］ 馮仲科．空間數(shù)據(jù)的最佳內(nèi)插法（kriging法）及其在GIS中應(yīng)用的模型［J］．測繪科學(xué)，1995（3）：22-26.

［10］ Kang-tsung Chang［美］．地理信息系統(tǒng)導(dǎo)論匯［M］．北京：科學(xué)出版社，2004：259-261.

［11］ 楊子凡，蔡立群，張仁陟，等．基于耕地地力評價結(jié)果的景泰縣農(nóng)業(yè)分區(qū)研究［J］．中國農(nóng)學(xué)通報，2012，28（26）：125-130.