-
土地作为人类进行各类开发建设活动的载体,对社会的发展起着举足轻重的作用.土地覆被变化是人类活动与自然环境相互作用最直接的表现形式,是全球气候和环境变化研究关注的重要内容[1].土地利用变化引起土地覆被变化,由于人类开发建设活动的干预,在自然、经济和政治等各种因素的作用下,土地利用不断改变着土地覆被变化[2],土地覆被的变化反过来又对土地利用产生重要影响[3],近年来,土地利用/覆被的变化及驱动力的研究受到专家学者的一致关注[4-9].土地利用/覆被变化可以形象地展示人类作用于自然环境的过程,景观格局变化是土地利用变化在景观生态学研究中的表达,区域景观格局的研究是揭示该区域土地利用/覆被变化以及生态状况和空间变异特征的有效手段[10].景观格局变化是各类景观要素在一定时空尺度内发生变化而引起景观空间结构的变化[11-12].人类在土地上进行的各类社会经济活动,特别是在城市化进程加快、建设用地需求持续增加的社会背景下的活动,必将对区域的景观格局产生巨大影响,因此有必要对景观格局的动态变化进行分析.关于土地利用/覆被和景观格局变化的研究成果较多,但研究区大多选择东部平原区或沿海沿江的城市群,对于西部山区城市虽有研究,但针对重庆市两江新区土地利用/覆被和景观格局变化研究成果较少,所用评价指数也较单一,本研究共采用动态度、破碎度等17项指数,以城镇化快速推进的国家级开发开放新区——两江新区为例进行研究,增加了分析的科学性和可靠性,以期引导决策部门制定合理的土地利用规划和生态保护政策,缓解经济发展与生态文明建设之间的矛盾,保障区域土地节约集约利用,促进两江新区经济、社会和生态协调可持续发展.
A GIS-Based Research of the Changes of Land Use/Land Cover and Landscape Pattern——A Case Study of Liangjiang New Area of Chongqing
-
摘要: 为了解国家级经济开发区土地利用现状和存在的问题,协调建设用地扩张和生态保护之间的关系,以重庆市两江新区为例,基于2009年、2016年的土地利用变更调查数据库,借助ArcGIS 10.2平台,共选取17项指数对两江新区的土地利用/覆被和景观格局进行分析.结果表明:2009年-2016年,建设用地面积上升,其占土地总面积比上升15.24%,农用地与未利用地面积下降且逐步向建设用地转化;两江新区景观多样性和均匀度均有所下降,空间异质性程度逐渐降低;耕地、园地、未利用地和交通水利用地空间分布不断地趋向离散化和破碎化,其他农用地和城乡建设用地破碎度大且不断下降;农用地和未利用地的分维数增加,空间分布趋于不规则,建设用地分维数降低,斑块形状趋于单一化和规则化.Abstract: In order to understand the present situation of land use in national economic development zones and the existing problems in it and to coordinate the relationship between the expansion of construction land and the protection of cultivated land, a case study was made of Liangjiang New Area in Chongqing, in which 17 descriptive indexes were selected to analyze the land use/land cover and landscape pattern of the area based on the database of land use change between 2009 and 2016 and using ArcGIS 10.2 technology as a platform. The results showed that from 2009 to 2016 the area of construction land increased, with a growth of 15.24% in the ratio of construction land area to total land area, and the area of agricultural land and unused land decreased and they were gradually converted into construction land. Both diversity and evenness of the landscape of Liangjiang New Area decreased, and its spatial heterogeneity also decreased. The spatial distribution of cultivated land, garden land, unused land and traffic and water conservancy land tended to be more discretized and fragmentized; and the fragmentation of other agricultural land and the urban and rural building land were big, and both landscapes decreased gradually. The fractal dimension of the agricultural land and unused land increased, and their spatial distribution tended to be irregularized, while that of the construction land decreased and the shape of its patches tended to be simplified and regularized.
-
Key words:
- land use/cover /
- landscape pattern /
- Liangjiang New Area /
- Chongqing .
-
表 1 土地利用/覆被变化特征描述指数
指数 表达式 表达意义 土地利用数量变化 单一土地利用动态度(D) $D=\left(\frac{U_{b}-U_{a}}{U_{a}}\right) \times \frac{1}{T} \times 100 \%$ 表示某一地类面积在单位时间内的变化情况,值越大,表明变化越剧烈. 综合土地利用动态度(F) $F=\left[\sum\limits_{i=1}^{n} J_{a b} / \sum\limits_{i=1}^{n} U_{a}\right] \times \frac{1}{T} \times 100 \%$ 反映研究区整体土地利用变化的速度,能够综合反映研究区土地利用年变化情况. 土地利用变化强度(Q) $Q=\left(\frac{U_{b}-U_{a}}{U}\right) \times \frac{1}{T} \times 100 \%$ 反映土地利用变化过程中,人类活动对自然环境的影响情况和各地类面积的变化情况. 土地利用开发度(K) $K=\frac{X_{a b}}{U_{a}} \times \frac{1}{T} \times 100 \%$ 反映某地类面积增加的情况,以单位时间内某地类新开发面积占研究期初面积的比值来表达. 土地利用耗减度(H) $H=\left(\frac{J_{a b}}{U_{a}}\right) \times \frac{1}{T} \times 100 \%$ 反映某地类面积减少的情况,以单位时间内某地类新损耗面积占研究期初面积比值来表达. 土地利用程度变化 土地利用程度综合指数(Z) $Z=100 \times \sum\limits_{i=1}^{n} A_{i} \times C_{i}$ 定量反映研究区土地利用程度的变化情况和集约程度. 土地利用变化率( $\Delta_{L_{b-a}}$ )$\begin{array}{l} {\Delta _{{L_{b - a}}}} = {L_b} - {L_a}\\ {\Delta _{{{{L}}_{b - a}}}} = 100 \times \left[{\sum\limits_{i = 1}^n {{{\rm{A}}_i}} \times {{\rm{C}}_{ib}} - \sum\limits_{i = 1}^n {{{\rm{A}}_i}} \times {{\rm{C}}_{ia}}} \right] \end{array}$ 反映区域土地利用的整体水平和变化方向.若值大于0,则表示该区域处于发展期,反之处于衰退期. 土地利用区域差异 土地利用相对变化率(R) $R=\left[\frac{K_{b}}{K_{a}}\right] /\left[\frac{C_{b}}{C_{a}}\right]$ 它能很好地反映各区域土地利用变化的不同,若值大于1,表示该区域变化高于整个研究区. 注:U为研究区总面积,Ua为研究起始时某地类的面积;Ub为研究期末某地类的面积;T为研究时段长;Xab是指研究期内某地类转变为其它地类的总面积;Jab表示研究期内由其它地类转变为该地类的总面积;Ai为研究区第i级土地利用程度分级指数;Ci为研究区第i级土地利用程度分级面积百分比;n为土地利用程度分级数;Cia,Cib分别为某区域第i级土地利用程度在研究期始末面积百分比;Ka,Kb分别代表研究区某区域某地类研究期始末的面积;Ca,Cb分别代表整个研究区某地类研究期始末的面积. 表 2 景观格局特征选取指数及其生态意义
指数 表达式 生态意义 斑块个数(NP) NP=ni 用来描述景观的空间分布差异,它与破碎程度呈同向变化关系,值越大,破碎程度越深. 斑块面积(PA) PA=ai 斑块面积是景观格局最基本的空间特征 平均斑块面积(MPA) MPA=ai 用景观的斑块面积与斑块个数的比值表示,可以反映景观的构成和多样性. 多样性(H) $H=-\sum\limits_{i=1}^{m}\left[P_{i} \times {ln} \left(P_{i}\right)\right]$ 表示各斑块类型的数量和其在空间分布上是否均匀,也描述斑块的复杂程度、多样性和种类丰富与否的状况. 均匀度(E) E=H/Hmax 表示不同景观类型在其个数或面积方面是否均匀 优势度(D) D=ln(m)-H 表示某个或某几个斑块类型在整个景观中的支配作用,值越大,表示受某个或几个斑块类型支配越明显. 分离度(S) $\begin{array}{c}{\mathrm{S}=\mathrm{B} / \mathrm{A}, \mathrm{A}=\mathrm{P} \mathrm{A}_{\mathrm{i}} / \mathrm{TA}} \\ {\mathrm{B}=0.5 \sqrt{n_{i} / \mathrm{TA}}}\end{array}$ 衡量斑块类型在空间分布上是否分散,值越大,表示此景观类型分布越分散. 破碎度(C) C=(ni)/MPA 表示整个景观或某一景观类型的破碎度,值越大,则景观空间分异程度越深,它能体现景观对人类开发建设活动的反映程度. 分维数(F) F=2ln(0.25Li)/lnai 表示景观中斑块形状的复杂程度,值越大,说明该斑块的差异性越大,形态越不均匀. 注:ni表中为i类景观类型的斑块个数;ai为i类景观类型的斑块总面积;Pi为各景观类型在总景观中所占的比例;m为景观类型数;Li为i类景观类型的斑块周长;Hmax=ln(m)为景观最大多样性;TA为景观类型总面积. 表 3 2009年-2016年两江新区各土地利用类型面积变化情况
地类名称 2009年 2016年 2009-2016年差值 一级地类 二级地类 面积/hm2 百分比/% 面积/hm2 百分比/% 面积/hm2 百分比/% 农用地 耕地 44 741.83 38.16 30 635.67 26.13 -14 106.16 -12.03 园地 8 627.56 7.36 6 492.06 5.54 -2 135.5 -1.82 林地 23 548.33 20.08 22 615.89 19.29 -932.44 -0.80 其他农用地 2 013.23 1.72 1 496.56 1.28 -516.67 -0.44 汇总 78 930.95 64.32 61 240.18 52.24 -17 690.77 -15.09 建设用地 城乡建设用地 30 876.88 26.33 46 123.95 39.34 15 247.07 13.00 交通水利用地 2 447.26 2.09 5072.61 4.33 2 625.35 2.24 汇总 33 324.14 28.42 51 196.56 43.67 17 872.42 15.24 未利用地 未利用地 5 001.64 4.27 4 819.99 4.11 -181.65 -0.15 表 4 2009年-2016年两江新区土地利用格局转移概率矩阵
% 2009年 2016年 耕地 园地 林地 其他农用地 城乡建设用地 交通水利用地 未利用地 耕地 68.40 0.03 2.15 0.13 25.96 3.33 0.00 园地 0.03 75.08 0.01 0.04 21.06 3.77 0.01 林地 0.01 0.00 91.49 0.04 7.09 1.38 0.00 其他农用地 0.01 0.00 0.02 70.37 25.01 4.59 0.00 城乡建设用地 0.08 0.00 0.33 0.02 98.35 1.19 0.01 交通水利用地 0.00 0.00 0.22 0.00 0.85 98.92 0.01 未利用地 0.00 0.00 0.00 0.05 2.62 1.08 96.25 表 5 2009年-2016年两江新区土地利用/覆被描述指数
% 指数 地类名称 耕地 园地 林地 其他农用地 城乡建设用地 交通水利用地 未利用地 土地利用数量变化 单一土地利用动态度 -4.50 -3.54 -0.57 -3.67 7.05 15.33 -0.52 综合土地利用动态度 2.39 土地利用变化强度 -1.72 -0.26 -0.11 -0.06 1.86 0.32 -0.02 土地利用开发度 0.01 0.02 0.65 0.57 7.29 15.48 0.02 土地利用耗减度 4.51 3.56 1.22 4.23 0.24 0.15 0.54 土地利用程度变化 土地利用变化率 16.79 表 6 两江新区土地覆被类型相对变化率
% 行政区 土地覆被类型 耕地 园地 林地 其他农用地 城乡建设用地 交通水利用地 未利用地 北碚区 0.86 1.08 1.04 0.77 1.42 0.83 1.01 渝北区 1.08 1.03 1.00 1.14 0.95 1.07 0.99 江北区 0.66 0.81 0.99 0.62 0.95 0.87 1.01 表 7 2009年-2016年两江新区景观格局指数
景观类型 斑块个数/个 斑块面积/hm2 平均斑块面积/hm2 分离度 破碎度 分维数 2009年 2016年 2009年 2016年 2009年 2016年 2009年 2016年 2009年 2016年 2009年 2016年 耕地 5 326 8 788 44 741.83 30 635.67 8.40 3.49 0.17 0.27 6.34 25.21 2.77 2.84 园地 4 275 4 725 8 627.56 6 492.06 2.02 1.37 0.35 0.43 21.18 34.39 3.03 3.09 林地 9 116 8 996 23 548.33 22 615.89 2.58 2.51 0.31 0.32 35.29 35.78 2.84 2.85 其他农用地 5 099 4 052 2 013.23 1 496.56 0.39 0.37 0.80 0.79 129.15 109.71 3.28 3.35 城乡建设用地 24 841 20 833 30 876.88 46 123.95 1.24 2.21 0.45 0.34 199.85 114.10 2.79 2.69 交通水利用地 189 534 2 447.26 5 072.61 12.95 9.50 0.14 0.16 0.15 0.56 3.08 2.95 未利用地 1 108 1 124 5 001.64 4 819.99 4.51 4.29 0.23 0.24 2.46 2.62 2.97 2.98 -
[1] 刘纪远, 匡文慧, 张增祥, 等. 20世纪80年代末以来中国土地利用变化的基本特征与空间格局[J].地理学报, 2014, 69(1):3-14. doi: http://d.old.wanfangdata.com.cn/Periodical/dlxb201401001 [2] 陈佑启, 杨鹏.国际上土地利用/覆被变化研究的新进展[J].经济地理, 2001, 21(1):95-100. doi: 10.3969/j.issn.1000-8462.2001.01.022 [3] 李秀彬.全球环境变化研究的核心领域:土地利用/土地覆被变化的国际研究动向[J].地理学报, 1996, 51(6):553-558. doi: 10.3321/j.issn:0375-5444.1996.06.010 [4] 熊侣英, 师学义.黄土山丘区土地利用变化对生态系统服务价值的影响——以长河流域为例[J].水土保持研究, 2018, 25(2):335-340. doi: http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=stbcyj201802052 [5] 陈百明, 刘新卫, 杨红. LUCC研究的最新进展评述[J].地理科学进展, 2003(1):22-29. doi: 10.3969/j.issn.1007-6301.2003.01.003 [6] 摆万奇, 柏书琴.土地利用和覆盖变化在全球变化研究中的地位与作用[J].地域研究与开发, 1999, 18(4):13-16. doi: 10.3969/j.issn.1003-2363.1999.04.004 [7] 叶宝莹, 黄方, 刘湘南, 等.土地利用/覆被变化的驱动力模型研究——以嫩江中上游地区为例[J].东北师大学报(自然科学版), 2002, 34(1):100-104. doi: 10.3321/j.issn:1000-1832.2002.01.020 [8] 于兴修, 杨桂山.中国土地利用/土地覆被变化研究的现状与问题[J].地理科学进展, 2002, 21(1):51-57. doi: 10.3969/j.issn.1007-6301.2002.01.007 [9] 刘成武, 黄利民.土地利用/土地覆盖变化的研究方法[J].地域研究与开发, 2004, 23(4):11-14. doi: 10.3969/j.issn.1003-2363.2004.04.003 [10] 潘如玉, 申丽娟, 谢德体, 等.重庆旧城改造区土地覆被与景观格局变化——以北碚区为例[J].经济地理, 2015, 35(1):156-157. doi: http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jjdl201501022 [11] 岳德鹏, 王计平, 刘永兵, 等.京郊西北地区农地利用与景观格局时空特征[J].农业工程学报, 2008, 24(6):89-95. doi: 10.3321/j.issn:1002-6819.2008.06.018 [12] 冯应斌, 慕卫东.国家级新区建设用地扩展的景观格局特征分析——以重庆市两江新区为例[J].江苏农业科学, 2017, 45(8):224-228. doi: http://d.old.wanfangdata.com.cn/Periodical/jsnykx201708062 [13] 王秀兰, 包玉海.土地利用动态变化研究方法探讨[J].地理科学进展, 1999, 18(1):81-87. doi: 10.3969/j.issn.1007-6301.1999.01.012 [14] 黄庆旭, 何春阳, 史培军, 等.城市扩展多尺度驱动机制分析-以北京为例[J].经济地理, 2009, 29(5):714-721. doi: http://www.cnki.com.cn/Article/CJFDTotal-JJDL200905004.htm [15] 赵健, 魏成阶, 黄丽芳, 等.土地利用动态变化的研究方法及其在海南岛的应用[J].地理研究, 2001, 20(6):723-730. doi: 10.3321/j.issn:1000-0585.2001.06.010 [16] 朱会义, 李秀彬.关于区域土地利用变化指数模型方法的讨论[J].地理学报, 2003, 58(5):643-650. doi: 10.3321/j.issn:0375-5444.2003.05.001 [17] 杨俊, 廖和平, 庄伟, 等.基于惩罚型变权的低丘缓坡开发建设适宜性评价——以重庆市两江新区为例[J].西南大学学报(自然科学版), 2016, 38(8):113-119. doi: http://xbgjxt.swu.edu.cn/jsuns/jsuns/ch/reader/view_abstract.aspx?file_no=201608018&flag=1 [18] 李月臣, 宫鹏, 何春阳, 等.中国北方13省土地利用景观格局变化分析(1989-1999)[J].水土保持学报, 2005, 19(5):143-146. doi: 10.3321/j.issn:1009-2242.2005.05.036 [19] 黄俊芳, 王让会, 师庆东.基于RS与GIS的三工河流域生态景观格局分析[J].干旱区研究, 2004, 21(1):33-37. doi: http://d.old.wanfangdata.com.cn/Periodical/ghqyj200401007 [20] 吴莉, 侯西勇, 徐新良, 等.山东沿海地区土地利用和景观格局变化[J].农业工程学报, 2013, 29(5):207-216. doi: http://d.old.wanfangdata.com.cn/Periodical/nygcxb201305030