汪楊駿，侯太平，張 韌,2，錢龍霞，王 鋒，龍 強

(1. 解放軍理工大學(xué) 氣象海洋學(xué)院，江蘇 南京 211101；2. 南京信息工程大學(xué) 氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心，江蘇 南京 210044；3. 河北省唐山市曹妃甸工業(yè)區(qū)氣象局，河北 唐山 063000)

基于動態(tài)極值理論和Copula函數(shù)的極端海平面高度預(yù)測建模

汪楊駿1，侯太平1，張 韌1,2，錢龍霞1，王 鋒3，龍 強3

(1. 解放軍理工大學(xué) 氣象海洋學(xué)院，江蘇 南京 211101；2. 南京信息工程大學(xué) 氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心，江蘇 南京 210044；3. 河北省唐山市曹妃甸工業(yè)區(qū)氣象局，河北 唐山 063000)

全球氣候變化背景下，海平面上升是一個潛在的重大風(fēng)險，為防范氣候災(zāi)害，應(yīng)對極端氣象海洋事件，需客觀、定量地對未來極端海平面變化進(jìn)行科學(xué)預(yù)測。為此，基于Copula函數(shù)和動態(tài)極值分析理論，綜合考慮平均海平面變化(包括垂直陸地運動和基準(zhǔn)的局地變化)與潮、涌、浪等其他氣候變化的增水對極端海平面高度的影響，采用DREAM方法改進(jìn)Bayes推斷對動態(tài)極值模型的參數(shù)空間估計問題，提出一種新的模型對未來極端海平面高度變化進(jìn)行預(yù)測，旨在改進(jìn)傳統(tǒng)模型存在的不確定性問題，并運用該模型對氣候變化背景下廈門地區(qū)未來35年的海平面變化情景進(jìn)行了模型應(yīng)用和實驗?zāi)M。

極端海平面高度；動態(tài)極值模型；DREAM方法；Bayes推斷；Copula函數(shù)

Abstract: In the context of global climate change, sea level rise is a potentially major risk. In order to prevent the climate disasters,and to deal with extreme weather events, it is necessary to objectively and quantitatively predict the future extreme sea level changes. To this end, based on the theory of Copula function and that the dynamic extreme value analysis, the authers of this paper have made comprehensive consideration of the influence of the mean sea level changes (including the vertical land movement and the baseline) and the changes of the water level of the tide, surge and wave,and used Bayes method to improve the parameter estimation of DREAM. A new method is proposed to predict the future extreme sea level changes, and the proposed method is used to improve the existing problem of the traditional method, and a model is used to simulate the sea level changes in Xiamen area in 35 years.

Keywords: extreme sea level； dynamic extreme value model； DREAM； Bayes inference； Copula function

極端海平面事件是指由惡劣天氣氣候(如大風(fēng))引起的沿岸洪水、風(fēng)暴潮、高水位事件等。不同于平均海平面高度緩慢持續(xù)的變化，極端海平面高度的變化對沿海城市和港口有著直接而顯著的影響。據(jù)國家海洋局公布的《海洋災(zāi)害公報》顯示，僅2014年我國沿海共發(fā)生風(fēng)暴潮過程9次，造成直接經(jīng)濟損失135.78億元，死亡6人(含失蹤)。隨著我國“一帶一路”計劃的展開，作為“海上絲綢之路”重要聯(lián)接節(jié)點的沿海港口城市將迎來巨大的發(fā)展機遇，但與此同時，更為聚集的人口、密集的基礎(chǔ)設(shè)施，使得這些港口城市在遭受極端事件時所面臨的災(zāi)害也更加嚴(yán)重。因此，科學(xué)、準(zhǔn)確地預(yù)測未來極端海平面高度變化成為十分迫切的需求。

研究表明極端海平面高度變化的趨勢與平均海平面高度變化趨勢基本一致[1-2]。也有學(xué)者指出極端海平面高度變化大于平均海平面高度[3-4]，并認(rèn)為這可能是由于局地的氣候變化引起的[5]。Menendez等[6]指出ENSO對太平洋地區(qū)和季風(fēng)區(qū)的極端海平面高度的年際變化有重要的影響；在歐洲南部，NAO對極端海平面高度變化也有顯著負(fù)反饋作用[1]；Ullmann等[4]指出在卡馬圭，受由歐洲中部海平面氣壓上升引起的偏南風(fēng)影響，年最大海平面上升速率是平均海平面變化的兩倍；Abeysirigunawardena等[5]指出自20世紀(jì)70年代中葉以來，受強PDO正位相影響，年海平面高度極值變化趨勢是年平均海平面的兩倍。除上述氣候尺度的海平面變化特征外，多數(shù)學(xué)者認(rèn)為風(fēng)暴潮增水對極端海平面高度有著更為顯著地影響[7-9]。

極端海平面高度變化有許多預(yù)測方法，如Hunter[10]采用耿貝爾極值分布模型和全球海平面的模式預(yù)估數(shù)據(jù)，分析了2100年全球平均極端海平面高度期望及其頻率變化，由于全球尺度上，無需考慮局地氣候因素影響，因此風(fēng)暴潮的影響也可被假定為固定不變。Tebaldi等[8]采用半經(jīng)驗公式和POT極值分布模型分別求得2050年美國沿岸海平面高度和風(fēng)暴潮增水變化，兩者求和得到美國沿岸極端海平面高度變化。Menendez等[11]采用動態(tài)廣義極值分布模型預(yù)估了2020年極端海平面高度季節(jié)變化情況。汪楊駿等[12]采用耿貝爾極值分布模型借助Arcgis平臺預(yù)估了2050年中國寧波地區(qū)極端海平面高度變化及其造成的損失。

上述研究從多空間、時間尺度開展了極端海平面高度預(yù)測，但仍存在一些問題亟待解決。第一，局地尺度上影響極端海平面高度的因素主要有平均海平面高度(包括陸地的垂向運動和基準(zhǔn)的變化)、浪、潮、涌的變化以及氣候變化引起的極端海平面高度變化，這些要素的數(shù)據(jù)來源不盡相同，數(shù)據(jù)存在較大不確定性，即便簡單相加也會造成誤差增大，且其不確定性也無法衡量。第二，相對于靜態(tài)極值分布，動態(tài)極值分布含有的參數(shù)空間維度更高，如何有效地對高緯度、高非線性的參數(shù)進(jìn)行最優(yōu)估計，尚未有效解決。針對上述問題，本文通過計算未來海平面高度平均分布和浪、潮、涌以及氣候變化引起的極端海平面高度的動態(tài)極值分布，采用Copula函數(shù)求得未來極端海平面高度的聯(lián)合分布，用積分期望作為海平面高度的預(yù)測值；針對傳統(tǒng)的最大似然估計方法參數(shù)估計時難以滿足正則條件的問題[13]，引入Bayes推斷方法進(jìn)行參數(shù)估計以更有效獲取綜合樣本和先驗信息，將Bayes推斷和DREAM方法結(jié)合以求解高維動態(tài)極值分布的最優(yōu)參數(shù)空間。

1 研究方法與技術(shù)途徑

IPCC AR5在其特別報告中指出未來極端海平面高度變化是由平均海平面高度變化與極端海平面增水(如風(fēng)暴潮、天文大潮、浪等)綜合作用而成[14]。

圖1 未來極端海平面高度預(yù)測流程圖Fig. 1 The flow chart of future extreme sea level forecast

方法步驟和技術(shù)流程如圖1所示。首先對原始數(shù)據(jù)(包括相對海平面高度的多模式集合預(yù)估值、陸地垂直運動的多模式集合預(yù)估值以及極端海平面增水年最大觀測資料)進(jìn)行預(yù)處理，通過Copula模型求得平均高海平面的預(yù)測高度及其概率分布，采用動態(tài)極值理論預(yù)測未來給定重現(xiàn)期的極端海平面增水及其概率分布(其中用Bayes推斷對廣義極值分布模型進(jìn)行參數(shù)估計，DREAM算法求得各參數(shù)后驗分布)；最后采用Copula模型聯(lián)立得到未來極端海平面高度及其概率分布從而對極端海平面高度的期望進(jìn)行預(yù)測。

1.1 數(shù)據(jù)處理

采用CMIP5多模式預(yù)估結(jié)果作為相對平均海平面高度預(yù)測值和陸地垂直運動高度預(yù)測值。對于浪、潮、涌及局地氣候變化引起的極端海平面增水，作如下處理：

1)選取研究區(qū)域的驗潮站站點的逐時和逐月海平面高度數(shù)據(jù)。

2)即使是統(tǒng)一驗潮站，其數(shù)據(jù)依然存在重復(fù)記錄、數(shù)據(jù)缺省、驗潮零點的變動、明顯的異常值和不同的采樣間隔等問題。因此要對驗潮站數(shù)據(jù)進(jìn)行誤差訂正。

3)訂正后的驗潮站記錄包含平均海平面高度的變化以及浪、潮、涌及局地氣候變化引起的海平面增水變化。這里采用逐時海平面高度數(shù)據(jù)減去年平均海平面高度數(shù)據(jù)，剔除平均海平面高度變化的影響，將其作為浪、潮、涌及局地氣候變化引起的海平面增水變化數(shù)據(jù)集。

4)從該數(shù)據(jù)集中選取歷年最大值，作為極端海平面增水的歷史觀測資料。

1.2 Copula理論

Copula函數(shù)描述的是時變量間的相關(guān)性，實際上是一類將聯(lián)合分布函數(shù)與它們各自的邊緣分布函數(shù)連接在一起的函數(shù)。Copula函數(shù)最早由Sklar提出，用不同邊緣分布的變量來構(gòu)造聯(lián)合分布[15]。其定義如下：令F為一個n維的邊緣分布函數(shù)集F1，F(xiàn)2，…，F(xiàn)n，對于所有的實數(shù)x，這些邊緣分布存在一個聯(lián)合分布函數(shù)C，即

若F1，F(xiàn)2，…，F(xiàn)n都連續(xù)，則存在唯一的函數(shù)C。

根據(jù)Sklar定理，可以自由地把任意n個邊緣分布函數(shù)構(gòu)成一個n元的聯(lián)合分布函數(shù)。選取9個Copula函數(shù)作為連接函數(shù)，這9個函數(shù)涵蓋了目前應(yīng)用中的大部分，其特性及定義如表1所示。前7個Copula函數(shù)(Clayton，Ali-Mikhai-Haq(AMH)，Gumbel，F(xiàn)rank，Joe，A12，A14)隸屬于阿基米德函數(shù)簇，其中A12, A14由Nelsen在1999提出[16]。Farlie-Gumbel-Morgenstern(FGM)和Gauss為非阿基米德函數(shù)簇。

表1 本文所采用的Copula函數(shù)定義及參數(shù)空間Tab. 1 Definition and parameter domain of the Copula used in this paper

采用最大似然估計法對上述Copula函數(shù)進(jìn)行參數(shù)估計，用AIC[17](Akaike information criterion)方法和Genest-Rivest方法[18]進(jìn)行擬合優(yōu)度檢驗，選取最優(yōu)的Copula函數(shù)。

1.3 平均海平面高度

1.4 浪、潮、涌及局地氣候變化引起的極端海平面增水

1.4.1 動態(tài)極值模型

其中，[a]+=max[a,0]，μ為位置參數(shù)，ψ>0為尺度參數(shù)，ξ為形狀參數(shù)。廣義極值分布模型根據(jù)ξ的不同可分為三種不同的分布形狀：ξ=0為Gumbel分布；ξ>0為Fréchet分布；ξ<0為Weibull分布。

值得注意的是，式(3)中位置參數(shù)μ(t)是時間的函數(shù)，這意味著該廣義極值分布模型有隨時間變化的特性，故可稱為動態(tài)極值分布模型；相比于靜態(tài)極值分布模型，動態(tài)極值分布模型加入了“時變”的因子，使得模型用于對未來狀態(tài)的預(yù)測時更為靈活、準(zhǔn)確。

未知參數(shù)μ(t)的選取，Renard等[21]提供了多種可供選擇的函數(shù)諸如多項式函數(shù)、分段函數(shù)等。水文氣象文獻(xiàn)中，當(dāng)研究極端事件發(fā)生頻率變化趨勢時，μ(t)常取線性或?qū)?shù)-線性模型[22]。因此，本文未知參數(shù)μ(t)定義如下：

通過動態(tài)極值分布模型，用如下公式計算特定重現(xiàn)期R的極端增水。

1.4.2 Bayes推斷模型

根據(jù)式(5)可以求得浪、潮、涌及局地氣候變化引起的極端增水。采用貝葉斯推斷方法求解式(5)中的參數(shù)空間。傳統(tǒng)的參數(shù)估計方法如矩估計法、最大似然估計法、概率加權(quán)矩估計法、L矩估計法等，以樣本為出發(fā)點，把未知參數(shù)θ為未知常數(shù)，在一定的統(tǒng)計模型下作出統(tǒng)計推斷。而Bayes方法首先采用一個先驗分布去描述未知參數(shù)θ，在得到樣本觀測值之后，由樣本觀測值與先驗分布提供的信息得到未知參數(shù)θ的后驗分布。相比于傳統(tǒng)估計方法，Bayes統(tǒng)計推斷方法有以下優(yōu)點：1)一般而言，極值具有較少的觀測數(shù)據(jù)，若直接用較少的數(shù)據(jù)對統(tǒng)計模型中參數(shù)進(jìn)行估計會有較大的偏差。在此情況下，先對未知參數(shù)的分布有一個較全面的認(rèn)識，再結(jié)合較少的極值觀測值，能得到較好的參數(shù)估計。2)由于未知量必然存在不確定性，Bayes推斷方法可以采用概率與概率分布客觀的描述未知量的不確定性。

假設(shè)各觀測值是相互獨立的，則用Bayes估計參數(shù)的函數(shù)關(guān)系可表示如下：

這個夜晚左小龍?zhí)貏e難熬，還有一只不懂事的蚊子在他的左手中指上咬了一口，那一口恰好咬在骨頭斷裂的位置，還不能撓，真是生不如死。有的時候疼好忍，但癢就不好忍了，忍還不能撓是最不能忍的。左小龍在這個時候想起了泥巴。他突然想，不知道這個小姑娘現(xiàn)在怎么樣了。明天應(yīng)該去找找她，告訴她自己受傷了，當(dāng)然，是在見完黃瑩的情況下。

在使用Bayes公式時(式(6)、(7))，遇到兩個問題：第一，關(guān)于先驗分布的選取問題；第二，邊緣密度函數(shù)積分計算的問題。為此，采用非信息的正態(tài)分布作為μ1，μ0，ψ的先驗分布，采用N(0，0.32)作為ξ的先驗分布[21]。而關(guān)于邊緣密度函數(shù)的積分運算，當(dāng)θ維度較高，即使使用數(shù)值積分方法，計算也比較困難，因此引入DREAM(diffe rential evolution adaptive metropolis algorithm)方法[23]解決該這一問題。

1.4.3 DREAM算法

對于復(fù)雜Bayes推斷的后驗分布求解，傳統(tǒng)方法主要采用MCMC(markov chain monte carloalgorithm)方法[24]。Braak等用DE-MC(differential evolution markov chain)方法[25]改進(jìn)了傳統(tǒng)的MCMC方法，這種算法加快了傳統(tǒng)MCMC方法的運算與收斂速度。而DREAM方法在DE-MC方法的基礎(chǔ)上，基于一個群進(jìn)化架構(gòu)，采用自適應(yīng)差分進(jìn)化學(xué)習(xí)策略對傳統(tǒng)的MCMC進(jìn)行改進(jìn)，可以同時運行多個不同的鏈開展全局搜索，且自動調(diào)整目標(biāo)函數(shù)隨機搜索的范圍和方向，大大增強了傳統(tǒng)MCMC方法對于復(fù)雜、多維問題搜索的能力。

DREAM算法的計算步驟如下：

1)采用先驗分布生成初始種群{xi,i=1，…，N}

2)計算概率密度函數(shù)p(xi)，i=1，…，N。

FORi←1,...,NDO (CHAIN EVOLUTION)

3)在每一條鏈i中產(chǎn)生一個候選點yi，

其中，δ意味著產(chǎn)生建議點的對數(shù)，而r1(j)，r2(n)∈{1，…，N}，且對于j，n=1，…，δ，r1(j)≠r2(n)。e和ε的值分別取自Ud(-b,b)和Nd(0,b*)，b和b*為目標(biāo)函數(shù)的一個細(xì)小偏差。而γ為躍遷系數(shù)，取決于δ和d'，d'為更新的維度。

4)采用二項分布式以1-CR的概率，將每一維度上的建議值yi用xi替換。其中CR為遺傳算法中的交換率。若CR=1，則所有維度都更新，即d'=d。

5)計算p(yi)和候選點上的α(yi，xi)，α(yi，xi)為Metropolis比值。

6)若α(yi，xi)可接受，則xi=yi，否則xi保持不變。

END FOR(CHAIN EVOLUTION)

7)利用IQR[26](Inter-Quartile-Range)統(tǒng)計方法剔除離群數(shù)據(jù)鏈。

通過DREAM，可以得到各個參數(shù)的概率密度函數(shù)，用MC方法進(jìn)行抽樣帶入式(5)，即可求得在特定重現(xiàn)期R下未來時間t時極端海平面增水高度。

1.5 極端海平面高度

2 案例分析

基于上述研究思想對廈門2050年極端海平面高度進(jìn)行預(yù)測。

2.1 2050年廈門平均海平面高度預(yù)估

2050年廈門相對海平面高度預(yù)估值如表2所示。由于廈門的陸地垂直運動不明顯，為計算方便假設(shè)由此引起的海平面高度變化為0。因此，可用表2結(jié)果擬合得到2050年廈門相對海平面高度符合正態(tài)分布(結(jié)果通過置信度為99%的k-s檢驗，p-value=0.939 4)，其概率密度函數(shù)表示如下：

表2 多模式對2050年廈門相對海平面高度預(yù)測Tab. 2 Multi-model prediction of relative sea level in Xiamen in 2050

注：① 溫室氣體排放情景，是對未來氣候變化預(yù)估的基礎(chǔ)。IPCC評估報告采用新一代情景稱為“典型濃度目標(biāo)”(representative concentration pathways，RCPs)。這里representative代表許多種可能性中的一種可能性，用 concentration 而不用輻射強迫是要強調(diào)以濃度為目標(biāo)，pathways則不僅僅指某一個量而且包括達(dá)到這個量的過程。4 種情景稱為RCP8.5、RCP6、RCP4.5及RCP2.6。其中，RCP8.5指輻射強迫上升至8.5 W/m2，2100年CO2當(dāng)量濃度達(dá)到約1 370×10-6；RCP6指輻射強迫上升至6 W/m2，2100年CO2當(dāng)量濃度達(dá)到約850×10-6；RCP4.5指輻射強迫上升至4.5 W/m2，2100年CO2當(dāng)量濃度達(dá)到約6 500×10-6；RCP2.6指輻射強迫上升至2.6 W/m2，2100年CO2當(dāng)量濃度達(dá)到約490×10-6。

圖2 樣本鏈?zhǔn)諗壳闆rFig. 2 Convergence of the sample chain

2.2 2050年廈門極端海平面增水預(yù)估

選取廈門驗潮站1954-1997年海平面高度逐時資料(資料源：夏威夷大學(xué)海平面高度中心(UHSLC)：http://ilikai.soest.hawaii.edu/uhslc/datai.html)，資料本身已經(jīng)經(jīng)過訂正，可以直接用于研究。將逐時資料與年平均值相比，剔除平均海平面高度變化的影響后，得到1954-1997年年最大海平面高度歷史觀測資料。將其代入式(3)、(4)的動態(tài)極值模型中，通過Bayes推斷和DREAM算法計算得到各參數(shù)的特性及分布情況。

圖3 各參數(shù)的后驗概率密度Fig. 3 The posterior probability density of the parameters

圖4 各參數(shù)的分布空間Fig. 4 Distribution space of each parameter sample

變量名平均值標(biāo)準(zhǔn)差μ02.43450.1306μ10.00300.0026ψ0.40210.0359ξ-0.01140.0300

2.3 2050年廈門極端海平面高度預(yù)估

將上述求得的參數(shù)樣本通過計算機輸入到式(5)中，并令重現(xiàn)期R=100a，t=2 050a，便可得到2050年廈門百年一遇的海平面極端增水樣本預(yù)估值。而通過式(10)可得到等量的平均海平面高度樣本。引入Copula模型來計算兩者的聯(lián)合概率密度。采用表1所示9種Copula函數(shù)代入計算，通過AIC準(zhǔn)則和Genest-Rivest方法選取了Frank Copula作為連接函數(shù)，其中Genest-Rivest檢驗結(jié)果如圖5所示。

圖5表征的是廈門地區(qū)2050年平均海平面高度與海平面極端增水聯(lián)合概率密度的Frank Copula函數(shù)的Genest-Rivest檢驗結(jié)果。如圖5所示，該FrankCopula函數(shù)理論估計與經(jīng)驗估計擬合較好。絕大部分點都分布在45°對角線附近，因此可以作為廈門2050年平均海平面高度與極端海平面增水概率密度函數(shù)的連接函數(shù)，將其代入式(9)，通過計算機運算得到2050年廈門百年一遇極端海平面高度分布情況如圖6所示。

圖5 Genest-Rivest方法檢驗結(jié)果(廈門)Fig. 5 Test result of Genest-Rivest method in Xiamen

圖6 2050年廈門百年一遇極端海平面高度概率密度Fig. 6 Probability density of a-hundred-year return period extreme sea level in Xiamen in 2050

圖6表征了2050年，廈門遭遇的百年一遇極端海平面高度的分布情況。從圖中可以得到，廈門遭遇百年一遇的極端海平面高度概率密度整體呈現(xiàn)出中間大，兩頭小的“鐘形結(jié)構(gòu)”，發(fā)生的期望為4.92 m，其95%置信區(qū)間為4.40～5.63 m。

廈門地區(qū)現(xiàn)采用的極端海平面預(yù)測模型，是在假設(shè)極端海平面增水不變的情況下，將當(dāng)前測得的百年一遇的極端海平面高度線性疊加上未來平均海平面高度的變化[27]，按照此方法得到2050年廈門地區(qū)極端海平面度應(yīng)為4.72 m。本文模型在考慮了極端海平面增水變化因素后，得到2050年廈門地區(qū)極端海平面高度期望為4.92 m，其95%置信區(qū)間為4.40～5.63 m。

在防波堤的建設(shè)過程中，若設(shè)計標(biāo)高過小則不能有效防范極端高水位事件，從而對當(dāng)?shù)氐慕?jīng)濟與人員安全造成威脅，而設(shè)計標(biāo)高過大則會提高建設(shè)成本。采用動態(tài)極值理論和Copula函數(shù)對極端海平面高度進(jìn)行建模預(yù)測，相比于其他預(yù)測模型，該模型得到的廈門市未來百年一遇的極端海平面高度是一條概率密度曲線。因此結(jié)合risk-cost-benefit的分析方法就可以對不同設(shè)計標(biāo)高下防波堤建設(shè)成本與其所能抵御的極端海平面事件的損失進(jìn)行量化比較，從而實現(xiàn)效益最大化。

3 結(jié) 語

針對傳統(tǒng)方法在分析極端海平面高度變化時存在的問題和局限，提出一種新的未來極端海平面高度預(yù)測模型。通過分析和案例計算得到如下結(jié)論：

1)通過引入動態(tài)極值模型來對極端海平面高度預(yù)測，較好體現(xiàn)了極端海平面高度隨時間變化的特性。

2)引入Bayes方法對模型參數(shù)空間進(jìn)行估計，可以綜合考慮局地水文經(jīng)驗來改善樣本量少帶來的不確定性。同時，得到的結(jié)果是各參數(shù)的分布，而非一個特定的值。能夠盡可能提取樣本的信息以刻畫參數(shù)的不確定性，為衡量誤差提供了一種簡捷的技術(shù)途徑。

3)DREAM算法所具有的高維全局搜索能力，可加速最優(yōu)參數(shù)空間的搜索效率，使構(gòu)造的動態(tài)極值模型可以更為復(fù)雜，所包含的物理統(tǒng)計信息可以更為完備。

4)引入Copula函數(shù)來求解兩個變量的聯(lián)合分布，用期望和概率分布代替單一加權(quán)平均。能夠有效減少誤差，更好地刻畫極端海平面變化的不確定性。

5)利用該預(yù)測模型對廈門地區(qū)2050年極端海平面高度情況進(jìn)行了預(yù)測。結(jié)果顯示，2050年，廈門遭受百年一遇的極端海平面高度的期望為4.92 m，95%的置信區(qū)間為4.40～5.63 m，可將此作為設(shè)計防波堤高度的輔助決策支持。

[1] MARCOS M, TSIMPLIS M N, SHAW A G P. Sea level extremes in Southern Europe[J]. Journal of Geophysical Research, 2009, 114(C1):362-370.

[2] HAIGH I, NICHOLLS R, WELLS N. Assessing changes in extreme sea levels: Application to the English Channel, 1900-2006[J]. Continental Shelf Research, 2010, 30(9):1042-1055.

[3] CHURCH J, HUNTER J, MCINNES K, et al. Sea-level rise around the Australian coastline and the changing frequency of extreme sea-level events[J]. Australian Meteorological Magazine, 2006, 55(4):253-260.

[4] ULLMANN A, PIRAZZOLI P A, TOMASIN A. 2007: Sea surges in Camargue: trends over the 20th century[J]. Continental Shelf Research, 2007, 27(7):922-934.

[5] ABEYISIRIGUNAWARDENA D S, WALKER I J. Sea level responses to climatic variability and change in Northern British Columbia[J]. Atmosphere-Ocean, 2010, 46(3):277-296.

[6] MENéNDEZ M, WOODWORTH P L. Changes in extreme high water levels based on a quasi-global tide-gauge data set[J]. Journal of Geophysical Research, 2010, 115(C10):234-244.

[7] CHURCH J A, WHITE N J. A 20th century acceleration in global sea-level rise[J]. Geophysical Research Letters, 2006, 33(1):313-324.

[8] TEBALDI C, STRAUSS B H, ZERVAS C E. Modelling sea level rise impacts on storm surges along US coasts[J]. Environmental Research Letters, 2012, 7(1):14032-14042.

[9] STERL A, BRINK H V D, VRIES H D, et al. An ensemble study of extreme storm surge related water levels in the North Sea in a changing climate[J]. Ocean Science, 2009, 5(3):369-378.

[10] HUNTER J. A simple technique for estimating an allowance for uncertain sea-level rise[J]. Climatic Change, 2012, 113(2):239-252.

[11] MENENDEZ M, MENDEZ F J, LOSADA I J. Forecasting seasonal to interannual variability in extreme sea levels[J]. Ices Journal of Marine Science, 2009, 66(7):1490-1496.

[12] 汪楊駿,張韌,錢龍霞,等. 海平面上升引發(fā)的極端高水位的頻率風(fēng)險評估模型及其應(yīng)用——以寧波為例[J]. 災(zāi)害學(xué),2016,31(1):213-218. (WANG Yangjun, ZHANG Ren,QIAN Longxia, et al. Model for probabilistic risk assessment in extreme high water level caused by rising sea level and its application———A case study in Ningbo[J]. Journal of Catastrophology,2016,31 (1): 213-218.(in Chinese))

[13] SMITH R L. Extreme value theory based on the r largest annual events[J]. Journal of Hydrology, 1986, 86(1-2):27-43.

[14] FIELD C B, BARROS V, STOCKER T F, et al. Managing the risks of extreme events and disasters to advance climate change adaptation [M]. Special Report of the Intergovernmental Panel on Climate Change, 2012.

[15] SKLAR M. Fonctions de répartitionn dimensions et leurs marges[J]. Publ. Inst. Statist. Univ. Pairs,1960,8:229-231.

[16] NELSEN R B.An introduction to copulas(Lecture notes in statistics)[M]. Springer,1998:1378597.

[17] AKAIKE H. A new look at the statistical model identification[J]. IEEE Trans. Autom: Control, 1974, 19(6):716-723.

[18] GENEST C,RIVEST L P.Statistical inference procedures for bivariate Archimedean copulas[J].Journal of the American Statistical Association,2012,88(423):1034-1043.

[19] BACKUS G A, LOERY T S, WARREN D E. The near-term risk of climate uncertainty among the U.S. states[J]. Climatic Change, 2013, 116(3-4):495-522.

[20] III J P. Statistical inference using extreme order statistics[J]. Annals of Statistics, 1975 (1):119-131.

[21] RENARD B, GARRETA V, LANG M. An application of bayesian analysis and markov chain monte carlo methods to the estimation of a regional trend in annual maxima[J]. Water Resources Research, 2006, 42(12):4407-4413.

[22] BEGUERíA S, ANGULO M M, VICENTE S, et al. Assessing trends in extreme precipitation events intensity and magnitude using non-stationary peaks-over-threshold analysis: a case study in northeast Spain from 1930 to 2006[J]. International Journal of Climatology, 2011, 31(14):2102-2114.

[23] VRUGT J A, DIKS C G H, HYMAN J M, et al. Accelerating markov chain monte carlo simulation by differential evolution with self-adaptive randomized subspace sampling[J]. International Journal of Nonlinear Sciences & Numerical Simulation, 2009, 10(3):273-290.

[24] METROPOLIS N, ROSENBLUTH A W, ROSENBLUTH M N, et al. Equation of state calculations by fast computing machines[J]. Journal of Chemical Physics, 1953, 21(6):1087-1092.

[25] BRAAK C J F T, VRUGT J A. Differential evolution markov chain with snooker updater and02fewer chains[J]. Statistics & Computing, 2008, 18(4):435-446.

[26] UPTON G, COOK I. Understanding statistics[M]. Oxford:Oxford University Press,1996:55.

[27] BERNARDO I J M, BERGER J O, DAWID A P, et al. Inference from iterative simulation using multiple sequences[J]. Statistical Science, 1992, 7(4):457-472.

[28] 李曉剛. 廈門市海平面上升規(guī)劃對策[J]. 現(xiàn)代城市研究, 2008(5):27-33.(LI Xiaogang. Plan countermeasueabout sea level rising of Xiamen[J].Journal of Modern Urban Research, 2008(5):27-33.(in Chinese))

Prediction modeling of extreme sea level based on dynamic extreme value theory and Copula function

WANG Yangjun1， HOU Taiping1， ZHANG Ren1,2， QIAN Longxia1， WANG Feng3， LONG Qiang3

(1. Institute of Meteorology and Ocean, PLA University of Science and Technology, Nanjing 211101, China; 2. Collaborative Innovation Center on Forecast Meteorological Disaster Warning and Assessment, Nanjing 210044, China;3. Meteorological Bureau of Caofeidian Industrial Zone, Tangshan 063000, China)

1005-9865(2016)04-0062-09

P731.23

A

10.16483/j.issn.1005-9865.2016.04.009

2015-08-04

氣象水文預(yù)先研究項目(407010602)；唐山市曹妃甸工業(yè)區(qū)專項(CQZ-2014001)

汪楊駿(1990-)，男，碩士研究生，從事海洋環(huán)境風(fēng)險評估研究。E-mail：492670449@qq.com

張 韌，教授，博士生導(dǎo)師，從事海洋環(huán)境影響評估與決策研究。E-mail: zrpaper@163.com