滨海滩涂湿地是指大潮高潮位与低潮位之间的潮浸地带, 包括潮间带泥滩、沙滩和咸水沼泽等, 是滨海湿地的重要组成部分。滨海滩涂湿地资源丰富, 是全球最有价值的生态系统之一, 在食物生产、水源涵养、水质净化、气候调节、保持土壤、营养循环、生物多样性维持、支撑区域经济社会发展等方面发挥着不可替代的作用[1—3]。但由于中国沿海地区人口众多, 对滩涂湿地的掠夺性开发利用, 导致滩涂湿地面积减少和功能退化, 在一定程度上影响了区域可持续发展[4—6]。滨海滩涂湿地形成、发育和演变以及分布格局, 受到海陆相互作用以及人为影响的共同塑造, 尤其国家关于滨海湿地保护管理政策是滨海湿地面积变化趋势的重要驱动。定量分析滨海滩涂湿地的格局演变过程, 有助于为滨海滩涂的资源评估、合理利用和保护管理提供科学依据[4, 7—11]。

近年来, 利用遥感卫星影像进行专题地图制作已经较为成熟, 利用遥感技术可以实现滨海滩涂湿地变化研究[12—16]。通过基于对象的影像分析来进行滨海湿地分类研究较为广泛, Google Earth Engine(GEE)、云计算和大功率计算经常用于量化滨海地区土地利用/覆盖变化的时间序列研究[17—19]。然而, 大多数研究时间间隔较大, 分类不够细致, 未明确区分光滩和盐沼面积动态变化[20—21], 导致滨海湿地保护管理措施政策与滨海湿地变化难以做到实时对应, 难以满足对滨海湿地高质量管理的需求。本研究以30年长时间序列下的中国滨海滩涂湿地空间分布特征为核心, 结合“陆域—滩涂—海域”之间的相互作用机制以及快速城市化背景, 从自然因素和人为因素等方面对典型区域的滨海滩涂时空演变进行驱动力分析, 为中国滨海滩涂湿地保护与合理利用提供支持。

中国滨海滩涂湿地北至鸭绿江口, 南达广西北仑河口, 涉及辽宁省、河北省、天津市、山东省、江苏省、上海市、浙江省、福建省、广东省、广西壮族自治区, 海南省和台湾省的海岸潮间带[3]。综合考虑滨海滩涂地区的地形地貌和海岸类型, 以从全球高分辨率地理数据集(GSHHG)中获取的海岸线为基准, 向陆延伸10 km, 向海延伸40 km的区域作为本研究滨海滩涂湿地信息提取的范围。

本研究采用的遥感数据源为Landsat卫星影像。使用了1990—2020年逐年覆盖研究区域所有可用的L1T正射校正Landsat图像, 传感器类型包括搭载于Landsat 5卫星的专题制图仪TM(Thematic Mapper)、Landsat 7卫星的增强型专题制图仪ETM+(Enhanced Thematic Mapper plus)和Landsat 8卫星的陆地成像仪OLI(Operational Land Imager), 信息如表 1所示。Landsat卫星影像来源于美国地质调查局(USGS), 在Google Earth Engine平台可通过在线代码编辑器编写代码调用经大气校正的地表反射率(Surface Reflectance)影像。辅助数据则包括来源于全球行政区域划分数据集(GADM, https://gadm.org/)的行政区域划分数据、全球高分辨率地理数据集(GSHHG)的海岸线数据。此外, 还有研究区域的地面高程及水深数据, 数据来源于全球地形水深数据(ETOPO1)。

参照全国第二次湿地资源调查分类标准, 结合滨海湿地的实际情况, 研究将滨海滩涂湿地分类为光滩和盐沼(含红树林滩涂等)。所有Landsat图像预处理任务均在GEE平台上进行, 利用Landsat时间序列影像选取地面参考样本, 利用Google Earth Pro软件识别已利用地、水体、光滩和盐沼等四种地物类型, 将样本点范围以*.kml格式从Google Earth Pro中导出, 在ArcMap 10.7中转化为shp格式, 导入GEE API中进行调用[22—23]。将获取的Landsat地表反射率影像集使用CFMask算法对云、云阴影进行识别, 并利用QA波段进行去除, 在GEE API中调用随机森林分类器, 以土地利用类型为目标变量, 将Landsat影像的原有光学波段、光谱指数和ETOPO 1数据作为随机森林分类器的特征变量, 决策树的数量设置为100, 其他参数保持默认设置, 实行大范围滨海滩涂湿地空间分布信息提取[24—26]。研究假定5a内训练样本保持5%的变化, 在目标年的高分辨率影像上选取700个随机点(滩涂类别200个, 滨海植被类别100个, 其他类别400个), 并以此为5a的逐年随机森林分类算法的训练样本点, 利用谷歌地球高分辨率影像, 通过分层随机采样的方法对验证数据集和训练数据集进行精度验证, 利用ArcGIS 10.7软件对1990—2020年制图产品, 每年选取随机点导入Google Earth Pro软件, 进行目视判读, 验证并生成研究中滨海滩涂湿地分类结果。

中国滨海滩涂湿地面积分布特征自东向西总体呈先下降后增长的趋势, 自南向北面积分布趋势变化较大(图 1)。其中, 光滩主要集中在中国沿海东部地区(118°E—123°E), 在地势低的海岸线与泥沙流速高的地区[27—28], 如黄河三角洲、辽河三角洲、江苏沿海岸段等。大范围盐沼主要集中在三个经度区间：21—23°N, 珠江三角洲所在地；30—34°N, 长江三角洲以及江苏盐城滩涂湿地所在地；37—41°N, 辽河及黄河三角洲所在地。盐沼多与光滩连接, 是在光滩的基础上演替而来, 所以其分布与光滩的区域较为一致。长江三角洲、黄河三角洲、江苏盐城等沿海地段, 这些地区的盐沼面积共计约占中国的45.4%。

根据滩涂湿地空间分布数据(表 2)可知, 江苏省滩涂湿地面积最大, 2020年江苏滩涂湿地总面积为24.47×104hm2, 其次为浙江省14.94×104hm2、山东省10.55×104hm2。光滩面积较大的省份为江苏省、浙江省、山东省, 面积及占比分别为21.99×104hm2(31.07%)、13.98×104hm2(19.74%)、9.37×104hm2(13.24%)。盐沼面积较大的省份则为江苏省、上海市、山东省、浙江省, 面积及占比分别为2.48×104hm2(25.77%)、2.42×104hm2(25.16%)、1.17×104hm2(12.21%)、0.96×104hm2(10.01%)。

从近30年沿海滩涂湿地面积年份变化情况(图 2)上来看, 1990年至2020年, 中国滩涂湿地面积变化趋势为明显减少, 总面积由1990年140.99×104hm2减少到2020年的80.39×104hm2, 减少了42.98%。其中光滩面积减少显著, 减少61.05×104hm2, 减少了46.30%；盐沼面积变化总体较为平稳并略有增加, 增加了0.45×104hm2, 增加了4.95%。以10年为时间段分析变化趋势显示, 1990年至2000年, 滩涂湿地面积呈现持续减少趋势, 减少了20.56×104hm2(14.58%)；2000年至2010年, 滩涂湿地面积呈现先动态平稳后减少趋势, 减少了15.69×104hm2(13.02%)；2010年至2020年, 滩涂湿地面积呈现先持续减少后保持平稳的趋势, 减少了24.35×104hm2(23.25%)。

中国沿海各省份滩涂湿地(光滩、盐沼)面积年份变化差异明显(见图 3)。江苏省滩涂湿地总面积在沿海省份中保持最大, 但损失高达16.85×104hm2, 减少了40.79%。其他滩涂湿地总面积减少较大的省份有山东省、辽宁省, 与1990年相比分别减少了16.25×104hm2(60.65%)、9.28×104hm2(58.39%)。此外还有河北省和天津市, 与1990年相比减少了5.05×104hm2(70.02%)和1.43×104hm2(65.09%)。2020年光滩主要集中分布在江苏省、浙江省、山东省、辽宁省、福建省等沿海地区。近30年来, 江苏省光滩面积损失15.99×104hm2(42.09%)、浙江省为2.71×104hm2(16.26%)、山东省为16.32×104hm2(63.51%)、辽宁省8.70×104hm2(57.65%)、福建省4.98×104hm2(45.31%)。此外减少比例较大的省份还有河北省和天津市, 与1990年相比减少了4.96×104hm2(69.91%)和1.43×104hm2(65.39%)。沿海各省份2020年盐沼面积较1990年变化差异明显。面积减少的省份有江苏省、广东省、辽宁省、河北省、天津市, 减少面积及占比分别为0.87×104hm2(25.97%)、0.61×104hm2(43.34%)、0.58×104hm2(72.2%)、0.09×104hm2(77.43%)、0.01×103hm2(9.18%)。其余沿海省份盐沼面积均增加, 其中广西增幅最大, 达到687.99%(0.68×104hm2)。增幅较为明显的省份还有台湾556.31%(0.02×104hm2), 福建365.79%(0.36×104hm2), 浙江159.62%(0.59×104hm2)。山东、江苏、福建、广东、广西、海南分阶段滩涂湿地面积变化趋势与中国滩涂面积变化趋势呈现高度的一致性, 即先减少后动态平稳之后再减少最后动态平稳；辽宁、天津在最初1990年至2000年先出现面积增加后与中国变化一致；河北则为持续性减少；上海2000年至2010年滩涂面积为增加趋势；浙江滩涂面积在2016—2020年持续增加；台湾滩涂面积的变化趋势则为先减少后动态平稳再增加最后减少。

江苏盐城分布有中国最典型的滩涂湿地, 占江苏省滩涂总面积近70%, 是由古黄河三角洲和长江三角洲的泥沙在黄海和东海潮汐作用下形成的淤泥质滩涂[29—30]。1990年至2020年, 盐城滨海滩涂湿地面积变化趋势明显减少, 总面积由1990年的24.55×104hm2减少到2020年的16.24×104hm2, 其中滩涂面积减少了8.23×104hm2, 盐沼面积减少了0.08×104hm2。盐城滩涂湿地面积变化趋势与中国滩涂面积变化趋势一致, 近30年内以2010年为时间节点的前后两个阶段均呈现先减少后动态平稳的规律。该时段盐城滩涂湿地中的光滩面积呈下降趋势, 盐沼面积基本保持稳定状态, 其变化趋势如图 4所示。

中国滨海地区共分布有4大三角洲类型的滩涂湿地, 包括长江三角洲、黄河三角洲、辽河三角洲和珠江三角洲湿地[31—33]。1990年至2020年的四大三角洲滩涂湿地面积整体呈下降趋势, 光滩和盐沼面积均随时间明显减少, 如图 5所示。长江、黄河、辽河、珠江三角洲2020年对比1990年光滩面积分别减少了3.82×104hm2(46.42%)、14.73×104hm2(66.78%)、3.03×104hm2(49.50%)、2.15×104hm2(66.33%)。长江、黄河、辽河、珠江三角洲2020年对比1990年盐沼面积分别减少了0.87×104hm2(52.53%)、0.02×104hm2(2.17%)、0.53×104hm2(77.29%)、0.61×104hm2(68.52%)。

滨海滩涂湿地是泥沙在入海口或海岸淤积而成, 其变化受到泥沙沉积和和水动力环境的影响。大面积滩涂湿地一般位于沿海地段、地势低平、排水不畅、泥流速度快的区域, 例如黄河三角洲、长江三角洲等。河流水沙及其造陆运动为三角洲滩涂湿地发展提供了源动力, 河口冲淤是三角洲滩涂湿地面积不断增长的重要原因。但是, 受气候变化影响流域降水发生变化, 以及黄土高原森林覆盖率的增加, 黄河入海泥沙量呈阶梯状减少[34], 2020年利津水文站年输沙量为3.14×108t[35], 为多年均值的49%；同时黄河故道处海岸侵蚀加快, 新黄河口处冲淤能力减弱, 导致黄河三角洲2020年对比1990年滩涂湿地面积减少了14.75×104hm2(63.76%)。此外, 黄河三角洲近几年互花米草面积增加, 也是光滩面积减少和盐沼面积增加的重要原因[36—37]。长江三角洲地形演变受流域来水来沙影响较大, 2020年对比1990年滩涂湿地面积减少了4.69×104hm2(47.45%), 2000年以后, 大通水文站的年输沙量出现了断崖式下降, 2020年大通水文站年输沙量为1.64×108t[38], 为多年均值的47%, 长江口来沙量锐减, 发生由淤积主导向侵蚀主导的转变[39—41], 长江三角洲滩涂湿地面积减少。

国家政策是主导滨海滩涂湿地面积变化的主要原因。在2000年至2010年, 《中华人民共和国海域使用管理法》(2001年)、《国务院关于进一步加强海洋管理工作若干问题的通知》(2004年)、《国务院办公厅关于加强湿地保护管理的通知》(2004年)先后出台, 提出了海域的合理开发、可持续利用以及严控制围填海, 中国滩涂湿地面积快速减少的趋势得到短期遏制, 但仍然持续减少, 总计减少15.69×104hm2(13.02%)。2010年至2020年, 滩涂湿地总面积呈现先持续减少后保持平稳的趋势, 总计减少了24.35×104hm2(23.25%)。在此期间, 《国务院关于印发全国主体功能区规划的通知》(2010年)和《全国海洋经济发展“十二五”规划》(2012年)发布, 提出在中国优先发展沿海区域的战略背景下, 沿海基础建设以及滨海产业的发展导致沿海省份用地紧张, 沿海城市通过“围填海”方式缓解区域发展的用地需求, 因此, 2010年至2015年期间滩涂湿地面积急剧减少。在此之后, 《全国海洋主体功能区规划》(2015年)、《湿地保护修复制度方案》(2016年)和《国务院关于加强滨海湿地保护严格管控围填海的通知》(2018年)等一系列有利于滨海滩涂湿地保护的政策出台, 进一步加强滨海湿地保护, 严格管控围填海活动, 尤其是在2019年, “全面禁止围填海”写进十九届四中全会公报, 对违法违规项目进行生态修复, 对严重破坏海洋生态环境的坚决予以拆除, 中国滩涂湿地面积再次归于平稳趋势, 2015年—2020年间, 除个别省份略有减少外, 到2020年中国滩涂湿地面积已略有增加, 国务院有关文件精神得到有效落实, 《中华人民共和国湿地保护法》于2022年6月起实施, 为滨海湿地保护提供了法律保障。

基建占用、渔业发展、滩涂围垦、油田开发等因素是造成滨海滩涂湿地减少的直接原因[9, 42]。围垦养殖、基建占用和油田开发等直接引起滨海湿地面积减少的热点区域主要分布在江苏盐城、辽河三角洲、黄河三角洲、长江三角洲和珠江三角洲等地区。盐城是国家实施沿海开发战略的重要节点, “开发沿海滩涂、建设海上盐城”战略促使大面积的滩涂围垦在相当长一段时间内持续进行, 除用于种植业、水产养殖业外, 还被用于工业开发和港口扩建[29]。辽河三角洲、黄河三角洲地区石油资源丰富, 油田开发是导致其滩涂湿地面积减少的主要原因[33, 43]。长江三角洲地区作为中国经济发展最活跃、人口最密集的地区之一, 其城镇化和经济快速发展造成用地需求的不断增长, 其变化驱动力由围垦后进行农业生产转为城市扩张, 成为滩涂面积减少的最主要原因[27]。珠江三角洲地区是我国第二大经济区, 存在大量的浅海养殖活动, 此外, 为实现产业结构转型, 在深圳、广州和珠海等城市的滨海滩涂区域, 通过围海造地建立产业园区。

此外, 在滩涂围垦热点区主要分布在渤海湾、江苏省和杭州湾[44], 滩涂湿地广泛受到人类活动干扰[45]。对近10年围垦热点区域的分析, 滩涂湿地空间变化分布情况见图 6。河北省人工围填海面积达2.80×104hm2, 转化为人工湿地0.39×104hm2, 转化为滩涂湿地0.33×104hm2；天津市人工围填海面积达1.51×104hm2, 转化为人工湿地0.01×102hm2, 转化为滩涂湿地0.25×102hm2, 这两个省份对滩涂的人工开垦主要用于工业发展, 此外, 天津滨海新区建设是天津地区围海造地, 滩涂面积损失的重要原因。江苏省人工围填海面积达4.63×104hm2, 转化为人工湿地3.94×104hm2, 转化为滩涂湿地2.09×104hm2, 江苏省滩涂湿地大部分转换为养殖鱼塘和农业用地及港口、道路等基础建设, 人工开垦面积较大、利用方式较为均衡。上海市人工围填海面积达0.70×104hm2, 转化为人工湿地3.57×104hm2, 转化为滩涂湿地0.01×104hm2, 滩涂围垦主要转化为养殖塘等人工湿地。浙江省人工围填海面积达3.41×104hm2, 转化为人工湿地11.84×104hm2, 转化为滩涂湿地3.90×104hm2, 其对滩涂的人工开垦相较于工业发展更多的是转化为农业(水稻田)、渔业(养殖塘)发展。

近30年来, 中国滨海滩涂湿地面积不断减少, 从1990年的140.99×104hm2减少到2020年的80.39×104hm2, 其中光滩面积损失46.30%、盐沼面积增加4.95%。中国滩涂湿地时空变化主要受自然条件和人为因素两方面影响, 自然因素在大时间尺度上控制着滩涂湿地演化, 人为因素在短期内影响滩涂湿地的动态变化。不同省份滩涂湿地面积变化总体趋势与全国一致, 国家关于土地利用、经济发展、环境保护的一系列相关政策是影响滩涂围垦的重要因素, 各省份因发展需求和压力不同, 滩涂湿地面积减小的变化程度有一定差异。随着人们湿地保护意识的增强, 在全面禁止围填海的前提下, 通过严守生态保护红线、加强滨海滩涂湿地保护恢复、合理利用滨海滩涂湿地资源、明确各方责任职责等一系列政策措施, 将有效促进中国滨海滩涂湿地的可持续发展。