周良勇,薛春汀,刘 健,梁 源
(1.国土资源部 海洋油气资源和环境地质重点实验室,山东 青岛266071;2.国土资源部 青岛海洋地质研究所,山东 青岛266071)
山东半岛东、北部主要为海滩海岸、障壁岛-潟湖海岸和基岩海岸,沙质海滩发育。了解海滩的地貌动力特征对海岸工程开发(如滨海公路、浴场修建)和旅游开发等有重要意义。根据2003—2006年的海滩调查资料,分析了山东半岛东北部海滩平面形态、剖面序列及剖面的年际和季节性变化,并利用沉积物-波浪参数(Dean参数),研究了山东半岛东、北部海滩的动力地貌的类型。
海滩倾向受到断层构造等影响[3]。海滩倾向不同,受到的波浪的作用也不同,这影响着海滩动力地貌的变化。我们的研究区域,包括牟平养马岛西侧至荣成市褚岛一带的一系列海滩(图1中养马岛至威海之间的M1~M7为以往监测剖面①王文海,吴桑云,夏东兴,等.典型岸段海岸侵蚀及其对策研究.1996-08.,P1~P16为我们的监测剖面;波高、潮差来自文献[4-5])。以成山头(或称成山角)为界,养马岛西侧至成山头部分为北部海岸,成山头至褚岛为东部海岸。两侧季节性波浪条件明显不同,为研究波浪动力因素对海滩动力地貌影响的提供了良好的条件。主要海滩的平面形态及倾向见图2。
据《山东省海岸带和海涂资源调查报告集》和《中国海湾志》第三分册[4-5],研究区波浪四季都是以风成浪为主,频率为94%~100%。年均H1/10为0.2~0.7m。年均波浪周期在1.6~3.4s。北部海岸波浪特征不同于东部海岸的。北部海岸最大月均波高在冬季,夏季较小,主要波浪方向为N和NNW;而东部海岸波浪特征正相反,冬季较小,夏季较大,主要波浪方向是N,NNW,SSW和SW。极端波高出现在成山头站,波高达8.0m。研究区潮汐小,年平均大潮潮差<2m,属于弱潮海岸(表1和图1)。烟台至威海以西(高科技开发区)地区为正规半日潮;威海市以西及东部海岸为不正规半日潮。平均大潮潮差为0.80~2.08m,属弱潮海岸。这里潮差小,属浪控海岸和波浪-潮汐共同作用的海岸[2]。浪潮作用指数K[6]和相对潮汐指数RTR值[7]反应了潮汐和波浪对海滩的相对影响程度。这里利用本地烟台芝罘岛至石岛区资料计算了浪潮作用指数K和相对潮汐指数RTR值。山东半岛易受气旋、冷锋和热带气旋等天气系统影响,形成风暴潮[7]。受冷锋影响,冬季经常出现大风天气(>17m/s)。海滩沉积物为富含石英和长石的沙质,主要海岸滩沉积物来自河流和海岸花岗岩等风化产物。较大的河流有辛安河、沽河、沁水河等。
研究区海滩的平面特征受到岬角、澙湖出口等的控制。这里使用沙质海滩开敞程度(岬角影响程度)和凹陷性比值来确定控制海滩发育的的海湾特征。沙质海滩开敞程度定义为海湾宽度(CL)和海湾岸线的长度(SL)的比值[8-10],比值(SL/CL)说明海滩暴露程度。此外,海滩的暴露性也可以用海湾凹陷值(a)和岬角间隔距离(R0)比率表示[11]。海滩的这些参数采用1996年1∶10 000航空影像图像在Mapinfo软件下进行测量。
2003年和2004年进行了取样。在16个剖面(图1)和调查点(图3a)上共获得287个海滩表层(0~10 cm)沉积物样品。调查点在高潮线附近取沉积物样;剖面上,沿测线每个地形测点取一个样。用筛析法进行粒度分析。调查点沉积物粒径取其沉积物中值粒径,剖面的粒径为各样品中值粒径的平均值。
Dean参数和海滩地貌特征被用于划分海滩动力地貌类型[12-13]。Dean参数(Ω)由波浪参数和沙粒的沉降速率决定,计算公式:
其中Hb为破波带波高,T为波浪周期,ω为沉积物沉降速率。
选用5个波浪站的波浪统计数据:芝罘岛站、旗杆嘴站(威海)、青矶岛站、成山角站和桑沟湾站。多数调查点的样品是2003-08取得,对应使用夏季或8月的波浪值。调查点26和28(图4)是4月取得,对应使用附近青矶岛站4月的波浪值。剖面上的样品都是2003-10取得,因此对应使用波浪站10月或秋季的波浪值。
由于收集的波浪波高为H1/10,需要通过2步计算获得Hb。首先通过海港水文规范估算有效波高(H1/3)(式2),然后利用 Kormar and Gauhan的半经验公式(式3)[14]计算Hb:
分别向空白猕猴桃和土壤样品中添加4个水平的螺虫乙酯及其4种代谢产物的混合标准溶液,添加水平分别为0.05、1.0和2.0 mg/kg。按1.3节方法进行前处理,按1.4节条件检测,计算添加回收率及相对标准偏差(RSD)。
式(3)中,H0为深水波高;L0=(gT2)/(2π),为深水波长。
沉降速率根据Ferguson和Church的公式[15],用中值粒径进行计算:
其中C1和C2为常数,C1=18,C2=1.0;D为中值粒径(D50),g为重力加速度,R为沉积物密度(根据石英密度计算),υ为运动黏度,υ=1.00×10-6m2·s-1(水温20℃时)。
根据Dean值可对海滩类型进行划分:Ω<2时为反射型海滩,2≤Ω≤5时为过渡型海滩,Ω>5时为消散型海滩;过渡型海滩进一步可分为4种亚类:沿岸坝-凹槽型、韵律沙坝型、横向沙坝裂流型和脊-槽型[8]。
2003—2006年进行了5次剖面重复测量。剖面一般选在海滩中部,使其更具带表性。起点有固定的标志,使测量结果可对比。沿海岸选则了16个剖面进行监测。使用水准仪进行剖面测量,用罗盘确定方向,测量范围为后滨(或海岸沙丘)至低潮线。并且在剖面上取表层5~10cm沉积样品,进行了粒度分析。对剖面进行比较,分析了空间上海滩序列变化和季节性变化。
山东半岛东北部的沙质海岸按岸滩形态组合特征和成因,可分岬湾海岸、沙坝-潟湖海岸和夷直海岸。岬湾海岸分布在威海北部、成山头附近和东部荣成诸小湾地区,岬湾海岸发育袋状海滩。沙坝-潟湖海岸广泛分布在山东半岛北部、东部地区,典型的如朝阳港外的海滩、月湖海滩等。夷直海岸为养马岛至双岛湾的一段海滩。
统计的各沙质海滩开敞程度(岬角影响程度,CL/SL)和凹陷性比值(a/R0)(表2和图2)表明,两者负相关,即CL/SL越大,a/R0则越小,相关系数r2=0.86。养马岛至双岛湾口海滩是平直海岸而非弧形海岸,为开敞型,其CL/SL值最大,接近1;成山头附近马兰湾海滩(图1中P10)为袋状海滩最小,为0.17。北部海滩开敞度好于东部海滩的。桑沟湾是半开敞海湾,CL/SL为0.16,尽管局部为弧形海岸,但整个海湾是封闭性强。
表2 山东半岛东、北部沙质海滩开敞程度(CL/SL)和凹陷性比值(a/R0)Table 2 The degree of headland(CL/SL)and indentation ratio(a/R0)of sandy beaches along the northern and eastern coasts of Shandong Peninsula
山东半岛东北部海滩沉积物组分在细沙至粗沙之间,主要矿物成分是石英和长石,此外还有少量钙质贝壳碎片。利用粒度数据和与之相同季节的5个站位的波浪数据(表3)[4-5],可计算出各海滩Dean参数。图3显示在各个海滩的取样位置、海滩粒度的中值粒径和海滩的Dean参数。图3中①~④代表不同波浪站范围。成山头附近(图3b中24~30)海滩粒度比较粗。这一带有些砾石滩是>2mm的砾石构成的。其他海滩的粒度为0.25~0.80mm,并有从成山头向西和向南两翼变小的趋势。Dean参数显示出研究区海滩差异,据此大致可分为2部分:1)北部海滩在威海西面Dean值较大(>5),属于消散型海滩;2)威海东侧和东部海岸的Dean值明显较小(<2),属于过渡类型或反射类型。
表3 山东半岛东、北部海岸夏秋波浪Table 3 Wave parameters in summer and autumn along the northern and eastern coasts of Shandong Peninsula
图3 海滩取样位置、中值粒径及Dean值分布Fig.3 Sampling locations,median grain sizes and Dean numbers of the beaches studied
重复测量显示海滩的形态在空间上和季节性的变化,以及其与Dean参数关系(图4水位线根据低潮位置和潮汐表推测)。从以下5个方面考察海滩的时空变化:1)海滩坡度和宽度;2)平滑海滩/有折的海滩;3)滩面形态,如上凸、下凹或平直;4)海滩坝数量变化;5)风成沙丘。
图4 海滩地形剖面对比(位置见图1)Fig.4 Contrast of beach topographic profiles(the locations referred Fig.1)
北部夷平海滩的坡度较小、宽度较大,一般坡度在1.1°~2.1°,最大宽度180m。有些海滩剖面潮间带有明显拐点,可将剖面分成陡的向陆海滩和缓的向海海滩。如2004-10剖面P4的海滩形态,陡的一段坡度达到6.4°。北部开敞弧形海岸海滩和东部类似海滩的宽度接近,都比夷直海岸坡度大,宽度变小。封闭性强的海滩坡度最大、宽度最小。北部海岸以下凹的形态为主,而东部海岸则以平直及上凸的海滩为主,如P12和P13。这反映北部海岸受到更强的风浪的作用。
海滩沙坝是沉积物静离岸输运在碎波带形成的潮间带凸起。海滩沙坝出现与坡度、宽度关系密切。研究区整个区域的海滩显示了这种相关性,即坡度小、宽度大时易出现沙坝。北部海岸最西侧靠近辛安河附近的海滩,坡度小,沉积物供应充足,发育海岸沙丘,则海滩上形成多个直的海滩沙坝,如烟台市牟平区(图5a)。而宽度窄的封闭海滩,没有海滩沙坝。
图5 山东半岛东、北部各类海滩照片Fig.5 Photographs of different beaches along the northern and eastern coasts of Shandong Peninsula
Dean参数对北部开敞的平直海滩剖面的状态和封闭海湾海滩状态的判断较为准确。北部开敞的平直海滩剖面的Dean参数较大,大于或接近5,而其形态为消散型为主,表现为坡度小、沉积物来源丰富,其中剖面P1为典型的有多个坝的消散型(图5a、图5b和图5c)。封闭海湾海滩的Dean参数为0.3~1.4,其形态为变化较小的反射型海滩,与Dean参数判断的一致(图4和图5e)。而开敞的岬湾海滩的Dean参数为0.5~2,根据海滩状态模型,其应为反射型,而实际海滩形态则包括了反射类型(图4中P12和P13)和过渡类型,而且以有坝的过渡类型为主(图5d和图5f)。因此,对开敞的岬湾海滩海滩,Dean参数不能准确反映其状态。
海滩剖面季节性变化也比较明显(图6)。北部海岸夏季(10月)容易形成沙坝,冬季(2月)则较少,如北部的平直海岸剖面和北部的岬湾海滩剖面的变化。北部海滩冬季海滩常变低数十厘米,或者海滩靠岸的上部沙坝受侵蚀,而下部堆积,使海滩变平滑。王文海等监测的1992年冬季和1993年夏季海滩剖面(图1中M1~M4和M7),也揭示海滩沙坝夏季出现,而冬季消失。而东部海岸海滩剖面则相反,海滩冬季出现沙坝,而夏季消失①。沙坝的变化显示北部和东部海岸季节性变化的不同。
图6 山东半岛东部、北部剖面的季节性变化及海滩坝的出现/消失Fig.6 The seasonal changes of beach profiles and the occurrence and disappearance of beach bars along the northern and eastern coasts of Shandong Peninsula
威海以西基本为消散型海滩,威海至成山头主要为过渡型海滩,成山头至楮岛以反射型为主,这与Dean参数反映的海滩动力地貌一些特征的判定结果基本是一致的。当然,控制海滩发育的还有其它地质控制因素,如构造对海滩倾向的控制、沉积物供给等[3,17-18]。尽管这个模型有一定的局限性,但数据显示,利用Dean参数可对山东半岛东、北部海滩进行初步分类。相关性分析显示海滩地形参数和动力参数的关系(图7)。
图7 海滩地形参数和动力参数Fig.7 Topographical and dynamic parameters of the beaches studied
海滩潮间坡度与Ω值有一定相关性(r2=0.54),而与岸滩平均沉积物中值粒度、潮差、相对潮差之间的相关性不强。这反映出地貌与单纯的沉积物或动力因素关联不强,而与综合因素有关。K值则与海滩平均沉积物粒度有较强的相关性(r2=0.648 8),但与海滩潮间坡度无关。这与华南海滩的特征[3]有一定的区别。
双岛湾以东至研究区东南端,是波浪和潮汐共同作用的海岸[6],季节性风浪对海滩地形发育影响非常明显。
图8 山东半岛东、北部波浪站波浪玫瑰图和各月平均波高变化[1]Fig.8 Wave rose diagrams and monthly averaged wave heights from the wave stations along the northern and eastern coasts of Shandong Peninsula[1]
北部海岸带波浪特征不同于东部。整体上,北部海岸年均浪高大,而东部年均浪高小(图1)。东部海域密集的海带养殖,这有减小波浪能量的作用[19]。从浪向玫瑰图看,在北部,海岸受到N,NE,和NW风引起的波浪作用,主浪向为北;而东部则同时受到北部和南部风浪的影响,至少有2个主浪向(图8)。月平均浪高图反映:北部海岸冬季浪高大,夏季浪高小;而东部则相反,冬季浪高小,夏季浪高大。由于冬季N和NW向风暴影响,因此北部冬季波高较大。与之相反,东部海岸受控于东风和东南风引起的波浪,因此夏季波高较大。这可以解释北部和东部海滩沙坝的出现和消失。这是由于海滩沙坝易出现在较小的涌浪条件下,而在大风浪时常受到破坏。因此北部和东部海滩沙坝的季节性变化正好相反:北部夏季出现潮间沙坝,冬季消失,而东部剖面夏季剖面较平滑,而冬季易出现潮间沙坝。构造决定着海滩的倾向,而倾向与波浪动力条件有关,这些因素共同影响着海滩的变化。
根据海滩剖面形态、平面形态和Dean参数,可以把山东半岛东北部海滩分为3种主要的海滩状态类型。北部海滩以消散型和过渡型为主,而东部海滩主要为过渡型和反射型的海滩。海滩的状态主要受到水动力因素的控制,因此,Dean参数在区域范围内是适用的指标性参数。
季节性变化显示,波浪条件控制山东半岛北部和东部海滩剖面变化。两段海岸季节性波浪的不同,使北部和东部海滩剖面变化趋势也不同。潮间沙坝清楚反映了这种不同。
(References):
[1]CHEN Z C.Response characteristics of beach morphologies in an arc-shaped coast to typhoon waves[J].Chinese Science Bulletin,1995,41(12):1015-1019.
[2]WANG W J,YANG X W.State and evolution of sandy beaches along the coast of South China[J].Journal of Tropical Ocean,1996,15(4):9-15.王文介,杨雪舞.华南沿海海滩状态和演变[J].热带海洋,1996,15(4):9-15.
[3]CAI F,CAO H,SU X,et al.Analysis on morphodyanmics of sandy beaches in South China[J].Journal of Coastal Research,2007,23(1):236-246.
[4]Scientific and Technological Committee of Shandong Province(STCSP).Integrated investigation report on coastal zone of Shandong Province[M].Beijing:China Science and Technology Press,1990:1-793.山东省科学技术委员会.山东省海岸带和海涂资源综合调查报告集(综合调查报告)[M].北京:中国科学技术出版社,1990:1-793.
[5]Compilation Committee of"Records of Bays in China"(CCRBC).Records of Bays in China(Volume 3:Bays in north and east of Shandong Peninsula)[M].Beijing:China Ocean Press,1991:179-469.中国海湾志编纂委员会.中国海湾志第3分册(山东半岛北部和东部海湾)[M].北京:海洋出版社,1991:179-469.
[6]CUI J R,XIA D X.The relationship between coastal morphology and the characteristics of waves and tides of Shandong peninsula[J].Journal of Oceanography of Huanghai &Bohai Sea,1992,10(3):20-25.崔金瑞,夏东兴.山东半岛海岸地貌与波浪、潮汐特征的关系[J].黄渤海海洋,1992,10(3):20-25.
[7]MASSELINK G,SHORT A D.The effect of tide range on beach morphodynamics and morphology:A conceptual beach model[J].Journal of Coastal Research,1993,9(3):785-800.
[8]ZHOU S L,YAN L F.Analysis and forecast of weather of storm tide in east Shandong[J].Meteorological science and technology,2002,30(1):37-42.周淑玲,闫丽凤.影响山东半岛风暴潮的分析和预报[J].气象科技,2002,30(1):37-42.
[9]SHORT A D,MASSELINK G.Embayed and structurally controlled beaches[C]∥SHORT A D.Handbook of beach.Chichester:Wiley,1999:142-161.
[10]KLEIN A H F,BENEDET F L,SCHUMACHER D H.Short-term beach rotation processes in distinct headland bay beach system[J].Journal of Coastal Research,2002,18(3):442-458.
[11]XIA D X,WANG W H.Formation and evolution of bays[C]∥CHEN Z S,WANG W H,WU S Y,et al.Introduction of Bays in China.Beijing:China Ocean Press,2007:179-469.夏东兴,王文海.海湾的形成与发育[C]∥陈则实,王文海,吴桑云,等编.中国海湾引论.北京:海洋出版社,2007:179-469.
[12]SILVESTER R,HSU J R C.Coastal Stabilization[M].Singapore:World Scientific,1997:578.
[13]WRIGHT L D,SHORT A D.Morphodynamic variability of surf zones and beaches:A synthesis[J].Marine Geology,1984,56:93-118.
[14]SHORT A D.Australian beach systems-nature and distribution[J].Journal of Coastal Research,2006,22(1):11-27.
[15]KORMAR P,GAUGHAN M K.Airy wave theory and breaker height prediction[C]∥Proceedings of the 13th conference on coastal engineering,ASCP,1972:405-418.
[16]FERGUSON R I,CHURCH M.A simple universal equation for grain settling velocity[J].Journal of Sedimentary Research,2004,74(6):933-937.
[17]ANTHONY E J.Sediment-wave parametric characterization of beaches[J].Journal of Coastal Research,1998,14(1):347-352.
[18]JACKSON D W T,COOPER J A G,DEL RIO L.Geological control of beach morphodynamic state[J].Marine Geology,2005,216:297-314.
[19]ZHANG X L,ZHU M Y,LI R X,et al.Simultaneous and consecutive multi-parameter monitoring of shellfish culture environment[J].Advances in Marine Science,2004,22(3),340-345.张学雷,朱明远,李瑞香,等.贝类养殖环境的多参数同步连续监测[J].海洋科学进展,2004,22(3):340-345.