邵亚婷, 卢 霞, 叶 慧, 张 森, 孙 敏, 林雅丽, 赵 倩

(江苏海洋大学, 江苏 连云港 222005)

滨海湿地是介于陆地和海洋生态系统之间复杂的过渡性区域[1-2], 具有稳定环境、保护物种基因等功能, 在全球和区域生态过程中发挥着极为重要的作用[3]。盐城滨海湿地属于典型的近海及海岸湿地, 是太平洋西海岸、亚洲大陆边缘最大的海岸型湿地, 拥有独特的淤泥质海岸带及丰富多样的滩涂湿地生态系统, 已列入世界重点湿地名录[4]。江苏盐城湿地珍禽国家级自然保护区(简称盐城保护区)是盐城滨海湿地中保护最完整的部分, 是我国最大的沿海湿地自然保护区, 是全球重要的湿地之一。保护区内河流众多, 沼泽湿地众多, 蕴藏着丰富的生物资源, 为鸟类的栖息和其他动物的生存提供了得天独厚的自然条件, 是各种动物的繁殖地与鸟类的重要越冬地[5]。同时, 盐城保护区作为生产力最高的天然生态系统, 还为其周边居民提供着各种宝贵资源。然而, 在人口和经济的压力之下, 保护区内的实验区已经大部分被开垦为鱼塘或农田, 缓冲区也很大程度上受到人为干扰, 目前只有核心区仍然保留着大片原生湿地; 同时, 过热的经济开发活动也在不断导致滩涂湿地的资源遭到破坏。

土地覆盖是指自然营造物和人工建筑物所覆盖的地表诸要素的综合体, 包括地表植被、土壤、湖泊、沼泽湿地及各种建筑物(如道路等), 具有特定的时间和空间属性, 其形态和状态可在多种时空尺度上变化[6]。几十年来, 盐城保护区的土地覆盖发生了巨大的变化[7], 土地覆盖变化能够体现人与自然的相互作用, 人类对地表景观的开发、利用而引起土地覆盖发生变化。因此,及时有效地研究盐城保护区土地覆盖的变化对保护江苏盐城滨海湿地生态系统的安全与稳定, 维持人与自然的和谐共存非常重要。利用遥感技术对盐城滨海湿地土地覆盖变化的研究已经成为国内外研究的热点[7-18]。左平等[7]基于1976年的地形测绘数据和1987、1992、1997、2002、2007年的遥感数据提取土地利用和景观的时空变化, 并分析了江苏盐城滨海湿地的景观变化和驱动力。贺秋华[8]以江苏滨海地区为研究对象, 以1987年利用现状图、1997年和2006年TM遥感影像为数据源, 运用3S技术, 采用定性与定量相结合的手段, 系统分析了该区域土地利用时空变化特征。张华兵等[9]将盐城国家级自然保护区核心区划分为人工管理区和自然湿地区, 根据1987、1997、2007年3个时相的景观资料, 运用遥感(RS)、地理信息系统(GIS)技术和景观生态学方法, 分析不同驱动力下湿地景观格局的变化差异。王凯等[10]以盐城国家级珍禽自然保护区核心区为例, 基于3期(2002、2010、2013年)遥感影像数据, 通过监督分类方法研究了近10年来核心区湿地土地利用类型的变化情况, 为保护区核心区进一步加强管理与保护提供科学依据。郝敬锋等[11]基于RS和GIS技术, 利用1987—2007年期间3个时相的遥感图像提取了江苏省滨海湿地信息, 定量分析了研究区湿地资源的演变趋势和驱动机制。杜培军等[12]基于景观生态学原理、RS和GIS技术, 分析了1992—2009年江苏滨海湿地景观格局动态变化。丁晶晶等[15]选用1992年、2002年和2006年的TM影像, 分析了江苏省盐城海岸带湿地的景观格局特征及其景观动态变化。尽管以上学者利用RS和GIS相结合的方法开展了盐城滨海湿地土地覆盖变化研究, 但在盐城保护区土地利用及景观格局变化监测方面的研究稍显集中, 而对盐城保护区土地覆盖分类算法及其长时间序列的变化监测研究偏少。为深入探讨盐城保护区核心区的土地覆盖分类算法及其变化监测, 本文选用1973—2018年期间6期遥感影像, 通过构建土地覆盖分类规则集, 应用面向对象的分类方法提取研究区土地覆盖分类信息; 基于土地利用转移矩阵分析土地覆盖变化特征。通过该研究可为人类活动、盐城保护区的保护和建设提供有利的依据和建设性意见。

1 数据和方法

1.1 研究区概况

江苏盐城湿地珍禽国家级自然保护区位于江苏省盐城市境内的沿海地带(119°53′45″E~121°18′12″E, 32°48′47″N~34°29′28″N), 地跨响水县、滨海县、射阳县、亭湖区、大丰区、东台市的沿海滩涂, 其海岸线长达约500 km, 区域面积为247 260 hm2。盐城保护区属于北亚热带向南暖温带过渡的气候类型, 具有季风盛行、四季分明、日照充足、无霜期、降水丰富、雨热同期的特征。主要受海洋性和大陆性气候影响, 年平均气温为13.7~14.6 , ℃ 光照充足, 无霜期210~224 d, 雨水丰沛, 年平均降水量980~1 070 mm。

盐城保护区内的滩涂主要是1128—1855年, 黄河从黄土高原携带大量的泥沙以及来自海底的部分泥沙在潮流等海洋动力作用下淤积而成的粉沙淤泥质滨海平原。1855年黄河北迁, 泥沙来源骤减, 导致射阳河以北开始侵蚀, 射阳河以南仍然淤积, 盐城保护区海岸线总体上呈现南於北蚀的状态。

盐城保护区成立于1983年, 成立时面积为453 300 hm2, 核心区面积为17 400 hm2。1992年11月, 江苏盐城沿海湿地珍禽省级保护区经国务院批准升级为江苏盐城国家级珍禽自然保护区。2006年国务院对保护区进行了调整, 调整后保护区总面积和核心区面积分别为284 179 hm2、21 889 hm2, 保护区主要保护对象为丹顶鹤等珍稀野生动物及其赖以生存的滩涂湿地生态系统。2012年, 保护区范围再次进行了调整, 调整后保护区总面积247 260 hm2, 其中核心区面积22 596 hm2, 缓冲区56 742 hm2, 实验区167 922 hm2, 核心区、缓冲区和实验区等五块区域组成保护区。

选 取 1992年 保 护 区 的 核 心 区(120°29′2″E~ 120°42′44″E, 33°28′40″N~33°39′6″N)范围作为研究区域, 核心区位于新洋港和斗龙港之间的滩涂湿地, 以西至陈李线(盐城市响水县陈家港镇至南通市海安县李堡镇的省级公路), 以东至黄河水深3 m的区域。核心区所在的滩涂不断向海域淤积延伸, 目前拥有江苏最大的原生滩涂湿地, 研究区位置如图1所示。

1.2 数据获取与预处理

为开展盐城保护区核心区土地覆盖分类, 利用地理空间数据云平台选取1973年、1985年、2001年、2008年、2013年和2018年等6期美国Landsat卫星影像, 数据信息如表1所示。对遥感影像, 主要进行波段组合、图像配准、图像裁剪等预处理, 图像配准精度控制在0.5个像素之内。图像裁剪的矢量边界主要依据1992年11月经国务院批准的盐城保护区核心区的空间范围。

图1 江苏盐城湿地珍禽国家级自然保护区及保护区核心区位置 Fig. 1 Map of the core area of the Rare Bird National Nature Reserve in Yancheng Wetland, Jiangsu Province

1.3 构建研究区土地覆盖分类体系

根据第三次全国土地调查土地分类体系和中国滨海湿地分类系统[19], 结合盐城保护区土地利用覆盖现状和开发利用特点, 同时依据前人的研究成果[20-21], 通过分析盐城保护区核心区遥感影像, 对核心区土地覆盖的空间分布进行了研究, 确立了盐城保护区核心区土地覆盖分类系统(表2)。从表中可以看出: 研究区土地覆盖方式分为自然湿地和人工湿地两个一级类, 其中每个一级类又包含不同的二级类别, 最终将研究区划分为8个类型, 分别是海水、光滩、河流、米草群落、碱蓬群落、芦苇群落、茅草群落、养殖区。

1.4 研究方法

面向对象分类方法(object-oriented classification method), 是指将分割得到的同质影像对象作为最小单元, 以实现较高层次的遥感图像分类和目标地物提取[22]。基于像元的分类方法无法充分利用影像中存在的形状和纹理等特征, 可能会造成分类精度不高。而面向对象分类方法可以充分利用目标地物的光谱、形状和纹理信息, 提高了分类精度和分类效率。面向对象影像分类方法包括影像的分割和分类。

表1 研究区遥感影像数据信息 Tab. 1 Data of remote sensing images in research areas

表2 盐城保护区核心区土地覆盖分类体系 Tab. 2 Land cover classification system in the core area of Yancheng Reserve

1.4.1 影像分割

影像分割是面向对象分类方法的核心, 影像分割首先将单个像素合并为较小的影像对象, 再将较小的影像对象合并成较大的多边形对象, 是从下到上、逐级合并的过程[22]。多尺度分割后, 影像的最小单元不再是单个的像素, 而是具有相同光谱特征和共同属性的一个个对象, 后续的影像分析和处理也都是根据这些对象进行的[23-25]。分割参数的选择直接影响分割效果, 多尺度分割参数包含各波段权重、分割尺度、光谱因子和形状因子(紧密度和光滑度)。

1.4.2 影像分类

影像分割后进行分类, 分类方法包括基于样本的影像分类和基于规则集的影像分类。其中基于样本的分类与常规的光谱监督分类有显著不同, 其样本的对比参数更多, 不仅仅是光谱信息, 还包括空间、纹理等信息; 基于规则集的分类也是根据影像对象的不同特征建立规则集进行分类[26]。选取合适的影像对象分类特征是分类最重要的一部分, 分别统计计算分割后影像对象的光谱特征、纹理特征以及形状几何特征, 并通过建立规则集对分割后的影像进行分类。表3列出了本文选取的影像对象的分类特征。其中, 比值植被指数(RVI)强化了植被在近红外和红光波段反射率的差异, 在植被覆盖浓密的情况下对植被十分敏感, 非常适合高覆盖度的植被监测。

表3 影像对象的分类特征描述 Tab. 3 Description of classification features of image objects

2 结果与分析

2.1 研究区土地覆盖分类信息提取

利用 eCognition 软件对Landsat卫星影像进行基于规则集的面向对象分类, 首先对影像数据进行多尺度分割, 其次构建规则集进行分类。

2.1.1 确定最优多尺度分割参数

根据各时相遥感影像邻近像素亮度、纹理、颜色等对影像进行分割。通过反复试验和对比, 对每一时期的遥感影像进行多尺度分割, 最终得到各时相遥感影像的最优多尺度分割参数如表4所示。

表4 多尺度分割参数 Tab. 4 Multiscale segmentation parameters

2.1.2 基于规则集的分类

规则集构建和分类的主要步骤是:

1) 区分水体与非水体

根据盐城保护区核心区土地覆盖分类体系先将初始层定义为水体(海水、河流)与非水体(光滩、米草、碱蓬、芦苇、茅草、养殖区), 水体在可见光波段处的光谱特征比较明显, 不断试验确定水体的特征阈值, 构建水体的分类规则, 通过TS将海水提取出来, 河流的形状是细长的, 因此利用IShape将河流提取出来。

2) 区分植被与非植被

根据分类体系将未分类的非水体分为植被(米 草、碱蓬、芦苇、茅草)和非植被(光滩、养殖区), 根据光滩和养殖区的光谱特征和影像信息构建分类规则。利用NDVI和Brightness提取光滩, 养殖区中含水量大, 形状比较规则, 利用NIR、Blue和长/宽将养殖区提取出来。

3) 区分植被类型

研究区中植被主要包括米草、碱蓬、芦苇和茅草。在eCognition中创建NDVI、NDMI、RVI等可以区分植被的光谱特征, 根据不同植被在不同波段的光谱特征差异, 经过反复比较并通过多次试验确定了米草、碱蓬、芦苇、茅草的分类阈值。1973、1985、2001、2008、2013和2018年盐城保护区核心区的土地覆盖分类规则集如图2—图7所示。

图2 1973年遥感影像分类规则集 Fig. 2 Classification rules set of remote sensing images in 1973

图3 1985年遥感影像分类规则集 Fig. 3 Classification rules set of remote sensing images in 1985

4) 分类及后处理

基于规则集进行土地覆盖分类, 并进行合并地类、去除小图斑、平滑等一系列后处理, 利用 eCognition软件随机产生500个分布随机且均匀的样本点, 并根据Google Earth地图对样本进行分类, 建立混淆矩阵, 得到最终分类精度评价结果如表5所示, 盐城保护区核心区土地覆盖分类结果(图8)。

图4 2001年遥感影像分类规则集 Fig. 4 Classification rules set of remote sensing images in 2001

图5 2008年遥感影像分类规则集 Fig. 5 Classification rules set of remote sensing images in 2008

图6 2013年遥感影像分类规则集 Fig. 6 Classification rules set of remote sensing images in 2013

图7 2018年遥感影像分类规则集 Fig. 7 Classification rules set of remote sensing images in 2013

由图8和表6可知: 核心区在2001年之前基本为自然湿地, 基本未受人为干扰。1973年, 盐城保护区核心区土地覆盖由海到陆依次是海水、光滩、河流、碱蓬、茅草、芦苇; 其中光滩覆盖面积最大, 占总面积的47.90%。1985年, 盐城保护区核心区土地覆盖由海到陆依次是海水、光滩、河流、米草、碱蓬、茅草、芦苇; 其中碱蓬覆盖面积最大, 占总面积的26.97%, 新出现的米草群落面积为455.99 hm2, 占总面积的2.04%。2001年, 盐城保护区核心区土地覆盖由海到陆依次是海水、光滩、河流、米草、碱蓬、芦苇、茅草、养殖区; 其中碱蓬覆盖面积最大, 占总面积的23.26%。2008年, 盐城保护区核心区土地覆盖由海到陆依次是海水、光滩、河流、米草、碱蓬、芦苇、养殖区; 其中芦苇的覆盖面积最大, 占总面积的19.65%。2013年, 盐城保护区核心区土地覆盖由海到陆依次是海水、光滩、河流、米草、碱蓬、芦苇、养殖区; 其中米草群落的覆盖面积最大, 占总面积的17.84%。2018年, 盐城保护区核心区土地覆盖由海到陆依次是海水、光滩、河流、米草、碱蓬、芦苇、养殖区; 其中米草群落覆盖面积最大, 占总面积的20.00%; 养殖区面积为4 425.32 hm2, 占总面积的19.84%。

1973—2018年期间, 盐城保护区核心区土地覆盖变化趋势: 海水呈现先增长、后下降、再增长、再下降的趋势; 光滩呈现先下降后增长的趋势; 河流呈现平缓的趋势; 米草于20世纪60年代引入中国[4], 从遥感影像上看出, 1973年以后开始出现米草, 呈现从无到有, 持续增长的趋势; 碱蓬呈现增长、下降、再增长、再下降的趋势; 芦苇呈现先增长后下降的趋势; 茅草呈现先增长后下降的趋势; 养殖区从1985年以后开始出现, 呈现从无到有, 先增长、后平缓下降、再增长的趋势。

表5 分类精度评价结果 Tab. 5 Classification accuracy

图8 1973—2018年盐城保护区核心区土地覆盖分类结果 Fig. 8 Classification of land covers in the core areas of the Yancheng Nature Reserve from 1973 to 2018

表6 盐城保护区核心区各年土地覆盖类型面积 Tab. 6 Areas of land cover types in the core areas of the Yancheng Nature Reserve in each year

2.2 研究区土地覆盖变化监测

2.2.1 土地转移矩阵

土地利用转移矩阵是马尔科夫模型在土地利用变化方面的应用, 土地利用转移矩阵可以定量地表明不同土地利用类型之间的转化情况。土地利用转移矩阵来源于系统分析中对系统状态与状态转移的定量描述。

2.2.2 土地覆盖类型变化及转移特征分析

利用 ArcGIS 软件对eCognition导出的土地覆盖分类矢量数据进行融合和叠置分析, 得到1973—1985、1985—2001、2001—2008、2008—2013、2013— 2018年的土地覆盖变化(图9), 并计算得到转移矩阵(表7—表11), 以定量分析盐城保护区核心区土地覆盖动态变化特征。

在1973—1985年期间, 土地覆盖转变方式主要为光滩向海水、碱蓬、茅草的转变, 碱蓬向茅草的转变, 茅草向芦苇的转变。其中光滩向海水、碱蓬、茅草的转变主要分布在东部、中部、北部区域; 碱蓬向茅草的转变主要分布在西南部区域; 茅草向芦苇的转变主要分布在西北部区域。海水覆盖面积增加了1 133.59 hm2光滩覆盖面积减少了4 437.45 hm2, 河流覆盖面积增加了44.71 hm2, 碱蓬覆盖面积增加了1 540.27 hm2, 茅草覆盖面积增加了777.91 hm2, 芦苇覆盖面积增加了484.98 hm2, 米草覆盖面积增加了455.98 hm2。其中, 光滩转变了5 122.59 hm2, 转变率为47.96%, 有1 194.77 hm2(11.18%)转变为海水, 487.65 hm2(4.57%)转变为河流, 1 982.64 hm2(18.56%)转变为碱蓬, 128.85 hm2(1.21%)转变为芦苇, 880.96 hm2(8.25%)转变为茅草, 447.73 hm2(4.19%)转变为米草。

在1985—2001年期间, 土地覆盖转变方式主要为光滩向米草的转变, 碱蓬向芦苇、茅草的转变, 茅草向芦苇、养殖区的转变。其中, 光滩向米草的转变主要分布在中东部区域; 碱蓬向芦苇、茅草的转变主要分布在西部区域; 茅草向芦苇、养殖区的转变主要分布在北部、西南部区域。海水覆盖面积减少了238.87 hm2, 光滩覆盖面积减少了2 315.18 hm2, 河流覆盖面积减少了197.07 hm2, 米草覆盖面积增加了1 972.85 hm2, 增加了4.33倍, 碱蓬覆盖面积减少了786.72 hm2, 茅草覆盖面积减少了1 985.21 hm2, 芦苇覆盖面积增加了1 127.31 hm2, 养殖区覆盖面积增加了2 395.87 hm2。其中, 碱蓬转变了2 138.52 hm2, 转变率为35.80%, 有83.29 hm2(1.39%)转变为河流, 842.09 hm2(14.10%)转变为芦苇, 807.77 hm2(13.52%)转变为茅草, 74.04 hm2(1.24%)转变为米草, 331.32 hm2(5.55%)转变为养殖区。

在2001—2008年期间, 土地覆盖转变方式主要为光滩向米草的转变, 碱蓬向芦苇、米草的转变, 芦苇向养殖区的转变, 茅草向芦苇、养殖区的转变。其中, 光滩向米草的转变主要分布在东部靠海区域; 碱蓬向芦苇、米草的转变主要分布在中西部、中东部区域; 芦苇向养殖区的转变主要分布在西南、西北部区域; 茅草向芦苇、养殖区的转变主要分布在西北部和西部区域。海水覆盖面积减少了377.14 hm2, 光滩覆盖面积减少了255.88 hm2, 河流覆盖面积减少了108.79 hm2, 米草覆盖面积增加了1 151.39 hm2, 碱蓬覆盖面积减少了1 790.46 hm2, 茅草覆盖面积减少了1 249.47 hm2, 芦苇覆盖面积增加了1 721.78 hm2, 养殖区覆盖面积增加了908.56 hm2。其中, 碱蓬转变了2 105.26 hm2, 转变率为40.58%, 有117.07 hm2(2.26%)转变为河流, 1 548.65 hm2(29.85%)转变为芦苇, 302.51 hm2(5.83%)转变为米草, 137.03 hm2(2.64%)转变为养殖区。

图9 1973—2018年核心区土地覆盖变化图 Fig. 9 Maps of land cover changes of the core areas from 1973 to 2018

表7 1973—1985年核心区土地覆盖类型面积转移矩阵 Tab. 7 Transfer matrix of the core areas of land cover types from 1973 to 1985

表8 1985—2001年核心区土地覆盖类型面积转移矩阵 Tab. 8 Transfer matrix of the core areas of land cover types from 1985 to 2001

表9 2001—2008年核心区土地覆盖类型面积转移矩阵 Tab. 9 Transfer matrix of the core areas of land cover types from 2001 to 2008

表10 2008—2013年核心区土地覆盖类型面积转移矩阵 Tab. 10 Transfer matrix of the core areas of land cover types from 2008 to 2013

表11 2013—2018年核心区土地覆盖类型面积转移矩阵 Tab. 11 Transfer matrix of the core areas of land cover types from 2013 to 2018

在2008—2013年期间, 土地覆盖转变方式主要为光滩向海水、米草的转变, 米草向碱蓬的转变, 碱蓬向芦苇、米草的转变, 芦苇向碱蓬、养殖区、米草的转变。其中光滩向海水、米草的转变主要分布在东部靠海区域; 米草向碱蓬的转变主要分布在东部区域; 碱蓬向芦苇、米草的转变主要分布在中部、东北部区域; 芦苇向碱蓬、养殖区、米草的转变主要分布在中部、西北部、东北部区域。海水覆盖面积增加了945.59 hm2, 光滩覆盖面积减少了1 148.49 hm2, 河流覆盖面积增加了75.29 hm2, 米草覆盖面积增加了398.68 hm2, 碱蓬覆盖面积增加了433.98 hm2, 芦苇覆盖面积减少了582.68 hm2, 养殖区覆盖面积减少了122.37 hm2。其中, 芦苇转变了1 668.38 hm2, 转变率为38.08%, 有54.48 hm2(1.24%)转变为河流, 796.71 hm2(18.18%)转变为碱蓬, 180.92 hm2(4.13%)转变为米草, 636.27 hm2(14.53%)转变为养殖区。

在2013—2018年期间, 土地覆盖转变方式主要为海水向光滩的转变, 碱蓬向芦苇、米草的转变, 芦苇向碱蓬、养殖区的转变。其中, 海水向光滩的转变主要分布在东部靠海区域; 碱蓬向芦苇、米草的转变主要分布在西部和中部区域; 芦苇向碱蓬、养殖区的转变主要分布在中部、西北部区域。海水覆盖面积减少了1 195.84 hm2, 光滩覆盖面积增加了1 331.03 hm2, 河流覆盖面积减少了39.98 hm2, 米草覆盖面积增加了480.81 hm2, 碱蓬覆盖面积减少了1 172.70 hm2, 芦苇覆盖面积减少了646.57 hm2, 养殖区覆盖面积增加了1 243.26 hm2。其中, 米草转变了394.16 hm2, 转变率为9.91%, 有179.51 hm2(4.51%)转变为光滩, 45.09 hm2(1.13%)转变为河流, 120.98 hm2(3.05%)转变为碱蓬, 1.26 hm2(0.03%)转变为芦苇, 47.33 hm2(1.19%)转变为养殖区。

盐城保护区核心区土地覆盖在1973—2018年期间变化显著, 其中光滩面积减少最大, 依次是茅草、碱蓬, 而米草、养殖区、芦苇面积显著增加。在1973— 2013年, 光滩持续减少8 156.99 hm2; 而在2013— 2018年, 增加了1 331.03 hm2, 整体上光滩减少了6 825.96 hm2。在1973—1985年, 茅草增加了777.91 hm2, 而从1985年开始持续减少, 直至2008年全部消失, 整体上茅草减少了2 429.77 hm2。在1973—1985年, 碱蓬增加了1 540.27 hm2, 而在1985—2018年持续减少, 整体上碱蓬减少了1 775.64 hm2。在1973—2018年, 米草增加了4 459.71 hm2。在1985—2018年, 养殖区增加了4 425.32 hm2。在1973—2008年, 芦苇持续增加了3 334.08 hm2, 而在2008—2018年持续减少了1 229.26 hm2, 整体上芦苇增加了2 104.82 hm2。

3 结论与讨论

1) 基于长时间序列的遥感影像数据采用面向对象方法对盐城保护区核心区土地覆盖类型进行分类, 并利用土地利用转移矩阵得到土地覆盖变化特征信息。研究表明, 1973—2018年, 盐城保护区核心区土地覆盖变化总体上表现为光滩大面积减少, 转变方式是光滩→碱蓬群落→米草群落; 茅草群落完全消失, 转变方式是茅草群落→芦苇群落→养殖区; 碱蓬群落明显减少, 集中在研究区中部地区, 西部靠近海堤的主要演替是碱蓬群落→芦苇群落, 东部靠近滩涂的主要演替是碱蓬群落→米草群落; 养殖区大面积增加, 集中在研究区西南部地区; 芦苇群落明显增加, 集中在研究区西部和西北部地区。

2) 1973—2013年, 植物群落的扩张导致光滩面积的减少; 2013—2018年, 由于海岸滩涂的淤蚀, 光滩面积增加, 海水面积减少。米草群落和芦苇群落的扩张, 导致碱蓬群落减少; 芦苇群落和人工养殖区的扩张, 导致茅草群落消失。因此, 盐城保护区核心区土地覆盖类型变化的主要原因包括自然因素和人为干扰两个方面。自然因素包括海岸滩涂的侵蚀和淤积以及滨海湿地的自然演替; 人为干扰包括为开发滨海地区资源的相关政策和措施相继出台和实施以及引种外来物种米草常年不断的快速扩张。

针对1973—2018年对盐城保护区核心区土地覆盖变化特征的研究表明: 米草群落的持续扩张导致盐城保护区核心区的先锋碱蓬植被群落的退化, 应实施科学有效的管理, 防止米草过度繁殖, 恢复碱蓬群落。人工养殖区覆盖面积的持续增加给生态系统带来了一定压力, 应采取措施防止人工湿地(养殖区) 不断向自然湿地方向延伸, 尽量减少人为干扰对核心区的影响, 维持核心区湿地生态系统的安全与稳定, 实现盐城滨海湿地的可持续发展, 达到人与自然和谐共存的目的。