山东青岛灵山湾南部海滩的时空演变及其影响因素*

吴 闯 张晓东 许淑梅 胡日军 姜胜辉 杨作升

1 海底科学与探测技术教育部重点实验室，中国海洋大学，山东青岛 2661002 中国海洋大学海洋地球科学学院，山东青岛 266100

1 概述

海滩是陆地和海洋之间的过渡地带，具有重要的社会、经济和生态价值(李广雪等，2013；Ranasinghe，2016；Vousdoukasetal.， 2020)。海滩不断调整其位置和形态以适应变化的天气、气候、海洋动力、海平面、陆源输入、地壳升降和人类活动等(Bruun，1988；陈子燊，2000；蔡锋等，2005；Jackson and Nordstrom，2020)。近几十年来，在相对海平面持续快速上升以及河流入海泥沙锐减背景下，全球海滩普遍发生侵蚀(Luijendijketal.， 2018；Vousdoukasetal.， 2020)；在中国，高强度的沿海开发导致海滩侵蚀问题更加突出(蔡锋等，2008；陈子燊等，2010；侯西勇等，2016；Zhangetal.， 2021)。

全面准确评估海滩的时空演变，是海滩侵蚀防护的前提。海滩实地测量数据可涵盖整个海滩，数据精度也高，是研究海滩时空演变的基本手段；但实测海滩费时费力，长期、连续的实测数据仅局限在为数不多的几个海滩(Turneretal.， 2016；黄祖明等，2021)。卫星图像时空范围广且有大量免费资源。使用卫星图像评估海滩的时空演变，对于缺乏长期连续实测数据的海滩的侵蚀防护具有重要意义(Zhangetal.， 2021)。

已有基于卫星图像的研究普遍选用几幅至几十幅卫星图像，先在图像上识别海滩岸线，然后使用DSAS软件(Thieleretal.， 2009)或其他软件获取海滩岸线在研究剖面上的位置变化，岸线提取过程相对繁琐以致于采用的卫星图像数量有限。使用卫星图像研究海滩演变，研究结果的精度除受制于卫星图像的分辨率外，还受海滩剖面形态和历史水位等诸多因素的制约，使用数量有限的卫星图像研究宽度和变化速率均相对较小的海滩的时空演变在准确性方面存在不足(Pardo-Pascualetal.， 2012；Zhangetal.， 2021)。

Zhang 等(2018，2019，2020)和张晓东等(2021)采用聚焦剖面方法直接提取海岸线在研究剖面上的位置，并开发了聚焦剖面的计算机辅助海岸线位置识别软件(CASPRS)。聚焦剖面避免了在卫星图像上解译和校正所有海岸线，有利于对大量卫星图像的处理；计算机辅助识别极大地放宽了卫星图像的筛选标准并确保了识别结果的正确性，研究结果的时间分辨率、精度和可靠性大幅提高(张晓东等，2021)。Zhang等(2021)进一步使用亚像素海岸线识别技术(Pardo-Pascualetal.， 2012，2018；Vosetal.， 2019)改进了CASPRS软件，并基于880幅Landsat图像综合使用多种海滩岸线指标研究了海阳海滩的时空演变，研究结果揭示了河流输沙减少、岸上工程和离岸工程等对海滩演变的影响和机制。

本研究综合使用692幅Landsat和188幅Sentinel图像，共880幅卫星图像，研究青岛灵山湾南部海滩(简称“灵南海滩”)在1984—2021年的演变，基于海滩实测结果评估利用卫星图像得出的海滩岸线的误差，目的在于深化和推进聚焦剖面、使用海量卫星图像评估海滩时空演变的研究，认识灵南海滩对自然变化和人类活动的响应。

2 研究区概况

灵山湾位于山东半岛东南沿岸、青岛市黄岛区境内。灵山湾沿岸发育多级地貌，高程50m以上为低山丘陵，随高程降低依次发育剥蚀面、基座阶地、冲积—洪积平原、海积平原、潮间浅滩和浪蚀平台(崔承琦和李家丰，1987)。研究海滩位于灵山湾南部。灵南海滩长约10km，向陆微凹，两侧有岬角保护(图1)。滩面沉积物为细砂和中砂(平均粒径1.2～2.6Φ)，分选好至较好(分选系数0.38～0.83)(徐方建等，2014)，是青岛市质量较好的海滩(王永红等，2017)。根据海滩地貌特征，崔承琦和李家丰(1987)认为灵南海滩的泥沙存在向南和向北2个输运方向，优势方向为自北向南。徐方建等(2014)和赵永芳等(2016)基于海滩沉积物粒度的研究结果佐证了上述观点。

底图描绘自2021/3/12的Sentinel图像，潮位相对于平均海平面为-2.03m，朝阳山嘴的形态描绘自早期的Landsat图像图 1 山东青岛灵南海滩及研究剖面位置Fig.1 Location of Lingnan Beach in Qingdao and study transects in Shandong Province

风河又名王戈庄河，是本区最大的入海河流，长35km，流域面积303km2，1976—1982年的年均输沙量为9×104t；本区其他河流还有两河、大卢河(又名隐珠河)和青草河等(崔承琦和李家丰，1987)。在早期(1980s前)，大量河流入海泥沙的沿岸输运导致灵南海滩向海淤进，淤进速率可达10～15m/a以上(崔承琦和李家丰，1987)。但近几十年来，由于大量水坝的建设，山东半岛河流入海泥沙几近断绝，加之海滩采沙和相对海平面上升的影响，山东半岛滨海沙滩普遍侵蚀(庄振业等，1989)。由于建坝，包括修建于1957—1967年的铁山大坝以及后期陆续在河道上修建的众多拦河坝，本区最大河流——风河近年来的入海泥沙量基本为零，风河河口段反而因为海相来沙发生淤积(匡翠萍等，2019)。杨鸣等(2005)的研究结果表明，灵南海滩在2000年前后的平均侵蚀速率约为7m/a。

本区潮汐属正规半日潮，平均高、低潮位分别为1.38、-1.40m。近岸潮流为顺岸往复流，涨潮流向SW、落潮流向NE。本区常浪向为S向，频率为48%；强浪向为SE向，最大波高7.6m。9216号台风(Polly)导致青岛沿岸在1992/8/31—9/1日出现风暴潮，致使胶南(灵南海滩当时属胶南)部分岸线后退12～14m(王文海等，1994；张晓慧等，2006)。

灵南海滩北邻1984年成立的国家级青岛经济技术开发区。在2014年，经国务院批复，青岛经济技术开发区的范围扩大至黄岛区全境，灵南海滩目前在开发区境内。近10年来，灵南海滩周边发展迅速。城市阳台和海军公园建于2012年，海滩北侧和南侧的堤坝分别建于2013和2016年，海滩北侧岬角——朝阳山嘴周边海域在2014年被圈围，海滩北侧的大型离岸人工岛——星光岛建于2014年。

3 研究资料与方法

3.1 卫星图像来源和预处理

本研究使用谷歌地球引擎(GoogleEarthEngine，https: //earthengine.google.com/)云平台下载灵南海滩1984—2021年云量小于80%的所有Landsat和Sentinel图像。人工删除50%海滩岸线不可见的图像，剩余880幅图像(图 2)。Landsat和Sentinel图像的分辨率分别为30m和10m。为方便后续处理，Landsat图像的分辨率采用双三次插值方法提升至10m。卫星图像利用近红外、红和绿波段合成，合成图像根据每个波段数值的频率分布进行了增强。在合成图像上，植被为红色，湿沙滩为褐色，水体为蓝色，干沙滩、浪花和云为灰白色(图 3)。

A—卫星图像的时间和类型；B—卫星图像每年的数量图 2 山东青岛灵南海滩1984—2021年的Landsat和Sentinel图像信息Fig.2 Information of the Landsat and Sentinel images of Ling-nan Beach in Qingdao of Shandong Province from 1984 to 2021

A—2021/3/26日的Landsat图像；B—2021/3/22日的Sentinel图像图 3 计算机辅助海岸线位置识别软件(CASPRS)界面以及干湿线和水边线在山东青岛灵南海滩剖面上的位置Fig.3 Operation interface of computer-aided shoreline position recognition software(CASPRS)and positions of dry/wet line and waterline on study transect in Lingnan Beach in Qingdao of Shandong Province

3.2 水边线和干湿线位置的提取

水边线是海滩和海水之间的瞬时分界线，干湿线则是上次高潮时海水在海滩上留下的痕迹线(Boak and Turner，2006；Pagnetal.， 2016；Zhangetal.， 2021)。本研究设置6条剖面T1～T6(图 1)并使用CASPRS软件提取水边线和干湿线在剖面上的位置(图 3)。CASPRS软件以亚像素的分辨率通过自动和交互2种方式识别海滩岸线在研究剖面上的位置，自动识别的分辨率为1m，交互识别的分辨率取决于图像的放大倍数。CASPRS软件可对图像进行1～30倍的缩放，在放大倍数为10时，交互识别的分辨率为1m。先用该软件自动提取水边线和干湿线在研究剖面上的位置，然后逐一检查和更正。

3.3 海滩剖面形态反演和坡度计算

水边线在海滩剖面上的位置受瞬时潮位控制。当潮位较高时，海滩出露较少，水边线位置靠近剖面起点，潮位较低时相反。干湿线在剖面上的位置主要受历史高潮水位的影响，变化规律与水边线一致。本研究采用潮汐模拟模型NAO.99b(Matsumotoetal.， 2000)计算本区的历史潮位，计算结果为逐时潮位。Zhang等(2018)研究表明，NAO.99b 模拟结果不包含潮位季节变化信息，因此，参照Zhang等(2018)的方法，本研究采用青岛港潮高的季节改正值(国家信息中心，2019)对模拟潮位进行校正。卫星图像拍摄时(精确到秒)的瞬时潮位根据逐时模拟潮位采用线性插值方法得出，卫星图像拍摄前的大潮高潮位采用模拟潮位在卫星图像拍摄前25h内的最大值。本研究还收集了2011/10/1—31日丁家咀码头(灵南海滩东北6km)的实测潮位，并对模拟潮位进行了误差评估。

3.4 平均高潮线和低潮线的计算方法

虽然干湿线位置受上次高潮位影响，但其变化幅度相对较小，再加上本研究采用的卫星图像较多，因此，本研究直接使用干湿线指标指代海滩的高潮线，并用其平均值或拟合值指示平均高潮线。受潮位波动的影响，水边线位置的变幅较大，本研究使用分段潮位校正方法计算低潮线的位置，进而使用其平均值或拟合值指示平均低潮线(Zhangetal.， 2021)。参照Zhang 等(2021)并根据灵南海滩的剖面形态，本研究选用历史潮位在-1.4～0.5m的水边线位置，根据公式(1)计算低潮线的位置：

P低潮线=P水边线-(平均低潮位-瞬时潮位)÷tan(坡度)

(1)

其中P低潮线和P水边线分别为低潮线和水边线在剖面上的位置，本区的平均低潮位是-1.4m，坡度的单位是度。

3.5 误差评估方法

在海滩岸线识别过程中，发现Sentinel图像存在明显的系统偏移。为评估该系统偏移并消除其影响，本研究选择平行于海滩岸线的4个线状地物(道路中心和水塘边缘)，平行于研究剖面设置了4条校准剖面C1～C4，利用CASPRS软件识别2015年后上述4个线状地物在校准剖面上的位置，进而比较和校准Landsat和Sentinel图像的系统偏移。本研究利用RTK实测结果评估根据卫星图像得出的高、低潮线的系统误差，与RTK实测值比较的是上述2个指标最近1年(2020/4/29—2021/4/28)的平均值。最近1年高、低潮线的标准偏差被用来评估其随机误差。

4 结果

4.1 卫星图像中的历史潮位

本区的模拟和实测潮位基本一致，平均绝对误差为0.14m，相关系数为0.99，但在实测潮位低于-0.5m时模拟潮位偏高约0.1m(图 4)。利用模拟潮位得出的灵南海滩卫星图像拍摄前的大潮高潮位在0.51～2.35m，其均值与本区平均高潮位一致；卫星图像拍摄时的瞬时潮位在-2.47～1.80m，均值(-0.36m)显著低于平均海平面(图 5-A)。

A—2011/10/1—31日青岛丁家咀码头(35°55′N，120°8′E)的实测潮位和NAO.99b(Matsumoto et al.， 2000)模拟潮位；B—上述实测潮位和模拟潮位的交汇图；C—模拟潮位的系统偏差图 4 山东青岛灵南海滩潮位数据分析Fig.4 Data analysis of tidal height of Lingnan Beach in Qingdao of Shandong Province

A—卫星图像拍摄前的大潮高潮位和拍摄时的瞬时潮位；B—1984—2021年日内逐时平均潮位；C—青岛港潮位的季节性波动图 5 山东青岛灵南海滩历史潮位及其季节和日内波动Fig.5 Historical water levels and the seasonal and daily tidal height deviations of Lingnan Beach in Qingdao of Shandong Province

4.2 干湿线和水边线的位置

图 6为干湿线和水边线在灵南海滩T1～T6剖面上的位置。水边线位置的波动范围较大，在150～300m，干湿线位置的波动范围较小，不到30m，干湿线位于水边线靠岸一侧。在时间变化方面，干湿线位置呈阶段性向岸蚀退的趋势，靠岸侧水边线位置的变化与干湿线位置的变化一致，但靠海侧水边线位置的变化趋势不明显。在所有880幅卫星图像上，水边线和干湿线位置的平均提取率分别为91%和71%。

4.3 海滩剖面形态和坡度

根据干湿线和水边线位置的时间变化(图 6)，选择变化较小的时段，结合模拟潮位，本研究分别绘制了T1～T6剖面早期(T1～T5剖面为1984—1995年，T6剖面为1984—1991年)和晚期(2010—2021年)的形态，并选取晚期的剖面形态与RTK实测结果进行了对比(图 7)。

图 6 山东青岛灵南海滩T1～T6剖面干湿线和水边线位置Fig.6 Dry/wet line and waterline positions on transects T1～T6 of Lingnan Beach in Qingdao of Shandong Province

图 7 山东青岛灵南海滩T1～T6剖面形态和坡度Fig.7 Shape and slope of transects T1～T6 of Lingnan Beach in Qingdao of Shandong Province

灵南海滩早期的剖面形态呈两段式，分界线的潮位在T1～T4和T5～T6剖面分别为0.5m和1m。上段较陡，坡度在3.4°～17°；下段较缓，坡度在0.8°～1.1°。在晚期，原上下段分界线的潮位未变，上段坡度也基本不变；下段的形态在中部T3～T4剖面变化较小，但在南北两侧的T1～T2和T5～T6剖面却发生了明显分化。T1～T2和T5～T6剖面下段分别在-1.7m和-1.5m潮位处发生转折，本研究进一步把其分为中段和下段。在晚期，T1～T2和T5～T6剖面中段侵蚀显著，但坡度仅在T5剖面变化较大；下段的侵蚀上大下小，坡度明显变缓至0.3°～0.4°。

4.4 平均高潮线和低潮线的演变

灵南海滩平均高、低潮线在1984—2021年发生了阶段性、基本同步的变化(图 8)。T1、T2和T6剖面的平均高潮线分别在1997—2006、2004和1992—1997年发生快速侵蚀，侵蚀速率为7.2m/a、约60m/a和11.9m/a，总侵蚀量为66m、30m和71m，持续时间分别为9.2年、0.5年和6.0年；平均低潮线快速侵蚀的起始时间与平均高潮线一致，但侵蚀速率和总侵蚀量略小，侵蚀结束时间滞后1～3年。在其他时期，平均高潮线、低潮线的侵蚀速率分别在-0.3～1.3m/a和-0.7～2.0m/a，平均低潮线侵蚀速率的变幅较大。

图 8 山东青岛灵南海滩平均高、低潮线在T1～T6剖面的演变Fig.8 Position changes of mean high- and low-water lines on transects T1～T6 of Lingnan Beach in Qingdao of Shandong Province

4.5 误差评估结果

上述4个用于校准的线状地物在Landsat和Sentinel图像上的平均位置之差为0.6～2.3m，说明从谷歌地球引擎下载的Landsat和Sentinel图像的地理位置基本一致、2个图像库可对比。4个线状地物在Sentinel图像上的位置存在较强的正相关关系，相关系数为0.64～0.75，而其在Landsat图像上的位置的相关性普遍较差，说明Sentinel图像存在系统性位置偏差而Landsat图像的系统性位置偏差较小。本研究采用上述4个线状地物相对其平均位置的偏离值的均值作为校正量，对干湿线和水边线在Sentinel图像上的位置进行了校正。4.2节至4.4节中的数据均为校正后的结果。

平均高、低潮线在Sentinel和Landsat图像上的系统和随机误差见表 1。总体上，平均高潮线的误差较小，平均低潮线的误差较大。平均高、低潮线在Sentinel图像上的误差普遍小于其在Landsat图像上的误差，但平均低潮线的系统误差例外。

表 1 山东青岛灵南海滩平均高、低潮线在Sentinel和Landsat图像上的系统和随机误差(m)Table 1 Systematic and random errors(m) of mean high- and low-water lines on Sentinel and Landsat images of Lingnan Beach in Qingdao of Shandong Province

5 讨论

5.1 海滩岸线的位置误差及其影响因素

灵南海滩平均高潮线的系统误差为3.7～8.2m(表 1)，与Pardo-Pascual等(2018)得出的3.1～5.5m基本相当；平均高潮线的随机误差为6.0～7.5m(表 1)，与Pardo-Pascual等(2018)得出的4.9～6.5m和Zhang等(2021)得出的7.8m也基本相当。然而Vos等(2019)得出的海滩岸线的系统误差仅为0.5～2.0m，但其随机误差较大，为7.3～12.7m。Vos等(2019)采用神经网络识别海滩岸线，在原理上，使用无偏数据集训练神经网络可以使系统误差为0。但是，由于潮位误差和海水的透光性等诸多因素的影响，利用海滩剖面形态和瞬时潮位得出的水边线与卫星图像上的水边线并不是同一条线，而且前者在卫星图像上的解译特征相对不够明确，这可能是Vos等(2019)识别结果的随机误差偏大的原因。

灵南海滩平均低潮线的随机误差为14.8～19.8m(表 1)，显著高于Pardo-Pascual等(2018)和Vos等(2019)的上述结果，这主要是由于灵南海滩在-1.4～0.5m区间(计算平均低潮线的潮位区间)的坡度较小所致。灵南海滩在-1.4～0.5m区间的坡度为0.8°～1.3°(图 7)，相当于1︰72～1︰44，显著小于Pardo-Pascual等(2018)和Vos等(2019)所研究海滩1︰23～1︰5的坡度。海滩坡度越小，由于潮位误差和海滩剖面的不规则形态引起的海滩岸线误差越大(Zhangetal.， 2021)。根据Zhang等(2021)在山东半岛海阳海滩得出的计算公式以及T1～T6剖面在-1.4～0.5m区间的坡度，灵南海滩随机误差的计算结果为11.4～14.1m，比本研究结果略小。区别于海阳海滩，灵南海滩的剖面形态起伏较大，例如在T3剖面发育有潮间带沙坝(图 7)，以致T3剖面的平均低潮线波动较大(图 8)，其随机误差也偏大(表 1)。

灵南海滩平均低潮线的系统误差为11.4～22.7m(表 1)，相对较大。本研究的模拟潮位在实测潮位低于-0.5m时偏高约0.1m，这会导致通过潮位校正得出的灵南海滩平均低潮线偏大4～7m。由于海水的透光性，卫星反演的低潮线相对实测的平均低潮线更靠海，进一步加大了平均低潮线的系统误差。另外，本研究依据单次海滩测量进行误差评估，评估结果还受海滩剖面季节变化的影响。

海滩岸线的系统误差对计算海滩变化速率影响较小(Luijendijketal.， 2018；Zhangetal.， 2021)，但对反演剖面形态影响较大，利用卫星反演得出的海滩坡度相对实测结果普遍偏小(图 7)。上述系统误差仍需通过更多次的海滩剖面重复测量进一步评估，并根据评估结果对卫星反演结果进行校正。海滩岸线的随机误差由于其内在的随机性不能被校正，但可通过大量卫星图像的应用减弱其影响。因此，研究宽度和变化速率均相对较小、研究结果对误差较为敏感的海滩的演变，采用的卫星图像越多越好。

5.2 影响灵南海滩演变的因素

如上所述，灵南海滩的T1、T2和T6剖面在历史上发生了阶段性的快速侵蚀(图 8)。T1剖面平均高潮线在1997—2006年的快速侵蚀与海滩上部挖沙有关。海滩上部挖沙留下的沙坑与海滩有沙埂间隔，在正常海况下，海滩岸线尚能维持，但在风暴潮来临之际，狭窄的沙埂易被摧毁，造成海滩岸线显著后退(王文海等，1994)。根据T2剖面附近养殖场改建前后的卫星图像可知，T2剖面平均高潮线在2004年的快速后退应与其周边养殖场的改建有关，改建后养殖场的宽度比原来显著减小。根据1992—1997年的卫星图像可知，T6剖面平均高潮线在1992—1997年的快速侵蚀与其周边养殖场靠海侧堤坝的破坏有关： 在1992/6/30，养殖场外侧海滩的干滩已经消失，经过9216号台风(1992/8/31—9/1)，在1992/9/18养殖场外侧海滩的堤坝出现决口，决口持续扩大直至1998年才恢复。另外，9216号台风后的海滩决口并不大，但一直到5年后才恢复说明该养殖场当时已废弃。

沙埂和堤坝损坏后，平均高潮线随即开始快速侵蚀，平均低潮线的侵蚀相对滞后。T1、T2和T6剖面下段(-1.7～-1.5m以下)的侵蚀相对中段较小(图 7)，这可能与海滩滩面泥沙在冲流带(swash zone)的上冲有关。冲流带位于瞬时水位附近，是潮汐和波浪等引起的上下冲流联合作用的地带(Masselink and Puleo， 2006；Puleo and Torres-Freyermuth，2016；黎树式等，2017)，其在海滩剖面上的位置目前尚无严格界定。本研究根据本区潮位的频率分布，选择潮位分布频率大于最大频率50%的区间作为冲流带在海滩剖面上的位置(图 9)。根据该方法，灵南海滩冲流带的下界位于-1.5m左右，与T1～T2和T5～T6剖面的侵蚀情况基本一致(图 7)。平均高潮线的快速蚀退给冲流带的上冲泥沙提供了容纳空间，进而造成冲流带的侵蚀以及平均低潮线的快速蚀退。

图 9 山东青岛灵南海滩历史潮位的频率分布和冲流带Fig.9 The frequency distribution of historical tidal heights and the swash zone of Lingnan Beach in Qingdao of Shandong Province

尽管已有研究揭示了台风对山东半岛海岸有较强的破坏作用(王文海等，1994；高伟等，2020)，而且本研究也显示9216号台风导致灵南海滩在T6剖面附近溃岸，但是，9216号台风在T6剖面处直接产生的海滩侵蚀量相对较小，在灵南海滩其他剖面处也未造成明显的海滩侵蚀，且其他时期的风暴潮(张晓慧等，2006)也未对灵南海滩的侵蚀产生明显影响(图 6，图 8)，这说明导致灵南海滩快速侵蚀的主要原因是沿岸挖沙、养殖场的废弃，风暴潮仅起次要的辅助作用。

除在上述3个时段发生快速侵蚀外，灵南海滩普遍发生慢速侵蚀(图 8)，这是相对海平面上升和河流入海泥沙锐减等的必然结果。河口附近的T3剖面的平均高、低潮线以及北部的T5～T6剖面的平均低潮线侵蚀速率相对较大，在1.3～2.0m/a。T3剖面位于河口附近，该处海滩的侵蚀应与河流入海泥沙减少有关。T5～T6剖面平均低潮线的侵蚀除与周边河流入海泥沙减少有关外，还与前述的本区泥沙优势输运方向为自北向南有关。截至目前，T5～T6剖面平均高、低潮线的长期慢速侵蚀趋势未发生明显改变(图 8)，本区北部2013—2014年的围填海工程(包括星光岛)对灵南海滩长期演变的影响尚不明确。

6 结论和建议

本研究基于880幅卫星图像，利用CASPRS软件识别干湿线和水边线在研究剖面上的位置；结合模拟潮位分析海滩剖面形态并计算海滩坡度，利用干湿线得出海滩平均高潮线，通过潮位校正得出平均低潮线，采用多种指标综合研究灵南海滩在1984—2021年的时空演变；结合历史资料分析影响因素，还利用海滩实测数据对研究结果进行误差评估。大量卫星图像的应用提高了研究结果的时间分辨率、精度和可靠性。

在灵南海滩，卫星提取的平均高潮线的系统误差为3.7～8.2m，随机误差为6.0～7.5m，与目前最精确的研究结果相当；平均低潮线的系统误差为11.4～22.7m，随机误差为14.8～19.8m，相对较大。灵南海滩中下部较缓的剖面形态是其平均低潮线误差较大的主要原因。除海滩坡度外，海滩剖面不规则的形态、模拟潮位的误差、海水的透光性以及海滩剖面的季节变化等也会影响海滩岸线的误差。

灵南海滩的部分岸段发生了阶段性快速侵蚀。该快速侵蚀发生在研究区北部的T1、T2剖面以及南部的T6剖面，发生时间分别为1997—2006年、2004年和1992—1997年；平均高潮线的侵蚀速率为7.2～60m/a，总侵蚀量为30～71m，持续时间为0.5～9.2年；平均低潮线的开始侵蚀时间与平均高潮线一致，总侵蚀量大致相当，但结束时间滞后1～3年，侵蚀速率稍小。该快速侵蚀主要是海滩挖沙、养殖场改建和废弃等人为因素造成的，风暴潮起次要的辅助作用。为避免海滩的快速侵蚀，海滩挖沙应杜绝，距海滩较近的养殖场“退养还滩”后，应对养殖场进行回填。在其他时段，灵南海滩普遍发生缓慢侵蚀，侵蚀速率一般小于2m/a，这是相对海平面上升和河流入海泥沙锐减等的结果。

致谢感谢审稿专家提出了宝贵的意见和建议，感谢中国海洋大学本科生研究发展计划(OUC-SRDP)的支持，还感谢檀夏伟、王佳星、张家宝、李胤辉和张亚荣同学对本研究的贡献。与本研究有关的野外调查照片和典型卫星图像请参见OSID，这些图片以及CASPRS软件也可以通过https: //github.com/ouczxd/CASPRS访问和下载。