基于数字海图及遥感的近60年崇明东滩湿地演变分析

郑宗生，周云轩，田 波，姜晓轶，刘志国

(1．上海海洋大学信息学院，上海 201306;2．华东师范大学河口海岸学国家重点实验室，上海 200062;3．国家海洋信息中心，天津 300171;4．国家海洋局东海信息中心，上海 200137)

0 引言

湿地被称为地球之肾，是地球3大生态系统之一，是人类最重要的生存环境之一，特别是在沿海经济发达地区，湿地是重要的后备土地资源［1］。对湿地资源的开发和保护进行有效合理地规划，需掌握湿地的演变规律，但是湿地面积宽广、表面泥泞、变化频繁，给常规测量带来很大困难。20世纪60年代，遥感技术作为新兴的探测技术迅速兴起，其大尺度、快速、同步且高频动态观测等优势有助于实现研究区域生态环境和资源的宏观、动态、同步的监测，大大弥补了常规观测方法的不足。

自从美国Landsat卫星发射以来，已有大量国内外学者尝试利用遥感技术监测潮滩和湿地的动态变化。监测方法可分为定性判读法［2－4］、基准线定量分析法［5－9］和水边线定量分析法［10－14］等。上述方法为监测潮滩、湿地的动态变化提供了多种手段，但是湿地不仅包括卫星图像能够监测到的陆面部分，还包括被高浊水体覆盖的水下沟槽，并且湿地水下沟槽的变化与陆面部分的变化相关联。滩面的地貌塑造过程与异常天气直接相关，单从遥感图像上观测潮滩变化，往往忽略滩/槽冲淤交互(滩冲槽淤)的过程，因此湿地演化研究需将2部分作为整体统一考虑，通过2部分的对比验证才能进一步加深对湿地动态变化的认识。特别是近年长江口入海泥沙逐年减少，几十年后将下降到2亿t以下［15］，这是否会直接导致潮滩、湿地的侵蚀，还需要将湿地作为一个整体加以研究。卫星数据可快速获取湿地出露部分的动态变化，而数字海图可分析湿地水下的地形演变。鉴于此，本文整合数字海图数据及多时相遥感图像，利用水边线法及等深线叠加对比分析法综合分析上海市崇明东滩湿地近60 a来的演化，为指导海岸带土地利用、海岸工程建设及预测东滩湿地未来发育趋势提供依据。

1 研究区概况

崇明东滩湿地1992年被列入《中国保护湿地名录》，1998年经上海市人民政府批准建立鸟类自然保护区，2002年2月2日被正式列入“拉姆萨国际湿地保护公约”的国际重要湿地名录。崇明东滩湿地位于长江口崇明岛的东端，E 121°47'～122°05'，N31°25'～31°38'(图1)，由长江下泻泥沙沉积形成，目前其自然淤涨速度仍可达到每年200～300 m［16］。人类对湿地的重要影响主要是围垦及放牧活动，1990年和2001年，大堤的2次修建表明崇明东滩海岸线从1990年位置最多向外推进达7．2 km;同时，过度的人工放牧也直接或间接的破坏了原有的湿地生态环境［17－18］。

图1 崇明东滩位置图Fig．1 Location of Chongming Dongtan

2 研究方法

利用多时相遥感图像及现场测量，分别提取水边线和海堤信息，并对潮滩类型进行分类，选取相近潮位水边线进行对比分析，同时整合不同年份海图数据进行潮滩面积及等深线变化分析，从定性及定量角度综合分析崇明东滩的演化情况，技术路线如图2所示。

图2 技术路线图Fig．2 Technique flow chart

2．1 现场测量

2005年9月9—14日在崇明东滩设立了3条垂直岸线的样带，分别位于潮滩北部(断面Ⅰ)、中部(断面Ⅱ)和南部(断面Ⅲ)，如图1所示。每隔50m对断面进行高程测量，高程基准统一到吴淞高程基准。高程及地理坐标使用Ashtech Z_XtremeTM GPS Receiver System测定，其空间定位精度达到cm级，高程测量精度达mm级。样带上每隔100 m选取一个样方进行植被生态因子测量，样方选取单一优势植被群落，大小为1 m×1 m，每个样方记录植被类型、植被平均高度并估算植被盖度。

2．2 数字海图生成

对收集的14幅由交通部安全监督局、海事局测量的1958年12月—2003年5月海图进行扫描，其具体参数见表1。

表1 海图数据Tab．1 Nautical charts

将相近年份的海图资料进行整合处理。以经纬度交叉点作为控制点，对海图进行几何纠正，每幅图均匀选取10～15个控制点，以提高拟合效果［19－20］;将映射坐标的偏差作为误差并求其均方根，几何纠正误差控制在5 m以内。对纠正后的海图选取50个水深点重复数字化10次，将各水深点横纵坐标与数字化后横纵坐标之差的均方根作为数字化的人为误差。结果表明50个采样点的绝对误差小于3．112 m，满足工程分析的精度。对数字化后的海图数据进行拓扑关系检查和水深异常值排查，排除数字化过程中等深线交叉、水深点异常等的错误录入。上述操作通过ArcInfo 10．0软件的相应模块完成。

2．3 遥感数据处理

选取1979—2005年之间的10个时相的遥感图像(天气状况均为无云)作为数据源，为了分析潮滩的动态变化，需要对不同时相遥感图像(表2)进行几何纠正等预处理，使所有数据统一到相同的坐标系统之下。以上海市1∶50 000地形图为基准，选取道路交叉点、桥梁、水系交叉点作为控制点进行几何纠正;同时引入2004年采集的长江口区82个地面控制点作为几何纠正坐标控制点。所有遥感图像的配准误差均小于0．5个像素，上述工作主要利用Erdas 9．1遥感处理软件完成。

表2 卫星图像数据Tab．2 Satellite images data

分别对经过几何纠正的MSS图像的4(R)，3(G)，2(B)波段，TM 及 ETM+图像的5(R)，4(G)，3(B)波段进行合成。同时为了提高崇明东滩区图像目视解译的精度，对图像进行了分段线性拉伸等标准化处理，使卫星图像层次感更强，细微结构纹理更清晰。

2．4 水边线信息提取

水边线遥感信息提取方法一般分为单波段法和多波段法。常用的多波段法又可分为比值法［21］和谱间关系阈值法［22］。长江口区影响水边线信息提取的主要因素是其附近的高浊度悬沙及浅滩表层的残余水体［23］。通过对长江口区卫星图像的波段谱及现场资料分析，认为TM5，TM6及TM7是提取水边线的有效波段，但由于陆源热水的注入，影响了TM6波段的提取效果［24］。针对研究区潮滩光谱特点，本文首先按卫星图像的潮时进行分类，分为涨潮图像及落潮图像。由于在涨潮前湿地较长时间暴露，表层的残余水体较少，TM5波段可作为此时的主要分类波段;落潮时刻潮水退去，表层浅沟中残留大量水体，TM5波段对水体敏感的特点使水边线提取结果偏向陆地一侧，此时TM7波段可作为主要分类波段。MSS图像的分析结果表明，MSS2波段为提取水边线最佳波段，而MSS4波段的提取效果最差。

依据上述波段选取原则，水边线提取采用人工参与及计算机自动解译结合的方式完成。具体操作为:①将选取的最优波段图像进行均值平滑，消除噪声并剔除陆地区的细节信息，利用ISODATA非监督分类方法对图像进行分类;②将分类结果进行再分类，归并不同植被类型及不同浊度水体的分类结果;③对再分类结果进行二值化处理;④利用区域增长算法从二值化图像中提取矢量水陆边界，邻域系统采用4邻域，阈值根据多次的试验结果及实际水边线位置进行设定;⑤生成的结果转换为ArcInfo 10．0识别的．shp矢量格式，导入Geodatabase数据库。

3 结果分析

3．1 长期趋势分析

崇明东滩自2001年大提修建以后，基本保持稳定的冲淤态势，本文选取崇明东滩1979—2005年遥感图像，分析其历史演化趋势(图3)。

图3 1979—2005年崇明东滩湿地演化过程Fig．3 Evolution process of Chongming Dongtan wetland from 1979 to 2005

图3清楚地再现了东滩近26 a来淤涨与围垦交替转换的情况:图像上比较亮的绿色调为植被，而水体在该波段有高吸收作用，使高、中潮滩为暗绿色，低潮滩为更深的暗绿色;1979—1989年10 a间最大的变化发生在东滩东南部团结沙的并陆，团结沙的分流作用使白港水道径流及泥沙逐年减少，1979年白港潮汐汊道上口筑坝堵汊，使团结沙面积迅速增大;在1989年图像上可以见到团结沙已经和陆地连为一体，原来的白港汊道退化为潮沟系统;目前东滩东南部最大的潮沟系统应是在此潮沟系统的基础上发展而来，东滩在围垦过程中，大堤外潮沟系统上部基本是原来白港水道的走向，在中部转向东南，与潮汐入射方向一致。

从图3还可定性地分析出:崇明东滩“南冲北淤”现象明显，造成南侧围堤以外的潮滩宽度较北侧窄;1979—1989年是东滩演化最快的10 a，人类工程的建设是导致其变化的主要原因;目前长江北支与白港水道表现相同的态势，北支上部与南支主航道成近直角相交，分流分沙比的减少直接导致北支累积性淤涨。

图4将数字海图数据0 m线作为分析基线，整合东滩湿地白港水道封堵前(1951年、1958年)及封堵后(1991年、2000年)的水下地形变化情况。可以看出团结沙在封堵前是一块孤立的沙体且“南冲北淤”，沙体向北侧移动。沙体0 m线逐渐与水道的走向一致，沙体的西南—东北方向淤积加快。从1991年0 m线看出，上口的封堵使大量的泥沙在团结沙东北部落淤，而0 m线继续向北侧移动，出现了侵蚀情况，这与卫星图像的分析结果一致。主要原因是汊道封堵加大了沙体靠海一侧的沿岸流，受外海潮汐的共同影响，使沙体继续向北侧移动。

图4 1951—2000年0 m等深线变化图Fig．4 0m depth contour change from 1951 to 2000

结合ArcInfo软件统计分析功能分别对数字海图0 m线及大堤之间的面积进行分析，并对比崇明东滩湿地的围垦面积。1968—1992年间，东滩湿地0 m线上面积增加了近125 km2，年均增加5．1 km2;1992—1998年间0 m线上面积增加7 km2，年均增加下降到1．1 km2;1998—2000年增加16 km2，年均增加又恢复到5 km2左右。结合同期的湿地围垦数据:1968—1992年间围垦面积为64 km2，1992—1998 年间为 46 km2，1998—2000 年间为 6 km2。3个时间段围垦面积与0 m线上增加面积的比分别为0．51，7．46 及0．36，其中 1992—1998 年间围垦面积是0 m线上增加面积的7倍多。

卫星图像显示的水边线受成像瞬间潮位的影响，在图像上直接解译湿地的高、中、低潮滩存在困难，根据实地地貌测量与植被调查，本文按标志性植被来划分潮间带:高潮滩以芦苇分布下界为限，中潮滩处于芦苇分布下限与海三棱藨草分布下限之间，低潮滩为海三棱藨草分布至0 m线之间。以此标准对东滩湿地潮间带重新进行面积统计(表3)。

表3 1991—2005年崇明东滩湿地潮间带面积变化Tab．3 Intertidal area variation of Chongming Dongtan wetland from 1991 to 2005

表3进一步说明在东滩湿地自然演化过程当中，虽然湿地总面积在逐年增加，但是人类围垦使高、中、低潮滩的结构发生了改变。自然状态下高潮滩面积要大于中、低潮滩面积，而崇明东滩刚刚相反。高潮滩面积所占比例逐年下降，即使在围垦活动结束后的2002—2005年的3 a间这个比例也从8．8%下降到7．4%。从现场观测剖面来看，东滩的高潮滩主要分布在北部、中部的互花米草带内，在东南部高潮滩相对较少，且呈不连续分布。

3．2 近期演化分析

采用水边线方法，对相近潮位的图像进行分析，可定量计算水边线位置的水平沉积速率。本文以东滩周围横沙、佘山、九段东、中浚、北槽中等5个潮位站的信息为依据选取4组潮位相近的卫星图像进行叠加分析(图5)。图5定量表明崇明东滩在中低潮滩仍然表现为“南冲北淤”的趋势，这与上述定性分析结果相一致。其中淤积最快的是东滩候鸟保护区岸段，1989—2002年的13 a间2 m(理论基面)水边线向前推进的距离最大达到3．9 km，年均的水平淤积速率超过300 m/a;北六滧港附近形成舌状淤积体，最宽处达到1．6 km;与此相反，南部奚家港至团结沙岸段2m水边线则表现为侵蚀后退，其中后退距离最大达500 m，年均水平侵蚀速率为27．68 m/a。图5(c)的高潮滩水边线也表现了相同的沉积趋势，但是高潮滩的最大沉积方向更靠南侧。从图5(a)—(c)中可以看出，1990年后，崇明东滩的淤积主要集中东北部、中部及东南部区域，北部岸段的淤积明显变缓;南部由于护岸工程的影响，侵蚀作用亦趋于停止，成为相对稳定的岸段。由图5(d)还可以看出，2003—2005年间，崇明东滩候鸟保护区北部甚至出现了轻微的侵蚀，表明此岸段冲淤交替现象的存在。

图5 不同年份同潮位水边线对比分析Fig．5 Waterline comparison with the same tidal level in different years

崇明东滩两侧北港、北支水流在此撞击形成低流速区，致使水流能量降低;同时，由于东滩滩面平缓纵比降小，河床断面扩大，促使潮流分散，造成潮流流速减慢:因此大量泥沙在此落淤，滩地得到迅速淤涨。另外，由于泥沙颗粒表面的电化学絮凝和盐淡水异重流产生的垂向环流共同作用，底层泥沙含量增加，在东滩缓流区内聚积形成高含沙量区，加速了崇明东滩的淤涨;而南侧的侵蚀主要是受波浪作用的影响。在1990年后，崇明东滩的淤积主要集中在候鸟保护区岸段，北部岸段的淤积速率明显减缓;南部沿岸护岸工程外的潮滩面积已经很小，由于工程的影响，侵蚀作用亦趋于停止，岸段相对稳定。2003—2005年间，东滩北部甚至出现了轻微的侵蚀(图5(d))，说明此岸段冲淤交替现象的存在。这主要与异常的天气条件有关，如风暴潮，台风等，这些都可能造成滩面的侵蚀。

4 结论

利用遥感图像及数字海图综合分析湿地出露滩体及水下部分的冲淤演变过程，能够从整体上反映1951—2005年间崇明东滩的发育进程及形态变化，更好地再现了东滩的地貌过程及人类活动的影响。主要结论如下:

1)崇明东滩湿地目前的演化过程已趋于稳定，在长江口河势稳定的情况下，将保持“南冲北淤、部分交替、中部淤涨”的态势。

2)人类活动对崇明东滩的影响主要集中在高潮滩部分。东滩中低潮滩及水下地形的变化主要是受该地区水动力条件控制，但是人类庞大的封堵工程也会诱发潮滩空间格局的较大变化。

3)人类活动是崇明东滩湿地自然演变的重要驱动力。白港潮汐通道的封堵使东滩湿地成陆过程大大加快。东滩湿地0 m线较为稳定的向东部扩展，在人类逐年围垦活动的干预下，湿地面积虽然逐年增加，但潮间带结构已偏离湿地自然状态，高潮滩面积所占比重减少，且呈不连续分布状态。海堤直接暴露于外部海洋环境中，势必降低了对灾害天气的防御能力。

综上所述，在整个东滩湿地的演化过程中，人类既有起正作用的水道封堵使湿地面积增加，又有起负作用的围垦造陆使湿地面积减少。如长江三峡工程及上游水库的修建，使长江口入海泥沙减少到2003年的2亿t，不足20世纪60年代的1/2，这必将对长江口潮滩整体冲淤态势产生影响。但潮滩的冲淤演变过程受流域来水、来沙、水动力条件、异常天气及人类活动等多种因素综合影响，特别是崇明东滩区的演化过程还需重点考虑局部水动力条件及人类活动的影响。当前，无论是哪方面考虑，人类都在以自身的利益出发，谋求湿地资源的无偿馈赠，崇明东滩湿地区也不例外;而长江上游流域来沙的减少及全球海平面的上升也必然直接或间接地影响着长江口湿地资源的自然演化态势，因此，如何对崇明东滩湿地进行合理规划，开发与保护并重，当前及长远利益兼顾，是人类值得思考的重要课题。

［1］ 周 鑫，左 平，滕厚峰，等．基于土地利用变化的生态系统服务价值核算——以江苏盐城滨海湿地为例［J］．海洋通报，2011，30(6):656－661．Zhou X，Zuo P，Teng H F，et al．Estimation of ecosystem service values based on land use changes:A case study of Yancheng coastal wetlands［J］．Marine Science Bulletin，2011，30(6):656－661．

［2］ 马仲荃，严蔚芸，曹 豪．彩色红外航空摄影在海涂调查中的应用——以杭州湾上海岸段遥感为例［J］．海洋通报，1989，6(4):64－69．Ma ZQ，Yan W Y，Cao H．The application of color infrared aerial photography to the investigation of tidal wetland in the Shanghai section of the Hangzhou Bay［J］．Marine Science Bulletin，1989，6(4):64－69．

［3］ 王贵明，董裕国．珠江韩江三角洲海岸变迁遥感解译对比研究［J］．海洋科学，1997(4):50－53．Wang GM，Dong YG．Interpretation and study of remote sensing of coastal－line migration in the Pearl River and the Hanjiang River Delta［J］．Marine Sciences，1997(4):50－53．

［4］ Blodget H W，Taylor P T，Roark JH．Shoreline changes along the Rosetta－Nile promontory:Monitoring with satellite observations［J］．Marine Geology，1991，99(1/2):67－77．

［5］ 冯永玖，刘 丹，韩 震．遥感和GIS支持下的九段沙岸线提取及变迁研究［J］．国土资源遥感，2012，24(1):65－69．Feng Y J，Liu D，Han Z．Shoreline extraction and change analysis of the Jiuduansha islands with the support of remote sensing and GIS technologies［J］．Remote Sensing for Land and Resources，2012，24(1):65－69．

［6］ 宋启帆，王少军，张 志，等．基于World ViewⅡ图像的钨矿区水体信息提取方法研究［J］．国土资源遥感，2012，24(2):33－37．Song Q F，Wang S J，Zhang Z，et al．A water information extraction method based on World View Ⅱ remote sensing image in tungsten ore districts［J］．Remote Sensing for Land and Resources，2012，24(2):33－37．

［7］ 陈 蕾，邓孺孺，陈启东，等．基于水质类型的TM图像水体信息提取［J］．国土资源遥感，2012，24(1):90－94．Chen L，Deng R R，Chen Q D，et al．The extraction of water body information from TM imagery based on water quality types［J］．Remote Sensing for Land and Resources，2012，24(1):90－94．

［8］ 庞家珍，姜明星，李福林．黄河口径流、泥沙、海岸线变化及其发展趋势［J］．海洋湖沼通报，2000(4):1－6．Pang JZ，Jiang M X，Li F L．Changes and development trend of runoff，sediment discharge and coastline of the Yellow River estuary［J］．Transaction of Oceanology and Limnology，2000(4):1－6．

［9］ 黄海军，樊 辉．1976年黄河改道以来三角洲近岸区变化遥感监测［J］．海洋与湖沼，2004，35(4):306－314．Huang H J，Fan H．Monitoring changes of nearshore zones in the Huanghe(Yellow River)delta since 1976［J］．Oceanologia Et Limnologia Sinica，2004，35(4):306－314．

［10］ Mason D C，Davenport I，Flather R A．Interpolation of an intertidal digital elevation model from heighted shorelines:A case study in the western wash［J］．Estuarine，Coastal and Shelf Science，1997，45(5):599－612．

［11］ Mason D C，Davenport I，Flather R A，et al．A digital elevation model of the intertidal areas of thewash produced by thewaterline method［J］．International Journal of Remote Sensing，1998，19(8):1455－1460．

［12］ Mason D C，Amin M，Davenport I J，et al．Measurement of recent intertidal sediment transport in Morecambe Bay using the waterline method［J］．Estuarine，Coastal and Shelf Science，1999，49(3):427－456．

［13］ Chen L C，Rau JY．Detection of shoreline changes for tideland areas using multi－temporal satellite images［J］．International Journal of Remote Sensing，1998，19(17):3383－3397．

［14］ 韩 震，恽才兴．伶仃洋大铲湾潮滩冲淤遥感反演研究［J］．海洋学报，2003，25(5):58－64．Han Z，Yun C X．Deposition and erosion remote－sensing reverse of Dachan Bay beach in Lingdingyang Estuary［J］．Acta Oceanologica Sinica，2003，25(5):58－64．

［15］ 杨世伦，朱 骏，赵庆英．长江供沙量减少对水下三角洲发育影响的初步研究［J］．海洋学报，2003，25(5):83－91．Yang SL，Zhu J，Zhao QY．A preliminary study on the influence of Changjiang River sediment supply on subaqueous delta［J］．Acta Oceanologica Sinica，2003，25(5):83－91．

［16］ 刘清玉，戴雪荣，何小勤．崇明东滩沉积环境探讨［J］．海洋地质动态，2003，19(12):1－4．Liu Q Y，Dai X R，He X Q．The sedimentary environment of the east tidal flat of the Chongming island［J］．Marine Geology Letters，2003，19(12):1－4．

［17］ 赵云龙，安传光，林 凌，等．放牧对滩涂底栖动物的影响［J］．应用生态学报，2007，18(5):1086－1090．Zhao Y L，An C G，Lin L，et al．Effects of grazing on zoobenthos community in beach［J］．Chinese Journal of Applied Ecology，2007，18(5):1086－1090．

［18］ 高 宇，王 卿，何美梅，等．滩涂作业对上海崇明东滩自然保护区的影响评价［J］．生态学报，2007，27(9):3752－3760．Gao Y，Wang Q，He M M．To evaluate the impact of economic activities on the mudflats of Chongming Dongtan birds nature reserve，Shanghai［J］．Acta Ecologica Sinica，2007，27(9):3752－3760．

［19］ 赵书河，冯学智，赵 锐．中巴资源一号卫星南京幅数据质量与几何纠正评价［J］．遥感技术与应用，2000，15(3):170－174．Zhao SH，Feng X Z，Zhao R．Evaluation on data quality and geometric correction of China－Brazil resources satellite No．1 data in Nanjing area［J］．Remote Sensing Technology and Application，2000，15(3):170－174．

［20］ 韩 涛．在1:25万电子地图辅助下的TM影像精校正研究［J］．遥感技术与应用，2003，18(1):6－9．Han T．Research for geometrical correction of TM image based on 1 ∶250 000 digital map［J］．Remote Sensing Technology and Application，2003，18(1):6－9．

［21］ 徐涵秋．利用改进的归一化差异水体指数(MNDWI)提取水体信息的研究［J］．遥感学报，2005，9(5):589－595．Xu H Q．A study on information extraction of water body with the modified normalized difference water index(MNDWI)［J］．Journal of Remote Sensing，2005，9(5):589－595．

［22］ 杨存建，周成虎．基于知识发现的TM图像居民地自动提取研究［J］．遥感技术与应用，2001，16(1):1－6．Yang C J，Zhou CH．Extracting residential area from TM image on the basis of knowledge discovered［J］．Remote Sensing Technology and Application，2001，16(1):1－6．

［23］ 郑宗生，周云轩，蒋雪中，等．崇明东滩水边线信息提取与潮滩DEM 的建立［J］．遥感技术与应用，2007，22(1):35－38．Zheng Z S，Zhou Y X，Jiang X Z，et al．Waterline extraction and DEM reconstruction in Chongming Dongtan［J］．Remote Sensing Technology and Application，2007，22(1):35－38．

［24］ 郑宗生，周云轩，刘志国，等．基于水动力模型及遥感水边线方法的潮滩高程反演［J］．长江流域资源环境，2008，17(5):763－767．Zheng Z S，Zhou Y X，Liu ZG，et al．Dem reconstruction based on hydrodynamic model and waterline method［J］．Resources and Environment in the Yangtze Basin，2008，17(5):763－767．