太倉市耕地土壤基本屬性空間分布特征研究

摘要 為掌握江蘇省太倉市耕地土壤基本屬性，于2021年秋季采集區域內88個土壤樣品，樣品經制備、實驗室檢測后，獲取土壤有機質（SOM）、全氮（TN）、全磷（TP）、全鉀（TK）、酸堿度（pH）5項屬性值。基于ArcGIS軟件地統計分析模塊中的普通克里格插值功能，得到土壤屬性的空間分布特征。結果顯示，太倉市SOM、TN、TP和TK分別為25.26、1.48、1.45和10.53 g/kg。在空間分布上，SOM和TN具有一定的相似性，整體上呈“西高東低”的趨勢；TP和TK分別為“北高南低”和“東高西低”。太倉市pH均值為7.51，呈中性偏弱堿性。

關鍵詞 土壤屬性；克里格插值；空間分布；變化趨勢

中圖分類號 S153.6 文獻標識碼 A 文章編號 0517-6611（2024）14-0113-04

doi：10.3969/j.issn.0517-6611.2024.14.025

Study on the Spatial Distribution Characteristics of Basic Properties of Cultivated Soil in Taicang City

ZHANG Hai-dong1，LI Hui-bin2，SHI Yun2 et al

（1.Institute of Agricultural Sciences in Taihu Area of Jiangsu/Suzhou Academy of Agricultural Sciences，Suzhou， Jiangsu 215105;2.Institute of Agricultural Resources and Regional Planning， CAAS，Beijing 100089）

Abstract In order to understand the basic soil properties of cultivated land in Taicang City，Jiangsu Province，this study collected 88 soil samples in the autumn of 2021. After preparation and laboratory testing，five attribute values of soil organic matter （SOM），total nitrogen （TN），total phosphorus （TP），total potassium （TK），and pH were obtained. Based on the ordinary Kriging interpolation function in the statistical analysis module of ArcGIS software，the spatial distribution characteristics of soil attributes were obtained. The results showed that the SOM，TN，TP，and TK in Taicang City were 25.26，1.48，1.45 and 10.53 g/kg，respectively. In terms of spatial distribution，SOM and TN have certain similarities，showing an overall trend of “high in the west and low in the east”;TP and TK are “high in the north and low in the south” and “high in the east and low in the west”，respectively. The average pH value in Taicang City is 7.51，which is neutral to slightly alkaline.

Key words Soil properties;Kriging interpolation;Spatial distribution;Changing trends

基金項目 蘇州市農業科技創新項目（SNG2020072）;2019年吉林聯合基金“玉米精準作業智能決策”；國家自然科學基金—區域創新發展聯合基金（吉林）（U19A2061）。

作者簡介 張海東（1984—），男，江蘇東臺人，副研究員，從事智慧農業研究。*通信作者，研究員，從事農業遙感與智慧農業研究。

收稿日期 2023-08-24；修回日期 2023-09-12

土壤有機質（soil organic matter，SOM）具有提供作物養分、土壤保水保肥等生物、物理和化學作用［1-3］，是土壤肥力的基礎，氮磷鉀是作物生長所必需的養分元素［4］，土壤酸堿度（pH）是影響土壤養分有效性的重要因素之一［5］。研究上述土壤屬性的空間特征，對于監測和管理土壤、制定農業政策等具有重要意義。Liu等［6］用地表動態反饋方法預測了江蘇省中部平原地區土壤有機質含量等屬性的空間分布。龐龍輝等［7］對青海省土壤全氮、有機碳、粉粒和pH進行了預測制圖。杜龍全等［8］基于青海省土系調查的 205個土壤剖面數據，利用隨機森林模型，分別建立了表層（0 ～ 20 cm）土壤全氮（TN）、全鉀（TK）、全磷（TP）與環境因素變量（地形、氣候、植被、遙感數據）之間的定量關系，對青海省土壤養分含量的空間分布進行了預測。筆者以太倉市耕地土壤有機質、全氮、全磷、全鉀、pH為研究對象，通過地統計方法獲取其空間分布特征，為把握太倉市耕地土壤屬性區域特征及相關管理輔助提供技術支撐。

1 材料與方法

1.1 研究區概況

太倉市（120°58′～121°20′E、31°20′～31°45′N）位于江蘇省東南部，長江口南岸，東瀕長江，南臨上海市，西連昆山市，北接常熟市，陸域面積665.96 km2。全境地勢平坦，屬長江三角洲沖積平原，自東北向西南略呈傾斜。屬北亞熱帶南部濕潤氣候區，光照充足，四季分明，年均氣溫16.2 ℃，年均降水量1 173 mm，其中5—9月雨量約占全年雨量的55.7%。《太倉市第三次國土調查主要數據公報》顯示，現有耕地面積為21 627.39 hm2（2019年12月31日）。

1.2 樣品布設與采集

以遙感解譯獲取的太倉地耕地分布圖為底圖，采用網格布點法布設并采集 88個有效土壤樣點（圖1），用 GPS導航定位至樣點，樣品采集深度為0～20 cm，

采用梅花形多點混合法進行取樣，四分法取1 kg樣品，采樣

時間為 2021年11月。采集的土樣經室內風干后，挑揀出土壤中混雜的動植物殘體等，研磨過篩待測。

1.3 檢測內容與方法

土壤有機質（SOM）含量測定采用重鉻酸鉀（K2Cr2O7）氧化滴定法，土壤全氮（TN）測定采用開氏法，土壤全磷（TP）含量測定采用酸溶-鉬銻抗比色法;土壤全鉀（TK）含量測定采用火焰光度法，pH測定采用電位法（水土比為 2.5∶1）［9］。

1.4 數據分析

數據的處理、分析和統計圖制作借助于SPSS 17.0和Origin 9.1軟件進行。考慮到空間插值方法要求數據呈正態分布，因此將數據取對數進行log變換，使得數據在一定程度上更符合正態分布的特征。

1.5 空間特征獲取

普通克里格（ordinary kriging，OK）是常用的插值方法之一，以5項土壤屬性數據為基礎，基于地理信息系統軟件ArcGIS 10.0的地統計分析模塊，采用普通克里格插值方法生成各種屬性的空間分布柵格圖。

2 結果與分析

2.1 耕地土壤屬性

由圖2可知，5項屬性值整體上均為正態分布，其中，TP和pH斜偏度較正態分布具有一定差距，分別呈右偏態和左偏態分布。經log變換后，各屬性值峰度均接近3，偏度均接近0（表1），表明變換后的屬性值更加接近正態分布。

太倉市耕地SOM含量平均值為25.26 g/kg，變幅為10.83～40.23 g/kg，表明在點位上存在較大的異質性。根據江蘇省耕地質量監測指標分級標準，處于2級（較高，20～30 g/kg）以上區間的樣點數為72個，占比高達81.8%，以2級為主，1級（高，>30 g/kg）區間的樣點占比為23.9%（圖2），沒有樣點處于4級（較低，6～10 g/kg）和5級（低，≤6 g/kg）。

TN含量為0.63～2.99 g/kg，平均含量為1.48 g/kg。86.4%的樣點分布于2級（較高，1.5～2.0 g/kg）和3級（中，1.0～1.5 g/kg），兩者樣點數均為38個，各占43.2%。TN和SOM數據擬合結果顯示，兩者間具有較好的線性相關性（圖3）。

TP含量為0.67～3.44 g/kg，平均含量為1.45 g/kg。從樣點屬性值來看，太倉市TP含量較高，處于1級（>1 g/kg）和2級（0.8～1.0>30 g/kg）的樣點數分別為59和25個，總占比達95.5%。

TK含量為6.56～15.68 g/kg，平均含量為10.53 g/kg。太倉市土壤TK含量處于較低水平，主要分布在4級（10～15 g/kg），占比達97.8%。

pH為6.10～8.27，均值為7.51，整體呈中性偏弱堿性。pH處于1級（6.5～7.5）的樣點數27個，2級（5.5～6.5）樣點數5個，3級（7.5～8.5）樣點數量最多達56個。

2.2 空間分布特征

由普通克里格插值結果（圖4）可知，太倉市SOM等5項土壤屬性值在空間上呈不同的分布特征。與已有研究結果類似，太倉市耕層土壤SOM與TN含量具有較好的相關性，在空間上也較為一致，均呈“西高東低”分布。TN含量分布呈“外高內低”的分布特征，縣域外部向內逐漸

降低，北部緊依長江的橫涇鎮最高，沙溪鎮和城廂鎮相對較低。TK含量則呈“東北西低”趨勢，由東部的瀏河鎮和浮橋鎮向西南呈梯度遞減的趨勢，位于太倉西部的沙溪鎮和雙鳳鎮部分區域TK含量高于周邊地區。太倉pH在空間上較為均勻，以中性偏弱堿性為主，最低值出現在中部，但與整體分布差異性較小。

2.3 變化趨勢

與1979年開始的第二次土壤普查數據相比，40多年來，太倉市SOM平均含量提升了2.46 g/kg，TN平均含量基本持平，SOM和TN含量總體上處于中等水平，且增加緩慢。pH整體下降了0.27，TP含量有小幅下降，貯量水平整體較高；TK下降幅度較大，減少近50%，鉀素較為缺乏，應適當增加鉀肥用量（表2）。

3 結論與討論

該研究基于采樣數據進行空間插值，在空間上分析了太倉市SOM、TN、TP、TK和pH的分布特征。結果表明，pH較40年前有所下降，表明在工業污染、長期施用化肥的情況下，耕地土壤存在酸化的風險，可能會降低一些養分的有效性，也會提高部分土壤重金屬元素的有效性，從而增加農產品中有害成分的含量。除pH外，其他4項土壤屬性均受自然因素如成土母質等［10-11］和人為因素如種植作物類型、管理方式等在空間變異較大。同時，通過該研究數據與第二次土壤普查數據的對比，發現經秸稈還田、綠肥還田等措施的實施，SOM含量呈緩慢上升趨勢。TK含量則降幅較大，一定程度上要加強關注。

參考文獻

［1］ 黎鈺鑫，趙小敏，郭熙，等.南昌市近郊耕層土壤有機質空間格局及其影響因素分析［J］.江西農業大學學報，2023，45（3）：749-758.

［2］ 顏祥照，姚艷敏，張霄羽.土壤有機質遙感制圖研究進展與展望［J］.中國農業信息，2019，31（3）：13-26.

［3］ 周榮飛，王海孝.喀什地區耕地土壤養分分析及肥力評價［J］.安徽農業科學，2024，52（6）：159-165，184.

［4］ 張妍娥，楊鵬，馬延龍，等.蘭州市耕地養分綜合評價及其影響因素研究［J］.土壤通報，2023，54（5）：1023-1031.

［5］ 王世豪，徐新良，黃麟，等.1980s—2010s東北農田土壤養分時空變化特征［J］.應用生態學報，2023，34（4）：865-875.

［6］ LIU F，ROSSITER D G，SONG X D，et al. An approach for broad-scale predictive soil properties mapping in low-relief areas based on responses to solar radiation［J］. Soil science society of America journal，2020，84（1）：144-162.

［7］ 龐龍輝，劉峰，趙霞，等. 青海省表層土壤屬性數字制圖［J］. 土壤通報，2019，50（3）：505-513.

［8］ 杜龍全，劉峰，史舟，等. 大面積高寒山區土壤養分空間預測與管理分區［J］.土壤，2022，54（6）：1273-1282.

［9］ 魯如坤.土壤農業化學分析方法［M］.北京：中國農業科技出版社，2000.

［10］ 郭宗祥，左其東，李梅，等.江蘇省太倉市耕地地力調查與質量評價：土壤pH、有機質、全氮、有效磷、速效鉀和CEC的變化［J］.土壤，2007，39（2）：318-321.

［11］ 沈文忠，張緒美，李梅，等.近30年太倉市耕層土壤碳氮含量的時空分布特征［J］.中國農技推廣，2022，38（9）：53-59.