孫 悅,蔡湘文,程 洋,蘇春田

(1.桂林理工大學(xué) 地球科學(xué)學(xué)院,廣西 桂林 541006;2.中國地質(zhì)科學(xué)院 巖溶地質(zhì)研究所,廣西 桂林 541004)

基于RS和GIS的巖溶石漠化與大氣降雨關(guān)系
——以貴州省典型巖溶石漠化地區(qū)為例

孫 悅1,蔡湘文1,程 洋2,蘇春田2

(1.桂林理工大學(xué) 地球科學(xué)學(xué)院,廣西 桂林 541006;2.中國地質(zhì)科學(xué)院 巖溶地質(zhì)研究所,廣西 桂林 541004)

以貴州省典型巖溶石漠化地區(qū)為例,采用不同時相的TM多光譜數(shù)據(jù)和19個基準(zhǔn)地面氣象觀測站及自動站1999—2013年的逐日降雨數(shù)據(jù)為數(shù)據(jù)源,以遙感(RS)、地理信息系統(tǒng)(GIS)為技術(shù)手段和平臺，1999—2013年的年均降雨量和年暴雨頻率為評價因子,運用地統(tǒng)計學(xué)和空間信息分析方法分析巖溶石漠化與大氣降雨的關(guān)系。結(jié)果表明:降雨(非暴雨)是加劇巖溶石漠化進(jìn)程的因素之一,也是減緩巖溶石漠化進(jìn)程的重要因素;暴雨加速了巖溶石漠化的演化進(jìn)程,對巖溶石漠化的恢復(fù)起反作用。

遙感;地理信息系統(tǒng);巖溶石漠化;降雨;暴雨

0 引 言

巖溶石漠化是指在南方熱帶、亞熱帶脆弱的巖溶環(huán)境中,在不合理的人類活動基礎(chǔ)上造成的植被退化、土壤侵蝕和大面積的基巖裸露，土地生產(chǎn)力嚴(yán)重下降,地表出現(xiàn)類似荒漠景觀的土地退化過程[1-4]。

當(dāng)前,巖溶石漠化已成為中國西南巖溶區(qū)最嚴(yán)重的地質(zhì)環(huán)境問題, 嚴(yán)重制約了當(dāng)?shù)厣鐣?jīng)濟(jì)的發(fā)展, 威脅當(dāng)?shù)鼐用竦纳鏃l件, 對長江、 珠江下游地區(qū)的防洪和生態(tài)環(huán)境是巨大的潛在威脅。 為此， 黨的十八大提出“要實施重大生態(tài)修復(fù)工程,增強(qiáng)生態(tài)產(chǎn)品生產(chǎn)能力,推進(jìn)石漠化綜合治理”。

巖溶石漠化的成因及其演變機(jī)理的研究是進(jìn)行有效的石漠化綜合治理的前提。近年來,許多研究者對巖溶石漠化成因及其演變機(jī)理的研究取得了不錯的成績[5-13],證明了地層巖性、地形坡度、水文地貌、土地利用類型和人口密度是巖溶石漠化發(fā)生和發(fā)展的主要驅(qū)動因素,為石漠化治理提供了理論依據(jù)。本次研究的目標(biāo)是定量分析巖溶石漠化與大氣降雨的關(guān)系,推動不同地形條件下石漠化演變規(guī)律研究,完善巖溶石漠化成因及其演變機(jī)理的理論體系,為制定有前瞻性和針對性的生態(tài)環(huán)境保護(hù)和恢復(fù)規(guī)劃提供科學(xué)依據(jù),推進(jìn)石漠化防治工作的進(jìn)程。

地處中國西南的貴州省以連續(xù)分布面積最大、發(fā)育類型最齊全、景色最秀麗的巖溶景觀著稱于世。但是,由于巖溶生態(tài)地質(zhì)環(huán)境的脆弱性和人類不合理社會經(jīng)濟(jì)活動,貴州省是我國巖溶石漠化面積最大、發(fā)育程度最嚴(yán)重的省份。以貴州省為研究區(qū)定量分析巖溶石漠化與大氣降雨的關(guān)系有重要的現(xiàn)實意義。

1 研究思路

基于RS技術(shù)分別提取1999年和2013年的巖溶石漠化信息,以兩期巖溶石漠化信息為基礎(chǔ),分析貴州省的石漠化演變情況;利用地統(tǒng)計學(xué)研究貴州省的年均降雨量和年均暴雨頻率;最后,基于GIS平臺和空間信息分析法研究大氣降雨與巖溶石漠化的關(guān)系。以“獲取數(shù)據(jù)→處理數(shù)據(jù)→分析數(shù)據(jù)→結(jié)論”為研究路線,詳細(xì)的技術(shù)路線如圖1所示。

2 巖溶石漠化信息提取

2.1 基于RS技術(shù)的巖溶石漠化信息自動提取

本次研究的遙感數(shù)據(jù)為1999年9月和2013年8月獲取的2個時相的TM數(shù)據(jù),2景遙感圖像的云量均小于5%,無明顯噪聲,影像質(zhì)量很好，滿足巖溶石漠化信息自動提取的需求(圖2、圖3)。

根據(jù)地物波譜特性, 植被在TM4(0.76～0.90 μm)波段的高反射率和絕大多數(shù)礦物在TM5(1.55～1.75 μm)波段的高反射率(對裸露基巖反映較好),故首先選取TM5/TM4運算為石漠化信息增強(qiáng)方式,采用二次圖像信息分析分類方法[2，12]提取石漠化信息;然后通過人機(jī)交互解譯提高石漠化信息的提取精度,減少“同物異譜”和“同譜異物”現(xiàn)象給計算機(jī)自動識別和信息提取帶來的誤差。這種巖溶石漠化信息遙感提取方法有較好的實踐結(jié)果,其成果精度優(yōu)于1∶25萬成果精度要求[14],符合本次研究的精度要求。

2.2 巖溶石漠化演變分析

根據(jù)石漠化的演變程度不同, 將石漠化演變分為5個等級:較好改善、 輕度改善、 基本未變、 輕度惡化、 嚴(yán)重惡化。 將1999和2013年的巖溶石漠化信息進(jìn)行空間疊加分析, 研究區(qū)巖溶石漠化演變?nèi)鐖D4所示， 各石漠化演變等級面積統(tǒng)計見表1。

圖1 技術(shù)路線圖

圖2 1999年巖溶石漠化圖

圖3 2013年巖溶石漠化圖

圖4 1999—2013年巖溶石漠化演變圖

表1 1999—2013年巖溶石漠化演變等級面積統(tǒng)計

Table 1 Statistics of karst rocky desertification evolution areas from 1999 to 2013

碳酸鹽巖輕度惡化重度惡化基本未變輕度改善較好改善面積/km21212161401304104192131602159比例/%-1.160.2585.9510.861.78

注：碳酸鹽巖面積是根據(jù)1∶25萬區(qū)域地質(zhì)圖統(tǒng)計。

可知,研究區(qū)石漠化以改善為主,總改善面積為15 319 km2,是總惡化面積的近9倍,占碳酸鹽巖總面積的12.64%;總惡化面積為1 705 km2,占碳酸鹽巖總面積的1.41%。

3 1999—2013年貴州省大氣降雨分析

3.1 1999—2013年貴州省大氣降雨概況

氣象數(shù)據(jù)來自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)中的《中國地面氣候資料日值數(shù)據(jù)集》(http://www.cma.gov.cn/2011qxfw/2011qsjgx/),選取了貴州省19個基本、基準(zhǔn)地面氣象觀測站及自動站1999—2013年的逐日降雨數(shù)據(jù),統(tǒng)計了1999—2013年貴州省的年均降雨量和年均暴雨(日降雨量大于50 mm)天數(shù)，詳見表2。

3.2 基于地統(tǒng)計學(xué)和空間差值研究貴州大氣降雨的空間分布規(guī)律

采用協(xié)同克里格法對年均降雨量和年均暴雨天數(shù)進(jìn)行空間差值,貴州省1999—2013年年均降雨量空間分布和年均暴雨天數(shù)空間分布如圖5、圖6所示。

表2 1999—2013年貴州省大氣降雨統(tǒng)計

圖5 1999—2013年年均降雨量空間分布圖

4 巖溶石漠化與大氣降雨關(guān)系分析

采用基于GIS平臺的空間信息分析方法分析巖溶石漠化與大氣降雨的相關(guān)關(guān)系。將1999—2013年的貴州省巖溶石漠化演變信息與貴州省降雨量空間分布圖和年均暴雨天數(shù)空間分布圖分別進(jìn)行空間疊加分析,相關(guān)數(shù)據(jù)見表3和圖7。

圖6 1999—2013年年均暴雨天數(shù)空間分布圖

表3 巖溶石漠化演變與年均降雨量相關(guān)性統(tǒng)計

Table 3 Statistics of evolution and annual mean rainfall in karst rocky desertification area %

年均降雨量/mm最小最大8009009001000100011001100120012001300輕度惡化率0.210.320.732.252.85重度惡化率0.030.060.140.590.61基本未變率1.964.047.7026.9727.92輕度改善率6.894.375.3916.5820.74較好改善率1.070.650.802.314.92

注: 惡化率=年均降雨量惡化面積/年均降雨量面積。

圖7 巖溶石漠化演變與年均降雨量

可見,巖溶石漠化惡化與降雨量呈明顯的正相關(guān),石漠化惡化隨降雨量的增加變嚴(yán)重,尤其是降雨量從1 000 mm遞增到1 200 mm時,石漠化惡化率增加了近3倍。大量的降雨導(dǎo)致水土流失,基巖裸露,最終加劇巖溶石漠化進(jìn)程。石漠化改善與降雨量沒有明顯的相關(guān)關(guān)系,大氣降雨在造成水土流失的同時也會促進(jìn)植被的生長,非暴雨降雨是石漠化恢復(fù)與重建的重要因素之一。

由表4和圖8可知,巖溶石漠化惡化與年均暴雨天數(shù)呈明顯的正相關(guān),石漠化惡化隨年均暴雨日數(shù)的增加變嚴(yán)重。暴雨對地表造成強(qiáng)烈的沖刷,加劇水土流失、基巖裸露,大大縮短巖溶石漠化演化進(jìn)程。石漠化改善與降雨量沒有明顯的相關(guān)關(guān)系,年均暴雨天數(shù)超過5.5天的地區(qū)石漠化改善率明顯低于年均暴雨天數(shù)小于5.5天的地區(qū),暴雨對石漠化恢復(fù)與重建起反作用。

表4 研究區(qū)巖溶石漠化演變信息與年均暴雨天數(shù)統(tǒng)計

注: 惡化率=年均暴雨次數(shù)惡化面積/年均暴雨天數(shù)面積。

圖8 巖溶石漠化演變與年均暴雨天數(shù)關(guān)系

5 結(jié) 論

RS和GIS是研究大空間尺度問題最有效的技術(shù)方法之一,利用RS技術(shù)快速、準(zhǔn)確地提取了研究區(qū)超過12萬km2的巖溶石漠化信息。GIS強(qiáng)大的數(shù)據(jù)分析功能為石漠化演變規(guī)律、年均降雨量和年均暴雨日數(shù)的空間分布的分析提供了基礎(chǔ)。大氣降雨是巖溶石漠化演變的重要因素,降雨(非暴雨)一方面加速了巖溶石漠化的惡化,石漠化惡化隨降雨量的增加變嚴(yán)重;另一方面,降雨(非暴雨)會促進(jìn)植被的生長,是石漠化改善的重要因素之一。暴雨是石漠化惡化的重要因素,巖溶石漠化惡化與年均暴雨天數(shù)呈明顯的正相關(guān),石漠化惡化隨年均暴雨日數(shù)的增加變嚴(yán)重,暴雨對石漠化改善起反作用。

通過上述研究,在石漠化治理工程中應(yīng)考慮緩解暴雨沖刷的工程措施,減輕暴雨對巖溶石漠化的不利影響。大尺度巖溶石漠化信息遙感自動提取的方法還可以改進(jìn),可以引入面向?qū)ο蟮倪b感分類方法和基于時間序列的遙感分類方法。氣象站點的數(shù)量還應(yīng)該進(jìn)一步增加,氣象數(shù)據(jù)的空間差值算法還有較大的提升空間,這兩者的進(jìn)步能顯著提升年均降雨量空間分布圖和年均暴雨天數(shù)空間分布圖的精度。

[1]袁道先.我國西南巖溶石山的環(huán)境地質(zhì)問題[J].世界科技研究與發(fā)展,1997,19(5):41-43.

[2]王世杰.喀斯特石漠化概念演繹及其科學(xué)內(nèi)涵的探討[J].中國巖溶,2002,21(2): 101-105.

[3]王德爐,朱守謙,黃寶龍.石漠化的概念及其內(nèi)涵[J].南京林業(yè)大學(xué)學(xué)報:自然科學(xué)版, 2004, 28(6): 87-90.

[4]謝雨萍，吳虹，劉澤東，等.恭城縣巖溶石漠化環(huán)境變化定量遙感研究[J].桂林工學(xué)院學(xué)報，2009, 29(1): 65-71.

[5]程洋,陳建平,皇甫江云,等.基于RS和GIS的巖溶石漠化惡化趨勢定量預(yù)測——以廣西都安瑤族自治縣典型巖溶石漠化地區(qū)為例[J].國土資源遙感, 2012, 24(3):135-139.

[6]張信寶,王世杰,白曉永,等.貴州石漠化空間分布與喀斯特地貌、巖性、降水和人口密度的關(guān)系[J].地球與環(huán)境, 2013, 41(1): 1-6.

[7]李陽兵,白曉永,周國富,等.中國典型石漠化地區(qū)土地利用與石漠化的關(guān)系[J].地理學(xué)報, 2006, 61(6): 624-632.

[8]胡寶清,廖赤眉,嚴(yán)志強(qiáng),等.基于RS和GIS的喀斯特石漠化驅(qū)動機(jī)制分析——以廣西都安瑤族自治縣為例[J].山地學(xué)報,2004, 22(5): 583-590.

[9]李陽兵,王世杰,周夢維,等.不同空間尺度下喀斯特石漠化與坡度的關(guān)系[J].水土保持研究, 2009,16(5): 70-72，77.

[10]王玉娟,楊勝天,溫志群，等.貴州典型喀斯特灌叢草坡類型區(qū)土壤水分及其影響因子研究[J].北京師范大學(xué)學(xué)報:自然科學(xué)版, 2008,44(5):529-532.

[11]張躍紅,安裕倫,馬良瑞,等. 1960—2010 年貴州省喀斯特山區(qū)陡坡土地利用變化[J].地理科學(xué)進(jìn)展, 2012,31(7): 878-884.

[12]童立強(qiáng).西南巖溶石山地區(qū)石漠化信息自動提取技術(shù)研究[J].國土資源遙感,2003, 15(4):35-38.

[13]孟小軍，莫源富.基干遙感技術(shù)的打狗河流域1999—2009年間石漠化及綠化研究[J].桂林理工大學(xué)學(xué)報, 2013, 33(4): 622-628.

[14]童立強(qiáng),劉春玲,聶洪峰.中國南方巖溶石山地區(qū)石漠化遙感調(diào)查與演變研究[M].北京:科學(xué)出版社,2013.

Relationship between karst rocky desertification and atmospheric precipitation based on RS and GIS —A typical karst rocky desertification area case from Guizhou

SUN Yue1, CAI Xiang-wen1,CHENG Yang2, SU Chun-tian2

(1.College of Earth Sciences,Guilin University of Technology, Guilin 541006,China;2.Institute of Karst Geology, Chinese Academy of Geological Science,Guilin 541004,China)

Taking the typical karst rocky desertification area in Guizhou Province as an example, TM multispectral data in different time phases and daily rainfall data from 1999 to 2013 in 19 basic ground meteorological observation stations and automatic stations as data source, remote sensing (RS) and geographical information system (GIS) as technical means and platform, the average annual rainfall and rainstorm frequency from 1999 to 2013 as evaluation factors,the relationship between karst rocky desertification and atmospheric precipitation is analyzed by the method of geo-statistics and spatial information analysis.The results show that rainfall (non storm) is one of the factors that aggravate the karst rocky desertification process. It is also an important factor in slowing the karst rocky desertification process. Rainstorm accelerated the evolution process of karst rocky desertification, and it was counterproductive to the recovery of karst rocky desertification.

RS;GIS;karst rocky desertification;rainfall;rainstorm

1674-9057(2015)04-0834-05

10.3969/j.issn.1674-9057.2015.04.023

2015-05-04

中國地質(zhì)調(diào)查局地質(zhì)調(diào)查項目(12120114069001-01；12120115046601);桂林理工大學(xué)自然地理學(xué)重點學(xué)科項目

孫 悅(1990—),女,碩士研究生,研究方向:環(huán)境遙感。

蔡湘文,博士,副教授, 2004042@glut.edu.cn。

孫悅,蔡湘文,程洋,等.基于RS和GIS的巖溶石漠化與大氣降雨關(guān)系——以貴州省典型巖溶石漠化地區(qū)為例[J].桂林理工大學(xué)學(xué)報，2015,35(4):834-838.

TP79；P642.254；P426.615

A