堆龍德慶縣農(nóng)田土壤養(yǎng)分空間變異研究

胡 俊

(拉薩市農(nóng)業(yè)技術推廣總站·西藏拉薩·850000)

堆龍德慶縣農(nóng)田土壤養(yǎng)分空間變異研究

胡 俊*

(拉薩市農(nóng)業(yè)技術推廣總站·西藏拉薩·850000)

為推進西藏測土配方施肥工作深入開展,提升農(nóng)民施肥技術,以堆龍德慶縣為例,采用GIS和地統(tǒng)計學相結合的方法,對其土壤養(yǎng)分的空間變異及分布特征進行了分析。研究結果表明,土壤有機質(zhì)含量的空間變異函數(shù)為高斯模型,全氮、有效磷和速效鉀均為球狀模型,PH為指數(shù)模型。所有指標具有中等的空間相關性,說明其受結構因素和隨機因素共同影響。Kriging插值圖較為直觀地描述了堆龍德慶縣土壤養(yǎng)分的空間分布規(guī)律,并且與目前農(nóng)民的施肥習慣、耕作制度等因素高度相符。

空間分析 GIS 地統(tǒng)計學 Kriging 土壤養(yǎng)分

土壤具有不均一變化特性和高度的空間變異性,對土壤空間變異量化過程的研究具有重要意義[1]。地統(tǒng)計學已經(jīng)被證明是分析土壤特性空間分布特征及其變異規(guī)律的最為有效的方法之一[2]。我國相關科研部門從20世紀80年代以來,運用地統(tǒng)計學結合GIS方法對農(nóng)田土壤養(yǎng)分空間變異特性進行了研究和應用。經(jīng)過研究表明:與經(jīng)典統(tǒng)計學相比,地統(tǒng)計學能更好地揭示隨機變量在空間上的分布、變異和相關特性[3]。西藏自治區(qū)從2006年開始測土配方施肥項目,進行了大量的野外土樣采集和化驗工作,掌握了大量的一手資料。本文擬采用堆龍德慶縣測土配方數(shù)據(jù)庫土樣化驗數(shù)據(jù),運用GIS軟件平臺和地統(tǒng)計學對土壤養(yǎng)分進行空間變異分析,有利于深刻認識農(nóng)業(yè)生產(chǎn)的復雜性和時空變異性,有利于構建循環(huán)高效的農(nóng)田生產(chǎn)系統(tǒng),推動農(nóng)業(yè)的健康、可持續(xù)發(fā)展[4]。

1 材料與方法

1.1 研究區(qū)域概況

堆龍德慶縣位于西藏自治區(qū)中南部,拉薩市的西北邊緣,其地理位置為北緯29°26'～30°39',東經(jīng)90°10'～90°27'之間,總土地面積2672.3平方公里。堆龍德慶縣全境均屬高原溫帶半干旱氣候區(qū),多年平均日照時數(shù)3000小時以上,全年太陽輻射總量186千卡/厘米2,日照百分率為68%。年平均氣溫7.5℃,六月最高氣溫15.4℃,一月平均最低氣溫-2.3℃,≥10℃積溫2124.3℃,年降水量444.8毫米,降水非常集中,水熱同期。年平均相對濕度45%,年蒸發(fā)量為降水量的5倍。種植業(yè)較單一,以春青稞和冬小麥為主,一年一熟制。

1.2 土壤采集和測定

按照農(nóng)業(yè)部《測土配方技術規(guī)范》關于采樣單元的規(guī)定,按照100畝為一個采樣單位的原則,將堆龍德慶縣8.27萬畝耕地劃分為827個采樣單元,在全縣范圍內(nèi)共采集混合土樣827個,采樣深度為0—20cm,并用GPS定位。土壤測定項目包括:有機質(zhì)、全氮、堿解氮、全磷、全鉀、有效磷、速效鉀、水溶態(tài)硼、PH。測試方法如下:PH用電位法、有機質(zhì)用重鉻酸鉀-硫酸溶液-油浴法、有效磷用碳酸氫鈉浸提-鉬銻抗比色法、速效鉀用乙酸銨浸提-火焰光度法、堿解氮用擴散法、全氮用氫氧化鈉熔融-鉬銻抗比色法、全鉀用氫氧化鈉熔融-火焰光度法、水溶態(tài)硼用姜黃素比色法。在Arcgis中導入采樣點坐標,生成采樣點分布圖如下:

1.3 研究方法

1.3.1 描述性統(tǒng)計。

運用SPSS17.0軟件,對有機質(zhì)、全氮、有效磷、速效鉀及PH進行描述性統(tǒng)計,掌握數(shù)據(jù)的基本特征。

1.3.2 正態(tài)性檢驗。

運用SPSS17.0非參數(shù)檢驗K-S檢驗,結合Arcgis地統(tǒng)計分析模塊中的Explore Date中的Histogram、Normal QQPlot及半方差/協(xié)方差云圖對數(shù)據(jù)進行正態(tài)性判斷。通過檢驗,所有指標均符合正態(tài)分布。

1.3.3 地統(tǒng)計分析。

運用Arcgis9.3中的地統(tǒng)計分析模塊,選擇合適的模型,擬合半方差函數(shù),求得相關參數(shù)。采用Kriging插值法對縣域內(nèi)所有指標進行插值,得到養(yǎng)分含量分布圖。

2 結果分析

2.1 土壤養(yǎng)分描述性統(tǒng)計分析

運用傳統(tǒng)統(tǒng)計方法,對堆龍德慶縣農(nóng)田土壤養(yǎng)分進行描述性統(tǒng)計,其統(tǒng)計特征值結果見表1。土壤有機質(zhì)和全氮平均值為23.62和1.45g/kg,變程為3.5-58.4、0.12-3.27g/kg;有效磷和速效救平均含量為10.6、69.63mg/kg,變程為1.8-106.6、20-187.0mg/kg;PH平均7.35,變程5.6-8.7。按照標準(CV＜10%為弱變異性,10%≤CV≤30%為中等變異性,CV＞30%為強變異性[7]),堆龍德慶縣農(nóng)田土壤的PH為弱變異,有機質(zhì)和全氮含量為中等變異,有效磷和速效鉀含量為強變異。其中尤以速效磷的變異程度最大,達到78%。

表1 堆龍德慶縣農(nóng)田土壤養(yǎng)分描述性統(tǒng)計特征值

2.2 土壤養(yǎng)分空間變異特征

塊金效應值可以反映土壤養(yǎng)分的空間相關程度。土壤養(yǎng)分空間變異是由結構性因素和隨機性因素共同作用的結果[3]。塊金效應值表達式[C0/(C0+C)],其中C0為塊金值,表示由試驗誤差和小于實際取樣尺度引起的變異,表示隨機部分的空間異質(zhì)性;(C0+C)為基臺值,表示系統(tǒng)內(nèi)的總變異;塊金效應值即塊金值與基臺值的比值,表示隨機部分引起的空間變異占總變異的比例。若塊金效應值大,說明隨機部分引起的空間異質(zhì)性程度高,而結構性因素作用小;相反,若塊金效應值小,說明隨機部分引起的空間異質(zhì)性程度低,而結構性因素作用大。在實際工作中,若塊金系數(shù)＜25%,則認為具有強烈的空間相關性;若塊金系數(shù)在25%-75%,則認為有中等的空間相關性;若＞75%,則認為空間相關性很弱[8]。

圖1 有機質(zhì)半方差圖

圖2 全氮半方差圖

圖3 有效磷半方差圖

圖4 速效鉀半方差圖

圖5 PH半方差圖

擬合最佳模型,首先用Geostatistical Analysis模塊進行分析,得到各種模型擬合的各種參數(shù)值,在眾多模型參數(shù)中選取Mean參數(shù)最接近零,Root-Mean-square參數(shù)和Average standard Error參數(shù)最小,Root -Mean-Square-standardized參數(shù)最接近1的模型類型,最后確定此模型即為最佳擬合模型。堆龍德慶縣土壤有機質(zhì)用高斯模型擬合,全氮、有效磷、速效鉀用球狀模型擬合,PH用指數(shù)模型擬合。見圖1-5,其擬合模型和相關參數(shù)見表2。

有機質(zhì)、全氮、有效磷、速效鉀、PH的塊金效應值均在27%-75%之間,具有中等的空間相關性,說明其變異受結構因素和施肥、田間管理等隨機因素的共同影響?？臻g自相關性總的情況是:全氮＞有機質(zhì)＞PH＞速效鉀＞有效磷。有效磷的塊金效應值最大,為68.67%,說明人為活動對有效磷的影響較大,削弱了結構性因素造成的強的空間自相關作用。變程的大小為PH＞速效鉀＞有機質(zhì)＞全氮＞有效磷,表明PH空間自相關距離最短,有效磷空間自相關距離最長。

表2 堆龍德慶縣土壤養(yǎng)分理論模型和相關參數(shù)

2.3 土壤養(yǎng)分空間分布

圖6-8是堆龍德慶縣土壤養(yǎng)分含量的空間布局插值圖,顯示了不同級別養(yǎng)分含量的分布狀況。從中可以直觀地看出,土壤有機質(zhì)和全氮含量的空間分布規(guī)律相似,其高值中心出現(xiàn)在北部和中西部地區(qū),低值中心出現(xiàn)在中東部地區(qū)。這種分布一方面說明有機質(zhì)與全氮含量相關程度較高;另一方面與農(nóng)民施肥習慣相符,北部和中西部地區(qū)的農(nóng)區(qū)畜牧業(yè)比重比中東部地區(qū)高,有機肥質(zhì)量和施用量普遍比中東部地區(qū)高,由于長期施用大量的優(yōu)質(zhì)有機肥,造成土壤有機質(zhì)和全氮含量普遍高于中東部地區(qū)的分布格局。有效磷的高值分布于東南部地區(qū),低值區(qū)為中部和北部地區(qū),這與東南部農(nóng)區(qū)農(nóng)民喜歡施用高磷二元復合肥磷酸二銨有關系,與實際相吻合。速效鉀的高值區(qū)分布在中北部地區(qū),低值區(qū)分布在中南部地區(qū),這與該縣的種植制度和施肥習慣緊密相關,一方面堆龍德慶縣長期以來不施用鉀肥,只施用氮肥和磷肥;另一方面中北部地區(qū)長期以來只種植春青稞(生育期為120天左右),中南部地區(qū)不僅種植春青稞,還種植大面積的冬小麥(生育期為300天左右),由于耕作制度的差異,顯然中南部地區(qū)耗鉀量比中北部地區(qū)大,而鉀肥有沒有得到及時補充,造成了目前北高南低的分布格局。PH總的分布規(guī)律是北高南低,南部也出現(xiàn)兩個高值帶,其與氣候條件、施肥習慣、耕作制度、田間灌溉等因素關系復雜,有待進一步分析。

3 結論

通過對堆龍德慶縣土壤養(yǎng)分空間變異及分布特征的研究,得出主要結論如下:

3.1 從經(jīng)典統(tǒng)計學分析結果來看,堆龍德慶縣農(nóng)田土壤的PH為弱變異,有機質(zhì)和全氮含量為中等變異,有效磷和速效鉀含量為強變異。其中尤以速效磷的變異程度最大,達到78%。

3.2 地統(tǒng)計學分析結果表明,堆龍德慶縣農(nóng)田土壤有機質(zhì)適合高斯模型、全氮、有效磷、速效鉀適合球狀模型、PH適合指數(shù)模型所有指標的塊金效應值均在27%-75%之間,具有中等的空間相關性,說明其變異受結構因素和隨機因素的共同影響。

3.3 Kriging插值圖直觀地描述了堆龍德慶縣土壤養(yǎng)分的分布格局。總的來看,土壤有機質(zhì)、全氮、有效磷、速效鉀及PH都呈現(xiàn)出一定的規(guī)律性分布,并且與目前農(nóng)民的施肥習慣、耕作制度等因素高度相符。

[1]張敏.基于GIS和地統(tǒng)計學的土壤養(yǎng)分空間變異研究.河南農(nóng)業(yè)大學碩士學位論文,2010:1.

[2]郭旭東,傅伯杰,馬克明.基于GIS和地統(tǒng)計學的土壤養(yǎng)分空間變異特征研究[J].應用學報,2000.11 (4):557-563.

[3]楊忠華,劉芳,趙澤英,馮延璽.基于GIS和地統(tǒng)計學的農(nóng)田土壤養(yǎng)分空間變異性研究.貴州農(nóng)業(yè)科學,2009,37(9):120-123.

[4]楊玉建.農(nóng)田土壤速效鉀含量的空間變異研究[C].農(nóng)業(yè)信息技術應用研究和科學發(fā)展研討會文集. 2008(7):54-58.

[5]王紹強,朱松麗,周成虎.中國土壤土層厚度的空間變異特征[J].地理研究,2001,20(2):161-169.

[6]李海濱,林忠輝,劉蘇峽.Kriging法在區(qū)域土壤水分估值中的應用[J].地理研究,200 1,20(4):446-452.

[7]薛正平,楊星衛(wèi),段項鎖等.土壤養(yǎng)分空間變異及合理取樣數(shù)研究[J].農(nóng)業(yè)工程學報,2002,(4):6-9.

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

Study on spatial variability of soil nutrients in Duilongdeqing county

Hu jun
(Lhasa City Agricultural Technology Extension Station 850000)

To promote Tibetan soil testing and fertilizer work in depth,fertilization technology to enhance farmers to Duilongdeqin County as an example,the use of GIS and geostatistics method of combining its spatial variability and distribution of soil nutrients were analyzed.The results showed that the spatial variability of soil organic matter content as a function of the Gaussian model,total nitrogen,phosphorus and potassium are spherical model,PH for the exponential model.All indicators with moderate spatial correlation,indicating that structural factors and random factors affected by both.Kriging interpolation maps more intuitive description of the distribution Duilongdeqin County soil nutrients in space,and with the current farmers'fertilizer habits,and other factors cropping system highly consistent.

spatial analysis;GIS;geostatistics;Kriging;soil nutrient

胡俊(1978-),男,高級農(nóng)藝師。主要從事農(nóng)業(yè)科研與技術推廣工作。項目來源:西藏測土配方施肥項目。