Landscape Pattern Variations of Forest, Wetland and Grassland Resources and Driving Forces in Liaoning Province, China
-
摘要:
目的 比较研究辽宁省2005年和2015年的森林、湿地与草地景观格局变化及驱动力,为3类自然资源的保护工作提供科学支持。 方法 基于2005、2015年2期辽宁省土地利用数据,利用ArcGIS和Fragstats软件计算景观转移矩阵和景观指数,并分析各类资源转化的空间分布,应用主成分分析法对辽宁省森林、湿地与草地景观格局变化的驱动力进行研究。 结果 1)辽宁省森林资源呈现东部多于西部的分布格局,林地景观面积增加了4284 km2,破碎化程度减小,景观转移变化主要发生在西部地区;湿地资源分布呈现沿海大于内陆、集中于河流沿岸的分布格局,近10年来湿地景观面积减小了0.55%,湿地景观形状趋于复杂,单元趋于离散,破碎化程度加剧,景观转移变化主要发生在河流沿岸;草地资源呈现西部多于东部的分布格局,变化率为-5.00%,破碎化程度加剧,斑块数量、面积、密度的3重减少表明草地景观退化程度严重,转移变化面积大且在全省均有分布。2)与2005年相比,2015年景观整体破碎程度加深,离散小单元增多,异质性减弱,原有的景观格局被破坏,主导斑块的控制能力下降。3)辽宁省景观变化的主要驱动因素为人口、经济、政策以及自然环境因素等。 结论 受人为因素和自然因素影响,辽宁省森林景观面积增加,破碎化程度得到改善,而湿地和草地景观存在破碎化程度加深的现象。结合辽宁省森林、湿地、草地景观格局的变化与实际,建议在现有基础上继续加强对森林和湿地资源的保护,草地资源的保护可以借鉴森林资源保护经验,加强立法与保护,从而实现资源的保护和可持续利用。 Abstract:Objective In this study, we compared the landscape pattern variations between 2005 and 2015 for forest, wetland and grassland and driving forces in Liaoning Province, so as to provide scientific support for natural resources conservation. Methods With the two-phases land use data in 2005 and 2015, the landscape transition matrix and landscape index were calculated and the spatial distribution of resources transformation was analyzed using ArcGIS and Fragstats, the driving forces of the change of landscape patterns of forest, wetland and grassland in Liaoning province were studied using principal component analysis. Results The results showed: 1)the distribution area of forest resources was larger in the east than the west of Liaoning Province. The forest landscape area had been increased by 4284 km2, when forest fragmentation had been reduced and the transition of forest landscape primarily was observed in the western areas. The existence of wetland resources was larger in coastal areas than inland and riverside areas . The wetland landscape had been reduced by 0.55%, while the status of the wetland landscape became more unexpected especially riverside areas, overall degraded amount, area and density of wetlands in the recent decade. The grassland resources showed larger distribution in the west than that in the east, which had been lost by 5.00% in Liaoning Province. Worse than wetland fragmentation, grassland fragmentation had been intensified at multiple scales within the whole province. 2) The landscape fragmentation in Liaoning Province had been severer from 2005 to 2015, in specific with more discrete small units and decreased and finally the lost control of dominant patches. 3) The landscape changes in Liaoning Province could be driven by population, economy, policy and natural environment factors, etc. Conclusion Affected by human and natural factors, the area of forest landscape had been increased, however overall fragmentation had been improved in the whole province, due to wetland and grassland losses. Based on our study, sustainable forest and wetland conservation is suggested. The grassland conservation should be conducted by following the successful experiences and strengthening the legislation. -
Key words:
- land use /
- landscape pattern /
- driving force /
- resource conservation /
- Liaoning province
-
表 1 景观指数及描述
Table 1. Landscape index and description
景观指数
Landscape index生态学意义
Ecological Significance应用水平
Level计算公式
Calculation Formula斑块数 NP 表示景观中斑块的总数,用来衡量景观的复杂程度,斑块数越多说明景观越复杂。Indicates the total number of patches in the landscape, which is used to measure the complexity of the landscape. When the landscape become more complex, the NP is larger. 类型/景观
Class/LandscapeNP=N 斑块密度 PD 表示单位面积的斑块数,斑块密度越大,表示景观越破碎。Indicates the number of patches per unit area. When the landscape become more fragmented, the PD is larger. 类型/景观
Class/LandscapePD=N/A
式中:A为斑块面积。A is the area of patch.景观形状指数LSI 表示景观中斑块形状的规则程度,当景观中斑块形状不规则或偏离正方形时,LSI增大。Indicates the regularity of the patch shape in the landscape. When the patch shape in the landscape is irregular or deviates from the square, the LSI increases. 类型
Class${L}{S}{I}=\dfrac{\text{0.25}E}{\sqrt{A_{ } } }$
式中:E为景观中所有斑块边界总长度。E is the sum of the lengths (m) of all edge segments involving the corresponding patch type.散布与并列指标 IJI 表示各个斑块类型间的散布与并列情况,IJI接近0时表明斑块类型仅与少数几种其他类型相邻接。Indicates the distribution and juxtaposition of various patch types, IJI approaches 0 when the corresponding patch type is adjacent to only one other patch type. 类型
Class${I}{J}{I} = \dfrac{ { - \displaystyle\sum\limits_{i = 1}^m {\sum\limits_{k = i + 1}^m {\left[ {\left({\dfrac{ { {e_{ik} } } }{E} } \right) \cdot \ln \left({\dfrac{ { {e_{ik} } } }{E} } \right)} \right]} } } }{ {\ln \left[ {0.5m\left({m - 1} \right)} \right]} }\left({100} \right)$
式中:eik为斑块类型i和k之间边缘总长度,m为景观中存在的斑块类型数量。eik is the total length of edge in landscape between patch types i and k; m is the number of patch types present in the landscape.平均斑块面积AREA_MN 表示各类型景观的平均面积,与斑块数一起可以反映景观的破碎程度,斑块数量越多,平均斑块面积越小,景观越破碎。Indicates the average area of various types of landscapes, which can reflect the degree of fragmentation of the landscape together with the number of patches. When the landscape become more fragmented, the NP is larger and the AREA_MN is smaller. 类型
Class${A}{R}{E}{A}\_{M}{N}=\dfrac{A}{N}1{0}^{-6}$ 聚集度指数AI 表示景观水平上不同斑块类型的聚合程度,景观由较多离散小斑块组成时,聚集度较小;景观以连续大斑块或斑块高度连接时,聚集度较大。
Indicates the degree of aggregation of different patch types at the landscape level. When the landscape is composed of more discrete small patches, the degree of aggregation is smaller; when the landscape is highly connected by continuous large patches, the degree of aggregation is larger.类型/景观
Class/Landscape${A}{I}= \left[ {\dfrac{ { {g_{ii} } } }{ {\max \to {g_{ii} } } }{P_i} } \right]\left({100} \right)$
式中:gii为斑块类型i像素之间的相似连接的数量;max→gii为斑块类型i像素之间的相似连接的最大数量,Pi为斑块类型i组成的景观所占比例。gii is the number of like adjacencies between pixels of patch type,. max-gii is the maximum number of like adjacencies between pixels of patch type, Pi is the proportion of landscape comprised of patch type i.Shannon多样性指数 SHDI 可以表示景观的异质性,在景观系统中,土地利用越丰富,破碎化程度越高,SHDI值越高。Indicates the heterogeneity of the landscape. In the landscape system, when the land use is more intensive, the degree of fragmentation is higher and the SHDI is larger. 景观
Landscape${S}{H}{D}{I} = - \displaystyle\sum\limits_{i = 1}^m {\left[ { {P_i} \times \ln \left({ {P_i} } \right)} \right]}$ Shannon均匀度指数 SHEI 描述景观不同时期的多样性变化,该值越小优势度越高,说明景观受到一种或几种优势斑块类型支配。Indicate the diversity of the landscape in different periods. When the landscape is dominated by one or more dominant patch types, the SHEI is smaller and the degree of dominance is higher 景观
Landscape${S}{H}{E}{I} = \dfrac{ { - \displaystyle\sum\limits_{i = 1}^m {\left[ { {P_i} \times \ln \left({ {P_i} } \right)} \right]} } }{ {\ln m} }$ 蔓延度指数CONTAG 描述景观不同斑块的团聚程度或延展形式,该值越高表示某种优势斑块类型连接性越好。
Indicates the degree of agglomeration or extension of different patches in the landscape. When the connectivity of a certain dominant patch type is better, the CONTAG is higher.景观
Landscape${C}{O}{N}{T}{A}{G} = \left[ {1 + \displaystyle\sum\limits_{i = 1}^m {\sum\limits_{j = 1}^n {\dfrac{ { {P_{ij} }\ln \left({ {P_{ij} } } \right)} }{ {2\ln \left(m \right)} } } } } \right]\left({100} \right)$
式中: Pij是随机选择的两个相邻栅格细胞属于类型i与j的概率。Pij is the probability that two adjacent grid cells randomly selected belong to types i and j.表 2 辽宁省景观转移矩阵
Table 2. Transition matrix of landscape type of Liaoning province
耕地
Agricultural
area/km2林地
Forest area/km2草地
Grassland area/km2湿地
Wetland area/km2建筑用地
Build-up area/km2未利用地
Unused area/km2转入面积
Transfer-in area/km2转出面积
Transfer-out area/km2净变化量
Net change/km2变化率
Change rate/%耕地
Agricultural area42 339 5 613 1 048 2 204 3 758 194 10 233 12 817 −2 584 −0.47 林地
Forest area4 596 50 069 876 384 831 79 11 049 6 765 4 284 0.75 草地
Grassland area1 446 4 848 2 536 96 238 101 2 094 6 730 −4 636 −5.00 湿地
Wetland area2 581 323 87 10 747 978 496 3 284 4 465 −1 182 −0.55 建筑用地
Build-up area1 506 240 61 307 6 471 29 5 857 2 144 3 713 4.63 未利用地
Unused area104 25 22 292 52 883 898 493 405 2.94 表 3 辽宁省2005和2015年的景观水平景观格局指数
Table 3. Landscape index at the landscape level of Liaoning province in 2005 and 2015
年份 Year 斑块数 NP 斑块密度 PD 聚集度指数 AI Shannon多样性指数 SHDI Shannon均匀度指数SHEI 蔓延度指数CONTAG 2005 70 902 0.48 96.277 5 1.355 9 0.682 2 39.051 6 2015 83 864 0.57 96.143 1 1.333 2 0.727 2 34.346 5 表 4 辽宁省2005和2015年的类型水平景观格局指数
Table 4. Landscape index at the landscape level of Liaoning province in 2005 and 2015
年份
Year景观类型
Landscape type斑块数
NP斑块密度
PD景观形状指数
LSI散布与并列指标
IJI平均斑块面积
AREA_MN聚集度指数
AI2005 林地Forestland 16 042 0.109 5 218.665 0 47.737 5 354.367 3 97.260 8 湿地Wetland 4 507 0.030 6 142.289 8 70.679 0 319.959 8 96.568 5 草地Grassland 8 178 0.055 8 170.988 3 56.923 8 113.322 9 94.700 9 2015 林地Forestland 15 048 0.095 9 216.209 4 45.946 2 406.175 6 97.388 2 湿地Wetland 5 963 0.040 7 146.500 4 74.889 6 272.861 7 96.364 1 草地Grassland 5 645 0.038 4 139.038 3 67.814 4 82.016 9 93.910 3 表 5 特征值和主成分贡献率
Table 5. Eigenvalues and principal component contribution rate
成分
Components特征值
Eigenvalues贡献率
Proportion of
Variance/%累计贡献率
Cumulative
proportion/%1 8.512 85.119 85.119 2 1.089 10.893 96.012 3 0.224 2.235 98.247 4 0.121 1.209 99.456 5 0.032 0.323 99.779 6 0.017 0.172 99.952 7 0.004 0.041 99.993 8 0 0.004 99.997 9 0 0.003 100 10 1.22E-05 0 100 表 6 主成分载荷矩阵
Table 6. Loading matrix of principal component
变量
Variable第一主成分
Component 1第一主成分
Component 2X1 0.916 0.134 X2 −0.973 0.006 X3 0.999 0.003 X4 0.997 0.025 X5 0.988 0.101 X6 0.98 −0.105 X7 0.992 0.071 X8 0.986 0.018 X9 0.031 0.993 X10 0.896 0.311 -
[1] 布仁仓, 胡远满, 常禹, 等, 2005. 景观指数之间的相关分析[J]. 生态学报, 25(10): 2764-2775. doi: 10.3321/j.issn:1000-0933.2005.10.044 [2] 蔡为民, 唐华俊, 吕钢, 等, 2006. 景观格局分析法与土地利用转换矩阵在土地利用特征研究中的应用[J]. 中国土地科学, 20(1): 75-80. [3] 陈利顶, 李秀珍, 傅伯杰, 等, 2014. 中国景观生态学发展历程与未来研究重点[J]. 生态学报, 34(12): 3129-3141. [4] 冯志贤, 张继贤, 侯伟, 等, 2017. 基于地表覆盖分类的生态环境人为干扰度分析:以北京市为例[J]. 生态学杂志, 28(8): 2611-2620. [5] 盖美, 吴慧歌, 曲本亮, 2016. 新一轮东北振兴背景下的辽宁省水资源利用效率及其空间关联格局研究[J]. 资源科学, 38(7): 1336-1349. [6] 国家林业局, 2018. 2016退耕还林工程生态效益监测国家报告 [M]. 北京: 中国林业出版社. [7] 郭少壮, 白红英, 孟清, 等, 2020. 秦岭地区林地与草地景观格局变化及其驱动因素[J]. 生态学报, 40(1): 130-140. [8] 辽宁省统计局, 2006. 辽宁统计年鉴 [M]. 北京: 中国统计出版社. [9] 辽宁省统计局, 2016. 辽宁统计年鉴 [M]. 北京: 中国统计出版社. [10] 任嘉衍, 刘慧敏, 丁圣彦, 等, 2017. 伊河流域景观格局变化及其驱动机制[J]. 应用生态学报, 28(8): 2611-2620. [11] 孙凤云, 刘淼, 胡远满, 等, 2013. 基于多种方法的景观格局动态变化综合分析:以辽宁省铁岭市为例[J]. 生态学杂志, 32(8): 2163-2171. [12] 王兵, 鲁绍伟, 尤文忠, 等, 2010. 辽宁省森林生态系统服务价值评估[J]. 应用生态学报, 21(7): 1792-1798. [13] 王美玲, 张继超, 王舶鉴, 等, 2017. 长白山区森林景观格局动态[J]. 生态学杂志, 36(11): 3138-3147. [14] 徐新良, 刘纪元, 庄大方, 2012. 国家尺度土地利用/覆被变化遥感监测方法[J]. 安徽农业科学, 40(4): 2365-2369. doi: 10.3969/j.issn.0517-6611.2012.04.161 [15] 杨振海, 2012. 基于Fragstats软件和TM影像的海南天然草地景观格局分析[J]. 安徽农业科学, 40(11): 6664,6675. [16] 邬建国, 2007. 景观生态学: 格局、过程、尺度与等级 [M]. 2版. 北京: 高等教育出版社. [17] 魏帆, 郭广轩, 张金萍, 等, 2018. 1985—2015年围填海活动影响下的环渤海湿地演变特征[J]. 生态学杂志, 37(5): 1527-1537. [18] 薛春燕, 2019. 基于森林资源调查数据的城市森林景观格局动态变化研究[J]. 陕西林业科技, 47(3): 36-42. doi: 10.3969/j.issn.1001-2117.2019.03.010 [19] 阳文锐, 2015. 北京城市景观格局时空变化及驱动力[J]. 生态学报, 35(13): 4357-4366. [20] 张林艳, 夏既胜, 叶万辉, 2008. 景观格局分析指数选取刍论[J]. 云南地理环境研究, 20(5): 38-43. doi: 10.3969/j.issn.1001-7852.2008.05.008 [21] 张莹莹, 蔡晓斌, 杨超, 等, 2019. 1974—2017 年洪湖湿地自然保护区景观格局演变及驱动力分析[J]. 湖泊科学, 31(1): 171-182. doi: 10.18307/2019.0116 [22] 赵同谦, 欧阳志云, 贾良清, 等, 2004. 中国草地生态系统服务功能间接价值评价[J]. 生态学报, 24(6): 1101-1110. doi: 10.3321/j.issn:1000-0933.2004.06.002 [23] 中华人民共和国中央人民政府, (2019-05-23)[2020-10-14]. 中共中央国务院关于建立国土空间规划体系并监督实施的若干意见[EB/OL]. http://www.gov.cn/zhengce/2019-05/23/content_5394187.htm. [24] 周华锋, 马克明, 傅伯杰, 1999. 人类活动对北京东灵山地区景观格局影响分析[J]. 自然资源学报, 14(2): 117-122. doi: 10.3321/j.issn:1000-3037.1999.02.004 [25] Costanza R, d’Arge R, de Groot R, et al, 1997. The value of the world’s ecosystem services and natural capital[J]. Nature, 387(6330): 253-260. [26] Enayat A, Fateme Z, Patrick M R, et al, 2020. Structuring and evaluating decision support processes to enhance the robustness of complex human–natural systems[J]. Environmental Modelling and Software, 123: 104551. doi: 10.1016/j.envsoft.2019.104551 [27] Liu J Y, Kuang W H, Zhang Z X, et al, 2014. Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s[J]. Journal of Geographical Sciences, 24(2): 195-210. doi: 10.1007/s11442-014-1082-6 [28] Mooney H A, Duraiappah A, Larigauderie A, 2013. Evolution of natural and social interactions in global change research programs[J]. PNAS, 110(Suppl 1): 3665-3672. [29] Ramalho C E, Laliberté E, Poot P, et al, 2016. Complex effects of fragmentation on remnant woodland plant communities of a rapidly urbanizing biodiversity hotspot[J]. Ecology, 95(9): 2466-2478. [30] Tian H, Xu H, Hall C A S, 1995. Pattern and change of a boreal forest landscape in northeastern China[J]. Water, Air, and Soil Pollution, 82(1/2): 465-476. doi: 10.1007/BF01182856 [31] Tian H Q, Chen G, Zhang C, et al, 2012. Century-scale response of ecosystem carbon storage to multifactorial global change in the Southern United States[J]. Ecosystems, 15(4): 674-694. doi: 10.1007/s10021-012-9539-x -