EN中文

资源科学 >

2017 , Vol. 39 >Issue 10: 1975 - 1988

DOI: https://doi.org/10.18402/resci.2017.10.16

基于地理加权回归的石漠化影响因子分布研究

  • 许尔琪
展开
  • 中国科学院地理科学与资源研究所陆地表层格局与模拟重点实验室, 北京 100101

作者简介：许尔琪,男,广东汕头人,博士,主要从事土地利用及空间格局、生态环境效应研究。E-mail：

收稿日期: 2017-08-23

  要求修回日期: 2017-09-14

  网络出版日期: 2017-10-20

基金资助

国家自然科学基金项目（41601095）

国家重点基础研究发展计划项目（973 计划）（2015CB452702）

收起

Spatial variation in drivers of karst rocky desertification based on geographically weighted regression model

  • XU Erqi
Expand
  • Key Laboratory of Land Surface Pattern and Simulation,Institute of Geographic Sciences and Natural Resources Research,Chinese Academy of Sciences,Beijing 100101,China

Received date: 2017-08-23

  Request revised date: 2017-09-14

  Online published: 2017-10-20

Copyright

《资源科学》编辑部
Fold

摘要

分析和识别石漠化的关键影响因子,有助于有效治理和恢复石漠化。以往研究对影响因子的空间局部差异关注较少,本文以黔桂喀斯特山地为研究区,选取自然及社会经济等12个影响因子,利用地理加权回归（GWR）模型,在普通线性回归的基础上嵌入空间因素,分析石漠化影响因子的空间分异。结果表明:①黔桂喀斯特山地的石漠化Moran’s I大于正态函数在99%显著水平,存在明显空间聚集现象;②GWR模型R2（0.508）明显高于传统统计模型的R2（0.156）,回归模型拟合效果显著提高;③12个影响因子与石漠化关系呈现不同数值大小、正负效应和线性组合关系的空间分布差异;④人类活动叠加在喀斯特特殊的岩性、土壤和植被组合上,显著影响石漠化分布;同时,局部区域高强度人类活动导致石漠化的急剧变化。GWR模型可揭示石漠化影响因子的空间分异规律和局部的关键影响因子,刻画多因子组合作用对石漠化的影响,有助于差别化的小流域石漠化治理。

关键词： 石漠化; 影响因子; 地理加权回归模型; 空间分异; 黔桂喀斯特山地

本文引用格式

许尔琪 . 基于地理加权回归的石漠化影响因子分布研究[J]. 资源科学, 2017 , 39(10) : 1975 -1988 . DOI: 10.18402/resci.2017.10.16

Abstract

Analysis and identification of key drivers of karst rocky desertification （KRD）can contribute to effective management and restoration. Ignoring heterogeneity leads to statistical bias and influences the specificity of desertification control planning. Taking Guizhou and Guangxi Karst Mountainous Regions as the study area,this paper chose twelve factors,including socioeconomic,spatial distance,topography,climate,soil,lithology and land use,as drivers of KRD. Geographically Weighted Regression modeling （GWR）and embedding spatial factors to the traditional ordinary linear regression （OLS）model were used to analyze the spatial distribution of influence on KRD. An obvious spatial agglomeration of KRD in the study area according to Moran’s I at the significant level of 99% was found. Coefficients of determination （R2）of GWR were much higher for OLS （0.508 vs. 0.156）,indicating a much better fit for the GWR model. Coefficients of GWR between twelve drivers and an obvious spatial distribution of value magnitudes,negative or positive effects and combined types. The specifics of the karst background create a fragile and vulnerable environment that is susceptible to human activities. Meanwhile,intense human activities lead to a sharp change in KRD status,which included the predatory sabotage for casing the severe KRD and the KRD restoration projects for reversing KRD to no KRD. The regression coefficients of the twelve drivers and their linear combined characteristics showed different spatial distributions based on GWR modeling. Using the GWR model revealed the spatial discrimination of the effects of driving forces on KRD and identified key drivers of KRD at the local area,revealing the spatial distribution of the joint effect of different driving forces helping provide a scientific reference to differential KRD control at small watershed scales.

Key words： Karst rocky desertification; driving force; geographically weighted regression model; spatial distribution; Guizhou and Guangxi Karst Mountainous Region

1 引言

“石漠化”是指热带、亚热带脆弱的岩溶环境中,由于不合理人类干扰和破坏,发生植被退化、土壤侵蚀和基岩裸露等的土地退化现象[1-4]。世界岩溶地貌集中分布区主要包括有欧洲中南部、北美东部和东亚地区等三大片区。前两个片区上百年前就有石漠化,如土耳其、法国、摩洛哥、意大利、克罗地亚等等都曾经有石漠化[5-7],但由于上述地区地质环境脆弱性较小、人口和经济压力相对较轻,石漠化多为生态地质环境因素影响[8,9]。然而,地处东亚片区中心地带的西南岩溶地区,其石漠化问题不仅受地质背景的影响,更是自然因素和人类活动叠加的结果[10-14]。因此,分析和识别石漠化的关键驱动因素及其贡献程度,有助于治理和恢复石漠化。
石漠化吞噬西南岩溶地区民众的生存空间,影响珠江、长江的生态安全,制约区域经济社会可持续发展[10,15,16]。研究表明,石漠化的影响因子众多[14,17-20],自然因素包括地质条件、气候因素、地形因子和地貌类型等[21-24],人为因子既包括毁林开荒、坡地种植、超载过牧和开山开矿等加剧石漠化的行为,也包括生态修复工程等治理石漠化的因素[9,25-28]。不同案例区影响石漠化的关键因素各不相同,贡献率也有显著差异[14,18,19],但这些研究却很少关注石漠化与影响因子关系在案例区内的空间差异。考虑到石漠化治理和恢复多是以小流域为单元1）( 1）《中华人民共和国国民经济和社会发展第十个五年计划纲要》提出“加快小流域治理,减少水土流失,推进黔桂滇岩溶地区石漠化综合治理”。),如果缺乏对因子影响程度空间分布的认识,将影响石漠化治理规划和工程制定的准确性,尤其随着研究区面积的增加,空间差异和不确定性将进一步增加。因此,本研究选定黔桂喀斯特山地,采用地理加权回归（GWR）模型,辨析石漠化与影响因子关系的空间差异,为喀斯特地区石漠化治理提供科学依据。
Table 1 Classification standard of karst rocky desertification

表1 石漠化等级分级标准

石漠化强度等级 基岩裸露/% 植被+土壤覆盖/% 分级
无石漠化 <20 >80 1
潜在石漠化 20~30 70~80 2
轻度石漠化 30~50 50~70 3
中度石漠化 50~70 30~50 4
强度石漠化 70~90 10~30 5
极强度石漠化 >90 <10 6

2 研究区概况、研究方法与数据来源

2.1 研究区概况

本文以位于贵州和广西的黔桂喀斯特山地为研究区[29],该区覆盖峰林、峰丛和洼地等主要喀斯特地貌类型,面积为21.41万km2,约位于22º8'54"N-28º12'27"N,104º18'27"E-110º20'40"E之间（图1）。地势为西北高,东南低,山地多,平地少。该区属亚热带季风湿润气候,年均降雨约为800~1900mm,年均温度为9℃~23℃。作为中国石漠化分布的最严重区域,各类型石漠化皆有分布,是石漠化与影响因子关系研究的典型区。
Figure 1 Boundary of Guizhou and Guangxi Karst Mountainous Region

图1 黔桂喀斯特山地范围

2.2 研究方法与数据来源

2.2.1 石漠化分级
因石漠化分级标准较多,尚无统一定论,本文采用国家林业局官方发布的最新一期石漠化空间分布数据（2011年）[30],并对部分结果进行人工修正,该数据主要依据基岩裸露、植被覆盖和土壤覆盖程度进行石漠化分级[31]（表1）,将石漠化强度等级分为无石漠化、潜在石漠化、轻度石漠化、中度石漠化、强度石漠化和极强度石漠化等6级,分别用1~6表示。
2.2.2 影响因子选取
对比以往研究对石漠化影响因子贡献程度的结果[14,17-20],考虑数据的可获得性,本文共选取12个影响因子（图2）,包括社会经济、空间距离、地形、气候、土壤和土地利用等因素,其数据来源及处理过程见表2（见第1978页）。以1km×1km栅格单元对石漠化和影响因子进行空间重采样,共计21.41万个单元,进行后续分析。
Figure 2 Driving forces of karst rocky desertification in Guizhou and Guangxi Karst Mountainous Region

图2 黔桂喀斯特山地石漠化影响因子空间差异

Table 2 Driving forces of karst rocky desertification and data sources and processing in Guizhou and Guangxi Karst Mountainous Region

表2 黔桂喀斯特山地石漠化影响因子及其数据来源、处理过程

类型 影响因子（代码） 描述 数据来源及处理
社会经济 人口密度（x1） 空间单元的人口密度数据（2010年） 中国科学院资源环境科学数据中心,人口密度数据集[32]
GDP（x2） 空间单元的国民生产总值（2010年） 中国科学院资源环境科学数据中心,GDP数据[33]
空间距离 到主要铁路距离
（x3）		 空间单元中心到线状铁路的最近距离
（2010年）		 中国科学院资源环境科学数据中心,道路数据集[34]和距离分析工具（ArcGIS 10.1）
到主要公路距离
（x4）		 空间单元中心到线状公路的最近距离
（2010年）		 中国科学院资源环境科学数据中心,公路数据集[35]和距离分析工具（ArcGIS 10.1）
地形 海拔（x5） 空间单元中心点的海拔数值 中国科学院计算机网络信息中心地理空间数据云[36]
坡度（x6） 空间单元切平面与水平面夹角 中国科学院计算机网络信息中心地理空间数据云）[36]和地形分析工具（ArcGIS 10.1）
气候 年均降水（x7） 空间单元平均降水数值（1980-2011年） 寒区旱区科学数据中心中国区域高时空分辨率地面气象要素驱动数据集[37]
年均气温（x8） 空间单元平均气温数值（1980-2011年） 寒区旱区科学数据中心中国区域高时空分辨率地面气象要素驱动数据集[37]
植被 植被指数（x9） 空间单元的植被指数数据（2011年） NASA的陆地过程分布式数据档案中心[38]
立地 土壤砂粘百分比（x10） 空间单元土壤砂粘百分比 寒区旱区科学数据中心,中国土壤特征数据集[39]
碳酸盐岩比例分级（x11） 空间单元的碳酸盐岩比例分级数据 喀斯特数据中心[40],1级为连续性灰岩、白云岩,2级为灰岩白云岩混合岩,3级为灰岩夹碎屑岩、白云岩夹碎屑岩,4级为灰岩白云岩碎屑岩互层、灰岩碎屑岩互层,5级为碎屑岩
土地利用 土地利用影响分级（x12） 空间单元土地利用对石漠化影响分级
（2010年）		 中国科学院资源环境科学数据中心,土地利用数据集[41],1级为建设用地、水域、水田和未利用地,2级为林地,3级为旱地,4级为草地
2.2.3 数据预处理
为避免指标之间由于量纲和数量级的影响,本文采取Z-score标准化方法将石漠化分级及影响因子进行标准化处理,转化函数如下:
x ′ = x - μ σ （1）
式中 x ′ 为标准化后数值; x 为原始数值; μ 为样本数据均值; σ 为样本数据标准差。
采用容差及方差膨胀因子（VIF）对影响因子与石漠化关系进行共线性检验,以避免由于因素之间高度共线影响回归分析结果。各影响因子容差均大于0.1,方差膨胀因子（VIF）均介于0~10之间（表3）,因素之间不存在高度共线性,可直接进行回归分析。
Table 3 Collinearity test results between karst rocky desertification and driving factors in Guizhou and Guangxi Karst Mountainous Region

表3 黔桂喀斯特山地石漠化与影响因子共线性检验结果

影响因子 共线性统计量 影响因子 共线性统计量
容差 VIF 容差 VIF
x1 0.682 1.466 x7 0.554 1.805
x2 0.773 1.294 x8 0.225 4.439
x3 0.880 1.137 x9 0.673 1.485
x4 0.900 1.110 x10 0.962 1.040
x5 0.244 4.102 x11 0.951 1.051
x6 0.814 1.228 x12 0.931 1.074
2.2.4 全局空间自相关检验
选取Moran’s I指数检验石漠化分布的空间自相关性,判断其在空间上集聚的特点和平均集聚程度[42],其计算公式如下:
I = n ∑ i = 1 n ∑ j = 1 n W ij ( X i - X ̅ ) ( X j - X ̅ ) ∑ i = 1 n ∑ j = 1 n W ij ∑ i = 1 n ( X i - X ̅ ) 2 （2）
式中 i ， j 为所在单元;n为空间单元总数; W ij 为空间单元 i ， j 之间的影响程度,若相邻则值为1,不相邻为0; X i , X j 为 X 在相应空间单元 i 和 j 上的石漠化分级数值; X ̅ 为 X 的平均值。I值区间为[-1,1],值越趋近于1或者-1,表明空间分布差异性越大,接近0则代表单元间不相关。应用ArcGIS 10.1软件中的Moran’s I指数工具进行计算。
2.2.5 地理加权回归模型
地理加权回归模型（GWR）是对普通线性回归-最小二乘法（OLS）的扩展,将数据地理位置镶嵌到回归模型之中[43,44],可体现空间分异规律,反映自变量与因变量在空间上的相互影响关系。其计算公式如下:
y i = β 0 u i ， v i + ∑ j = 1 p β j u i ， v i x ij + ε i （3）
式中i为样本点;yi为样本点i的因变量; （ui,vi）为样本点i的空间位置; β j u i ， v i 为连续函数 β j u ， v 在样本点i的值;P为样本点i独立变量的个数; ε i 为随机误差项。
利用ArcGIS10.1 软件中的GWR工具来实现模型构建,以高斯函数确定权重,并选取AIC方法确定最有效带宽。

3 结果及分析

3.1 黔桂喀斯特山地石漠化分布

黔桂喀斯特山地2011年在各石漠化等级均有所分布（图3）,其中,无石漠化和潜在石漠化面积最大,分别为13.59万km2和3.35万km2,各占总面积的63.48%和15.66%;其余类型石漠化面积较低,重度、中度和轻度石漠化分列三到五位,比例分别为8.23%、6.73%和4.84%,极重度石漠化面积最少,比重仅为1.06%。Moran’s I指数为0.56,标准化统计量Z（I）值为514.97,大于正态函数在99%显著水平时的临界值2.58,表明石漠化存在着显著空间集聚现象。
Figure 3 Spatial distribution of karst rocky desertification in Guizhou and Guangxi Karst Mountainous Region in 2011

图3 2011年黔桂喀斯特山地石漠化空间分布

无石漠化是全区基底,其余等级石漠化在贵州和广西的空间分布有一定差异。广西范围石漠化的分布相对集中且等级差别相对明显,广西东南部和西部石漠化等级相对较低,主要以无石漠化等级分布为主,而中部石漠化等级相对较高,重度石漠化和极重度石漠化等级集中分布在该区域。贵州范围内石漠化的空间分布格局更为分散,以潜在石漠化和轻度石漠化两种等级分布为主,在各个区域皆有所分布,中度以上等级石漠化则多分布在西部。

3.2 石漠化GWR模型拟合结果

GWR模型与OLS模型相比,其拟合效果显著提高（表4）,GWR的拟合优度R2（0.508）明显高于OLS中的R2（0.156）,GWR的AIC（42 551.165）远小于比OLS的AIC（564 901.811）。同时,GWR模型的标准化残差值的范围在[-3.37,3.65],其中约99.76%的范围在[-2.58,2.58],表明标准化残差值在99%的显著性水平下是随机分布。
Table 4 Parameter estimation and test results of the GWR model

表4 GWR 模型参数估计及检验结果

模型参数 OLS模型 GWR模型
Bandwidth - 6 928.466
Residual Squares - 9 823.062
Effective Number - 1 678.221
Sigma - 0.779
AIC 564 901.811 42 551.165
R2 0.156 0.554
Adjusted R2 0.156 0.508
GWR模型对自变量在每个空间位置的回归系数进行局部统计,分别得到回归系数的统计特征（表5）以及直方图累积分布（图4）。人口密度（x1）、到主要公路的距离（x4）、海拔（x5）、坡度（x6）、年均气温（x8）、植被指数（x9）、土壤砂粘百分比（x10）、碳酸盐比例分级（x11）和土地利用分级（x12）等因素对石漠化影响程度的波动区间较小,GDP（x2）、到主要铁路的距离（x3）、年均降水（x7）等因素对石漠化的影响程度波动区间较大。因子回归系数均呈现出正负变动幅度,主要集中在[-1,1]范围内,围绕在0上下浮动。除人口密度、GDP、海拔、坡度和植被指数等5个因子有一定偏态分布,其余因子呈较明显正态分布。GDP、植被指数、碳酸盐岩百分比等因子回归系数小于0的比例较大,其他因素与石漠化的回归系数大于0的比重较大。
Table 5 Statistical characteristics of GWR regression coefficients for driving forces of karst rocky desertification in Guizhou and Guangxi Karst Mountainous Region

表5 黔桂喀斯特山地石漠化影响因子的GWR模型回归系数的描述性统计

影响因子 平均值 标准差 最大值 最小值 下四分位值 中位值 上四分位值
x1 -0.064 0.191 1.230 -1.734 -0.148 0.033 0.089
x2 -0.018 0.596 14.950 -12.900 -0.078 -0.012 0.018
x3 0.022 1.235 14.126 -18.324 -0.489 0.001 0.534
x4 0.024 0.197 1.049 -1.112 -0.068 0.001 0.115
x5 0.179 0.585 5.539 -2.585 -0.020 0.112 0.557
x6 0.082 0.106 0.835 -0.351 0.007 0.058 0.135
x7 -0.028 1.310 22.305 -18.087 -0.781 0.007 0.817
x8 0.040 0.565 3.670 -6.051 -0.187 0.002 0.233
x9 -0.136 0.134 0.349 -0.850 -0.098 -0.053 0.011
x10 0.008 0.163 1.105 -1.962 -0.062 0.003 0.082
x11 -0.098 0.193 1.417 -1.727 -0.183 -0.049 0.002
x12 0.010 0.080 0.527 -0.471 -0.027 0.001 0.044
Figure 4 Histogram statistics of regression coefficients of GWR for different driving forces in Guizhou and Guangxi Karst Mountainous Region

图4 黔桂喀斯特山地石漠化影响因子回归系数统计

3.3 石漠化影响因子的空间分异规律

利用GWR模型,分析石漠化影响因子回归系数的空间分布（图5）,结果表明,每一因子回归系数呈现不同的正负和数值的空间变化,且各因子系数的空间分布各不相同。一方面,不合理人类活动可加剧石漠化,有效保护的行为和工程则可治理石漠化,这使得量化人类活动强弱的驱动因子与石漠化在不同空间位置可能正负效应和影响大小的差异。另一方面,因子多以联合作用影响石漠化的分布,由于某个关键因子在局部的影响程度以及因子间的相关性,使得其他影响因子在回归方程中系数发生变化。因此,影响因子与石漠化关系,在空间上呈现出影响程度、正负效应和因子组合的显著差异。
Figure 5 Spatial distribution of GWR regression coefficients for driving forces of karst rocky desertification in Guizhou and Guangxi Karst Mountainous Region

图5 黔桂喀斯特山地石漠化影响因子回归系数

3.3.1 社会经济因子对石漠化影响的空间分异规律
人口密度回归系数呈现“中间低南北高”分布,与石漠化呈明显正相关区域主要分布在黔西、黔北、桂西和桂南,人口集聚加大了对土地的压力,加剧石漠化。呈明显负相关区域,集中分布在黔南及桂中西的天等、兴宾和上林等县（市）,研究表明,很多石漠化严重等级区域的人口较少,生态环境较差[45],石漠化受自然因子（如岩性、土壤和坡度）的影响更大,使得人口密度回归系数为负值。
GDP对石漠化的影响程度变化幅度较大（图5）,但仍主要集中在[-1,1]之间（图4）,回归系数总体呈现“西低中东高”格局。黔东,桂东以及黔桂交界的GDP与石漠化呈明显正相关,区域经济发展的提高以资源掠夺和破坏环境为代价,如桂西北红水河流域建立高能耗的林板厂[46],大量砍伐木材,严重破坏植被,加剧石漠化;呈较大负相关区域主要分布在黔西、桂中和桂西南,生态建设与产业发展和农民增收相结合,如黔西南的“晴隆模式”、“坪上模式”和“顶坛模式”等[47],既增加区域GDP,又增大对石漠化治理投入,从而,有效减轻对土地压力和治理石漠化,使得GDP与石漠化关系以负相关为主。
3.3.2 空间距离因子对石漠化影响的空间分异规律
两类空间距离因子对石漠化影响无明显空间分布规律,距离因子影响了人类活动的可达性,距离因子小,表明人为干扰活动范围和强度大,同时,石漠化治理工程实施可能性也增加,使得因子与石漠化正负相关关系并存[48];随着距离因子增大,人为活动减弱,以距离因子与其他因子复杂组合影响为主。到主要铁路距离回归系数在贵州极值分布少,而广西波动较大。在贵州修文、清镇以及广西东兰、巴马、罗城和环江等县（市）内,因子与石漠化正相关较大;负相关较大区域主要集中广西河池、都安、上林、靖西、凌云和那坡等县（市）内。
各县（市）主要公路可达性较好,道路的影响较为相似,因此回归系数波动较小（表5）,正负交错分布。到主要公路距离与石漠化正相关较大区域（回归系数为0.329~1.049）主要集中在贵州凯里、织金、康宁、紫云,以及广西河池、宜州、南丹、忻城、那坡和大新等县（市）内（图5d）;负相关较大区域则零散分布在黔西北以及桂中、桂南,这些区域道路可达性相对较差,以土壤和岩性因子的影响为主。
3.3.3 地形因子对石漠化影响的空间分异规律
海拔和坡度回归系数变化均相对较小,分别在
-2.585~5.539和-0.351~0.835范围内（表5）,呈“东南高西北低”的空间分布（图5e,图5f）。海拔和坡度影响植被和土层的分布,随着海拔和坡度增加,植被生境变差,土层变薄,容易加剧石漠化。因此,海拔与石漠化关系以正相关为主,正相关较大区域主要集中在贵州六盘水、西秀以及广西中东部;零星区域因子与石漠化负相关较大,包括贵州大方、兴仁、安龙和广西都安、德保、武宜、天等、兴宾等县（市）,随着海拔升高,人类活动受到限制,有利于石漠化改善。
坡度回归系数以正值为主,负值多分布在贵州。黔中和桂中东部该因子与石漠化正相关较大,负相关较大区域则主要分布在黔西、黔北和桂西南。坡度较低区域,坡度增加,发生水土流失风险越高,容易加剧石漠化,坡度与石漠化呈正相关关系;当坡度增加到一定程度,限制了人类活动,随着坡度增加,人类活动减弱,有助于石漠化恢复,坡度与石漠化为负相关。如清镇市王家寨流域研究表明,石漠化严重程度随坡度增加而升高,超过27°后石漠化程度则逐渐降低[49]。
3.3.4 气候因子对石漠化影响的空间分异规律
年均降水和年均气温回归系数皆分布较为分散。贵州和广西中东部年均降水量与石漠化正相关较大,降水量增加,增强了对土层的冲刷动力,加剧水土流失和石漠化;负相关较大区域（回归系数为-18.087~-3.831）则分布在黔东北和黔桂交界（图5g）,区域降雨增加有利于植被生长,植被盖度的增加又有利于石漠化改善,使得两者呈负相关关系。
年均气温与石漠化正相关关系较大区域主要集中在贵州中部,以及广西河池、宜州、南丹、忻城、那坡等县（市）内,气温升高,加速喀斯特发育和土层干化,进而影响植被生长,产生石漠化;负相关关系较大区域主要集中在贵州南部以及广西中南部,上述区域降雨大、植被盖度高,三者相互影响值得气温与石漠化呈负相关。
3.3.5 植被因子对石漠化影响的空间分异规律
植被指数越大,表明植被生长状况越好,越有利于水土保持,因此,植被指数与石漠化的关系以负相关为主。回归系数的波动相对较小,变化范围分布在-0.850~0.349（表5）,呈现“东高西低”的空间分布格局（图5i）。与石漠化负相关较大区域主要集中在贵州的六盘水市辖区、大方、纳雍、关岭、晴隆、安龙、隆林以及广西德保、平果、田东、凌云等县（市）;正相关较大区域零散分布在贵州东部和广西东北部,这些区域植被指数数值较高且变化很小,包括人口密度、GDP等人为因子对石漠化的影响更大,使得回归关系呈现一定的弱正相关。
3.3.6 立地因子对石漠化影响的空间分异规律
土壤质地对石漠化影响的变化范围相对较小,土壤砂粘百分比回归系数分别在-1.962~1.105（表5）,分布较为零散,正负交错。研究区东部该因子与石漠化正相关较大;负相关较大区域则零散分布在贵州大方、长顺、惠水、独山和广西柳江等县（市）内,土壤质地的变化与石漠化存在一定非线性关系[50],使得该因子无明显空间分布特征。
岩性因子回归系数变化幅度也较小,在-1.727~1.417之间,呈现“东高西低”格局（图5k）。碳酸盐岩比例分级与石漠化负相关较大区域主要集中在黔西南和桂西北,该因子级别越低,则碳酸盐岩比例越大,有利于岩溶发育,加剧石漠化;与石漠化正相关较大区域主要集中在黔东和桂东北,研究表明植被在特殊的喀斯特环境有很强适应能力,裸岩暴露夹缝中仍能生长[51],因此,上述区域人为因子贡献大,通过合理的保护和修复有效治理石漠化,促进植被恢复,使得该因子系数呈正值。
3.3.7 土地利用分级对石漠化影响的空间分异规律
土地利用分级回归系数范围为-0.471~0.527（表5）,对石漠化影响程度变化较小,呈“桂中高、其余低”分布格局（图5l）。该因子与石漠化在贵州晴隆、施秉以及广西中部多个县市呈较大正相关,土地利用强度大,容易导致水土流失,加剧石漠化;负相关较大区域分布零散,包括有贵州清镇、石阡、余庆、纳雍、紫云、独山,以及广西环江、罗城、凤山、东林、凌云和德保等县（市）内,脆弱的喀斯特是石漠化形成的主因,自然因子贡献大,使得该指标为负值。

4 结论

黔桂喀斯特山地石漠化Moran’s I指数为0.56,标准化统计量Z（I）值为514.97,大于正态函数在99%显著水平的数值,存在明显的空间集聚分布。GWR模型分析加入空间要素,克服了经典回归模型统计条件假设的缺陷,其R2（0.508）明显比OLS模型的R2（0.156）高,拟合效果显著提高。
GWR计算的各石漠化影响因子回归系数在正负区间皆有分布,主要集中在[-1,1]范围内,呈现不同空间分布:人口呈“中间低南北高”,GDP为“西低中东高”,植被指数和岩性为“西低东高”,海拔和坡度则是“东南高西北低”,土地利用分级呈“桂中高、其余低”,到主要铁路的距离、到主要公路的距离、年均降水、年均气温和土壤质地对石漠化的影响空间分布较为零散。
喀斯特地区地形破碎,岩性、土壤和植被等类型多样,加之复杂人类活动,对石漠化分布产生显著影响,体现出单因子的正负效应、多种因子的组合类型和贡献程度的空间明显差异,如社会经济发展对石漠化的正负效应,人口与海拔、坡度的组合关系,降雨与植被对石漠化的交互作用等。人类活动叠加在喀斯特特殊的岩性、土壤和植被构成上,更容易发生石漠化,同时,在局部范围,高强度破坏和资源掠夺可导致严重石漠化发生,高投入石漠化治理措施又能够有效治理石漠化。
因此,石漠化治理工程应考虑不同影响因子组合对石漠化的影响及其空间差异。应用GWR模型,分析各影响因子对石漠化影响程度和正负效应,揭示影响因子在不同空间位置的组合,辨析其中影响石漠化的关键影响因子,从而为差别化小流域治理提供科学依据,以有效改善和控制石漠化。

The authors have declared that no competing interests exist.

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
熊平生,袁道先,谢世友. 我国南方岩溶山区石漠化基本问题研究进展[J]. 中国岩溶,2010,29(4):355-362.

[Xiong P S,Yuan D X,Xie S Y.Progress of research on rocky desertification in South China Karst Mountain[J]. Carsologica Sinica,2010,29(4):355-362.]

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

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

[3]
袁道先,蔡桂鸿. 岩溶环境学[M]. 重庆:重庆出版社,1988.

[Yuan D X,Cai G H.Karst Environment[M]. Chongqing:Chongqing Publishing House,1988.]

[4]
李阳兵,王世杰,容丽. 关于喀斯特石漠和石漠化概念的讨论[J]. 中国沙漠,2004,24(6):689-695.

[Li Y B,Wang S J,Rong L.Discuss on concepts of rock desertification and rock desert in karst region[J]. Journal of Desert Research,2004,24(6):689-695.]

[5]
Parise M,Pascali V.Surface and subsurface environmental de-gradation in the karst of Apulia(southern Italy)[J]. Environ-mental Geology,2003,44(3):247-256.

[6]
Gams I.Origin of the term “karst,” and the transformation of the classical karst(kras)[J]. Environmental Geology,1993,21(3):110-114.

[7]
Sauro U.Human impact on the karst of the Venetian Fore-Alps,Italy[J]. Environmental Geology,1993,21(3):115-121.

[8]
李阳兵,罗光杰,白晓永,等. 典型峰丛洼地耕地,聚落及其与喀斯特石漠化的相互关系-案例研究[J]. 生态学报,2014,34(9):2195-2207.

[Li Y B,Luo G J,Bai X Y,et al.The correl-ations among arable land,settlement and karst rocky desertifi-cationcases study based on typical peak-cluster depression[J]. Acta Ecologica Sinica,2014,34(9):2195-2207.]

[9]
袁道先. 岩溶石漠化问题的全球视野和我国的治理对策与经验[J]. 草业科学,2008,25(9):19-25.

[Yuan D X.Global view on karst rock desertification and integrating control measures and experiences of China[J]. Pratacultural Science,2008,25(9):19-25.]

[10]
Wang S J,Li R L,Sun C X,et al.How types of carbonate rock assemblages constrain the distribution of karst rocky desertified land in Guizhou Province,PR China:phenomena and mechani-sms[J]. Land Degradation & Development,2004,15(2):123-131.

[11]
李阳兵,王世杰,容丽. 关于中国西南石漠化的若干问题[J]. 长江流域资源与环境,2003,12(6):593-598.

[Li Y B,Wang S J,Rong L.Problems of karst rocky desertification in southwest China[J]. Resources and Environment in the Yangtze Basin,2003,12(6):593-598.]

[12]
苏维词. 中国西南岩溶山区石漠化的现状成因及治理的优化模式[J]. 水土保持学报,2002,16(2):29-32.

[Su W C.Cont-rolling model for rocky desertification of karst mountainous region and its preventing strategy in Southwest,China[J]. Journal of Soil Water Conservation,2002,16(2):29-32.]

[13]
Yan X,Cai Y L.Multi-scale anthropogenic driving forces of karst rocky desertification in Southwest China[J]. Land Degrad-ation & Development,2015,26(2):193-200.

[14]
Xu E,Zhang H.Characterization and interaction of driving factors in karst rocky desertification:a case study from Chang-shun,China[J]. Solid Earth,2014,5(2):1329-1340.

[15]
Bai X Y,Wang S J,Xiong K N.Assessing spatial-temporal evolution processes of karst rocky desertification land:indi-cations for restoration strategies[J]. Land Degradation & Development,2011,24(1):47-56.

[16]
王世杰. 喀斯特石漠化-中国西南最严重的生态地质环境问题[J]. 矿物岩石地球化学通报,2003,22(2):120-126.

[Wang S J.The most serious eco-geologically environmental problem in southwestern China-karst rocky desertification[J]. Bulletin of Mineralogy Petrology and Geochemistry,2003,22(2):120-126.]

[17]
Jiang Y,Li L,Groves C,et al.Relationships between rocky desertification and spatial pattern of land use in typical karst area,Southwest China[J]. Environmental Earth Sciences,2009,59(4):881-890.

[18]
Liu Y S,Wang J Y,Deng X Z.Rocky land desertification and its driving forces in the karst areas of rural Guangxi,Southwest China[J]. Journal of Mountain Science,2008,5(4):350-357.

[19]
Yang Q Q,Wang K L,Zhang C H,et al.Spatio-temporal evolu-tion of rocky desertification and its driving forces in karst areas of Northwestern Guangxi,China[J]. Environmental Earth Sci-ences,2011,64(2):383-393.

[20]
Li Y,Xie J,Luo G,et al.The evolution of a karst rocky desertification land ecosystem and its driving forces in the Houzhaihe area,China[J]. Open Journal of Ecology,2015,5(10):501-512.

[21]
李瑞玲,王世杰,周德全,等. 贵州岩溶地区岩性与土地石漠化的相关分析[J]. 地理学报,2003,58(2):314-320.

[Li R L,Wang S J,Zhou D Q,et al.The correlation between rock deserti-fication and lithology in karst area of Guizhou[J]. Acta Geo-graphica Sinica,2003,58(2):314-320.]

[22]
李瑞玲,王世杰,熊康宁,等. 贵州省岩溶地区坡度与土地石漠化空间相关分析[J]. 水土保持通报,2006,26(4):82-86.

[Li R L,Wang S J,Xiong K N,et al.Correlation between rocky desertification and slope degree in karst area of Guizhou[J]. Bulletin of Soil and Water Conservation,2006,26(4):82-86.]

[23]
王世杰,李阳兵,李瑞玲. 喀斯特石漠化的形成背景,演化与治理[J]. 第四纪研究,2003,23(6):657-666.

[Wang S J,Li Y B,Li R L.Karst rocky desertification:Formation background,evolution and comprehensive taming[J]. Quaternary Sciences,2003,23(6):657-666.]

[24]
张信宝,王世杰,白晓永,等. 贵州石漠化空间分布与喀斯特地貌,岩性,降水和人口密度的关系[J]. 地球与环境,2013,41(1):1-6.

[Zhang X B,Wang S J,Bai X Y,et al.Relationships between the spatial distribution of karst land desertification and geomorphology,lithology,precipitation,and population density in Guizhou province[J]. Earth and Environment,2013,41(1):1-6.]

[25]
Qi X K,Wang K L,Zhang C H.Effectiveness of ecological restoration projects in a karst region of southwest China assessed using vegetation succession mapping[J]. Ecological Engi-neering,2013,54:245-253.

[26]
Zhang M,Wang K,Liu H,et al.How ecological restoration alters ecosystem services:an analysis of vegetation carbon sequestration in the karst area of northwest Guangxi,China[J]. Environmental Earth Sciences,2015,74(6):5307-5317.

[27]
Tong X,Wang K,Yue Y,et al.Quantifying the effectiveness of ecological restoration projects on long-term vegetation dynamics in the karst regions of Southwest China[J]. International Journal of Applied Earth Observation and Geoinformation,2017,54:105-113.

[28]
李阳兵,罗光杰,邵景安,等. 岩溶山地聚落人口空间分布与演化模式[J]. 地理学报,2013,67(12):1666-1674.

[Li Y B,Luo G J,Shao J A,et al.The evolving models of rural settlements and population distribution in karst mountains[J]. Acta Geographica Sinica,2013,67(12):1666-1674.]

[29]
许尔琪,张红旗. 喀斯特山地土地利用变化的垂直分布特征[J]. 中国生态农业学报,2016,24(12):1693-1702.

[Xu E Q,Zhang H Q.Vertical distribution of land use in karst moun-tainous region[J]. Chinese Journal of Eco-Agriculture,2016,24(12):1693-1702.]

[30]
中华人民共和国国家林业局. 中国石漠化状况公报[EB/OL].(2012-06-18)[2017-08-23].

[State Forestry Administr-ation of People 's Republic of China. Bulletin of China 's Karst Rock Desertification[EB/OL].(2012-06-18)[2017-08-23].

[31]
Li Y B,Shao J A,Yang H, et al. The relations between land useand karst rocky desertification in a typical karst area,China[J]. Environmental Geology,57(3):621-627.

[32]
中国科学院资源环境科学数据中心.人口密度数据(2010)[EB/OL].(2014-06-01)[2017-03-10].

[Data Center for Resources and Environmental Sciences,Chinese Academy of Sciences. Population Density Dataset (2010)[EB/OL]. (2014-06-01)[2017-03-10].

[33]
中国科学院资源环境科学数据中心. GDP数据[EB/OL]. (2014-06-01)[2017-03-10].

[Data Center for Resources and Environmental Sciences,Chinese Academy of Sciences. GDP Dataset (2010)[EB/OL]. (2014-06-01)[2017-03-10].

[34]
中国科学院资源环境科学数据中心. 道路数据集[EB/OL]. (2017-03-01)[2017-03-10].

[Data Center for Resources and Environmental Sciences,Chinese Academy of Sciences. Road Dataset [EB/OL]. (2017-03-01)[2017-03-10].

[35]
中国科学院资源环境科学数据中心. 公路数据集[EB/OL]. (2017-03-01)[2017-03-10].

[Data Center for Resources and Environmental Sciences,Chinese Academy of Sciences. Road Dataset [EB/OL]. (2017-03-01)[2017-03-10].

[36]
中国科学院计算机网络信息中心地理空间数据云. SRTM DEM数据 [EB/OL]. (2017-01-01)[2017-03-10].

[International Scientific & Technical Data Mirror Site,Computer Network Information Center,Chinese Academy of Sciences. SRTM DEM Dataset [EB/OL]. (2017-01-01)[2017-03-10].

[37]
寒区旱区科学数据中心.中国区域高时空分辨率地面气象要素驱动数据集(2011)[EB/OL]. (2016-06-06)[2017-03-10].

[Cold and Arid Regions Science Data Center at Lanzhou. China Meteorological Forcing Dataset (2011)[EB/OL]. (2016-06-06)[2017-03-10].

[38]
The Land Processes Distributed Active Archive Center,LP DAAC/NASA.MODIS/Aqua Vegetation Indices[EB/OL].(2017-03-01)[2017-03-10].

[39]
寒区旱区科学数据中心. 中国土壤特征数据集 [EB/OL]. (2013-06-07)[2017-03-10].

[Cold and Arid Regions Science Data Center at Lanzhou. Soil Particle-size Distribution Dataset (2011)[EB/OL]. (2013-06-07)[2017-03-10].

[40]
喀斯特数据中心. 岩性数据[EB/OL]. (2011-10-30)[2017-03-10].

[Karst Data Center. Lithology Dataset[EB/OL]. (2011-10-30)[2017-03-10].

[41]
中国科学院资源环境科学数据中心. 土地利用数据(2010)[EB/OL]. (2012-02-20)[2017-03-10].

[Data Center for Resources and Environ-mental Sciences,Chinese Academy of Sciences. Land Use Dataset(2010)[EB/OL]. (2012-02-20)[2017-03-10].

[42]
Sawada M. Global Spatial Autocorrelation Indices-Moran's I,Geary's C and the General Cross-Product Statistic[EB/OL].(2009-11-23)[2017-08-23].

[43]
Brunsdon C,Fotheringham A S,Charlton M E.Geographically weighted regression:a method for exploring spatial nonstation-arity[J]. Geographical Analysis,1996,28(4):281-298.

[44]
Fotheringham A S,Brunsdon C,Charlton M E.Geographically Weighted Regression:The Analysis of Spatially Varying Re-lationships[M]. Chichester:John Wiley & Sons,2003.

[45]
白晓永,熊康宁,李阳兵,等. 喀斯特山区不同强度石漠化与人口因素空间差异性的定量研究[J]. 山地学报,2006,24(2):242-248.

[Bai X Y,Xiong K N,Li Y B,et al.Quantitative research of the space dependence of different grades rocky desertification and the factor of population in karst mountain area[J]. Journal of Mountain Science,2006,24(2):242-248.]

[46]
覃勇荣,蓝崇钰,束文圣,等. 广西石灰岩植被破坏造成的石漠化问题[J]. 热带林业,2007,35(S1):48-51.

[Qin Y R,Lan C Y,Shu W S,et al.Problem of karst desertification caused by destru-ction of limestone vegetation in Guangxi,China[J]. Tropical Forestry,2007,35(S1):48-51.]

[47]
邓家富. 黔西南州石漠化治理的主要做法及成功模式[J]. 中国水土保持,2014,34(1):4-7.

[Deng J F.Main practice and successful mode of karst rocky desertification control in Qianxinan[J]. Soil and Water Conservation in China,2014,34(1):4-7. ]

[48]
Xu E Q,Zhang H Q,Li M X. Mining spatial information to investigate the evolution of karst rocky desertification and its human driving forces in Changshun,China[J]. Science of the Total Environment,2013,458-460:419-426.

[49]
李阳兵,王世杰,周梦维,等. 不同空间尺度下喀斯特石漠化与坡度的关系[J]. 水土保持研究,2009,16(5):70-77.

[Li Y B,Wang S J,Zhou M W,et al.Relationship between karst rocky desertification distribution and slope under different spatial scales[J]. Research of Soil and Water Conservation,2009,16(5):70-77.]

[50]
盛茂银,刘洋,熊康宁. 中国南方喀斯特石漠化演替过程中土壤理化性质的响应[J]. 生态学报,2013,33(19):6303-6313.

[Sheng M Y,Liu Y,Xiong K N.Response of soil physical-chemical properties to rocky desertification succession in South China Karst[J]. Acta Ecologica Sinica,2013,33(19):6303-6313.]

[51]
谭成江,冉景丞,莫家伟,等. 茂兰保护区石漠化现状,成因及治理对策[J]. 山地农业生物学报,2011,30(5):440-442.

[Tan C J,Rang J C,Mo J W,et al.Status,causes and control measures of stony desertificaion in Maolan Nature Reserve[J]. Journal of Moutain Agriculture and Biology,2011,30(5):440-442.]

Options
文章导航

/

〈 〉