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

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

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

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

中圖分類(lèi)號(hào)：Ｓ１５２．４文獻(xiàn)標(biāo)識(shí)碼：Ａ文章編號(hào)：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)分的直接來(lái)源，因此作物產(chǎn)量與土壤速效養(yǎng)分含量密切相關(guān)［２］。隨著人口壓力的增大和糧食危機(jī)的加劇，獲得更多的產(chǎn)量成為農(nóng)業(yè)生產(chǎn)的主要目的，而施肥是調(diào)控作物產(chǎn)量及產(chǎn)量構(gòu)成動(dòng)態(tài)轉(zhuǎn)化的重要手段［３］，土壤養(yǎng)分狀況又是合理施肥的直接依據(jù)［４］。

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

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

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

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

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

１．２研究方法

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

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

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

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

２結(jié)果與分析

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

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

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

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

２．３插值分析

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

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

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

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

２．４影響因素分析

２．４．１地震壓實(shí)研究區(qū)域?yàn)椤埃担保病钡卣鹬貫?zāi)區(qū)，在災(zāi)民安置的過(guò)程中造成了大面積的土壤壓實(shí)，利用顯著性測(cè)驗(yàn)中的ｔ測(cè)驗(yàn)法［２０］，測(cè)驗(yàn)該區(qū)壓實(shí)與未壓實(shí)區(qū)域的土壤養(yǎng)分狀況，結(jié)果見(jiàn)表３。由表３可知，水田在壓實(shí)后，土壤堿解氮和有效磷與未壓實(shí)比較存在極顯著差異，速效鉀呈顯著差異；旱地在壓實(shí)后，除速效鉀與未壓實(shí)存在極顯著差異外，堿解氮和有效磷都與未壓實(shí)差異不顯著。

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

３結(jié)論

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

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

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

參考文獻(xiàn)：

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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