Marine Meteorology

Reconstruction of the tropical cyclones activity in the Nansha Islands since the Little Ice Age from the atoll lagoon sediments

  • YANG Hongqiang , 1, 2, 3 ,
  • TAN Fei 1, 2, 4 ,
  • XU Huilong 1, 2 ,
  • ZHANG Xiyang 1, 2 ,
  • SHI Qi 1, 2 ,
  • TAO Shichen 1, 2
Expand
  • 1. Key Laboratory of Ocean and Marginal Sea Geology, South China Sea Institute of Oceanology, Innovation Academy of South China Sea Ecology and Environmental Engineering, Chinese Academy of Sciences, Guangzhou 510301, China
  • 2. Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511458, China
  • 3. Nansha Marine Ecological and Environmental Research Station, Chinese Academy of Sciences, Sansha 573199, China
  • 4. University of Chinese Academy of Sciences, Beijing 100049, China
YANG Hongqiang. email:

Copy editor: LIN Qiang

Received date: 2022-01-27

  Revised date: 2022-04-14

  Online published: 2022-04-18

Supported by

Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou)(GML2019ZD0206)

Innovation Academy of South China Sea Ecology and Environmental Engineering, Chinese Academy of Sciences(ISEE2018PY01)

Strategic Priority Research Program of the Chinese Academy of Sciences(XDA13010103)

Abstract

Tropical cyclones (TCs) and the resultant storm surges and heavy rainfall have posed huge socio-economic threats to the coastal areas of the South China Sea (SCS) and its surrounding areas. Holocene TCs studies before instrumental record help to accurately predict trends in TCs activity in the context of global warming. Here, we presented a new sub-annually resolved paleostorm record from the Anle reef lagoon sediments of the Nansha Islands in the southern SCS since the relative cold Little Ice Age (LIA) and a total of 28 storm event layers were identified. The new reconstruction indicates that the site has witnessed frequent TCs variations on a decadal to centennial-scale since the LIA, and two main storm active periods have been developed. The most intense storm activity period occurred in the early LIA (AD 1471~1620), and the other relatively active period was witnessed in the Current Warm Period (CWP) of AD 1930~1960. Although the storm activity in CWP was strengthened, it was significantly lower than that in the early LIA. Regional comparison with the storm reconstruction record of the Yongshu Reef during the same period reveals distinct spatio-temporal heterogeneity in the southern SCS. More high-resolution storm records from similar or adjacent regions could effectively reduce the uncertainty of storm reconstructions and increase the accuracy of reconstructed records.

Cite this article

YANG Hongqiang , TAN Fei , XU Huilong , ZHANG Xiyang , SHI Qi , TAO Shichen . Reconstruction of the tropical cyclones activity in the Nansha Islands since the Little Ice Age from the atoll lagoon sediments[J]. Journal of Tropical Oceanography, 2022 , 41(6) : 171 -182 . DOI: 10.11978/2022016

极端热带气旋活动对世界各地沿海生态系统和城市基础设施构成了重大威胁(Peduzzi et al, 2012)。西北太平洋海盆是地球上遭受热带气旋灾害最为频繁和剧烈的区域, 每年热带气旋个数约占全球气旋活动总量的30%左右(Mei et al, 2015; 廖菲 等, 2019; 陶丽 等, 2020), 给中国大陆东南沿海地带及南海周边带来巨大的经济和渔业损失(Yang et al, 2015; Soria et al, 2016; 张敏 等, 2019)。现代观测记录和气候模型预测结果均表明在全球变暖背景下热带气旋活动可能会向更强和更频繁的趋势发展(Knutson et al, 2015; Bacmeister et al, 2018)。相对短暂的器测记录阻碍了我们对热带气旋长期变化趋势的认识(Landsea et al, 2006; Toomey et al, 2013; Chen et al, 2019b; 高抒 等, 2019)。利用沉积物重建的长尺度古风暴记录, 可为全球热带气旋的预测提供重要参照(Knutson et al, 2015; 周亮 等, 2015; Liu et al, 2016; 高抒 等, 2019; 徐笑梅 等, 2019; Ishizawa et al, 2020)。大洋型珊瑚环礁潟湖中的粗粒沉积物可用于古风暴或海啸事件的判定(Yu et al, 2009; Toomey et al, 2013; Yue et al, 2019)。环礁潟湖沉积物具有沉积速率高、成岩蚀变小、可精确定年的特点, 可发育连续的高分辨率风暴记录。南海珊瑚礁区是研究百年~千年尺度古风暴的理想区域(Yu et al, 2009)。研究表明小冰期南海地区遭受了强烈且频繁的热带气旋活动(Liu et al, 2001; Chen et al, 2019a; Sun et al, 2020)。然而, 南海南部热带地区的高分辨率古风暴活动记录目前还较少, 对南海小冰期热带气旋活动的时空差异性、活动变率和驱动机制的认识不足。
本研究利用南沙群岛安乐环礁礁后潟湖沉积物中的粗粒组分, 重建该区小冰期以来的热带气旋活动, 进一步证明小冰期早期南海南部处于强烈的热带气旋活动期, 并与附近的永暑礁潟湖风暴记录进行对比, 探讨了区域热带气旋活动的时空差异及可能机制。

1 材料及方法

1.1 研究区位置和样品采集

南海位于亚洲—澳大利亚季风的中心转换地带, 同时受到太平洋和印度洋海气作用的相互影响(Contreras-Rosales et al, 2019)。南海不仅受到西太暖池的热带气旋的影响, 而且也是西太平洋热带气旋的重要起源区(Zuki et al, 2008)。研究区位于南沙群岛的九章环礁中的安乐礁, 为大洋型珊瑚环礁(图1)。据美国国家海洋和大气局 (National Oceanic and Atmospheric Administration, NOAA)的 IBTrACS Hurricane Dataset的热带气旋路径数据库, 在1945—2020年, 共有75次热带气旋活动登录南沙群岛附近区域, 其中大于Cat 1级别强度的风暴事件达17次(图1b)。安乐环礁位于九章群礁东北部, 很少受到人为活动的干扰。安乐礁形似耳状, 从西北部延伸到东南部, 是一个几乎完全封闭的珊瑚环礁, 仅由南部的几个潮汐通道与外海和九章环礁的深潟湖相连通。环礁内为障壁礁后潟湖, 水深6~30m。在天气温和晴朗时, 潟湖沉积物很少受到礁坪外围波浪及潮汐作用的影响(图1c)。卫星观测数据和南沙永暑礁验潮站仪器观测资料显示, 研究区多年平均潮差为0.5~1.4m, 潮汐对礁后潟湖沉积物的影响较小。强热带风暴或台风引起的波浪或风暴潮超过2m, 风暴潮可将粗粒沉积物从周围的宽缓礁坪输送到礁后潟湖。因此, 安乐礁潟湖内部的粗颗粒沉积物能够记录热带气旋活动, 可以用来重建热带气旋活动规律和变化趋势。
图1 研究区位置及气象和地貌学特征

a. 南海和西北太平洋地理位置及热带气旋历史记录, 该图依据审图号GS(2016)2937号底图制作, 黄色圆圈标注为本研究中的采样位置安乐礁潟湖, 永暑礁NY-4风暴记录以蓝色圆圈表示; b. IBTrACS Hurricane Datase数据库中提供的南沙群岛附近历史热带气旋移动路径(该图来源于 https://coast.noaa.gov/hurricanes/, 图中为热带风暴强度及以上级别热带气旋路径, 白色虚线圆圈代表以安乐礁为中心半径为300km范围内的热带气旋活动路径); c. 安乐礁位置和地貌分区(该图根据google卫星遥感图像制作, 黄色圆圈代表岩心AL-1采集地点, 黄色虚线表示礁坪和潟湖之间的边界)

Fig. 1 Study location and the regional atmosphere-geomorphology context

2017年9月至10月期间, 研究团队在安乐礁礁后潟湖中(9°55'56.75"N, 114°31'19.55"E)利用内径9cm的冲击不锈钢管采集了长度约5m的潟湖沉积物岩心(图1c), 岩心采收率为100%, 确保了垂向地层的完整性。取心地水深12m, 为砂质的潟湖底环境, 该环境不利于原生大型造礁石珊瑚丛体生长, 因此该地的珊瑚丛体断枝或大的珊瑚砾块沉积物通常为高能风暴浪异地搬运形成。

1.2 实验室粒度分析

在野外提取潟湖沉积物岩心后, 岩心随后被运往中国科学院大洋与边缘海地质重点实验室进一步处理分析。在室内首先对岩心进行纵向剖分、分离、拍照和沉积物定性描述。根据目视初步观察, 安乐礁潟湖岩心的背景沉积物为粒径小于0.25mm的细粉砂—细砂质细粒沉积, 基本无陆源沉积物输入; 顶部和底部含多层可能由风暴事件作用引起的粗粒沉积层, 包含极粗砂、珊瑚断枝或砾块, 分布杂乱, 可见清晰的棱角和完整生物贝壳。为了确定沉积物粗颗粒的百分比, 参考Bramante等(2020)等的粒度分析流程, 以1cm等间距对岩心进行了连续取样, 总共采集了496个样品用于粒度分析。然后使用250、500、1000和2000μm孔径的样品筛对分散的干样进行筛分并称量不同粒级的干重, 再使用63μm和125μm孔径的样品筛对250μm以内的样品进行湿筛, 并对各部分筛后样品烘干称重, 以获取每个粒级的重量百分比。
本研究将潟湖中的粗粒沉积物作为异常强水流引起的沉积物搬运的识别标志。沉积物颗粒大小与沉降速度呈正相关, 地表水流的强度决定了能移动的最大颗粒大小和距离。在正常平缓的水流下主要是沉积细粉砂, 而风暴驱动的快速水流可以搬运和沉积厚的粗砂层。一般情况下沉积颗粒沉降速度随粒径的增大而单调增加, 颗粒粒径可作为异常水流搬运的替代指标(Bramante et al, 2020)。对于球形和近球形颗粒, 颗粒直径和沉降速度之间存在非常强的关系(Dietrich, 1982)。但在环礁潟湖中产生的许多粗粒碳酸盐生物碎屑沉积物, 特别是对于大于2mm的生物碎屑, 通常形状不规则, 且可能含有大量的孔隙, 导致搬运阈值和沉降速度低于实际粒径的预测值。如果这类颗粒构成沉积物粗粒组分的大部分, 那么沉积物颗粒与沉积物搬运、沉降和水流强度关系的推断可能无效, 导致生物碎屑沉积物粒度和沉降速度之间的关系存在偏差(Maiklem, 1968; Braithwaite, 1973; Kench et al, 1996)。参考Bramante 等(2020)在太平洋贾鲁特环礁蓝洞沉积物风暴事件重建中使用的方法, 本文根据潟湖沉积物粒径的总体分布特征, 将沉积物粒度分为0.25~2mm的粗组分和小于0.25mm的细组分。最后根据粗组分颗粒(0.25~2mm)的百分比来识别风暴事件, 大于2mm的极其粗颗粒并未纳入粗组分百分比计算。本处理方法可以有效降低一些异常的粗粒峰。在这些沉积层中通常含有2个粗组分峰, 孤立和异常粗大的生物碎屑, 如大贝壳或者珊瑚砾块产生了大部分未校正的一个粗粒峰, 另一个粗峰为自身沉积物构成, 但常被>2mm的粗峰给掩盖掉。在实际处理中, 我们仅将这种>2mm的粗粒峰作为风暴层识别的参考和补充。

1.3 年代地层框架建立

风暴事件精确年龄的确定对古风暴活动重建至关重要。同一层位中纯净的鹿角珊瑚(Acropora spp.)断枝的U/Th年龄代表该层最小年龄, 是建立潟湖沉积物年代框架最为可靠的材料(Yu et al, 2009; Chen et al, 2019b)。本研究共选取21个无成岩改造的鹿角珊瑚断枝进行U/Th定年, 对潜在的风暴粗粒层进行加密采样测试。
U/Th测年处理流程和分析方法遵循Clark等(2014)提出的方法。首先将每个珊瑚断枝样品的表面抛光, 去除表面的藻类和生物有机质膜。然后用玛瑙杵或刻刀将珊瑚骨骼粉碎到颗粒大小为~0.5mm的颗粒。将破碎后的样品浸泡于预先清洗干净的含有15% H2O2的玻璃烧杯中48h以去除有机残留物。然后用MQ水超声清洗样品多次, 直到水变澄清。然后将多余的水倒出, 样品在40℃下充分干燥。在双目显微镜下检查每个样品的纯净程度, 最终选择约150mg的纯净珊瑚骨骼进行年代测定。在年代测定之前, 对前处理后的样品进行XRD射线衍射分析, 以检测文石质珊瑚骨骼是否发生成岩蚀变。在澳大利亚昆士兰大学地球与环境科学学院的放射性同位素实验室的 Nu Plasma MC-ICP-MS上进行珊瑚断枝的 U/Th年龄测定。U/Th比值是利用 Isoplot/EX 3.0软件程序, 基于Cheng等(2000)提出的衰变常数计算出的活度比值。所有年龄数值扣除实验室平均空白, 测试误差均表示为2σ标准偏差。岩心顶部12cm沉积物由于在取样过程中受到搅动, 未进行粒度分析。使用BACON建模软件建立了高精度的年龄框架, 并利用贝叶斯统计方法计算了岩心每1cm间隔区间的加权平均年龄和年龄不确定性(Blaauw et al, 2011)。

1.4 古风暴事件层阈值和频率

本文通过沉积学观察和粒度分析确定风暴事件层, 并将两种方法识别的风暴层进行交叉验证。安乐礁潟湖岩心的风暴层与上下的正常沉积层突变接触, 颜色较深, 发育更粗粒的颗粒或珊瑚断枝和贝壳等, 厚度1~20cm不等。因为事件层几乎是瞬间形成的, 在年龄模型中, 本文将厚度大于2cm的事件层校正为1cm。采用Bramante等(2020)和Wallace等(2019)提出的方法确定风暴事件层的粒度阈值, 以区分事件层和正常的背景沉积。首先, 对整体粗组分序列进行11点滑动平均滤波处理, 粗组分超过平均粗组分含量1σ的样品定义为异常值(Outliers)(>10%), 较厚的事件层对应于连续的异常值, 在2σ时间误差范围对其进行深度校正, 代表一次强烈的风暴事件层。然后, 将排除异常值后的粗组分数据再次进行11点滑动平均滤波处理。最后, 我们将超过滑动平均滤波后粗组分百分比1.5σ(2.35%)的粗峰定义为风暴事件层。从分析中去除异常值可以防止特别大的事件层掩盖附近较小事件峰的存在, 使用11点滑动窗口可以消除背景沉积中粗组分信号的年代际变化的影响(Donnelly et al, 2015; Bramante et al, 2020)。
为评估该方法对研究区热带气旋的响应能力, 岩心顶部的事件沉积与现代热带气旋观测数据进行了校准。从IBTrACS version 4收集了AD 1945—2020期间所有过境研究区附近的热带气旋活动信息, 将超过风暴阈值的粗粒层与可能影响安乐礁潟湖的现代热带气旋记录进行比较, 评估事件层阈值的可信度。在最终确定了风暴事件层阈值和岩心中全部的事件层数量后, 参考Lane等(2011)的方法, 使用100年窗口滑动重新计算了研究区每百年风暴事件的发生频率, 以评估百年尺度的热带气旋变化特征。50年滑动窗口的分析用于揭示年代际周期内高频或低频的热带气旋活动。
使用1945—2020年IBTrACS的热带气旋最佳路径数据集经过安乐礁潟湖115km范围内的Cat 1级别或大于Cat 1热带气旋级别的频率来计算研究区附近局地热带气旋估计频率(热带风暴和Cat 1~Cat 5不同热带气旋级别强度的详细定义见 https://www.nhc.noaa.gov/aboutsshws.php)。研究区局地热带气旋范围(115km)划定是参考Rodysill等(2020)根据巴哈马现代飓风事件和碳酸盐沉积记录及水动力模拟得出的结果。以安乐潟湖为中心, 统计研究区300km半径范围内的Cat 1级别及更强热带气旋的频率(基本涵盖了整个南沙群岛), 获得南海南部区域热带气旋估计频率。比如, 根据IBTrACS观测记录, 安乐礁附近的热带气旋记录始于AD 1945, 在AD 1945—2020期间, 在半径 115km的范围内登陆或者直接影响安乐礁强度大于或等于Cat 1级别热带气旋的数目为5个, 据此每百年的当地热带气旋频率λ=5÷75×100=6.7(小数点后保留1位有效数字)。本文使用局地热带气旋估计频率来界定安乐礁潟湖的风暴相对活跃期和平静期, 当某个时期100年滑动窗口风暴事件超过该区的局地热带气旋频率λ界限即定义为风暴活跃期。

2 研究结果

2.1 年代框架

基于21个纯净的鹿角珊瑚断枝高精度U/Th年龄建立了AL-1的年代地层框架(图2c)。整体上, U/Th年龄随着沉积物深度增加而增加, 表明潟湖沉积序列的连续性。79cm和187cm处的铀系年龄出现倒转, 被剔除(图2c)。根据年代模型, AL-1记录在小冰期以来的记录中, 其沉积速率为0.86cm·a-1, 即1cm沉积物的分辨率为~1.37a, 是目前南海区域分辨率较高的古风暴记录之一。根据沉积速率的变化, AL-1可大致分为下部(496~234cm)、中部(233~153cm)、上部(152~0cm)三部分。下部234cm(校准深度后为209cm)为AL-1初始沉积速率较快的部分, 沉积速率为1.15cm·a-1。234cm到153cm(校准深度为128cm)的中部, 沉积速率降低至0.33cm·a-1。而岩心上部沉积速率增加到1.42cm·a-1。岩心中粗组分含量与沉积速率的变化一致, 岩心底部和顶部的粗组分含量相对较高, 中间层段粗组分含量相对较低(图2a)。事件层粗组分含量变化范围为6.1%~52%, 平均为11%。事件层整体以灰白色中—粗砂为主, 且常含粒径>2mm的珊瑚断枝和生物壳体, 可与正常沉积层的浅灰色泥质—细砂质碳酸盐沉积区别。
图2 安乐礁潟湖AL-1岩心沉积物剖面和U/Th年龄模型

a. AL-1粗组分含量剖面, 红色圆圈代表超过粗组分含量百分比11点滑动平均值的1.5σ的风暴事件层; b. 原始粗组分含量>10%的异常值(Outliers); c. 基于鹿角珊瑚U/Th年龄建立的岩心AL-1的BACON年代模型

Fig. 2 Profiles of sediment core AL-1 collected from the Anle lagoon and U/Th age model

2.2 历史风暴事件频率

研究区热带气旋活动频率在年代际到百年尺度上变化很大。AL-1岩心中共识别出了28层事件沉积层, 其中底部含有17个事件层, 中部3个, 顶部8个(图3)。本文使用经过安乐礁附近115km范围内的大于或等于Cat 1级别的每百年区域热带气旋频率(λ=6.7)来定义风暴活动期。风暴事件频率的100年滑动窗口分析结果显示, AL-1岩心记录了两个主要风暴活跃时期。持续时间最长的风暴活跃期(A1)发生于小冰期(Little Ice Age, LIA)早期(AD 1421— 1610), 平均每百年风暴事件数为9.2个 (图3), 尤以AD 1510—1550期间热带气旋活动最为频繁, 每百年平均为12个, 几乎为现代区域热带气旋频率的2倍。随后从AD 1610—1890, 安乐礁潟湖热带气旋活动趋于平静, 热带气旋频率呈下降趋势, 最平静时期为AD 1810—1890期间, 每百年热带气旋频率为0, 没有任何沉积记录超过风暴事件阈值。另一个相对风暴活跃期(A2)位于AD 1930—1960, 平均每百年热带气旋数为7.5个。50年滑动窗口风暴事件频率揭示了年代际的古风暴极高活跃期, AD 1496—1529和 AD 1534—1557期间热带气旋事件频率分别达到每50年6.1个和7.0个。整体上, 50年滑动窗口分析显示了和100年滑动窗口分析类似的热带气旋活动变化规律, 即在小冰期早期的大约2个世纪期间(AD 1400—1600), 研究区热带气旋活动频繁且强烈, 此后研究区热带气旋活动频率和强度显著减弱。在不同的滑动窗口中, 小冰期晚期是热带气旋活动最不活跃的时期。AD 1931—1966在50年滑动窗口中则表现为现代热带气旋最活跃期, 基本和百年滑动窗口同步。
图3 AL-1岩心粗组分异常和风暴事件层数百年滑动窗口

红线代表粗组分含量的11点滑动平均值(去掉>10%的异常值), 绿线代表1.5σ的阈值, 超过此阈值则代表一次风暴事件层; 黑色五角星代表一次风暴事件; 蓝线代表AL-1热带气旋事件频数的100年滑动窗口数, 黑色虚线代表基于IBTrACS v4观测数据得到的局域热带气旋估计频(λ=6.8)

Fig. 3 The Coarse sediment anomalies and 100-year sliding window on the number of storm event layers recorded by core AL-1

3 讨论

3.1 风暴事件识别和现代记录校准

评估地质代用指标适用性和确定风暴事件阈值有助于准确重建长期的古风暴变化历史(Donnelly et al, 2015; Winkler et al, 2020)。礁后潟湖是主要堆积区, 其沉积物主要来源于礁坪和礁脊, 潟湖中的原生点礁亦可作为物质来源。潟湖水动力从礁脊到潟湖盆地逐渐减弱, 粒度逐渐变细, 沉积物以粉砂—细砂质沉积为主。正常天气条件下, 波浪能量通常在珊瑚礁100~500m消散(Kench et al, 2006), 限制了粗颗粒珊瑚礁沉积物进入潟湖。风暴期间的大浪和风暴潮会冲刷 浅的礁脊边缘和礁坪, 将珊瑚碎片、断枝、珊瑚藻和底栖有孔虫等粗粒骨骼碎屑搬运到富含砂质沉积的礁后潟湖(Toomey et al, 2013; Klostermann et al, 2014)。对安乐环礁的实地调查发现活体Acropora spp. 是礁坪的优势物种之一, 其骨骼相对脆弱, 易受到风暴潮的影响, 破碎后被迅速搬运和堆积在潟湖。岩心AL-1距离礁坪约600m(水深约12m), 没有活珊瑚生存。堆积在粗粒事件层中的珊瑚枝大多呈白色和灰白色, 保存有明显的棱角和新鲜断面, 表明这些珊瑚枝是被迅速搬运和沉积的。因此, 低能潟湖环境中较粗沉积物可作为识别热带古风暴活动的标志(Yu et al, 2009; Yue et al, 2019)。不同粗组分的百分含量可以代表古风暴强度。
热带气旋强度、路径、平移速度等参数均可影响风暴沉积(Winkler et al, 2020; Rodysill et al, 2020)。在佛罗里达州的Apalachee海湾, 水动力学模拟表明, 尽管Cat 4级别和Cat 5级别的飓风产生了最极端的强风暴潮, 但中等强度Cat 1~ Cat 3级别的飓风也能够产生强风暴潮, 这与风场、热带气旋路径、移动速度或更接近的路径等因素有关(Lin et al, 2014; Winkler et al, 2020)。风暴沉积可能发生在风暴频繁的时期, 一个粗粒层可能记录了多个风暴。我们将安乐礁潟湖附近现代热带气旋移动路径的观测数据与岩心中的粗颗粒事件层进行比较, 以评估岩心对不同热带气旋记录的潜力, 并验证事件阈值的可靠性。
在1945—2020年期间, 一共有2个Cat 1级别及以上强度的热带气旋经过安乐礁潟湖60km半径范围内, 而热带风暴强度级别及以上的热带气旋数为10个。其中, 在1962年12月28和29日, Cat 3级别的热带气旋Lucy在115km的范围内, 以平均风速大于46.3m·s-1, 最大风速超过51.7m·s-1的强度持续影响安乐礁超过48h, 最近距离研究区不到30km, 最有可能在研究区内留下事件沉积的证据。虽然热带风暴相对于Cat1~Cat5级别的热带气旋强度更低, 但受热带风暴移动路径与距研究区的距离、风暴持续时间等因素影响, 也可能形成风暴事件层, 比如热带风暴Tess (1988)移动路径最近距安乐潟湖不到10km, 最大风速超过30.9m·s-1。在半径115km的范围内过境安乐潟湖的热带风暴及以上级别的热带气旋共有16个, Cat 1级别及以上的热带气旋为5个, Cat 2级别及以上有3个(表1)。
表1 1945—2020年登陆南沙群岛安乐礁和永暑礁不同半径范围内的不同强度热带气旋的频率

Tab. 1 The historical tropical cyclones passing within different distances from the Anle Reef and the Yongshu Reef in the Nansha Islands, southern SCS between 1945~2020

珊瑚礁区 半径范围 热带风暴及以上级别热带气旋 每百年频率 ≥cat 1级别热带气旋 每百年
频率
≥cat 2级别热带气旋 每百年
频率
≥cat 3级别热带气旋 每百年
频率
安乐礁 r<300km 75.0 100.0 17.0 22.7 9.0 12.0 3.0 4.0
r=115km 16.0 21.3 5.0 6.7 3.0 4.0 1.0 1.3
r=60km 10.0 13.3 2.0 2.7 1.0 1.3 1.0 1.3
永暑礁 r<300km 67.0 89.3 15.0 20.0 8.0 10.7 1.0 1.3
r=115km 18.0 24.0 6.0 8.0 2.0 2.7 1.0 1.3
r=60km 6.0 8.0 2.0 2.7 1.0 1.3 1.0 1.3
本文根据岩心顶部沉积物的粗组分含量和事件阈值建立了AD 1900以来的现代风暴沉积序列, 并将其进一步与IBTrACS观测值进行校准。从AD 1900—2000, 一共识别出11个事件层(图4)。考虑到沉积特征、粗组分含量和多个事件层的年龄位于2σ年代误差范围内, 将事件层E8附近的4个粗粒峰合并为一个事件层, 该事件层的年龄定义为粗组分含量最高层对应的年龄。最终识别出以下8个事件层。
图4 安乐礁现代器测热带气旋移动路径和沉积物记录敏感性校准

a. 经过南海南部安乐礁115km半径范围内的热带气旋活动历史路径; b. 60km范围内近距离影响安乐礁的热带气旋活动路径; c. 岩心顶部1900—2001 AD的粗组分含量百分比与现代热带气旋观测记录的校准和敏感性评估。蓝色五角星代表超过事件阈值的粗粒沉积峰, 其中事件层E8为附近几个小峰的合并; 不同颜色和大小的圆分别代表了某一热带气旋过境安乐礁附近115km范围内的强度和距离安乐礁的最短直线距离, 圆半径越大代表距离安乐礁越近。每个事件层标注了基于贝叶斯模型的2σ年龄误差。现代热带气旋移动路径观测记录来源于 https://coast.noaa.gov/hurricanes/在线数据库网站, 统计时间区间为1945—2020年

Fig. 4 Sensitivity calibration of the sediment records to the modern instrumented tropical cyclone tracks passed near the Anle Reef

事件层E1: 沉积于AD 1971—2004(中值年龄为AD 1987), 最可能由Cat 2级别的热带气旋Faith造成, Faith在1998年12月12日在~55km的范围内自东北向西南向袭击了安乐礁。鉴于沉积物记录无法区分同一校正日历年误差范围内的事件, 热带风暴Tess无法被排除, 其在1988年11月5日登陆安乐礁东北方向, 最短距离安乐潟湖<10km。在此时间范围内的其他热带风暴的强度和风暴距离不足以产生风暴事件层。
事件层E2: 沉积于AD 1969—2001(中值年龄为AD 1984), 鉴于非常接近的年龄和较近的风暴距离, E2很可能是热带风暴Tess造成的, 其在1988年11月5日, 以距离安乐潟湖小于10km的距离经过安乐礁东北方向。其他热带风暴如Mamic(1982年3月22日, 距离安乐礁北110km)强度太小和距离太远不足以形成事件层。
事件层E3: 沉积于AD 1962—1989(中值年龄为AD 1973), 则可能是距离相对较远的Cat 1级别的热带气旋Freda(1967年11月7—8日, 距安乐礁~90km)或者近端(<30km)的热带风暴Thelma导致(1973年11月14日, 距离安乐礁南~27km)。
事件层E4: 沉积于AD 1954—1972(中值年龄为AD 1961), 是岩心顶部最为明显的粗粒事件层, 该事件层包含了41~61cm, 厚度达20cm的连续沉积。E4的粗组分含量平均为27.80%, 最大为52%, 并含有珊瑚砾石, 与粉砂质碳酸盐背景沉积物混杂沉积。Cat 1级别的热带气旋Freda和Cat 3级别热带气旋Lucy(1962年11月28—29日, 距离安乐礁北~29km)在此期间都曾登陆安乐礁。但是, 根据事件层的强度和厚度, 过境风速超过51.7m·s-1的热带气旋Lucy是导致此次厚层风暴事件沉积的最可能原因。距离热带气旋移动路径短、强度大、登陆过境时间长等因素最终在安乐礁产生了最强烈的风暴事件。
事件层E5(沉积于AD 1949—1961, 中值年龄为AD 1955)和事件层E6(沉积于AD 1936—1956, 中值年龄为AD 1947)很可能由Cat 2级别的热带气旋Amy导致, Amy于1951年12月16—17日过境安乐礁潟湖北部, 最近距离安乐礁~48km, 最大风速为43.7m·s-1。考虑到事件层E6仅稍微超过风暴阈值, 我们也可以合理地认为它们与近距离的, 但强度相对较弱热带风暴有关, 但是热带风暴Betty(1949年12月7日, 距安乐礁北52km)的强度和过境持续时间不大可能造成事件沉积。AD 1970之前的热带气旋轨迹和强度记录存在较大不确定性, 加之定年误差, 使特定风暴和粗异常事件层之间的相关性降低(Winkler et al, 2020)。事件层E7和E8分别沉积于AD 1941和AD 1934, 而IBTrACS数据库中南沙群岛最早的观测记录始于AD 1945, 无法确定这两个事件层和特定热带气旋的对应关系。
综上, 风暴事件层的年龄大体对应过境研究区115km范围内的Cat 1级别热带气旋和更强级别的强风暴事件。因此, 粗粒沉积层可以表征风暴事件, 并且达到相似风暴事件阈值的沉积记录能可靠地捕获过去对应的风暴记录。考虑到热带气旋最佳路径数据最大风速估计的不确定性, 早期热带气旋移动路径和影响可靠半径的缺失, 本文粗略的估计还需要潟湖风暴流体动力学模型的验证。

3.2 南海南部区域古风暴频率变化及比较

对比南海不同地区的研究有助于深入了解过去整个区域或泛盆地范围的热带古风暴活动。本文将安乐礁潟湖的古风暴记录与南海南部的其他风暴重建记录进行比较。南海南部长度超过千年的古风暴重建记录比较稀缺, 目前仅有的记录来自研究区西南200km的永暑礁潟湖沉积岩心NY-4(Yu et al, 2009)。NY-4利用潟湖沉积物中大于1mm的粗颗粒峰重建了4000余年的风暴记录。参照安乐礁潟湖附近的局地热带气旋估计频率, 本文计算了以永暑礁为中心, 115km半径范围内的热带气旋估计频率(表1)。在和AL-1记录同区间内, NY-4一共记录了32个风暴事件, 和AL-1风暴数量大致相当(28次)(图5)。NY-4古风暴频率100年和50年滑动窗口显示出明显的差异(图5c、d)。在100年滑动窗口中, NY-4记录了4个持续10~20年短周期的风暴活跃期(分别位于AD 1694—1703、1707—1739、1741—1754和1823—1832)和一个位于AD1768— 1818期间持续长50年的风暴活跃期。考虑到NY-4 整个岩心的2σ分辨率误差为18a·cm-1, 无法准确地将低于此年代误差的事件层单独分离出来, 大体认为NY-4记录了在AD 1664—1832期间持续168年的百年尺度风暴活跃期。50年滑动窗口结果则揭示了两个分别位于AD 1665—1723和 AD 1758—1837的风暴活跃期, 这两个活跃期和100年滑动窗口可以很好的对应。但是, 与100年滑动窗口相反, 在AD 1714—1757之间表现为一段风暴平静期。鉴于整个NY-4岩心年龄的不确定性, 无法获悉50年滑动窗口尺度上明显的风暴短暂减弱是由较短的年代际外部气候驱动, 或者表征风暴活动的随机性。
图5 小冰期以来南海南部安乐礁和永暑礁潟湖的古风暴事件频率对比

a. 安乐礁古风暴事件数量100年滑动窗口; b. 安乐礁古风暴事件数量50年滑动窗口; c. 永暑礁古风暴事件100年滑动窗口数; d. 永暑礁古风暴事件50年滑动窗口数(Yu et al, 2009)。滑动窗口数中的黑色虚线代表该经过该区域115km范围内的局地热带气旋估计频率。50年滑动窗口中的估计频率按照100年滑动窗口频率的50%计算

Fig. 5 Comparison of the paleostorms events recorded in the Anle and Yongshu Reef Lagoon, southern South China Sea since the LIA

值得注意的是, 与AL-1岩心记录不同, NY-4在小冰期中晚期(AD 1664—1832)发育多个风暴事件层, 而在AL-1记录中, 风暴最活跃时期位于小冰期早期。记录的不一致性可能反映了影响两个地点风暴数的不同以及特定地点对风暴活动记录的敏感性差异。例如, 从1945—2020年, 一共6次Cat 1级别或以上的热带气旋过境永暑礁115km的半径范围, 然而NY-4记录中仅记录到了AD 1975的一次风暴事件层沉积(表1)。与之对比, AL-1记录了至少5次风暴事件层。在1962年, Cat 3级别的热带气旋Lucy(风速>51.7m·s-1)过境安乐礁潟湖和永暑礁115km范围, 在安乐礁沉积了厚粗粒事件层, 永暑礁却没有留下明显的风暴沉积证据。类似的现象也见于Cat 3级别的热带气旋Faith(1998年)和热带风暴Tess(1988年)。
热带气旋活动的随机性也是影响其移动路径和登陆强度的重要因素之一(Wallace et al, 2021a, b)。由于研究点的沉积记录只捕捉到近距离经过的风暴, 因此这些风暴可能仅影响过境岛礁, 却对其他区域的影响有限。来自巴哈马Long岛(Wallace et al, 2021b)和佛罗里达西北部Apalachee 湾(Lin et al, 2014)的热带气旋研究证实了上述论断。研究指出, 热带气旋强度、热带气旋性质(包括活动方向和最大风速半径)在决定是否会产生足够的风暴潮作用, 进而留下粗粒沉积物方面也发挥了重要作用。在相似的古风暴气候背景值下, 随机过程可以部分解释古风暴活动的时空差异性。特定区域的某一时期的古风暴活动频率可能是只是随机过程而非气候因素控制(Wallace et al, 2021a)。
不同地点的沉积过程也存在显著差异, NY-4和AL-1的沉积物来源、水体深度、潟湖的封闭程度等局部水动力环境都可能造成事件记录的差异。此外, 相对较低的分辨率也增加了NY-4年代模型的不确定性。风暴识别和统计方法的差异也可能是造成两个记录差异的原因。NY-4基于经验原则, 以>1mm粗颗粒峰作为事件层识别标准, 这可能包含了>2mm生物碎屑或者珊瑚断枝导致的异常峰, 造成风暴层的识别偏差。NY-4记录中>1mm粗颗粒峰也并未完全捕捉到现代器测记录中的风暴事件, 表明这一识别方法中粗粒峰的风暴指示意义可能存在偏差, 风暴事件识别的阈值条件或划分标准需进一步校准。遗憾的是, 受远海采样条件的制约, AL-1和NY-4均只提供了单一钻孔的记录, 缺乏对照孔, 增加了利用该记录推断南海南部历史风暴活动的不确定性。因此, 更多来自同一研究区域的高分辨率风暴重建记录将有助于更好地约束风暴活动的随机变化特征, 更完整地还原区域风暴活动历史, 从而更好探讨南海南部安乐礁和永暑礁记录之间的时空差异及原因。

4 结论

通过对南沙群岛安乐环礁礁后潟湖粗粒沉积物的研究, 本文获得了小冰期以来的准年分辨率热带气旋活动记录。采样分辨率和测年精度的提高, 增进了我们对小冰期南海南部热带气旋活动的理解。研究表明潟湖沉积可靠地记录了经过安乐礁115km范围内的Cat 1级别和更强级别的热带气旋活动。在AD 1945—2020, AL-1岩心顶部沉积物记录了6个可追溯的风暴事件。小冰期以来, 南海南部安乐礁热带气旋活动在年代际到百年尺度上频繁变化, 发育两个主要的风暴活跃期。在小冰期早期(AD 1471—1620)经历了最为强烈的风暴活跃期, 另一个风暴活跃期位于现代暖期的AD 1930— 1960, 热带气旋活动虽有所加强, 但明显低于小冰期早期。南海南部区域重建记录的对比表明热带气旋活动具有较大的时空差异性, 这可能与局部沉积环境、热带气旋活动的随机性变化、识别方法差异等有关。来自同一研究区域的多条高分辨率古风暴重建记录将有助于更好地约束热带气旋活动的随机变化特征, 增加热带气旋活动重建的准确度。
[1]
高抒, 贾建军, 杨阳, 等, 2019. 陆架海岸台风沉积记录及信息提取[J]. 海洋学报, 41(10): 141-160.

GAO SHU, JIA JIANJUN, et al, 2019. Obtaining typhoon information from sedimentary records in coastal-shelf waters[J]. Haiyang Xuebao, 41(10): 141-160. (in Chinese with English abstract)

[2]
廖菲, 李文婷, 张子然, 等, 2019. 1949-2017年南海海域热带气旋强度和路径快速变化统计特征[J]. 海洋学报, 41(9): 126-135.

LIAO FEI, LI WENTING, ZHANG ZIRAN, et al, 2019. Analysis of rapid changes of tropical cyclones over the South China Sea for 1949-2017[J]. Haiyang Xuebao, 41(9): 126-135. (in Chinese with English abstract)

[3]
陶丽, 张艺帆, 王学兵, 2020. 南海与西北太平洋地区夏季热带气旋潜在生成指数的改进[J]. 大气科学学报, 43(4): 603-616.

TAO LI, ZHANG YIFAN, WANG XUEBING, 2020. Improvement of genesis potential index for western North Pacific tropical cyclones[J]. Transactions of Atmospheric Sciences, 43(4): 603-616. (in Chinese with English abstract)

[4]
徐笑梅, 高抒, 周亮, 等, 2019. 海南岛东北部海岸极端波浪事件沉积记录[J]. 海洋学报, 41(6): 48-63.

XU XIAOMEI, GAO SHU, ZHOU LIANG, et al, 2019. Sedimentary records of extreme wave events on the northeastern Hainan Island coast, southern China[J]. Haiyang Xuebao, 41(6): 48-63. (in Chinese with English abstract)

[5]
张敏, 罗军, 胡金磊, 等, 2019. 雷州市沿海风暴潮淹没危险性评估[J]. 热带海洋学报, 38(2): 1-12.

DOI

ZHANG MIN, LUO JUN, HU JINLEI, et al, 2019. Inundation risk assessment of storm surge along Lei Zhou coastal areas[J]. Journal of Tropical Oceanography, 38(2): 1-12. (in Chinese with English abstract)

DOI

[6]
周亮, 高抒, 杨阳, 等, 2015. 海南岛东南部海湾350年古风暴事件沉积与历史文献记录对比[J]. 海洋学报, 37(9): 84-94.

ZHOU LIANG, GAO SHU, et al, 2015. Comparison of paleostorm events between sedimentary and historical archives: A 350 year record from southeastern Hainan Island coastal embayments[J]. Haiyang Xuebao, 37(9): 84-94. (in Chinese with English abstract)

[7]
BACMEISTER J T, REED K A, HANNAY C, et al, 2018. Projected changes in tropical cyclone activity under future warming scenarios using a high-resolution climate model[J]. Climatic Change, 146(3-4): 547-560.

DOI

[8]
BLAAUW M, CHRISTEN J A, 2011. Flexible paleoclimate age-depth models using an autoregressive gamma process[J]. Bayesian Analysis, 6(3): 457-474.

DOI

[9]
BRAITHWAITE C J R, 1973. Settling behaviour related to sieve analysis of skeletal sands[J]. Sedimentology, 20(2): 251-262.

DOI

[10]
BRAMANTE J F, FORD M R, KENCH P S, et al, 2020. Increased typhoon activity in the Pacific deep tropics driven by Little Ice Age circulation changes[J]. Nature Geoscience, 13(12): 806-811.

DOI

[11]
CHEN H-F, LIU Y-C, CHIANG C-W, et al, 2019a. China's historical record when searching for tropical cyclones corresponding to Intertropical Convergence Zone (ITCZ) shifts over the past 2 kyr[J]. Climate of the Past, 15(1): 279-289.

DOI

[12]
CHEN TIANRAN, ROFF G, FENG YUEXING, et al, 2019b. Tropical sand cays as natural Paleocyclone archives[J]. Geophysical Research Letters, 46(16): 9796-9803.

DOI

[13]
CHENG H, EDWARDS R L, HOFF J, et al, 2000. The half-lives of uranium-234 and thorium-230[J]. Chemical Geology, 169(1-2): 17-33.

DOI

[14]
CLARK T R, ROFF G, ZHAO JIANXIN, et al, 2014. Testing the precision and accuracy of the U-Th chronometer for dating coral mortality events in the last 100 years[J]. Quaternary Geochronology, 23: 35-45.

DOI

[15]
CONTRERAS-ROSALES L A, JENNERJAHN T, STEINKE S, et al, 2019. Holocene changes in biome size and tropical cyclone activity around the Northern South China Sea[J]. Quaternary Science Reviews, 215: 45-63.

DOI

[16]
DIETRICH W E, 1982. Settling velocity of natural particles[J]. Water Resources Research, 18(6): 1615-1626.

DOI

[17]
DONNELLY J P, HAWKES A D, LANE P, et al, 2015. Climate forcing of unprecedented intense-hurricane activity in the last 2000 years[J]. Earth’s Future, 3(2): 49-65.

DOI

[18]
ISHIZAWA T, GOTO K, YOKOYAMA Y, et al, 2020. Dating tsunami deposits: Present knowledge and challenges[J]. Earth-Science Reviews, 200: 102971.

DOI

[19]
KENCH P S, MCLEAN R F, 1996. Hydraulic characteristics of bioclastic deposits: New possibilities for environmental interpretation using settling velocity fractions[J]. Sedimentology, 43(3): 561-570.

DOI

[20]
KENCH P S, BRANDER R W, 2006. Response of reef island shorelines to seasonal climate oscillations: South Maalhosmadulu atoll, Maldives[J]. Journal of Geophysical Research: Earth Surface, 111(F1): F01001.

[21]
KLOSTERMANN L, GISCHLER E, STORZ D, et al, 2014. Sedimentary record of late Holocene event beds in a mid-ocean atoll lagoon, Maldives, Indian Ocean: Potential for deposition by tsunamis[J]. Marine Geology, 348: 37-43.

DOI

[22]
KNUTSON T, CAMARGO S J, CHAN J C L, et al, 2020. Tropical cyclones and climate change assessment: Part II: Projected response to anthropogenic warming[J]. Bulletin of the American Meteorological Society, 101(3): E303-E322.

DOI

[23]
KNUTSON T R, SIRUTIS J J, ZHAO Ming, et al, 2015. Global projections of intense tropical cyclone activity for the late twenty-first century from dynamical downscaling of CMIP5/RCP4.5 scenarios[J]. Journal of Climate, 28(18): 7203-7224.

DOI

[24]
LANDSEA C W, HARPER B A, HOARAU K, et al, 2006. Can we detect trends in extreme tropical cyclones?[J]. Science, 313(5786): 452-454.

DOI

[25]
LANE P, DONNELLY J P, WOODRUFF J D, et al, 2011. A decadally-resolved paleohurricane record archived in the late Holocene sediments of a Florida sinkhole[J]. Marine Geology, 287(1-4): 14-30.

DOI

[26]
LIN NING, LANE P, EMANUEL K A, et al, 2014. Heightened hurricane surge risk in northwest Florida revealed from climatological-hydrodynamic modeling and paleorecord reconstruction[J]. Journal of Geophysical Research: Atmospheres, 119(14): 8606-8623.

DOI

[27]
LIU ENTAO, ZHAO JIANXIN, FENG YUEXING, et al, 2016. U-Th age distribution of coral fragments from multiple rubble ridges within the Frankland Islands, Great Barrier Reef: implications for past storminess history[J]. Quaternary Science Reviews, 143: 51-68.

DOI

[28]
LIU K B, SHEN CAIMING, LOUIE K S, 2001. A 1,000-year history of typhoon landfalls in Guangdong, southern China, reconstructed from Chinese historical documentary records[J]. Annals of the Association of American Geographers, 91(3): 453-464.

DOI

[29]
MAIKLEM W R, 1968. Some hydraulic properties of bioclastic carbonate grains[J]. Sedimentology, 10(2): 101-109.

DOI

[30]
MEI WEI, XIE SHANGPING, PRIMEAU F, et al, 2015. Northwestern pacific typhoon intensity controlled by changes in ocean temperatures[J]. Science Advances, 1(4): e1500014.

DOI

[31]
PEDUZZI P, CHATENOUX B, DAO H, et al, 2012. Global trends in tropical cyclone risk[J]. Nature Climate Change, 2(4): 289-294.

DOI

[32]
RODYSILL J R, DONNELLY J P, SULLIVAN R, et al, 2020. Historically unprecedented Northern Gulf of Mexico hurricane activity from 650 to 1250 CE[J]. Scientific Reports, 10(1): 19092.

DOI PMID

[33]
SORIA J L A, SWITZER A D, VILLANOY C L, et al, 2016. Repeat storm surge disasters of typhoon Haiyan and its 1897 predecessor in the Philippines[J]. Bulletin of the American Meteorological Society, 97(1): 31-48.

DOI

[34]
SUN LI, BI SHUOBEN, CHEN CHANGCHUN, et al, 2020. Typhoon frequency sequence reconstruction and characteristics analysis in the Southeast coastal area over China during the Ming and Qing dynasties[J]. Natural Hazards, 100(3): 1105-1116.

DOI

[35]
TOOMEY M R, DONNELLY J P, WOODRUFF J D, 2013. Reconstructing mid-late Holocene cyclone variability in the Central Pacific using sedimentary records from Tahaa, French Polynesia[J]. Quaternary Science Reviews, 77: 181-189.

DOI

[36]
WALLACE E J, DONNELLY J P, HENGSTUM P J, et al, 2019. Intense Hurricane activity over the past 1500 years at South Andros Island, the Bahamas[J]. Paleoceanography and Paleoclimatology, 34(11): 1761-1783.

DOI

[37]
WALLACE E J, COATS S, EMANUEL K, et al, 2021a. Centennial-scale shifts in storm frequency captured in paleohurricane records from the Bahamas arise predominantly from random variability[J]. Geophysical Research Letters, 48(1): e2020GL091145.

[38]
WALLACE E J, DONNELLY J P, VAN HENGSTUM P J, et al, 2021b. 1,050 years of hurricane strikes on long island in the Bahamas[J]. Paleoceanography and Paleoclimatology, 36(3): e2020PA004156.

[39]
WINKLER T S, VAN HENGSTUM P J, DONNELLY J P, et al, 2020. Revising evidence of hurricane strikes on Abaco Island (The Bahamas) over the last 700 years[J]. Scientific Reports, 10(1): 16556.

DOI

[40]
YANG HONGQIANG, YU KEFU, ZHAO MEIXIA, et al, 2015. Impact on the coral reefs at Yongle Atoll, Xisha Islands, South China Sea from a strong typhoon direct sweep: Wutip, September 2013[J]. Journal of Asian Earth Sciences, 114: 457-466.

DOI

[41]
YU KEFU, ZHAO JIANXIN, SHI QI, et al, 2009. Reconstruction of storm/tsunami records over the last 4000 years using transported coral blocks and lagoon sediments in the southern South China Sea[J]. Quaternary International, 195(1-2): 128-137.

DOI

[42]
YUE YUANFU, YU KEFU, TAO SHICHEN, et al, 2019. 3500-year western Pacific storm record warns of additional storm activity in a warming warm pool[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 521: 57-71.

DOI

[43]
ZUKI Z M, LUPO A R, 2008. Interannual variability of tropical cyclone activity in the southern South China Sea[J]. Journal of Geophysical Research: Atmosphere, 113(D6): D06106.

Outlines

/