印度尼西亚海及周边海域层结的时空变化特征分析

梁浩然，谢玲玲,3*，周磊，郑全安，李明明,3

( 1. 广东海洋大学 海洋与气象学院 近海海洋变化与灾害预警实验室，广东 湛江 524088；2. 广东省高等学校陆架及深远海气候资源与环境重点实验室，广东 湛江 524088；3. 南方海洋科学与工程广东省实验室（湛江），广东 湛江 524025；4. 上海交通大学 海洋研究院，上海 200240；5. 南方海洋科学与工程广东省实验室（珠海），广东 珠海 519080；6. 美国马里兰大学 大气与海洋科学系，马里兰州 大学市 20742)

1 引言

印度尼西亚海(Indonesian seas，简称印尼海)横跨赤道, 位于亚-澳大陆之间，具有复杂的岸线、狭窄的水道、众多的岛屿以及崎岖的海底地形[1]，是太平洋和印度洋在低纬度海域进行物质和能量交换的主要通道[2]和全球热盐输送带的关键环节[3-4]。该海域还是太平洋-印度洋大尺度海洋波动的波导汇聚区[5-6]和大气季节内振荡活跃区，在海盆季风和气候模态变化、全球物质能量平衡，以及气候系统演变中具有重要作用[3,7-8]。同时，作为连接南海与印度洋的重要通道，印尼海还是我国“海上丝绸之路”的要地。认识印尼海水文特性和动力过程，是开展海上丝绸之路等海洋活动的重要基础。

海水层结和密度跃层是海洋基本水文特征，是海洋环流、内波等动力过程、海洋生物地球化学过程以及气候系统变化等科学研究的重要参数，直接影响着海洋能量循环和物质垂向输送[9-10]。海洋层化在海洋军事、航海输运以及渔业养殖[11]等海洋应用方面也有着重要作用。因此，开展印尼海海水层结特征研究，具有重要的科学意义和实用价值。

目前，印尼海海水层结时空变化还鲜有研究，前人多关注了该海域印尼贯穿流（ITF）和潮汐潮流特征。研究显示，ITF 共分为西、中、东3 条路径，西部路径经卡里马塔海峡流入，水体输运量约为2×106～4×106m3/s[12-15]；中部路径通过望加锡海峡流入，输运量为9×106～11.6×106m3/s，是ITF 的主要路径[16-18]；东部路径通过马鲁古海和哈马黑拉海流入，估计输运量为1×106～3×106m3/s[19-20]。ITF 具有显著的多尺度变化特征[21]，在季节内尺度上，ITF 受大气季节内震荡（MJO）影响明显[22]；季节尺度上，ITF 具有明显的年周期和半年周期变化，其最大通量出现在8-9 月[23]；年际上ITF 主要受厄尔尼诺-南方涛动（ENSO）和IOD事件影响，并反作用于ENSO 事件[24-26]。印尼海域内潮汐、潮流复杂多变，正压潮与内潮之间、不同海域内潮之间以及不同周期内潮分潮之间均会发生相互作用[27-28]。该海域正压潮以半日潮为主[29]，海峡处M2分潮振幅可达2 m[30]；并且半日分潮和全日分潮的传播方向存在明显差别，半日分潮由东印度洋传入，全日分潮由西北太平洋传入[27,29]；大部分区域内潮流速大于正压潮，在苏拉威西海内潮最大流速可达0.2 m/s[31]。前人还对该海域的水团来源进行了分析，Nof[32]、Gordon 和Fine[33]证实盐度较低的北太平洋水通过望加锡海峡温跃层进入印尼海，盐度较高的南太平洋水通过印尼海东部温跃层下层进行输运。张艳慧等[34]发现印尼海东部中层水体可能是由南太平洋南极中层水（AAIW）和南太平洋热带水（SPTW）垂向混合形成的，也可能是印度洋的AAIW 向北延伸的结果。Wang等[35]发现，马鲁古海次表层水团来自于北太平洋热带水（NPTW），中层受来自于南太平洋的AAIW 控制；哈马黑拉海次表层水体来自于SPTW，中层来自AAIW；水团来源存在季节和年际变化。张晶等[7]研究表明从卡里马塔海峡输入的高温低盐的南海表层水会对印尼海产生重要影响。

在南海、西北太平洋等印尼海周边海域，部分研究分析了浮力频率和跃层的变化。在南海研究方面，刘桂芳[36]简要分析了粤东陆架区夏季密度跃层强度和厚度的分布特征，得出在远海密度跃层强度逐渐减弱、厚度逐渐增加的规律；张绪东等[37]和刘金芳等[38]详细阐述了台湾周围海域密度跃层的空间分布和季节变化规律；廖光洪等[39]发现南海北部浮力频率存在显著的季节变化特征，而海盆深水区浮力频率最大值季节性变化较弱，但空间分布存在明显的季节差异；郑曼立等[40]得出粤西和琼东海区冬季近岸区域混合层较深，外海密度跃层深度位于60～120 m 且层结较强，浮力频率大于10-2s-2；袁瑞等[41]发现在南海范围内，拉尼娜年海洋层化现象较弱，厄尔尼诺年层结较强。在西北太平洋研究方面，Emery 等[42]发现北太平洋浮力频率分布存在明显的地域差异和季节性差异；蒲书箴等[43]发现热带西太平洋浮力频率具有明显的年际变化，与ENSO 有密切关系；江伟等[44]利用模式同化的温盐数据产品对西北太平洋密度跃层的类型、上下界以及强度的空间分布特征和季节变化进行分析，得出西北太平洋不同纬度的密度跃层在季节上存在明显差异。近期，Yamaguchi 和Suga[45]发现随着全球变暖，自1960 年以来全球约40%的海洋密度层结正在加强，垂向物质输运变弱，其中热带区域更加显著。

印尼海海域上层温度受局地太阳辐射以及赤道纬向风引起的自由开尔文波和罗斯贝波影响[6]，同时ITF 各路径将南海表层大量淡水和热量以及南、北太平洋次表层高温高盐水和中层低盐水输入印尼海[7,32-35,46]，经蒸发、降雨以及潮汐混合等多种因素的作用，使印尼海形成了独特的层结和密度跃层分布特征。但在前人研究中，针对印尼海密度跃层和层结强度三维空间特征的研究还较少，尤其对ITF 不同路径上的差异变化认识较少。此外，印尼海区跃层和层结强度的季节和年际变化规律也有待进一步研究[47]。因此，本文利用气候态WOA13（World Ocean Atlas 2013 Version 3.1.1）数据和SODA v3.3.1（Simple Ocean Data Assimilation）数据，研究印尼海及周围海域（20°～25°N，95°～141°E），尤其是ITF 3 条路径上的层结强度和密度跃层的空间特征和时间变化规律。

2 数据与方法

2.1 数据来源

本文所用数据取自美国国家海洋和大气管理局（NOAA）国家海洋数据中心的WOA13 数据集产品[48-49]和美国马里兰大学的SODA v3.3.1 海洋数据集[50-51]，下载网址分别为https://www.nodc.noaa.gov/OC5/woa13/woa13data.html 和https://www.atmos.umd.edu/～ocean/index_files/soda3.3.1_mn_download.htm。WOA13 为全球海洋气候学数据，包含1955-2012 年气候态月平均、季平均和年平均的温度、盐度、溶解氧等多种海洋要素，水平空间分辨率有5°、1°和0.25° 3 种，垂向从表层到最大深度层5 500 m 共分102 层。SODA v3.3.1为海洋资料同化数据，包含1993-2015 年全球月平均温度、盐度、流场、风场和海面高度等多种全球海洋要素，空间分辨率为0.5°×0.5°，垂向方向共50 层，深度范围为5～5 390 m。

选取WOA13 数据集中空间分辨率为0.25°×0.25°的气候态年平均和月平均温盐数据，数据纬度范围为20°～25°N，经度范围为95°～141°E，共33 485 个网格点，垂向上取上层500 m、共37 层。选取SODA v3.3.1数据集1993-2015 年的月平均温盐数据，数据空间分辨率为0.5°×0.5°，纬度范围为20°～25°N，经度范围为95°～141°E，共8 463 个网格点，垂向上取上层500 m、共24 层。

2.2 ITF 路径选取

在印尼海众多海峡中，望加锡海峡是ITF 的主要流入通道，其次是利法马托拉海峡和卡里马塔海峡[33]；ITF 3 个主要流出通道分别是龙目海峡、翁拜海峡和帝汶海峡[21]，在INSTANT 观测中这3 个海峡的平均流量分别为2.6×106m3/s、4.9×106m3/s 和7.5×106m3/s[52]。由于卡里马塔海峡水深较浅，早期对南海海水是否通过卡里马塔海峡进入印尼海存在争议[53]。后续观测和数值模拟显示卡里马塔海峡是南海南部海区的主要流出通道，流量达1.6×106～3.6×106m3/s[54-58]。Qu等[12,59-60]给出了太平洋-印度洋贯穿流南海分支示意图，指出北赤道流在菲律宾海岸分叉出向北运动的黑潮，到达吕宋海峡附近后，一支分流穿过吕宋海峡进入南海腹地，途径西沙群岛、卡里马塔海峡、爪哇海、龙目海峡最后进入印度洋。由于卡里马塔海峡水深小于50 m，通过这一海峡主要将南海表层的高温低盐水输入印尼海域，因此本文选取这一路径上表层（25 m 处）气候态流场（图2a）流速最大值为ITF 西部路径断面（图1 黄色实线）。参照Gordon[61]和Wang等[35]给出的印尼贯穿流分支示意图，在中部路径，次表层的NPTW 通过棉兰老流，向南进入苏拉威西海，途径望加锡海峡和弗洛里斯海，从帝汶海峡进入印度洋。东部路径主要有两条通道，一是次表层NPTW通过棉兰老流分支进入马鲁古海，后进入班达海；另一条路径是次表层SPTW经新几内亚沿岸流进入哈马黑拉海和塞兰海，途径班达海通过帝汶海峡进入印度洋；本文主要研究SPTW进入哈马黑拉海的路径断面。为进一步确定中部和东部路径，选取中部和东部路径上150 m 深度的气候态流场（图2b）流速最大值为ITF 中部路径断面（图1 红色实线）和ITF 东部路径断面（图1 粉色实线）

图1 印尼海及周边海域地形图Fig. 1 Map of Indonesia seas and surrounding waters

图2 印尼海及周边海域气候态流场分布Fig. 2 2D distribution (x-y) of climatological velocities in the Indonesian seas and surrounding waters

2.3 数据处理

以图1 中A-F 点为例，对比WOA13 气候态年平均数据和SODA 逐月平均数据插值前后的密度和层结剖面，可以看出各海域插值前后上500 m 的密度剖面和N2剖面几乎完全重合（图未显示），因插值引起误差可忽略。同时，对比WOA13 和SODA 数据，二者剖面几乎完全相同，说明两种数据具有统一性。

本文采用垂向梯度法进行密度跃层判定[44,63]，在水深小于200 m 的近岸海域，跃层密度梯度大于0.1 kg/m4；在水深大于200 m 的海域，跃层密度梯度大于0.015 kg/m4。

3 气候态分布特征

3.1 沿ITF 路径断面分布

图3 给出ITF 3 条路径500 m 以上温盐断面分布。在西部路径（图3a、图3d），太平洋表层水温度较高，可达29℃，温度梯度较小，盐度较低；次表层为高温高盐的NPTW[64]，盐度最高可达35，温度梯度较大。南海表层温度较高，盐度较低；次表层南海热带水（SCSTW）温度略低于NPTW，等温线与太平洋相比有明显抬升，温度梯度最大值约在50～80 m 深度处。通过卡里马塔海峡，西部路径将南海表层相对高温低盐的海水输入印尼海域[7]，受赤道强太阳辐射、降雨、地表径流以及较强的潮致混合影响，使卡里马塔海峡和爪哇海整体呈现高温低盐的特征。印度洋海域表层为高温低盐的印度洋东部表层水（EIW），在次表层由东向西海水从低温低盐的亚-澳中间水（AAMW）逐渐转变为低温盐度较高的北印度洋中央水（NICW），随深度增加在100 m 附近海水温度梯度最大。

图3 ITF 路径上气候态温度（a-c）和盐度（d-f）断面分布Fig. 3 Distribution of climatological temperature（a-c）and salinity (d-f) along ITF pathways

中部路径（图3b、图3e），太平洋海域温盐的垂向分布特征与西部路径基本相同。在次表层，高盐的NPTW 通过棉兰老流进入印尼海域后与中层低盐水混合，盐度逐渐降低，特别是望加锡海峡处海水受地形影响产生较强的垂向混合作用，使这一区域次表层海水盐度明显低于太平洋；弗洛里斯海和望加锡海峡表层低盐水主要来自于通过卡里马塔海峡输入的南海水；最后混合较为均匀的低盐海水通过帝汶海峡流入印度洋，海水盐度逐渐升高。沿中部路径，表层海水的温度呈现逐渐降低的变化趋势，次表层温度等值线呈现抬升趋势，约在150 m 处，温度梯度最大。

在东部路径（图3c、图3f），太平洋区域表层为高温低盐的西部赤道-热带水（WET），最大温度可达30℃；次表层为高盐的SPTW，盐度最大值可达35.5，在150 m附近温度梯度最大。在次表层，高盐的SPTW向哈马黑拉海流入并与中层低盐水相互混合，盐度逐渐降低。塞兰海和班达海的温盐垂向分布特征基本相同，上层为高温低盐水，随深度增加温度逐渐降低，盐度逐渐增大，温度梯度在110 m 附近达到最大值，总体来看沿东部路径温度等值线呈现逐渐抬升的趋势。前人分析显示，印尼海中东部海域次表层主要由SPTW 和NPTW 控制，中层主要受北太平洋中层水（NPIW）和AAIW 的影响[33-35]；西部海域主要受南海表层水和SCSTW 影响[7]。

图4a 表示ITF 3 条路径上500 m 浮力频率断面分布。在西部路径（图4a）海水分层显著，其中太平洋海域混合层厚度约为50 m，跃层较深，分布在50～250 m，N2最大值为10-3.6s-2；从吕宋海峡进入南海，混合层厚度和跃层深度沿路径逐渐变浅，层结强度逐渐增加，在卡里马塔海峡以北N2最大值达到10-3.2s-2，混合层厚度小于20 m；爪哇海海水水深较浅，N2垂向分布较为均匀，强度最弱，N2多小于10-3.9s-2；进入东印度洋，混合层厚度增大到30 m，跃层约分布在30～170 m 深度层，N2达到10-3.6s-2。整个路径总体来看，南海南部海域的层结强度最大，卡里马塔海峡和爪哇海海域层结强度最小，太平洋和印度洋介于二者之间。

沿中部路径（图4b），太平洋上混合层厚度约为40 m，跃层分布深度约为40～220 m，在120 m 左右浮力频率最大，N2最大值约为10-3.4s-2；进入苏拉威西海，海水混合层变浅，跃层抬升至10～200 m 深度，在望加锡海峡和弗洛里斯海，上混合层厚度几乎为0；由帝汶海峡进入印度洋，上混合层厚度约为30 m，跃层分布深度约为30～220 m。与印尼海相比，太平洋和印度洋的上混合层厚度以及跃层深度明显较深。

沿东部路径（图4c），混合层厚度和跃层深度分布较为均匀，混合层约为20 m，跃层分布深度约为20～220 m，在100 m 上下深度层N2达到最大值约为10-3.4s-2，在太平洋和哈马黑拉海北部部分海域可以看出明显的双跃层结构。

对比图3 和图4 可以看出，在ITF 3 条路径上，海水表层和次表层的层结强度主要由海水温度控制，温度梯度最大值所在深度与层结强度最大值所在深度几乎重合，在表层和深层海水温度梯度较小的区域，海水层结强度较弱。而盐度对海水层结强度的贡献较小，东部和西部路径太平洋海域次表层盐度梯度很大，但层结强度却略小于其他区域。

图4 ITF 路径上气候态浮力频率（N2）断面分布（N2 的单位：s-2）Fig. 4 Distribution of climatological squared buoyancy frequency (N2) along ITF pathways (unit of N2 is s-2)

表1 列出气候态下ITF 3 条路径上跃层平均浮力频率、平均深度和厚度结果。可见，各路径跃层内平均N2差异较小，其中中部路径平均最大，为10-3.68s-2，东部路径最小，为10-3.71s-2；各路径上N2变化不大，其中西部路径偏差最大，为10-0.13s-2。各路径跃层深度和厚度平均值存在明显差异，东部路径跃层深度和厚度最大，分别为124 m 和192 m，中部次之，西部最小为99 m 和143 m。不同海域对比来看，印尼海域3 条路径的跃层深度平均值均要小于其他海域，跃层平均厚度分布特征与深度类似。

表1 ITF 路径平均跃层强度、深度和厚度Table 1 Mean depth and thickness of pycnocline along ITF pathways

表1 ITF 路径平均跃层强度、深度和厚度Table 1 Mean depth and thickness of pycnocline along ITF pathways

注：-表示无数据。

跃层平均强度log10N2跃层深度Hpyc/m跃层厚度Dpyc/m西部中部东部西部中部东部西部中部东部太平洋-3.78-3.68-3.72142139131180185202南海-3.64--79--126--印尼海区-3.86-3.66-3.714510011857179185印度洋-3.64-3.70-112123-127191-平均±偏差-3.70±0.13-3.68±0.04-3.71±0.0299±55120±27124±9143±38185±22192±18

3.2 平面分布特征

图5a 为垂向最大浮力频率（N2max）的平面分布。可见，在整个研究区，印度洋东侧苏门答腊岛附近海域N2max最高，约为10-2.9s-2，南海海盆东南部次之，达到10-3.1s-2，卡里马塔海峡和爪哇海N2max最小，为10-4.1s-2，随着远离海岸呈现逐渐减小趋势。ITF 中部路径途经海域N2max较高，其中热带太平洋、苏拉威西海北部和弗洛勒斯海西部N2max达到10-3.3s-2；东部路径N2max分布较为均匀，多在10-3.4s-2左右。总体而言，印尼海中部和南部海域N2max分布较为均匀。太平洋除棉兰老岛西部，其他区域N2max均小于10-3.5s-2，尤其是北赤道流和黑潮影响区域仅为10-3.7s-2。受太平洋影响，南海东北部陆架N2max较低为10-3.5s-2，海盆中部以及南部N2max逐渐增大，整体呈现从北向南逐渐增加趋势。东印度洋N2max整体呈现出随纬度增加的变化趋势。

图5 印尼海及周围海域气候态层结特征平面分布Fig. 5 2D distribution (x-y) of climatological stratification characteristics in the Indonesian seas and surrounding waters

结合图4、图5 以及表1 可以看出，从太平洋到印度洋ITF 3 条路径层结均存在由强变弱再变强的趋势，跃层深度和混合层深度则存在由深变浅再变深的趋势。由图3 可以看出，3 条路径从太平洋进入印尼海域海水温度等值线均存在较为明显的抬升现象，中部和东部路径次表层高盐水（SPTW 和NPTW）的盐度迅速降低，在苏拉威西海和哈马黑拉海海水盐度分别从35 和35.5 降至34.6 和34.7 左右，在垂向上次表层高盐水与中层的低盐水逐渐趋于均匀，水团性质发生显著变化；西部路径进入卡里马塔海峡后海水的盐度也有显著降低。已有研究表明，潮致混合是ITF 区域水团转化的重要成因[61,65-66]。印尼海域正压潮向斜压潮的能量转移十分剧烈，平均能量转移量是全球海洋转移量的20 倍[67]。同时，印尼海是一个半封闭海域，内潮生成后仍然被局限在这个区域内，几乎所有的内潮能量在印尼海域内被耗散[68]。该海域平均的垂直扩散率可达到1.5 cm2/s，在塞兰海、哈马黑拉海和望加锡海峡等区域扩散率最大值可达几十个cm2/s[65,69]。印尼海域内较强且分布不均匀的内潮混合对海域内的海表面温度（SST）分布特征也有较大影响[70-71]。耦合模型表明，如果没有内潮混合，海域内SST 会增加0.5℃[66]。由于较强的潮汐混合，使ITF 路径上印尼海域内海水的层结强度低于其他海域。进入印度洋后，沿西部和中部路径海水和N2均逐渐增大，并且从图4 断面中也可以看出跃层内N2有逐渐增大的变化趋势，其主要原因是ITF 将印尼海域混合较为均匀的低温低盐海水带入印度洋，并逐渐与印度洋海水进行混合同化，从而形成了随着远离ITF 流出海峡海水层结强度逐渐增大的变化过程。

4 季节分布特征

研究海区位于低纬度热带海域，受亚澳季风影响，冬季11 月至翌年4 月主要为西北风，夏季5-10 月主要为东南风[35]。选取W0A13 数据集中2 月和8 月气候态月平均数据分析研究海区层结季节变化。

图6 示出ITF 3 条路径夏、冬两季浮力频率断面分布。可见，沿西部路径（图6a 和图6b），太平洋夏季混合层厚度明显小于冬季，约为冬季一半，跃层深度夏冬季基本相同，太平洋夏季出现双跃层，跃层强度相当；南海夏季的层结强度则明显高于冬季约1 个量级，混合层深度和跃层深度明显小于冬季；卡里马塔海峡以北，夏季底层出现强层化，N2最大值达10-3.2s-2，而冬季则混合均匀，跃层现象不明显，而卡里马塔海峡南部，冬季层结强度大于夏季。爪哇海夏季层结较弱，冬季层结较强；印度洋夏季混合层厚度和跃层深度要明显大于冬季，跃层厚度和层结强度小于冬季。跃层深度夏冬季差异明显，夏季由太平洋到南海显著抬升，冬季太平洋和南海变化不大，而印度洋海域出现跃层抬升现象。夏、冬季层结差ΔN2如图6c 所示，可见太平洋和南海100 m 以上ΔN2为正值，夏冬季最大相差0.8 个量级，印尼和印度洋海域ΔN2负值，N2最大相差1 个量级。

图6 ITF 路径上夏、冬季浮力频率（N2）断面Fig. 6 Distribution of summer and winter squared buoyancy frequencies (N2) along ITF pathways

沿中部路径（图6d 和图6e），太平洋和印度洋夏冬两季跃层变化与西部路径跃层变化特征类似。苏拉威西海和弗洛里斯海夏季混合层可达50 m，跃层集中在50～200 m，浮力频率最大值出现在苏拉威西海，达到10-3s-2；冬季，印尼海混合层几乎消失，跃层分布在20～180 m 深度，平均层结强度大于夏季，N2最大值出现在弗洛里斯海，约为10-3s-2，与夏季相当。在100 m 以上深度层，太平洋海域层结强度大于冬季，最大差值达到1.5 个量级（图6f），印尼海和印度洋夏季层结强度小于冬季，最大差值可达到两个量级，100～500 m 深度之间的大部分海区夏季浮力频率要小于冬季。

东部路径（图6g 和图6h），太平洋夏季混合层厚度可达100 m，跃层集中在100～200 m，冬季混合层厚度约50 m，在50～200 m 处层结强度较大，夏季跃层厚度小于冬季，但平均层结强度大于冬季。哈马黑拉海和塞兰海夏季混合层厚度和跃层厚度略大于冬季，跃层内平均层结强度小与冬季。班达海夏季混合层厚度和跃层深度较大，跃层厚度和跃层内平均层结强度小于冬季。从夏冬两季ΔN2的断面分布来看，太平洋和班达海100 m 以上深度夏季浮力频率小于冬季，而100～200 m 深度层夏季浮力频率略大于冬季；哈马黑拉海和塞兰海100 m 以上深度层夏冬季N2基本相同，100～200 m 水深夏季小于冬季（图6i）。

印度尼西亚海层结强度的季节性差异由多种因素共同作用产生。该海域冬季盛行西北季风，夏季盛行东南季风，印尼大部分海域夏季风强度大于冬季风强度，因此由风应力引起的海水上层混合夏季要大于冬季。在西部路径，冬季南海大量的高温低盐水进入印尼海域[7]，使爪哇海处层结强度明显增大（图6b）；夏季受东南季风强迫，卡里马塔海峡上层海水向北流动[58,72-73]，混合均匀的爪哇海海水流回南海，使南海南部海域上层混合层厚度和底层的层结强度均明显增加（图6a）。在中部路径，望加锡海峡处夏季海水南向流速大于冬季，使该海域夏季混合强度大于冬季[2]；同时，冬季南海高温低盐水经过卡里马塔海峡输入望加锡海峡，使该海域层结强度增强（图6e）。此外，印尼海域的季节变化还受到太平洋斜压罗斯贝（Rossby）波和印度洋开尔文（Kelvin）波的影响[74-75]。在印度洋季风转换期间激发的开尔文波沿赤道向东传播至苏门答腊沿岸，然后沿岸界向东传播，最远可穿过龙目海峡抵达望加锡海峡[52]。Susanto 等[18]研究表明，5 月和11 月印度洋开尔文波对望加锡海峡南向海流速度影响最大。开尔文波抵达望加锡海峡后会减弱该海峡处海水的南向流动，进而影响该海域的层结强度。

5 年际变化分析

图8 1993 年1 月至2015 年12 月沿ITF 路径断面最大浮力频率（N2max）逐月变化Fig. 8 Distance-time distribution of monthly maximum squared buoyancy frequency along ITF pathways from January 1993 to December 2015

关注N2max的时间变化，可见明显年周期信号，这与印尼海域的季节性变化有关。同时，N2max变化在不同年份也有显著不同：1997 年5 月至1998 年5 月强El Niño 期间，ITF 各路径的最大浮力频率均明显大于其他年份；1998 年7 月至2001 年3 月强La Niña 期间，ITF 各路径N2max均明显小于其他年份。由此推测研究区层结年际变化可能与ENSO 事件有关。

为进一步研究该海域层结特征逐年变化趋势以及周期变化，在ITF 3 条路径上分别选取3 个特征较为明显的站点，按西部、中部、东部路线顺序分别为：A 点（18°N，112°E）、B 点（5°N，107°E）、C 点（10°S，104°E）、D 点（11°N，132°E）、E 点（4°S，120°E）、F 点（12°S，120°E）、G 点（0°N，133°E）、H 点（2.5°S，128°E）、I 点（6°S，139°E），如图1 红色“*”所示位置。

图9 示出A-I 点1993-2015 年N2max的变化趋势（纵轴（y轴）为N2max距平值），其中A、I 点的长期趋势通过95%置信度检验，C、D、E、F、H 点通过99%置信度检验。由图可见，在选取的9 个点中C（印度洋西部）、D（太平洋中部）两点在23 a 中N2max的变化趋势是减小的，其余7 点在23 a 中N2max均是增加趋势，其中F（印度洋中部）、H（哈马黑拉海）两点的增大趋势最为明显，在23 a 中N2max增加了近0.1 个量级。总体而言，在1993-2015 年的23 a 中，研究海区大多数海域的层结强度是逐渐增强的，各海域层结增强的原因主要与气候变暖引起的SST 增加有关。对于层结减弱的C 和D 点，进一步分析发现C 点附近海表面风速在23 a 中平均增长了约0.25 m/s，N2max所在深度的海水温度增长了0.67℃，D 点N2max所在深度的海水温度增长了2℃（图未示出）。由此我们推测D 点层结强度减弱主要是由下层海水升温引起的；C 点层结强度减弱是由上层海水混合增强以及下层海水升温共同作用产生。

图9 图1 中A-I 点1993-2015 年N2max 变化趋势Fig. 9 Variation of the maximum squared buoyancy frequency (N2max) at points A-I marked in Fig. 1 from 1993 to 2015

图10 和图11 分别示出A-I 点最大浮力频率和跃层平均深度的功率谱。由图10 可见，除A 点（南海）N2max多年周期较弱外，其余各点均有较为明显的4～7 a 以及10 a 以上的多年周期变化。D（太平洋中部）、H（哈马黑拉海）、I（班达海）点N2max的年周期信号较弱，其余各点处N2max均有显著的1 a（12 个月）周期变化。从图11 中可以看出，A、C（印度洋）两点Dpcline的多年周期较弱，年周期较强，D、I（班达海）两点Dpcline的年周期变化较弱，多年周期信号较强，其余各点的Dpcline均有较为显著的1 a 周期以及4～7 a 和10 a 以上周期变化。在这些周期中，1 a 周期信号主要与研究海域的季节变化特征有关，4～7 a的周期信号，与前人提出ITF 受ENSO 影响呈现的主周期相符，10 a周期可能与太平洋年代际振荡（PDO）有关。关于ENSO和PDO 对印尼海层结的调制作用有待更多研究。

图10 图1 中A-I 点最大浮力频率功率谱Fig. 10 Power spectra of maximum squared buoyancy frequency at points A-I marked in Fig. 1

图11 图1 中A-I 点跃层平均深度功率谱Fig. 11 Power spectra of the mean depth of pycnocline at points A-I marked in Fig. 1

6 结论

层结强度和跃层厚度是海洋水文动力和水声传播等领域的重要参数。印度尼西亚海尤其是ITF 不同路径上海水层结的时空特征不仅影响局地生态动力过程，也影响着西太平洋和东印度洋的水体交换。本文利用WOA13 气候态温盐数据和SODA 月平均温盐数据，给出了研究海域层结和跃层特征的空间三维分布及其季节和年际变化，对比分析了太平洋到印度洋的ITF 不同路径上海水层结特征的演化，还特别分析了层结的年际变化特征，并从机制上进行了探讨。

（1）在气候态平均下，ITF 各路径之间跃层内平均层结强度差异较小，其中东部路径平均值最小，为10-3.71s-2，西部路径次之，为10-3.70s-2，中部路径最大，为10-3.68s-2，同一路径上不同区域间平均层结强度差异也较小。各路径间跃层深度和厚度平均值存在明显差异，西部路径平均跃层深度和厚度最小，中部次之，东部最大。在印尼海海域，海水的层结强度、跃层深度和跃层厚度存在显著的空间分布不均，如卡里马塔海峡和爪哇海N2max小于10-4.1s-2，Dpcline小于50 m，明显小于印尼海其他海域。受强潮致混合作用影响，太平洋次表层高盐水进入印尼海后，水团性质发生显著变化，海水盐度明显降低，层结强度减弱，印尼海跃层平均深度和厚度均小于其他海域。

（2）受季风和ITF 流量季节性差异影响，研究海域层结的季节变化特征显著。西部路径夏季印尼海和印度洋海域上100 m 层的层结强度小于冬季，其余区域均大于冬季，其中印尼海层结季节性差异主要与卡里马塔海峡海水的输运方向有关。中部路径除太平洋100 m 以上区域，其余区域夏季层结强度均小于冬季；东部路径哈马黑拉海和塞兰海100 m 以上深度层夏冬季N2基本相同，其余海域夏季小于冬季。总体来看，在印尼海域内，爪哇海和卡里马塔海峡夏冬季N2max差异最为明显，最大可相差两个量级，其他大部分海域夏季N2max略小于冬季，差值多在0.5 量级以内；同时，夏季印尼海跃层深度大于冬季，夏季平均比冬季深30 m。

（3）对ITF 3 条路径上9 个特征点的N2max长期趋势分析显示，印尼海及周边大部分海域的层结强度呈现增大趋势，在23 a 内N2max最大增量接近0.1 个量级；这与Yamaguchi 和Suga[45]得出的结论基本相同。对特征点N2max和Dpcline的谱分析显示，各点除有较明显的1 a 周期变化，多数区域还有较为明显的4～7 a以及10 a 以上的周期信号，其中4～7 a 的周期变化可能主要与ENSO 事件有关。

本文主要分析了ITF 3 条路径及所经印尼海及周边海域层结的时空特征和变化，关于西太平洋和东印度洋通过印尼海和ITF 连接交换的过程和机制还有待进一步研究。