基于ERA5再分析數(shù)據(jù)的中國鄰近海域極端波高特征分析*

杜文彥 張旭日 張麗麗 尤再進(jìn) 石洪源

基于ERA5再分析數(shù)據(jù)的中國鄰近海域極端波高特征分析*

杜文彥1張旭日1張麗麗1尤再進(jìn)2石洪源1①

(1. 魯東大學(xué)水利工程學(xué)院 山東煙臺(tái) 264000; 2. 大連海事大學(xué)港口與航運(yùn)安全協(xié)同創(chuàng)新中心 遼寧大連 116000)

極端波浪對(duì)沿海地區(qū)基礎(chǔ)設(shè)施有著深遠(yuǎn)的影響, 了解它們的變化規(guī)律是進(jìn)行海岸帶風(fēng)險(xiǎn)分析和災(zāi)害預(yù)防的基礎(chǔ)。文章基于歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasting, ECMWF)第五代再分析數(shù)據(jù)集(ECMWF reanalysis v5, ERA5),對(duì)中國鄰近海域1979～2018年間極端波高展開時(shí)空特征分析, 并統(tǒng)計(jì)了40 a厄爾尼諾系數(shù), 利用廣義極值分布(generalized extreme value, GEV)函數(shù)探究了近40 a厄爾尼諾現(xiàn)象對(duì)我國海域極端波高的影響, 結(jié)果顯示: 統(tǒng)計(jì)1979～2018年整個(gè)研究區(qū)域前2%極端波浪年均值在6～10 m間浮動(dòng), 且整體趨勢(shì)遞增, 在四季趨勢(shì)變化中, 春夏極端波高增長(zhǎng)趨勢(shì)較秋冬高, 且波動(dòng)明顯, 在年際極端波高變化趨勢(shì)中有較大波動(dòng)時(shí)大多伴隨著厄爾尼諾或是拉尼娜現(xiàn)象的發(fā)生, 結(jié)合厄爾尼諾現(xiàn)象對(duì)GEV分布中位置參數(shù)的影響分布圖和極端波高年、季節(jié)際趨勢(shì)變化分布圖, 厄爾尼諾現(xiàn)象影響大的地區(qū)極端波高大多呈現(xiàn)增長(zhǎng)趨勢(shì), 表明厄爾尼諾現(xiàn)象對(duì)極端波高有較高影響。

極端波高; 厄爾尼諾; 中國海; ERA5 (ECMWF第五代再分析數(shù)據(jù)集)

海浪會(huì)對(duì)海岸結(jié)構(gòu)物、近海泥沙運(yùn)輸、海岸侵蝕等造成影響。因此, 它是沿海災(zāi)害的一個(gè)重要因素。大浪疊加潮汐將進(jìn)一步增加沿海地區(qū)的脆弱性(Wang, 2012)。海浪由風(fēng)產(chǎn)生, 很多學(xué)者已經(jīng)研究了全球波高的歷史變化趨勢(shì)(Patra, 2016; Patra, 2017)。此外, 年際氣候變化模式如厄爾尼諾-南方濤動(dòng)(El Ni?o-Southern Oscillation, ENSO)、北大西洋濤動(dòng)(North Atlantic Oscillation, NAO)等將造成大氣環(huán)流的改變, 最終影響整個(gè)全球海洋表面的風(fēng)應(yīng)力(Dodet, 2010)。因此, 氣候變化通過影響海洋-大氣相互作用、風(fēng)和臺(tái)風(fēng)活動(dòng)來影響全球波浪的變化(Hemer, 2010; Bromirski, 2013)。

極端波浪對(duì)沿海地區(qū)基礎(chǔ)設(shè)施有著深遠(yuǎn)的影響, 了解它們的變化規(guī)律是進(jìn)行海岸帶風(fēng)險(xiǎn)分析和災(zāi)害預(yù)防的基礎(chǔ)。然而, 關(guān)于氣候變化對(duì)全球范圍內(nèi)極端波浪影響的研究十分有限。眾所周知, 氣候變化對(duì)平均波高具有影響, 但由于短時(shí)間尺度的相關(guān)分析噪聲較大, 它們對(duì)極端波浪的影響尚未被深刻認(rèn)識(shí)。此外, 由于極端波浪的非正態(tài)性, 簡(jiǎn)單線性回歸分析不能應(yīng)用于其研究。因此, 學(xué)者們開始使用極端函數(shù)理論研究氣候變化與極端波浪的相關(guān)性。Izaguirre等(2011)對(duì)衛(wèi)星高度計(jì)數(shù)據(jù)集的月最大顯著波高(significant wave heights, SWH)進(jìn)行了非平穩(wěn)廣義極值(generalized extreme value, GEV)分析, 將氣候變化指數(shù)作為GEV位置參數(shù)的附加協(xié)變量。分析了ENSO、NAO等10種大尺度氣候變化對(duì)全球年極端波高的影響, 但尚未評(píng)估對(duì)極端波高的季節(jié)貢獻(xiàn)。Kumar等(2016)基于歐洲中期天氣預(yù)報(bào)中心再分析數(shù)據(jù)集ERA-Interim和ERA-20C再分析數(shù)據(jù), 利用GEV分析研究了ENSO和PDO對(duì)北方冬季(11月到次年2月)期間的極端波高的影響。Kumar等(2019)利用相同的方法及ERA-20C再分析數(shù)據(jù), 研究了印度洋上的極端SWH對(duì)ENSO、SAM的季節(jié)性反應(yīng), 發(fā)現(xiàn)了同樣的物理機(jī)制。Patra等(2020)等利用同樣的方法及衛(wèi)星和歐洲中期天氣預(yù)報(bào)中心第五代再分析數(shù)據(jù)集(ECMWF reanalysis v5, ERA5), 分析了ENSO等10余種氣候變化對(duì)全球極端波高的影響。

目前, 針對(duì)中國近海波浪特征的研究很多, 有的基于實(shí)測(cè)數(shù)據(jù)(吳迪茜等, 2021), 有的基于衛(wèi)星遙感資料(李文博等, 2017)和數(shù)值模型(鄭崇偉等, 2013), 也有的基于再分析數(shù)據(jù)(喬守文等, 2020; 邱文博等, 2021), 研究結(jié)果指出, 我國海域波浪呈現(xiàn)出明顯的季節(jié)特征, 其中秋冬季節(jié)波高大, 春夏季節(jié)波高小。因全球變暖等因素, 厄爾尼諾現(xiàn)象發(fā)生頻率及影響程度不斷增強(qiáng), 目前有關(guān)我國海域極端波高的變化特征研究較少, 更鮮有研究定量指出厄爾尼諾現(xiàn)象與極端波高之間的關(guān)聯(lián)性。因此, 本文采用ERA5數(shù)據(jù)分析我國極端波高的變化特征, 并利用GEV函數(shù)探究厄爾尼諾現(xiàn)象對(duì)我國海域極端波高的影響。

1 數(shù)據(jù)資料

1.1 研究區(qū)域

中國近海海域包括渤海、黃海、東海和南海以及臺(tái)灣東側(cè)的太平洋部分海域, 位于北太平洋西部的邊緣地區(qū)南北跨越44個(gè)緯度, 本文研究的中國臨近海范圍為(7°～41°N, 105.5°～126.5°E)(圖1)。

1.2 數(shù)據(jù)來源

ECMWF是全球幾家最主要的再分析數(shù)據(jù)中心之一。ERA5是ECMWF第五代再分析數(shù)據(jù)集。它是歐洲中期天氣預(yù)報(bào)中心最新的再分析產(chǎn)品。ERA5涵蓋了1979年到最近的再分析數(shù)據(jù), 目前還在不斷更新中。相比ERA-Interim, ERA5同化系統(tǒng)使用了若干專為再分析而開發(fā)的綜合預(yù)報(bào)系統(tǒng)(integrated forecasting system, IFS cycle 41r2), 數(shù)據(jù)準(zhǔn)確性得到了顯著改善。ERA5數(shù)據(jù)具有高分辨率: 水平分辨率為31 km (～0.25°), 垂向?yàn)閺牡乇黹_始137層。以往ERA數(shù)據(jù)產(chǎn)品在波浪波高模擬中普遍存在偏低的現(xiàn)象, 而ERA5數(shù)據(jù)已經(jīng)有較大的改進(jìn)(李朝等, 2021), ERA5提供的變量從ERA-Interim的100個(gè)增加到240個(gè), 其中包括耦合波浪模型提供的波高和波向, 使用戶能夠更準(zhǔn)確地分析過去的大氣和海洋狀態(tài)。許多研究人員將其用于大氣和海洋研究(Stopa, 2013)。ERA5在我國的適用性已有學(xué)者(Shi, 2021)研究, 因此, 本文選擇ERA5再分析數(shù)據(jù)用于本文相關(guān)研究。

圖1 研究區(qū)域地形水深圖

本文按照世界氣象學(xué)會(huì)的定義, 四季劃分標(biāo)準(zhǔn)為: 春季(MAM), 包括3、4和5月; 夏季(JJA), 包括6、7和8月; 秋季(SON), 包括9、10和11月; 冬季(DJF), 包括12、1和2月。

2 研究區(qū)域極端波高時(shí)空變化特征分析

2.1 極端波高的年特征分析

早期的文章對(duì)極端波高的定義是, 在總體波高數(shù)據(jù)中選取一個(gè)最大值數(shù)據(jù), 但由于臺(tái)風(fēng)等極端天氣對(duì)這種極值的影響極大, 故這種方法可能會(huì)提供不可靠的極端樣本值的估計(jì)。為克服這一缺點(diǎn), 本文采用了波高大于整體數(shù)據(jù)閾值的波高均值來定義極端波高, 所選閾值為每個(gè)數(shù)據(jù)點(diǎn)的所有波高數(shù)據(jù)的第98個(gè)百分位(Menéndez, 2008)。

2.1.1 極端波高年趨勢(shì)變化空間分布 通過對(duì)1979～2018年前2%極端波高統(tǒng)計(jì)和推算, 得出極端波高年均圖、波高年趨勢(shì)變化分布圖和50年一遇及100年一遇(圖2)。由圖2可見, 渤海及黃海海域的極端波高年均值約在3 m, 東海南部、臺(tái)灣海峽周邊及南海東北部海域極端波高年均值較大, 約在5～6 m。

由圖2可見, 臺(tái)灣島東北部地區(qū)及東海南部大部分地區(qū)極端波高呈明顯增長(zhǎng)趨勢(shì), 增長(zhǎng)幅度在0.025 m/a左右, 南海大部分海域也處于增長(zhǎng)趨勢(shì), 增長(zhǎng)幅度基本在0.01 m/a左右浮動(dòng), 渤海、黃海大部分海域及海南島周邊海域呈減小趨勢(shì), 其削減幅度在0.005 ～0.010 m/a左右。在50年一遇和100年一遇極值波高圖中, 渤海和黃海大部分海域的波高在4～6 m, 可以看出波高大值區(qū)域集中在東海海域, 南海大部分海域波高在10 m左右, 在東沙群島和巴士海峽附近波高較大, 在11～14 m左右。

圖2 1979～2018年極端波高年均值圖(a)、極端波高年趨勢(shì)變化空間分布圖(b)、50年一遇極值波高推算圖(c)和100年一遇極值波高推算圖(d)

2.1.2 極端波高年際變化分布 圖3為1979～2018年間本文研究區(qū)域2%極端年均波浪高度的變化情況, 整個(gè)區(qū)域逐小時(shí)的有效波高累加求平均, 可以看出本文研究海域2%極端年均波浪在3～4 m間浮動(dòng), 且整體趨勢(shì)遞增, 其中1988年和1996年波高值出現(xiàn)較大波動(dòng)。

圖3 1979～2018年極端波高季節(jié)變化趨勢(shì)圖

2.2 極端波高的季節(jié)特征分析

2.2.1 極端波高季節(jié)趨勢(shì)變化空間分布 通過對(duì)40年季均前2%極端波高統(tǒng)計(jì), 得出極端波高季均分布圖及波高季節(jié)趨勢(shì)變化分布圖(圖4～5)。由圖4可見, 春季和夏季的季節(jié)均值極端波高較秋季和冬季小, 春季大部分海域極端波高在3.5～4 m左右, 臺(tái)灣海峽西側(cè)極端波高較高在4.5 m左右; 夏季東海東南部極端波高值在5 m以上在四個(gè)海域中最為顯著; 秋季研究海域較大值集中在東海南部和南海中北部; 冬季整個(gè)研究區(qū)域極端波高值普遍較高, 大部分海域極端波高在4～4.5 m左右。

1979～2018年極端波高各季節(jié)趨勢(shì)變化分布圖顯示, 春季渤海及黃海呈下降, 其下降幅度約在0.005 m/a, 太平洋西部地區(qū)則為增長(zhǎng)趨勢(shì), 增長(zhǎng)幅度在0.035 m/a左右, 南海大部分地區(qū)增長(zhǎng)幅度在0.02 m/a; 在夏季和秋季, 南海西南部為下降趨勢(shì), 在0.02 m/a左右, 而東海南部及西部則為增長(zhǎng)趨勢(shì), 夏季在臺(tái)灣島東部增長(zhǎng)趨勢(shì)明顯, 約在0.04 m/a, 秋季臺(tái)灣島東北部及東海中南部增長(zhǎng)趨勢(shì)幅度約為0.03 m/a; 冬季期間, 渤海地區(qū)呈現(xiàn)下降趨勢(shì), 幅度在0.01 m/a左右浮動(dòng), 臺(tái)灣海峽南側(cè)及西沙群島周邊海域則為明顯的上升趨勢(shì), 上升幅度在0.025 m/a左右。

圖4 1979～2018年極端波高四季均值分布圖

圖5 1979～2018年極端波高四季趨勢(shì)變化空間分布

2.2.2 極端波高季節(jié)變化分布 為了體現(xiàn)變化趨勢(shì)的季節(jié)性差異, 本文還分析了近40 a期間研究區(qū)域海域逐春季、逐夏季、逐秋季、逐冬季2%極端波浪的變化趨勢(shì)(圖6)。由圖6可見, 大部分海域的大浪頻率表現(xiàn)出顯著性遞增趨勢(shì); 春季波高波動(dòng)幅度最大, 幅度最顯著時(shí)能達(dá)到1 m, 在1.9～3.5 m之間且整體波高較小, 其中在1991、1996、2005、2006和2011年波動(dòng)較為明顯; 夏季波浪高度波動(dòng)明顯, 極端波高在1.5～3.9 m左右, 其中1987、1988、1998、2010和2015年波動(dòng)值顯著, 最大可達(dá)到2 m左右; 夏季和秋季極端波高變化趨勢(shì)較春季和冬季相對(duì)平緩, 秋季極端波高較大在3～4.5 m之間浮動(dòng), 其中在1983、1993、2007、2009和2015年出現(xiàn)突變值, 最大波動(dòng)在1983年出現(xiàn)可達(dá)到1.4 m左右; 冬季極端波高變化趨勢(shì)為上升趨勢(shì), 波高在3～4.2 m左右浮動(dòng), 在1981、1988、1991、1993和1999年出現(xiàn)明顯波動(dòng), 在1993年發(fā)生波動(dòng)最顯著, 達(dá)到0.8 m。

3 厄爾尼諾現(xiàn)象對(duì)研究區(qū)域極端波高的影響

3.1 厄爾尼諾系數(shù)(Ni?o3.4)統(tǒng)計(jì)分析

本文使用Ni?o3.4區(qū)域平均的海表面溫度異常作為厄爾尼諾指數(shù)。根據(jù)中國氣象局頒布的《厄爾尼諾/拉尼娜事件判定方法》, 定義Ni?o3.4指數(shù)≥0.5 °C并持續(xù)5個(gè)月以上為一次厄爾尼諾事件, Ni?o3.4指數(shù)≤-0.5 °C并持續(xù)5個(gè)月以上為一次拉尼娜事件。另外, 定義事件峰值強(qiáng)度絕對(duì)值達(dá)到或超過2.5 °C的事件為超強(qiáng)事件。為研究ENSO的影響, 本文利用美國國家海洋和大氣管理局(National Oceanic and Atmospheric Administration, NOAA)-氣候預(yù)測(cè)中心(Climate Prediction Center, CPC)(https://www.cpc.ncep. noaa.gov/data/indices/ersst5.nino.mth.81-10.ascii)每月的Ni?o3.4指數(shù)(海表面溫度在5°S～5°N, 170°～120°W上的平均異常, 基于ERST.v5計(jì)算)進(jìn)行研究分析。從圖7可看出1979～2018年的厄爾尼諾指數(shù)。根據(jù)厄爾尼諾事件的定義, 1979～2018年共發(fā)生了10次厄爾尼諾事件和10次拉尼娜事件, 其中1982/1983、1986/1987、1991/1992、1994/1995、1997/1998、2002/2003、2004/2005、2006/2007、2009/2010、2015/2016為厄爾尼諾事件; 1983/1984、1984/1985、1988/1989、1995/1996、1998/2000、2005/2006、2007/2008、2010/2011、2011/2012、2017/2018為拉尼娜事件。其中, 2015/2016年厄爾尼諾事件峰值強(qiáng)度為2.8 °C, 是一次超強(qiáng)厄爾尼諾事件。

圖6 1979～2018年極端波高季節(jié)變化趨勢(shì)圖

圖7 1979～2018年厄爾尼諾系數(shù)(Ni?o3.4)統(tǒng)計(jì)

3.2 極端波高特征與厄爾尼諾系數(shù)相關(guān)性分析

3.2.1 研究方法 極值理論(extreme value theory, EVT)是統(tǒng)計(jì)學(xué)的重要分支之一, 在氣象和水文極端中得到了廣泛的應(yīng)用(Shinyie, 2013; Vasiliades, 2015)。極值統(tǒng)計(jì)中有兩個(gè)最基本的分布, 廣義極值分布(generalized extreme value, GEV)和廣義帕累托分布(generalized Pareto distribution, GPD)。GEV分布用統(tǒng)一的形式來表示Gumbel分布、Frechet分布和Weibull分布三種極值分布類型, 通常選取數(shù)據(jù)中的“年最大值”或“區(qū)組最大值”數(shù)據(jù)來建立模型。GEV累積分布函數(shù)表示為

3.2.2 研究區(qū)域極端波高與厄爾尼諾系數(shù)相關(guān)性分析

圖8 1979～2018年相關(guān)性年空間分布

Fig.8 Spatial distribution of correlation years from 1979 to 2018

3.2.2.2 相關(guān)性季節(jié)空間分布 Ni?o3.4極端波響應(yīng)模式的季節(jié)變化如圖9所示。GEV分析使用位置參數(shù)和Ni?o3.4的回歸系數(shù)季節(jié)平均指數(shù)。顯著的厄爾尼諾事件影響全年持續(xù), 并且在不同季節(jié)的空間范圍有所不同。春季中國渤海、黃海以及太平洋西部地區(qū)受氣候指數(shù)影響大; 在夏季, 氣候指數(shù)對(duì)南海東南部、東海東部及南部影響較大; 在秋季, 南海北部地區(qū)可以看到強(qiáng)烈的氣候指數(shù)影響; 冬季期間, 黃海東南部、東海部分海域、渤海地區(qū)、臺(tái)灣海峽西側(cè)及西沙群島周邊海域受氣候指數(shù)影響較大。結(jié)合2.2.1中極端波高季節(jié)趨勢(shì)變化分布, 可以看出在厄爾尼諾現(xiàn)象影響大的地區(qū)極端波高普遍具有較大幅度波動(dòng)。

圖9 1979～2018年相關(guān)性四季空間分布

4 結(jié)論

基于研究海域ERA5后報(bào)資料, 本文詳細(xì)分析了2%極端波浪的空間分布情況、長(zhǎng)期年均變化趨勢(shì)及變化趨勢(shì)的空間分布情況, 極端波浪季節(jié)分布、趨勢(shì)變化空間分布以及極端波浪的季均變化趨勢(shì)和厄爾尼諾系數(shù)變化, 并通過GEV極值波高與厄爾尼諾系數(shù)回歸方程結(jié)合分析得出厄爾尼諾現(xiàn)象對(duì)極端波高的影響分布, 結(jié)果表明:

(1) 由2%極端波高的數(shù)值年際變化可以看出, 在長(zhǎng)期變化趨勢(shì)中, 本文研究海域2%極端波浪年均值在3～4 m間浮動(dòng), 且整體趨勢(shì)遞增, 其中1988和1996年波高值出現(xiàn)的波動(dòng)最為顯著。渤海及黃海海域的極端波高年均值約在3 m, 東海南部、臺(tái)灣海峽周邊及南海東北部海域極端波高年均值較大, 約在5～6 m。

在四季趨勢(shì)變化中, 春夏極端波高增長(zhǎng)趨勢(shì)較秋冬高, 且波動(dòng)明顯。春季在1.9～3.5 m之間且整體波高較小, 其中在1991、1996、2005、2006和2011年波動(dòng)明顯, 幅度最顯著時(shí)能達(dá)到1 m; 夏季波浪高度波動(dòng)明顯, 極端波高在1.5～3.9 m左右, 其中1987、1988、1998、2010和2015年波動(dòng)值顯著, 最大可達(dá)到2 m左右。春季和夏季的季節(jié)均值極端波高較秋季和冬季小, 春季大部分海域極端波高在3.5～4 m左右; 夏季和秋季極端波高變化趨勢(shì)較春季和冬季相對(duì)平緩; 秋季研究海域較大值集中在東海南部和南海中北部; 冬季整個(gè)研究區(qū)域極端波高值普遍較高, 大部分海域極端波高在4～4.5 m左右。

結(jié)合厄爾尼諾系數(shù)變化趨勢(shì)和極端波高的數(shù)值年際變化, 極端波高出現(xiàn)顯著波動(dòng)值的1988和1995年可以看到, 1994和1995年都發(fā)生了厄爾尼諾和拉尼娜現(xiàn)象, 且1987年整年都為厄爾尼諾現(xiàn)象; 而在四季極端波高有較大波動(dòng)時(shí)大多伴隨著厄爾尼諾或是拉尼娜現(xiàn)象的發(fā)生。

(2) 在2%極端波高年際趨勢(shì)變化分布中, 南海與東海為增長(zhǎng)趨勢(shì), 渤海海域及黃海部分海域呈下降趨勢(shì)。在極端波高季節(jié)際趨勢(shì)變化分布中, 春季和冬季渤海及黃海呈下降趨勢(shì), 太平洋西部地區(qū)和南海大部分地區(qū)則為增長(zhǎng)趨勢(shì); 在夏季, 南海西南部為下降趨勢(shì), 而臺(tái)灣島東部為明顯增長(zhǎng); 秋季臺(tái)灣島東北部及東海中南部增長(zhǎng)趨勢(shì), 南海西南部為下降趨勢(shì)。

(3) 在極端波高的GEV位置參數(shù)對(duì)全年Ni?o3.4指數(shù)回歸系數(shù)的空間分布中, 東海大部分海域及南海海域受氣候變量影響較大; 在不同季節(jié)回歸系數(shù)分布中, 夏秋冬季氣候變量對(duì)南海的影響較大, 冬春對(duì)東海海域影響較大。結(jié)合厄爾尼諾現(xiàn)象對(duì)GEV分布中位置參數(shù)的影響分布圖和極端波高年、季節(jié)際趨勢(shì)變化分布圖, 可以看出在厄爾尼諾現(xiàn)象影響大的地區(qū)極端波高普遍具有較大幅度波動(dòng)。

喬守文, 孫嘉誠, 石洪源, 等, 2020. 基于ERA5的黃渤海附近海域波浪能資源時(shí)空特征分析[J]. 海洋與湖沼, 51(6): 1350-1358.

李文博, 張坤蘭, 吳克儉, 2017. 基于衛(wèi)星資料的南海波候特征研究[J]. 海洋湖沼通報(bào)(6): 11-20, doi: 10.13984/j.cnki. cn37-1141.2017.06.002.

李朝, 侯一筠, 李水清, 等, 2021. 兩類典型臺(tái)風(fēng)路徑影響下的黃、渤海海浪場(chǎng)特征研究[J]. 海洋與湖沼, 52(1): 51-65.

吳迪茜, 陳智杰, 2021. 平潭島東部海域波浪特征及波浪能分析[J]. 應(yīng)用海洋學(xué)學(xué)報(bào), 40(2): 293-302.

邱文博, 李冠宇, 徐俊臣, 等, 2021. 黃渤海海域波浪時(shí)空變化特征分析[J]. 海洋科學(xué), 45(7): 1-8.

鄭崇偉, 林剛, 邵龍?zhí)? 2013. 1988—2010年中國海大浪頻率及其長(zhǎng)期變化趨勢(shì)[J]. 廈門大學(xué)學(xué)報(bào)(自然科學(xué)版), 52(3): 395-399.

BROMIRSKI P D, CAYAN D R, HELLY J,, 2013. Wave power variability and trends across the North Pacific [J]. Journal of Geophysical Research: Oceans, 118(12): 6329- 6348, doi: 10.1002/2013JC009189.

DODET G, BERTIN X, TABORDA R, 2010. Wave climate variability in the North-East Atlantic Ocean over the last six decades [J]. Ocean Modelling, 31(3/4): 120-131, doi: 10. 1016/j.ocemod.2009.10.010.

HEMER M A, CHURCH J A, HUNTER J R, 2010. Variability and trends in the directional wave climate of the Southern Hemisphere [J]. International Journal of Climatology, 30(4): 475-491, doi: 10.1002/joc.1900.

IZAGUIRRE C, MéNDEZ F J, MENéNDEZ M,, 2011. Global extreme wave height variability based on satellite data [J]. Geophysical Research Letters, 38(10): L10607, doi: 10.1029/2011GL047302.

KUMAR P, KAUR S, WELLER E,, 2019. Influence of natural climate variability on the extreme ocean surface wave heights over the Indian Ocean [J]. Journal of Geophysical Research: Oceans, 124(8): 6176-6199, doi: 10.1029/2019JC015391.

KUMAR P, MIN S K, WELLER E,, 2016. Influence of climate variability on extreme ocean surface wave heights assessed from ERA-interim and ERA-20C [J]. Journal of Climate, 29(11): 4031-4046, doi: 10.1175/JCLI-D-15-0580.1.

MENéNDEZ M, MéNDEZ F J, LOSADA I J,, 2008. Variability of extreme wave heights in the northeast Pacific Ocean based on buoy measurements [J]. Geophysical Research Letters, 35(22): L22607, doi: 10.1029/2008GL035394.

PATRA A, BHASKARAN P K, 2016. Trends in wind-wave climate over the head Bay of Bengal region [J]. International Journal of Climatology, 36(13): 4222-4240, doi: 10.1002/ joc.4627.

PATRA A, BHASKARAN P K, 2017. Temporal variability in wind–wave climate and its validation with ESSO-NIOT wave atlas for the head Bay of Bengal [J]. Climate Dynamics, 49(4): 1271-1288, doi: 10.1007/s00382-016-3385-z.

PATRA A, MIN S K, SEONG M G, 2020. Climate variability impacts on global extreme wave heights: seasonal assessment using satellite data and ERA5 reanalysis [J]. Journal of Geophysical Research: Oceans, 125(12): e2020JC016754, doi: 10.1029/2020JC016754.

SHI H Y, CAO X F, LI Q J,, 2021. Evaluating the accuracy of ERA5 wave reanalysis in the water around China [J]. Journal of Ocean University of China, 20(1): 1-9.

SHINYIE W L, ISMAIL N, JEMAIN A A, 2013. Semi-parametric estimation for selecting optimal threshold of extreme rainfall events [J]. Water Resources Management, 27(7): 2325-2352.

STOPA J E, CHEUNG K F, TOLMAN H L,, 2013. Patterns and cycles in the Climate Forecast System Reanalysis wind and wave data [J]. Ocean Modelling, 70: 207-220, doi: 10.1016/j.ocemod.2012.10.005.

VASILIADES L, GALIATSATOU P, LOUKAS A, 2015. Nonstationary frequency analysis of annual maximum rainfall using climate covariates [J]. Water Resources Management, 29(2): 339-358.

WANG X L, FENG Y, SWAIL V R, 2012. North Atlantic wave height trends as reconstructed from the 20th century reanalysis [J]. Geophysical Research Letters, 39(18): L18705, doi: 10.1029/2012GL053381.

CHARACTERISTICS OF EXTREME WAVE HEIGHT IN CHINA’S MARGINAL SEAS BASED ON ERA5 REANALYSIS DATA

DU Wen-Yan1, ZHANG Xu-Ri1, ZHANG Li-Li1, YOU Zai-Jin2, SHI Hong-Yuan1

(1. School of Hydraulic Engineering, Ludong University, Yantai 264000, China; 2. Centre for Ports and Maritime Safety, Dalian Maritime University, Dalian 116000, China)

Extreme waves have a profound impact on coastal infrastructure, and it is the basis of understanding of their variation law and prevention of disasters in coastal zones. Based on the ECMWF reanalysis the fifth generation ERA5 data set, the spatio-temporal characteristics of extreme wave heights in China’s marginal seas from 1979 to 2018 were analyzed and the 40-year El Ni?o coefficients calculated. The GEV (generalized extreme value) function was used to explore the influence of El Ni?o event on extreme wave height in the study area in recent 40 years. Results show that the annual mean value of extreme waves in the top 2% of the whole research area fluctuated between 6 m and 10 m from 1979 to 2018 and the fluctuation is increasing in overall, which is more obvious in spring and summer for having greater extreme wave height than in autumn and winter duo reasonably to the El Ni?o or La Nina works. Combining the influence of El Ni?o on the location of GEV distribution and the interannual and inter-seasonal trends of extreme wave heights, we found that most of the extreme wave heights in areas where El Ni?o has a great influence show an increasing trend, which indicates that El Ni?o has a greater impact on extreme wave heights.

extreme wave height; El Ni?o; China Sea; ERA5 (ECMWF reanalysis v5)

* 基于我國資源特性的海洋能高效利用創(chuàng)新技術(shù)研發(fā), 2018YFB1501901號(hào); 山東省濱海沙灘防護(hù)工程環(huán)境災(zāi)害及防災(zāi)減災(zāi)對(duì)策研究, U1806227號(hào); 山東省粉沙質(zhì)海岸港口工程驟沖驟於災(zāi)害及防災(zāi)減災(zāi)措施研究(課題1), U1906231號(hào)。杜文彥, 碩士研究生, E-mail: 910718958@qq.com

石洪源, 碩士生導(dǎo)師, 講師, E-mail: shihongyuan1234@163.com

2022-01-03,

2022-03-10

P576

10.11693/hyhz20220100003