海 洋 学 报
ACTAOCEANOLOGICASINICA
Vol131,No11January2009
太平洋海平面变化特征及影响因素分析
顾小丽1,李培良1*
(11中国海洋大学物理海洋实验室,山东青岛266100)
摘要:采用经验模态分解法对太平洋沿岸验潮站的月平均海平面资料进行处理,结合T/P高度计
资料、Church(2004)重构SSH资料、Ishii(2005)月均海温资料,研究太平洋海平面年际、年代际变化以及资料长度内海平面变化趋势。太平洋沿岸海平面总体呈上升趋势,平均上升速率为114mm/a,趋势项分布有明显的区域性和纬度特征。ENSO对太平洋地区海平面年际变化有显著影响,海平面年际变化与Nino3指数在西(东)太平洋为负(正)相关,海平面年际变化与Nino3指数的相关性在热带太平洋最大,并随着纬度升高相关性减弱,且不同地区年际变化有滞后ENSO时间不等的最大相关。海平面年际变化与PDO指数在西(东)太平洋为负(正)相关,海平面与PDO的相关性分布有区域性和随时间演变特征。年代际变化对目前使用卫星高度计资料分析海平面长期趋势项的预测有直接影响,可能完全掩盖海平面长期变化趋势。关键词:海平面变化;年际变化;年代际震荡;ENSO
中图分类号:P731 文献标识码:A
文章编号:0253-4193(2009)01-0028-09
1 引言
海平面是上层海水特征和海洋环流的重要指示,海平面的变化对海洋大气系统以及人类的生产
生活都有着重要影响。近年来由于温室效应加剧,导致的海平面加速上升,使海平面变化成为社会关注的焦点问题。太平洋海平面变化受其特有的地理环境、海气相互作用、全球变暖大背景等因素的影响,表现出特殊的规律和特征。
近年来中外学者利用验潮站资料、卫星高度计资料对海平面长期变化趋势进行了深入的研究。研究结果表明,太平洋海平面在过去十几到几十年间总体趋势是上升的,且在近十几年有加速上升的现象,海平面长期变化趋势有显著的区域性分布特征[5],而海平面的年际和年代际变化对使用短期资料(特别是卫星高度计资料)分析海平面的
收稿日期:2008-08-18;修订日期:2008-12-30。
[1-4]
长期变化趋势有显著的影响。
Daniel等[6],Wu等[7],Albert等[8]的研究表明太平洋地区的气候演变有明显年际、年代际变化特征。这种年际和年代际变化在海温场、大气压强、大气环流等方面都有显著表现,对季风和降水有重大影响[9]。研究表明太平洋年际变率与ENSO相关[10-12],年代际变化与PDO有关[8-9]。Smith[13],Han等[14],Yoshi等[15]的研究表明热带太平洋海平面也存在显著的年际和年代际变化,而ENSO和PDO是海平面年际和年代际变化的重要影响因素。作为高频变化的重要背景,海平面年代际变化对高频变化有重要的调制作用。而年际和年代际变化影响着海平面低频变化和长期趋势的预测。吕俊梅等、Albert等认为PDO在热带太平洋的叠加使ENSO循环表现出年代际的低频震荡特征,该调制作用是近年来ENSO现象频发、强度较大的
[17]
[8]
[16]
资助项目:教育部科学技术重点项目(108159);国家自然科学基金资助项目(40506006;40676013);国家重点基础研究发展计划资助项目(2007CB411807);国家高技术研究发展计划(/8630计划)资助项目(2007AA0912012)。
作者简介:顾小丽(1984)),女,江苏省宿迁市人,在读硕士,从事海平面变化研究。E-mail:ccho.dgl@gmail.com*通讯作者:李培良(1975)),博士,副教授。E-mail:lpliang@ouc.edu.cn
1期 顾小丽等:太平洋海平面变化特征及影响因素分析29
重要因素。Cazenave等[18]利用最近15a的高度计资料得到热带西太平洋海平面有显著的上升趋势,而热带东太平洋下降明显,该结果与过去50a的热比容高度分布相反,太平洋海平面年代际震荡的空间分布[19]是影响这一结果的重要因素。因而,海平面年代际变化对海平面长期趋势项的准确预测有非常大的影响,对其深入的研究有重要意义。
本文利用经验模态分解(EMD)[20-21]方法处理PSMSL(PermanentServiceforMeanSeaLevel)的验潮站资料,研究太平洋地区海平面长期变化趋势,海平面年际、年代际变化的分布特征、主要影响因素以及年代际变化对海平面变化趋势的影响。
W),西太平洋(10bS)10bN,140bE)180b)(图中蓝色长方形区)作为计算重构海面高度和比容海平面在热带太平洋东、西区域平均的范围。为研究海平面的年
际和年代际变化,对海平面时间序列取距平,然后进行35个月或10a的低通滤波。计算比容海面高度的温度场为1945-2003年Ishii(2005)月均三维海温资料,空间分辨率为1b@1b。采用Church(2004)CSIRO提供的海表高度重构资料,空间分辨率为1b@1b。海面高度资料为1993-2007年共15a的T/P和Jason1联合高度计资料。取Nino3指数代表ENSO现象的海洋效应(http://www.cpc.ncep.noaa.gov/data/indices/),PDO指数来源于http://jisao.wash-ington.edu/pdo/PDO.latest。
212 EMD方法
EMD方法的基本思路是:如果一个原始的数据序列x(t)中极大值或极小值点的数目比上跨零点或下跨零点的数目多两个或两个以上,那就对该数列进行平稳化处理,产生一系列具有不同特征尺度的数据序列,称固有模态函数(intrinsicmodefunc-tion)。EMD分解为一个从高频到低频的极值筛选过程,分解得原始数据的一系列IMF分量c1,c2,c3,,,cn与原数据序列的关系如下: x(t)=
t=1
2 资料来源及数据处理
211 资料来源
验潮站资料为PSMSL提供的月均、年均海面RLR(revisedlocalreference)数据,文中选取了太平洋海区27a以上的240个站点的资料,站点分布见图1。计算比容海面高度的温度场为1945)2003年Ishii(2005)月均三维海温资料,空间分辨率为1b@1b。采用Church(2004)CSIRO(theCommonwealthScien-tificandIndustrialResearchOrganisation)提供的海表高度重构资料,空间分辨率为1b@1b。海面高度资料为1993)2007年共15a的T/P和Jasonl联合高度计资料。取Nino3指数代表ENSO现象的海洋效应(http://www.cpc.ncep.noaa.gov/data/indices/),PDO指数来源于http://jisao.washington.edu/pdo/PDO.latest。选取东太平洋(10bS)10bN,80b~120b
Ec
n
t
+rn,(1)
式中,rn(t)为原始数据序列整个时间范围内的固有适定单调函数,或是最多一个极值的函数,称为总体自适定趋势,以下文中提到的趋势项都指该总体
自适定趋势。
[21]
图1 太平洋验潮站海平面变化趋势分布
蓝色长方形区域为热带东、西太平洋区域平均范围;暖色系点为海平面上升站点,冷色系点为海平面下降站点30海洋学报 31卷
3 太平洋海平面长期趋势变化特征
EMD方法可以很好的提取时间序列的非线性变化趋势,本文对太平洋沿岸的240个验潮站资料作EMD分解,得到各点海平面变化趋势。太平洋地区海平面总体呈上升趋势,其速率平均值为114mm/a,与Church等[1]全球海平面在过去50~100a间上升速率约1~2mm/a的结论较一致。各个站点的长期变化趋势空间分布极不平衡,呈现出西部上升快、东部上升不显著的区域分布特征(见图1)。一半以上(约57%)的站点海平面变化趋势在-2和2mm/a之间,而少数站点表现出极大的上升或下降趋势。
西太平洋区域平均趋势的线性上升值为1190mm/a。中高纬(25bN以北)日韩附近海平面平均趋势为2117mm/a,中低纬(25bS)25bN)海平面平均趋势为2133mm/a,两者都大于西太的总体平均,上升站多于下降站点。朝鲜半岛、日本群岛附近海平面升降趋势较复杂,大部分站点海平面上升,日本东部由于受黑潮及其延伸体以及陆地沉降等因素的影响海平面上升很快,个别站点受地壳运动影响表现为下降。中低纬暖池附近海域海平面上升趋势比较显著,这与赤道西太平洋存在的海面高度变化速率的高值中心相吻合。澳大利亚附近海区下降站点居多,且由于站点稀疏,海面变化趋势代表性不强,很难反应总体特征。
东太平洋海域海平面上升趋势不明显,其区域平均约0120mm/a。北部中高纬地区总体趋势在下降,平均年变率约-0143mm/a,这与该区域地壳上升有关[23]。处于高纬的阿拉斯加和阿留申地区海平面下降尤为显著,区域平均达-2113mm/a,是验潮站海平面下降最大地区,太平洋验潮站海平面下降的极大值点几乎都分布于该区。东太平洋中低纬地区平均海平面变化速率为1191mm/a,由于此区站点较稀疏、区域范围较大,平均变化趋势代表性不佳。南美沿岸验潮站较少,但大都表现为下降趋势。
[22]
在年际变化特征(图2,3中蓝色实线,其他未画出站点也有类似的特征),年际变化的振幅可达150mm。网格点上和区域平均的重构海表高度也表现出年际变化特征(图略),由此可见,太平洋海平面有明显的年际变化特征。海平面的年际变化与ENSO的发生发展关系密切。
图2为热带太平洋验潮站海平面变化与Nino3指数的关系,其中图2a-d为热带西太平洋站点,图2e-h为热带东太平洋站点,图中数据均经35个月的低通滤波处理。由图可见热带西太平洋海平面与Nino3指数负相关,二者的相关系数可达-0191。热带东太平洋海平面与Nino3指数同步相关,不同站点二者相关性大小不等(在0142~0187之间)。验潮站海平面资料与Nino3指数呈现出东(西)太平洋二者正(负)相关的关系,东、西区域平均(图1蓝线区域)的Church等[24]重构海表面高度资料与Nino3指数也有相似的相关性。Church等重构海面高度资料(见图4中蓝色和黑色线)在热带西(东)太平洋与ENSO负(正)相关为0178(-0178)。图2)5都表现出在强ENSO年海平面与Nino3指数相关更为明显、海平面变化振幅增大的特征,其最大振幅差可达300mm。由此可见,ENSO对热带太平洋海平面的年际变化影响显著。
图3为北太平洋中高纬度验潮站海平面年际变化与Nino3指数的关系,其中图3a-d为西太平洋站点,图3e-h为东太平洋站点。与热带太平洋相似,西(东)太平洋中高纬地区海平面年际变化与Nino3指数负(正)相关,然而二者的相关性远小于热带地区,相关系数在-012到-0154(0160~0182)之间,且海平面年际变化表现出滞后Nino3指数1个月到1a不等的相关。比较图2、图3不难看出,中高纬度地区验潮站海平面的年际变化强度小于热带低纬度地区,热带太平洋地区年际变化振幅可达140mm,而中高纬站点的年际变化振幅不到100mm。由此可见,海平面的年际变化在热带地区更为显著,ENSO对低纬度地区的影响较大,随纬度升高,影响减弱。在中高纬度海区,海平面变化与ENSO的发生发展表现出一定的超前滞后关系。
[24]
[16]
[24]
4 太平洋海平面年际变化特征
由于南太平洋站点较少,将以20bS以北海区作进一步研究,另外由于验潮站较多难以一一画出,本文将取部分有代表性的站点绘图。将验潮站资料进行35个月的低通滤波,可见太平洋验潮站海平面存5 太平洋海平面年代际变化特征
由图3验潮站资料可看出在年际变化曲线上叠加有年代际变化的信号,其中以图3a,c,f,g,h尤为明显,振幅可达到50mm。1期 顾小丽等:太平洋海平面变化特征及影响因素分析31
众多观测和研究结果显示1977年前后北太平洋气候出现一次显著突变,该现象引起了学者对年代际变化的认识。由于太平洋年代际涛动对气候的影响、对ENSO循环的调制作用以及对海平面长期趋势项的影响[6,8],年代际变化受到越来越多的关注。吕俊梅等的研究表明PDO是太平洋地区年代际变化的主导模态(图5,7中的蓝色实线),ENSO现象因受PDO的调制作用,表现出年代际变化特征(图5,7中的青色虚线)。太平洋海平面变化有明显的年代际特征,该年代际变化在验潮站资料(图5,图6中红色实线)、重构的海表高度资料(图7蓝色实线)和比容海面高度(图7红色实线)中均有反映。对于验潮站资料,年代际变化振幅可达50mm,重构场区域平均的年代际振幅也达25mm。
海平面的年代际变化与PDO的影响有关,图5
[17]
为热带太平洋验潮站海平面年代际变化与PDO指数的关系,图中a-d站点位于西太平洋沿岸,e-h站位于东太平洋沿岸,数据均经十年低通滤波处理。由图可见,热带西太平洋海平面年代际变化与PDO指数表现为负相关,且部分地区海平面有滞后PDO指数约1a的最大负相关。热带东太平洋海平面年代际变化与PDO指数为显著正相关(图5e-h),部分地区出现前者滞后后者一年左右的最大正相关。热带太平洋地区重构的海面高度场和比容海平面也有与PDO相关的年代际变化现象(见图7),且对于重构场,海平面年代际变化滞后PDO一年左右,在西(东)太平洋达到最大负(正)相关。这说明太平洋年代际涛动对热带太平洋海平面有显著影响,海平面的响应变化滞后PDO一年左右。由于卫星高度计资料时间长度有限,难以分辨年代际变化,在此不作分析。
图2 热带太平洋(25bS)25bN)验潮站海平面年际变化asl与Nino3指数(ini)
蓝色实线为验潮站海平面变化,红色虚线为Nino3指数。图中a,b,c、d位于热带西太平洋(站点经纬度依次为13109bN,123145bE;10118bN,123155bE;9131bN,138108bE;8144bN,167144bE);e,f,g,h位于热带东太平洋(站
点经纬度依次为8158bN,79134bW;2112bN,80155bW;0145bN,90119bW;23138bN,70124bW)
32海洋学报 31卷
图6为北太平洋中高纬地区验潮站海平面年代际变化与PDO的关系,图中a-d站点位于西北太平洋沿岸,e-h站位于东北太平洋沿岸。由图可见,西(东)北太平洋海平面年代际变化与PDO指数表现为显著的负(正)相关,二者的相关系数可达-0150(0190)。验潮站海平面的年代际变化与PDO的相关性有区域性分布和随时间演变的特征:不同
地区二者相关性大小不同,不同时期海平面超前、滞后PDO的关系不同。总体上,西太平洋海平面年代
际变化与PDO指数呈反相变化;然而在中国沿海和日本海附近站点二者负相关不显著,年代际变化强度较小(振幅在40mm以下),这些地区海平面变化可能还受其它因素的影响;在韩国和北海道附近站点甚至表现为前者滞后后者一年左右的正相关。日
1期 顾小丽等:太平洋海平面变化特征及影响因素分析33
34海洋学报 31卷
本太平洋沿岸站点海平面的年代际变化随纬度变化,35bN以北(南)海平面年代际变化超前(滞后)PDO指数1a左右达最大正(负)相关。东太平洋的阿拉斯加地区,海平面年代际变化振幅20~40mm,主要表现为超前PDO指数一年左右的最大正相关,向南二者总体为正相关,振幅增大,于2000年之前前者超前后者5个月左右的达最大正相关,2000年之后为同步相关。1977年到2000年间,PDO变化有两次极大值(图6中蓝色实线),西(东)太平洋中高纬度站点海平面年代际变化也表现出相对应的两个极小(大)值,由此可见太平洋海平面年代际变化受PDO影响显著。比较图5,7可见太平洋中高纬地区的验潮站海平面与PDO的相关性大于低纬度地区,由此可见北太平洋中纬度地区海平面变化受北太平洋年代际震荡的影响更为显著。海平面年代际变化的这种空间分布和时间演变特征与PDO模的时空演变有关。
[25]
升、热带东太平洋下降,1993)2003年的热容海平面趋势在热带太平洋有相似的空间分布。然而,1953)2003年50a的热容海平面变化趋势的空间分布与上面的结论相反,即1953)2003年的50a间,热容海平面变化趋势在热带西太平洋下降,热带东太平洋上升。为探求这种截然相反结果的原因,本文分析了比容海平面和重构场的年代际变化在1993年之后的10a多时间内的变化。
图7为热带东西太平洋重构的海面高度与比容海面高度的年代际变化于1993)2003年间的变化趋势。在热带西太平洋,比容海面高度和重构场的年代际变化在后10a间都表现为上升趋势(上升速率分别为0154和014mm/a),而东太平洋则显著下降(下降速率分别为-0135和-0146mm/a)。正是这种年代际变化对Cazenave[18]15a的海平面趋势的预测产生重大影响,导致了完全相反的海平面变化趋势,从而影响结果的准确性。因此年代际变化对海平面长期趋势项的预测有直接而重大的影响,有可能导致完全相反的趋势项结果。
Cazenave等[18]的研究表明,1993)2008年卫星高度计观测的海平面变化趋势在热带西太平洋上
图7 重构的海平面年代际变化(蓝色实线)于1993)2002年间的趋势项变化和比容
海面高度(红色实线)于1993)2003年间的趋势项变化
a为热带西太平洋,b为热带东太平洋。蓝色虚线为重构海平面后10a趋势,红色虚线为比容海
平面后11a趋势。图中资料经10a低通滤波
1期 顾小丽等:太平洋海平面变化特征及影响因素分析35
6 小结
太平洋地区验潮站海平面总体呈上升趋势,其平均速率为114mm/a,升降速率的大小在空间上分布极不平衡,呈现出西部上升快、东北部主要为下降的区域分布特征,少数站点表现出极大的上升或下降趋势。
太平洋海平面有明显的年际和年代际变化特征。ENSO是影响海平面年际变化的主要因素,二者在西(东)太平洋表现为负(正)相关。ENSO对低纬度地区海平面的影响较大,随着纬度的升高影响有所减弱。海平面年代际变化与PDO指数在西(东)太平洋为负(正)相关,且由低纬到中高纬,二者
的相关性增大。PDO对北太平洋中高纬地区海平面的影响大于热带地区。海平面与PDO的相关性有区域分布特征和随时间演变特征:不同地区二者相关性大小不同,不同时期海平面超前、滞后PDO的关系不同。中国和日本海附近海平面年代际变化可能还受其他因素影响,日本太平洋沿岸海平面年代际变化在35bN以南(北)海平面年代际变化滞后(超前)PDO指数1a左右达最大负(正)相关。由南向北,东太平洋海平面年代际变化超前PDO达到最大正相关的时间由5个月增大到1a左右。年代际变化对目前使用卫星高度计资料分析海平面长期趋势项的预测有直接影响,可能完全掩盖海平面长期变化趋势。
参考文献:
[1] CHURCHJA,GREGORYJM,HUYBRECHTSP,etal1Changesinsealevel[R]MClimateChangetheScientificBasis,Contribution
ofWorkingGroupItotheThirdAssessmentReportoftheIntergovernmentalPanelonClimateChange.Cambridge:CambridgeUnivPress,20011
[2] CAZENAVEA,NEREMRS1Present-daysealevelchange:observationsandcauses[J]1RevGeophys,2004,42;doi:1011029/
2003RG0001391
[3] LOMBARDA,CAZENAVEA,LETRAONPY,etal.Contributionofthermalexpansiontopresent-daysealevelriserevisited[J]1
GlobalPlanetChange,2005,47:1-161
[4] ANTONOVJ,LEVITUSS,BOYERTP1Thermostericsealevelrise:1955-2003[J].GeophysResLett,2005,32;doi:1011029/
2005GL0231121
[5] CABANESC,CAZENAVEA,LEPROVOSTC1Sealevelriseduringpast40yearsdeterminedfromsatelliteandinsituobservations
[J]1Science,2001,294:840)8421
[6] DANIELMW,THOMSONRE1InterannualtomultidecadaltimescaleclimatevariationsinthenortheastPacific[J]1JournalofClimate,
2000,13:3209)32211
[7] WUL,LIUZ,GALLIMORER,etal1Pacificdecadalvariability:thetropicalPacificmodeandtheNorthPacificmode[J]1Journalof
Climate,2003,16:1101)11201
[8] ALBERTMM,ARTHURJM.InterdecadalvariabilityandclimatechangeintheeasterntropicalPacific:Areview[J].Progressin
Oceanography,2006,69:267)2841
[9] 李崇银,朱锦红,孙照渤1年代际气候变化研究[J]1气候与环境研究,2002,7(2):209)2181[10] 莫如平.中国降水、气温的年际振荡及其与厄尔尼诺的联系[J].海洋学报,1989,11(2):16)241
[11] 荣增瑞,刘玉光,陈满春,等1全球和南海海平面变化及其与厄尔尼诺的关系[J]1海洋通报,2008,27(1):1)81[12] 杨修群,朱益民,谢倩,等1太平洋年代际振荡的研究进展[J]1大气科学,2004,28(6):1)51
[13] SMITHTM1TropicalPacificsealevelvariations(1948-1998)[J]1JClimate,2000,13(15):2757)27691
[14] HANW,WEBSTERPJ1Forcingmechanismsofsea-levelinterannualvariabilityintheBayofBengal[J]1JPhysOceanogr,2002,32:
216)2391
[15] YOSHINS,SHOSHIROM,NIKLASS,etal1Decadalsealevelvariabilityinthesouthpacificinaglobaleddy-resolvingoceanmodel
hindcast[J]1AmericanMeteorologicalSociety,2008,38:1731)17461
[16] TIWARIAVM,CABABNESBC,DOMINHBK,etal.CorrelationofinterannualsealevelvariationsintheIndianOceanfromTopex/
Poseidonaltimetry,temperaturedataandtidegaugeswithENSO[J]1GlobalandPlanetaryChange,2004,43:183)1961
[17] 吕俊梅,琚建华,张庆云,等1太平洋海温场两种不同时间尺度气候模态的分析[J]1海洋学报,2005,27(5):30)361
[18] CAZENAVEA,ALIXL,WILLIAML1Present-daysealevelrise:asynthesis[J]1CRGeoscience,2008,doi:1011016/
1crte12008j10710081
[19] LOMBARDA,CAZENAVEA,LETRAONPY,etal1Contributionofthermalexpansiontopresent-daysealevelriserevisited[J]136
GlobalPlanetChange,2005,47:1)161
海洋学报 31卷
[20] HUANGNE,SHENZ,LONGSR,eta1.TheempiricalmodedecompositionandtheHilbertspectrumfornonlinearandnon-stationary
timeseriesanalysis[J].ProcRSocLond,1998,454A:903)995.
[21] WUZ,HUANGNE,LONGSR,etal1Onthetrend,detrending,andvariabilityofnonlinearandnonstationarytimeseries[J].
PNAS,2007,104(38):14889)148941
[22] 申辉,郭佩芳,钱成春11993)2001年全球海平面高度变化特征[J]1海洋与湖沼,2002,33(3):225)2321[23] 马继瑞,田素珍,郑文振,等1太平洋水位站相对海平面升降趋势分析[J]1海洋学报,1996,18(5):14)201
[24] CHURCHJA,WHITENJ,COLEMANR,etal1Estimatesoftheregionaldistributionofsealevelriseoverthe1950to2000period
[J],JClimate,2004,17(13):2609)26251
[25] 咸鹏,李崇银1北太平洋海温变化的年代际模及其演变特征[J]1大气科学,2003,27(5):861-8671
Pacificsealevelvariationsanditsfactors
GUXiao-li,LIPe-iliang
1
1
(11PhysicalOceanographyLaboratory,OceanUniversityofChina,Qingdao266100,China)
Abstract:DecomposingmonthlysealeveldataattidegaugesalongthecoastofthePacificOceanintheEMDmethod,andcombiningwithotherdata,thePacificOceansealevelpsinterannualanddecadaloscilla-tionwereinvestigatel,aswellasinformationonthesealeveltrends1Overall,thePacificOceancoastalsealevelrises,withanaverageratein114mm/a,andshowshighlynon-uniformpropertiesingeography1ThereareobviousinterannualanddecadalvariabilityatthePacificsealevel1TheinterannualvariationsofthePacificOceansealeveloscillationsaresignificantlyaffectedbytheENSO1IntheeastPacificOceantheinterannualsealevelvariabilityandNino3indexarepositivecorrelation,whilenegativecorrelationsinthewestPacificOcean1ThecorrelationbetweenbothissignificantinthetropicalPacificOcean,andthecorre-lationisweakenedwithlatitudeincreasing1Atdifferentareasthiscorrelationcoefficientisdifferent1Inthewest(east)PacificOceansealeveldecadalvariationperformssignificantnegative(positive)correlationwiththePDOindexoverall1Thecorrelationbetweenthemisnon-uniformgeographicallyandevolveswithtime1Sealeveldecadalchangescouldleadtoexactlytheoppositesealeveltrend,duetodecadalvariationdirectimpactonthetrends1
Keywords:sealevelvariations;interannualvariation;decadaloscillation;ENSO
因篇幅问题不能全部显示,请点此查看更多更全内容