滨海生态脆弱区土地景观格局动态变化分析

刘洪柱+毛晓曦+王树涛

摘要：为探讨滨海生态脆弱区土地利用景观的时空动态变化过程，以河北省黄骅市为研究区，在1990、2000、2014年3期LandSat遥感影像的基础上，利用动态度模型、平面重心模型以及Arcgis 9.3、Fragststs 3.3等软件分析了该地区1990—2014年间土地景观格局的动态变化特征。结果表明，1990—2014年年间，从景观类型变化上看，建设用地和盐田面积增多，耕地、草地、湿地和盐碱地面积减少；园地和水域面积整体上较为稳定。从变化速度上看，盐田的变化速度最大，其次是湿地。从景观类型分布重心上看，各类景观都有不同程度的迁移，其中建设用地迁移的距离最远，其次是草地和盐碱地。从景观格局特征上看，景观整体破碎化程度增加，异质性增高。从景观格局梯度上看，水平样带上，一定范围内距离中心城区越远，土地景观类型越复杂，聚合度越差；但距海岸线30 km范围内，景观多样性降低，聚合度较大幅度上升；垂直样带上，由南到北，Shannon多样性指数整体呈上升趋势，中心城区向南，聚合度逐渐增加，向北则呈现出聚合度先小幅上升后又下降的趋势。研究结果可为该地区制定土地利用与生态环境建设规划，实现可持续发展与生态城市建设提供理论依据。

关键词：景观格局；动态变化；黄骅市；滨海生态脆弱区

中图分类号： F323.211文献标志码：A

文章编号：1002-1302（2017）05-0245-07

景观格局一般指景观的空间结构特征，是由一系列大小、形状各异的斑块，在自然因素及人为因素的驱动下互相作用、互相干扰，形成的具有空间异质性的景观区域[1-3]。土地景观格局的形成，不仅反映了区域内生态系统的演变，更反映出人与自然之间相互作用、相互依存的关系，景观格局指数能够高度浓缩景观格局信息，反映其结构组成和空间配置某些方面特征的简单定量指标[4-5]。景观格局分析是景观生态学研究的重要方法，研究土地景观格局的整体特征，分析其局部演变规律，有助于帮助人们认识土地资源变化过程，从而为土地利用规划的制定提供依据[6-8]。滨海位于海陆过渡带，其独特的地形地貌、海岸带、滩涂、湿地和野生动植物群落等自然景观是沿海特色的标志，也是保障滨海城市未来生活品质的基础。但自然因子急剧的梯度变化和脉冲式冲积使滨海生态系统处于脆弱状态，加之人类活动的强大干扰，滨海地区人口、资源、环境问题越发严重，土地景观格局也发生了较为明显的变化。因此，探讨滨海生态脆弱区土地景观格局的动态变化规律，可以为其土地生态治理提供理论支持。

本研究以滨海生态脆弱区河北省黄骅市为例展开土地景观格局动态变化的研究，首先研究了土地景观类型面积变化特征、变化速度以及空间重心的迁移，然后选择指标从类型水平和景观水平上阐述了该地区的土地景观格局变化特征，最后从水平和垂直2个梯度上对土地景观格局的动态变化进行了分析。本研究对当地制定土地利用与生态环境建设规划、实现可持续发展与生态城市建设具有重要的现实意义和指导意义。

1资料与方法

1.1研究区概况

黄骅市位于河北省东南部，是河北省沧州市管辖的县级市，东临渤海、西依沧州、南靠山东、北倚京津。位于 117.08°～117.82°E和38.15°～38.65°N之间，总面积为 1 717.50 km2，处在华北冲积平原黑龙港流域的最东端，地势低洼平坦，海拔高度一般为3～5 m[9]。地处暖温带半湿润季风气候区，因为靠近渤海而略具有海洋气候特征，季风显著，四季分明。历年平均日照时数为2 755 h，年平均蒸发量为 1 908.7 mm，年平均降水量为627 mm，降水年际变化大，时空分布不均。处于“环渤海”及“环京津”的枢纽地帶，不仅是沧州渤海新区的重要组成部分，也是“现代化沿海经济隆起带”的重要窗口。

1.2数据来源与处理

1.2.1数据来源

遥感数据来自于美国地质调查局网站和我国地理空间数据云。非遥感数据包括2014年黄骅市土地利用现状调查图、1 ∶[KG-*3]50万土壤类型图、DEM高程图、行政区划图，黄骅市国民经济和社会发展统计年鉴、黄骅市水利志等。

1.2.2数据处理

本研究使用ENVI 4.7软件进行遥感影像的预处理，预处理包括几何校正和影像裁剪。几何校正的目的是纠正系统及非系统性因素对遥感影像造成的影响，达到与标准图像或地图相吻合的目的[10-12]。影像裁剪的目的是将研究之外的区域去除。在进行了上述2步数据预处理操作后，采用监督分类法对1990、2000、2014年3期的影像图进行解译，最终得到黄骅市土地景观类型图

2结果与分析

2.1土地景观类型时空变化分析

2.1.1土地景观类型转移矩阵分析

为了探求研究区不同类型景观的数量转入转出情况，分别叠加1990年与2000年、2000年与2014年遥感影像解译图，得到转移面积矩阵，结果见表1、表2。横向表示了不同景观类型的转出面积，纵向表示了不同景观类型的补给来源。

2.1.1.1黄骅市1990年和2000年景观类型转移矩阵分析

2.2土地景观格局特征变化分析

考虑到黄骅市景观格局空间分布的实际情况，本研究选取的景观指数包括平均斑块面积（AREA_MN）、斑块所占景观面积的比例（PLAND）、斑块个数（NP）、斑块密度（PD）、景观形状指数（LSI）、丛聚指数（CLUMPY）、聚合度（AI）、蔓延度（CONTAG）、斑块内聚力指数（COHESION）、Shannon多样性指数（SHDI）、Shannon均匀度指数（SHEI）[20-21]。

2.2.1类型水平变化特征

在Fragstats 3.3软件支持下，计算得出各景观类型在斑块类型水平上的指数变化结果（表6）。从表6可以看出，耕地为研究区的基质景观，斑块所占景观面积的比例最大，1990年、2000年和2014年分别为 36.84%、36.42%、35.85%，斑块密度逐年增加，景观形状指数先增后减，表明耕地斑块的空间破碎度增强，但形状逐步趋向规则；园地、草地斑块所占景观面积的比例与平均斑块面积均较小，其中园地斑块密度逐年降低，丛聚指数和斑块内聚力指数均随时间的推移而升高，这说明园地景观破碎度逐渐小，且斑块愈来愈聚集，连通度提升，呈集中趋势发展；草地的景观形状指数增加幅度较大，由初期的63.83%增加到末期的99.41%，这说明草地斑块形状变异较大；建设用地斑块所占景观的面积比例逐渐增加，斑块密度逐渐增大，平均斑块面积增加明显，由8.45 km2增加到17.43 km2，表明其分布比较分散，破碎化程度较高，空间异质性程度也大；湿地的斑块密度、景观形状指数逐渐增加，平均斑块面积逐年减少，说明湿地面积有所减少，景观破碎度增加，且斑块形状愈来愈不规则；盐碱地面积减少较为明显，其面积占全部景观的比例由1990年的22.60%减少到2014年的21.59%，反映出盐碱地综合治理取得了一定的效果，但破碎度逐渐增加；盐田景观面积逐年增加，破碎化程度加重；水域景观在类型水平上变化不大，基本保持初期的状态。

2.2.2景观水平变化特征

经Fragstats 3.3景观格局指数软件计算得出景观水平指数的结果（表7）以及各指标的空间分布图。从表7可以看出，25年间，黄骅市景观斑块数量（NP）、斑块密度（PD）、景观形状指数（LSI）随年份的变化均呈现增长趋势，说明景观整体破碎化程度逐渐加重；蔓延度（CONTAG）均小于50%，表明研究区景观的破碎程度大于集聚程度，也就是说研究区景观的连通性较差；并且，该指数随年份变化呈现减少的趋势，也同样说明了景观趋于破碎化的状态。

1990—2014年，Shannon多样性指数（SHDI）逐年增大，分别为1.722 6、1.725 4、1.786 6，其中2000—2014年间年变化显著，说明土地利用类型越来越丰富，破碎化程度越来越高，其不定性的信息含量也越来越大，景观格局的异质性增高；从Shannon均匀度指数（SHEI）变化来看，SHEI值逐年增大，并且越来越趋近于1，这表明各景观类型分布趋于均匀，某种或某几种景观在研究区的主导地位有所弱化，景观优势度降低，各景观类型在空间上的分布出现均衡化。聚集度（AI）则呈现逐年减少的趋势，其变化说明景观的聚集程度逐年减弱，破碎化加重。[FL）]

2.3土地景观格局梯度变化分析

景观的空间梯度，指的是沿着某一特定的方向景观呈现出逐渐变化的、有规律的空间特征。本研究选择了关联度较小的2个指数——Shannon多样性指数（SHDI）和聚合度（AI）指数。

本研究选择了水平和垂直2条样带TR1和TR2。TR1基本与研究区内的河流水系平行，TR2垂直于TR1，2条样带均穿过中心城区，梯度变化更具代表性和典型性。样带宽为5 km，以中心城区的几何中心为原点，向左、向南为负值；在Fragstats 3.3中设置移动窗口通过“5 km×5 km”的正方形从中心点开始移动，最后输出结果。

Shannon多样性指数（SHDI）的梯度变化如图6所示。从水平方向上看，随着时间的推移，Shannon多样性指数整体上逐年增加；由中心城区向西，SHDI值随距离的增加不断增大，向东则先升高后降低。这表明水平方向上，中心城区连片的建设用地斑块占主导地位，导致Shannon多样性指数较低；一定范围内距离中心城区越远，土地景观类型就趋向复杂多样；但距海岸线30 km范围内又开始下降，这主要由于成片盐田的存在使得景观多样性降低。从垂直方向来看，由南到北，SHDI整体呈现上升趋势，分析其原因，主要是南部土地景观

类型较为单一，主要为耕地；北部广泛地分布着草地、园地、盐碱地等景观类型，导致景观的复杂化和丰富化。从时间变化上观察，发现中心城区以北的土地景观多样性逐年增加，中心城区向南，Shannon多样性指数则随年份的变化先升高后降低。

聚合度（AI）的梯度变化如图7所示。从水平方向上看，随着时间的推移，聚合度整体上逐年降低；中心城区向西，AI值随距离的增加不断降低，向东则先降低后升高，与SHDI的变化趋势大致相反。這表明在水平方向，一定程度上距离中心城区越远，土地景观类型的聚合度就越差，越趋向于离散分布；而海岸线附近由于大片的盐田而出现较大幅度的上升。从垂直方向来看，中心城区向南聚合度逐渐增加，时间上呈现先降低后升高的规律；中心城区向北，聚合度先增加后降低，并且年份上逐年下降。这主要是由于南部大片耕地在土地整理、盐碱地治理的政策下日趋连片、聚集分布；中心城区以北则受到人类活动的干扰较剧烈，各类型景观随时间推移日趋破碎，5 km范围内景观类型较为均一，聚合度小幅上升，而 5～25 km范围内草地、盐碱地、农村居民点显著增多，景观类型较为复杂，聚合度下降。

3结论与讨论

1990—2014年间，建设用地和盐田面积明显增多，耕地、草地、湿地和盐碱地面积明显减少，其中湿地是研究期间面积减少最多的地类；园地和水域面积整体上较为稳定，变化不大。从变化速度上看，1990—2000年盐田的变化速度最大，其次是湿地；从各景观类型的分布重心来看，1990—2000年和2000—2014年这2个研究时段内，各类景观都有不同程度的迁移。其中建设用地在1990—2000年间迁移的距离最远，其次是草地和盐碱地。

从景观类型水平特征变化来看，耕地斑块形状在2000—2014年间逐步趋向规则；草地景观形状指数增加较为明显，形状变异较大；建设用地扩张严重，斑块密度逐渐增大，景观破碎度增加；湿地的斑块密度、景观形状指数有所增加，平均斑块面积减少，说明湿地景观破碎度增加，斑块形状愈来愈不规则；盐碱地面积减少较为明显，占全部景观的面积明显减少，但破碎度逐渐增加；盐田景观面积逐年增加，破碎化程度加重；水域景观斑块类型水平变化较小，基本保持原来水平。从景观总体水平特征变化来看，黄骅市景观斑块数量、斑块密度、景观形状指数随年份增加均呈现增长趋势，表明景观整体破碎化程度增加，异质性增高。

从景观梯度观察发现，水平样带上，一定范围内距离中心城区越远，土地景观类型就越趋向复杂多样，其聚合度就越差；但距海岸线30 km范围内，由于成片盐田的存在，景观多样性降低，聚合度则出现较大幅度的上升。垂直样带上，由于南部主要为耕地景观，北部土地受人类活动的干扰较剧烈，并且广泛分布着草地、园地、盐碱地等景观类型，导致由南到北，Shannon多样性指数整体呈上升趋势；中心城区向南，聚合度逐渐增加，向北则因景观类型由均一到复杂而呈现出聚合度先小幅上升后又下降的趋势。

参考文献：

[1]赵阳，余新晓，贾剑波，等. 红门川流域土地利用景观动态演变及驱动力分析[J]. 农业工程学报，2013，29（9）：239-248.

[2]段翰晨，王涛，薛娴，等. 科尔沁沙地沙漠化时空演变及其景观格局——以内蒙古自治区奈曼旗为例[J]. 地理学报，2012，67（7）：917-928.

[3]李文杰，乌铁红，李晓佳，等. 内蒙古希拉穆仁草原旅游地景观格局动态变化[J]. 地理科学，2013，33（3）：307-313.

[4]刘吉平，赵丹丹，田学智，等. 1954—2010年三江平原土地利用景观格局动态变化及驱动力[J]. 生态学报，2014，34（12）：3234-3244.

[5]Wu J G，Hobbs R. Key issues and research priorities in landscape ecology：an idiosyncratic synthesis[J]. Landscape Ecology，2002，17（4）：355-365.

[6]Geist H J，Lambin E F. Proximate causes and underlying driving forces of tropical deforestation[J]. Bioscience，2002，52（2）：143-150.

[7]Moser D，Zechmeister H G，Plutzar C，et al. Landscape patch shape complexity as an effective measure for plant species richness in rural landscapes[J]. Landscape Ecology，2002，17（7）：657-669.

[8]Costanza R，DArge R，de Groot R，et al. The value of the worlds ecosystem services and natural capital[J]. Nature，1998，25（1）：3-15.

[9]李楠. 黄骅市水土资源空间优化耦合研究[D]. 石家庄：河北师范大学，2012.[ZK）]

[10]刘红玉，吕宪国，张世奎. 三江平原流域湿地景观多样性及其50年变化研究[J]. 生态学报，2004，24（7）：1472-1479.

[11]曹雪，羅平，李满春. 基于扩展CA模型的土地利用变化时空模拟研究：以深圳市为例[J]. 资源科学，2011，33（1）：127-133.

[12]龙腾飞，焦伟利，王威. 基于面特征的遥感图像几何校正模型[J]. 测绘学报，2013，42（4）：540-545.

[13][JP2]刘纪远，匡文慧，张增祥，等. 20世纪80年代末以来中国土地利[JP3]用变化的基本特征与空间格局[J]. 地理学报，2014，69（1）：3-14.[JP]

[14]马士彬，张勇荣，安裕伦. 山区城市土地利用动态空间分布特征——以贵州省六盘水市为例[J]. 自然资源学报，2012，27（3）：489-496.

[15]常成，刘霞，张光灿，等. 利用马尔柯夫过程预测蒙阴县土地利用/覆被格局变化[J]. 土壤，2010，42（2）：309-313.

[16]李保杰，顾和和，纪亚洲. 矿区土地复垦景观格局变化和生态效应[J]. 农业工程学报，2012，28（3）：251-256.

[17]徐丽华，王欢欢，张结存，等. 近15年来杭州市土地利用结构的时空演变[J]. 经济地理，2014，34（7）：135-142.

[18]亢庆，张增祥，王长有，等. 艾比湖绿洲农业区土地利用动态与盐碱化影响的遥感应用研究[J]. 农业工程学报，2006，22（2）：73-78.

[19]赵锐锋，陈亚宁，李卫红，等. 塔里木河干流区土地覆被变化与景观格局分析[J]. 地理学报，2009，64（1）：95-106.

[20]林娜，徐涵秋，何慧. 南方红壤水土流失区土地利用动态变化——以长汀河田盆地区为例[J]. 生态学报，2013，33（10）：2983-2991.

[21]崔晓伟，张磊，朱亮，等. 三峡库区开县蓄水前后景观格局变化特征[J]. 农业工程学报，2012，28（4）：227-234.