班达海(Banda Sea)位于印度尼西亚东部，西靠苏拉威西岛，北、东、南三面为班达弧所环绕，通过海峡与周围的萨武海、帝汶海、阿拉弗拉海、塞兰海、马鲁古海和弗洛勒斯海等连通(图 1)。班达海总面积约47万平方公里，平均水深3064m，最深处位于班达海东部的韦伯海槽约7440m(陈史坚，1989)。印度尼西亚贯穿流(Indonesian Throughflow，ITF)通过印度尼西亚海内众多的海峡从太平洋进入印度洋，是全球海洋热盐环流传送带的关键环节，维持全球大洋间的物质、动量和能量平衡(杜岩等，2011)。班达海作为ITF流经的一个重要海区，接收由ITF运输的来自热带西太平洋的大部分水体和热量，经几个重要的海峡通道流入印度洋。ITF运输总量的约80%通过望加锡海峡(Makassar Strait)，剩余的20%流经东部的哈马黑拉海(Halmahera Sea)或利法马托拉通道(Lifamatola Passage)(Gordon et al，1996)。根据INSTANT(International Nusantara Stratification and Transport)计划的结果，2004—2006年平均ITF入流约为12.7Sv(1Sv=10⁶m³/s)，其中通过望加锡海峡的流入量约为11.6Sv，通过利法马托拉通道的流入量约为1.1Sv; 年平均ITF流出量约为15.0Sv，除了其中约2.6Sv通过龙目海峡(Lombok Strait)流出外，剩余的约12.4Sv通过班达海流入印度洋(Gordon et al，1996，1999，2001，2010)。

观测表明，班达海表层存在明显的季节和年际变化。Gordon等(2001)通过对班达海海表面温度(SST)、海面高度异常(SLA)和局地风场数据的分析，证实了Wyrtki(1961)提出的班达海海表面温度(SST)的季节变化主要受季风影响的机制，东班达海和阿拉弗拉海(Arafura Sea)存在很强的由局地风应力旋度造成的Ekman抽吸，从每年的4月持续到11月，平均上升流量约1.27Sv，其将深层较冷的水输送到表层，从而影响SST。在正常年份下，Ekman上升流导致了班达海较低的海表面温度以及较浅的温跃层，而在ENSO事件发生的年份，这一现象有所改变: EI Niño年东南季风盛行时，Ekman上升流减弱而海表面温度与正常年份相比较冷; La Niña年东南季风盛行时，则对应较强的Ekman上升流以及较暖的海表面温度。Gordon等(2001)认为这是由于受ENSO等大尺度运动的影响，班达海温跃层发生抬升或下降，从而影响海表面温度。Qu等(2008)利用一个高分辨率海洋环流模式证实了印尼海域存在“缓存现象”，具体表现为印尼海上层在每年11月到翌年6月西北季风盛行时表现为“储水期”，储存了更多的水，而在7月至10月东南季风盛行时为“排水期”，他们认为这一现象主要受局地风场导致的Ekman流的辐聚辐散调制。Kida等(2009)发现该海域SST的季节变化大多主要受季风控制，也有一部分受表层热通量的影响。Gordon等(2010)分析INSTANT计划资料时发现，入流海域垂向剖面和出流区域存在明显差异，具体表现为出流区域等温线比入流区域浅15m左右，这说明在印尼贯穿流流经海域存在显著的上升流，并且分为两个部分: 500m以浅的温跃层上升流和600m以深的深层水上升流。Syamsudin等(2010)在对班达海南边界流流速和流量的研究中发现其存在丰富的季节内、季节和年际变化，他们认为这些信号主要受局地风场的季节变化、来自赤道印度洋激发的Kelvin波和ENSO的影响。Yuan等(2013)通过对ITF观测资料的滞后相关分析，证实了Yuan等(2011)提出的ITF是IOD(Indian Ocean Dipole)事件在年际尺度上影响热带太平洋的重要信号通道，在IOD事件发生的年份，ITF自赤道太平洋向印度洋输送更多的暖水，使太平洋温跃层抬升，其激发的上升Kelvin波向东传播，有利于ENSO事件的产生。王健等(2014)发现望加锡海峡中ITF的年际变化主要受与ENSO和DMI(Dipole Mode Index)等大尺度过程相关的Kelvin波、Rossby波调制，局地风场的作用并不是很重要。

本文采用法国空间局AVISO提供的1993年1月—2013年12月的网格化绝对动力地形(MADT: maps of absolute dynamic topography)资料，分析班达海区(图 1黑色框内区域，经纬度范围为3.5°—7.5°S，123.5°—134.5°E)海面高度的变化。第二部分描述数据和方法; 第三部分对MADT资料进行了气候态月平均处理，分析了班达海面高度的季节变化特征; 第四部分计算了MADT年际异常场，并对该异常场进行EOF分解，将得到的第一模态时间序列与Niño3.4和DMI指数进行相关分析，讨论该海区海平面高度的年际变化特征及其与ENSO、IOD事件的关系; 第五部分为总结与讨论。

1 资料与方法

本文采用的卫星测高资料为法国空间局AVSIO下SSALTO/DUACS系统发布的融合了Topex/ Poseidon(T/P)、Janson-1/2、ERS-1/2等多颗测高卫星得到的格点化绝对动力地形资料，时间跨度为1993年1月—2013年12月，时间分辨率为1天，空间分辨率为0.25°×0.25°。绝对动力地形(ADT)为海平面异常(SLA)与平均动力地形(MDT)之和。相比SLA数据，ADT能保留长时间尺度定常流影响的海平面高度信号。风场数据采用欧洲中期天气预报中心ECMWF提供的ERA-Interim数值预报再分析月平均10m风场。时间跨度为1993年1月—2013年12月，空间分辨率为0.75°×0.75°。

MADT数据下载网址: http://www.aviso.altimetry. fr/en/data/products/sea-surface-height-products/global/madt.html;

风场数据下载网址: http://apps.ecmwf.int/datasets/ data/interim-full-moda/levtype=sfc/。

本文采用的海洋200m以浅的月平均温盐场数据来自日本气象厅提供的Ishii V6.13版再分析资料。该资料主要基于WOD05/WOA05(World Ocean Database/ Atlas 2005)、全球温盐剖面计划IRD(L′institut de recherché pour le development，France)采集的温盐资料以及ARGO剖面数据等，时间跨度为1993年1月—2012年12月，空间分辨率为1°×1°。

Ishii温盐场数据下载网址: https://amaterasu.ees.hokudai.ac.jp/~ism/pub/ProjD/。

本文首先计算了班达海MADT气候态月平均空间分布，研究班达海区域海平面的季节变化特征及其空间分布; 随后做出MADT区域月平均的距平场以及对应比容海面高度距平，从整体角度分析班达海面高度季节变化和比容海面高度变化的关系，其中比容海面高度的计算参照Gill(1982)和Shao等(2015):

其中D为水深，ρ₀为标准海水密度(取参照温度T=0，盐度S=35)，ρ为海水密度，ρ₀和ρ的计算均使用基于TEOS-10(Thermodynamic Equation of Seawater 2010)的GSW工具包(下载网址: http://www.teos-10.org/ software.html)。然后对班达海10m风场进行了同样的空间分布处理，并计算了气候态月平均风应力旋度的空间分布，将其与班达海海平面的空间分布进行对比，其中风应力旋度参照Fang等(2006)，Lian等(2015)文章中的计算方法

其中$\text{curl}\vec{\tau }$为风应力旋度，$\vec{\tau }$为风应力，$\nabla =\left( \frac{\partial }{\partial x},\frac{\partial }{\partial y} \right)$为哈密顿算子。风应力计算方法如下

其中为风应力，ρ为大气密度(本文取1.3kg/m³)，为海洋上10m高度风场矢量，U为风速标量，C_d为拖曳系数。对于拖曳系数C_d的取法，本文采用Yelland等(1996)的方法

为研究其年际变化，本文首先对班达海区域逐日MADT资料进行了月平均处理，基于月平均MADT资料计算了其年际异常场，即减去其气候态月平均，得到MADTA(MADT Anomaly)，再对1993年1月—2013年12月的MADTA资料进行了经验模态分解(EOF，empirical orthogonal function)，将其第一模态时间序列与Niño3.4指数、DMI指数作超前滞后相关，分析该海域海平面年际变化与ENSO事件以及IOD事件之间的关系。其中Niño3.4指数定义为太平洋5°S—5°N，120°—170°W区域的海表面温度异常(SSTA)区域平均; DMI指数定义为印度洋10°S— 10°N，50°—70°E区域与0°—10°S，90°—110°E区域的海表面温度异常(SSTA)区域平均之差(Saji et al，1999)。之后对比容海面高度也做了相同的处理，分析班达海海平面年际变化特征与比容海平面变化的关系。

Niño3.4指数与DMI指数下载网址: http://climexp.knmi.nl/selectindex.cgi?id=someone@somewhere

2 海平面的季节变化 2.1 季节变化特征

海平面的变化主要受两个因素影响: 海水质量的变化和海水密度的变化导致的海面高度变化(Gilson et al，1998)，其中后者被称为比容海面高度(Steric height，SH)。Gill等(1973)提出海平面的季节变化主要由大气压强、局地风场和热膨胀造成的比容海平面这三个因素的季节变化所决定。图 2为班达海气候态月平均海面高度距平(MADTA)(减去1993—2013年多年年平均海面高度)和基于Ishii温盐场得到的200m以浅比容海面高度距平(steric height anomaly，SHA)气候态月平均时间序列，从图中可以看出，班达海海面高度和比容海面高度均表现出明显且一致的季节变化特征，其海面高度最大值均出现在3月，最低值出现在8、9月。每年1—4月，班达海海面高度略高于比容海面高度，而每年6—9月，班达海海面高度略低于比容海面高度，前者周年变化幅度约为16.5cm，后者约为14cm，比容海面高度约占海平面高度变化的84%。

图 3为班达海MADTA和SHA空间分布。可见，班达海海面高度3月最高，8月最低，其中班达海东部近阿拉弗拉海区域与西南部近弗洛勒斯海区域海面高度存在着一个“跷跷板”形态。具体表现为: 班达海东部近阿拉弗拉海海面高度周年变化幅度约为27cm，其海平面从每年11月开始升高，到翌年3月达到峰值，随后开始降低，8月达到最低值; 而班达海西南区域近弗洛勒斯海周年变化振幅约为12cm，每年11月海平面开始下降，到翌年2、3月达到最低，随后海平面开始升高，到5、6月达到最高，这与班达海东部3月高、8月低的变化态势相反。SHA的空间分布(图中等值线)与MADTA基本相似，11月到翌年5月较高，6—10月较低，并且在1月、2月、3月、6月、7月、11月也都出现了较为明显的“跷跷板”形态。这一结果符合Gordon等(2001)对班达海SST的研究结果。他们计算了1982—2000年气候态月平均SST以及1993—2000年月平均SLA，发现每年班达海海面高度低值出现在8月底9月初，对应较低的SST以及较浅的温跃层; 最高值出现在4、5月，对应较高的SST; 空间分布上班达海东部近阿拉弗拉海区SST变化幅度最大。然而在图 3中的部分区域，MADTA和SHA表现出了明显差异: 2月、12月的班达海东北部MADTA均明显高于SHA，5月、6月、8月的班达海东北部MADTA明显低于SHA。结合图 2曲线中两者表现出的季节性差异，说明班达海海平面的变化还存在其他较为重要的动力机制。Gordon等(2001)对班达海SST的分析表明，班达海温跃层较浅，且海平面在西北季风盛行时升高，东南季风盛行时下降，这与本地季节变化的Ekman上升流相对应。因此每年1—3月和6—9月两者的差异可能是局地季风影响下，海平面的辐聚辐散导致的。

2.2 季节变化机制

当班达海盛行西北季风时(11月到翌年3月)，弗洛勒斯海大部分表层海水流入班达海后下沉，只有一小部分表层海水通过班达海向北输送至马鲁古海和哈马黑拉海，向南汇入帝汶海; 相反，当东南季风盛行时(5—9月)，班达海和阿拉弗拉海出现上升流来补偿通过帝汶海峡流入印度洋的海水(Wyrtki， 1958)。图 4为班达海各月(1993—2013年平均)海表风场以及风应力旋度分布。从图中可以看出，在5—9月，班达海受东南季风控制，此时班达海东部风场风应力旋度为负，西南部风应力旋度为正; 而12月到翌年3月盛行西北季风，风应力旋度与5—9月相反，即东部北部为正旋度，西部南部为负旋度，而4月和10、11月为季风转换期。结合图 3海面高度距平来看，西北季风盛行时，班达海平均风应力旋度为正，海水辐聚为主，整体海平面较高，其中西南部近弗洛勒斯海的风应力负旋度导致海水辐散，海平面较低，东部近阿拉弗拉海的风应力正旋度导致海水辐聚，海平面较高; 而6月到10月平均风应力旋度为负，整体海平面较低，其中西南部近弗洛勒斯海的风应力正旋度导致海水辐聚，海平面较高，东部近阿拉弗拉海的风应力负旋度导致海水辐散，海平面较低。这一点和图 2中比容海面高度与总体海面高度之间的差异相符合。除了向海洋输入旋度场影响海面高度以及上层热含量外，局地风场激发的Ekman输运也会对班达海海平面造成一定影响。

图 5为局地风场激发的班达海Ekman输运和MADTA与200m以浅SHA之差，如图中所示，在12月到翌年3月，西北季风驱动班达海上层东北向的Ekman输运，使部分海水在班达海东北部的塞兰岛南部和新几内亚岛西南部堆积，从而造成了海水质量上的“南低北高”; 而每年5—9月班达海盛行东南季风激发了班达海上层西南向的Ekman输运，导致海水向西南方向堆积，东北陆地边界附近海水质量减少，造成了质量上的“南高北低”。这一现象也解释了MADTA与SHA空间上的差异，12月和2月东北部的MADTA高于SHA原因是局地风场激发的东北向Ekman输运在东北陆地边界堆积，海水质量增加; 同样5月、6月、8月西南向Ekman输运使东北陆地边界海水质量减少，且班达海东北部水深较浅(图 1)，海水辐散引发的上升流对上层热含量的影响较小，从而导致东北部的MADTA低于SHA。

班达海在1—3月盛行西北季风，驱动上层东北向的Ekman输运，使海水向班达海东北部堆积，并且正的风应力旋度导致班达海上层海水辐聚，总体海面高度略高于比容海面高度，另一方面上层海水堆积引起的上层热含量增加会导致比容海平面的升高; 而6—9月东南季风盛行，驱动上层西南向的Ekman输运，使海水从班达海西南部流出，并且负的风应力旋度导致班达海上层海水辐散，质量减少，总体海面高度略小于比容海面高度，而海水辐散引起的Ekman上升流带来的低温高盐水会使上层热含量减少，比容海面高度较低。

3 班达海海平面高度年际变化

为研究班达海区域海平面的年际变化，本文首先计算了MADT的年际异常场MADTA(MADT Anomaly)，然后对MADTA进行EOF分解。EOF分解后得到前三模态方差贡献率分别为96.10%、1.56%、0.90%，其中第一模态方差贡献率占到了95%以上，而后几个模态所占方差很小，因此本文着重对第一模态进行分析讨论。图 6为EOF分解第一模态(EOF1)的空间分布，从第一模态空间分布中可以看出，班达海区整体海平面变化表现出同相位，其中班达海东部海区海平面变化幅度略高于西南海区近弗洛勒斯海的变化幅度。

图 7为第一模态时间序列PC1(负)与Niño3.4指数和DMI指数，其中Niño3.4指数和DMI指数分别代表太平洋ENSO事件和印度洋IOD事件。据此分析班达海面高度的年际变化与ENSO和IOD之间的关系。从图 7中可以看出，PC1与Niño3.4指数有很强的相关性，尤其在1997—1998、2002—2003、2004—2005、2006—2007和2009—2010的El Niño年以及1995—1996、1999—2000、2011—2012、2007— 2008和2010—2011的La Niña年中均表现出了较为明显的相位变化，在正负IOD事件的年份，PC1也表现出一定相关性，具体表现为EI Niño年和正IOD事件对应较高的海平面，La Niña年和负IOD事件对应较低的海平面。

图 8为PC1与Niño3.4指数、DMI指数的超前滞后相关。从图中可以看出，PC1与Niño3.4指数和DMI指数均呈负相关，其中PC1滞后Niño3.4指数1个月的相关系数为0.76并超过95%置信度水平。班达海海平面年际尺度上的变化与ENSO间的关系，许多前人的工作中也有总结(Wyrtki，1961; Bray et al，1996; Gordon et al，2001; Iskandar，2010; Syamsudin et al，2010)。其中Gordon等(2001)用1993年1月—1999年12月的SLA数据与Niño3指数对比，证实了Bray等(1996)提出的班达海海平面La Niña年较高而EI Niño年较低的变化特征。本文利用更长时间资料，发现2000—2013年班达海平面与ENSO间依然显著相关，尤其在El Niño年与La Niña年中，海平面的起伏有明显的相位变化。

班达海SST的年际变化与ENSO和IOD事件有密切的联系，具体表现为La Niña年对应较高的SST，EI Niño年和IOD事件对应较低的SST(Gordon et al，2001; Wijffels et al，2004; McClean et al，2005; Iskandar，2010)。而SST的变化与海洋上层热含量的变化密不可分，图 9为减去气候态月平均的MADT和SH的区域平均时间序列，与MADT相似，比容海面高度也表现出明显的年际变化，并且变化趋势与MADT基本一致，在El Niño年比容海面降低，La Niña年升高，两者同期相关系数为0.93，并通过99%置信度检验; MADT年际变化幅度约为36cm，比容海平面年际变化幅度约为17cm，两者在极端ENSO事件发生的月份存在显著的高度差，这说明ENSO等大尺度运动发生时，不仅伴随着由于班达海上层热含量的变化导致的海平面高度变化，还存在着明显局地海水质量的变化。

El Niño建立初期，赤道西太平洋西风异常，激发暖的下沉Kelvin波和冷的上升Rossby波，其中下沉Kelvin波向东传播，使赤道东太平洋海区SST上升，促进El Niño的形成; 而上升Rossby波西传至印度尼西亚群岛，其中一部分以上升沿岸陷波(coastally trapped waves，CTW)形式传入班达海，使班达海温跃层抬升，SST降低; La Niña年则情况相反(Wajsowicz et al，2003; Iskandar，2010)。因此El Niño发生的年份班达海比容海面较低，La Niña年则对应较高的比容海面。ENSO事件发生时，通常伴有大尺度的天气过程，例如El Niño现象发生时，沃克环流减弱，西太平洋大气上升区东移，班达海表面出现西风异常; La Niña现象发生时，沃克环流增强，西太平洋大气上升区向西移动，班达海表面出现东风异常。这一大尺度天气过程势必造成班达海海水的堆积和扩散，从而影响班达海海平面高度。

4 总结与讨论

本文利用1993年1月—2013年12月的卫星测高数据，分析了班达海海平面的季节和年际变化特征。在季节尺度上，班达海海平面12月—翌年4月较高，7—10月较低，最低出现在8月中下旬，其中区域平均海平面高度3月中上旬最高，随后一直下降，到8月中下旬达到最低，全年变化幅度为16.5cm左右。班达海海平面季节变化具有明显的空间差异，东部近阿拉弗拉海全年海平面变化幅度较大，12月—翌年4月较高，7—10月较低; 而西部近弗洛勒斯海海平面变化幅度较小，12月—翌年4月较低，7—10月较高，两者全年表现出偶极子形态。海平面的变化受比容海平面变化和海水质量变化共同调制。其中班达海比容海平面变化与总体海平面变化相关密切，全年变化幅度为14cm 左右，占海平面变化的82%。受季风影响，班达海上层的辐聚(辐散)是导致海水质量变化的主要原因: 每年12月—翌年3月班达海盛行西北季风，正的风应力旋度造成了上层海洋的辐聚，温跃层下降，海平面上升; 5—9月东南季风盛行，负的风应力旋度造成了上层海洋的辐散，温跃层上升，海平面下降。同时激发的下降(上升)Ekman流，会影响海洋上层热含量的变化，从而影响比容海平面。而局地风场激发的水平方向Ekman输运会造成海水质量的季节性南北空间分布差异。因此班达海海面高度的季节变化主要由上层热含量变化导致的比容海平面变化和局地风场的季节变化共同调制，其中局地风场不仅通过海水质量的堆积与分散影响海平面高度变化，还能通过引起的上升、下降流影响海洋上层热含量的分布，从而调控比容海平面的变化。

年际尺度上，班达海海平面呈典型的ENSO形态，其中MADTA的EOF分解后第一模态方差贡献率高达96%，第一模态时间序列(PC1)与Niño3.4指数和DMI指数均呈负相关，其中PC1滞后Niño3.4指数1个月的相关系数为0.76，超过95%置信度水平。这说明班达海海平面的年际变化与大尺度的ENSO事件有密切的关系。在El Niño、La Niña的年份，PC1均表现出明显的相位变化。班达海比容海平面也存在明显的年际变化特征，其年际异常序列与MADT年际异常序列有较强的相关性，同期相关系数为0.93。MADT年际变化幅度约为36cm，比容海平面变化幅度约为17cm，两者在极端ENSO事件发生的月份存在明显的高度差，说明ENSO对班达海海平面的影响不仅体现在对其上层热盐含量的影响上，还体现在海水质量的变化上。

本文工作进行时，曾尝试使用重力卫星GRACE资料，定量分析班达海海水质量的变化。然而班达海所处的印尼海域岛屿众多，地形复杂，空间分辨率1°×1°的GRACE资料在班达海的准确性有待商榷。本文下一步工作将利用ISTANT计划在ITF主要入流和出流通道中观测得到的流量资料，结合班达海海平面的变化，分析两者间的关系。并且通过班达海入流、出流的流量差，推算海水质量变化对班达海海平面变化的影响。