南海-北印度洋波高的長期變化：II.趨勢的區域性、季節性差異

高占勝，鄭崇偉，，李訓強，劉文鵬

（1.海軍大連艦艇學院航海系，遼寧大連116018；2.解放軍理工大學氣象海洋學院，江蘇南京211101）

?

南海-北印度洋波高的長期變化：II.趨勢的區域性、季節性差異

高占勝1，鄭崇偉1，2，李訓強2，劉文鵬1

（1.海軍大連艦艇學院航海系，遼寧大連116018；2.解放軍理工大學氣象海洋學院，江蘇南京211101）

摘要：利用ERA-40海表10 m風場驅動WAVEWATCH-III（WW3）海浪模式，模擬得到南海-北印度洋1957年9月—2002年8月的海浪資料，采用一元線性回歸方法，分析了該海域有效波高的長期趨勢，以期為研究全球氣候變化、波浪能資源開發、防災減災等提供參考。結果表明：（1）1958—2001年期間，該海域的SWH有線性遞增趨勢，遞增率為0.0017 m/a，且變化趨勢表現出很大的區域性差異：僅部分小范圍海域呈顯著性遞減，其余大部分海域的SWH呈顯著性遞增；（2）SWH的變化趨勢存在較大的季節性差異：各個季節呈顯著性遞減的區域范圍都較小；低緯度的遞增趨勢主要體現在春季和冬季，尤其冬季幾乎整個南海-北印度洋的SWH均呈顯著性遞增趨勢；索馬里以東一近似圓形海域的遞增趨勢主要體現在夏季。

關鍵詞：WAVEWATCH-II海浪模式；線性趨勢；區域性差異；季節性差異

1　引言

海洋對人類的生產生活［1-3］、以及軍事行動［4］都有著重要影響，隨著“一帶一路”戰略正式展開，南海-北印度洋的重要性愈發凸顯［5-7］。有效波高（Significant Wave Height，SWH）的長期變化趨勢對全球氣候變化、波浪能資源開發、防災減災等都有著重要的影響。前人對在資料匱乏的情況下對SWH的長期變化趨勢做了較多工作，Gulev等［8］利用綜合海洋大氣數據集（International Comprehensive Ocean-Atmosphere Data Set，ICOADS）資料，分析發現1958—2002年期間，北大西洋和北太平洋中緯度海域的SWH在冬季呈顯著性遞增，趨勢為10—40 cm/decade。Dodet等［9］發現在1953—2009年期間，大西洋東北部海域的SWH以20 cm/ decade的速度顯著性遞增。Semedo等［10］研究發現1957—2002年期間，北大西洋和北太平洋大部分海域的SWH表現出顯著的遞增趨勢。鄭崇偉等［11-12］研究發現1988—2009年期間，中國海的SWH和波浪能流密度均呈顯著性逐年線性遞增。鄭崇偉等［13］在國內首次對南海-北印度洋SWH的長期變化趨勢展開研究，發現在1979—2014年期間，整個南海SWH的遞增趨勢強勁，大部分海域在0.4 cm/（s·a）以上，大值區分布于南海北部、臺灣東南部近海，在0.6 cm/（s·a）以上。北部灣北部、泰國灣北部近海的SWH沒有顯著的變化趨勢。孟加拉灣、印度半島西部大范圍海域的SWH沒有表現出顯著的變化趨勢。在阿拉伯海西部海域、印度洋15°S—0°大范圍海域表現出顯著性遞增，趨勢為0.1—0.4 cm/（s·a）。僅部分零星海域的SWH表現出顯著性遞減趨勢。

南海-北印度洋對我國至關重要，它不僅僅是重要國際航線，也是我國亞丁灣護航的必經之地、重要的能源通道之一。本文利用ERA-40海表10 m風場驅動WAVEWATCH-Ⅲ（WW3）海浪模式，模擬得到南海-北印度洋1957年9月—2002年8月的海浪資料，采用一元線性回歸方法，分析了該海域SWH的線性趨勢、以及變化趨勢的區域性和季節性差異，以期可為研究全球氣候變化、波浪能資源開發、防災減災等提供參考。

2　資料與方法

本文利用模擬得到的近45 a（1957年9月—2002年8月）的海浪數據，計算了南海-北印度洋SWH的長期變化趨勢，主要計算了該海域SWH的整體變化趨勢、變化趨勢的區域性差異。數據來源：以ERA-40海表10 m風場驅動第三代海浪數值模式WW3，模擬得到南海-北印度洋近45 a逐3 h的海浪場數據。計算范圍為10°S—30°N，40°—140°E，地形分辨率取0.5°×0.5°，計算時間從1957年9月1日00:00時（世界時，下同）—2002年8月31日18:00時。積分步長900 s，每3 h輸出一次結果，模擬的海浪資料的已經過與ICOADS資料對比檢驗［7］，精度較高。ERA-40海表10 m風場來自歐洲中期天氣預報中心（European Centre for Medium-Range Weather Forecasts，ECMWF），其時間分辨率為6 h，空間分辨率為2.5°×2.5°。該數據的空間范圍覆蓋全球，時間范圍從1957年9月1日00:00時—2002 年8月31日18:00時。

3　南海-北印度洋SWH的長期變化趨勢

3.1 SWH的整體變化趨勢

將整個南海-北印度洋的SWH從1958—2001年進行逐年區域平均，分析近44 a期間該海域SWH的整體變化趨勢，見圖1。長期變化趨勢的具體計算方法如下：

回歸方程

由圖1可見，變化趨勢的相關系數|r|=0.63＞r0.05= 0.29，通過了95%的信度檢驗，表明該海域SWH的線性遞增趨勢顯著，回歸系數為0.0017，即1958—2001年期間，南海-北印度洋海域的SWH整體上以0.0017 m/a的速度顯著性逐年線性遞增。在1958—1969年的10余年期間，SWH的波動較為劇烈，在1969—1978年的10a期間表現出較為明顯的遞增趨勢（0.0076 m/a），1979—1986年期間的變化較為平緩，1987—1990年期間為第二波遞增趨勢顯著的時間段（0.035 m/a），1990—2001年期間便顯出緩慢的遞減趨勢。從圖1不難發現，SWH的長期遞增趨勢較為緩慢，梅勇等［14］對該海域的海表風速、有效波高進行過EOF分析，發現二者的第一模態時間序列都呈遞增趨勢，增速緩慢，本文的發現與梅勇等［14］吻合。鄭崇偉等曾對南海重要島礁的波候特征、海浪發電的可行性展開過論證［15-16］，并率先對第一島鏈的波候（海浪氣候態）特征展開了系統性研究［17］，發現該島礁的SWH呈遞增趨勢。郭隨平等［18］曾發現南海的海浪場變化與厄爾尼諾存在相關，在以后的工作中，可借鑒該方法，深入分析南海-北印度洋SWH長期變化趨勢的物理機制。

圖1 1958—2001年南海-北印度洋的逐年平均有效波高及其變化趨勢

3.2 SWH變化趨勢的區域性差異

近年來，對大洋SWH變化趨勢的研究已有不少，但多是區域平均后分析整體變化趨勢，很難體現出不同海域變化趨勢的區域性差異，這種分析方法往往會將某些重點海域的變化特征忽略，本文將南海-北印度洋0.5°×0.5°逐網格點上的SWH從1958—2001年取逐年平均，分析每個網格點上SWH的逐年線性趨勢，分析變化趨勢的區域性差異，見圖2。

近44 a期間，南海-北印度洋大部分海域的SWH表現出顯著性遞增趨勢，主要有蘇門答臘島-索馬里一帶的大范圍海域、呂宋海峽-中南半島東南海域一帶海域、呂宋島東北部海域；僅部分小范圍海域呈顯著性遞減趨勢，其余海域則無明顯的變化趨勢。遞增趨勢較強的海域主要分布于：呂宋島東北部近海（0.003—0.005 m/a）、呂宋海峽（0.002—0.004 m/a）、中南半島東南海域（即傳統的南海大風區，趨勢為0.003—0.004 m/a）、斯里蘭卡至蘇門答臘島一帶的大范圍海域（0.003—0.006 m/a）、索馬里以東的一橢圓形海域（0.002—0.005 m/a）；SWH在泰國灣、北部灣、孟加拉灣的遞減趨勢基本在-0.002 m/a左右，爪哇海的遞減趨勢較為強勁，約-0.003—-0.005 m/a。

3.3 SWH變化趨勢的季節性差異

本文還將SWH從1958—2001年進行逐季節平均，分析逐春季、逐夏季、逐秋季、逐冬季的變化趨勢，比較SWH變化趨勢的季節性差異，見圖3。

逐春季：近44 a期間，印度洋低緯海域、阿拉伯半島東部近海、印度半島15°N西部近海、琉球群島附近海域的SWH表現出顯著性遞增趨勢，尤其是在印度洋80°—100°E低緯海域，遞增趨勢較為強勁，達到0.004—0.007 m/a；SWH呈顯著性遞減的范圍較小，主要分布于：北部灣（約-0.002 m/a）、泰國灣（-0.003—-0.001 m/a）、加里曼丹島西部近海（約-0.001 m/a）、加里曼丹島北部近海（-0.004—-0.002 m/a）、亞丁灣（約-0.002 m/a），以及一些零星海域，詳見圖3a。

逐夏季：SWH呈顯著性遞增的區域主要分布于蘇門答臘島-索馬里一帶、呂宋島東部大范圍海域，以及一些散亂的小范圍海域，遞增趨勢較強的海域分布于：蘇門答臘島-斯里蘭卡之間一橢圓形海域（0.004—0.006 m/a）、索馬里東部一圓形海域（0.004—0.008 m/a）；呈顯著性遞減的區域主要分布于：幾乎整個孟加拉灣（-0.006—-0.004 m/a）、加里曼丹島南部和西南部海域（-0.008—-0.01 m/a），以及一些零星海域，尤其是加里曼丹島南部的爪哇海，遞減趨勢尤為強勁，達到-0.012 m/a左右，詳見圖3b。

圖2 1958—2001年南海-北印度洋有效波高的變化趨勢（單位：10-3m/a），彩色區表示通過了95%的信度檢驗

圖3 1958—2001年期間不同季節有效波高變化趨勢（單位：10-3m/a），彩色區表示通過了95%的信度檢驗

逐秋季：SWH呈顯著性遞增的區域主要分布于蘇門答臘島-索馬里一帶、呂宋海峽及其東部近海，以及一些零星海域，遞增趨勢較強的海域分布于：斯里蘭卡南部一橢圓形海域（約0.004 m/a）、呂宋海峽及其東部近海（約0.006 m/a）；呈顯著性遞減的區域主要分布于：孟加拉灣中部和東北部海域（約-0.004 m/a）、加里曼丹島南部和西南部海域（-0.006— -0.002 m/a），以及一些零星海域，尤其是加里曼丹島南部的爪哇海，遞減趨勢較為強勁，達到-0.006 m/a左右，詳見圖3c。

逐冬季：幾乎整個南海-北印度洋海域的SWH均呈顯著性遞增，基本在0.002 m/a以上，遞增趨勢較為強勁的海域分布于：印度洋低緯度中東部海域（約0.008 m/a）、呂宋海峽（0.008—0.01 m/a）、呂宋島東北部近海（0.006—0.01 m/a）；僅部分小范圍零星海域的SWH表現出顯著性遞減趨勢，詳見圖3 d。

對比不難發現不同海域SWH的變化趨勢由不同季節占主導。

4　結論

本文利用ERA-40海表10 m風場驅動WW3海浪模式，模擬得到南海-北印度洋1957年9月—2002年8月的海浪資料，采用一元線性回歸方法，分析了該海域SWH的長期趨勢，主要計算了各個區域的逐年變化趨勢，以及各個季節的變化趨勢，結果表明：

（1）1958—2001年期間，該海域的SWH有線性遞增趨勢，遞增率為0.0017 m/a，且變化趨勢表現出很大的區域性差異：僅部分小范圍海域呈顯著性遞減趨勢，其余大部分海域的SWH呈顯著性遞增，其中遞增趨勢較強的海域集中分布于呂宋島東北部近海（0.003—0.005 m/a）、呂宋海峽（0.002—0.004 m/a）、南海大風區（0.003—0.004 m/a）、斯里蘭卡至蘇門答臘島一帶的大范圍海域（0.003—0.006 m/a）、索馬里以東的一橢圓形海域（0.002—0.005 m/a）；

（2）近44 a期間，南海-北印度洋海域SWH的變化趨勢存在較大的季節性差異，不同海域SWH的變化趨勢由不同季節主導。各個季節呈顯著性遞減的區域范圍都較小，僅夏季的范圍略大于其余季節、遞減的強度也稍強于其余季節。低緯度的遞增趨勢主要體現在春季和冬季，尤其冬季幾乎整個南海-北印度洋的SWH均呈顯著性遞增趨勢；索馬里以東一近似圓形海域的遞增趨勢主要體現在夏季。

參考文獻：

［1］高占勝.臺風“菲特”期間T639預報風場在東中國海的有效性檢驗［J］.海洋學研究，2015，33（1）:9-15.

［2］高占勝.Grads在中國近海大浪頻率統計中的運用［J］.海洋科學前沿，2014，1:50-54.

［3］許富祥.海浪預報知識講座.海洋預報，2001，18（1）:82-84.

［4］鄭崇偉，潘靜，黃剛.利用WW3模式實現中國海擊水概率數值預報［J］.北京航空航天大學學報，2014，40（3）:314-320.

［5］鄭崇偉，潘靜，孫威，等.經略21世紀海上絲路之海洋環境特征系列研究［J］.海洋開發與管理，2015，32（7）:4-9.

［6］鄭崇偉，李訓強，高占勝，等.經略21世紀海上絲路之海洋環境特征：風候統計分析［J］.海洋開發與管理，2015，32（8）:4-11.

［7］鄭崇偉，高占勝，張雨，等.經略21世紀海上絲路之海洋環境特征：極值風速和極值波高［J］.海洋開發與管理，2015，32（11）:4-11.

［8］Gulev S K，Grigorieva V.Variability of the Winter Wind Waves and Swell in the North Atlantic and North Pacific as Revealed by the Voluntary Observing Ship Data［J］.Journal of Climate，2006，19 （21）:5667-5685.

［9］Dodet G，Bertin X，Taborda R.Wave Climate Variability in the North-East Atlantic Ocean over the Last Six Decades［J］.Ocean Modelling，2010，31（3-4）:120-131.

［10］Semedo A，Su?elj K，Rutgersson A，et al.A Global View on the Wind Sea and Swell Climate and Variability from ERA-40［J］. Journal of Climate，2011，24（5）:1461-1479.

［11］鄭崇偉，李訓強.基于WAVEWATCH-III模式的近22年中國海波浪能資源評估［J］.中國海洋大學學報（自然科學版），2011，41 （11）：5-12.

［12］Zheng C W，Zhuang H，Li X，et al.Wind Energy and Wave Energy Resources Assessment in the East China Sea and South China Sea［J］.Science China Technological Sciences，2012，55 （1）:163-173.

［13］鄭崇偉，付敏，芮震峰，等.經略21世紀海上絲路之海洋環境特征：波候統計分析［J］.海洋開發與管理，2015，32（10）:1-7.

［14］梅勇，宋帥，周林.北印度洋—南海海域海浪場、風場的年際變化特征分析［J］.海洋預報，2010，27（5）:27-33.

［15］鄭崇偉，李崇銀.中國南海島礁建設：重點島礁的風候、波候特征分析［J］.中國海洋大學學報：自然科學版，2015，45（9）:1-6.

［16］鄭崇偉，李崇銀.中國南海島礁建設：風力發電、海浪發電［J］.中國海洋大學學報：自然科學版，2015，45（9）:7-14.

［17］鄭崇偉，陳璇，李崇銀.朝鮮半島周邊海域波候觀測分析［J］.中國海洋大學學報：自然科學版，2015，45（9）:21-27.

［18］郭隨平，莊卉，鄭崇偉，等.南海海浪場與厄爾尼諾的相關性分析［J］.海洋預報，2012，29（6）:37-43.

中圖分類號：P731.22

文獻標識碼：A

文章編號：1003-0239（2016）02-0039-00

DOI:10.11737/j.issn.1003-0239.2016.02.006

收稿日期：2015-07-12

基金項目：國家自然科學基金項目（41205073，41275099）

作者簡介：高占勝（1965-），男，副教授，碩士，主要從事船舶操縱的理論和方法及軍事海洋環境方面的研究。E-mail:1310615779@qq.com

Long-term trend of SWH from year 1958 to 2001 in the South China Sea and the North Indian Ocean

GAO Zhan-sheng1，ZHENG Chong-wei1，2，LI Xun-qiang2，LIU Wen-peng1
（1.People's Liberation Army Dalian Naval Academy，Dalian 116018 China；2.College of Meteorology and Oceanography，People's Liberation Army University of Science and Technology，Nanjing 211101 China）

Abstract：In this study，the wave in the South China Sea and North Indian Ocean from was simulated from December 1957 to August 2002 by using WAVEWATCH-Ⅲ（WW3）wave model with ERA-40 wind data as the driving field.The long-term trend of SWH was analyzed based on the linear regression method to provide references for studying the global climate change，wave energy resource development and disaster prevention and reduction.Results show that:（1）During 1958 to 2001，the averaged SWH in the study area exhibits a significant increasing trend，with a rate of 0.0017m/a.The change trend of SWH also has obvious regional difference，increasing in most areas，while decreasing in a small-scale area.（2）The variation trend of SWH has a significant seasonal difference.The scale of waters with decreasing trend of SWH is small in each season.The increasing trend of SWH in lower latitudes mainly reflected in the spring and winter，especially in the winter almost the entire South China Sea-North Indian Ocean SWH showed a statistically significant increasing trend.

Key words：WAVEWATCH-Ⅲwave model；linear trend；regional difference；seasonal difference