舒 田,熊康寧,陳麗莎,肖 杰

(1.貴州師范大學(xué) 喀斯特研究院/國(guó)家喀斯特石漠化防治工程技術(shù)研究中心，貴州 貴陽(yáng) 550001；2.貴州省農(nóng)業(yè)科學(xué)院 科技信息研究所，貴州 貴陽(yáng) 550006)

【研究意義】喀斯特石漠化是指在熱帶亞熱帶暖溫帶濕潤(rùn)半濕潤(rùn)氣候條件的喀斯特環(huán)境背景下，受人類(lèi)不合理社會(huì)經(jīng)濟(jì)活動(dòng)的干擾破壞所造成的土壤嚴(yán)重侵蝕，基巖大面積出露，土地生產(chǎn)力嚴(yán)重下降，地表出現(xiàn)類(lèi)似荒漠景觀的土地退化過(guò)程[1]。土地生產(chǎn)力下降最直接的因素就是土壤肥力的下降及養(yǎng)分的缺失。土壤養(yǎng)分狀況與石漠化過(guò)程關(guān)系緊密，不同石漠化等級(jí)下土壤的養(yǎng)分含量不同[2]。【前人研究進(jìn)展】土壤有機(jī)質(zhì)是指存在于土壤中的含碳有機(jī)物質(zhì)，是植物生長(zhǎng)所需營(yíng)養(yǎng)元素的重要來(lái)源。全氮是土壤各種形態(tài)氮素含量之和，是植物生長(zhǎng)發(fā)育的必需營(yíng)養(yǎng)元素。揭示土壤有機(jī)質(zhì)和全氮的空間變異規(guī)律和掌握其分布狀況，對(duì)農(nóng)業(yè)生產(chǎn)和環(huán)境管理具有重要意義[3]。因此，了解石漠化地區(qū)土壤有機(jī)質(zhì)和全氮分布狀況，不僅可以反映植被恢復(fù)模式對(duì)土壤的響應(yīng)，也可為土壤培肥提供參考依據(jù)[4]，同時(shí)能夠因地制宜地選擇作物，發(fā)展地方特色生態(tài)產(chǎn)業(yè)，為石漠化治理和生態(tài)恢復(fù)提供指導(dǎo)。封山育林、退耕還林還草、林草種植和自然恢復(fù)等生態(tài)治理措施，對(duì)石漠化地區(qū)土壤養(yǎng)分提升具有重要作用。【本研究切入點(diǎn)】通過(guò)對(duì)2010和2015年關(guān)嶺花江示范區(qū)的石漠化等級(jí)遙感解譯，利用土地石漠化過(guò)程的定點(diǎn)采樣監(jiān)測(cè)數(shù)據(jù)并進(jìn)行Kriging插值?！緮M解決的關(guān)鍵問(wèn)題】分析土壤有機(jī)質(zhì)和全氮變異特征和不同石漠化等級(jí)土壤有機(jī)質(zhì)和全氮變異規(guī)律，從而確定土壤有機(jī)質(zhì)和全氮在石漠化作用過(guò)程和生態(tài)治理過(guò)程中的虧缺及變異，為貴州石漠化山區(qū)土壤資源有效利用、水肥保持、植被重建及生態(tài)修復(fù)提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

選擇關(guān)嶺-貞豐花江喀斯特高原峽谷中-強(qiáng)度石漠化綜合治理示范區(qū)(以下簡(jiǎn)稱(chēng)“花江示范區(qū)”)為研究區(qū)。花江示范區(qū)位于貴州省安順市關(guān)嶺縣與黔西南州貞豐縣交界處的北盤(pán)江峽谷花江段，位置105°36′30″～105°46′30″E、25°39′13″～25°41′00″N，屬亞熱帶濕潤(rùn)季風(fēng)氣候，土壤以黃壤、黃色石灰土為主，土地利用以林地和耕地為主，植被類(lèi)型多以次生灌木林為主，區(qū)域內(nèi)水土流失以微度和輕度為主，多年平均降雨量1100 mm，年平均氣溫18.4 ℃，海拔430～1370 m，整體地勢(shì)陡峭，是典型的喀斯特峽谷強(qiáng)度石漠化地區(qū)，其石漠化等級(jí)基本情況見(jiàn)表1。示范區(qū)共涉及3個(gè)鄉(xiāng)鎮(zhèn)、8個(gè)行政村，50個(gè)村民組，包括貞豐縣北盤(pán)江鎮(zhèn)的水淹壩村、查爾巖村、銀洞灣村，關(guān)嶺縣板貴鄉(xiāng)的木工村、壩山村、三家寨村、孔落箐村和花江鎮(zhèn)的五里村[5]，示范區(qū)總面積5164.07 hm2，喀斯特面積占示范區(qū)總面積87.93 %。區(qū)域內(nèi)涉及法郎、木工和頂壇3個(gè)小流域，2010年總?cè)丝?612人，其中農(nóng)業(yè)人口占89.7 %，農(nóng)民人均純收入5390元，人口密度150人/km2，產(chǎn)業(yè)結(jié)構(gòu)仍以第一產(chǎn)業(yè)為主。

表1 研究區(qū)不同石漠化等級(jí)的基本情況Table 1 Basic information of different grade rocky desertification in study area

表2 花江示范區(qū)不同石漠化等級(jí)的分類(lèi)標(biāo)準(zhǔn)Table 2 Grading standard of different rocky desertification in Huajiang demonstration area

1.2 數(shù)據(jù)來(lái)源

(1) 遙感影像?；ń痉秴^(qū)石漠化解譯遙感影像，采用成像時(shí)間為2010年ALOS PRISM(全色波段空間分辨率2.5 m)與ALOS AVNIR-2影像(多光譜波段空間分辨率10 m)，成像時(shí)間為2015年的資源三號(hào)衛(wèi)星數(shù)據(jù)(全色波段空間分辨率2.1 m，多光譜波段空間分辨率5.8 m)，分別對(duì)全色波段和多光譜波段進(jìn)行融合，從而保證示范區(qū)石漠化解譯影像數(shù)據(jù)空間分辨率的一致性。

(2) 有機(jī)質(zhì)和全氮含量?；ń痉秴^(qū)土壤有機(jī)質(zhì)和全氮含量分別采用2010年貞豐和關(guān)嶺縣測(cè)土配方施肥數(shù)據(jù)和2015年7月針對(duì)4種不同石漠化等級(jí)實(shí)施的治理恢復(fù)模式地取樣后的室內(nèi)土壤養(yǎng)分測(cè)定數(shù)據(jù)進(jìn)行Kriging插值得到，采樣深度為0～20 cm表層土，自然風(fēng)干后剔除石塊和雜質(zhì)后研磨。土壤有機(jī)質(zhì)和全氮分別采用油浴加熱重鉻酸鉀氧化容量法和半微量凱氏蒸餾法測(cè)定，各指標(biāo)進(jìn)行3次平行測(cè)定，取平均值。

1.3 研究方法

1.3.1 花江示范區(qū)石漠化等級(jí)分類(lèi)標(biāo)準(zhǔn)的確定及其等級(jí)面積統(tǒng)計(jì) 將遙感數(shù)據(jù)在ENVI5.3中對(duì)ALOS和資源三號(hào)衛(wèi)星數(shù)據(jù)進(jìn)行全色和多光譜波段的融合，基于“3S”技術(shù)，并結(jié)合花江示范區(qū)坡度、巖性條件、基巖裸露率、土壤信息、植被覆蓋度等作為主要判讀因子，輔以示范區(qū)內(nèi)平均降雨量、土層厚度、水土流失狀況、土地利用現(xiàn)狀等因子，參照文獻(xiàn)[6]的分級(jí)方法確定示范區(qū)石漠化等級(jí)的分類(lèi)標(biāo)準(zhǔn)(表2)，基于ArcGIS的Spatial Analyst和疊置對(duì)研究區(qū)2010-2015年石漠化等級(jí)進(jìn)行圖斑判別、定級(jí)和不同石漠化等級(jí)面積統(tǒng)計(jì)。通過(guò)空間分析，研究監(jiān)測(cè)時(shí)段各個(gè)圖斑石漠化等級(jí)的演變情況，分析石漠化的內(nèi)部變化及區(qū)域變化，統(tǒng)計(jì)分析各圖斑石漠化的轉(zhuǎn)移規(guī)律。

1.3.2 花江示范區(qū)土壤等級(jí)的確定 土壤有機(jī)質(zhì)和全氮的變異受成土母質(zhì)、氣候、土地利用方式/覆蓋變化、耕作與施肥方式、栽培歷史等自然和人為因素的影響，土壤養(yǎng)分具有高度的空間異質(zhì)性和變異性[7]。研究參照貴州省第二次土壤普查分級(jí)標(biāo)準(zhǔn)(表3)[8]確定花江示范區(qū)的土壤等級(jí)。

表3 貴州省第二次土壤普查的分級(jí)標(biāo)準(zhǔn)Table 3 Grading standard of second soil general survey in Guizhou

注：各養(yǎng)分含量范圍不含最大值。
Note: Each nutrient content range does not contain the maximum value.

圖1 2010-2015年花江示范區(qū)喀斯特石漠化的等級(jí)分布Fig.1 Grades distribution of karst rocky desertification of Huajiang demonstration area during 2010-2015

1.3.3 石漠化動(dòng)態(tài)度 采用土地利用動(dòng)態(tài)度方法[9]研究區(qū)域內(nèi)石漠化變化的速度，計(jì)算花江示范區(qū)不同等級(jí)石漠化動(dòng)態(tài)度。動(dòng)態(tài)度指數(shù)綜合考慮了研究時(shí)段區(qū)域石漠化等級(jí)類(lèi)型間的轉(zhuǎn)移，用以反映區(qū)域石漠化變化的劇烈程度和變化速度，便于在不同空間尺度上找出石漠化變化的熱點(diǎn)區(qū)域。單一石漠化動(dòng)態(tài)度表達(dá)某一研究區(qū)一定時(shí)間范圍內(nèi)某種石漠化等級(jí)類(lèi)型的數(shù)量變化情況，可用以下分式表示：

式中，K來(lái)描述研究區(qū)某等級(jí)石漠化的變化即石漠化動(dòng)態(tài)度，Ua、Ub分別為各等級(jí)研究期初及研究期末的石漠化面積，T為研究時(shí)段長(zhǎng)，當(dāng)T的時(shí)段設(shè)定為年時(shí)，K值就是研究區(qū)某等級(jí)石漠化年動(dòng)態(tài)度。

1.4 數(shù)據(jù)處理

采用SPSS19.0和Excel 2017對(duì)數(shù)據(jù)進(jìn)行整理，計(jì)算石漠化土壤有機(jī)質(zhì)和全氮含量的平均值、標(biāo)準(zhǔn)偏差和變異系數(shù)，通過(guò)制作圖表對(duì)不同石漠化等級(jí)與土壤有機(jī)質(zhì)和全氮數(shù)據(jù)進(jìn)行相關(guān)性分析。

2 結(jié)果與分析

2.1 2010-2015年花江示范區(qū)石漠化等級(jí)的面積變化

從表4和圖1可知，花江示范區(qū)內(nèi)喀斯特面積4540.52 hm2，占示范區(qū)總面積的87.93 %。經(jīng)過(guò)幾年的石漠化治理，無(wú)石漠化面積2010年為859.57 hm2，2015年增至1066.30 hm2，5年時(shí)間增加206.73 hm2，年動(dòng)態(tài)度為4.81 %，變化貢獻(xiàn)主要來(lái)自示范區(qū)西北部的五里村以及南部的水淹壩村。潛在石漠化面積2015年較2010年增加42.40 hm2，年動(dòng)態(tài)度為0.94 %，變化量不大。輕度石漠化面積2010年為1490.14 hm2，2015年增至1497.96 hm2，期間增加7.83 hm2。強(qiáng)度石漠化面積2010年為628.09 hm2，減少到2015年為500.81 hm2，期間減少127.29 hm2，平均以25.46 hm2/a的速度遞減，年動(dòng)態(tài)度為4.05 %。中度石漠化2015年較2010年減少129.66 hm2，與強(qiáng)度石漠化變化量相當(dāng)，其年動(dòng)態(tài)度為3.93 %?？梢?jiàn)，示范區(qū)石漠化變化的主要類(lèi)型為無(wú)石漠化、強(qiáng)度石漠化和中度石漠化。

表4 2010-2015年花江示范區(qū)喀斯特石漠化等級(jí)的面積年動(dòng)態(tài)變化Table 4 Area dynamic change of karst rocky desertification grades of Huajiang demonstration area during 2010-2015

注：-表示無(wú)數(shù)據(jù)(下同)；變化量的正數(shù)表示增加，負(fù)數(shù)表示減少。
Note:- indicates no data (the same below); a positive number of changes indicates an increase, and a negative number indicates a decrease.

表5 2010-2015年花江示范區(qū)石漠化等級(jí)的轉(zhuǎn)移矩陣Table 5 Transition matrix of rocky desertification of Huajiang demonstration area during 2010-2015 (hm2)

注：表頭為2015年石漠化等級(jí)，表列為2010年石漠化等級(jí)。
Note: The first row and column are rocky desertification grades in 2015 and 2010 respectively.

表6 2010-2015年花江示范區(qū)土壤的有機(jī)質(zhì)及全氮含量的變化Table 6 Content variance of soil organic matter and total nitrogen of Huajiang demonstration area during 2010-2015

注：標(biāo)準(zhǔn)偏差無(wú)單位，變異系數(shù)單位為 %。
Note: No unit for standard deviation, the unit of variation coefficient is %.

從表5可知，2010-2015年花江示范區(qū)轉(zhuǎn)換面積最大的類(lèi)型為輕度石漠化，其次是潛在石漠化和無(wú)石漠化，石漠化轉(zhuǎn)換主要集中在程度較低的等級(jí)，中度、強(qiáng)度石漠化轉(zhuǎn)換面積相對(duì)小，示范區(qū)石漠化治理的重點(diǎn)為石漠化較輕的等級(jí)。強(qiáng)度和中度石漠化主要向中度、輕度、潛在和無(wú)石漠化等級(jí)趨勢(shì)變化，而部分輕度石漠化區(qū)域變成中度石漠化，潛在石漠化變成輕度和中度石漠化，但變化面積較小，呈石漠化等級(jí)趨于兩級(jí)變化的態(tài)勢(shì)?？傮w而言，石漠化治理效果得明顯，其惡化趨勢(shì)得到有效控制。

2.2 土壤有機(jī)質(zhì)和全氮含量的變化

2.2.1 平均含量 從表6看出，石漠化治理以來(lái)，花江示范區(qū)土壤有機(jī)質(zhì)平均含量由2010年的31.28 g/kg升至2015年的41.83 g/kg，屬于上等水平的豐富向極豐富變化，期間變異系數(shù)由23.79 %升至27.47 %。土壤全氮平均含量由2010年的2.22 g/kg降至2015年的1.52 g/kg，其含量由上等水平的極豐富向豐富水平變化，期間變異系數(shù)由21.17 %升至26.32 %。變異系數(shù)(c.v.%)的大小表示土壤特性空間變異性的大小，通常變異系數(shù)c.v.%≤10 %時(shí)為弱變異性，10 %

2.2.2 不同石漠化等級(jí)含量 從圖2看出，2010-2015年潛在石漠化至強(qiáng)度石漠化土壤有機(jī)質(zhì)含量均呈先升后降趨勢(shì)，輕度和中度石漠化土壤有機(jī)質(zhì)含量明顯高于潛在和強(qiáng)度石漠化。有機(jī)質(zhì)含量：2010年不同石漠化等級(jí)土壤有機(jī)質(zhì)含量變幅較小，其中，輕度和中度石漠化土壤有機(jī)質(zhì)含量較高，分別為33.20和32.79 g/kg；強(qiáng)度石漠化含量最低，為31.22 g/kg。2015年各石漠化等級(jí)土壤有機(jī)質(zhì)含量變幅較大，其中，輕度石漠化土壤有機(jī)質(zhì)含量最高，為51.19 g/kg；中度石漠化其次，為47.85 g/kg；強(qiáng)度石漠化最小，為31.54 g/kg。全氮含量：2010年不同石漠化等級(jí)土壤全氮含量接近，變幅很小，為2.18～2.25 g/kg。2015年土壤全氮含量變幅較大，其中，中度石漠化土壤全氮含量最高，為1.69 g/kg；潛在石漠化其次，為1.63 g/kg；輕度石漠化最低，為1.40 g/kg?？傮w而言，石漠化治理前(2010年)輕度和中度石漠化有機(jī)質(zhì)和全氮顯著高于強(qiáng)度石漠化，石漠化治理后(2015年)潛在、輕度和中度石漠化有機(jī)質(zhì)含量顯著高于強(qiáng)度石漠化，不同等級(jí)石漠化間差異不顯著。

圖2 2010-2015年花江示范區(qū)不同石漠化等級(jí)土壤養(yǎng)分含量的變化Fig.2 Soil nutrients variation in different rocky desertification grades in Huajiang demonstration area during 2010-2015

圖3 2010-2015年不同石漠化等級(jí)土壤養(yǎng)分的變異系數(shù)Fig.3 Variable coefficients of soil nutrients in different rocky desertification grades during 2010-2015

從圖3看出，2010年潛在石漠化至強(qiáng)度石漠化土壤有機(jī)質(zhì)和土壤全氮變異系數(shù)均呈先降后升趨勢(shì)，2015年潛在石漠化至強(qiáng)度石漠化土壤有機(jī)質(zhì)和土壤全氮變異系數(shù)呈先升后降趨勢(shì)，均屬中等變異，且變異性逐年增大。有機(jī)質(zhì)變異系數(shù)：2010年強(qiáng)度石漠化土壤有機(jī)質(zhì)變異系數(shù)最大，為24.15 %；潛在石漠化其次，為23.28 %；中度石漠化最小，為20.98 %。2015年中度石漠化土壤有機(jī)質(zhì)變異系數(shù)最大，為26.83 %；輕度石漠化其次，為24.01 %；強(qiáng)度石漠化最小，為19.21 %。全氮變異系數(shù)：2010年潛在石漠化土壤有機(jī)質(zhì)變異系數(shù)最大，為21.08 %；輕度石漠化其次，為20.89 %；中度石漠化最小，為18.47 %。2015年輕度石漠化土壤有機(jī)質(zhì)變異系數(shù)最大，為37.14 %；潛在石漠化其次，為25.15 %；強(qiáng)度石漠化最小，為17.36 %。

3 討 論

隨著花江示范區(qū)石漠化治理及生態(tài)修復(fù)工程的實(shí)施，特別是區(qū)域內(nèi)以頂壇小流域?yàn)楹诵牡氖鷳B(tài)工程的布局，石漠化趨勢(shì)得到有效控制，中度和強(qiáng)度石漠化嚴(yán)重等級(jí)面積逐漸減少，無(wú)石漠化土地面積逐漸增加，土壤有機(jī)質(zhì)和全氮含量隨著石漠化等級(jí)的變化也產(chǎn)生相應(yīng)變化。

土壤有機(jī)質(zhì)對(duì)生態(tài)修復(fù)具有明顯的指示作用，對(duì)于生態(tài)修復(fù)的響應(yīng)較快。因此，有機(jī)質(zhì)可作為石漠化治理過(guò)程中植被恢復(fù)的物種選擇、恢復(fù)階段的治理效果的評(píng)價(jià)指標(biāo)[11]。但有機(jī)質(zhì)含量受土壤質(zhì)地、植被類(lèi)型、土壤歸還量和人為種植制度等多因素綜合影響，彭佳佳等[12]研究表明，在沙化草地上種植植被后土壤肥力提高，土壤有機(jī)質(zhì)含量顯著增加，并且隨著植被恢復(fù)時(shí)間的延長(zhǎng)而逐漸提高。花江示范區(qū)廣泛種植花椒和金銀花等經(jīng)濟(jì)林，構(gòu)樹(shù)、鐵線(xiàn)蓮和白三葉等灌草叢，以及火龍果等攀援肉質(zhì)灌木，一定程度上起到水土保持作用，從而減輕了土壤養(yǎng)分的流失，同時(shí)枯落物的蓄積和腐蝕分解對(duì)增加有機(jī)質(zhì)含量具有很大促進(jìn)作用。潛在石漠化土壤受除草、施肥及種植作物的影響導(dǎo)致土壤有機(jī)質(zhì)含量低于輕度、中度石漠化退耕還林還草區(qū)。

土壤全氮含量變化趨勢(shì)與有機(jī)質(zhì)變化一致，在無(wú)外源氮肥施入的情況下，土壤氮素主要來(lái)源于有機(jī)質(zhì)的分解、生物固氮和大氣氮沉降。石漠化治理恢復(fù)過(guò)程中，花江示范區(qū)地表植被逐漸恢復(fù)，生物多樣性明顯增加，故土壤全氮含量也趨于增加。然而，2015年土壤全氮含量較2010年略有降低。可能原因：石漠化治理過(guò)程中引種如花椒和金銀花等均為喜氮植物所致。因此，在生態(tài)恢復(fù)治理和農(nóng)事活動(dòng)中應(yīng)增加氮肥的施用量。

4 結(jié) 論

經(jīng)過(guò)石漠化治理，無(wú)石漠化面積從2010年859.57 hm2增至2015年的1066.30 hm2，年動(dòng)態(tài)度為4.81 %；潛在石漠化面積增加42.40 hm2，年動(dòng)態(tài)度為0.94 %；輕度石漠化面積增加7.83 hm2；強(qiáng)度石漠化面積由2010年的628.09 hm2減至2015年的500.81 hm2，平均以25.46 hm2/a的速度遞減，年動(dòng)態(tài)度為4.05 %；中度石漠化減少129.66 hm2，年動(dòng)態(tài)度為3.93 %。強(qiáng)度、中度石漠化主要向中度、輕度、潛在和無(wú)石漠化等級(jí)趨勢(shì)變化?？傮w而言，石漠化治理效果得明顯，其惡化趨勢(shì)得到有效控制?；ń痉秴^(qū)土壤有機(jī)質(zhì)平均含量由2010年的31.28 g/kg升至2015年的41.83 g/kg，土壤全氮平均含量由2010年的2.22 g/kg降至2015年的1.52 g/kg，土壤有機(jī)質(zhì)和全氮含量均為中等變異，且變異性逐年增大。

從土壤有機(jī)質(zhì)和全氮含量看，2010-2015年潛在石漠化至強(qiáng)度石漠化土壤有機(jī)質(zhì)含量均呈先升后降趨勢(shì)，輕度和中度石漠化土壤有機(jī)質(zhì)含量明顯高于潛在和強(qiáng)度石漠化。2010年不同石漠化等級(jí)土壤有機(jī)質(zhì)含量變幅較小，為31.22～33.20 g/kg。2015年各石漠化等級(jí)土壤有機(jī)質(zhì)含量變幅較大，為31.54～51.19 g/kg。2010年不同石漠化等級(jí)土壤全氮含量接近，變幅很小，為2.18～2.25 g/kg。2015年土壤全氮含量變幅較大，為1.40～1.69 g/kg，與吳育忠[11]的研究結(jié)果相近。總之，石漠化治理前(2010年)輕度和中度石漠化有機(jī)質(zhì)和全氮顯著高于強(qiáng)度石漠化，石漠化治理后(2015年)潛在、輕度和中度石漠化有機(jī)質(zhì)含量顯著高于強(qiáng)度石漠化，不同等級(jí)石漠化間差異不顯著。

從變異系數(shù)看，2010年潛在石漠化至強(qiáng)度石漠化土壤有機(jī)質(zhì)和土壤全氮變異系數(shù)均呈先降后升趨勢(shì)，2015年潛在石漠化至強(qiáng)度石漠化土壤有機(jī)質(zhì)和土壤全氮變異系數(shù)呈先升后降趨勢(shì)，均屬中等變異，且變異性逐年增大。有機(jī)質(zhì)變異系數(shù)和全氮變異系數(shù)2015年的變幅均大于2010年。