海平面上升、强台风和风暴潮对厦门海域极值水位的影响及危险性预估

许炜宏，蔡榕硕*

( 1. 自然资源部第三海洋研究所，福建 厦门 361005)

1 引言

在气候变暖背景下，全球海平面加速上升，强台（飓）风趋于频繁，强台（飓）风和风暴潮导致低海拔（小于10 m）沿海地区的洪涝灾害日趋严重[1]。这是由于海平面的持续上升，极端海面高度（极值水位）有升高的现象，从而引起低海拔沿海地区极易被海水淹没或发生洪涝，而未来极值水位事件的重现期还有明显缩短的趋势[2]。极值水位主要由海平面高度、风暴增水和潮汐等构成，而海平面高度的小幅上升也会显著增加沿海地区洪水灾害发生的频率和强度，这是因为海平面的上升抬高了风暴潮、潮汐和波浪的基础水位，并且这种变化在陆架海区特别明显，尤其是当风暴潮、波浪和天文大潮叠加到上升的海平面后产生的极值水位对沿海地区低洼地造成的洪涝灾害尤其严重[3]。

近几十年来，台（飓）风和风暴潮造成全球沿海地区的洪水灾害及社会经济的损失与日俱增。1975-2016 年，全球80.2%因洪水死亡的人口位于距海岸线100 km 的地区内[4]。其中，2005 年的“卡特尼娜”飓风袭击美国新奥尔良，造成新奥尔良80%地区被洪水淹没，部分地段水深高达6 m，并导致约1 500 人死亡，至少1 000 亿美元的经济损失[5]。1949-2009 年，我国沿海地区共发生了220 余次较为严重的台风风暴潮灾害[6]。《中国海洋灾害公报》（1989-2018 年）显示，台风风暴潮是造成我国直接经济损失最严重的海洋灾害[7]；并且，自2000 年以来，台风风暴潮的出现频率较之前半个世纪有显著的增加[8]。

基于《中国海洋灾害公报》（1989-2018 年）公布的数据，本文统计表明，我国东南沿海地区特别是福建省和广东省是受到台风风暴潮影响最大的地区（表1）。其中，厦门市位于福建省南部沿海，台湾海峡西岸，由厦门岛、鼓浪屿和位于大陆沿海的同安区、翔安区、集美区和海沧区等组成，面积约为1 699 km2，有234 km的海岸线，常住人口411 万，是我国东南沿海重要的中心城市。而厦门岛南部和东南部地势低平，呈现朝向东南开口的汤匙形态。历史上5903 号、9914 号、1614号等台风都曾引发了厦门市尤其是厦门岛低洼地较为严重的海岸洪水[9]。研究表明，气候变化背景下未来东亚地区热带气旋的路径将逐渐偏北，并且达到强度峰值时的纬度越来越高[10]；到2050 年，厦门市将可能是全球沿海地区受到洪水影响最大的城市之一，届时，海岸洪水每年给厦门市造成的经济损失将达到全市GDP 的0.22%[11]。因此，气候变化背景下海平面上升、台风和风暴潮对我国东南沿海尤其是厦门地区的影响与灾害风险，以及应对策略是值得高度重视的科学与社会问题。

表1 1989-2018 年我国沿海各省因风暴潮造成损失情况Table 1 The losses caused by storm surges in each coastal provinces of China from 1989 to 2018

1980 年以来，我国沿海海平面呈现明显的上升趋势，上升速率达到3.3 mm/a，高于同期全球平均速率[12]，并且自2012 年以来沿海海平面不断达到新的高度。其中，1993-2017 年，沿海海平面上升速率和幅度分别为3.9 mm/a、100 mm[13]，也高于同期全球平均的3.0～3.1 mm/a、80 mm[14-15]。其中，1993-2012 年，福建和广东沿海海平面上升速率甚至分别达到3～4 mm/a、(3.6±0.7) mm/a[16-17]。厦门地区潮位站的观测结果显示，1993-2012 年厦门沿海海平面以3.1 mm/a 的速度上升[16]。2001 年之后影响我国东南沿海地区的台风、强台风频率增加[18]，风暴潮灾害的频次和强度也随之提高。由此可见，未来在海平面持续上升背景下，台风和风暴潮将引起极值水位的进一步上升，并导致其危险性（致灾强度和频率）的增加，从而对我国东南沿海地区特别是厦门市沿海地区带来更大影响和灾害风险。

然而，目前关于台风和风暴潮对沿海城市洪水灾害的影响研究主要侧重于单个致灾因子如风暴潮[19]、海平面上升[20]对海岸洪水产生的影响，而对于台风、风暴潮、天文潮、潮流和极端降水等多致灾因子对沿海地区极值水位的叠加影响研究还相对较少。其中有关多致灾因子对我国东南沿海的厦门等地区的综合影响研究尚未见专门报道，而这对于沿海地区低洼地的防灾减灾有重要的意义和迫切性。为此，本文选择了历史上对厦门沿海地区造成严重影响的9914 号“丹恩”和1614 号“莫兰蒂”台风事件，分析了这两次台风影响期间的风暴增水、天文潮位、潮流场和降水过程，以及大气环流形势的变化，评估了海平面上升、风暴潮增水、天文潮和降水等多致灾因子对厦门沿海地区极值水位的叠加影响，预估了不同气候情景下未来海平面的变化和台风等极端事件对厦门沿海地区极值水位危险性的影响，以期为我国东南沿海地区应对气候和减灾防灾提供科学依据。

2 数据与方法

2.1 资料

2.1.1 大气数据

本文选用了9914 号“丹恩”和1614 号“莫兰蒂”两个台风案例，其中，台风过程的最大风速、逐时降水站点数据源用中国气象局（CMA）厦门基本站（24.48°N，118.07°E）数据资料（http://data.cma.cn/data/detail/dataCode/A.0012.0001.html），时间覆盖范围为1951 年1 月至2016 年12 月。台风路径资料源用中国气象局的热带气旋最佳路径数据集（CMA-BST）（http://tcdata.typhoon.org.cn/zjljsjj_zlhq.html），时间覆盖范围为1949-2018 年。台风的大气环境背景场数据源用欧洲中期数值预报中心ERA-5 全球逐小时再分析格点资料，地面10 m 风场以及海平面气压场（https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/）的空间分辨率为0.25°×0.25°，时间覆盖范围为1979 年1 月至2016 年12 月。

2.1.2 海洋数据

本文以厦门海洋环境监测站(24.44°N，118.08°E，简称厦门验潮站)潮位数据代表厦门沿海的潮位特征。潮位资料来自厦门海洋环境监测站的逐时观测资料，以及夏威夷大学海平面中心（UHSLC）的厦门验潮站逐时潮位数据资料，时间覆盖范围为1954-1999 年（http://uhslc.soest.hawaii.edu/）。海洋环境背景场数据选用欧洲中期数值预报中心的ERA-5 全球逐小时再分析格点资料（分辨率为0.5°×0.5°）中的波浪方向及有效波高数据。未来海平面高度数据选取国际耦合模式比较计划第五阶段（CMIP5）的模式结果，主要为CMCC-CM、CNRM-CM5、MIROC-ESM 等29个模式中2030-2100 年的模拟试验数据结果，并采用多年加权平均获得未来海平面上升预估值（可能变化范围）[21]。

2.2 方法

本文应用了数值模拟方法、水文概率统计模型和Kendall 相关系数统计方法等。具体细节如下。

2.2.1 潮流场模拟

平面二维潮流基本方程组包括连续方程和动量方程如下：

式中，h表示水位；H表示水深；u、v分别表示x、y方向的流速分量；f为科氏力系数；C表示谢才系数；t表示时间；g为重力加速度。

初始条件为

具体的计算参数和条件如下所述：

（1）边界条件：岸边界法向量取0，水边界取自于全球近海潮波模型TPOX 7.2。

（2）计算区域：24.35°～24.68°N，117.96°～118.47°E范围的海域。

（3）岸线数据：采用美国国家海洋和大气管理局（NOAA）的全球自洽分层高分辨率地理数据库(GSHHG v2.3.7)。

（4）水深地形数据：取自NOAA 的ETOPO1（1-Minute Gridded Global Relief Data Collection）全球地形起伏模型，空间分辨率为1′，选取范围为24.3°～24.7°N，117.8°～118.5°E。

（5）风力强迫项数值：取自欧洲中期数值预报中心ERA-5 全球逐小时再分析格点资料（分辨率为0.25°×0.25°）中的地面10 m 风场。

2.2.2 水文概率统计模型

2.2.3 肯德尔（Kendall）相关系数

肯德尔相关系数是一个用来测量两个随机变量相关性的统计值。假设两个随机变量分别为X、Y（也可以看做两个集合），它们的元素个数均为N，X与Y中的对应元素组成一个元素对集合XY。肯德尔相关系数τ的取值范围在-1 到1 之间，当τ为1 时，表示两个随机变量拥有一致的等级相关性；当 τ为-1 时，表示两个随机变量拥有完全相反的等级相关性；当 τ为0 时，表示两个随机变量是相互独立的。肯德尔相关系数表达式为

本文中，X为极值水位高度，Y为对应时刻的降水量，C表示集合XY中拥有一致性的元素对数（两个元素为一对）；D表示集合XY中拥有不一致性的元素对数[24]。

2.2.4 未来极值水位计算方法

未来极值水位高度（ESL）主要通过海平面上升（SLR）的情况和未来风暴潮事件的水位（ ST）叠加估算而得[25]

3 结果与分析

3.1 9914 号“丹恩”和1614 号“莫兰蒂”台风概况

1999 年10 月2 日，9914 号台风“丹恩”在菲律宾以东洋面生成，经菲律宾进入南海后，突然北折，正面袭击我国闽南地区，并于北京时间1999 年10 月9 日上午10 时许，在漳州镇海角附近登陆，登陆时中心最大风力35 m/s，中心最低气压970 hPa（图1）。10 月9 日中午12 时，台风再次登陆厦门市海沧区。期间，恰逢农历九月初一的天文大潮，台风风暴增水叠加天文大潮，引起潮位异常增高。与此同时，台风带来的短时强降水，造成厦门沿海地区严重洪水和内涝。此次台风造成厦门地区多处堤防损毁、12 人死亡、4 人失踪、727 人受伤，直接经济损失约18.6 亿元。

2016 年9 月10 日14 时1614 号台风“莫兰蒂”在菲律宾以东洋面生成后，向西北方向行进，穿过巴士海峡后，直接扑向厦门，成为自“5903”号台风后，登陆闽南的最强台风（图1）。2016 年9 月15 日凌晨3 时许在厦门市翔安区登陆，登陆时中心最大风力52 m/s，中心最低气压935 hPa。此次台风造成厦门地区因灾死亡1 人、重伤2 人，直接经济损失约102亿元[26]。

3.2 9914 号台风期间海洋和大气过程对极值水位的影响

3.2.1 9914 号台风期间天文潮位、风暴增水和降水过程对极值水位的影响

图2 为1999 年10 月9 日01 时至20 时，9914 号台风过程的降水、风暴增水和厦门验潮站潮位的变化。由图2 可见，在9 日（农历九月初一）9914 号台风登陆前（03 时至12 时），厦门验潮站的潮位逐步升高，到12 时达到天文大潮最高潮位585 cm，并出现了实测潮位为732 cm 的极值水位（厦门验潮站警戒水位700 cm）；期间，风暴增水基本保持着1.2 m 以上的增水幅度，最大增水达1.71 m，尤其是从08 时至11 时，风暴增水基本处于最高值。台风过程累计降水量达208 mm，每小时降水量逐渐增大，至12 时达到过程降水量的峰值（29 mm/h），即每小时最大降水量。与此同时，厦门验潮站水位达到峰值，即出现732 cm 的极值水位。即天文大潮位、风暴增水和降水量最大同时出现（所谓“三碰头”），三者的共同作用造成了异常偏高的超警戒极值水位的出现，从而导致厦门岛沿岸低洼地发生严重的洪涝灾害。有研究表明，当强风暴潮与强降水相遇时，出现洪水的概率比二者单独出现要高得多[27-28]。

图2 9914 号台风“丹恩”期间厦门海域的天文潮位、风暴增水和降水过程Fig. 2 The astronomical tide, storm surge and precipitation in coastal waters of Xiamen during the period of No.9914 Typhoon Dan

3.2.2 9914 号台风期间大气环流对极值水位的影响

图3 为9914 号台风“丹恩”登陆前后的大气背景场变化。该图是采用经实测数据修正的ERA-5 逐时海平面风场和气压场再分析资料，以及CMA 台风最佳路径集资料绘制的。研究表明，风力和强降水均是风暴潮增水的重要因素。随着台风向岸抵近，大量海水涌入厦门海域的同时，由于风场和低气压的作用，并且到台风登陆时（12 时），降水量达到最大，为29 mm/h，厦门海域内出现持续且较强的风暴增水。由图3 可见，本次台风在登陆前后的移动速度偏慢。观测资料显示，厦门气象基本站在9 日07 时至当天15 时出现持续的12 级（32.7 m/s）大风。台风登陆前的持续大风时间基本与厦门验潮站的大幅度增水时间重叠。例如，在台风登陆前的一段时间（9 日06 时前）（图3a，图3b），福建沿海以偏北大风为主。在台风即将登陆前的时段内（9 日10-12 时），厦门海域的风向逐步转为东南风、南风和西南风向，大气的低层风场有利于驱动外海的海水向厦门海域方向输运堆积；同时,该海域处于台风低压区，这又有利于台风对海水的抽吸作用。值得注意的是，台风强度在临近登陆前还出现了一次增强，最大风速达到40 m/s，中心最低气压为965 hPa。简言之，台风登陆前（12 时之前），在持续向岸强风的驱动和台风中心低压的抽吸作用下，外海海水不断向厦门海域内堆积，造成厦门验潮站从9 日凌晨至中午观测到较长时间的显著增水现象。而台风登陆后（12 时之后），随着台风移动方向和位置的变化，台风不断北移、风力减弱以及风向转变，并且低压中心的气压逐步升高，因此，有利于风暴增水的条件随之消失。

图3 1999 年10 月9 日9914 号台风“丹恩”登陆前后厦门海域海平面气压和地面10 m 风场的变化（资料源自ERA-5）Fig. 3 The change of sea level pressure and the 10 m wind field in coastal waters of Xiamen during the landing time of No.9914 Typhoon Dan on October 9, 1999 (data from ERA-5)

3.2.3 9914 号台风期间潮流和波浪对极值水位的影响

图4 为模拟的9914 号台风过程中厦门岛周边海域（简称厦门海域）4 个典型时刻的表层潮流场变化。本文关于厦门海域潮流场的潮位和潮流模拟验证与文献[29]验证材料与方法一致，故省略。由图可见，在涨潮阶段，潮波由开边界传入，被金门岛阻断分为南北两支。北支经厦门岛北部进入同安湾海域。南支达到厦门岛南部海域时，又被分为两支，一支沿岸线北上，而另一支则进入九龙江河口区。

图4 1999 年10 月9 日9914 号台风“丹恩”期间厦门岛周边海域的表层潮流场Fig. 4 Surface flow field in coastal waters of Xiamen during the period of No.9914 Typhoon Dan on October 9, 1999

9914 号台风登陆前（06-10 时），厦门海域正处于从低平潮向高平潮过渡、落潮流向涨潮流转变过程中。至10 时，由外海向厦门海域的涨潮流达到涨急阶段，外海的海水大量涌入厦门海域内；至12 时，涨潮流流速趋于最小，潮位转为高平潮，来自外海的海水在厦门海域内堆积，使得潮位达到最高位。台风期间，向岸的持续风场有利于外海表层水向厦门海域输送，增强了涨潮流对厦门岛东南部沿岸形成超警戒极值水位的作用。在12 时的高平潮之后，潮流场逐步转为落潮流，至15 时，潮流达到落急阶段；随后，台风北移、风向改变，气压场削弱，潮位迅速回落。

图5 为9914 号台风期间厦门附近海域的平均波浪方向和有效波高的分布。厦门海域的波浪方向与地面风向出现了同步的变化。台风登陆前（06-11 时），在向岸大风的驱动下，厦门海域掀起了巨浪，06 时波浪由外海向厦门海域传播（图5b），至12 时台风中心附近海域的有效波高达到了4 m 以上（图5c）。

图5 1999 年10 月9 日9914 号台风“丹恩”期间厦门附近海域平均波浪方向和有效波高（资料源自ERA-5）Fig. 5 Wave direction and significant wave height in coastal waters of Xiamen during the period of No.9914 Typhoon Dan on October 10,1999 (data from ERA-5)

台风大风、暴雨峰值与天文高潮相遇，厦门验潮站增水的幅度处于高位，出现风（风暴增水）、雨（暴雨）和潮（天文高潮）最大“三碰头”的情景。在三者叠加的作用下，厦门验潮站实测潮位达到过程峰值，厦门沿海出现了40 年一遇的超警戒极值水位。13 时后，台风登陆后减弱并北上，海水和潮波逐渐转为沿岸线方向传播，降水强度骤减，潮水退去，水位回落。

台风期间，1999 年10 月9 日11-12 时，除了上述“三碰头”的重要影响之外，向岸的强风、急流和波浪也同时发挥了显著的作用。调查表明，1999 年10 月9 日，在台风正面的袭击下，恰遇天文大潮，伴随狂风和暴雨，风暴潮挟卷着大量海水，巨浪冲击厦门南部尤其是东南部海岸并越过堤岸，造成了长达1 160 m护岸的损毁和沿岸街道的淹没，沙滩冲刷剧烈，岸线明显后退，防护林木折毁严重[30]。

综上所述，9914 号台风期间，厦门沿海地区发生的洪涝灾害，是海平面上升背景下风（向岸强风）、雨（强降水）、浪（巨浪）、潮（高潮位）、流（急流）等多致灾因子共同作用的结果。

3.3 1614 号台风期间海洋和大气过程对极值水位的影响

图6 为2016 年9 月14 日18 时至9 月15 日13 时，1614 号台风过程的降水、风暴增水和厦门验潮站潮位的变化。在1614 号台风登陆前（14 日18—23 时），厦门验潮站潮位不断升高，至23 时达到最高潮位，但期间的增水基本相同；15 日02 时风暴增水开始加强，台风于03 时登陆，风暴增水在05 时达到最大值。1614 号台风过程的明显降水始于14 日22 时，并逐步增大，台风登陆后3 h（03-05 时）的每小时降水量达到最强阶段，在04 时达到每小时最大降水量（58.9 mm）。

图6 1614 号台风“莫兰蒂”影响期间厦门海域天文潮位、风暴增水和降水过程Fig. 6 The astronomical tide, storm surge and precipitation in coastal waters of Xiamen during the period of No.1614 Typhoon Meranti

对比1614 号和9914 号台风期间的潮位、风暴增水和降水过程，可以发现有以下较明显的差异。一是1614 号台风过程的强降水时间短，降水强度大。1614 号台风的降水主要集中在台风登陆后的03-05 时，该时段厦门地区降水量激增，累计降水量将近120 mm。二是此次台风过程虽然也出现了明显的增水，但风暴增水的强度相对较小，且由于此次台风过程不处于天文大潮期，因此台风期间厦门验潮站并未出现超警戒的极值水位。值得注意的是，在台风登陆后2 h，厦门海域处于天文低潮期。前人的研究认为在落潮时段，天文-风暴潮耦合非线性项产生的是负增水效果[31]。但是，此时却有相反的现象，台风登陆后1 h（04 时）出现台风过程的最大增水值，1 h 后（05 时）风暴增水峰值出现。由此看来，强降水也是本次台风期间出现风暴最大增水的重要因素之一。

本文进一步计算分析了近20 年来影响厦门的风暴潮增水事件以及对应的降水数据的Kendall 相关系数。结果表明，Kendall 系数τ值为0.15。这表明，在厦门海域的风暴潮增水事件中，降水与风暴潮增水存在一定程度的正相关关系。结合本文分析的两个台风事件的具体降水、增水过程以及历史数据统计分析，可以推测强降水和风暴潮的叠加也是造成增水量突然增加的重要因素。

图7 为1614 号台风“莫兰蒂”登陆前后厦门附近海域海平面气压和风场的变化。相对于9914 号台风“丹恩”而言，1614 号台风“莫兰蒂”有着更大的影响范围、更大的强度，对海水水体有更强的驱动力和抽吸作用。在14 日20 时，厦门海域进入了台风的10 级风圈，海平面气压低于950 hPa，台风最大风力超过了50 m/s。在台风登陆的前两小时，相关海域风场情况与9914 号台风的情况几乎一致。厦门附近海域由原先偏北大风为主，逐渐向东南风向转变。在向岸强风以及台风抽吸作用下，验潮站已有50 cm 以上幅度的增水。而在风暴增水达到峰值的05 时前后，风场和气压场的配置并未出现明显变化。09 时过后，台风北移深入内陆，台风的风力明显减弱，低压中心填塞。

图7 2016 年9 月14-15 日1614 号台风“莫兰蒂”登陆前后厦门附近海域海平面气压和10 m 风场的变化（资料源自ERA-5）Fig. 7 The change of sea level pressure and the 10 m wind field in coastal waters of Xiamen during the landing time of No.1614 Typhoon Meranti from September 14, 2016 to September 15, 2016 (data from ERA-5)

另外，值得注意的是，1614 号台风“莫兰蒂”登陆前后虽未引起厦门沿海出现超警戒极值水位，但却造成厦门地区空前的森林和城市行道树的大量损毁[32-33]。

图8 为1614 号台风期间厦门附近海域的平均波浪方向和有效波高的分布。在1614 号台风的影响下，厦门附近海域的海洋环境场与风场也有着较为一致的变化关系。1614 号台风过程的潮流场与9914号台风过程的情况基本一致，但是在此次过程的台风登陆时刻，厦门海域的天文潮过程处于落潮阶段，因此，并未出现超高潮位的情况。由于此次台风的强度空前，近岸出现了较强的增水，厦门验潮站实测的最大风暴增水幅度达到1 m。15 日03 时台风登陆前后，在大风的驱动下，巨浪袭击岸线，厦门附近海域有效波高都在3 m 左右。在出现波浪越堤效应的同时，叠加有超强台风带来的短时强降水，这是造成此处沿岸低洼地淹没的主要原因。

图8 2016 年9 月14-15 日1614 号台风“莫兰蒂”期间厦门岛周边海域平均波浪方向和有效波高（资料源自ERA-5）Fig. 8 Wave direction and significant wave height in coastal waters of Xiamen during the period of No.1614 Typhoon Meranti from September 14, 2016 to September 15, 2016 (data from ERA-5)

3.4 RCP4.5、RCP8.5 情景下厦门验潮站极值水位和重现期预估

图9 为在温室气体中等和高排放（RCP4.5、RCP8.5）情景下厦门验潮站极值水位和重现期（相对于1954-1999 年）的变化。未来随着气候的持续变暖，全球海平面将进一步上升，许多地区沿海地区极值水位危险性将明显的增加[2-3]。为了预估未来厦门沿海地区极值水位危险性的变化，本文应用了1954-1997 年间厦门历史验潮站的观测数据，以及RCP4.5 和RCP8.5 情景下厦门海域海平面变化的预估值[21]，通过水文概率模型P-Ⅲ型线进行适线分析，获得了年极值水位的频率曲线（图9）。其中，极值水位的高度与其重现期呈对数线性关系。

图9 在RCP4.5(a)和RCP8.5(b)情景下未来厦门验潮站的极值水位事件和重现期（相对于1954-1999 年平均海平面）的变化Fig. 9 The relation between expected extreme sea level events and return period at Xiamen tidal gauge station, references to 1954-1999 mean sea level and future conditions for RCP 4.5 (a) and RCP 8.5 (b) scenarios

基于历史观测数据的概率曲线分析表明，在叠加不同情景下未来厦门沿海海平面变化的预估值（取中值）后，由图9（不同颜色实线）可以看出，在RCP4.5情景下，9914 号台风风暴潮引起厦门沿海地区发生40 年一遇的极值水位（732 cm）事件，到2030 年，732 cm水位的出现频率将提高到18 年一遇，在2050 年达6 年一遇，而到2100 年，将成为低于1 年一遇的事件。而在RCP8.5 的情景下，同期的频率比RCP4.5 情景更高，但差距不是很大。

表2 和表3 为不同气候（RCP4.5、RCP8.5）情景下当前厦门验潮站50 年一遇和100 年一遇的极值水位和重现期（相对于1954-1999 年）的变化。由表2可见，当前50 年一遇的极值水位为737 cm，在RCP4.5情景下，到2030 年，50 年一遇的极值水位将变为750 cm；到2050 年将变为762 cm；到2100 年将变为798 cm（取海平面变化预估值的中值）。在RCP8.5 情景下，到2030 年，当前50 年一遇的极值水位高度将变为751 cm；到2050 年将变为765 cm；到2100 年将变为822 cm。

表2 不同气候（RCP4.5、RCP8.5）情景下厦门验潮站50 年一遇和100 年一遇的极值水位的变化（相对于1954-1999 年）Table 2 1-in-50-year and 1-in-100-year extreme sea level events at Xiamen tidal gauge station, referenced to 1954-1999 for RCP(4.5, 8.5) scenarios

表3 不同气候（RCP4.5、RCP8.5）情景下厦门验潮站当前50 年一遇和100 年一遇极值水位的重现期的变化（相对于1954-1999 年）Table 3 The return period of recent 1-in- 50-year and 1-in-100-year extreme sea level events at Xiamen tidal gauge station, referenced to 1954-1999 for RCP (4.5,8.5) scenarios

另外，值得注意的是，当前厦门验潮站100 年一遇极值水位在RCP4.5 情景下，相对于1954-1999 年平均海平面，到2030 年将变为55 年一遇；至2050 年，将变为30 年一遇；到2100 年，将成为2 年一遇。而在RCP8.5 情景下，到2030 年，将变为53 年一遇；至2050 年，将变为25 年一遇；至2100 年，将低于1 年一遇。由此可见，未来厦门地区极值水位的危险性将大幅上升。

4 结论与讨论

4.1 结论

基于历史上9914 号、1614 号两次台风过程的观测和再分析数据，本文分析了台风登陆厦门地区前后，厦门海域的海洋和大气环境的变化，包括天文潮位、风暴增水和降水、海平面气压和风场等变化，影响厦门海域极值水位的各种因素和过程，以及不同气候情景下未来厦门海域极值水位和重现期的变化，得到以下几点结论：

（1）在1999 年10 月9 日9914 号台风“丹恩”登陆前的数小时内，厦门附近海域持续有向岸12 级大风（32.7 m/s），最大风速达到40 m/s，中心最低气压为965 hPa，厦门验潮站持续出现1.2 m 的风暴增水；至台风登陆时（12 时），每小时降水量达到最大值（29 mm），天文大潮达到最高潮位585 cm，厦门验潮站出现732 cm的极值水位，超过了700 cm 的警戒水位。

（2）9914 号台风影响期间，在持续向岸强风和台风中心低压的驱动、抽吸以及涨潮流的作用下，厦门附近外海海域的海水不断向厦门海域内输运堆积，强风、低压和涨潮流对风暴增水起到了正面的驱动作用。天文大潮最高潮位、风暴最大增水和最大降水量三者的同时出现与共同作用造成了异常偏高极值水位（732 cm）的出现。台风正面袭击厦门地区时，恰遇天文大潮，风（向岸强风）、雨（暴雨）、浪（巨浪）、潮（高潮位）、流（急流）造成的风暴潮是厦门沿海地区发生人员伤亡、护岸损毁、街道淹没，以及大量经济损失等严重灾情的重要原因。

（3）在不同气候（RCP4.5、RCP8.5）情景下，到21 世纪中叶和末期，当前厦门验潮站50 年和100年一遇极值水位将变得很频繁，其危险性将显著上升。例如，在RCP4.5 情景下，到2050 年和2100 年，当前百年一遇的极值水位将变为30 年一遇和2 年一遇；在RCP8.5 情景下，到2050 年和2100 年，当前百年一遇的极值水位将变为25 年一遇和低于1 年一遇。

（4）在不同气候（RCP4.5、RCP8.5）情景下，未来全球和我国沿海海平面将持续上升，强台风将趋多，这将抬高极值水位的基础水位[3]，因此，我国沿海极值水位将有显著增高的趋势，极值水位的重现期将缩短，但具体的变化将有待我们进一步深入预估。可以预计的是，我国东南沿海包括厦门地区如果仅维持现有的防洪排涝标准和基础设施，海岸洪涝灾害风险将极可能显著增加。因此，为降低未来的沿海地区洪涝灾害风险，应加强研究并采取充分的适应气候变化措施。

4.2 讨论

Emanuel[34]模拟研究表明，到21 世纪末，在登陆前迅速增强的热带气旋可能将会大幅增加。Kossin[35]的研究发现，在全球平均温度上升0.5℃, 全球热带气旋的移动速度减缓约10%，尤其是在西北太平洋地区，减幅达到了21%。这表明未来热带气旋（台风）登陆期间的风暴增水和降水可能比现在的预估更高。本文利用了验潮站的历史观测资料和未来极值水位重现期的计算方法，分析的未来海平面变化情景下极值水位和重现期的变化，在一定程度上反映了不同气候情景下未来极值水位事件危险性（强度和重现期）的变化，但尚未更充分考虑未来强热带气旋（强台风）变化情况，因此，还有待将此情况一并纳入预估分析。

致谢：衷心感谢王慧研究员和王彥明博士对本文工作的大力支持和帮助。