西秀區(qū)土壤有效磷分布狀況及變異特征分析

王秋媛

(貴州省安順市西秀區(qū)農業(yè)局土壤肥料站，貴州 安順　561000)

?

■研究報告

西秀區(qū)土壤有效磷分布狀況及變異特征分析

王秋媛

(貴州省安順市西秀區(qū)農業(yè)局土壤肥料站，貴州 安順561000)

以安順市西秀區(qū)2014年采集的669個耕地土樣有效磷測試數據為研究對象，運用地統(tǒng)計學與GIS技術相結合的方法分析土壤有效磷分布狀況及變異特征。結果表明，土壤有效磷含量(mg/kg)的平均值26.4，最小值5.2，最大值89.2，極差84。半方差分析模型擬合比較，指數模型擬合效果最好，塊基比為0.48，說明土壤有效磷的分布狀況是結構性因素(自然因素)和隨機性因素(人為因素)共同作用的結果。普通克里格法內插得到土壤有效磷分布狀況圖，參照第2次土壤普查土壤養(yǎng)分分級標準劃分西秀區(qū)土壤有效磷為4個等級，經面積統(tǒng)計匯總，26.32%的耕地有效磷處于中等水平，68.12%的耕地有效磷處于豐富狀態(tài)，4.85%的耕地有效磷偏高，有效磷缺乏的僅占0.76%。有效磷中等的耕地應穩(wěn)定磷肥施用量，確保土壤有效磷不下降；有效磷豐富的耕地可在配方施肥設計上適量降低磷肥施用量；有效磷偏高的耕地采取隔季施磷的措施，充分發(fā)揮利用土壤中的殘磷。

土壤有效磷；豐缺狀況

磷素(P2O5)是作物生長發(fā)育不可缺少的營養(yǎng)元素之一，作物對磷的需求僅次于氮和鉀，磷對作物吸收養(yǎng)分有很大的促進作用[1]。磷素充足能促進作物體內的物質合成和代謝，作物的產量和品質也得到提高和改善[2]。禾谷類作物缺磷時，會延遲植株的分蘗、抽穗、開花和成熟，在田間可見到水稻“僵苗”、“坐蔸”，玉米的穗禿頂增長，油菜脫莢等現(xiàn)象，但磷肥施用過量也容易引起作物早衰，影響產量和效益[2]。在安順市西秀區(qū)耕地施肥過程中，幾乎所有耕地都存在長期連續(xù)施用磷肥現(xiàn)象。由于對耕地土壤磷素豐缺狀況的認識多停留在單純的土樣測試數據和傳統(tǒng)的施肥經驗上，導致宏觀施肥決策的盲目性和不準確性，使配方施肥方案制定較為粗糙，針對性不強，影響土壤養(yǎng)分管理的進一步提高，極易造成浪費和污染。為此，特運用地統(tǒng)計學與GIS技術相結合的方法對西秀區(qū)土壤有效磷分布狀況和變異特征進行研究,為指導科學合理施用磷肥提供參考。

1　材料與方法

1.1地理概況

安順市西秀區(qū)位于貴州省中部偏西，地處東經105°44′32″～106°21′58″，北緯25°56′30″～26°26′42″之間，東接長順縣，南鄰鎮(zhèn)寧、紫云縣，西連普定縣，北接平壩、織金縣。全區(qū)東西寬61km，南北長59.5km，土地總面積1546平方公里，平均海拔1431m，相對高差593m，年平均氣溫14.2℃，年日照時數1261h，年平均降水量1121mm，無霜期278d，屬亞熱帶季風濕潤氣候，是典型的農業(yè)區(qū)，主要糧油農作物有水稻、玉米、油菜，總播種面積63583.13hm2。全區(qū)耕地土壤分7個土類、25個亞類、86個土屬、143個土種，以黃壤、水稻土、石灰土為主，占全區(qū)耕地總面積的95.5%以上[3]。

1.2土樣采集與測試方法

首先在室內將土地利用現(xiàn)狀圖與土壤圖套合，根據套合單元，在保證采樣點具有典型性和代表性，同時兼顧空間分布均勻性的基礎上，確定采樣點位置，繪制采樣布點圖。在實地采樣點，選擇有代表性的一塊田塊，采用“S”形布點采集7～20個耕層土樣，將土樣混勻，用四分法處理，取1公斤混合土樣，用手持式GPS接收機定位，記錄采樣點的經緯度。及時風干土樣，按照分析要求研磨過篩，裝入樣品瓶。土樣分析按照《測土配方施肥技術規(guī)范》[4]的要求，采用碳酸氫鈉浸提—鉬銻抗比色法測定。采用2014年西秀區(qū)測土配方施肥項目實施中采集的669個土壤樣品有效磷測試數據為研究對象。樣點分布見圖1。

圖1　樣點分布圖Fig.1　Soil sample sites distribution map

1.3分析方法

用Excel的數據分析工具進行描述性統(tǒng)計分析，用統(tǒng)計軟件GS+7.0進行方差分析，用ArcGIS10.0軟件中的統(tǒng)計分析擴展模塊進行普通克里格插值和繪制有效磷分布狀況圖，用空間分析模塊下的區(qū)域統(tǒng)計功能統(tǒng)計不同有效磷含量等級耕地面積。

2　結果與分析

2.1有效磷描述性統(tǒng)計分析

有效磷描述性統(tǒng)計分析(表1)結果看出，有效磷的平均值25.49mg/kg，中位數21.8，眾數9.5，標準差16.56，方差274.28，有效磷含量范圍在5.2～89.2 mg/kg之間，最小值5.2mg /kg，最大值89.2mg /kg，極差84mg /kg。

表1土壤有效磷描述性統(tǒng)計

Table 1　Descriptive statistics of soil effective phosphorus

2.2有效磷正態(tài)分布檢驗

用GS+7.0軟件對數據進行正態(tài)分布檢驗，原始數據趨向正態(tài)分布(圖2)，滿足半方差函數的計算和克里格方法進行土壤養(yǎng)分特性空間分析的前提條件。

圖2　有效磷正態(tài)分布檢驗圖Fig.2　The effective phosphorus normal distribution test

2.3有效磷的半方差分析

土壤性質是一個典型的區(qū)域化變量，土壤性質變化不是隨機的，而是具有高度的空間相關性，對于具有空間結構變量的空間變異分析，采用經典統(tǒng)計學來分析是遠遠不夠的，目前，國內外對土壤空間變異研究最多的是采用地統(tǒng)計學方法[5]。地統(tǒng)計學方法以半方差函數為基本工具，可以對既具有隨機性又具有結構性的各種變量在空間上的分布進行研究，常用的理論模型有線性(Linear)、球狀(Spherical)、指數(Exponential)、高斯(Gussian)4種模型[6]。對數據進行半方差分析，用4種模型分別進行計算擬合,擬合結果見圖3，相關參數見表2。按照決定系數(R2)最大、殘差(RSS)最小的原則選擇最優(yōu)半方差函數模型[7]，經過圖3和表2比較，指數模型擬合效果最好。半方差函數,也稱空間變異函數,是用來描述區(qū)域化變量結構性和隨機性并存這一空間特征而提出的,其中塊金值(C0)、基臺值(C0+C)、變程作為半方差函數的重要參數,用來度量區(qū)域化變量在一定尺度上的空間變異和相關程度[8]。塊金值表示隨機性因素引起的空間異質性；基臺值表示系統(tǒng)內的總變異；塊金值與基臺值的比值(塊基比)表示隨機性因素引起的空間變異性占系統(tǒng)總變異的比例，如果塊基比高,說明隨機性因素引起的空間異質性程度較高；從結構性因素的角度來看,塊基比可以表明系統(tǒng)變量的空間相關性的程度，如果比值<0.25,說明系統(tǒng)具有強烈的空間相關性，如果比值在0.25～0.75之間表明系統(tǒng)具有中等的空間相關性，若>0.75說明系統(tǒng)空間相關性很弱[9]。土壤養(yǎng)分分布是由結構性和隨機性因素共同作用的結果，結構性因素，如土壤類型、母質、地形、氣候等可以導致土壤養(yǎng)分強烈的空間相關性,而隨機性因素如施肥、耕作、灌溉等各種人為活動使土壤養(yǎng)分的空間相關性減弱，朝均一化發(fā)展[10]。從表2指數模型(擬合效果最好)參數可見，有效磷的塊基比為0.48，表明土壤有效磷含量的變化是結構性因素(自然因素)和隨機性因素(人為因素)共同作用的結果。

圖3　有效磷4種模型擬合半方差圖Fig.3　Four semi-variance exponent curve model of soil A-P map

表2有效磷半方差函數理論模型及相關參數

Table 2　The semi-variance theoretical model and correlation parameter of soil A-P

2.4有效磷的空間分布格局

ArcGIS10.0軟件地統(tǒng)計分析擴展模塊提供的空間插值方法——克里格法是利用原始數據和半方差函數的結構性，對未采樣點的區(qū)域化變量進行無偏最優(yōu)估值的一種插值方法，在各種內插方法中普通克里格法效果較好[11，12]。參照第二次土壤普查土壤養(yǎng)分分級標準[13]，結合西秀區(qū)土壤有效磷含量的實際情況，劃分4個豐缺等級：≤10 mg /kg為缺乏，10～20 mg /kg為中等，20～40 mg /kg為豐富，≥40 mg /kg為偏高。對空間插值分布圖分級后與行政區(qū)劃圖疊置，得到以鄉(xiāng)鎮(zhèn)為單位的土壤有效磷含量等值區(qū)圖(圖4)，以定量指標和定位方式直觀反映了轄區(qū)土壤有效磷分布狀況和分布規(guī)律。

用ArcGIS10.0軟件中的空間分析模塊下的區(qū)域統(tǒng)計功能，以耕地資源管理單元圖提取各等值區(qū)面積，并按含量等級匯總地塊面積(表3)。

表3西秀區(qū)土壤有效磷分級評價標準及面積匯總

Table 3　The grading evaluation standard and area statistical of soil A-P in Xixiu

圖4　土壤有效磷分布狀況圖Fig.4　The spatial interpolation distribution map of A-P in soil

土壤養(yǎng)分分布圖是反映土壤養(yǎng)分狀況的應用圖件，是數字農業(yè)的基礎。它能以定量指標和定位方式，直觀地反映調查區(qū)土壤養(yǎng)分含量狀況和分布規(guī)律[14]。從圖4看出，西秀區(qū)耕地土壤有效磷的分布具有明顯的差異，有效磷處于豐富狀態(tài)的耕地占大多數，主要分布在北部、東部、和東南部；有效磷處于中等狀態(tài)的耕地占四分之一，主要分布在轄區(qū)西南部；有效磷偏高的耕地不足5%，主要分布在城區(qū)的東南角及轄區(qū)東部、東北部3個較小區(qū)域；有效磷極缺或缺乏的耕地極少。

3　結論

3.1有效磷含量描述性統(tǒng)計分析結果表明，有效磷(mg/kg)的平均值26.4，中位數22.3，眾數9.5，最小值5.2，最大值89.2，極差84。

3.2應用地統(tǒng)計學方法進行半方差分析，四種模型擬合比較，球狀模型擬合效果最好，塊基比為0.48，說明土壤有效磷的分布狀況是結構性因素(自然因素)和隨機性因素(人為因素)共同作用的結果。

3.3利用普通克里格法進行內插,繪制土壤有效磷分布狀況圖，以定量指標和定位方式直觀表達耕地土壤有效磷豐缺狀況。26.32%的耕地有效磷處于中等水平，68.12%的耕地有效磷處于豐富狀態(tài)，4.85%的耕地有效磷偏高，有效磷缺乏的僅占0.76%。

3.4有效磷中等的耕地應穩(wěn)定磷肥施用量，確保土壤有效磷不下降；有效磷豐富的耕地可在配方施肥設計上適量降低1/4～1/3施用量；有效磷偏高的耕地采取隔季施磷的措施，充分發(fā)揮利用土壤中的殘磷，同時每年或隔年取土分析，監(jiān)測有效磷變化情況，避免盲目施磷造成污染。

[1]周寶庫,張秀英,張喜林,等.化肥施用技術問答[M].第三版.北京:化學工業(yè)出版社,2011.5.60-60.

[2]褚天鐸等.化肥科學使用指南[M].修訂版.北京:金盾出版社,2002.6.66-67.

[3]貴州省土壤普查辦公室.貴州省土壤[M].貴陽:貴州科技出版社,1994.

[4]王凱.測土配方施肥技術規(guī)范(農業(yè)行業(yè)標準)[M].北京:中國農業(yè)出版社,2006.

[5]方元,吳志峰,等.空間分析方法在土壤變異研究中的應用[J].農機化研究，2006,(4): 170-171.

[6]王政權.地統(tǒng)計學及在生態(tài)學中的應用[M].北京:科學出版社,1999.

[7]李中元.貴州喀斯特山區(qū)村級尺度下農田土壤養(yǎng)分空間變異與采樣設計研究[D].貴州大學碩士畢業(yè)論文,2008.

[8]譚萬能,李志安,鄒碧,等.地統(tǒng)計學方法在土壤學中的應用[J].熱帶地理, 2005年12月，第25卷第4期. 307-311.

[9]史利江,鄭麗波,柳云龍.農田土壤養(yǎng)分空間變異特征研究[J].河南農業(yè)大學學報2008, 42(1): 51-56.

[10]張一揚,等.基于地統(tǒng)計與GIS的土壤有機質空間變異性研究—以湖南省主產煙區(qū)為例[J].青島農業(yè)大學學報(自然科學版)2008,25(1): 48-52.

[11]王煒明. 基于GIS的地統(tǒng)計學方法在土壤科學中的應用[J]. 中國農學通報，2007年5月，第23卷第5期. 404-408.

[12]丁園圓,沈潤平,劉莤.基于地統(tǒng)計學的儀征市土壤養(yǎng)分空間分布研究[D].中國氣象學會2007年年會，生態(tài)氣象業(yè)務建設與農業(yè)氣象災害預警分會場論文集 , 2007年.

[13]全國土壤普查辦公室.全國第二次土壤普查暫行技術規(guī)程[M].北京:中國農業(yè)出版社,1979 .59-62.

[14]周敏,王秋兵,賈艷萍. 基于鄉(xiāng)鎮(zhèn)級耕層土壤養(yǎng)分空間變異研究[J].中國土壤與肥料，2009，2:11-16.

Distribution and Variation Characteristic of Available Phosphate Fertilizer in Xixiu District

WANG Qiu-yuan

(XixiuStationofSoilandFertilizer,Xixiu,Guizhou, 561000,China)

A method incorporated Geo statistics and Geo-information system (GIS) was applied to analysis the distribution and variation characters of available phosphate fertilizer (P) of 669 soil samples collected from Xixiu district, Anshun city in 2014. The results indicated that the content of available phosphate range from 5.2mg/kg to 89.2mg/kg, with the average of 26.4mg/kg, and the range was 84; Semi-variance analysis was carried out and the result showed that the exponent curve model was perfect, and the nugget/sill ratio was 0.48, which indicating that the distribution of efficient phosphate was determined by the structural factors (nature factors) and random factors (Artificial factor). The spatial interpolation distribution map of efficient P was analyzed by ordinary Kriging and the soil nutrients in Xixiu district could be divide into four classes by reference the results of the 2ndgeneral survey of soil nutrients, and the medium, abundance, higher and deficient of available phosphate content were 26.32%, 68.12%, 4.85% and 0.76%, respectively among the total area of cultivated land. So, appropriate fertilizer application methods should be emphasized in the areas with different efficient P content to make use of the variability residual P in the soil.

soil；situation

2016-04-25

農業(yè)部、財政部測土配方施肥補貼資金項目

王秋媛(1989-)，女，助理農藝師，從事土壤肥料工作。