近40年阜康市耕層土壤養分變化特征

高 潔，武紅旗*，李新梅，范燕敏，侯艷娜，楊強軍

(1.新疆農業大學草業與環境科學學院，新疆 烏魯木齊 830052；2.阜康市農業技術推廣中心，新疆 阜康 831500)

【研究意義】土壤養分是土壤肥力的重要標識[1]，是農作物生長的必要條件，其含量直接影響農作物生長和產量，對農業可持續發展有重要意義[2-4]。【前人研究進展】國內外學者針對不同尺度[5-9]、不同地貌類型[10-12]、不同插值方法[13-16]相繼開展了大量研究，取得了豐碩的研究成果。其中，土壤養分的空間變異性已成為土壤科學研究的重要領域，諸多學者利用經典統計與地統計相結合的方法，對土壤養分變化進行研究，該方法彌補了經典統計學空間信息上的不足，具有較好的應用前景[17-19]。同時，對阜康市土壤養分含量變化及空間變異性的研究已有學者進行，馬木提阿吉·烏買爾等[20]對阜康市1982、2012年土壤養分含量變化進行了分析。馬媛[21]等對阜康市土壤微量元素的空間變異性進行了討論。喬娟峰等[22]對阜康荒地土壤有機質的空間變異性進行了研究。但發現研究的時間均較早且跨度較小。隨著糧食需求的增長，耕地不斷擴張、化肥施用不合理等現象造成環境問題日益加劇。土壤養分變化已成為研究的熱點之一。【本研究切入點】然而，現有土壤養分的研究方法多注重定量描述，通過表格進行簡單描述性統計，不能很好的地對長時期土壤養分變化狀況進行可視化表達。而地學信息圖譜，能夠將復雜的現象簡化為空間單元格局圖和時間序列變化過程譜[23]。因此，本文將地學信息圖譜與土壤養分變化研究相結合，對阜康市耕層土壤養分進行整體系統的分析。結合經典統計與地統計的方法，確定最優半變異函數模型，最優搜索方向，提高了插值精度，借鑒地學信息圖譜理論，基于時段首期、中期、末期構建土壤養分變化圖譜，對阜康市近40年土壤養分主要等級變化進行分析，豐富了土壤養分變化研究方法體系，深入挖掘土壤養分變化規律及轉換過程。該方法有助于及時掌握田間養分變化信息、降低生產成本，從而實現農田養分精準管理。【擬解決的關鍵問題】本研究旨在為阜康市農田養分管理、科學施肥、耕地的可持續利用提供理論依據。

1 材料與方法

1.1 研究區概況

阜康市位于天山北麓，準噶爾盆地南緣，地處87°46′～88°44′E，43°45′～45°30′N，海拔450～5445 m，年均氣溫6.7 ℃，年均降雨量205 mm[24]。屬于典型的干旱區綠洲荒漠帶，溫帶大陸性干旱氣候，具有從山地到荒漠的完整垂直地帶性，地貌特征主要分為南部山區、中部平原、北部荒漠3個單元。土壤類型受地形影響也分為3大類型，即山地土壤、平原土壤、沙漠土壤。轄區內共有7條河流，均發源于天山北麓，水資源較豐富。阜康市下轄3個街道，4個鄉。研究區位于阜康市中部耕地，屬平原地區，耕地總面積為55 253.24 hm2[25]，糧食作物以小麥為主，土壤類型以潮土、灰漠土、灌漠土為主。

1.2 資料收集

本研究收集了阜康市3個時期(1982、2010、2018年)耕地土壤有機質、堿解氮、有效磷、有效鉀養分數據。共獲取1982年樣點數據477個，主要來源于阜康市第二次土壤普查數據及土壤類型圖，通過矢量化獲得；2010年通過GPS定點獲得樣點數據275個，主要來源于阜康市測土配方施肥項目；2018年樣點數據共90個，數據來源于農業農村部阜康市耕地質量變更調查項目，點位數據通過GPS定點獲得。影像數據來源于美國地質調查局(United States Geological Survey，USGS)2018年的Landsat-8 OLI影像，時間為9月。其它數據包括阜康市行政區劃圖、1949-2008年《昌吉六十年》[26]、1989-2018年《新疆統計年鑒》[25]。

圖1 農田采樣點分布Fig.1 The farmland sampling point distribution

1.3 數據處理

以阜康市2018年土地利用現狀圖為基礎，獲取耕地邊界。由于采樣數據往往存在異常值，在對數據進行空間結構分析前需對異常值進行剔除，本研究采用經典統計與Voronoi圖結合的方法，對異常值以3倍標準差原則進行剔除，利用SPSS 19軟件對數據進行K-S檢驗、描述性統計、相關性分析等，對不符合正態分布的數據采用Minitab 17軟件進行Box-Cox轉換，獲取最優擬合參數，GS+9.0軟件對數據進行半方差函數分析，獲取最優擬合模型，運用ArcGIS10.2對采樣數據進行空間自相關分析及普通克里格插值，并采用交叉驗證(Cross-Validation)的方法驗證模型精度，依據“全國第二次土壤普查養分分級標準”對阜康市各時期土壤養分進行分級。

1.3.1 土壤養分變化模式圖譜 借鑒鮑文東[27]、呂曉[28]等土地利用變化研究方法對土壤養分進行研究，探討近40年阜康市不同等級土壤養分間相互轉換，利用圖譜融合代碼對土壤養分進行圖譜分析，公式如下：

D=100A+10B+C

式中：D代表新的圖譜單元代碼；A代表1982年阜康市土壤養分的圖譜單元代碼；B代表2010年阜康市土壤養分的圖譜單元代碼；C代表2018年阜康市土壤養分的圖譜單元代碼。

表2 土壤養分數據正態分布檢驗

1.3.2 圖譜單元取值、排序 本研究土壤養分等級跨度為6，故土壤養分圖譜單元取值為1～6。圖譜單元是地學信息圖譜的基本單位，記錄了該單元在某一時期的土壤養分等級，通過對各圖譜單元面積進行排序，篩選出面積最大的前3種圖譜類型進行分析。

2 結果與分析

2.1 土壤養分基本統計特征

對剔除異常值后的數據進行統計(表1)，從均值上看，有機質3個時期含量變化表現為先減少后增加；堿解氮先增加后減少，變化幅度較穩定，在49～50 mg/kg之間上下浮動；有效磷和有效鉀均有不同程度的增加。1982-2018年養分增長率依次為9.33 %、1.43 %、2.40 %、12.71 %。變異系數(Coefficient of Variation，c.v.)用來描述數據的離散、變異等特點，c.v.<10 %變異性較弱，10 %100 %說明數據具有強變異性，本研究3個時期阜康市土壤養分均呈中等強度的空間變異。利用K-S檢驗法對數據進行正態分析，選取符合正態分布的數據，進行地統計分析。

表1 1982-2018年阜康市耕層土壤養分含量描述性統計特征

2.2 土壤養分空間變異性分析

采用半方差函數、Morans’I自相關指數對土壤養分空間結構進行分析。由表3可知，3個時期土壤養分Z值均大于1.96，養分在全局具有空間相關性，Morans’I指數大于0，空間正相關。半方差函數主要擬合模型為指數模型和球狀模型，其中有效磷各時期模型擬合精度最高，R2均大于0.5。塊金系數2010年有機質為9.13 %，屬于強空間變異，說明養分主要受地形、氣候、土壤性狀等結構性因子影響，其余各時期養分塊金系數均在25 %～75 %之間，屬于中等強度的空間變異，表明結構性因子對養分的影響減弱，隨機性因子如施肥、種植結構、灌溉方式等影響增強，變程反映土壤養分空間相關性的范圍大小。

表3 土壤養分空間結構特征

堿解氮、有效磷3個時期塊金系數相差不大，說明各時期養分在一定范圍內變異程度較小。各變量各時期半方差函數結果與Morans’I指數一致，有機質空間相關性先增加后減弱，有效鉀先減弱后增加，堿解氮、有效磷各時期空間相關性逐漸減弱，隨著時間推移，地塊差異性逐漸增加。

2.3 土壤養分插值分析

確定最優插值參數，即選擇最優步長，最優搜索鄰域進行插值。采用交叉驗證方法檢驗精度，選取標準：平均誤差(Mean，ME)和標準預測誤差(Mean Standardized，MSE)2個參數的絕對值越接近0，標準均方根誤差(Root-Mean-Square Standardized，RMSSE)越接近1越好，均方根誤差(Root-Mean-Square，RMSE)和平均標準誤差(Average Standard Error，ASE)值越接近，差值越小，精度越高。通常以均方根誤差(RMSE)和平均標準誤差(ASE)的最小差值作為第一驗證精度，標準預測誤差(MSE)作為第二驗證精度[29]。

插值中步長大小使用ArcGIS 10.2工具中平均最鄰近法求得，再根據實際變程對步長大小進行調整。通過對比不同扇區類型的插值精度，確定最優搜索鄰域。由表4可知，各時期變量最優扇區類型以8扇區、4扇區為主，且交叉驗證精度符合插值要求：均方根誤差與平均標準誤差差值最小，平均誤差和標準預測誤差接近0，標準均方根誤差接近1。

表4 土壤養分插值參數及精度驗證

2.4 各鄉鎮土壤養分含量狀況分析

選擇最優參數對數據進行插值統計，獲取各時期各鄉鎮耕地土壤養分平均值。由表5可知，3個時期有機質含量最高的均為九運街鎮，最低為1982年準東辦事處，含量為15.53 g/kg，3個時期整體有機質均含量大于10 g/kg，處于中等偏上水平。整體趨勢先減后增。堿解氮以滋泥泉子鎮2010年含量最高，為70.00 mg/kg，近40年堿解氮的最低含量變化不大，最高含量呈先增后減趨勢，根據全國第二次土壤養分分級標準，堿解氮近40年土壤養分含量整體小于60 mg/kg，比較缺乏，應在耕作過程中適量補充氮肥。

表5 各鄉鎮土壤養分含量

有效磷1982-2010年含量變化相差不大，都處于中等水平，2018年有效磷含量最高達到32.65 mg/kg，相比之前顯著增加，過多的磷素不僅浪費肥料、對土壤產生有害影響，當磷含量超過一定程度時還會對水體產生危害，應注意磷的減施。有效鉀含量以2010年準東辦事處最高，為315.23 mg/kg，總體呈增長趨勢，含量較豐富。

2.5 土壤養分變化

2.5.1 土壤養分含量分級 依據“全國第二次土壤普查養分分級標準”對阜康市各時期養分含量進行分級，從表6～9可知，3個時期有效鉀含量均集中在一、二、三、四級，含量豐富，能滿足作物生長的需要。堿解氮3個時期含量主要集中在四、五、六級，含量缺乏，應采取相應措施來提高氮肥含量，以確保耕地肥力。有效磷前2個時期含量主要集中在三、四、五級，到了2018年顯著增加，主要集中在二級，磷肥的投入顯著增大。有機質含量變化主要集中在三、四級，總體含量較豐富。

表6 有機質含量分級

2.5.2 養分等級變化 由表10可知，有機質變化主要在三、四級之間轉換，面積最大為17 162.14 hm2，

表7 堿解氮含量分級

表8 有效磷含量分級

表9 有效鉀含量分級

表10 1982-2018年阜康市土壤養分等級変化圖譜特征

變化比率為42.77 %，總體處于中等偏上水平。堿解氮變換主要集中在四、五級之間，面積最大為15 296.38 hm2，變化率為53.54 %，由于阜康市本身堿解氮含量屬中等偏下水平，導致四、五級之間相互轉換較為集中。有效磷的變化1982-2010年間的變化主要集中在三與四級之間的相互轉換，2018年含量猛增至二級，這與描述性統計得出的結論相符。有效鉀在近40年變化主要集中在一、二級，總體來說含量豐富。

為了更好描述土壤養分變化模式，研究土壤養分等級變化轉入轉出，將研究區時段劃分為前期變化型、后期變化型、反復變化型、持續變化型。前期變化型指1982-2010年養分等級發生變化，后期變化型指在2010-2018年養分等級發生變化，反復變化型指養分等級只在2010年發生變化，持續變化型指養分等級在1982-2018持續發生變化。

由圖2可知，有機質以反復變化型為主，主要集中在中部地區，后期變化型次之，持續變化和前期變化所占比例較小，堿解氮同樣以反復變化型為主，主要集中在東部地區，其余所占比率較小，有效磷主要以持續變化和后期變化為主，有效鉀各變化類型相差不大。

圖2 土壤養分含量等級變化模式圖譜Fig.2 Map of soil nutrient content gradation pattern

3 討 論

由變異性、相關性分析可知，研究區養分變化主要受人為因素影響。描述性統計表明：有機質先減后增，堿解氮先增后減，有效磷、有效鉀均增加。這與養分等級變化趨勢相符，和馬木提阿吉·烏買爾等[21]研究結論一致。

20世紀80年代初，改革開放拉開序幕，提出了家庭聯產承包責任制，農民對于養地的觀念不強，結合1982-2018年阜康市化肥施用量統計(圖3)，在全國第二次土壤普查時期(1979-1992年)，灌溉方式主要是大水漫灌，施肥方式落后且施肥量較少，阜康市主要以種植油料作物為主，種植結構單一，有機質含量“減少”。經過幾十年的發展，棉花種植面積有所增加，2005年測土配方施肥政策推行、秸稈還田政策實施，化肥施用量不斷增加，導致有機質含量“增加”。而氮肥是決定產量的重要元素，農民廣施尿素，堿解氮“升高”。隨著測土配方施肥政策的推行，人們對土壤養分逐漸有了正確的認識，氮含量“降低”。新疆測土配方施肥以后，發現磷含量較低，導致磷肥的施用量不斷增大，所以磷在2018年出現“猛增”。由圖3可知，測土配方施肥政策實施后，大大提高了化肥的生產效率，十一五期間，阜康市新開墾的荒地增多，需要投入更多肥料來提高地力，農作物產量不斷提高，導致作物需肥量增加。新開墾的耕地主要在中下部地區，主要以灌漠土、潮土為主，土壤中黏粒含量較高，鉀更容易被吸附，同時，新開墾的耕地本身鉀含量就很高，導致整體鉀肥含量持續“增加”，除此之外，農業機械化水平、灌溉方式、牲畜數量等都是影響土壤養分含量變化的重要因素。

圖3 1982-2018年阜康市化肥施用量Fig.3 Fertilizer use in Fukang form 1982 to 2018

地學信息圖譜與土壤養分變化研究相結合的方法，是在傳統土壤養分研究基礎上，將長時期復雜的土壤養分變化狀況，以圖譜形式直觀展現出來，可以掌握養分在一定時期內的動態轉化過程，更加完整的表達土壤養分變化，有助于了解田間養分變化趨勢，為阜康市耕地土壤養分精準管理提供依據。

4 結 論

本研究以阜康市為例，基于多時期土壤養分數據，揭示了近40年阜康市土壤養分變化規律。

(1)1982-2018年阜康市養分含量總體表現為有機質含量先減少后增加，堿解氮先增加后減少、有效磷、有效鉀含量均呈上升趨勢。

(2)從變異性分析和相關性分析來看，除有機質2010年具有強相關性，其余堿解氮、有效磷、有效鉀各時期都表現為中等強度的相關性，主要受人為因素影響。

(3)經近40年土地利用后，研究區除堿解氮含量較低，其余養分均處于中等或中上水平。養分等級變化：有機質、堿解氮以反復型變化為主，有效磷以持續性變化為主，有效鉀各變化類型相差不大。

(4)研究區應注意增氮，減磷，穩鉀，提高土壤供氮能力；結合作物所需肥力，因地制宜，增加有機肥的投入。