1 研究方法与数据来源
1.1 研究区概况与数据来源
1.1.1 研究区概况
1.1.2 数据来源
表1 所用数据类型及数据来源Tab.1 Data types and data sources used |
| 数据类型 | 数据来源(年份) | 数据描述 | 分辨率 |
|---|---|---|---|
| 土地利用数据 | 中国科学院资源环境科学数据中心(2020年) | 土地利用类型 | 30 m |
| 气象数据 | 国家地球系统科学数据中心(2020年) | 年均降水、年均温 | 1 km |
| 土壤数据 | HWSD官网(2023年) | 土壤类型、深度等 | 1 km |
| DEM数据 | 地理空间数据云(2015年) | 高程 | 30 m |
| 人口密度数据 | Worldpop官网(2020年) | 人口密度 | 100 m |
| 道路数据 | 国家基础地理信息中心(2020年) | 道路类型、道路长度等 |
1.2 研究方法
1.2.1 道路网络结构特征指标构建
1.2.2 景观格局破碎化特征指数计算
表2 景观格局指数及含义描述Tab.2 Landscape pattern indices and their descriptions |
| 维度 | 指数 | 描述 | 计算公式 |
|---|---|---|---|
| 数量 | 斑块密度(DP) | 度量单位面积某种景观类型斑块的数量 | DP= |
| 面积 | 最大斑块指数(ILP) | 度量景观内最大斑块占景观总面积的比例 | ILP= |
| 形状 | 景观形状指数(ILS) | 反映整个景观内斑块形状特点 | ILS= |
| 聚合 | 聚合度(IA) | 度量区域景观格局中不同景观组分斑块的聚集程度 | IA= |
| 连通 | 凝聚度指数(ICOHESION) | 度量生态景观类型的自然连通度 | ICOHESION= |
注:N为某景观类型的斑块总数;A为景观总面积;aij为斑块ij的面积;E为景观中两斑块类型之间的边缘总长度;Pij为斑块ij的周长;m为景观中斑块类型总数。 |
1.2.3 回归模型及变量筛选
表3 自变量维度及名称描述Tab.3 Dimensions and names of the independent variable |
| 指标维度 | 指标描述 | 指标名称 |
|---|---|---|
| 自然基底 | 地形条件 | 高程(x1) |
| 坡度(x2) | ||
| 气象条件 | 年均气温(x3) | |
| 年均降雨量(x4) | ||
| 土壤条件 | 土壤保持量(x5) | |
| 道路特征 | 道路密度 | 高速道路密度(x6) |
| 高等级道路密度(x7) | ||
| 低等级道路密度(x8) | ||
| 距道路距离 | 距公路距离(x9) | |
| 距高速道路距离(x10) | ||
| 距高等级道路距离(x11) | ||
| 距低等级道路距离(x12) | ||
| 道路形态 | 道路曲率(x13) | |
| 道路网络连通性 | 度中心性(x14) | |
| 中介中心性(x15) | ||
| 平均最短路径长度(x16) | ||
| 人口分布 | 分布范围 | 建设面积占比(x17) |
| 分布强度 | 平均人口密度(x18) |
2 结果分析
2.1 陕南秦巴山区道路网络空间分布
2.2 陕南秦巴山区道路网络结构特征
2.3 陕南秦巴山区景观格局破碎化多维特征空间分布
2.4 道路结构特征对景观破碎化的差异性影响
2.4.1 回归模型构建
表4 五维度景观破碎化特征探索性回归结果Tab.4 Exploratory regression results for five-dimensional landscape fragmentation characteristics |
| 维度 | 因变量 | 校正R2 | AIC | JB | VIF | SA | 自变量 |
|---|---|---|---|---|---|---|---|
| 数量 | 斑块密度(DP) | 0.64 | 5 789.14 | 0 | 1.29 | 0 | -x1-x2+x6+x7+x8-x13+x14-x9 |
| 形状 | 景观形状指数(ILS) | 0.48 | 3 101.42 | 0 | 1.33 | 0 | x3-x2-x6-x13+x14-x12 |
| 面积 | 最大斑块面积(ILP) | 0.39 | 19 970.25 | 0 | 1.25 | 0 | x1+x2-x8-x14+x9 |
| 聚合 | 聚合度(IA) | 0.56 | 2 077.39 | 0 | 1.17 | 0 | -x3+x2-x7-x8+x11+x12 |
| 连通 | 凝聚度(ICOHESION) | 0.52 | -4 919.63 | 0 | 1.82 | 0 | x1+x2-x6-x7-x8 |
表5 五维度景观破碎化特征不同回归模型拟合结果对比Tab.5 Comparison of fit results for different regression models of five-dimensional landscape fragmentation characteristics |
| 维度 | 普通最小二乘回归(OLS) | 地理加权回归(GWR) | 多尺度地理加权回归(MGWR) | ||||||
|---|---|---|---|---|---|---|---|---|---|
| R2 | 校正R2 | AIC | R2 | 校正R2 | AIC | R2 | 校正R2 | AIC | |
| 数量 | 0.66 | 0.66 | 5 821.98 | 0.80 | 0.75 | 5 138.29 | 0.82 | 0.80 | 3 082.52 |
| 形状 | 0.46 | 0.46 | 3 133.64 | 0.64 | 0.62 | 2 309.55 | 0.72 | 0.70 | 4 058.22 |
| 面积 | 0.40 | 0.40 | 20 442.90 | 0.56 | 0.51 | 5 152.31 | 0.55 | 0.53 | 5 076.14 |
| 聚合 | 0.58 | 0.58 | 2 091.24 | 0.74 | 0.73 | 1 010.54 | 0.83 | 0.82 | 2 886.67 |
| 连通 | 0.50 | 0.50 | -4 961.87 | 0.67 | 0.64 | -5 681.49 | 0.71 | 0.68 | 4 220.71 |
2.4.2 关键影响因素分析
表6 各维度景观破碎化特征多尺度地理加权回归结果Tab.6 MGWR results for quantitative features of landscape fragmentation |
| 维度 | 解释变量 | 带宽/m | 回归系数平均值 | 变量解释度 |
|---|---|---|---|---|
| 数量 | 高程 | 43 924.38 | -0.42 | 82.41 |
| 距公路距离 | 51 422.74 | -0.35 | 76.83 | |
| 低等级道路密度 | 43 924.38 | 0.27 | 69.56 | |
| 坡度 | 201 529.08 | -0.18 | 100 | |
| 高等级道路密度 | 126 274.20 | 0.13 | 79.92 | |
| 形状 | 年均温 | 43 924.38 | 0.28 | 75.40 |
| 坡度 | 43 924.38 | -0.22 | 45.15 | |
| 道路曲率 | 43 924.38 | -0.15 | 40.82 | |
| 度中心性 | 425 241.26 | 0.18 | 100 | |
| 面积 | 距公路距离 | 58 824.76 | 0.65 | 98.47 |
| 坡度 | 74 623.51 | 0.16 | 45.22 | |
| 低等级道路密度 | 96 715.41 | -0.11 | 27.68 | |
| 度中心性 | 425 241.26 | -0.14 | 100 | |
| 聚合 | 年均温 | 43 924.38 | -0.47 | 77.21 |
| 坡度 | 43 924.38 | 0.35 | 75.14 | |
| 连通 | 高程 | 43 924.38 | 0.21 | 33.20 |
| 坡度 | 43 924.38 | 0.27 | 40.52 | |
| 高等级道路密度 | 43 924.38 | -0.22 | 38.33 | |
| 低等级道路密度 | 43 924.38 | -0.32 | 80.25 |