渤海湾M2 分潮的季节变化：增强调和分析的应用

余鹰，王道胜,2*

( 1. 中国地质大学（武汉） 海洋学院 海洋地质资源湖北省重点实验室，湖北 武汉 430074；2. 中国地质大学 深圳研究院，广东 深圳 518057)

1 引言

由于海洋物理性质和海岸形态的变化，潮汐存在多个时间尺度的变化[1]，尤其是沿海和极地地区的M2分潮[2]。前人研究表明，英国海岸[3]、西北欧大陆架[4]、黄海和东海[5]、哈得孙河口[6]和东南亚海岸[1]等沿海地区，M2分潮都存在着季节变化。此外，相关研究表明极地地区M2分潮也存在季节变化[2,7-8]。

为了分析M2分潮的季节变化，传统上将观测的水位数据分段，并对每一段进行经典调和分析（Classical Harmonic Analysis, CHA），该方法可以通过MATLAB 中的T_TIDE 工具包来实现[9]。上述过程称为分段调和分析（Segmented Harmonic Analysis, SHA），该方法在前人研究中得到广泛使用[2,5,10-13]。当直接使用SHA 时，得到的振幅和迟角是离散的，需要对其进行插值以得到连续的结果。如果采用时间窗滑动的方式，一些峰值和快速变化的信号会被平滑，若某些分段中存在较多缺测，则可能导致计算结果不准确。除此之外，在使用SHA 进行月度分析时，一些无法分辨的潮汐成分也需谨慎处理。Jin 等[10]提出了增强调和分析（Enhanced Harmonic Analysis, EHA）方法，可以直接得到内潮不同潮汐成分时变的振幅和迟角。Pan等[14]基于EHA 方法开发了S_TIDE 工具包，并利用S_TIDE 研究了哥伦比亚河口的潮汐-河流相互作用。根据Pan 等[14]所述，在S_TIDE 中如何针对分析的潮汐成分选择合适的独立点数是一项相对复杂的工作，应通过一系列敏感性试验来测试，以避免产生虚假的结果。

在SHA 和EHA 中，所有潮汐成分的调和常数均被认为是随时间变化的。然而，潮汐季节变化的原因错综复杂[15]，在大多数情况下，只有主要分潮（如M2分潮）的时间变化是显著并需要关注的，而其他分潮的时间变化比较微弱且可以忽略不计。EHA 认为所有分潮的振幅和迟角都是随时间变化的，一方面增加了计算的复杂度和需要的计算机内存，另一方面显著减少了能够分辨的分潮数量。Wang 等[16]对EHA 进行了改进，使其不仅能得到主要分潮时变的振幅和迟角，同时也能得到其他分潮不变的振幅和迟角。为便于对EHA 改进前后的效果进行对比分析，改进前的EHA，即：假定所有潮汐成分的调和常数均随时间变化的EHA，被记为OEHA（Original EHA），而改进后的EHA 则简记为EHA。本文使用渤海湾的水位观测数据，通过设计理想实验，构建人造“水位数据”验证EHA 的可行性和有效性，并利用EHA 得到了渤海湾M2分潮振幅的季节变化，同时讨论了渤海湾M2分潮振幅季节变化的原因。

2 方法与数据

2.1 EHA

海平面高度由多个因素决定，可表达为[16-18]

式中， ζ(t)代表海平面高度(即：水位)，右边3 项分别是平均海平面高度(Mean Sea Level, MSL)、潮汐引起的海表面高度、气象因素导致的海表面高度。在不考虑节点调制的情况下，潮汐的计算通常是k个分潮之和，即

式中， ω，a和φ分别代表每个分潮的角速度（单位：rad/s），振幅（单位：m）和迟角（单位：rad）。在OEHA 中， ζ0，a和φ均被认为是随时间变化的[10,14]，即

在OEHA 的基础上，Wang 等[16]进行了改进。在改进后的EHA 中，只有几个主要分潮的振幅和迟角是时变的，而其他分潮的振幅和迟角仍然是常数[16]。因此，海平面高度的表达式为

式中，具有时变振幅和迟角的分潮数量为i, 而振幅和迟角是常数的分潮数量为j，ζ0(t)包含气象因素导致的海表面高度[16]。类似于Jin 等[10]和Pan 等[14]，改进后的EHA 在求解时变的平均海平面高度、振幅和迟角时，引入了独立点方案和三次样条插值法，而上述的时变量和等号右边第三项中的常数振幅和迟角可以用最小二乘法同时求解。此外，误差评估算法与Jin 等[10]和Pan 等[14]相同。改进后的EHA 已更新到S_TIDE 中， 可以从以下网站获取： https://www.researchgate.net/project/A-non-stationary-tidal-analysistoolbox-S-TIDE。

2.2 数据来源

本文采用的数据是中国渤海湾的两个站点E1 和E2 的逐时水位数据(已减去平均水位)，观测数据的时间跨度为2013 年11 月1 日(世界时，下同)到2014 年11 月1 日，图1 标注了E1 和E2 站点的位置。如图2所示，E1 站点有大约1/4 的数据缺测，而E2 站点的数据是连续的，E1 和E2 的海平面振幅大多小于2 m。

图1 渤海湾观测站E1 和E2 的位置（红星），以及10 m 和40 m 水深线Fig. 1 The locations of the observation stations (E1 and E2,red stars) in the Bohai Bay, and the isobathymetric lines at 10 m and 40 m

图2 E1（a）和E2（b）站的水位观测数据Fig. 2 The observed sea level at E1(a) and E2(b) stations

2.3 EHA 的验证

在分析实际水位数据时，分潮振幅和迟角的真实值并不清楚，无法验证方法的准确性。因此，本文首先通过理想实验来验证EHA 方法的可行性和有效性。

2.3.1 人造“水位数据”

根据E2 站的水位数据来构造理想实验中的人造“水位数据”。图3 呈现了E2 站点水位数据的功率谱密度，很明显，水位的波动主要受半日潮(如M2和S2)影响；同时，M2分潮的能量最强，全日潮(如K1和O1)也很显著。由于半日潮和全日潮的非线性相互作用和平流作用[19]，MK3和M4分潮的频率同样出现谱峰。低频部分也出现了一些谱峰，不失一般性，选择0.235 cpd和0.392 cpd 两个频段的水位信号作为低频代表，分别标记为LF1 和LF2。

图3 E2 站水位数据和人造“水位数据”的功率谱密度Fig. 3 The power spectral densities of the observed sea level at E2 Station and the artificial sea level

假定MSL 为0，故人造“水位数据”将由潮频信号(包括M2、S2、K1、O1、MK3和M4)和低频信号(包括LF1 和LF2)组成，如下式所示：

在人造“水位数据”中，假设M2分潮的振幅和迟角是随时间变化的，而其他信号的振幅和迟角均是常数。另外，MK3、M4、LF1 和LF2 被视为噪声，不估计其振幅和迟角。为了获得M2分潮振幅和迟角较为真实的时间变化，在不考虑节点调制的情况下，通过SHA分析了E2 的真实水位数据，如图4 所示。结果表明，M2分潮的振幅和迟角存在明显的时间变化。为了获得更加平滑的时间变化特征，利用正弦函数对SHA方法所得M2分潮振幅和迟角进行了拟合。图4 中的红线展示了拟合结果，拟合公式如下：

图4 使用SHA 分析E2 站水位观测数据所获得M2 分潮时变的振幅(a)和迟角(b)及其拟合值Fig. 4 The temporally varying amplitude (a) and phase lag (b) of the M2 constituent obtained by analyzing the observed sea level at E2 Station using SHA and their fitting values

式中，t表示从世界时0 年1 月0 日起经过的时间（即MATLAB 中的“datenum”值），单位：d；此外，t的范围为2013 年11 月1 日至2014 年11 月1 日，与E2 站实际水位观测数据的时间范围相同；a′M2是M2分潮的时变振幅，单位：m； φ′M2为M2分潮的时变迟角，单位：（°）。参考CHA[9]，由M2分潮引起的水位变化的表达式为下

式中，tmid是t的中值。

为了获得其他潮汐信号较为真实的振幅和迟角，在不考虑节点调制的情况下，使用CHA 分析了E2 处的海平面数据，得到的常数形式的振幅和迟角如表1所示。LF1 和LF2 的迟角设置为4 个主要分潮的平均迟角，并根据功率谱密度的结果，通过反复试验得到了振幅值，见表1。由上述信号引起的水位变化与M2分潮引起的水位变化具有类似的形式(即式(8))，此处不再赘述。

表1 理想实验中的潮汐成分及其常数形式的振幅和迟角Table 1 The prescribed constant amplitudes and phase lags of the constituents in ideal twin experiments

从人造“水位数据”的功率谱密度（图3）可以看出，人造“水位数据”成功地再现了E2 站M2、S2、K1、O1、MK3、M4、LF1、LF2 等信号，且不包含其他频段的信号，表明由已知的振幅和迟角构造的“水位数据”可以作为E2 站真实水位观测数据的替代品。

2.3.2 实验和结果

在理想实验中，不考虑节点调制，采用CHA、SHA、OEHA 和EHA 分析人造“水位数据”，以获得M2、S2、K1和O14 个主要分潮的振幅和迟角，除此之外，人造“水位数据”中的其他信号，被视作噪声[20]。在使用SHA 时，每个月的数据作为一段进行分析，根据瑞利准则，区分上述4 个主要分潮的最小周期为14.8 d[21]，因此，只有当该段数据长度大于15 d 时，SHA 才对该段数据进行分析。此外，采用3 次样条插值法对每段所得离散的振幅和迟角进行插值以获得连续的结果。Pan 等[14]指出，在不同独立点数下，EHA 得到的时变振幅和迟角代表了不同时间尺度上的变化。根据Jin 等[10]的研究结果，在下面的理想实验中，OEHA和EHA 的独立点数被设为5，比季节数（春、夏、秋、冬）多1。

如图5 所示，SHA、OEHA 和EHA 得到的M2分潮时变的振幅和迟角几乎等于给定值；CHA 得到的结果是恒定的，不能捕捉到M2分潮的时间变化，而SHA、OEHA 和EHA 可以得到M2分潮的时间变化特征。

图5 理想实验中给定的M2 分潮振幅（a）和迟角（b），以及使用SHA、采用5 个独立点的OEHA、采用5 个独立点的EHA 和CHA 得到的结果Fig. 5 The amplitude (a) and phase lag (b) of the M2 constituent prescribed in ideal experiments and the estimated values using SHA,OEHA with 5 independent points, EHA with 5 independent points, and CHA

对于S2分潮而言，EHA 分析所得振幅和迟角比CHA 所得结果更接近给定值（图6a 和图6b），而SHA的结果与给定值相差甚远，尤其是迟角（图6b），这主要是因为M2分潮和S2分潮的频率很接近，当同时被假定为随时间变化时，分析结果可能会互相影响。使用OEHA 时，除2014 年10 月和11 月外，S2分潮的迟角几乎都等于给定值；此外，得到的S2分潮的振幅是时变的，与给定值有一定的偏差。对于K1分潮（图6c和图6d）和O1分潮（图6e 和图6f），CHA 和EHA 所得结果几乎等于给定值，而OEHA 得到的振幅和迟角整体上接近给定值，但在端点处差距较大，SHA 分析所得振幅和迟角具有明显的时间变化特征，与给定值相差较大。总的来说，虽然SHA、OEHA 和EHA 均能成功估计M2分潮的时间变化特征，但EHA 所得其他3 个主要分潮的振幅和迟角与给定值吻合最好。

图6 理想实验中给定的S2 分潮（a，b）、K1 分潮（c，d）和O1 分潮（e，f）的振幅和迟角，以及使用SHA、采用5 个独立点的OEHA、采用5 个独立点的EHA 和CHA 得到的结果Fig. 6 The amplitude and phase lag of the S2 (a, b), K1 (c, d) and O1 (e, f) constituent prescribed in ideal experiments and the estimated values using SHA, OEHA with 5 independent points, EHA with 5 independent points, and CHA

需要注意的是，E1 站的实际水位观测数据存在部分缺测，这可能会对分析结果造成影响。在理想敏感实验中，根据E1 站的缺测时段，剔除了人造“水位数据”中相应的数据，进而利用CHA、SHA、OEHA、EHA 对其进行了分析。虽然大约1/4 的时段没有“水位数据”，但是OEHA 和EHA 所得结果很好地捕获了M2分潮给定振幅和迟角的时间变化特征（图7a和图7b）。相反，SHA 所得M2分潮的振幅和迟角与给定值略有偏差，特别是在2013 年11 月至2014 年2 月期间，这主要是因为2013 年12 月和2014 年1 月的缺测较多，导致这两个月的“水位数据”没有被分析而无法获得M2分潮振幅和迟角值。EHA 分析得到的M2分潮所引起的水位变化与给定的M2分潮所导致的水位变化之间的绝均差为0.19 cm；而SHA 所得绝均差为0.28 cm，约为没有缺测时的4.5 倍。对于S2、K1和O1分潮，CHA 分析所得振幅和迟角值与给定值之间的差距比理想实验中的差距更大，如图8所示；此外，“水位数据”的缺测对SHA 和OEHA 的结果也有很大的影响，OEHA 分析所得3 个分潮的振幅和迟角与给定值之间的误差变得更大。相反，EHA得到的S2和K1分潮的振幅和迟角仍然接近给定值。虽然EHA 对O1分潮的分析结果与给定值略有偏差，但差异明显小于CHA、SHA 和OEHA 的结果。由此可见，EHA 总体上优于CHA、SHA 和OEHA。

图7 理想敏感实验中给定的M2 分潮振幅（a）和迟角（b），以及使用SHA、采用5 个独立点的OEHA、采用5 个独立点的EHA 和CHA 得到的结果Fig. 7 The amplitude (a) and phase lag (b) of the M2 constituent prescribed in ideal sensitivity experiments and the estimated values using SHA, OEHA with 5 independent points, EHA with 5 independent points, and CHA

图8 理想敏感实验中给定的S2 分潮（a, b）、K1 分潮（c, d）和O1 分潮（e, f）的振幅和迟角，以及使用SHA、采用5 个独立点的OEHA、采用5 个独立点的EHA 和CHA 得到的结果Fig. 8 The amplitude and phase lag of the S2 (a, b), K1 (c, d) and O1 (e, f) constituent prescribed in ideal sensitivity experiments and the estimated values using SHA, OEHA with 5 independent points, EHA with 5 independent points, and CHA

总体而言，无论“水位数据”中是否存在缺测，在利用EHA 进行分析时，M2分潮时变振幅和迟角以及其他分潮调和常数的给定值都能被准确反演，说明EHA 能有效分离具有不同变化特征的潮汐信号，可用于实际实验中M2分潮的季节变化研究。

3 EHA 在渤海湾M2 分潮季节变化研究中的应用

3.1 基准实验

为了确定基准实验中所需分析的潮汐成分，利用CHA 对E1 和E2 站的水位观测数据进行了调和分析，根据信噪比的大小[9,22]选择了M2、K1、S2和O14 个最重要的成分作为进行分析的分潮。利用EHA进行数据分析时，独立点数设置为5，仅将M2分潮的振幅和迟角假定为随时间变化。

CHA 和EHA 所得E2 站M2分潮振幅的分析结果如图9b 所示。从图中可以看出，EHA 分析所得M2分潮振幅具有明显的季节变化特征，其平均值与CHA分析所得常数形式的M2分潮振幅几乎相等，说明了EHA 所得分析结果的合理性。

图9 E2 站基准实验的实验结果Fig. 9 The results of benchmark experiment at E2 Station

EHA 所得随时间变化的M2振幅在夏季（大约7-8 月）达到最大值，在冬季（大约1-2 月）达到最小值，夏季M2分潮振幅比年平均值大8.42 cm（约11.01%），而冬季值比年平均值小8.76 cm（约11.46%），其季节变化幅度远大于在观测期内由18.6 a 的节点调制引起的振幅变化（不超过1%）。表2 列出了CHA 和EHA在E2 站分析所得其他分潮的振幅值，可以看出，EHA和CHA 分析所得其他3 个主要分潮的常数形式的振幅值之间的绝均差最大仅为0.07 cm，进一步表明EHA 不仅可以准确提取M2分潮振幅的季节变化特征，也能准确得到其他分潮的常数形式的振幅。

CHA 和EHA 分析所得E2 站MSL 如图9a 所示。可以看出，EHA 分析所得随时间变化的MSL 的平均值几乎等于CHA 分析所得常数形式的MSL，进一步表明EHA 的结果是合理的。但是需要注意的是，EHA分析所得E2 站处MSL 在1 年内的变化约为0.4 m，远大于传统上MSL 的全球平均变化（约为0.001 m/a）[23-24]，这主要是因为EHA 中随时间变化的MSL 包含气象因素导致的水位变化（式（3））。如Lü等[18]所述，风速和气压等气象因素均会影响到低频水位的变化，因此包含了气象因素影响的MSL 在2013 年11 月和2014 年11 月存在一定的差异，该现象在E1 站同样存在。EHA分析所得时变MSL 的最大值出现在夏季，而最小值出现在冬季。与M2分潮振幅的时间变化特征相似，EHA 分析所得MSL 的时间变化比较平稳，且不存在季节变化之外的信号。

CHA 和EHA 分析所得E1 站MSL 和M2分潮振幅的结果如图10 所示。可以看出，EHA 分析所得随时间变化的M2分潮振幅比较平滑，且与E2 站具有相同的变化趋势。在夏季，M2分潮振幅增加了约6.05%，而在冬季则减少了3.85%，这比在观测期内由18.6 a节点调制引起的不超过1%的振幅变化大得多。从表2 中可以看出，EHA 分析所得其他3 个分潮的常数形式的振幅与CHA 分析所得结果非常接近。此外，EHA 分析所得随时间变化的MSL 在夏季达到最大值，在冬季达到最小值，与E1 站M2分潮的季节变化特征相同。需要注意的是，E1 站和E2 站相距较近，但对比图9a 和图10a 发现，E1 站和E2 站的平均海平面的年变化并不完全一致，E2 站MSL 的最小值出现在11 月，而E1 站MSL 的最小值出现在1 月，这是由于E1 站2014年1 月和2 月的数据存在缺测，会给该时间段的平均海平面的计算带来偏差，从而影响最小MSL 出现的月份，但并没有显著影响两站的M2分潮振幅和MSL 的变化趋势，表明EHA 在即使存在缺测的E1 站，也可以得到相对合理的结果。

图10 E1 站基准实验的实验结果Fig. 10 The results of benchmark experiment at E1 Station

表2 基准实验中CHA 和EHA 分析所得4 个主要分潮的振幅（单位：cm）Table 2 The amplitudes of four main tidal constituents obtained by CHA and EHA in the benchmark experiment (unit:cm)

综上所述，EHA 在E1 站和E2 站分析所得M2分潮振幅具有相同的季节变化特征：夏季大，冬季小。E2站M2分潮振幅的相对变化约为11%，在E1 站约为5%，这与Foreman 等[11]、Huess 和Andersen[4]的结果（6%），以及Müller 等[2]根据模型和观测数据所分析的结果（5%～10%）相近。此外，EHA 在E1 站和E2 站分析所得时变MSL 具有与M2分潮振幅相同的季节性变化，这表明MSL 可能是导致M2分潮振幅季节变化的因素之一，因为潮波的传播受到水深的影响[13]。

3.2 对分潮选择的敏感性实验

在基准实验中，选择了4 个主要分潮进行评估，并且假定只有M2分潮的振幅和迟角随时间变化。当其他分潮也被分析或其他分潮振幅和迟角也被假定为随时间变化时，M2分潮的分析结果可能会受到影响。因此，对E2 站的数据进行了一些新的实验来测试EHA 分析所得的M2分潮振幅季节变化特征对分潮选择的敏感性。如表3 所示，在敏感性实验SE11-SE13中，对M2，K1，S2和O14 个分潮进行分析，这与基准实验相同，但具有时变振幅和迟角的分潮不同；在敏感性实验SE21-SE24 中，与基准实验不同的是，将对M2、K1、S2、O1、N2、K2、P1和Q18 个分潮进行分析，且振幅和迟角随时间变化的分潮也发生了变化。

表3 实际实验中的详细实验设置Table 3 The detailed experimental settings of the practical experiments

上述敏感性实验中，EHA 分析所得随时间变化的MSL 和M2分潮振幅如图11 所示。可以看出，当具有时变特征的分潮发生改变时（SE11-SE13），EHA分析所得随时间变化的M2分潮振幅和MSL 几乎与基准实验的结果完全一致。当分析的分潮数目发生变化时（SE21-SE24），EHA 分析所得随时间变化的MSL与基准实验的结果几乎完全一致（图11c）；此外，虽然随时间变化的M2分潮振幅与基准实验结果不完全相同，但整体变化趋势是一致的（图11d）。值得注意的是，SE21-SE24 中随时间变化的M2分潮振幅几乎彼此相等，这与SE11-SE13 所得结论一致。上述结果表明，EHA 中被分析的分潮和具有时变特征的分潮的选取并不会影响分析所得M2分潮振幅和MSL 的时间变化趋势和季节变化特征。

图11 分潮选择的敏感性实验结果Fig. 11 The results on the sensitivity of tidal constituent selection

3.3 对数据缺失的敏感性实验

时变M2分潮振幅和MSL 的分析结果也可能会受到数据缺测的影响。因此，利用E2 站的数据进行了一系列实验以测试EHA 分析结果对数据缺测的敏感性。如表3 所示，在敏感性实验SE31 中，根据E1 点缺测情况，相应地剔除了E2 站水位观测数据。E1 站水位数据中存在大约1/4 的缺测，因此，为了进一步检验数据缺失的影响，将E2 站水位观测数据平均分为4 份，并将这4 份水位观测数据分别作为SE32-SE35中的缺测数据。其他实验设置与基准实验完全相同。

敏感性实验SE31-SE35 的结果如图12 所示。可以看出，由于SE31-SE35 的水位观测数据存在约1/4 的数据缺失，因此EHA 分析所得M2分潮振幅和MSL 均与基准实验中的值不完全相同。对比实验SE31 和基准实验中分析所得E2 站MSL 的结果，可以看到当E2 站数据在2014 年1 月和2014 年2 月出现缺测时，虽然MSL 的最小值出现的月份在2013 年11 月，但是这种现象并没有数据无缺测时那么显著，并且2013 年11 月至2014 年2 月的分析结果与无缺测时的分析结果存在明显的差异，说明了数据缺失是导致E1 站和E2 站MSL 最小值出现月份不同的原因。但是，SE31-SE35 中所得M2分潮振幅和MSL 的变化趋势与基准实验中的相同：夏季大，冬季小。实验结果表明，数据在某一月份的集中缺测会给该月的计算结果带来影响，但不会显著影响EHA 分析所得M2分潮振幅和MSL 的季节变化特征，并进一步论证了EHA所得E1 站的M2分潮振幅季节变化是合理、可信的。

图12 数据缺失和独立点数的敏感性实验结果Fig. 12 The results on the sensitivity of missing data and independent points

3.4 对独立点数的敏感性实验

Jin 等[10]指出，如何选择EHA 中的独立点是一个需要进一步研究的问题。而Pan 等[14]则指出，EHA在不同独立点数时所得到的分析结果代表了不同时间尺度上的振荡。同样地，选择合适的独立点数对EHA 的表现也非常重要，尤其是对于E1 站非连续的水位数据。因此在新的敏感性实验中尝试了4 个、6 个、7 个和8 个等不同的独立点数。

如图12d 所示，当独立点数取为季节数的两倍（即：8）时，EHA 分析所得E2 站的M2 分潮振幅的时间变化趋势与基准实验结果一致，但存在除季节变化之外的震荡，结果表明：EHA 中，独立点越多并不一定会得到更真实的结果。

当独立点数大于3 且小于8 时，EHA 分析所得M2振幅和MSL 的结果并不完全与5 个独立点时的结果等同，但其季节性变化趋势相同。如Pan 等[14]所述，不同独立点数得到的结果不同，可能是因为采用不同独立点数的EHA 得到的结果包含了不同的信号。以上所有敏感性实验的结果表明：EHA 分析所得E2 站M2分潮振幅的季节变化特征是鲁棒的，不受实验设置的影响。

4 讨论

平均海平面、层化和涡动黏性系数的季节变化是造成沿海地区M2分潮振幅呈现季节性变化的主要原因[4,13,15]。渤海湾的风场不仅对平均海平面有重大影响[18]，因受东亚季风影响而具有明显的季节变化特征[25]。此外，季风还可以影响层化涡动黏性系数[26-29]。参照Devlin等[1]，采用北太平洋季风指数(WNPMI)来表征东亚季风[30]。WNPMI 下载网址为：http://apdrc.soest.hawaii.edu/projects/monsoon/daily-data.html#mon。如图13 所示，WNPMI 与E1 站和E2 站M2分潮振幅具有显著的正相关，相关系数分别为0.86 和0.87，表明：东亚季风可能是造成渤海湾M2分潮振幅季节性变化的内在原因。

图13 北太平洋季风指数（WNPMI）与E1 站和E2 站M2 分潮振幅的相关性Fig. 13 The correlation between WNPMI and M2 constituent amplitude of E1 and E2 stations

根据前人研究结果，我们大胆猜测，东亚季风可能通过以下3 个途径来影响渤海湾M2分潮振幅的季节变化：

（1）影响平均海平面的季节变化。亚洲季风导致了径流的季节性变化[31]，使得超过60%的黄河水在夏季入海[32]，也使得海河夏季的入海量是冬季的两倍左右[33]；与此同时，夏季东亚季风导致的渤海上空的东南风使一部分黄河水向北输送到渤海湾[28]，加上直接汇入渤海湾的海河入海量也增大，使得夏季渤海湾的淡水汇入增多[34]。另一方面，受东亚季风影响，夏季通过渤海海峡从渤海进入黄海的通量比冬季更小[27]。综上，受东亚季风影响，夏季从黄河和海河汇入渤海湾的淡水增加，与此同时，从渤海流出到黄海的海水更少，共同导致夏季渤海湾平均海平面更高。而潮波的传播受到水位的影响[13]，因此，受东亚季风季节变化影响的渤海湾平均海平面的季节变化导致了M2分潮振幅的季节变化。

（2）影响层化的季节变化。根据Yanagi 等[35]和Jeon 等[26]的研究，对于渤海湾而言，冬季从陆地到海洋的强西北风会导致海水降温，破坏之前形成的层化并导致垂向均一的层化状态。春末夏初，随着海水表层升温速率的增加以及黄河和海河入海量的增加，海面风应力和底部摩擦力引起的混合不足以消耗输入的势能，使得渤海湾层化逐渐加强[36-37]，并在夏季达到最大值[37]。因此，东亚季风能够导致层化的季节变化，进而引起M2分潮的季节变化[2]。

（3）影响涡动黏性系数。Pohlmann[38]发现北海的垂向涡黏性具有显著的年周期并与海面风应力的年变化有关；此外，垂直涡动黏性系数依赖于层化。Maas 和Van Haren[39]证实了层化情况下涡动黏性系数降低的可能性。Howarth[40]认为温跃层的存在可以

在理想实验中，根据E2 站的水位观测数据构建了人造“水位数据”，其中，M2分潮具有时变的振幅和迟角，而S2、K1、O1、MK3、M4分潮及LF1 和LF2 的振幅或迟角是常数。实验结果表明：无论人造“水位数据”是否存在缺测，EHA 分析所得随时间变化的M2振幅和迟角以及不随时间变化的S2、K1和O1分潮的振幅和迟角整体比CHA、SHA、OEHA 等方法的分析结果更接近给定值，表明EHA 是一种行之有效的方法。

利用EHA 分析实际水位观测数据时，分析结果发现，E1 站和E2 站时变M2分潮振幅具有相同的变化趋势：夏季大、冬季小。此外，M2分潮振幅的相对变化在E2 站约为11%，在E1 站约为5%。E1 站和E2 站随时间变化的平均海平面与M2分潮振幅随时间变化的趋势相同。EHA 分析所得其他潮汐成分（包括S2，K1和O1）的常数形式的振幅值与CHA 所得结果吻合较好。更重要的是，敏感实验的结果表明，M2分潮振幅和平均海平面的变化趋势是鲁棒的，不受实验设置的影响。根据前人研究[2-3,11,13,41]猜测，东亚季风的季节变化引起了平均海平面、层化和涡动黏性系数的季节变化，进而诱导了渤海湾M2分潮振幅的季节变化。降低涡动黏性系数，除此之外，Müller 等[2]指出，层化可以使得水体更加稳定，进而导致夏季的涡动黏性系数比混合均匀的夏季小几个数量级。因此，受东亚季风影响的层结的季节变化会导致涡动黏性系数的季节变化，进而诱发M2分潮的季节变化。

除了东亚季风可能影响渤海湾M2振幅的季节变化之外，水深、近岸与否等原因也有可能会对M2分潮振幅的季节变化造成影响。对比图9 和图10 可知，E1 站和E2 站的MSL 和M2分潮振幅的年变化并不完全一致，这种现象一方面是由数据缺测造成的，另一方面也可能与E1 和E2 两站的水深、近岸与否的差异有关。Müller 等[2]认为，相比远离海岸的地区，近岸地区的M2分潮的季节变化更显著。而在渤海地区，不同的水深和底摩擦系数也会对M2分潮的振幅带来影响[42]。

5 结论

本文利用EHA 分析了渤海湾E1 站和E2 站一年时间段的水位观测数据，同时估计了时变的和不随时间变化的分潮振幅和迟角，来探究渤海湾M2分潮振幅的季节变化。