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

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

關(guān)鍵詞：速效養(yǎng)分；土壤；空間分布；地震

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

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

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

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

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

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

土壤養(yǎng)分是土壤生產(chǎn)力的基礎(chǔ)，也是土壤重要的屬性之一。土壤養(yǎng)分的高低直接影響著農(nóng)業(yè)生產(chǎn)和環(huán)境安全［１］。土壤速效養(yǎng)分是植物養(yǎng)分的直接來源，因此作物產(chǎn)量與土壤速效養(yǎng)分含量密切相關(guān)［２］。隨著人口壓力的增大和糧食危機的加劇，獲得更多的產(chǎn)量成為農(nóng)業(yè)生產(chǎn)的主要目的，而施肥是調(diào)控作物產(chǎn)量及產(chǎn)量構(gòu)成動態(tài)轉(zhuǎn)化的重要手段［３］，土壤養(yǎng)分狀況又是合理施肥的直接依據(jù)［４］。

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

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

１研究區(qū)域概況與研究方法

１．１研究區(qū)域概況

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

１．２研究方法

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

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

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

１．２．４數(shù)據(jù)處理方法采用ＳＰＳＳ和Ｅｘｃｅｌ對數(shù)據(jù)進行基礎(chǔ)處理，采用ＡｒｃＧＩＳ ９．2進行地統(tǒng)計學(xué)和插值分析。

２結(jié)果與分析

２．１常規(guī)統(tǒng)計分析

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

２．２地統(tǒng)計學(xué)分析

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

２．３插值分析

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

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

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

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

２．４影響因素分析

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

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

３結(jié)論

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

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

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

參考文獻：

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

［２］ 黃國勤，王興祥．施用化肥對農(nóng)業(yè)生態(tài)環(huán)境的負面影響及對策［Ｊ］．生態(tài)環(huán)境，２００４，１３（４）：６５６－６６０．

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

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

［５］ 陸景陵．植物營養(yǎng)學(xué)（上冊）［Ｍ］．北京：中國農(nóng)業(yè)出版社，２００８． ２３－２４．

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

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

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

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

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

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

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

［１３］ 王政權(quán)．地統(tǒng)計學(xué)及在生態(tài)學(xué)中的應(yīng)用［Ｍ］．北京：科學(xué)出版社，１９９９．１－１４９．

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

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

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

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

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

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

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