震后土壤有效氮、磷、鉀空間分布特征

張毅，蔡艷，張世熔，楊喜珍

摘要：運用地理信息系統（ＧＩＳ）和地統計學相結合的方法研究了彭州市新黃村震后土壤堿解氮（ＡＮ）、有效磷（ＡＰ）和速效鉀（ＡＫ）的空間分布特征。結果表明，該區土壤堿解氮、有效磷、速效鉀的含量分別為（１２６．１９±６７．８９） ｍｇ／ｋｇ、（３７．１８±１７．１８） ｍｇ／ｋｇ、（１６５．６９±７９．８０） ｍｇ／ｋｇ；導致堿解氮空間變異的因素以隨機性因素為主，而有效磷和速效鉀為隨機性因素和結構性因素共同影響的結果，但影響程度有所差異。基于ＡｒｃＧＩＳ ９．２平臺的克里格插值分析結果顯示，該區堿解氮含量總體上呈現由北向南逐漸增加的趨勢，且呈條帶狀分布，高值區（＞１６５ ｍｇ／ｋｇ）主要分布在該區的中部和南部，低值區（＜９０ ｍｇ／ｋｇ）零星分布在研究區域內；有效磷含量在總體上呈現以低值區為中心向四周環狀遞增趨勢，且呈兩端高中間低的斑塊狀分布，高值區（＞４４ ｍｇ／ｋｇ）主要分布在北部山坡林地的中部，低值區（＜２０ ｍｇ／ｋｇ）主要分布在低山丘陵與平壩的過渡地帶；速效鉀含量總體上呈斑塊狀分布，且呈現兩端高中間低，西面高東面低的分布特征，以北部山坡林地和西南平壩一角的高值區（＞２２０ ｍｇ／ｋｇ）為中心，分別向南部和東部呈帶狀遞減，并在低山丘陵區東部的旱地以及丘陵向平壩過渡的東部邊緣區出現相對低值（＜１１０ ｍｇ／ｋｇ）分布。顯著性檢驗結果表明，震后壓實對土壤速效養分有顯著或極顯著影響，且水田明顯于旱地；不同土地利用方式對土壤速效養分含量也有顯著影響。

關鍵詞：速效養分；土壤；空間分布；地震

中圖分類號：Ｓ１５２．４文獻標識碼：Ａ文章編號：0439－８114（２０12）14－2948-05

Ｓｐａｔｉａｌ Ｄｉｓｔｒｉｂｕｔｉｏｎ Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ Ｓｏｉｌ Ａｖａｉｌａｂｌｅ Ｎｉｔｒｏｇｅｎ， Ａｖａｉｌａｂｌｅ Ｐｈｏｓｐｈｏｒｕｓ ａｎｄ Ａｖａｉｌａｂｌｅ Ｐｏｔａｓｓｉｕｍ ａｆｔｅｒ Ｅａｒｔｈｑｕａｋｅ：Taking New Huang Village of Pengzhou City as An Example

ＺＨＡＮＧ Ｙｉ，ＣＡＩ Ｙａｎ，ＺＨＡＮＧ Ｓｈｉ－ｒｏｎｇ，ＹＡＮＧ Ｘｉ－ｚｈｅｎ

（Ｃｏｌｌｅｇｅ ｏｆ Ｒｅｓｏｕｒｃｅｓ ａｎｄ Ｅｎｖｉｒｏｎｍｅｎｔ，Ｓｉｃｈｕａｎ Ａｇｒｉｃｕｌｔｕｒａｌ Ｕｎｉｖｅｒｓｉｔｙ， Ｃｈｅｎｇｄｕ ６２５０１４，Ｃｈｉｎａ）

Ａｂｓｔｒａｃｔ： Ｔｈｅ ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｓｏｉｌ ａｖａｉｌａｂｌｅ ｎｉｔｒｏｇｅｎ（ＡＮ）， ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ（ＡＰ） ａｎｄ ａｖａｉｌａｂｌｅ ｐｏｔａｓｓｉｕｍ（ＡＫ） ｉｎ ｎｅｗ ｈｕａｎｇ ｖｉｌｌａｇｅ ｉｎ Ｐｅｎｇｚｈｏｕ ａｆｔｅｒ ｅａｒｔｈｑｕａｋｅ ｗｅｒｅ ａｎａｌｙｚｅｄ ｕｓｉｎｇ ｇｅｏｓｔａｔｉｓｔｉｃｓ ａｎｄ ＧＩＳ． Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｔｈｅ ｃｏｎｔｅｎｔｓ ｏｆ ＡＮ， ＡＰ ａｎｄ ＡＫ ｗｅｒｅ １２６．１９±６７．８９ ｍｇ／ｋｇ， ３７．１８±１７．１８ ｍｇ／ｋｇ ａｎｄ １６５．６９±７９．８０ ｍｇ／ｋｇ ｒｅｓｐｅｃｔｉｖｅｌｙ． Ｆａｃｔｏｒｓ ｃａｕｓｉｎｇ ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ＡＮ ｗｅｒｅ ｍａｉｎｌｙ ｒａｎｄｏｍ ｆａｃｔｏｒｓ， ｗｈｉｃｈ ｏｆ ＡＰ ａｎｄ ＡＫ ｗｅｒｅ ｒａｎｄｏｍ ｆａｃｔｏｒｓ ａｎｄ sｔｒｕｃｔｕｒａｌ ｆａｃｔｏｒｓ ｊｏｉｎｔｌｙ． Ｋｒｉｇｉｎｇ ａｎａｌｙｓｉｓ ｂａｓｅｄ ｏｎ ＡｒｃＧＩＳ ９．２ ｓｈｏｗｅｄ ｔｈａｔ ＡＮ ｉｎ ｔｈｅ ａｒｅａ ｇｅｎｅｒａｌｌｙ ｓｈｏｗｅｄ ａ ｇｒａｄｕａｌ ｉｎｃｒｅａｓｉｎｇ ｔｒｅｎｄ ｆｒｏｍ ｎｏｒｔｈ ｔｏ ｓｏｕｔｈ， ａｎｄ ｓｈｏｗｅｄ ｔｈｅ ｓｔｒｉｐ ｄｉｓｔｒｉｂｕｔｉｏｎ， ａｎｄ ｗｈｏｓｅ ｈｉｇｈ－ｖａｌｕｅ ａｒｅａｓ （＞１６５ ｍｇ／ｋｇ） ｄｉｓｔｒｉｂｕｔｅｄ ｍａｉｎｌｙ ｉｎ ｔｈｅ ｃｅｎｔｒａｌ ａｎｄ ｓｏｕｔｈｅｒｎ， ｌｏｗ－ｖａｌｕｅ ａｒｅａｓ （＜９０ ｍｇ／ｋｇ） ｓｃａｔｔｅｒｅｄ ｉｎ ｔｈｅ ｓｔｕｄｉｅｄ ａｒｅａ； ＡＰ ｐｒｅｓｅｎｔｅｄ ｉｎｃｒｅａｓｉｎｇ ｔｒｅｎｄ ｗｉｔｈ ｌｏｗ ｖａｌｕｅ ａｒｅａ ｃｅｎｔｅｒｅｄ ｉｎ ｇｅｎｅｒａｌ， ａｎｄ ｅｎｄ ｏｆ ｔｈｅ ｓｔｕｄｉｅｄ ａｒｅａｓ ｗｅｒｅ ｈｉｇｈ， ａｎｄ ｈｉｇｈ－ｖａｌｕｅ ａｒｅａｓ （＞４４ ｍｇ／ｋｇ） ｗｅｒｅ ｍａｉｎｌｙ ｄｉｓｔｒｉｂｕｔｅｄ ｉｎ ｔｈｅ ｃｅｎｔｒａｌ ｐａｒｔ ｏｆ ｔｈｅ ｎｏｒｔｈｅｒｎ ｓｌｏｐｅｓ ｏｆ ｗｏｏｄｌａｎｄ， ａｎｄ ｌｏｗ－ｖａｌｕｅ ａｒｅａｓ （＜２０ ｍｇ／ｋｇ） ｗｅｒｅ ｍａｉｎｌｙ ｄｉｓｔｒｉｂｕｔｅｄ ｉｎ ｔｒａｎｓｉｔｉｏｎａｌ ｚｏｎｅ ｆｒｏｍ ｈｉｌｌｙ ｔｏ ｐｌａｉｍ； ＡＫ ｗｅｒｅ ｇｅｎｅｒａｌｌｙ ｐａｔｃｈｙ ｄｉｓｔｒｉｂｕｔｉｏｎ， ａｎｄ ｐｒｅｓｅｎｔｅｄ ｈｉｇｈ ｅｎｄ ａｎｄ ｌｏｗ ｉｎｔｅｒｍｅｄｉａｔｅ ｏｒ ｗｅｓｔ ｈｉｇｈ ａｎｄ ｅａｓｔ ｌｏｗ， ａｎｄ ｒｅｄｕｃｅｄ ｓｔｒｉｐ ｔｏ ｔｈｅ ｓｏｕｔｈ ａｎｄ ｔｈｅ ｅａｓｔ ｒｅｓｐｅｃｔｉｖｅｌｙ ｗｉｔｈ ｎｏｒｔｈｅｒｎ ｗｏｏｄｌａｎｄ ａｎｄ ｓｏｕｔｈｗｅｓｔｅｒｎ ｓｌｏｐｅｓ ａｓ ｈｉｇｈ ｖａｌｕｅ ｃｅｎｔｅｒ（＞２２０ ｍｇ／ｋｇ）． Ｓｉｇｎｉｆｉｃａｎｃｅ ｔｅｓｔ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｔｈｅ ｐｏｓｔ－ｅａｒｔｈｑｕａｋｅ ｃｏｍｐａｃｔｉｏｎ ｈａｄ ｎｏｔａｂｌｅ ｏｒ ｓｉｇｎｉｆｉｃａｎｔ ｉｎｆｌｕｅｎｃｅ ｏｎ ｓｏｉｌ ａｖａｉｌａｂｌｅ ｎｕｔｒｉｅｎｔｓ， ａｎｄ ｔｈｅ ｉｎｆｌｕｅｎｃｅ ｏｆ ｐａｄｄｙ ｆｉｅｌｄ ｗａｓ ｍｏｒｅ ｏｂｖｉｏｕｓ ｔｈａｎ ｔｈａｔ ｏｆ ｄｒｙ ｌａｎｄ； ａｎｄ ｄｉｆｆｅｒｅｎｔ ｌａｎｄ ｕｓｅ ｐａｔｔｅｒｎｓ ａｌｓｏ ｈａｄ ｎｏｔａｂｌｅ ｉｎｆｌｕｅｎｃｅ ｏｎ ｓｏｉｌ ａｖａｉｌａｂｌｅ ｎｕｔｒｉｅｎｔｓ．

Ｋｅｙ ｗｏｒｄｓ： ａｖａｉｌａｂｌｅ ｎｕｔｒｉｅｎｔｓ； ｓｏｉｌ； ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ； ｅａｒｔｈｑｕａｋｅ

土壤養分是土壤生產力的基礎，也是土壤重要的屬性之一。土壤養分的高低直接影響著農業生產和環境安全［１］。土壤速效養分是植物養分的直接來源，因此作物產量與土壤速效養分含量密切相關［２］。隨著人口壓力的增大和糧食危機的加劇，獲得更多的產量成為農業生產的主要目的，而施肥是調控作物產量及產量構成動態轉化的重要手段［３］，土壤養分狀況又是合理施肥的直接依據［４］。

氮、磷、鉀有植物營養三要素之稱。氮和磷均為植物的生命元素［５，６］，鉀為植物的品質元素和抗逆元素。在我國人口迅猛增長，耕地后備資源不足的現實下，只能依靠提高單產來滿足我國人口對糧食和其他農作物的需求，作為植物營養三要素的氮、磷、鉀必然成為眾多學者的研究對象［７－１２］。

土壤具有高度的空間異質性［１３］。２０世紀７０年代后，地統計學被引入土壤科學，加之地理信息系統的發展，許多土壤學家都采用地統計學和ＧＩＳ相結合的方法研究土壤氮、磷、鉀的空間分布特征［１４－１６］。但對地震后土壤速效養分空間分布特征的研究還鮮見報道。因此，以“５.１２”汶川大地震的典型地震災區為研究對象，運用地統計學和ＧＩＳ相結合的方法，研究該區土壤有效氮、磷、鉀空間分布特征，旨在為該區的災后重建工作及農業生產提供參考。

１研究區域概況與研究方法

１．１研究區域概況

彭州市新黃村，距成都主城區５５ ｋｍ，彭州市城區２３ ｋｍ，交通方便，地處東經１０３°４７′２８．９４″－１０３°４８′４６．３９″，北緯３１°９′５７．４３″－３１°１１′２４．７１″，土地面積２．７３５ ｋｍ２，屬湔江流域，地下水位不穩定。受龍門山斷帶及湔江沖積平原影響而具有山、坡、壩復合地貌特征。同時該區也是“５．１２”汶川大地震重災區，具有典型的地震災區代表性。

１．２研究方法

１．２．１土壤樣點設置根據研究區域的實際范圍大小，充分考慮不同土地利用方式，臨時安置點的位置，以及地震中倒塌房屋對農田土壤的影響，結合３Ｓ技術，在該區域內布設了２６個采樣點。采樣點分布如圖１所示。

１．２．２土壤樣品采集與處理結合ＧＰＳ技術在研究區域采樣后，帶回室內，自然風干后，除去非土壤物質，采用四分法取出大約２００ｇ風干樣品，磨細使其全部通過２０目篩，裝袋并做好標記，備用。

１．２．３土壤有效氮、磷、鉀含量的測定均采用常規方法進行測定，即堿解氮采用堿解擴散法測定，有效磷采用０．５ｍｏｌ／Ｌ ＮａＨＣＯ３浸提－鉬銻抗比色法測定，速效鉀采用１．０ ｍｏｌ／Ｌ 中性ＮＨ４ＯＡｃ浸提－火焰光度計法測定［１７］。

１．２．４數據處理方法采用ＳＰＳＳ和Ｅｘｃｅｌ對數據進行基礎處理，采用ＡｒｃＧＩＳ ９．2進行地統計學和插值分析。

２結果與分析

２．１常規統計分析

借助ＳＰＳＳ Ｖ１３．０軟件平臺，對分析結果進行常規統計分析，分析結果見表１。表１結果表明，研究區域土壤堿解氮、有效磷、速效鉀含量均呈正態分布，堿解氮含量為（１２６．１９±６７．８９） ｍｇ／ｋｇ，有效磷含量為（３７．１８±１７．80） ｍｇ／ｋｇ，速效鉀含量為（１６５．６９±７９．８０） ｍｇ／ｋｇ；從偏度上看，有效磷和速效鉀都呈左偏態分布，而堿解氮則呈右偏態分布；從峰度上看，３種有效養分都為負。

２．２地統計學分析

常規統計分析只能說明土壤速效養分的樣點情況，而對其空間分布特征卻難以反映。為更好地反映出該區土壤３種有效養分含量的空間變異特征，借助ＡｒｃＧＩＳ ９．２的地統計學相關組件對區域土壤有效氮、磷、鉀含量進行地統計學分析，結果見表２。從結構性因素的角度來看，模型塊金效應即Ｃ０／（Ｃ０＋Ｃ）反映變量空間變異的來源，比值越大表明人為因素（灌溉、施肥、耕作等）的影響越明顯；反之，結構性因素（母質、氣候、生物、地形等）的影響占主要地位［１８］。如果比值＜２５％，說明系統具有強的空間相關性；如果比值為２５％～７５％，表明系統具有中等的空間相關性；如果比值＞７５％，說明系統空間相關性很弱［１］。土壤養分分布是由結構性因素和隨機性因素共同作用的結果。一般研究認為，結構性因素，如氣候、母質、地形、土壤類型等自然因素可以導致土壤養分具有強的空間相關性，而隨機性因素，如施肥、耕作措施、種植制度等各種人為活動使得土壤養分的空間相關性減弱，朝均一化方向發展［１９］。由表２可知，研究區域土壤堿解氮含量的變程為１ ６７４．７ ｍ，長軸方位角為２８９．２°，約為東南—西北方向，堿解氮空間變異中塊金值與基臺值之比為０．７９（７９％），表明其空間異質性主要由隨機性因素引起，其空間分布的相關性很弱。該區域土壤有效磷的變程為１ ０７９．５ ｍ，長軸方位角為７６．３°，約為東北—西北方向，有效磷含量空間變異中塊金值與基臺值之比為０．３９（３９％），表明其空間異質性主要是隨機性因素和結構性因素共同引起的，其空間分布體現為中等空間相關性。與有效氮、磷比較，該區域土壤速效鉀的變程最小，為９９８．７ ｍ，長軸方位角為２８３．５°，約為東南—西北方向，速效鉀的空間變異中塊金值與基臺值之比為０．３６（３６％），表明其空間異質性主要是隨機性因素和結構性因素共同作用的結果，其空間分布體現為中等空間相關性。

２．３插值分析

為進一步反映出研究區域土壤有效養分的空間分布特征，基于ＡｒｃＧＩＳ ９．２軟件平臺對研究區域土壤有效氮、磷、鉀含量進行插值分析。

通過計算不同間距的半方差，得出該區域土壤堿解氮含量變異函數以高斯模型擬合效果最佳。趨勢分析表明區域內堿解氮含量在北—南方向上的變化趨勢不明顯，呈略微增長趨勢；在東南—西北方向上的變化趨勢同樣不明顯，為略微的遞增趨勢（圖２）。由以上所得半方差理論模型，在ＡｒｃＧＩＳ ９．２軟件平臺上采用普通克里格法（Ｏｒｄｉｎａｒｙ Ｋｒｉｇｉｎｇ）進行最優內插值獲得土壤堿解氮含量的空間分布圖（圖３）。空間分析結果表明，該區域土壤堿解氮含量呈條帶狀分布，高值區（＞１６５ ｍｇ／ｋｇ）主要分布在該區的中部和南部，低值區（＜９０ ｍｇ／ｋｇ）零星分布在研究區域內。全區堿解氮以含量為１１０～１４５ ｍｇ／ｋｇ的區域為主，占全區總面積的４１．５３％；含量為９０～１１０ ｍｇ／ｋｇ的區域次之，占全區總面積的３１．３８％；含量＜９０ ｍｇ／ｋｇ的區域最少，僅占全區總面積的１０．５０％。該區域堿解氮總體上的分布特征為由北向南逐漸增加。

通過計算不同間距的半方差，得出該區域土壤有效磷含量變異函數以指數模型擬合效果最佳。趨勢分析表明區域內有效磷含量分布在北—南方向上呈明顯正“Ｕ”型分布，在東—西方向上呈現明顯淺倒“Ｕ”型分布（圖４）。由以上所得半方差理論模型，在ＡｒｃＧＩＳ ９．２軟件平臺上采用普通克里格法（Ｏｒｄｉｎａｒｙ ｋｒｉｇｉｎｇ）進行最優內插值獲得土壤有效磷含量的空間分布圖（圖５）。圖５結果表明，該區域土壤有效磷含量呈大斑塊狀分布，高值區（＞４４ ｍｇ／ｋｇ）主要分布在北部山坡林地的中部和東南平壩一角，低值區（＜２０ ｍｇ／ｋｇ）主要分布在低山丘陵與平壩的過渡地帶。全區有效磷以含量３５～４４ ｍｇ／ｋｇ的區域為主，占全區總面積的４７．２９％，其次為含量２５～３５ ｍｇ／ｋｇ的區域，占全區總面積的３１．４６％；含量＜２０ ｍｇ／ｋｇ的區域最少，僅占全區總面積的３．９１％。該區有效磷含量在總體上呈現以低值區為中心向四周遞增的趨勢，呈兩端高中間低的特點。

通過計算不同間距的半方差，得出該區土壤速效鉀含量變異函數以指數模型擬合效果最佳。趨勢分析表明區域內土壤速效鉀含量分布在南—北方向呈淺度正“Ｕ”型分布，在東—西方向呈現明顯遞增趨勢（圖６）。由以上所得半方差理論模型，在ＡｒｃＧＩＳ ９．２軟件平臺上采用普通克里格法（Ｏｒｄｉｎａｒｙ Ｋｒｉｇｉｎｇ）進行最優內插值獲得土壤速效鉀含量的空間分布圖（圖７）。結果表明，該區域土壤速效鉀含量以北部山坡林地和西南平壩一角為高值區（＞２２０ ｍｇ／ｋｇ），分別向南部和東部呈條帶狀遞減，并在低山丘陵區東部的旱地以及丘陵向平壩過渡的東部邊緣區出現相對低值（＜１１０ ｍｇ／ｋｇ）分布。全區速效鉀以含量１４０～１９０ ｍｇ／ｋｇ的區域為主，占總面積的４１．０７％；其次為含量＞２２０ ｍｇ／ｋｇ的區域，占全區面積的２０．９４％；含量＜１１０ ｍｇ／ｋｇ的區域最少，僅占全區總面積的約４．７０％。研究區域速效鉀含量在總體上呈現兩端高中間低，西面高東面低的特點。

２．４影響因素分析

２．４．１地震壓實研究區域為“５．１２”地震重災區，在災民安置的過程中造成了大面積的土壤壓實，利用顯著性測驗中的ｔ測驗法［２０］，測驗該區壓實與未壓實區域的土壤養分狀況，結果見表３。由表３可知，水田在壓實后，土壤堿解氮和有效磷與未壓實比較存在極顯著差異，速效鉀呈顯著差異；旱地在壓實后，除速效鉀與未壓實存在極顯著差異外，堿解氮和有效磷都與未壓實差異不顯著。

２．４．２土地利用方式研究區域的土地利用方式主要有水田和旱地。利用顯著性測驗中的ｔ測驗法對該區的旱地和水田進行顯著性測驗。結果表明，堿解氮、有效磷、速效鉀的ｔ值分別為１．８３８、３．５８６**、７．６１０**，表明不同土地利用方式之間，有效磷和速效鉀存在極顯著差異，堿解氮的差異不顯著。

３結論

震后研究區域土壤堿解氮的空間分布特征主要受隨機性因素影響，有效磷和速效鉀是由隨機性因素和結構性因素共同作用的結果。堿解氮的塊金效應大于有效磷，有效磷與速效鉀相當，說明人為因素對堿解氮空間變異性的影響較有效磷與速效鉀強，堿解氮的變程大于速效磷，速效鉀最小。

該區域速效養分的低值區主要分布在中部地區即旱地利用區，其中，堿解氮和有效磷的高值區主要分布在南部即樣板房區，而速效鉀的高值區主要分布在北部即林地。堿解氮以條帶狀分布于研究區域內，由北向南呈增加趨勢；而有效磷和速效鉀則由區域中部的低值區向高值區逐漸過渡。

顯著性測驗顯示，地震壓實對土壤速效養分產生了一定的影響，且壓實對水田的影響較旱地明顯。壓實與未壓實比較，在水田中堿解氮和有效磷呈極顯著差異，速效鉀呈顯著差異，但旱地中除速效鉀呈極顯著差異外，堿解氮和有效磷差異不顯著。不同的土地利用方式對該區域有效磷和速效鉀的影響極顯著，但對堿解氮的影響不顯著。

參考文獻：

［１］ 趙軍，張久明，孟凱，等．地統計學及ＧＩＳ在黑土區域土壤養分空間異質性分析中的應用——以海倫市為例［Ｊ］．水土保持通報，２００４，２４（６）：５３－５７．

［２］ 黃國勤，王興祥．施用化肥對農業生態環境的負面影響及對策［Ｊ］．生態環境，２００４，１３（４）：６５６－６６０．

［３］ 楊志彬，陳兵林，周治國．花鈴期棉田速效養分時空變異特征及對棉花產量品質的影響［Ｊ］．作物學報，２００８，３４（８）：１３９３－１４０２．

［４］ 張樹清，孫小鳳．甘肅農田土壤氮磷鉀養分變化特征［Ｊ］．土壤通報，２００６，３７（１）：１３－１８．

［５］ 陸景陵．植物營養學（上冊）［Ｍ］．北京：中國農業出版社，２００８． ２３－２４．

［６］ 張甘霖，趙玉國，楊金玲，等． 城市土壤環境問題及其研究進展［Ｊ］．土壤學報，２００７，４４（５）：９２５－９３３．

［７］ 高祥照，馬文奇，崔勇，等．我國耕地土壤養分變化與肥料投入狀況［Ｊ］．植物營養與肥料學報，２０００，６（４）：３６３－３６９．

［８］ 肖軍，秦志偉，趙景波．農田土壤化肥污染及對策［Ｊ］．環境保護科學，２００５，31（5）：３２－３４．

［９］ 王伯仁，李冬初，黃晶．紅壤長期肥料定位試驗中土壤磷素肥力的演變［Ｊ］．水土保持學報，２００８，２２（５）：９６－101，141．

［１０］ 秦鐘立，秦松，劉洪斌．貴州省植煙區土壤ｐＨ值和養分空間變異特征的研究［Ｊ］．土壤通報，２００７，38（6）：１０４６－１０５１．

［１１］ 廖育林，鄭圣先，黃建余，等．施鉀對缺鉀稻田土壤鉀肥效應及土壤鉀素狀況的影響［Ｊ］．中國農學通報，２００８，２4（２）：２５５－２６０．

［１２］ 汪濤，楊元合，馬文紅． 中國土壤磷庫的大小、分布及其影響因素［Ｊ］．北京大學學報（自然科學版），２００８，44（6）：９４５－９５２．

［１３］ 王政權．地統計學及在生態學中的應用［Ｍ］．北京：科學出版社，１９９９．１－１４９．

［１４］ 蘇永中，趙哈林，崔建垣．農田沙漠化演變中土壤性狀特征及其空間變異性分析［Ｊ］．土壤學報， ２００４，４１（２）：２１０－２１７．

［１５］ 秦魚生，涂仕華，馮文強，等．成都平原水旱輪作種植下土壤養分特性空間變異研究［Ｊ］．土壤學報，２００８，４５（２）：３５５－３５９．

［１６］ 張玉銘，毛任釗，胡春勝，等．華北太行山前平原農田土壤養分的空間變異性研究［Ｊ］．應用生態學報，２００４，１５（１１）：２０４９－２０５４．

［１７］ 鮑士旦．土壤農化分析［Ｍ］．北京：中國農業出版社，２０００．５６－１０７．

［１８］ 盛建東，肖華，武紅旗．不同取樣尺度農田土壤速效養分空間變異特征初步研究［Ｊ］．干旱地區農業研究，２００５，２３（２）：６３－６７．

［１９］ 郭旭東，傅伯杰，馬克明，等．基于ＧＩＳ和地統計學的土壤養分空間變異特征研究——以河北省遵化市為例［Ｊ］．應用生態學報，２０００，１１（４）：５５７－５６３．

［２０］ 明道緒．田間試驗與統計分析［Ｍ］．北京：科學出版社，２００８．７８－９０．