基于“3S”的廣西典型巖溶區土壤厚度空間格局研究

尹 輝, 李 暉, 蔣忠誠, 楊奇勇

(1.惠州學院 旅游系, 廣東 惠州 516007; 2.廣西師范大學 環境與資源學院,廣西 桂林 541004; 3.中國地質科學院巖溶地質研究所, 廣西 桂林 541004)

基于“3S”的廣西典型巖溶區土壤厚度空間格局研究

尹 輝1, 李 暉2,3, 蔣忠誠3, 楊奇勇3

(1.惠州學院 旅游系, 廣東 惠州 516007; 2.廣西師范大學 環境與資源學院,廣西 桂林 541004; 3.中國地質科學院巖溶地質研究所, 廣西 桂林 541004)

運用地統計學與GIS技術，將室內布點和實地監測相結合，對廣西典型巖溶區果化示范區土壤厚度的空間格局進行了研究。結果表明：果化示范區土壤厚度的空間分布存在正向的全局空間自相關，其空間自相關尺度為2 200 m。土壤厚度以“低—低”型空間集聚為主要分布類型，其空間結構在1～2 200 m區間內為空間集聚；而在2 200～2 800 m區間內為空間孤立。土壤厚度的半方差函數模型符合球狀模型，最小空間變程為500.0 m，該變程可作為區內土壤厚度樣點間隔的布設依據。土壤厚度空間分布不均勻，主要呈條帶狀、弧狀和少量斑狀分布，分布梯度較為明顯。本研究可為果化示范區及相似巖溶峰叢洼地的農業生產、石漠化治理和生態重建提供參考。

土壤; 厚度; 空間自相關; 空間變異; 典型巖溶區

巖溶區土壤具有成土慢、易流失、土層薄且分布不連續等特點[1]。受自然環境和人類活動的影響，我國西南巖溶區的石漠化現象較為嚴重[2-4]。在脆弱的巖溶生態環境中，土層厚度是巖溶區石漠化的主要評估指標之一[5]。近年來，土壤厚度空間格局的相關研究已受到越來越多學者的高度重視。王紹強等[6]基于地統計學和GIS技術，對我國土壤土層厚度的空間變異特征進行了分析；劉憲春[7]以野外采樣調查為依據，輔以統計學和GIS技術,從不同角度探討了東北漫崗黑土層厚度、土體厚度等土壤指標的空間變化規律；謝運杰等[8]運用GIS技術，對黑龍江典型土壤有效土壤厚度進行了實地調研；曾憲勤等[9]運用插釬法，對北京市密云縣山區坡面土壤厚度的空間分布規律進行了調查分析，并對土壤厚度與海拔、坡向等的關系進行了探索；尹亮[10]運用經典統計學與地統計學相結合的方法，對貴州省金沙縣烏箐河小流域土壤厚度的空間變異及其分布格局進行了研究。縱觀土壤厚度空間格局的相關研究成果，在研究方法上已由傳統的量化統計轉向基于“3S”技術的地統計學與經典統計學相結合的空間研究，在技術手段與監測精度上有了質的飛躍。而較其它土壤理化指標的空間格局研究而言，土壤厚度空間格局的相關研究仍較為薄弱，特別是對于土壤淺薄、作物單產低、生態環境脆弱的巖溶區，缺乏從土壤厚度空間關聯與空間異質到樣點布設尺度的系統性研究。而作為巖溶區重要生態環境因子之一的土壤厚度，相關科學問題的提出與實地監測工作的開展，可為區域性農業生產、石漠化治理和生態重建提供科學依據。據此，本文以典型巖溶峰叢洼地——果化示范區為例，從巖溶區土層厚度的空間自相關和空間變異兩個角度，系統性地探討土壤厚度的空間格局，希冀為相似研究區的土壤厚度空間格局評估工作提供借鑒與參考。

1 材料與方法

1.1 研究區概況

果化示范區是我國“十五”科技攻關課題“喀斯特峰叢洼地生態重建技術與示范”和“十一五”科技支撐計劃課題“喀斯特峰叢洼地脆弱生態系統重建技術研究”的石漠化治理與生態重建，以及立體農業生態模式綜合示范區。果化示范區始建于2001年，位于廣西百色市平果縣果化鎮，主要為距果化鎮5 km處的龍何屯、布堯屯和隴堯屯及其各自的耕作區，地理坐標為23°22′28.7″—23°23′40.7″N，107°22′40.8″—107°23′56.9″E，總監測面積約為2.71 km2。示范區屬典型峰叢洼地，海拔為176.4～535.0 m，地處區域構造龍味背斜北翼，北西向斷裂發育，出露的地層主要為石炭系上統(C3)和二疊系下統(P1)。示范區屬于亞熱帶季風性氣候，年均氣溫約為21℃，區內總體降水量豐富，年均降水量多達1 322.3 mm，但時間分配不均，降水年內變化較大。示范區土壤以巖溶石灰土為主，主要發育棕色石灰土和黃棕色石灰土。示范區幾乎不存在園地、牧草地和地表水域，土地利用類型以種植玉米、黃豆的旱作坡耕地為主，植被類型主要為稀疏林地、灌叢和低矮灌草，絕大多數植物是陽生灌木和草本種類，如紅背山麻桿(Alchorneatrewioides)、黃荊條(VitexnegundoLinn)、灰毛漿果楝(Cipadessacinerascens)、類蘆(Neyraudiareynaudiana)、青香茅(Gramineae)、蔓生秀竹(Microsteniumvagans)等。

1.2 圖件來源與處理

平果果化示范區的主要圖件有：果化自然屯區劃圖、1∶1萬5 m等高距地形圖、地質圖、2009年9月的alos高分辨率遙感影像(全色2.5 m分辨率，多光譜10 m分辨率)。以ArcGIS 9.3為操作平臺，完成果化自然屯區劃圖和1∶1萬地形圖的空間配準和矢量化，生成數字化高程圖(DEM)、坡度圖等圖件；運用Envi 4.8遙感軟件，對alos影像進行輻射定標、大氣校正、正射校正、影像融合、影像裁剪等處理，為果化示范區土壤厚度的空間布點作準備。

1.3 樣點選取與實地監測

將上述專題數據提取到ArcGIS 9.3空間處理平臺，運用投影轉換工具，將各圖層的空間參考轉換為高斯—克呂格投影坐標系下的“Beijing_1954_3_Degree_GK_Zone_36”北京54投影坐標系，按照標準格網進行均勻布點，并將準采樣點的坐標輸入GPS，于2011年9—10月進行樣點定位和監測。實際監測過程中，對于待測樣點落在難以到達的險峰、深谷、陡密樹林深處等區域，在待測樣點5～50 m內進行調整，并用GPS記下實際監測樣點的坐標，如不能調整，則該監測樣點剔除。遵照這個原則共得到圖1所示的有效監測樣點191個(各樣點的平均間距約為120 m，樣點密度約為70個/km2)。

野外實地監測應用GPS的精確定位技術，依據土壤樣點的地理坐標，采取野外實地定位方法進行土壤剖面的調查。由于巖溶區土壤較為淺薄，每個土壤剖面均挖取至基巖可見為止。對于土壤厚度不足0.5 m的土壤樣點，用卷尺量取土層厚度，對于土壤厚度大于0.5 m的土壤樣點，在已挖取土壤剖面的基礎上，將鐵釬插入與基巖接觸的土層底部，用皮尺零刻度按壓在與土壤表層接觸的鐵釬部位，抽出鐵釬，將皮尺的另一端拉直，量測插入土壤的鐵釬長度值，作為該樣點的土壤厚度。同一土壤樣點，重復3次測量，將3次測量所得的土壤厚度均值填入土壤樣點剖面調查表，作為土壤厚度空間分析的基礎數據源。

圖1 果化示范區土壤樣點布設圖

1.4數據處理與分析

首先，將野外實地監測的每一個土壤樣點的地理坐標、土壤厚度值和空間環境信息輸入Excel中，導入ArcGIS 9.3并定義為WGS 1984投影坐標系，生成shp格式的土壤厚度樣點布設圖(圖1)，并運用ArcGIS 9.3的投影轉換工具，將shp文件的WGS 1984坐標系轉為“Beijing_1954_3_Degree_GK_Zone_36”北京54投影坐標系。然后，將投影轉換后的shp文件，導入空間自相關軟件Geoda 9.2完成土壤厚度空間權重矩陣、Moran′I指數和空間自相關性的計算，得到土壤厚度Moran散點圖、土壤厚度Moran′I指數圖、土壤厚度空間自相關和Lisa顯著性水平分布圖；再將shp文件的屬性信息導入地統計學軟件GS+7.0，得到土壤厚度半方差函數圖，并擬合出最優變異函數模型，得到相應的塊金值(C0)、基臺值(Sill)、變程(m)等主要參數；最后，打開ArcGIS 9.3的地統計模塊，參照GS+7.0已模擬的模型參數，完成shp格式土壤厚度的空間插值，并根據示范區邊界，在ArcGIS 9.3出圖界面輸出土壤厚空間分布圖。

2 結果與分析

2.1 土壤厚度的空間自相關格局

2.1.1 土壤厚度的全局空間自相關性 由圖2可知，示范區土壤厚度的Moran指數為正向的空間關聯，表明土壤厚度的空間分布存在正向的全局空間自相關。土壤厚度樣點多數位于第三象限，屬于“低—低”型空間集聚類型，表明多數土壤厚度的低值區為附近的低值區所包圍。從土壤厚度的高、低值空間集聚結構來看，“低—低”型空間集聚較為集中，而“高—高”型空間集聚則較為分散。

圖2 土壤厚度Moran散點圖

圖3表明，示范區土壤厚度的Moran′sI標準化指數在1～2 200 m區間內為正值，表明土壤厚度的空間結構在1～2 200 m區間內為空間集聚；而土壤厚度的Moran′s I標準化指數在2 200～2 800 m區間內為負值，表明土壤厚度的空間結構在2 200～2 800 m區間內為空間孤立，從而得到示范區土壤厚度的空間自相關尺度為2 200 m。

圖3 土壤厚度間隔距離與Moran′s I指數

2.1.2 土壤厚度的局部空間自相關性 圖4可見，示范區土壤厚度的空間分布存在明顯的空間集聚區：龍何耕作區北部的土壤厚度呈現出“高—高”型空間集聚區；隴堯西部的土壤厚度呈現出“低—低”型的空間集聚區；土壤厚度空間分布的“高—低”型和“低—高”型空間孤立區較少，其中，土壤厚度的“高—低”型空間孤立區主要分布在“低—低”型空間集聚區附近，而“低─高”型空間孤立區則主要伴隨“高—高”型空間集聚區出現。

2.2 土壤厚度的空間分布格局

表1可見，示范區土壤厚度的半方差函數模型符合球狀模型，其決定系數R2為0.641，擬合殘差僅為0.009，說明土壤厚度半方差函數模型的擬合效果較好，能夠反映示范區土壤厚度的空間結構特征。示范區土壤厚度的C0/Sill值為26.84之間(表1)，表明示范區土壤厚度約25%的變異性是由空間隨機部分引起的，而約75%的變異性是由土壤厚度的結構性變異引起的，其空間結構具有中等強度的空間自相關。

圖4 土壤厚度區域空間關聯和Lisa顯著性水平分布

土壤指標模型塊金值C0基臺值SillC0/Sill變程/mR2殘差RSS土壤厚度球狀0.0510.19026.84490.00.6410.009

從圖5所示的土壤厚度空間插值結果可知，示范區土壤厚度的空間分布不均，主要呈條帶狀、弧狀和少量斑狀分布，梯度分布較為明顯。土壤厚度的空間集聚結構主要表現為：龍何耕作區的北部呈現出高值區空間集聚；隴堯屯的中、西部地區則呈現出低值區空間集聚。從屯落尺度的土壤厚度集聚中心來看，龍何耕作區北部高值區的空間集聚較為明顯，但因龍何屯以北的區域被道路和峰叢隔斷，難以形成完整的閉合型集聚中心；龍何屯北部、布堯屯南部和隴堯屯東南部呈現出次一級的高值集聚中心，而隴堯中偏北的山峰區域則呈現出低值集聚中心。

圖5 土壤厚度空間插值

從土壤厚度分布來看，龍何屯及其耕作區的土壤厚度略高于布堯屯和隴堯屯，主要與龍何屯開展過較大規模的土地整理有關，龍何屯的土地整理工作從坡改梯、平整土地、土壤改造以及相應的農業措施著手，實現了龍何屯保土增厚的預期效果。從布堯屯和隴堯屯的實地調查來看，布堯屯和隴堯屯的土壤厚度主要與地形條件直接相關：地勢較為低平的屯落附近，以及各屯落耕作區的平整區域和低洼地段，土壤厚度相對較大；而山體高度、坡度較大的地段和巖石裸露率較高的區塊，土壤厚度相對較薄。

2.3 土壤厚度的合理取樣間隔

與經典統計學方法相比，地統計學方法在樣點的空間布局、取樣形狀、布局方式和取樣設計的有效性等方面更具優勢。因而，在確保土壤厚度研究精度的前提下，確定示范區土壤厚度的合理取樣間隔，除考慮土壤厚度的統計學特征以外還要考慮其空間結構特性。一般而言，空間變程決定了取樣間隔的大小，其變程越大，取樣間隔可相應加大，其取樣密度可漸趨稀疏，取樣數目可相應減少，從而降低了土壤取樣所需的人力、物力和財力損耗。由半方差函數模擬的變異參數可知，若進行較大尺度的土壤厚度空間監測，土壤厚度的最小空間變程即為土壤空間監測樣點的最適取樣間隔。由表1可見，示范區土壤厚度的最小空間變程490.0 m(約500.0 m)，可作為區內土壤樣點間隔的布設依據。

3 結論與討論

空間自相關分析表明，果化示范區土壤厚度的空間分布存在正向的全局空間自相關，其空間自相關尺度為2 200 m；半方差函數分析表明：果化示范區土壤厚度變異性為26.84%的結構性變異，屬于中等強度的空間自相關。空間插值分析表明，果化示范區土壤厚度空間分布不均，主要呈條帶狀、弧狀和少量斑狀分布，分布梯度較為明顯。其中，龍何耕作區的北部呈現出高值區空間集聚；隴堯屯的中、西部地區則呈現出低值區空間集聚。

相關研究表明，典型巖溶區的土壤理化指標空間變異較大[11-13]，若直接進行大尺度的區域高密度取樣，勢必會產生極大的人力、物力和財力消耗；而大尺度的低密度取樣，又難免會因監測樣點數量的不夠而影響其研究精度[14]。可見，探討土壤樣點的布設尺度是急待系統研究的科學問題。針對土壤理化指標的合理取樣，在考慮土壤指標統計特征的同時，還應考慮其空間結構性，特別應結合地統計學中的空間變程來決定土壤樣點的取樣間隔，以明確取樣形狀和布局方式，制定出合理高效的取樣方案。本文在土壤厚度空間尺度的研究方面，將微尺度、高密度的室內布點和小區域實地監測相結合，應用地統計學與GIS技術，探索空間異質性較大的巖溶區土壤厚度最適取樣間隔：本研究表明，示范區土壤厚度的最小空間變程約為500.0 m，該變程可作為區內土壤厚度樣點間隔的布設依據，相關研究方法可為相似研究區土壤理化指標樣點的空間布設、空間自相關和空間變異研究提供科學借鑒。

[1] 黃金國,魏興琥,王兮之.粵北巖溶山區土地石漠化成因及其生態經濟治理模式[J].水土保持研究,2013,20(4):105-109.

[2] 賴長鴻,覃家作,張文,等.四川省石漠化敏感性評價及其空間分布特征[J].水土保持研究,2013,20(4):99-104.

[3] 于晨曦,陳滸,熊康寧.喀斯特峽谷區干旱災害下的水資源工程性保障研究:以花江石漠化綜合示范區頂壇小流域為例[J].水土保持研究,2013,20(3):305-309.

[4] 張鳳太,張軍以,蘇維詞.喀斯特峰叢洼地農業生產活動的景觀效應初探[J].水土保持研究,2012,19(4):191-195.

[5] 尹輝.我國西南典型喀斯特峰叢洼地土壤理化特征研究[D].北京:中國地質科學院,2012.

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

[7] 劉憲春.東北漫崗黑土壤厚度空間分布規律研究[D].北京：北京師范大學,2005.

[8] 解運杰,劉鳳飛,白建宏,等.基于GIS技術的黑龍江省典型土壤有效土層厚度調查研究[J].水土保持研究,2005,12(6):251-253.

[9] 曾憲勤,劉寶元,劉瑛娜,等.北方石質山區坡面土壤厚度分布特征:以北京市密云縣為例[J].地理研究,2008,27(6):1281-1289.

[10] 尹亮.西南巖溶高原石漠化地區小流域土壤厚度空間分異規律研究[D].長沙:湖南大學,2012.

[11] 尹輝,李暉,蔣忠誠,等.典型巖溶區土壤水分的空間異質性研究[J].江蘇農業科學,2013,41(7):332-336.

[12] 劉艷,尹輝,劉煬,等.典型巖溶坡面土壤水分的空間變異研究[J].廣東農業科學,2013,40(13):156-160.

[13] 尹輝,李暉,蔣忠誠,等.廣西果化巖溶區土壤電導率的空間異質性[J].中國沙漠,2014,34(3):786-794.

[14] 尹輝,蔣忠誠,楊奇勇,等.典型巖溶區土壤pH值的空間異質性研究[J].廣東農業科學,2013,40(11):48-51.

StudyonSpatialPatternofSoilThicknessinTypicalKarstAreainGuangxiBasedon‘3S’Technology

YIN Hui1, LI Hui2,3, JIANG Zhong-cheng3, YANG Qi-yong3

(1.TourismDepartment,HuizhouUniversity,Huizhou,Guangdong516007,China; 2.CollegeofEnvironment&Resources,GuangxiNormalUniversity,Guilin,Guangxi541004,China; 3.InstituteofKarstGeology,CAGS,Guilin,Guangxi541004,China)

By using geo-statistics and GIS, combining indoor sampling point distribution and field monitoring of soil thickness in typical karst area of Guohua demonstration zone in Guangxi, the spatial pattern characteristics of soil thickness were studied. The results showed that the spatial distribution of soil thickness showed obviously positive global spatial autocorrelation. The spatial autocorrelation distance of soil thickness was 2 200 m. The ‘low-low’ spatial agglomeration was the main spatial distribution type, showing two types of spatial structures: the one was spatial agglomeration between 1 m to 2 200 m in spatial distances, the other was spatial isolation between 2 200 m to 2 800 m in spatial distances. Soil thickness of the semivariance function models was well in accord with spherical model. The minimum spatial variable scale of soil thickness was nearly 500.0 m, which could be served as the best interval of soil thickness for sampling point distribution in Guohua demonstration zone. The soil thickness showed uneven banded, arc-shaped and small amounts of porphyritic distribution, and the gradient distribution was obvious. This paper is of great theory and reference significance to guide agricultural production, to manage desertification, and to reconstruct the ecology in typical karst peak cluster-depressions of Guohua demonstration zone and other similar karst areas.

soil; thickness; spatial autocorrelation; spatial variability; typical karst area

2013-12-02

：2014-02-04

國家科技支撐計劃(2010BAE00739-02);廣西自然科學資助項目(2013GXNSFBA019222);巖溶所所控項目(121237128100212);桂林市科學研究與技術開發項目(20110120-6,20120122-2,20120125-1);廣西師范大學教育發展基金會教師成長基金(JS2012009，JS2012012);廣西教育廳科研項目(201106LX070、201203YB028、2013YB040);惠州學院博士科研啟動項目(C513.0206)

尹輝(1983—),男,廣西桂林人,講師,博士,主要研究方向為土壤、生態和“3S技術”。E-mail:yinhui741852963@163.com

李暉(1981—),男,河南南陽人,副教授,博士,主要從事水文水資源與生態環境等方面的研究。E-mail:lh1029@126.com

S159

：A

：1005-3409(2014)06-0025-05