海岸是陆海直接作用的地带，能量、物质输运、转化极其活跃，自然环境系统极其敏感和脆弱，极易受风暴潮和海平面上升等影响，导致海岸侵蚀及海水入侵(Nicholls et al，2006; Splinter et al，2011)。最近几十年来，日益强烈的土地利用、水库修建等人类活动，导致河流入海输沙及径流显著减少，海岸蚀退在全球各地普遍发生，且程度有日益提高之趋势(Yang et al，2003，2011a，b; Ericson，2006; Dallas et al，2011;Gao et al，2011)。海岸侵蚀除导致土地资源、人类财产直接损失外，还引发严峻的环境生态问题，制约、威胁着海岸带可持续发展，成为影响海洋经济发展的重要不确定性因素。

目前，海岸侵蚀或淤积放缓及其环境效应研究多集中于大河三角洲、平原海岸(Fanos，1995; Yang et al，2003; 王兆印等，2005，2006; Ericson et al，2006)，但对广布的山地海岸关注和研究较少。从引发海岸蚀退及环境变迁的人为驱动力看，海水养殖产业的海岸环境效应近年来受到国内外学者重视和研究(Anh et al，2010; e Silva et al，2010; Berlanga-Robles et al，2011; Brojo et al，2011; Swapan et al，2011; 战超等，2013)。但是，关于潮上带海水养殖对海岸动力地貌过程，尤其是地貌冲淤变化影响的研究还少见报道。莱州湾东部为典型的温带季风气候区砂岸，是我国工厂化海水养殖业最发达地区之一(仲少云等，2015)，同时也是我国海岸侵蚀发生暨研究较早(庄振业等，1989)，而且目前蚀退仍很严重的岸段。因此，本文以莱州湾东岸为例，研究潮上带海水养殖设施建设对海岸地貌冲淤和蚀退的影响。

1 区域背景

莱州湾系公元1855 年以来形成的弱潮海湾，为现代黄河三角洲将古渤海湾南部分隔而形成。海湾东部濒临胶东半岛西北部山地，沿其山前堆积平原发育砂质海岸，总体呈NE-SW 展布(图 1)。全新世中期该岸段为沙坝—泻湖海岸，坝外有屺坶岛、石虎嘴等基岩小岛，后经岛屿连陆、泻湖淤塞而成现代海岸(庄振业等，1994)。此砂岸系由NE-SW 泥沙流塑造而成(蔡爱智，1980)，而该泥沙流又系半岛西北部界河、王河等山地河流输沙、NE和NNE 向波浪与NE-SW 向海岸相互作用的结果(王琦等，1982)。

20 世纪50 年代末至60 年代初，胶东半岛西北部山地修建了大量水库，界河、王河等河流中上游流域基本被水库控制(王庆等，2003)。至20 世纪80 年代，开始出现该段海岸侵蚀后退的报道，是我国较早发现并研究的蚀退砂岸(庄振业等，1989)。有研究认为，水库拦截导致河流入海输沙减少是该段海岸蚀退的主要原因(庄振业等，1989; 常瑞芳等，1993)。近期野外调查发现，从1990s 年代初起该岸段蚀退范围扩大，强度有越来越大的趋势。

莱州湾东部属不规则半日潮，平均潮差1.0m;强浪和常浪向NE，平均波高1.06m，频率22%。该海域为风暴潮多发区，1952 年12 月到2013 年5 月有21次灾害性风暴潮，平均不到3 年1 次，羊角沟站记录的最高潮位达6.74m。在风暴潮和其它波浪交替作用下，海岸线(大潮平均高潮线)与风暴潮时高潮线之间的潮上带宽阔，20 世纪80 年代前为宽250—450m 的平坦风沙地，风暴潮时可被海水淹没(图 1)。80 年代初起利用潮上带修建养殖棚(池、塘)开展海水养殖，外侧建有防浪堤。

调查表明，从20 世纪90 年代初以来，莱州湾东部潮上带风沙地大部为养殖设施覆盖，形成了平行岸线、宽达数百米，基本连续展布的养殖设施带(仲少云等，2015)。土地覆被的高强度变化，必然导致风暴潮消能空间及边界条件变化。但是，关于其海岸地貌效应及其对海岸侵蚀的影响尚缺乏专门研究，而其结果有助于解释最近二十年随着水库修建时间的日渐久远，该段海岸侵蚀不仅未因动力地貌趋于平衡而减弱，反而越来越严重的地貌事实。

2 研究方法

在野外海岸地貌、沉积和蚀退状况调查的基础上，以界河口至石虎嘴岸段为典型岸段，运用水动力数值模拟方法，针对养殖设施带所占潮上带宽度的比例0%(即未利用)、30%和70%三种情境，分别对潮上带土地利用/土地覆被变化后的海岸沉积动力进行模拟计算，分析土地利用对海岸带波浪、潮流等的影响，揭示潮上带土地利用/土地覆被变化对海岸动力地貌、冲淤状态和侵蚀后退的影响。

2.1 模型搭建

根据该岸段水深、动力、沉积和地貌特征，基于二维浅水方程建立潮流场模型，采用丹麦水力学研究所(DHI)开发的MIKE21 软件进行计算。所用动量方程和连续方程分别如下:

式中: t 为时间; h 为水深; x、y 为笛卡尔坐标系中位移;为自由水面高度; $\bar{u}$，$\bar{v}$ 为流速在 x、y 方向分量; s 为源汇项排放量; p_a 为当地大气压; f 为科氏力系数(f = 2$\Omega \sin \varphi $，$\Omega $ 为地球自转角速度，$\varphi $ 为地理纬度);$\rho $ 为水密度; ${{\rho }_{0}}$ 为标准水密度; u_s 、v_s为源汇项流速在 x、y 方向分量; s_xx 、s_xy 、s_yx、s_yy为辐射应力分量; T_xx、T_xy、T_yx、T_yy为水平黏滞力项; ${{\tau }_{bx}}$、${{\tau }_{by}}$为底摩擦力在 x、y 方向分量; ${{\tau }_{sx}}$、${{\tau }_{sy}}$ 为表面摩擦力在x、y 方向分量; g 为重力加速度。

基于波作用守恒方程建立波浪场模型，如方程(4)所示:

式中: $\frac{\partial A}{\partial t}$为波能; ${{\nabla }_{x}}$ (c_g A) 为波能随地理空间的变化; ${{\nabla }_{\theta }}\left( {{c}_{\theta }}A \right)$为波能随方向的变化; ${{\nabla }_{\sigma }}\left( {{c}_{\sigma }}A \right)$为波能随频率的变化; S 为源函数。

基于对流扩散方程建立泥沙输运数值模型，如方程(5)所示:

式中: u，v 分别为速度在x，y 方向分量; $\bar{C}$ 为泥沙浓度;D_x，D_y 为扩散系数在x，y 方向分量; S 为对应泥沙沉降和悬起的源汇项; n 表示源汇总数; i 为源汇数，其余量的含义与上述相同。

根据莱州湾及渤海的潮汐及潮流场基本特征，模 型 计 算 海 域 的 范 围 确 定 为 点 Ⅰ(120°1.848′E，37°31.452′N)、点Ⅱ(120°10.068′E，37°37.050′N)、点Ⅲ(120°14.622′E，37°32.916′N)、 点 Ⅳ(120°6.006′E，37°27.258′N)连线所构成的矩形区域(图 2)，计算海域的初始流速和潮位均取值为0，东边界、西边界和北边界三个开边界(图 2)的水位通过潮位预报方法获得，计算方程如公式(6)所示:

式中: A₀为平均海面; F_i和(v₀ + u)_i为天文要素; H_i和g_i为调和常数。其中，调和常数选用 4 个分潮，包括O1和K1 两个日分潮、M2和S2 两个半日分潮。

计算海域的东边界设为外海来浪方向(图 2)。根据2003 年10 月10—12 日莱州湾“10.11”特大温带风暴潮资料，有效波高为4.5 m，波向为NE。水边线位置随潮位涨落而迁移，其最南位置为特大潮高潮线即潮上带南界。因此，模型考虑了动边界内网格节点的干湿变化，养殖设施带前缘假定为直墙边界。根据1985 年测量结果编绘的1∶10000 地形图，计算岸段潮上带南边界大致与岸线平行，其宽度统一确定为450 m。

计算海域的水深按1985 年测量的1∶250000 海图数据确定，西、北和东边界的开边界处水位根据点Ⅰ、Ⅱ、Ⅲ、Ⅳ的潮汐预报水位经线性插值得到(图 2a)。采用三角形网格法对计算海域进行单元划分，其中，外海区域的网格空间步长约为700 m，近岸区域的步长约为25 m(图 2b)。

根据潮上带养殖始于潮上带向陆侧，然后自陆向海展宽、平行海岸延伸的实际情况(仲少云等，2015)，通过调整潮上带南界至养殖设施带北侧前缘距离，即养殖设施覆盖占潮上带宽度的比例即潮上带土地利用程度，设定了三种不同的情境分别模拟计算: 潮上带未利用(即土地利用0%)、150 m 宽养殖设施带(即土地利用30%)和300 m 宽养殖设施带(即土地利用70%)。

2.2 模型验证

由于计算海域内没有水文测站，采用距离最近的龙口站1984 年9 月24 日4 时至25 日4 时的潮位资料进行验证。结果显示，该海域属于不规则半日潮，一日潮位变化包括两个涨潮和落潮过程，存在显著的高低潮不等现象(图 3)。根据验证结果，数值模拟计算潮位与实测潮位吻合较好，模拟结果符合动力地貌演变研究的精度要求。据此确定潮流场模型的涡扩散系数(Cs)为0.28，曼宁糙率系数(n)为32，干湿动边界相关的干水深为0.005 m，淹没水深为0.05 m，湿水深为0.1 m; 科氏力设置为随纬度不同而不同。

3 研究结果3.1 对波高的影响

潮上带养殖设施覆盖对近岸波高影响显著。对于选取的剖面1-1、2-2和3-3，在风暴潮作用12h和24h后，养殖设施占潮上带比例30%和70%时近岸波高较0%(即没有修建养殖池)时均显著增大，且越靠近岸线位置，增大幅度越大。波高随着潮上带土地利用程度的增强呈现非单调变化的趋势，养殖设施覆盖比例30%时最大。在不同的土地利用/土地覆被情境下，三个剖面波高的沿程波动性差异显著，但总体上自海向陆波动均减小。在岸线两侧存在自海向陆先增大、岸线附近最大、再减小的趋势。波浪自外海向近岸传播过程中，受底部摩阻等导致的能量损耗引起波高沿程衰减，在接近岸线时因发生破碎而迅速减小，在岸线附近完全破碎(图 4)。

由于模型所设置的直墙边界位置不同，其所产生的反射波与入射波相互叠加后的结果也有所不同，从而导致不同土地利用程度下波高大小的不同。考虑到海底地形的粗糙程度不同，不同剖面的底摩阻均不同，波浪在不同剖面传播过程中损耗的能量也有差别，因此不同剖面波高的沿程衰减变化也不同。再者，由于土地利用改变了海岸动力边界条件，致使同一剖面不同土地利用程度下的波浪沿程变化特征也不尽相同。此外，风暴潮作用时间不同，也会影响近岸波高大小和分布。

3.2 对水流的影响

受NE 向波浪入射方位的影响，在东边界附近海域的水流流速比其他位置处的流速大。涨急时刻，水流主流向在近岸海域为NE 向，在靠近北边界海域转为NW 向，在靠近西边界海域则转为NNE—E 向。落急时刻，西边界附近海域水流流向为SW 向，东边界附近由NE 向转为EES 向，共同流出北边界。潮上带养殖设施修建后，莱州湾东部海域的流场特征发生了显著改变。随着养殖设施覆盖比例的增大，无论落急时刻还是涨急时刻，东边界附近流速矢量分布更加密集，流速显著增大。靠近养殖设施带附近海域流速呈递增趋势，流向自NE 向SW 方向。此外，涨急时刻西边界海域水流主流向，由初始未建养殖设施时的NE 向逐渐转为E 向(图 5)。

潮上带养殖设施覆盖比例的提高显著增大了近岸海域流速。不同土地利用/土地覆被情境下，剖面1-1，2-2和3-3 的流速总体上自海向陆呈现先减小—后增大—再减小的变化趋势。在岸线两侧，自海向陆先剧增、再剧减的趋势更加显著，岸线附近流速最大(图 6)。

对于剖面1-1，养殖设施覆盖比例为30%时，除距岸约0.8—1.7 km 范围外，计算海域流速较未建养殖设施(即土地利用0%)时显著增大，岸线附近流速由 0.29 m/s 增至0.37 m/s; 覆盖比例70%时，除距岸约 5.0 km 之外范围外，计算海域流速较覆盖比例30%时显著增大，岸线附近流速达0.41 m/s(图 6a和图 6b)。

对于剖面2-2，养殖设施覆盖比例为30%时，潮上带和距岸约3.0 km 范围内流速较未建养殖设施(即土地利用0%)时显著增大，最大流速由0.26 m/s 增至 0.34 m/s; 覆盖比例70%时，整个计算范围内流速较覆盖比例30%时显著增大，岸线附近流速达到0.41 m/s(图 6c和图 6d)。

对于剖面3-3，养殖设施覆盖比例为30%时，潮上带和距岸约0.60 km 范围内流速较未建养殖设施(即土地利用0%)时显著增大，而且越靠近岸线增幅越大，岸线附近流速由0.12 m/s 增至0.30 m/s; 覆盖比例70%时，岸线附近流速较覆盖比例30%时显著增大，达0.36 m/s(图 6e和图 6f)。

3.3 海岸地貌冲淤及其分布

潮上带土地利用对海岸地貌冲淤及其分布影响显著。总体上来看，随着潮上带养殖设施覆盖比例的提高，海岸地貌冲淤强度和分布范围均显著增大。侵蚀区域不断向外海扩展，向下游延伸，最远至10 m等深线附近。其中，岸线附近侵蚀最为显著，形成NE-SW 向侵蚀条带，岸线显著后退(图 7)。

风暴潮作用12h 后，潮上带未建养殖设施(即土地利用0%)时，东北岸段较西南岸段侵蚀显著，在 5.0 m 水深以外有NE-SW 向侵蚀条带，侵蚀深度 0.04—0.08 m，主要在8.0 m等深线附近及以浅，岸线附近局部侵蚀显著; 而西南岸段侵蚀深度多在0.04 m以下，岸线附近及潮上带侵蚀深度0.04—0.08 m(图 7a)。养殖设施覆盖比例30%时，东北岸段地貌冲淤发生调整，8.0 m等深线附近侵蚀带向外海扩展，局部达10.0 m等深线，5.0 m 水深以外NE-SW 向侵蚀条带向外海扩展，侵蚀深度0.06—1.0 m; 而西南岸段侵蚀深度在0.04 m 左右，岸线附近及潮上带侵蚀深度达0.06—1.0 m(图 7c)。养殖设施覆盖比例70%时，东北岸段侵蚀范围显著增大，8.0 m等深线附近侵蚀条带局部到达10.0 m等深线以外，5.0 m 水深以外NE-SW向侵蚀条带也向外海延展，侵蚀深度0.08—1.50 m;而西南岸段侵蚀深度多在0.04—0.08 m，岸线附近及潮上带侵蚀深度达0.08—1.50 m(图 7e)。

风暴潮作用24h 后，潮上带未建养殖设施(即土地利用0%)时，5.0 m水深以外NE-SW向侵蚀条带有三分之一展宽至10.0 m等深线，侵蚀深度0.04—0.08 m; 在东北岸段5.0 m 水深到岸线附近，新形成了NE-SW向平行岸线的侵蚀条带，侵蚀深度达1.5—2.5m(图 7b)。养殖设施覆盖比例30%时，海岸地貌冲淤调整显著，5.0 m 水深以外NE-SW 向侵蚀条带有三分之二展宽至10.0 m等深线外，侵蚀深度达1.0—1.5 m; 东北段岸线附近侵蚀条带进一步向外海和陆侧展宽，侵蚀深度达2.0—3.0 m(图 7d)。养殖设施覆盖比例70%时，5.0 m 水深以外NE-SW向侵蚀条带整体延至10.0 m等深线附近，侵蚀深度达1.5—2.0m; 岸线附近侵蚀条带继续向西南延伸、向外海和陆侧展宽，侵蚀深度达到2.5—3.0 m，局部3.0 m 以上(图 7f)。

3.4 海岸横剖面变化

潮上带土地利用对海岸横剖面冲淤变化影响显著。随着潮上带养殖设施覆盖比例的提高，各个剖面的侵蚀深度均随之增强，侵蚀范围也随之增大，闭合深度则随之变深至约10.0 m，闭合位置外推至距岸约5.0—6.0 km(图 8)。

风暴潮作用12h 后，对于剖面1-1，潮上带未建养殖设施(即土地利用0%)时，岸线附近侵蚀明显，养殖设施带前缘侵蚀微小; 养殖设施覆盖比例30%和70%时侵蚀强度显著增大，集中在距岸约4.0 km范围内(图 8a)，养殖设施带前缘侵蚀略有增强。对于剖面2-2，未建养殖设施(即土地利用0%)时，岸线附近呈微侵蚀; 覆盖比例30%时，岸线附近仍呈微侵蚀，但比未建养殖设施时略有增强; 覆盖比例为70%时，岸线附近侵蚀较显著，集中在距岸约0.60 km 范围内;不同覆盖比例时养殖设施带前缘侵蚀均很微弱(图 8c)。对于剖面3-3，未建养殖设施(即土地利用0%)时，基本处于冲淤平衡状态; 覆盖比例30%时，岸线附近和养殖设施带前缘出现微弱侵蚀，其余仍处于冲淤平衡状态; 覆盖比例70%时，岸线附近和养殖设施带前缘侵蚀较显著(图 8e)。

风暴潮作用24h 后，对于剖面1-1，潮上带未建养殖设施时，岸线附近和距岸约1.0 km 内水下岸坡侵蚀显著，养殖设施带前缘侵蚀非常微弱。覆盖比例30%时，岸线附近和养殖设施带前缘侵蚀均明显增强，距岸约4.0 km 范围内水下岸坡侵蚀增强。覆盖比例70%时，岸线附近、养殖设施带前缘和距岸约4.0 km 范围内水下岸坡侵蚀均显著(图 8b)。对于剖面2-2，未建养殖设施时，岸线附近和养殖设施带前缘侵蚀显著，范围集中在距岸约0.60 km 内。覆盖比例30%时，岸线附近和养殖设施带前缘侵蚀略有增强，范围集中在距岸约0.60 km 内。覆盖比例70%时，岸线附近和养殖设施带前缘侵蚀显著，距岸约6.0 km 范围内水下岸坡出现轻微侵蚀(图 8d)。对于剖面3-3，未建养殖设施时，处于冲淤平衡状态，岸线附近和养殖设施带前缘微弱侵蚀。覆盖比例30%时，岸线附近和养殖设施带前缘侵蚀略有增强。覆盖比例70%时，岸线附近和养殖设施带前缘侵蚀显著，距岸约4.0 km范围内水下岸坡出现轻微侵蚀(图 8f)。

4 讨论

从20 世纪90 年代以来，莱州湾东岸潮上带风沙地大部分已为工厂化海水养殖棚(池、塘)及相关设施所覆盖，在海岸线与防护林之间形成了一条平行于岸线，并基本连续展布的养殖设施带(仲少云等，2015)。野外调查结果显示，随着养殖设施带的修建，最近20 多年来莱州湾东岸蚀退范围扩展、程度加重，与本文水动力数模计算结果相吻合，由此可见，本文所建立的数值模型是切实可行的。

如上述，三个剖面的流速总体上自海向陆先减小—后增大—再减小，岸线附近的流速最大。潮上带未建养殖设施时，三个剖面的最大流速均不超过0.29m/s(图 6)。据有关研究，该流速仅能起动部分细砂(窦国仁，1960，1999; 武汉水利电力学院河流动力学及河道整治教研组等，1961; 唐存本，1964)。因此，除岸线附近的侵蚀相对显著外，仅在剖面1-1和2-2距岸1km 范围内有微弱侵蚀，剖面3-3 则基本处于冲淤平衡状态(图 8)。此外，计算区域东北岸段的岸线附近，形成了NE-SW 向平行岸线的侵蚀条带，侵蚀强度较其它区域显著(图 7)。由此可见，流速最大的岸线附近也是地貌侵蚀发生最显著的区域。

随着养殖设施持续扩建，覆盖比例达70%时，剖面1-1和2-2 岸线附近和距岸5km 范围内水下岸坡流速显著增大，最大值0.41m/s，剖面1-1 增大更显著。剖面3-3 岸线附近流速较其它海域增大显著，最大值0.36m/s(图 6)。除砾石和部分粗砂外，上述流速几乎能起动所有粒级的泥沙(窦国仁，1960，1999; 张瑞瑾等，1961; 唐存本，1964)。因此，剖面1-1和2-2 岸线附近和距岸5km 范围内水下岸坡侵蚀均显著，剖面1-1 更为显著，剖面3-3 岸线附近侵蚀显著(图 8)。此外，在计算区域东北岸段，岸线附近平行海岸的侵蚀条带向下游延伸至西边界附近，10m等深线附近及其以浅侵蚀深度也显著增大(图 7)。由此可见，侵蚀显著部位与流速增大部位对应，且流速增大越显著处侵蚀也越显著。

计算海域属于宽浅型区域，水平尺度远大于垂直尺度。因此，可将实际三维浪流运动简化为二维模型，但在预测地貌冲淤演变时略显不足，无法同时直观地呈现剖面地貌形态和平面岸线的进退。此外，本研究仅考虑了一种风暴工况，在今后研究中可考虑多种工况来全面深入地探讨潮上带土地利用的地貌冲淤效应。

5 结论

自从20 世纪80 年代中期开始，在潮上带大规模地修建海水养殖设施，到90 年代初期在岸线与防护林之间形成了基本连续的养殖设施带，深刻地改变了海岸水动力边界条件，致使风暴潮时岸线附近波高和流速变化显著。其中，波高沿横剖面自海向陆波动减小，岸线两侧先增大、岸线附近达到最大后再减小; 流速自海向陆先减小后增大再减小，岸线两侧先剧增、再剧减。潮上带养殖设施修建显著改变了原有的地貌演变趋势，侵蚀强度和范围均显著增大。随着养殖设施覆盖比例的提高，在平面分布上，沿岸侵蚀条带向下游延伸，向外海展宽; 在剖面分布上，岸线附近的侵蚀深度和范围均明显增大，水下岸坡闭合水深加大，闭合位置向外海推移。因此，20 世纪80 年代中期以来的潮上带土地利用对海岸边界改变，导致风暴潮能量强度及分布发生深刻变化，是本地区最近二十年来海岸蚀退的重要原因。