海南陵水沙坝潟湖海岸近26年地貌类型空间格局演化

第３６卷第１期第四纪研究Ｖ０１．３６，Ｎｏ．１２０１６年１月ＱＵＡＴＥＲＮＡＲＹＳＣＩＥＮＣＥＳＪａｎｕａｒｙ，２０１６ｄｏｉ：１０．１１９２８／ｊ．ｉｓｓｎ．１００１—７４１０．２０１６．０１．１０文章编号１００１－７４１０（２０１６）Ｏ卜１０３—１０海南陵水沙坝溻湖海岸近２６年地貌类型空间格局演化术贾培宏①②③徐伟①③李功成④⑤殷勇①②③李鑫④（①南京大学地理与海洋科学学院，南京２１００２３；②南京大学中国南海研究协同创新中心，南京２１００２３；③南京大学海岸与海岛开发教育部重点实验室，南京２１００２３；④海南国际旅游岛先行试验区，陵水５７２４００；⑤同济大学建筑与城市规划学院，上海２０００９２）摘要通过对海南陵水县黎安、新村两个沙坝渴湖海岸地貌的野外踏勘，构建研究区地貌类型分类标准，据此采用德国智能影像分析软件ｅＣｏｇｎｉｔｉｏｎ提供的随机森林法模块，对１９８７—２０１３年２６年间６期Ｌａｎｄｓａｔ卫星遥感影进行地貌类型分类与空间分布位置图斑提取，结合ＡｒｃＧＩＳ软件据实测地形数据生成数字高程模型以及地貌类型转移矩阵数理分析技术，对研究区地貌格局演化特征及地貌类型转化趋势展开研究。结果表明，上述野外测量、随机森林遥感影像分类提取、ＧＩＳ空间分析及矩阵转换数理统计等多技术方法结合用于沙坝涡湖海岸地貌长周期演化研究是可行、有效的，此方法弥补了传统采用钻孔及沉积柱状样采集分析中以点代面、难以全域覆盖的海岸带沉积地貌分析方法的不足。研究发现，近２６年间，陵水双溻湖海岸地貌发生了显著变化，地貌类型转化趋势为海积平原、水下浅滩、海积阶地向养殖池及盐田转化，自然地貌间的相互转化不明显；大范围的潮滩围垦及海积平原挖掘改造，形成了一条人工地貌带，打破了原有自然地貌空间格局，改变了海积平原与水下浅滩之间的沉积动力作用方式以及物质能量交换模式，不断加强的人类活动已成为研究区海岸地貌改造的主要驱动力。主题词沙坝渴湖海岸地貌空间格局遥感地貌解译地貌转移矩阵海南陵水中圈分类号Ｐ７３６．１４，Ｐ７３７．１，ＴＰ７９文献标识码Ａ１引言物圈多计划（ＩＧＢＰ）及全球环境变化人文因素计划（ＩＨＤＰ）中海岸带关键研究区的典型代表∽卅１｜。海岸带处于大气圈、岩石圈、水圈和生物圈的我国沙坝溻湖海岸主要分布于华北、华南地交互地带，是全球变化中极为敏感的部位，如今人区¨２｜，在海南岛沿岸分布广泛。新村与黎安双渴类大量高强度的开发活动都集中于此¨ｑ］。沙坝漓湖位于海南岛东南部陵水县境内，新村渴湖在解放湖海岸是沉积海岸的重要组成部分，主要由沙坝、前已开始有码头建设，１９６７年两个漏湖开始发展水溺湖和潮汐汉道等地貌单元构成，在世界海岸线上产养殖，８０年代初开始沿岸开发养殖池及盐田¨引，广泛分布，约占全世界大陆岸线总长的１３％［４］，在双湖及沿岸自然环境开始受到人类活动的持续影中国的砂质海岸中，约６０％的海岸线属此种类型。响。２００９年，海南国际旅游岛建设发展意见正式提渴湖具有蓄水防洪作用，其外侧沙坝能够阻挡波浪出，２０１２年１月，新村与黎安片区被划为海南国际及风暴潮，防止岸线受侵蚀冲刷，湖内平时风浪平旅游岛先行试验区，预计未来１０年后，建造一座滨静，营养物质含量丰富，适宜港口及养殖场开发。海新城。近百年来，人类潮滩围垦、盐田与养殖池开发以及本研究利用１９８７—２０１３年２６年间６期多时相航道疏浚、防波堤修建等强活动干扰，导致沙坝渴Ｌａｎｄｓａｔ卫星遥感影像、实测海陆地形图、历史海图湖海岸消亡数目增加¨“１。这一体系对人类活动的等多数据源，运用专业遥感影像处理软件Ｅｎｖｉ及德响应具有极强的敏感性和脆弱性，是国际地圈与生国智能影像分析软件ｅＣｏｇｎｉｔｉｏｎ的随机森林法模块，第一作者简介：贾培宏女４３岁讲师海洋信息科学及海岸带规划管理Ｅ－ｍａｉｌ：ｊｉａｐｈ＠ｎｊｕ．ｅｄｕ。ｃｎ＋海南省政府项目“海南国际旅游岛先行试验区潮汐汉道海湾动力地貌及双湖连通工程选址研究”（批准号：２０１３０５８４）资助２０１５—１０—０２收稿，２０１５－１０—２８收修改稿万方数据第四纪研究借助ＧＩＳ数字高程模型三维模拟分析方法，解译提取陵水沙坝双漓湖海岸不同时期沙坝溻湖海岸地貌类型及其时空分布特征；运用地貌类型转移矩阵数理统计方法定量分析沙坝溺湖海岸地貌类型时空格局演化规律。探索一种借助长周期遥感影像资料和随机森林遥感影像分类提取、ＧＩＳ空间分析及矩阵转换数据理统计等多种方法结合用于沙坝淌湖海岸地貌长周期演化研究的有效技术方法，以弥补传统采用钻孔及沉积柱状样采集分析中以点代面、难以全域覆盖的海岸带沉积地貌演化分析方法的不足。根据研究结果，进一步探讨人类活动对双湖海岸地貌演化的作用过程、强度、方式，为陵水沙坝溻湖海岸资源的合理开发及可持续利用提供科学依据。回、沉积环境演变和现代海洋作用¨６’１１川。研究区位于海南岛的东南部沿海陵水黎族自治县，东濒南海，南与三亚市毗邻。新村、黎安双渴湖位于陵水县东南角，是两个典型的口小腹大的潮控天然渴湖，双湖腹背相依中间隔一沙坝，南部为南湾猴岛及牛白山形成天然屏障将之与南海隔离（图１）。新村港溻湖有曲港河、曲沟河两条流域面积都很小的淡水河注入，淡水河向该海湾输沙有限¨８。。海湾西部有一约２５０ｍ宽口门与陵水湾相通，涨、落潮三角洲都比较发育，湾内落潮流速大于涨潮流速，使该渴湖能够维持近于封闭海湾的稳定性¨９｜。该溻湖东西长５．１ｋｍ，南北宽４．５ｋｉｎ，岸线总长２６．４１ｋｍ，总面积１９．４３ｋｉｎ２，其中０ｍ等深线以深水域面积１２．５ｋｉｎ２，口门平均水深５．７ｍ，最大水深２研究区概况及数据来源２．１研究区概况海南岛属中国的第二大岛屿，位于南海西北部，属热带季风型岛屿，冬季盛行东北风，夏季盛行西南风，全年气温较高，热量充足，无霜、无寒露，雨季长，旱季短，属湿润地区，多年最低气温大于５。Ｃ，１月份平均气温１９℃以上¨…；雨量丰富，据海南省水资源公报（２０１３年）公布海南省多年平均降水量为１７５０ｍｍ。海南岛属于台风多发地带，平均每年受到８—９个台风影响，素有“台风走廊”之称，是我国乃至西太平洋地区受台风灾害最严重的地区之一¨５’”’。晚更新世以来，海南岛东部经历了复杂的海陆相互作用，包括海平面升降旋１１．２ｍ。其海底地貌类型单一，为一较平坦的浅海淌湖盆地，周围主要为侵蚀残丘、海积阶地和沙堤环绕，北岸为新村码头堤岸，南部有南湾猴岛和牛白山形成的天然屏障。湾内水域面积较大，水深条件良好，水深分布状况大致为南深北浅，底质以砂为主，深水处多为泥，有珊瑚礁零星分布。淌湖口门的渔港湾内ｔｌ然环境优美，渔排养殖大面积分布，漓湖内有３个水上蛋家渔村，是“蛋家文化”的发源地之一¨“。黎安溻湖最宽处约２．８２ｋｉｎ，长约４．３９ｋｍ，岸线总长２１．１８ｋｉｎ，面积约７．９２ｋｍ２，口门位于渴湖东南部，开口向南，东面为一南北走向隔开外海的巨型沙坝，西面为沙嘴，内侧距离口门约２００ｍ处有一长约７０ｍ的涨潮流三角洲。北倚东高岭，周围台图１研究区位置Ｆｉｇ．１Ｌｏｃａｔｉｏｎｏｆｔｈｅｓｔｕｄｙａｒｅａ万方数据１期贾培宏等：海南陵水沙坝渴湖海岸近２６年地貌类型空间格局演化地孤山相间，东高岭和涡湖间有小面积渴湖淤积平原，整个溻湖为海积阶地、丘陵和沙坝所环绕，内部相对平坦，外围为水下浅滩¨引，向内为大面积水下平原，湾内渔塘、渔排养殖业发达。双涡湖主要动力是波浪和潮汐，附近海域潮汐属不正规全Ｅｔ混合潮，据龚文平等¨州观测结果表明，新村溻湖潮流表现出时间一流速不对称特征，实测口门平均涨潮历时为１３小时２０分钟，平均落潮历时为１１小时２０分钟，涨潮历时大于落潮历时。据蒋增杰等¨副观测，黎安溻湖潮流也表现出时间一流速不对称特征，实测口门平均涨潮流历时为１２小时５０分钟，平均落潮流历时为１１小时２０分钟，涨潮历时大于落潮历时。本研究组在双漓湖内通过对南京大学野外获取柱状样的粒度、地球化学等多项指标研究，发现有３５０年来的古风暴记录¨“。研究区内有自然村镇３个：新村镇、黎安镇和三才镇。２００９年１２月３１日国务院印发了《国务院关于推进海南国际旅游岛建设发展的若干意见》；２０１２年１月海南省委省委省政府设立海南国际旅游岛先行试验区（以下简称“试验区”），管理新村及黎安片区的国际滨海旅游新城建设。２．２数据来源本研究数据来源包括近２６年多时相遥感影像及实测陆地、海底地形图、陵水县多年期土地利用数据（２００５—２０１４年）及ＧｏｏｇｌｅＥａｃｈ高分辨率ＳＰＯＴ影像。遥感影像采用Ｌａｎｄｓａｔ系列卫星１９８７—２０１３年间６期影像，如表１。实测地形图为１：１０００ＡｕｔｏＣＡＤＤＷＧ格式数字化地形图，由试验区２０１１～２０１３年委托专业测量人员实测获得。表１遥感影像及其参数Ｔａｂｌｅ１Ｐａｒａｍｅｔｅｒｓｏｆｒｅｍｏｔｅ‘ｓｅｎｓｉｎｇｉｍａｇｅｓ美国陆地探测卫星系统第一颗卫星Ｌａｎｄｓａｔ１于１９７２年发射，目前最新的Ｌａｎｄｓａｔ８于２０１３年２月发射，是目前观测周期最长的影像资料，非常适合做长周期海岸地貌演变分析。本研究收集了研究万方数据区所有可获得Ｌａｎｄｓａｔ影像资料３０多期，经过波段合成处理后，选择以下６期质量好、能反映研究区地貌转化节点的影像作为主要数据源（表１）。３数据处理过程及其方法３．１地貌类型划分通过野外地貌调查，结合《中国海湾志》中新村湾地貌图¨４１以及《海南国际旅游岛先行试验区潮汐汊道海湾动力地貌及双湖连通工程选址研究》报告对双湖沿岸地貌类型的划分方法旧…，根据遥感地貌类型分布特点，将双湖沿岸地貌类型划分为自然地貌类型６种：构造侵蚀中高丘、海积阶地、海积平原、水下浅滩、水下平原和海滩沙堤；人工地貌类型３种：养殖池、盐田和人为堆积高地（表２）。表２陵水沙坝渴湖海岸地貌类型分类标准Ｔａｂｌｅ２ＣｌａｓｓｉｆｙｉｎｇｃｒｉｔｅｒｉｏｎｏｆｓａｎｄｂａｒｒｉｅｒｌａｇｏｏｎｓｇｅｏｍｏｒｐｈｏｌｏｇｉｃｐａｔｔｅｒｎｓｉｎＬｉｎｇｓｈｕｉＣｏｕｎｔｙ３．２地貌类型解译地貌类型遥感影像解译根据Ｌａｎｄｓａｔ影像的７个波段光谱特征、数字高程模型（ＤＥＭ）以及影像纹理的几何特征相结合，以对解译精度有很好的控制。首先，在专业ＡｒｃＧＩＳ１０．１软件中，采用３ＤＡｎａｌｙｓｉｓ模块，对实测地形数据创建研究区数字高程模型（ＤｉｇｉｔａｌＥｌｅｖａｔｉｏｎＭｏｄｅｌ，简称ＤＥＭ），并实现遥感影像与ＤＥＭ的叠加，建立三维虚拟环境，以区分处于不同高程的地貌类型，如：侵蚀高丘高程一般在２５ｍ以上，海积阶地高程在１０～２５ｍ范围内，ｌＯｍ以下则为大面积平坦的海积平原。第四纪研究２０１６拄使用专业遥感影像处理软件ＥＮＶＩ５．１中的ＦＬＡＡＳＨ模块对遥感影像进行大气校正，消除大气散射、吸收及反射引起的误差，获取真实的地表反射率，以分析沙坝渴湖海岸不同地貌类型的光谱特征。本文用于区分空间位置联系紧密的有３种地貌类型：海积平原、水下浅滩及水下平原。海积平原属典型陆地地貌，水下平原属典型海洋地貌，水下程就是对于每个随机产生的决策树分类器，输入特征向量，森林中每棵树对样本进行分类，根据每个树的权重得到最后的分类结果∽２｜。近年来，随机森林法在多源数据和随机小样本的遥感影像分类中获得广泛应用旧３｜。在包括ＬａｎｄｓａｔＥＴＭ＋、多光谱、高光谱和ＧＩＳ数据等多源数据组合，以及遥感土地覆盖等分类应用方面都有优秀的表现旧１‘，其分类计算速度快、对噪声与过训练不敏感等特点弥补了人工神经网络、支持向量机等分类方法的局限性旧…。本研究选择随机森林法，综合多光谱影像、ＤＥＭ和几何纹理等多种分析解译方法，将多源数据纳入到影像分类决策树中，在提高分类精度的同时保证数据处理效率。分类使用的训练样本通过实地考察、陵水县多年期土地利用数据（２００５～２０１４年）及ＧｏｏｇｌｅＥａｒｔｈ高分辨率影像解译获得，每类地物３０个样本，９种地类共２７０个样本。德国ＤｅｆｉｎｉｅｎｓＩｍａｇｉｎｇ公司开发的智能化影像分析软件ｅＣｏｇｎｉｔｉｏｎ提供了实用的随机森林法模块，实现了多数据源融合、多尺度分割的自动分类功能。本研究在该软件中，结合３种决策特征将分类器的大小设定为１００棵决策树，获得１９８７至２０１３年间６期遥感影像的地貌分类结果。为了对ｅＣｏｇｎｉｔｉｏｎ自动提到的海岸地貌分类结浅滩则介于渴湖平原及水下平原之间。因此，水下浅滩与相邻两种地貌仅靠形状、纹理、颜色等几何及可见光色彩特征进行区分难以实现定量化反演。结合野外实地考察资料，以及２０１２年实测地形数据构建的ＤＥＭ，统计分析研究区遥感影像光谱特征后发现，水下浅滩与海积平原及水下平原在光谱特征上有明显差异（图２）。水下浅滩在可见光波段的反射率高于水下平原和海积平原，其红外波段特征介于水下平原及海积平原之间，因此可将其作为随机森林决策树的一个重要决策条件旧“”１。ＴＭｌＴＭ２ＴＭ３ＴＭ４ＴＭ５ＴＭ７波段果进行有效的评估，混淆矩阵常用的几项评价指标有：总体分类精度、Ｋａｐｐａ系数、错分误差、漏分误差等心４｜，本研究采用最常用的总体分类精度和Ｋａｐｐａ系数进行评估。总体分类精度等于被正确分类的像元总和除以总像元数，本研究在开展野外实测年份的２０１３年影像上随机选取２００个真实地貌样本数据，对随机森林分类结果的混淆矩阵进行检验，如表３所示，Ｃ１一Ｃ９为本次研究的分类（说明见表２），被正确分类的像元数目沿着混淆矩阵的对角线分布（总和为１６０７），总像元数等于所有真表３地貌分类结果混淆矩阵Ｔａｂｌｅ３ＣｏｎｆｕｓｉｏｎｍａｔｒｉｘｆｏｒｔｈｅｇｅｏｍｏｒｐｈｏｌｏｇｉｃｃｌａｓｓｉｆｉｃａｔｉｏｎＣ２０１７７８３０５Ｏ３４２００Ｃ３７０１７６０４０１００３２００Ｃ４０１１１８５０３４０６２００Ｃ５Ｏ０６０１７９１２０３０２００Ｃ６０３０６３１８３０５０２００Ｃ７００１０１０１８２１５１２００Ｃ８３０５２５０３１７８４２００Ｃ９Ｏ０１３３１２０２１７９２００Ｈ２水Ｆ、Ｆ原、水下浅滩、海积甲原的光谱特ｆＦｉｇ．２Ｓｐｅｃｔｒａｌｓｉｇｎａｔｕｒｅｓｏｆｕｎｄｅｒｗａｔｅｒｐｌａｉｎｓ，ｍａｒｉｎｅｕｎｄｅｒｗａｔｅｒｓｈｏａｌｓａｎｄｓｅｄｉｍｅｎｔａｒｙｐｌａｉｎｓ养殖池、盐田的光谱特征与水下浅滩、水下平原相近，为避免异物同谱现象对地貌类型解译的干扰，本研究引入几何特征来辨别养殖池与盐田。养殖池与盐田为人工地貌，呈现出规则的格网状，受蓄水量及其他各种人为因素影响，图像纹理粗糙（图３）。类别Ｃ１Ｃ１１６８ＯＯ３０６２０１２合计１７８１８１２１０２０２１９３２ｌｌ２１９２０７１９９１８００（ａ）养殖池（ｂ）水下浅滩（ｃ）水下平原Ｃ２Ｃ３图３养殖池、水下浅滩、水下、Ｆ原纹理特祉Ｆｉｇ．３Ｔｅｘｔｕｒｅｆｅａｔｕｒｅｏｆａｒｔｉｆｉｃｉａｌｐｏｎｄ，Ｃ４Ｃ５Ｃ６Ｃ７ｕｎｄｅｒｗａｔｅｒｓｈｏａｌａｎｄｕｎｄｅｒｗａｔｅｒｐｌａｉｎ３．３随机森林法随机森林是一种比较新的机器学习模型，其过Ｃ８Ｃ９合计２００万方数据１期贾培宏等：海南陵水沙坝渴湖海岸近２６年地貌类型空间格局演化实参考源的像元总数（为１８００），本次总体分类精度＝１６０７／１８００＝８７．５９３３％。Ｋａｐｐａ系数是通过把所有真实参考的像元总数乘以混淆矩阵对角线的和，再减去某一类中真实参考像元数与该类中被分类像元总数之积之后，再除以像元总数的平方减去某一类中真实参考像元总数与该类中被分类像元总数之积对所有类别求和的结果，为０．８８１，可见分类结果精度很高。２０１３年，整个溻湖已基本被养殖池地貌环绕，海积平原及水下浅滩的面积大幅缩小；２００１年始，涡湖东部海积平原与海积阶地上出现零散分布的盐田，后来，并没有大规模连片发展，而是部分废弃退回原始地貌状态。４．２地貌类型面积变化１９８７—２０１３年期间，双湖沿岸不同地貌类型均发生变化，最为显著的是水下浅滩大面积减少，养殖池、盐田和人为堆积高地等人工地貌大量增加；海滩沙堤总面积较小，变化趋势不明显；水下平原面积处于不断波动之中，构造侵蚀中高丘及海积阶４研究结果４．１地貌空间格局演化从多源数据提取地貌类型信息，获得近２６年陵水沙坝渴湖海岸６个年份（见表１）的地貌类型与空间格局演化图（图４）。图４表明，２０世纪８０年代末，研究区主要是以渴湖水下平原为中心的环形地貌格局，自水下平原向外依次是水下浅滩、海积平原、海积阶地、构造侵蚀中高丘，同时水下浅滩与海积平原在空间上相互连接。新村溻湖沿岸在海积平原与水下浅滩之间分布着少量养殖池与盐田，黎安溻湖沿岸基本呈自然地貌状态。９０年代开始，双湖沿岸人工地貌面积迅速增长打破了原有自然地貌格局，并且，由于双湖沿岸地貌资源、地形条件差异较大，两个渴湖虽然空间距离接近，人文环境类似，人工地貌在空间格局特征演化方面却存在着较大差异。桐栖港西侧海积阶地将新村溻湖北岸海积平原分为东西两部分。１９８７年，西部平原上已有较大规模盐田存在，东部平原上零星分布有小面积养殖池，此后一直呈现东西分异格局；至２０１３年东部已经全部被养殖池包围；西部向盐田演化并一直向南延伸，至２００７年，趋于稳定，同时，盐田南部开始出现人为堆积高地。新村涡湖东岸是连岛沙坝，海积平原狭窄，１９９６年大墩村附近出现零星分布的狭长养殖池，随后环湖向北、向南延伸，至２００７年，受海草保护区限制，南部沿岸人工地貌停止扩展，向北延伸的养殖池２０１３年与桐栖港东岸养殖池连接成片，形成一条环湖５．７ｋｍ的养殖池地貌带，将沿岸海积平原与水下浅滩彻底隔离。黎安淌湖沿岸环形地貌格局明显，四周皆有海积阶地及构造侵蚀中高丘阻挡，海积平原与水下浅滩呈狭窄的圈层状，环绕水下平原展布。自１９８７年，湾顶的海积平原开始向养殖池地貌转变；至地近２６年变化不大。自然地貌中水下浅滩变化显著，近２６年间，其面积从１９８７年的１９．１８ｋｍ２缩减到２０１３年的１５．０６ｋｍ２，减少近２１．５％（图５），有持续缩减趋势。其间，水下浅滩面积稍有回升，是因为新村溻湖落潮主水道东侧沙坝自１９８７年来不断向西延伸，至２００１年共向西延伸了近２００ｍ，与水道西侧沙坝连接形成了新村口门外面积达１．７ｋｍ２的水下浅滩¨５｜。双湖沿岸人工地貌面积均有较大幅度变化（图６和表４）。１９８７年双湖沿岸已经开始有养殖池及盐田分布，面积分别为０．４ｋｉｎ２及１．６ｋｉｎ２；至２０１３年，养殖池面积增长到４．６４ｋｉｎ２，是１９８７年的１１．６倍；盐田面积增长到２．９６ｋｉｎ２，是１９８７年的１．８５倍。４．３地貌类型转化通过地貌类型图斑面积提取及地貌类型转移矩阵分析，发现近２６年来，双湖沿岸地貌类型转化趋势主要表现为海积平原、水下浅滩、海积阶地向养殖池及盐田转化（表５）。１９８７—２０１３年，海积平原有１．６３ｋｍ２转化为养殖池、１．１０ｋｍ２转化为盐田；水下浅滩有２．８５ｋｍ２转化为养殖池、０．１ｌｋｍ２转化为盐田，海积阶地有０．０２ｋｍ２转化为养殖池、０．１２ｋｍ２净转化为盐田，共有５．８３ｋｍ２的自然地貌净转化为人工地貌。相反，自然地貌间的相互转化规模较小，主要表现为水下浅滩、海滩沙堤及水下平原向海积平原转化，其转化面积分别为１．０７ｋｍ２、０．６２ｋｍ２和０．１５ｋｍ２，总共１．８４ｋｍ２，远小于自然向人工地貌的转化。可见，人类活动正在作为一种重要的地貌营力，对双湖沿岸地貌加以改造。将地貌转移矩阵分为１９８７～２００１年（表６）和２００１—２０１３年（表７）两个时间段，其转化特征稍有不同，２００１年是地貌类型转换的一个明显时间节万方数据１０８第四纪研究２０１６往图４Ｆｉｇ．４１９８７—２０１３年的近２６年陵水双渴湖海岸地貌类型及空间格局（与图１经纬度范围相同）ｃｏａｓｔＧｅｏｍｏｒｐｈｏｌｏｇｉｃｔｙｐｅｓａｎｄｔｈｅｉｒｓｐａｔｉａｌｐａｔｔｅｒｎｓｏｆｓａｎｄ—ｂａｒｒｉｅｒｌａｇｏｏｎｉｎＬｉｎｓｈｕｉＣｏｕｎｔｙｉｎｔｈｅｓｅ２６ｙｅａｒｓ（１９８７～２０１３）８０５ｌ４６０——３４０２００海积平原水下平原震熊飘圉水下浅滩龟１／醛喧２０养殖池盐田人为堆积高地图５Ｆｉｇ．５自然地貌面积变化图６人｜［地貌面积变化Ａｒｅａｃｈａｎｇｅｓｏｆｎａｔｕｒａｌｇｅｏｍｏｒｐｈｏｌｏｇｙ万方数据１期贾培宏等：海南陵水沙坝渴湖海岸近２６年地貌类型空间格局演化１０９点。１９８７～２００１年间，海积平原倾向于向养殖池转表４人工地貌面积变化Ｔａｂｌｅ化，转化面积为１．４７ｋｍ２，向盐田转化的面积仅为４Ａｒｅａｃｈａｎｇｅｓｏｆａｒｔｉｆｉｃｉａｌｇｅｏｍｏｒｐｈｏｌｏｇｙ０．８９ｋｍ２；２００１—２０１３年间，海积平原向养殖池及盐田的转化面积均下降至０．２ｋｍ２左右；海积阶地的转化主要发生在１９８７—２００１年间向盐田转化０．０９ｋｍ２，占近２６年总转化面积的８１．２％；水下浅滩向养殖池及盐田转化的高峰时期在２００１—２０１３年间，向养殖池转化１．５９ｋｍ２、向盐田转化０．１１ｋｍ２。２００１年以前沿岸居民主要选择陆地地貌加以人工改造，人工地貌主要由海积平原与海积表５１９８７—２０１３年地貌类型转移矩阵（面积／ｋｉｎ２）Ｔａｂｌｅ５Ｔｈｅｔｒａｎｓｆｅｒｍａｔｒｉｘｏｆｔｈｅｇｅｏｍｏｒｐｈｏｌｏｇｉｃｔｙｐｅｓｆｒｏｍ１９８７ｔｏ２０１３（ａｒｅａ／ｋｍ２）海积阶地１５．１ｌ０．０００．０００．０００．０００．０２０．Ｏｌ０．００１５．１４海积平原０．００７３．６００．１００．０１０．２００．９３１．４８０．００７６．３２海滩沙堤０．０００．３３１．１４０．０００．１３０．０００．０００．００１．６０水下平原０．０００．０００．０６４７．１０２．６７０．０００．０００．００４９．８３水下浅滩０．０００．８２０．１９１．１１１６．０９０．０１０．９６０．００１９．１８盐田０．０００．０４０．０００．０００．０ｌ１．５５０．０００．００１．６０养殖池０．０００．０１０．０００．０００．０００．０００．３８０．０００．４０构造侵蚀中高丘Ｏ．００Ｏ．００Ｏ．ｏｏＯ．００Ｏ．０００．ｏｏＯ．００２２．０７２２．０７合计１５．１１７４．８１１．４８４８．２２１９．１０２．５２２．８４２２．０７１８６．１５万方数据１１０第四纪研究２０１６正图７Ｆｉｇ．７１９８７～２０１３年人上地貌类型空问分布ｔｏＴｈｅｓｐａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆａｒｔｉｆｉｃｉａｌｇｅｏｍｏｒｐｈｏｌｏｇｙｆｒｏｍ１９８７２０１３阶地转入；２００１年以后水下浅滩成为人类活动的主要改造对象。以１９８７年海积平原与水下浅滩分界线为界，将２０１３年环绕双湖的人工地貌带划分为海积平原转入区域（向陆）及水下浅滩转入区域（向湖）（图７），计算出１９８７～２０１３年双湖人工地貌的转入来源（表８和表９），分析人工地貌演化特征在空间上的差异性。表８新村渴湖沿岸人工地貌转入来源（面积／ｋｍ２）Ｔａｂｌｅ８Ｔｈｅｓｏｕｒｃｅ２５．６２％；水下浅滩转入黎安占６８．２２％，新村占３７．４５％。与地貌格局演化的差异性相似，人工地貌转入也因双湖沿岸地形差异造成。新村溻湖北侧与西侧大面积海积平原，为人工地貌的演化提供了优质资源，桐栖港北部曲沟河的泾流注入，为沿河扩张的养殖池提供了排水便利。黎安溻湖东西两侧均为狭长阶地，南北两侧有侵蚀高丘，海积平原狭窄，养殖池多通过对水下浅滩挖掘形成。５结论（１）研究结果表明，随机森林法在海岸带ｏｆａｒｔｉｆｉｃｉａｌｇｅｏｍｏｒｐｈｏｌｏｇｙａｒｏｕｎｄｔｈｅＸｉｎｃｕｎＬａｇｏｏｎ（ａｒｅａ／ｋｍ２）Ｌａｎｄｓａｔ多光谱遥感影像地貌分类中解译精度很高，结合ＧＩＳ数字高程模型时空分析技术，及地貌类型转移矩阵数理分析方法，对沙坝渴湖海岸地貌类型演化展开科学研究，是一种可行、有效的新方法。该方法弥补了传统采用钻孔及沉积柱状样采集分析中以点代面、难以全域覆盖的海岸带沉积地貌演化表９黎安漓湖沿岸人工地貌转入来源（面积／ｋｍ２）Ｔａｂｌｅ９Ｔｈｅｓｏｕｒｃｅ分析方法的不足。（２）分析发现，１９８７～２０１３年近２６年间，陵水沙坝溻湖海岸地貌发生了显著变化，地貌类型由海积平原、水下浅滩、海积阶地向养殖池、盐田转化；自然地貌间相互转化不明显；大范围潮滩围垦和海积平原挖掘改造形成了一条人工地貌带，打破了原来的地貌空间格局，割断了海积平原与水下浅滩之间的沉积动力作用方式和物质能源输运模式；人类活动已成为研究区海岸地貌演变的主要驱动力。ｏｆａｒｔｉｆｉｃｉａｌｇｅｏｍｏｒｐｈｏｌｏｇｙａｒｏｕｎｄｔｈｅＬｉ—ＡｎＬａｇｏｏｎ（ａｒｅａ／ｋｍ２）计算结果显示，双湖沿岸人工地貌由海积平原和水下浅滩转入的总比例新村占９９．８４％，黎安占９３．８４％；海积平原转入新村占６２．２９％，黎安占（３）人工地貌的转入来源在时空格局上具有明显差异。新村渴湖周边大面积的海积平原地貌是养殖池与盐田的主要转化来源；黎安渴湖四周的侵蚀万方数据１期贾培宏等：海南陵水沙坝涡湖海岸近２６年地貌类型空间格局演化１１１中高丘与海积阶地，并不适合养殖池与盐田的开发，海积平原面积有限，可开发空间小，因此水下浅滩是其人工地貌的主要转入源。无论从时间顺序还是地貌空间格局上，海积平原总是研究区人类活动的优先改造对象，当其开发殆尽、资源不足时，水下浅滩成为下一个开发目标，其中，２００１年是地貌类型转换的一个明显时间节点，将地貌转移矩阵分为１９８７～２００１年（表６）和２００１～２０１３年（表７）两个时间段，其转化特征稍有不同。２００１年以前沿岸居民主要选择陆地地貌加以人工改造，人工地貌主要由海积平原与海积阶地转入；２００１年以后水下浅滩成为人类活动的主要改造对象。参考文献（Ｒｅｆｅｒｅｎｃｅｓ）１恽才兴，韩震，蒋雪中等．海岸带及近海卫星遥感综合应用技术．北京：海洋出版社，２００５．１—８ＹｕｎＣａｉｘｉｎｇ，ＨａｎＺｈｅｎ，ＪｉａｎｇＸｕｅｚｈｏｎｇｅｔａ１．ＩｎｔｅｇｒａｔｅｄＡｐｐｌｉｃａｔｉｏｎｏｆＣｏａｓｔａｌＺｏｎｅａｎｄＯｆｆｓｈｏｒｅＳａｔｅｌｌｉｔｅＲｅｍｏｔｅＳｅｎｓｉｎｇＴｅｃｈｎｏｌｏｇｙ．Ｂｅｉｊｉｎｇ：ＣｈｉｎａＯｃｅａｎＰｒｅｓｓ．２００５．１～８２谢志仁，袁林旺．略论全新世海面变化的波动性及其环境意义．第四纪研究，２０１２，３２（６）：１０６５—１０７７ＸｉｅＺｈｉｒｅｎ，ＹｕａｎＬｉｎｗａｎｇ．ＦｌｕｃｔｕａｔｉｏｎｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆＨｏｌｏｃｅｎｅｓｅａ－ｌｅｖｅｌｃｈａｎｇｅａｎｄｉｔｓｅｎｖｉｒｏｎｍｅｎｔａｌｉｍｐｌｉｃａｔｉｏｎｓ．ＱｕａｔｅｒｎａｒｙＳｃｉｅｎｃｅｓ，２０１２，３２（６）：１０６５—１０７７３程和琴，陈祖军，阮仁良等．海平面变化与城市安全——以上海市为例．第四纪研究，２０１５，３５（２）：３６３～３７３ＣｈｅｎｇＨｅｑｉｎ，ＣｈｅｎＺｕｊｕｎ，ＲｕａｎＲｅｎｌｉａｎｇｅｔａ１．Ｓｅａｌｅｖｅｌｃｈａｎｇｅａｎｄｃｉｔｙｓａｆｅｔｙ——ＴｈｅＳｈａｎｇｈａｉａｓａｎｅｘａｍｐｌｅ．ＱｕａｔｅｒｎａｒｙＳｃｉｅｎｃｅｓ，２０１５，３５（２）：３６３～３７３４戴志军，施伟勇，陈浩．沙坝一滴湖海岸研究进展与展望．上海国土资源，２０１１，３２（３）：１２一１７ＤａｉＺｈｉｊｕｎ，ＳｈｉＷｅｉｙｏｎｇ，ＣｈｅｎＨａｏ．Ｒｅｓｅａｒｃｈｐｒｏｇｒｅｓｓａｎｄｐｒｏｓｐｅｃｔｓｏｆｂａｒｒｉｅｒ－ｌａｇｏｏｎｃｏａｓｔ．ＳｈａｎｇｈａｉＬａｎｄａｎｄＲｅｓｏｕｒｃｅｓ，２０１１，３２（３）：１２～１７５张乔民，陈欣树，王文介等．华南海岸沙坝漏湖型潮汐汉道口门地貌演变．海洋学报，１９９５，１７（２）：６９～７７ＺｈａｎｇＱｉａｏｍｉｎ，ＣｈｅｎＸｉｎｓｈｕ，ＷａｎｇＷｅｎｊｉｅｅｔａ１．Ｇｅｏｍｏｒｐｈｏｌｏｇｉｃａｌｅｖｏｌｕｔｉｏｎｏｆｔｈｅｅｎｔｒａｎｃｅｓｏｆｔｈｅｓａｎｄｂａｒ－ｌａｇｏｏｎｔｙｐｅｔｉｄａｌｉｎｌｅｔｓａｌｏｎｇｔｈｅＳｏｕｔｈｅｒｎＣｈｉｎａｃｏａｓｔｌｉｎｅｓ．ＡｃｔｓＯｃｅａｎｏｌｏｇｉｅａＳｉｎｉｃａ，１９９５，１７（２）：６９～７７６孙伟富．我国海岸溜湖遥感监测与典型渴湖分析．青岛：中国海洋大学博士学位论文，２０１３．２～２６ＳｕｎＷｅｉｆｕ．ＭｏｎｉｔｏｒｉｎｇａｎｄＡｎａｌｙｚｉｎｇｔｈｅＬａｇｏｏｎＤｙｎａｍｉｃｓｉｎＣｈｉｎａＵｓｉｎｇＲｅｍｏｔｅＳｅｎｓｉｎｇＩｍａｇｅｒｙ．Ｑｉｎｇｄａｏ：ＴｈｅＤｏｃｔｏｒＤｉｓｓｅｒｔａｔｉｏｎｏｆＯｃｅａｎＵｎｉｖｅｒｓｉｔｙｏｆＣｈｉｎａ．２０１３．２—２６７ＥｌｉａｓＥＰＬ。ＳｐｅｋＡＪＦ．Ｌｏｎｇ—ｔｅｒｍｍｏｒｐｈｏｄｙｎａｍｉｃｅｖｏｌｕｔｉｏｎｏｆｔｅｘｅｌｉｎｌｅｔａｎｄｉｔｓｅｂｂ－ｔｉｄａｌｄｅｌｔａ（ｔｈｅＮｅｔｈｅｒｌａｎｄｓ）．ＭａｒｉｎｅＧｅｏｌｏｇｙ，２００６，２２５（１～４）：５—２ｌ８ＷｏｏｄｒｏｆｆｅＣＤ．Ｃｏａｓｔｓ：Ｆｏｒｍ，ＰｒｏｃｅｓｓａｎｄＥｖｏｌｕｔｉｏｎ．Ｌｏｎｄｏｎ：ＣａｍｂｒｉｄｇｅＵｎｉｖｅｒｓｉｔｙＰｒｅｓｓ．２００２．４５～５２９杨红强，余克服．微环礁的高分辨率海平面指示意义．第四纪研万方数据究，２０１５，３５（２）：３５４—３６２ＹａｎｇＨｏｎｇｑｉａｎｇ，ＹｕＫｅｆｎ．Ｈｉｇｈ－ｒｅｓｏｌｕｔｉｏｎｓｅａｌｅｖｅｌｃｈａｎｇｅｓｉｎｄｉｃａｔｉｖｅｉｍｐｌｉｃａｔｉｏｎｓｏｆｍｉｃｒｏａｔｏｌｌ．ＱｕａｔｅｒｎａｒｙＳｃｉｅｎｃｅｓ，２０１５，３５（２）：３５４—３６２１０高抒．潮汐汉道开发中的水环境问题．水资源保护，２００２，１８（３）：１８—２１ＧａｏＳｈｕ．Ｗａｔｅｒｅｎｖｉｒｏｎｍｅｎｔｐｒｏｂｌｅｍｓｉｎｄｅｖｅｌｏｐｍｅｎｔｏｆｔｉｄａｌｉｎｌｅｔｓｙｓｔｅｍｓ．ＷａｔｅｒＲｅｓｏｕｒｃｅｓＰｒｏｔｅｃｔｉｏｎ，２００２，ｌＳ（３）：１８—２１１１ＤａｖｉｓＲＡ，ＹａｌｅＫＥ，ＰｅｋａｌａＪＭｅｔａ１．ＢａｒｒｉｅｒｉｓｌａｎｄｓｔｒａｔｉｇｒａｐｈｙａｎｄＨｏｌｏｃｅｎｅｈｉｓｔｏｒｙｏｆｗｅｓｔ－ｃｅｎｔｒａｌＦｌｏｒｉｄａ．ＭａｒｉｎｅＧｅｏｌｏｇｙ，２００３，２００（１～４）：１０３～１２３１２李春初，罗宪林，张镇元等．粤西水东沙坝渴湖海岸体系的形成演化．科学通报，１９８５，３１（２０）：１５７９～１５８２ＬｉＣｈｕｎｃｈｕ，ＬｕｏＸｉａｎｌｉｎ，ＺｈａｎｇＺｈｅｎｙｕａｎｅｔａ１．Ｆｏｒｍａｔｉｏｎａｎｄｅｖｏｌｕｔｉｏｎｏｆｔｈｅｅａｓｔｗａｔｅｒｂａｒｒｉｅｒ－ｌａｇｏｏｎｓｙｓｔｅｍｓｉｎｗｅｓｔｅｒｎＧｕａｎｇｄｏｎｇ．ＣｈｉｎｅｓｅＳｃｉｅｎｃｅＢｕｌｌｅｔｉｎ，１９８５，３１（２０）：１５７９～１５８２１３陵水黎族自治县地方志编纂委员会．陵水县志．北京：方志出版社．２００７．１～５０ＬｏｃａｌＣｈｒｏｎｉｃｌｅｓＣｏｍｐｉｌａｔｉｏｎＣｏｍｍｉｔｔｅｅｏｆＬｉｎｓｈｕｉＣｏｕｎｔｙ．ＣｏｕｎｔｙＣｈｒｏｎｉｃｌｅｏｆＬｉｎｓｈｕｉ．Ｂｅｉｊｉｎｇ：ＰｕｂｌｉｓｈｉｎｇＨｏｕｓｅｏｆＬｏｃａｌＲｅｃｏｒｄｓ，２００７．１～５０１４中国海湾志编纂委员会．中国海湾志・第十一分册・海南省海湾．北京：海洋出版社，１９９９．１—４２６ＣｈｉｎｅｓｅＣｏｍｐｉｌａｔｉｏｎＣｏｍｍｉｔｔｅｅｏｆＣｈｉｎａ’ｓＣｏａｓｔＥｍｂａｙｍｅｎｔｓ．Ｃｈｉｎａ’ｓＣｏａｓｔＥｍｂａｙｍｅｎｔｓ（Ｖ０１．１１：ＨａｉｎａｎＣｏａｓｔＥｍｂａｙｍｅｎｔｓ）．Ｂｅｉｊｉｎｇ：ＣｈｉｎａＯｃｅａｎＰｒｅｓｓ。１９９９．１—４２６１５周亮，高抒，杨阳等．海南岛东南部３５０年古风暴事件沉积与历史文献记录对比．海洋学报，２０１５，３７（９）：８４～９４ＺｈｏｕＬｉａｎｇ，ＧａｏＳｈｕ，ＹａｎｇＹａｎｇｅｔａ１．Ｃｏｍｐａｒｉｓｏｎｏｆｐａｌｅｏｓｔｏｒｍｅｖｅｎｔｓｂｅｔｗｅｅｎｓｅｄｉｍｅｎｔａｒｙａｎｄｈｉｓｔｏｒｉｃａｌａｒｃｈｉｖｅｓ：Ａ３５０ｙｅａｒｒｅｃｏｒｄｆｒｏｍｓｏｕｔｈｅａｓｔｅｒｎＨａｉｎａｎＩｓｌａｎｄｃｏａｓｔａｌｅｍｂａｙｍｅｎｔｓ．ＡｃｔａＯｃｅａｎｏｌｏｇｉｃａＳｉｎｉｃａ，２０１５，３７（９）：８４－９４１６刘艳霞，祁雅莉，黄海军等．渤海西南岸６０００年古海岸线重建及影响因素解析．第四纪研究，２０１５，３５（２）：２６５—２７４ＬｉｕＹａｎｘｉａ，ＯｉＹａｌｉ，ＨｕａｎｇＨａｉｊｕｎｅｔａ１．ＲｅｃｏｎｓｔｒｕｃｔｉｏｎａｎｄｄｉｆｆｅｒｅｎｔｉａｔｉｏｎｆａｃｔｏｒｏｆｐａｌｅｏｓｈｏｒｅｌｉｎｅｆｏｒｔｈｅｓｏｕｔｈｗｅｓｔｅｒｎｏｆｔｈｅＢｏｈａｉＳｅａｉｎ６ｋａＢ．Ｐ．ＱｕａｔｅｒｎａｒｙＳｃｉｅｎｃｅｓ，２０１５，３５（２）：２６５—２７４１７陈中原，刘演，陈静等．大河三角洲全新世海平面研究中的若干重要问题——以长江、尼罗河和密西西比河三角洲为例．第四纪研究，２０１５，３５（２）：２７５—２８０ＣｈｅｎＺｈｏｎｇｙｕａｎ，ＬｉｕＹａｎ，ＣｈｅｎＪｉｎｇｅｔａ１．ＴｈｅｋｅｙｉｓｓｕｅｓｏｆＨｏｌｏｃｅｎｅｓｅａ・ｌｅｖｅｌｓｔｕｄｙｏｆｍｅｇａ－ｄｅｌｔａｓ：ＥｘａｍｐｌｅｓｆｒｏｍＣｈａｎｇｊｉａｎｇＲｉｖｅｒ，ｔｈｅＮｉｌｅＲｉｖｅｒａｎｄｔｈｅＭｉｓｓｉｓｓｉｐｐｉＲｉｖｅｒ．ＱｕａｔｅｒｎａｒｙＳｃｉｅｎｃｅｓ，２０１５，３５（２）：２７５—２８０１８蒋增杰，方建光，张继红等．海南黎安港纳潮量及海水交换规律研究．海南大学学报自然科学版，２００９，２７（３）：２６１～２６４ＪｉａｎｇＺｅｎｇｊｉｅ，ＦａｎｇＪｉａｎｇｕａｎｇ，ＺｈａｎｇＪｉｈｏｎｇｅｔａ１．ＴｈｅｃａｌｃｕｌａｔｉｏｎｏｆｔｉｄａｌｗａｔｅｒｃａｐａｃｉｔｙａｎｄｗａｔｅｒｅｘｃｈａｎｇｅｃｈａｒａｃｔｅｒｉｓｔｉｃｏｆＬｉａｎＬａｇｏｏｎ．ＮａｔｕｒａｌＳｃｉｅｎｃｅＪｏｕｒｎａｌｏｆＨａｉｎａｎＵｎｉｖｅｒｓｉｔｙ，２００９，２７（３）：２６１—２６４１９龚文平，ＳｈｅｎＪｉａｎ，陈斌．用一维水力学方程求取泻湖水位及潮汐汉道断面流速——以海南陵水新村港为例．台湾海峡，２００７，２６（３）：３０１—３１３１１２第Ｂｉｎ．Ｏｂｔａｉｎｉｎｇｗａｔｅｒ四纪研究２０ｔ６拄ＧｏｎｇＷｅｎｐｉｎｇ，ＳｈｅｎＪｉａｎ，Ｃｈｅｎｅｌｅｖａｔｉｏｎｉｎ２３雷震．随机森林及其在遥感影像处理中应用研究．上海：上海交通大学博士学位论文，２０１２．２—１６０ＬｅｉＺｈｅｎ．ＲａｎｄｏｍＦｏｒｅｓｔａｎｄＩｔｓＡｐｐｌｉｃａｔｉｏｎｉｎＲｅｍｏｔｅＳｅｎｓｉｎｇ．ｔｈｅｌａｇｏｏｎａｎｄｃｒｏｓｓ—ｓｅｃｔｉｏｎａｌｌｙａｖｅｒａｇｅｄｖｅｌｏｃｉｔｙｏｆｔｈｅｔｉｄａｌｉｎｌｅｔｂｙｕｓｉｎｇｏｎｅ－ｄｉｍｅｎｓｉｏｎａｌｈｙｄｒａｕｌｉｃｅｑｕａｔｉｏｎ——ＡｃａｓｅｓｔｕｄｙｉｎＸｉｎｃｕｎｉｎｌｅｔＴａｉｗａｎＬｉｎｓｈｕｉ，Ｈａｉｎａｎ，Ｃｈｉｎａ．ＪｏｕｒｎａｌｏｆＯｃｅａｎｏｇｒａｐｈｙｉｎＳｈａｎｇｈａｉ：ＴｈｅＤｏｃｔｏｒＤｉｓｓｅｒｔａｔｉｏｎｏｆ２０１２．２—１６０ＳｈａｎｇｈａｉＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｓｔｒａｉｔ，２００７，２６（３）：３０１～３１３２０高抒，贾培宏，殷勇等．海南国际旅游岛先行试验区潮汐汉道海湾动力地貌及双湖连通工程选址研究，２０１４．１—２００ＧａｏＳｈｕ，ＪｉａＰｅｉｈｏｎｇ，ＹｉｎＹｏｎｇｅｔ２４张晓东，李德仁，龚健雅等．遥感影像与ＧＩＳ分析相结合的变化检测方法．武汉大学学报（信息科学版），２００６，３１（３）：２６６～２６９ＺｈａｎｇＸｉａｏｄｏｎｇ，ＬｉＤｅｒｅｎ，ＧｏｎｇＪｉａｎｙａｍｅｔｈｏｄｏｆｉｎｔｅｇｒａｔｉｎｇＳｃｉｅｎｃｅｒｅｍｏｔｅｓｅｎｓｉｎｇｅｔａ１．ＳｔｕｄｙｏｎＧｅｏｍｏｒｐｈｏ—ｄｙｎａｍｉｃｏｆｔｗｏＴｉｄａｌＩｎｌｅｔｓａｎｄＬｏｃａｔｅｔｈｅＣｏｎｎｅｃｔｉｎｇＳｉｔｅｏｆＩｄｌｇｏｏｎｓｉｎＨａｉｎａｎａ１．ＡｃｈａｎｇｅｄｅｔｅｃｔｉｏｎａｎｄＩｎｔｅｒｎａｔｉｏｎａｌＴｏｕｒｉｓｍＩｓｌａｎｄＰｉｌｏｔＺｏｎｅ，２０１４．１－２００２１ＣｈａｎＣａｎｄＧＩＳ．ＧｅｏｍａｔｉｃｓＷ．ＰａｅｌｉｎｃｋｘＤ．ＥｖａｌｕａｔｉｏｎｏｆｒａｎｄｏｍｆｏｒｅｓｔａｎｄａｄａｂｏｏｓｔｃｌａｓｓｉｆｉｃａｔｉｏｎａｉｒｂｏｒｎｅａｎｄｓｐｅｃｔｒａｌｂａｎｄｓｅｌｅｃｔｉｏｎｆｏｒＩｎｆｏｒｍａｔｉｏｎｏｆＷｕｈａｎＵｎｉｖｅｒｓｉｔｙ，２００６，３１（３）：２６６—２６９ｔｒｅｅ—ｂａｓｅｄｅｎｓｅｍｂｌｅｅｅｏｔｏｐｅＳｅｎｓｉｎｇ２２ｍａｐｐｉｎｇ２５龚文平，陈明和，温晓骥等．海南陵水新村港潮汐汉道演变及其稳定性分析．热带海洋学报，２００４，２３（４）：２５—３２Ｇｏｎｇｕｓｉｎｇｈｙｐｅｒｓｐｅｅｔｒａｌｉｍａｇｅｒｙ．ＲｅｍｏｔｅｏｆＥｎｖｉｒｏｎｍｅｎｔ，２００８．１１２（６）：２９９９—３０１１ＴＣ，ＢｅａｒｄＫＨｅｔＷｅｎｐｉｎｇ，ＣｈｅｎＭｉｎｇｈｅ，ＷｅｎＸｉａｏｊｉＸｉｎｃｕｎｅｔａ１．ＥｖｏｌｕｔｉｏｎａｎｄＣｕｔｌｅｒＤＲ，Ｅｄｗａｒｄｓｃｌａｓｓｉｆｉｃａｔｉｏｎｉｎａ１．ＲａｎｄｏｍｆｏｒｅｓｔｓｆｏｒｓｔａｂｉｌｉｔｙｏｆＪｏｕｒｎａｌｔｉｄａｌｉｎｌｅｔ，ＬｉｎｓｈｕｉＣｏｕｎｔｙ，ＨａｉｎａｎＰｒｏｖｉｎｃｅ．ｅｃｏｌｏｇｙ．Ｅｃｏｌｏｇｙ，２００７，船（１１）：２７８３—２７９２ｏｆＴｒｏｐｉｃａｌＯｃｅａｎｏｇｒａｐｈｙ，２００４，２３（４）：２５—３２ＴＨＥＥＶｏＬＵＴＩｏＮＯＦＣｏＡＳＴＡＬＡＮＤＩＴＳＧＥＯＭＯＲＰＨｏＬＯＧＹＰＡＴＴＥＲＮＩＮＲＥＣＥＮＴ２６ＹＥＡＲＳＳＰＡＴＩＡＬ－ＴＥＭＰｏＲＡＬＤＩＳＴＲＩＢＵＴＩＯＮＡＲｏＵＮＤＴＨＥＴＷＯＬＡＧＯＯＮＡＬＡＲＥＡＳＩＮＬＩＮＧＳＨＵＩＣＯＵＮＴＹ．ＨＡＩＮＡＮＩＳＬＡＮＤＪｉａＰｅｉｈｏｎｇ①②③ＸｕＷｅｉ①③ＬｉＧｏｎｇｃｈｅｎｇ＠⑤（①ＳｃｈｏｏｌｏｆＧｅｏｇｒａｐｈｗａｎｄＯｃｅａｎｏｇｒａｐｈｉｃＹｉｎＹｏｎｇ①②③ＬｉＸｉｎ④Ｓｃｉｅｎｃｅｓ，ＮａｎｊｉｎｇＵｎｉｖｅｍｉ移，Ｎａ蛳ｎｇ２１００２３；②ＴｈｅＫｅｙＬａｂｏｒａｔｏｒｙｏｆＣｏａｓｔ＆ｂｌａｎｄＤｅｖｅｌｏｐｍｅｎｔＩｎｎｏｖａｔｉｏｎＣｅｎｔｅｒｏｆＭｉｎｉｓｔｒｙｏｆＥｄｕｃａｔｉｏｎ，ＮａｎｊｉｎｇＵｎｉｔ＊ｒｓｉｔｙ，Ｎａｎｊｉｎｇ２１００２３；③ＣｏｌｌａｂｏｒａｔｉｔｅⅣ口彬ｎｇ２１００２３；④ＨａｉｎａｎＩｎｔｅｒｎａｔｉｏｎａｌＴｏｕｒｉｓｍｂｌａｎｄｏｆＳｏｕｔｈＣｈｉｎａＳｅａＳｔｕｄｉｅｓ，ＮａｎｊｉｎｇＵｎｉｔｍｉｔｙ，Ｐｆ‰Ｚｏｎｅ，Ｌｉｎｇｓｈｕｉ５７２４００；⑤ＣｏｌｌｅｇｅｏｆＡｒｃｈｉｔｅｃｔｕｒｅａｎｄＵｒｂａｎＰｌａｎｎｉｎｇ，ＴｏｎａｌＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ２０００９２）ＡｂｓｔｒａｃｔＴｈｅｒｅｓｅａｒｃｈｅｓｔａｂｌｉｓｈｅｓｔｗｏａｐｒｏｆｅｓｓｉｏｎａｌｃｒｉｔｅｒｉｏｎｔｏｃｌａｓｓｉｆｙｔｈｅｃｏａｓｔａｌｇｅｏｍｏｒｐｈｏｌｏｇｙｐａｔｔｅｒｎｓａｆｔｅｒｓｕｒｖｅｙｉｎｇａｔｓａｎｄ—ｂａｒｒｉｅｒｌａｇｏｏｎｓ，ｕｓｉｎｇｌａｎｄ－ｃｏｖｅｒｄａｔａｄｕｒｉｎｇ２００５—２０１４，ｂａｔｈｙｍｅｔｒｙｄａｔａａｎｄｄｉｇｉｔａｌｔｅｒｒａｉｎｍｏｄｅｌｄａｔａａｓｓｃａｌｅｄ１：１０００ａｎｄｆａｍａｔｅｄａｓＡｕｔｏＣＡＤＤＷＧ，ｉｎＬｉｎｇｓｈｕｉＣｏｕｎｔｙ，ＨａｉｎａｎＩｓｌａｎｄ．Ａｃｃｏｒｄｉｎｇａｓｔｏｔｈｉｓｃｒｉｔｅｒｉｏｎ，ｔｈｅｇｅｏｍｏｒｐｈｏｌｏｇｙｐａｔｔｅｒｎｓｈａｖｅｂｅｅｎｉｄｅｎｔｉｆｉｅｄａｎｄｃｌａｓｓｉｆｉｅｄｏｃｅａｎ６ｔｙｐｅｓｎａｔｕｒａｌｇｅｏｍｏｒｐｈｏｌｏｇｙ：ｅｒｏｓｉｏｎａｌｈｉｌｌｔｅｃｔｏｎｉｃｓ，ｓｈｏａｌｓａｎｄｍａｒｉｎｅｓｅｄｉｍｅｎｔａｒｙｔｅｒｒａｃｅ，ｓｅｄｉｍｅｎｔａｒｙｐｌａｉｎ，ｕｎｄｅｒｗａｔｅｒｐｌａｉｎ，ａｎｄ３ｔｙｐｅｓａｒｔｉｆｉｃｉａｌｇｅｏｍｏｒｐｈｏｌｏｇｙ：ａｑｕｅ－ｆａｒｍｉｎｇｐｏｎｄ，ｓａｌｔｅｒｎ，ｍａｎ－ｍａｄｅｓｔａｃｋｉｎｇｈｉｇｈｌａｎｄ；ｔｈｅｎｔｈｅｉｒｄｉｓｔｒｉｂｕｔｉｎｇｂｏｕｎｄａｒｉｅｓｈａｖｅｂｅｅｎｅｘｔｒａｃｔｅｄｏｎｔｈｅ６一ｓｔａｇｅＬａｎｄｓａｔＲＳｉｍａｇｅｓｆｒｏｍ１９８７ａｎｔｏ２０１３ｔｏｒｅｓｅａｒｃｈｔｈｅｅｎｏｖａｔｉｏｎｃｈａｒａｃｔｅｒｓａｎｄｉｔｓｓｈｉｆｔｉｎｇｄｉｒｅｃｔｉｏｎｂｙｅＣｏｇｎｉｔｉｏｎｓｏｆｔｅｒｗａｒｅ，ｗｈｉｃｈｉｓｉｎｔｅｌｌｉｇｅｎｔｉｍａｇｅａｎａｌｙｚｉｎｇｓｏｆｔｅｒｗａｒｅ，ｐｒｏｄｕｃｅｄｉｎＧｅｒｍａｎ，ｃｏｍｂｉｎｅｄｔｈｅＧＩＳａｎａｌｙｚｉｎｇｔｅｃｈｎｏｌｏｇｙａｎｄｔｈｅｔｒａｎｓｆｅｒｉｎｇｍａｔｒｉｘｍｅｔｈｏｄｓｂｙｐｕｔｔｉｎｇｄｅｃｉｓｉｔｉｏｎｃｌａｓｓｆｉｅｒｔｒｅｅｓｍｕｌｔｉ・ｓｏｕｒｃｅｄａｔａｉｎｔｏｔｈｅｔｏｏｆｉｍａｇｅｃｌａｓｓｉｆｉｃａｔｉｏｎｔｒｅｅｓｔｏｓｅｌｅｃｔ２７０ｓａｍｐｌｅｓｗｉｔｈ９ｔｙｐｅｓａｎｄ３０ｓａｍｐｌｅｓｐｅｒｔｙｐｅａｎｄｄｅｃｉｓｉｏｎｆｅａｔｕｒｅｓ．Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｏｎｄｅｓｉｇｎｔｈｅａｓ１００ｄｅｃｉｓｉｏｎｗｉｔｈ３一ｔｙｐｅｔｈａｔｔｈｅｍｅｔｈｏｄ．ｉｎｔｅｒｇｒａｔｅｄｔｈｅａｎｄｔｒａｎｓｆｏｆｉｎｇｍａｔｒｉｘｍａｋｅｓｕｐｔｈｅｓｕｒｖｅｙｉｎｇａｔｓｔａｔｉｓｔｉｃｓ，ｉｓｆｉｅｌｄ，ｒａｎｄｏｍｆｏｒｅｓｔｅｘｔｒａｃｔｉｏｎａｎＲＳｉｍａｇｅｓ，ＧＩＳｏｎ３Ｄａｎａｌｙｚｉｎｇｔｅｃｈｎｏｌｏｇｙｌｏｎｇ－ｔｅｒｍｏｕｔｅｆｆｅｃｔｉｖｅａｎｄｃｏｒｅｃａｐｔｉｂｌｅｔｏｒｅｓｅａｒｃｈｔｈｅｇｅｏｇｒａｐｈｏｌｏｇｉｃｅｖｏｌｕｔｉｏｎ．Ｉｔｒｅｓｅａｒｃｈｉｎｇｍｅｔｈｏｄｏｆａｎｄｓｅｎｄｉｍｅｎｔａｒｙｃｏｌｕｍｎａｎａｌｙｚｉｎｇ．Ｔｈｅｒｅｓｅａｒｃｈｆｉｎｄｓｔｈｅｃｏａｓｔａｌｇｅｏｍｏｒｐｈｏｌｏｇｙｈａｓｔｅｒｒａｃｅｃｈａｎｇｅｄｒｅｍａｒｋａｂｌｙｉｎｔｈｅｓｅ２６ｙｅａｒｓ．Ｔｈｅｓｅｄｉｍｅｎｔａｒｙｔｒａｎｓｆｅｒｒｉｎｇｔｏｐｌａｉｎ，ｕｎｄｅｒｗａｔｅｒｎｏｔｂａｎｋａｎｄｓｅｄｉｍｅｎｔａｒｙｈａｖｅｂｅｅｎｆａｒｍｉｎｇｐｏｎｄａｎｄｓａｌｔｆｉｅｌｄｔｏｃｒｅａｔｅａｎａｒｔｉｆｉｃｉａｌｇｅｏｍｏｒｐｈｏｌｏｇｙｂａｎｄａｒｏｕｎｄｔｈｅｔｗｏｌａｇｏｏｎｓａｆｔｅｒｔｈｅｃｏａｓｔｔｉｄａｌｓｈｏａｌｓｒｅｃｌｉｍａｔｅｄａｎｄｔｈｅｓｅｄｉｍｅｎｔａｒｙｐｌａｉｎｓｄｒｅｄｇｅｄ，ｂｕｔｓｈｉｆｔｉｎｇｗｉｔｈｉｎｔｈｅｍｓｅｌｖｅｓ．Ｔｈｅｏｒｉｇｉｎａｌｇｅｏｍｏｒｐｈｏｌｏｇｙｐａｔｔｅｒｎｓｈａｖｅｂｅｅｎｄｅｍａｇｅｄ；ｔｈｅｗａｙｓｏｆｓｅｄｉｍｅｎｔａｒｙｄｙｎａｍｉｃａｎｄｅｎｅｒｇｙｔｒａｎｓｆｏｒｉｎｇｈａｖｅｂｅｅｎｃｈａｎｇｅｄｔｈｏｒｏｕｇｈｌｙ．Ｔｈｅｉｎｃｒｅａｓｉｎｇａｎｔｈｒｏｐｏｌｏｇｉｃａｌａｃｔｉｖｉｔｉｅｓａｒｅｂｅｃｏｍｉｎｇａｍａｉｎｄｒｉｖｉｎｇｆｏｒｃｅｏｆｔｈｅｃｏａｓｔａｌｅｖｏｌｕｔｉｏｎ．Ｋｅｙｗｏｒｄｓｓａｎｄ－ｂａｒｒｉｅｒｌａｇｏｏｎ，ｇｅｏｍｏｒｐｈｏｌｏｇｙｐａｔｔｅｒｎｓ，ｒｅｍｏｔｅｓｅｎｓｉｎｇｉｎｔｅｒｐｒｅｔａｔｉｏｎ，ｇｅｏｍｏｒｐｈｏｌｏｇｉｃｔｒａｎｓｆｅｒｍａｔｒｉｘ，ＬｉｎｇｓｈｕｉＣｏｕｎｔｙｉｎＨａｉｎａｎＩｓｌａｎｄ万方数据