桂西北喀斯特区域植被变化趋势及其对气候和地形的响应

童晓伟，王克林，岳跃民，廖楚杰，徐艳芳，朱海涛

(1.中国科学院亚热带农业生态研究所亚热带农业生态过程重点实验室， 长沙 410125；2.中国科学院环江喀斯特生态系统观测研究站， 环江 547100；3.中国科学院大学， 北京 100049；4.环境保护部卫星环境应用中心， 北京 100094)

桂西北喀斯特区域植被变化趋势及其对气候和地形的响应

童晓伟1,2，3，王克林1，2，*，岳跃民1，2，廖楚杰1,2，3，徐艳芳1,2，3，朱海涛4

(1.中国科学院亚热带农业生态研究所亚热带农业生态过程重点实验室， 长沙 410125；2.中国科学院环江喀斯特生态系统观测研究站， 环江 547100；3.中国科学院大学， 北京 100049；4.环境保护部卫星环境应用中心， 北京 100094)

基于1999—2010年的SPOT NDVI数据，分析了河池市植被变化趋势及空间差异，并结合气象和地形数据分析了植被与气候、地形的关系。结果表明：(1)桂西北喀斯特地区植被变化总体上呈恢复趋势，年均气候因子对植被变化的作用不明显；(2)200—500m的海拔范围内植被恢复显著，但400—500m的海拔范围内有小面积植被退化现象，随着海拔增加，植被变化趋于稳定；(3)6—15°的坡度范围内植被恢复最显著，而2—6°及大于25°坡度范围存在植被退化的现象；(4)不同坡向上的植被恢复差异不明显，但随着坡向由阴坡转阳坡，植被总体恢复呈减小趋势。喀斯特地区人类生态建设取得一定成效，但由于人类活动的负面影响，在海拔400—500m、坡度大于25°的阳坡区域仍存在植被减少的现象。

SPOT-VGT NDVI；气候变化；地形面积差异修正系数；生态建设工程; 喀斯特

脆弱的地质背景加上长期以来人类的不合理干扰，使我国西南喀斯特地区自然植被不断遭到破坏，土壤侵蚀严重，基岩大面积裸露，土地生产力下降，呈现出类似荒漠的石漠化景观[1- 2]。植被保护与恢复对遏制石漠化发生、有效控制水土流失、改善喀斯特生态环境起着至关重要的作用。国家先后在西南喀斯特地区实施了退耕还林还草、保护、天然草地植被恢复与建设等一系列生态建设工程，开展生态恢复工作，尤其在2008年国务院批复了《岩溶地区石漠化综合治理大纲(2006—2015)》以来，喀斯特地区植被发生了较大变化。因此，植被动态监测与分析成为喀斯特地区生态建设工程成效监测与评估的关键内容之一。

我国西南喀斯特地区作为生态系统退化比较严重的典型生态脆弱区，地貌类型多样，生境具有高度的景观异质性，植被变化不仅受气候变化影响和人类活动的干扰，也受喀斯特特殊地形地貌的影响。目前，喀斯特地区植被长时间序列变化的研究较少，且多集中在大尺度植被与气候关系的探讨上，较少考虑区域尺度人为因素及喀斯特特殊地形地貌对喀斯特植被变化的影响[3- 5]。复杂的喀斯特山区地形，影响太阳辐射和降水的空间再分配，不同喀斯特地区土壤质地、水分和养分的差异造成植被生境的高度异质性，对植被变化产生重要影响。然而，在研究地形因子对植被变化的影响时，往往忽略了不同地形因子绝对面积差异的影响，如某一地形因子条件下，植被发生变化的面积可能较小，但该地形因子条件下植被发生变化的比例可能占所有不同地形条件下植被变化的比例较大，从而造成评价地形因子对植被变化影响的不确定性[6]。因此，本研究综合考虑气候因子、喀斯特地形因子对喀斯特区域植被变化的影响，提出地形面积差异修正方法，揭示不同地形因子条件下喀斯特植被恢复趋势差异，明确植被恢复成效显著的区域及特别需要加强生态治理的区域，为喀斯特地区石漠化综合治理、退耕还林等生态建设工程的生态效益评估提供科学依据。

1 研究地区与研究方法

1.1 研究区概况

河池市(106°34′—109°09′E，23°41′—25°37′N)地处广西西北边陲、云贵高原南麓，喀斯特地貌广泛发育，占地区总面积的66%，地形破碎、地貌类型多样，是举世闻名的中国西南喀斯特地区之一。现辖宜州、环江、都安、大化、凤山等11个县(区、市)(图1)，全市土地面积33500km2，总人口409万人，是广西壮族自治区少数民族聚居最多的地区之一。该地区属于亚热带季风气候区，热量丰富，雨量充沛，且雨热同期；无霜期长，历年平均无霜期约为329d。河池市石漠化严重，广西有12个县属于国家石漠化治理首批试点县，而河池就占有4个，2012年公布的石漠化监测结果显示，全市石漠化土地面积8520km2，潜在石漠化土地面积11470 km2，占全市土地总面积的近60%。

图1 河池市行政区划图Fig.1 The administrative map of Hechi

1.2 数据来源及处理

1.2.1 SPOT-NDVI和气象数据

归一化植被指数(NDVI)是反映植被覆盖状况的良好指标[7- 8]，是公认的用来表征植被变化的有效遥感参数[9- 10]。本文遥感数据选用的是 SPOT VEGETATION NDVI 10 天合成的产品(http://free.vgt.vito.be/)，每年36期数据，空间分辨率为1km，时间范围是1999年1月至2010年12月，共计432期影像。获得的SPOT-VGT遥感数据已经过MVC(maximum value composite)方法处理，为了进一步减少云、大气等噪声对数据的影响，利用Chen等[11]提出的基于Savitzky-Golay滤波器的平滑方法，对NDVI时间序列数据进行平滑滤波处理，在此基础上再次利用最大合成法获得月NDVI值，最后取12个月NDVI的平均值得到年均NDVI，用以表征植被年均生长状况。气候数据从中国气象局获得(http://cdc.cma.gov.cn/home.do)，包括河池地区及周边14个气象站点。数据内容为1999年1月—2010年12月逐月的月平均气温，降水量以及各站点的经度和纬度。在Arcgis 9.3软件的空间分析模块中对14个站点的气象数据进行Kriging插值，利用各项交叉验证指标，选取合适的插值模型，使插值结果精度满足数据分析要求。

1.2.2 土地覆盖和DEM数据

土地覆盖数据从寒区旱区科学数据中心获取(http://westdc.westgis.ac.cn/)，该数据是在中国科学院2000年1∶10万按县分幅的土地资源调查成果的基础上进行合并、矢栅转换(面积最大法)后得到全国幅的土地利用数据产品，采用中国科学院资源环境分类系统，是目前国家尺度上精度最高的土地利用数据产品，已被国内外学者广泛应用于相关研究中，其空间分辨率为1km。DEM数据从国际科学数据服务平台获得(http://datamirror.csdb.cn/)，空间分辨率90m，利用Arcgis9.3软件提取研究区海拔、坡度和坡向。本文采用100m等间距划分法将研究区海拔划分为12个等级，其中海拔1100m以上区域所占面积比例较小，故将其合并为一个等级；将坡向划分为东、西、南、北4个坡向，分别表示半阴坡，半阳坡，阳坡，阴坡[12]；根据中国农业区划委员会颁发的《土地利用现状调查技术规程》，将坡度划分为5级，即≤2°、2°—6°、6°—15°、15°—25°和>25°，其中<6°为平耕地，6—25°为缓坡， >25°为陡坡耕地。将以上3种地形因子分级栅格数据进行矢量化处理，生成分级矢量数据，统一坐标投影为WGS 1984 UTM Zone49 N。

1.3 研究方法

1.3.1 趋势分析

趋势线是对一组随时间变化的序列数据进行回归分析，预测其长期趋势的变化态势。趋势线包括直线和曲线两种，采用线性趋势线来模拟植被覆盖的年际变化，即最小二乘拟合直线,其斜率计算公式如：

(1)

式中,变量i为1—12的年序号；NDVIi表示第i年的年均NDVI值。变化趋势图反映了在12a的时间序列中，河池市植被NDVI的变化趋势。某像元的趋势线是这个像元12a的年均NDVI值用一元线性回归模拟出来的一个总的变化趋势，θslope即这条趋势线的斜率，θslope> 0说明植被在12a间的变化趋势是增加的，反之则是减少，θslope=0表示植被无变化趋势。对斜率的显著性进行F检验，并依据F值将植被变化趋势划分为不显著 (F<4.96)和显著 (F>4.96)两种。利用植被变化趋势和显著性检验的叠加分析，将研究区植被变化趋势分为减少，不变和增加3种类型。

对研究区年均降水和年均气温的时间序列变化趋势采用Mann-Kendall分析法， Mann-Kendall检验法是世界气象组织推荐的非参数统计检验法，广泛应用于分析降水、气温和径流等要素时间序列的趋势变化[13]。其检验统计量公式为：

(2)

式中， sign()为符号函数。当Xi-Xj小于、等于或大于零时，sign=(Xi-Xj)分别为-1、0或1；M-K统计量公式S大于、等于、小于零时分别为：

(3)

式中，Z为正值表示增加趋势，负值表示减少趋势。Z的绝对值在大于等于1.28、1.64、2.32时表示分别通过了信度90%、95%、99%显著性检验。

1.3.2 偏相关分析

偏相关分析是在对其他变量的影响进行控制的条件下，衡量多个变量中某两个变量之间的线性相关程度，在确定两个变量之间的内在线性联系时会更真实、更可靠。本研究利用Matlab 6.5对植被年均NDVI和两种主要气候因子进行偏相关分析，计算河池市各像元NDVI与气候因子的偏相关系数，分析各量之间的偏相关程度及空间分布特征。

1.3.3 地形面积差异修正

不同地形条件下植被变化特征不同，如某一植被变化趋势类型在特定地形条件下的分布面积可能较小，然而与整个研究区此植被变化趋势类型所占的面积比例相比，此地形条件下该植被变化趋势类型分布却较大，从而造成评价不同地形因子对植被变化影响的不确定性。因此，需要消除不同地形因子绝对面积不同的影响，揭示特定地形因子对植被变化的影响，明确不同地形条件对植被变化类型分布的影响及其演变趋势。本研究提出地形面积差异修正系数k，用来消除各地形因子条件下因地形绝对面积不同而引起的植被恢复评价的不确定性，其计算公式为：

(4)

2 结果与分析

2.1 植被变化总体趋势分析

SPOT NDVI时间序列分析显示，1999—2010年河池市69.9%的区域植被增加，29.9%的区域植被不变，植被减少的区域仅占0.2%(图2)。植被变化趋势和土地覆盖叠加分析表明，植被增加的区域，其土地覆盖类型以有林地及高、中覆盖度草地为主；西北及南部地区植被基本不变，且存在极少数植被减少的区域，其土地覆盖类型以中覆盖度草地及灌木疏林地为主。说明在1km空间尺度、年际尺度上，河池市近12年来植被以增加为主，表明研究区总体生态状况呈恢复改善的趋势，但由于毁林开垦、陡坡耕种、过度放牧等现象的存在，导致极小部分地区存在植被减少的现象。

图2 河池市植被变化趋势空间分布Fig.2 The distribution of vegetation change trends in Hechi

2.2 植被变化与气候因子的关系

近12年来河池市年均降水为1193 mm，年均温为20.4℃，1999—2010年，河池市年均降水在波动变化中呈较小的下降趋势，年均气温呈微弱的上升趋势(图3)。Mann-Kendall分析表明，河池市约63.4%的区域年均降水呈减少趋势，广泛分布在中部及东部，其中约3.1%的区域在置信度为 90%水平上呈显著减少趋势，主要分布在罗城、宜州及都安3县；研究区西部地区降水呈微弱增加趋势，仅3.2%的区域在置信度为 90%水平上呈显著增加趋势，主要分布在巴马和凤山县，其中巴马县有207km2的区域在置信度为 95%水平上存在显著增加趋势。对年均温而言，仅有1.54%的区域在置信度为 90%水平上呈显著增加趋势，主要分布在天峨、河池、大化和都安四县。由此可见，近12年来河池市年均降水呈小幅下降趋势，年均气温呈微弱的上升趋势，但总体变化趋势不显著。

图3 河池地区1999—2010年均降水量、气温变化趋势Fig.3 The change trends of temperature and precipitation in Hechi during 1999—2010

植被变化与气候因子(年均温、年均降水)的偏相关分析可知，河池市植被对降水、气温变化的响应存在明显差异(图4)，其偏相关显著区域可分为：负相关不显著，相关显著，正相关不显著3种。研究区植被NDVI与气温基本呈正相关，比例达95.4%，但相关性不显著；成负相关的区域主要分布在大化、凤山及都安3县，仅5.1%的区域其植被NDVI与气温通过了α=0.1的显著性检验水平。对降水而言，以东北—西南方向连线划分研究区，此线以南植被NDVI与降水呈负相关，以北呈正相关，两者所占比例相差不大，但显著性水平达到0.1的区域仅占4%，且以负相关为主，集中分布在都安县。由此表明，河池市年均气温变化对植被NDVI的影响略高于年均降水变化对其的影响，但整体表现为植被NDVI与气候因子的相关性不显著，说明在区域尺度上，近12年的年均气候变化对研究区植被变化作用不明显。

图4 河池市植被年均NDVI与气温、降水偏相关显著性检验空间分布Fig.4 The significant test of partial correlation between mean annual NDVI and temperature and precipitation in Hechi

2.3 植被变化对不同地形因子的响应

2.3.1 海拔对植被变化的影响

植被变化趋势类型在不同海拔范围内的分布差异较大，从不同植被变化类型随海拔的变化趋势来看(图5)，植被减少类型在高程小于500m的范围内呈趋于增强的分布趋势(1.45

从不同海拔范围内植被变化类型所占的面积百分比来看(图6)，不同海拔范围内，植被变化均以不变和增加为主，植被减少类型所占比例极小。植被增加类型的比例除在海拔小于100m的范围内较小外(约33%)，其余各海拔区段所占比例均大于54%，表明除海拔较低(小于100m)区域外，研究区植被总体在恢复，海拔200—500m范围内的植被恢复显著(植被增加类型所占比例均大于71%)，其中，200—300m海拔范围内植被恢复最显著，其植被增加类型所占比例约78%。然而， 200—500m的海拔区域也是植被覆盖变化易于波动的区域，易于恢复，同时也易于退化(400—500m)，这主要是两方面的原因：一方面耕地主要分布在此海拔范围内，受农作物等植被类型的生长周期变化的影响较大；另一方面退耕还林等生态建设工程的实施促进了植被恢复，而喀斯特地区毁林开垦、樵采薪材等人为逆向干扰活动导致了植被退化。随着海拔的升高，特别是海拔大于800m的区域，植被变化趋于稳定，这主要是因为海拔高的地区坡度往往也更陡峭，受人类活动干扰少，同时，封山育林等政策的实施也有利于高海拔地区的植被结构稳定。

图5 地形面积差异修正系数k随海拔的变化趋势Fig.5 Topographic factor k applied to different vegetation change trends depending on elevation

图6 不同植被变化类型在不同海拔范围内的面积比Fig.6 The area rate of vegetation change types at different elevations

2.3.2 坡度对植被变化的影响

从不同坡度上研究区植被变化趋势来看，坡度对植被增加和不变两种类型分布的影响程度差异较小，对植被减少类型的分布影响较大(图7)。植被减少类型随坡度增大表现为先增加后减少再趋于增加的分布趋势，在2—6°的坡度范围内植被减少类型分布最明显(k=1.87)，其次为>25°的坡度范围，6—25°坡度范围内植被减少类型的变化趋于减少。而植被增加类型随坡度增大差异不明显，仅在6—15°的坡度范围内植被增加类型分布趋于小幅度的增加(k=1.04)。植被不变的类型随坡度增大表现为先减后增的分布趋势，在小于2°及大于15°的坡度范围内植被不变的类型分布明显。从不同坡度上各植被类型所占的面积百分比来看(图8)，不同坡度范围内，植被变化均以增加和不变为主，且植被增加类型所占比例均在68%以上，植被减少类型所占比例极小，表明不同坡度上植被总体在恢复。综合k值及不同坡度上各植被变化趋势类型的面积比例可知，6—15°的坡度范围内植被减少类型分布不明显、植被增加类型分布最明显，且植被增加类型的面积比也最大(约73%)，因此，河池地区在6—15°坡度范围内植被恢复最显著。退耕还林还草政策和石漠化治理等一系列生态建设工程的实施，使喀斯特地区不同坡度上植被总体在恢复。然而，由于陡坡耕种等人类活动的负面影响，导致坡度大于25°的区域仍然存在极少数的植被减少；同时由于耕地主要分布于2—6°等坡度较缓的区域，农作物生长周期的变化使该坡度范围内的植被减少类型分布较明显。

图7 地形面积差异修正系数k随坡度的变化趋势Fig.7 Topographic factor k applied to different vegetation change trends depending on slope

图8 不同植被变化类型在不同坡度上的面积比Fig.8 The area rate of vegetation change types at different slopes

2.3.3 坡向对植被变化的影响

由不同坡向上植被变化类型的分布趋势可知(图9)，植被减少类型随着坡向由阴坡转阳坡(北坡→东坡、西坡→南坡)呈增加的分布趋势，北坡和东坡上此植被变化类型的分布趋于减少、南坡和西坡上此植被变化类型趋于增加；而植被增加类型在不同坡向上的分布差异不明显，仅在北坡趋于较小幅度的增加(k=1.06)、南坡趋于较小幅度的减少(k=0.96)。植被不变的类型随坡向由阴转阳呈先增加后减少的趋势，但变化幅度较弱(0.86

图9 地形面积差异修正系数k随坡向的变化趋势Fig.9 Topographic factor k applied to different vegetation change trends depending on aspect

图10 不同植被变化类型在不同坡向上的面积比Fig.10 The area rate of vegetation change types at different aspects

3 讨论

本研究利用12a的气象和植被指数数据，分析气候变化对喀斯特区域植被变化的影响，结果表明，河池市在近12年里气候变化趋势不显著，但降水已呈较小的下降趋势、气温呈微弱的上升趋势，说明研究区气候因子在一定程度上已呈现出暖干化的趋势；同时，NDVI与年均气温、降水存在一定相关性，但显著相关区域所占比例很小，据此可初步认为在区域尺度上，非气候因素是引起近12年来喀斯特区域植被增加的主要原因，这与其他地区区域尺度上的研究结论类似[15- 19]，说明在区域尺度上，喀斯特区域植被的增加虽不能完全排除气候变化的影响，但人类活动已成为影响植被变化的重要因素，近年来的植树造林、退耕还林还草等一系列生态建设工程所带来的生态效益正在呈现。

为了明确不同地形因子条件下植被恢复的差异，提出了地形面积差异修正系数k来消除由于不同地形因子绝对面积差异而导致的恢复成效评价的不确定性。该修正系数k以某一植被变化趋势类型占特定地形面积的比例作为基准，并归一化到以k=1作为植被演变的平衡变化点，不但有效消除了地形绝对面积差异的影响，而且使判断特定地形因子条件下植被的变化趋势成为可能。研究结果表明，不同喀斯特地形条件下植被恢复差异明显，从植被变化在单一地形因子的表现看，海拔200—500m范围内的植被恢复显著(200—300m恢复最显著)，但此海拔范围区域也是植被覆盖变化易于波动的区域，同时，海拔大于800m的区域，植被变化趋于稳定。参考现有相关研究[20]，并结合桂西北地区植被恢复野外调查，发现：广西人工造林主要选择在山坡中下部适宜人工造林的地段，这对200—500m海拔区段的植被增加起到了重要作用;同时，退耕还林、封山育林等生态建设工程对中上坡位、坡度较大的灌丛地和裸岩地采取全封的植被恢复策略，减少了人类活动对高海拔区域的干扰，而高海拔喀斯特地区植被主要以石山灌丛和次生林为主，植被恢复缓慢，从而使高海拔区域植被结构趋于稳定。坡度6—15°范围内的植被恢复最显著，而不同坡向间植被均恢复但差异不明显。由此说明，自2001年以来，国家相继在喀斯特地区实施的退耕还林工程、珠江防护林工程、石漠化治理试点工程等重点生态工程取得明显成效，喀斯特地区植被总体在恢复，生态环境趋于好转。另一方面，综合分析不同地形条件下喀斯特区域植被变化情况，发现在海拔400—500m、坡度大于25°的阳坡区域仍然存在植被减少的现象，这在一定程度上说明了喀斯特区域人类活动的干扰空间在扩大，趋向于陡坡地区[21- 22]。同时，国家林业局2012年公布的石漠化监测公报也表明喀斯特地区毁林开垦、樵采薪材的现象还较严重，陡坡耕种、过度放牧等现象还大量存在，使得部分区域仍然存在植被退化的现象。因此，后续喀斯特地区生态恢复与重建过程中要注意海拔200—500m范围内生态恢复成果的巩固，同时要特别加强海拔400—500m、坡度大于25°的阳坡区域的生态治理措施的实施力度。

目前为止，SPOT VGT NDVI资料的时间尺度还较短，同时，基于国家气象站点的气象数据插值存在一定的不确定性，针对喀斯特特殊地貌的气象插值方法也需要考虑更多因素才能获得精度更高的插值数据，因此，对于植被变化的影响因素(如气候变化、人类活动等)的分析还需要更长时间的遥感资料积累和气象插值方法上的改进。另一方面，本研究主要考虑了以生态建设工程为主的人类活动对植被变化的正面影响，没有考虑其负面影响，在今后的研究中应增加深入探讨人类生态建设对植被变化的负面影响，从而有利于更全面的评价喀斯特地区植被恢复成效。

4 结论

基于SPOT VGT NDVI时间序列数据和气候资料(温度和降水)，利用提出的地形面积差异修正系数，分析了桂西北1999—2010年来植被变化趋势及其与气候、地形因子的关系，揭示了不同地形因子条件下喀斯特植被恢复趋势差异，明确了植被恢复成效显著的区域及特别需要加强生态治理的区域，主要结论如下：

(1)近12年来，桂西北喀斯特地区植被变化总体上呈恢复趋势，约70%的区域植被呈增加趋势；区域尺度上，年均气候因子对桂西北喀斯特地区植被变化的作用不明显。

(2)海拔200—500m范围内的植被恢复显著(200—300m恢复最显著)，但此海拔范围区域也是植被覆盖变化易于波动的区域，而海拔大于800m的区域，植被覆盖变化趋于稳定；坡度6—15°范围内的植被恢复最显著，不同坡向上的植被均在恢复但差异不明显，且随着坡向由阴坡转阳坡，总体恢复呈减小趋势。

(3)由于人类干扰的负面影响，海拔400—500m、坡度大于25°的阳坡区域仍存在植被退化的现象，后续喀斯特地区植被恢复与重建过程中要注意海拔200—500m范围内生态恢复成果的巩固，同时要特别加强海拔400—500m、坡度大于25°的阳坡区域的生态治理措施的实施力度。

致谢：感谢中国西部环境与生态科学数据中心及中国气象科学数据共享服务网对本研究的数据支持。

[1] Wang S J.Concept deduction and its connotation of karst rocky desertification.Carsologica Sinica, 2002, 21(2): 101- 105.

[2] Liu C Q, Lang Y C, Li S L, Piao H C, Tu C L, Liu T Z, Zhang W, Zhu S F.Researches on biogeochemical processes and nutrient cycling in karstic ecological systems, southwest China: a review.Earth Science Frontiers, 2009, 16(6): 1- 12.

[3] Meng J J, Wang J.The response of vegetation dynamics to climate change in the southwestern karst region of China since the early 1980s.Geographical Research, 2007, 26(5): 857- 866.

[4] Zhang Y D, Zhang X H, Liu S R.Correlation analysis on normalized difference vegetation index (NDVI) of different vegetations and climatic factors in Southwest China.Chinese Journal of Applied Ecology, 2011, 22(2): 323- 330.

[5] Li H, Cai Y L, Chen R S, Chen Q, Yan X.Effect assessment of the project of grain for green in the karst region in Southwestern China: a case study of Bijie Prefecture.Acta Ecologica Sinica, 2011, 31(12): 3255- 3264.

[6] Liu T, Kinouchi T, Ledezma F.Characterization of recent glacier decline in the Cordillera Real by LANDSAT, ALOS, and ASTER data.Remote Sensing of Environment, 2013, 137: 158- 172.

[7] Jobbagy E G, Sala O E, Paruelo J M.Patterns and controls of primary production in the Patagonian steppe: a remote sensing approach.Ecology, 2002, 83(2): 307- 319.

[8] Boelman N T, Stieglitz M, Griffin K L, Shaver G R.Inter-annual variability of NDVI in response to long-term warming and fertilization in wet sedge and tussock tundra.Oecologia, 2005, 143(4): 588- 597.

[9] Schultz P A, Halpert M S.Global analysis of the relationships among a vegetation index, precipitation and land surface temperature.International Journal of Remote Sensing, 1995, 16(15): 2755- 2777.

[10] Piao S L, Fan J Y.Dynamic vegetation cover change over the last 18 years in China.Quaternary Sciences, 2001, 21(4): 294- 302.

[11] Chen J, Jonsson P, Tamura M, Gu Z H, Eklundh L, Matsushita B.A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter.Remote Sensing of Environment, 2004, 91(3/4): 332- 344.

[12] Yan Y H, He Z X, Yuan H, Xing F W.The ecological response of fern diversity to different slopes in Gudoushan Nature Reserve, Guangdong.Biodiversity Science, 2011, 19(1): 41- 47.

[13] Burn D H, Hag Elnur M A.Detection of hydrologic trends and variability.Journal of Hydrology, 2002, 261(1/4): 107- 122.

[14] Song T Q, Peng W X, Zeng F P, Wang K L, Qin W G, Tan W N, Liu L, Du H, Lu S Y.Spatial pattern of forest communities and environmental interpretation in Mulun National Nature Reserve, karst cluster-peak depression region.Chinese Journal of Plant Ecology, 2010, 34(3): 298- 308.

[15] Yang Y Z, Zhao P X, Hao H K, Chang M.Spatiotemporal variation of vegetation in northern Shaanxi of Northwest China based on SPOT-VGT NDVI.Chinese Journal of Applied Ecology, 2012, 23(7): 1897- 1903.

[16] Zhang B Q, Wu P T, Zhao X N.Detecting and analysis of spatial and temporal variation of vegetation cover in the Loess Plateau during 1982—2009.Transactions of the CSAE, 2011, 27(4): 287- 293.

[17] Li A, Wu J G, Huang J H.Distinguishing between human-induced and climate-driven vegetation changes: a critical application of RESTREND in inner Mongolia.Landscape Ecology, 2012, 27(7): 969- 982.

[18] Sun X P, Wang T M, Kou X J, Ge J P.Normalized difference vegetation index dynamic change and its driving factor analysis with long time series in the Jinghe River watershed on the Loess Plateau of China.Chinese Journal of Plant Ecology, 2012, 36(6): 511- 521.

[19] Li D K.Vegetation change and its response to climate in the hill and ravine region of the Loess Plateau in Northern Shaanxi.Acta Botanica Boreali-Occidentalia Sinica, 2009, 29(5): 1007- 1015.

[20] Jiang Z C.The study on rocky desertification and its comprehensive control of karst mountain area in Guangxi.Beijing: Science Press, 2011: 138- 140.

[21] Xu Y Q, Luo D, Feng Y, Peng J.A study on land use and land cover change in karst mountainous areas in southwestern China: a case study of the Maotiao River Watershed, Guizhou Province.Resources Science, 2010, 32(9): 1752- 1760.

[22] Yang Q Q, Wang K L, Chen H S, Zhang W, Tian R C.Effects of geology and landform on karst rock desertification: a case study in Dahua County of Guangxi, China.Journal of Mountain Science, 2009, 27(3): 311- 318.

参考文献：

[1] 王世杰.喀斯特石漠化概念演绎及其科学内涵的探讨.中国岩溶, 2002, 21(2): 101- 105.

[2] 刘丛强, 郎赟超, 李思亮, 朴何春, 涂成龙, 刘涛泽, 张伟, 朱书法.喀斯特生态系统生物地球化学过程与物质循环研究: 重要性、现状与趋势.地学前缘, 2009, 16(6): 1- 12.

[3] 蒙吉军, 王钧.20世纪80年代以来西南喀斯特地区植被变化对气候变化的响应.地理研究, 2007, 26(5): 857- 866.

[4] 张远东, 张笑鹤, 刘世荣.西南地区不同植被类型归一化植被指数与气候因子的相关分析.应用生态学报, 2011, 22(2): 323- 330.

[5] 李昊, 蔡运龙, 陈睿山, 陈琼, 严祥.基于植被遥感的西南喀斯特退耕还林工程效果评价- 以贵州省毕节地区为例.生态学报, 2011, 31(12): 3255- 3264.

[10] 朴世龙, 方精云.最近18年来中国植被覆盖的动态变化.第四纪研究, 2001, 21(4): 294- 302.

[12] 严岳鸿, 何祖霞, 苑虎, 邢福武.坡向差异对广东古兜山自然保护区蕨类植物多样性的生态影响.生物多样性, 2011, 19(1): 41- 47.

[14] 宋同清, 彭晚霞, 曾馥平, 王克林, 覃文更, 谭卫宁, 刘璐, 杜虎, 鹿士杨.木论喀斯特峰丛洼地森林群落空间格局及环境解释.植物生态学报, 2010, 34(3): 298- 308.

[15] 杨延征, 赵鹏祥, 郝红科, 常鸣.基于SPOT-VGT NDVI的陕北植被覆盖时空变化.应用生态学报, 2012, 23(7): 1897- 1903.

[16] 张宝庆, 吴普特, 赵西宁.近30a黄土高原植被覆盖时空演变监测与分析.农业工程学报, 2011, 27(4): 287- 293.

[18] 孙晓鹏, 王天明, 寇晓军, 葛剑平.黄土高原泾河流域长时间序列的归一化植被指数动态变化及其驱动因素分析.植物生态学报, 2012, 36(6): 511- 521.

[19] 李登科.陕北黄土高原丘陵沟壑区植被覆盖变化及其对气候的响应.西北植物学报, 2009, 29(5): 1007- 1015.

[20] 蒋忠诚.广西岩溶山区石漠化及其综合治理研究.北京: 科学出版社, 2011: 138- 140.

[21] 许月卿, 罗鼎, 冯艳, 彭建.西南喀斯特山区土地利用/覆被变化研究——以贵州省猫跳河流域为例.资源科学, 2010, 32(9): 1752- 1760.

[22] 杨青青, 王克林, 陈洪松, 张伟, 田日昌.地质地貌因素对喀斯特石漠化的影响——以广西大化县为例.山地学报, 2009, 27(3): 311- 318.

TrendsinvegetationandtheirresponsestoclimateandtopographyinnorthwestGuangxi

TONG Xiaowei1,2,3, WANG Kelin1,2,*,YUE Yuemin1,2, LIAO Chujie1,2,3, XU Yanfang1,2,3, ZHU Haitao4

1KeyLaboratoryofAgro-ecologicalProcessesinSubtropicalRegion,InstituteofSubtropicalAgriculture,theChineseAcademyofSciences,Changsha410125,China2HuanjiangObservationandResearchStationforKarstEco-systems,Huanjiang547100,China3UniversityofChineseAcademyofScience,Beijing100049,China4SatelliteEnvironmentCent,MinistryofEnvironmentProtection,Beijing100094,China

The annual dynamics and spatial variability of vegetation in northwest Guangxi were studied based on the SPOT NDVI datasets during 1999—2010.The response of vegetation to climate and topography was analyzed.The results showed that: (1) Vegetation had a large proportion of increasing trends in northwest Guangxi from 1999 to 2010.Partial correlation analysis between NDVI and climate parameters showed that climate had little effect on vegetation.(2) Elevation affected vegetation trend.It increased between 200 and 500 m, however, there were a small areas between 400 m and 500 m where vegetation decreased.Vegetation tended to be stable as the elevation increasing.(3) Vegetation recuperated on slopes between 6°and 15°, but decreased on slopes between 2°and 6° and above 25°.(4) The differences of vegetation restoration on aspects were not obvious.However, vegetation recovery trended to be decreasing as the aspect turning from shady to sunny slopes.Our study indicated that the implementation of environmental conservation policies and ecological construction had resulted in the improvement of vegetation in karst regions.However, there were some regions with vegetation degradation caused by human disturbance.In the further implementation of ecological restoration, the government should pay more attention to the sunny regions where the elevation between 400m and 500m and the steep slope that larger than 25°.

SPOT-VGT NDVI; climate change; correction coefficient for topographic areas; ecological construction projects；karst.

中国科学院西部行动计划项目(KZCX2-XB3- 10); 国家科技支撑计划课题(2010BAE00739- 02); 国家自然科学基金项目(41371418)

2013- 10- 16;

2014- 04- 23

10.5846/stxb201310162503

*通讯作者Corresponding author.E-mail: kelin@isa.ac.cn

童晓伟，王克林，岳跃民，廖楚杰，徐艳芳，朱海涛.桂西北喀斯特区域植被变化趋势及其对气候和地形的响应.生态学报,2014,34(12):3425- 3434.

Tong X W, Wang K L,Yue Y M, Liao C J, Xu Y F, Zhu H T.Trends in vegetation and their responses to climate and topography in northwest Guangxi.Acta Ecologica Sinica,2014,34(12):3425- 3434.