检索
E-mail
RSS

作者投稿 主编审稿 专家审稿 远程编辑

热带海洋学报 >

2019 , Vol. 38 >Issue 4: 20 - 29

DOI: https://doi.org/10.11978/2018113

海洋水文学

利用高分辨率水体反射地震资料研究吕宋海峡以东黑潮区混合

  • 张哲 , 1, 3 ,
  • 经志友 1 ,
  • 唐群署 2
展开
  • 1. 热带海洋环境国家重点实验室(中国科学院南海海洋研究所), 广东 广州 510301
  • 2. 中国科学院边缘海与大洋地质重点实验室(南海海洋研究所), 广东 广州 510301
  • 3. 中国科学院大学, 北京 100049;
张哲。E-mail:

张哲(1994—), 男, 山西省忻州市人, 硕士研究生, 主要从事反射地震方法计算混合参数研究。

Copy editor: 孙淑杰

收稿日期: 2018-11-04

  要求修回日期: 2019-01-17

  网络出版日期: 2019-07-21

基金资助

国家自然科学基金项目((41776040、41830538、41576070);)

热带海洋环境国家重点实验室自主研究项目((LTOZZ1701).)

版权

版权所有，未经授权，不得转载、摘编本刊文章，不得使用本刊的版式设计。
收起

Study on Kuroshio mixing located in the east of Luzon Strait by using high- resolution seawater seismic reflection data

  • ZHANG Zhe , 1, 3 ,
  • JING Zhiyou 1 ,
  • TANG Qunshu 2
Expand
  • 1. State Key Laboratory of Tropical Oceanography (South China Sea Institute of Oceanology, Chinese Academy of Sciences), Guangzhou 510301, China
  • 2. CAS Key Laboratory of Ocean and Marginal Sea Geology (South China Sea Institute of Oceanology), Guangzhou 510301, China
  • 3. University of Chinese Academy of Sciences, Beijing 100049, China;
ZHANG Zhe. E-mail:

Copy editor: SUN Shu-jie

Received date: 2018-11-04

  Request revised date: 2019-01-17

  Online published: 2019-07-21

Supported by

National Natural Science Foundation of China((41776040、41830538、41576070);)

Foundation of State Key Laboratory of Tropical Oceanography((LTOZZ1701).)

Copyright

Copyright reserved © 2019. Office of Acta Agronomica Sinica All articles published represent the opinions of the authors, and do not reflect the official policy of the Chinese Medical Association or the Editorial Board, unless this is clearly specified.
Fold

摘要

混合过程是海洋中普遍存在的一种形式, 对气候变化、物质分布等起到了重要作用。地震海洋学是近十多年发展起来的一门新兴学科, 被广泛应用到物理海洋学问题的研究中, 具有高空间分辨率的突出优点。文章利用反射地震资料, 通过斜率谱方法, 分别获得了吕宋海峡以东黑潮区湍流段与内波段的耗散率及扩散率。结果显示, 在剖面深度200~800m的平均耗散率为10 -7.0W·kg -1, 平均扩散率为10 -3.3m 2·s -1, 比大洋统计均值10 -5.0m 2·s -1高约1~2个量级, 与前人在吕宋海峡的观测结果相一致。湍流段和内波段的扩散率空间分布差异较大: 湍流段扩散率高值区对应强流区域, 推测这里是中尺度涡边缘, 其次中尺度不稳定过程引起扰动增强, 进而引起湍流混合的加强; 内波段扩散率高值区出现在吕宋岛弧附近, 推测是内波遇到岛弧地形发生破碎, 进而引起强的内波混合。

关键词： 黑潮区; 反射地震; 斜率谱; 混合

本文引用格式

张哲 , 经志友 , 唐群署 . 利用高分辨率水体反射地震资料研究吕宋海峡以东黑潮区混合[J]. 热带海洋学报, 2019 , 38(4) : 20 -29 . DOI: 10.11978/2018113

Abstract

Mixing is a ubiquitous motion in the ocean, which plays an important role in climate change and matter distribution. Seismic Oceanography is a new discipline in recent 10 years; it has been widely applied to Physical Oceanography, because it has the advantage of high spatial resolution. Seismic reflection data were used in this study, along with slope spectral method. We obtained the dissipation rate and diffusivity of turbulence interval and internal wave interval respectively in the Kuroshio (located in the east of the Luzon Strait). The results show that the average dissipation rate from 200- to 800-m depth in this profile is 10 -7.0W·kg -1, and the average diffusivity is 10 -3.3m 2·s -1, which is about 1-2 order of magnitude higher than the average value in the ocean 10 -5.0m 2·s -1. This result is consistent with previous mixing observations in the Luzon Strait. There is a big difference in diffusivity spatial distribution between turbulence interval and internal wave interval: the region where strong current occurs is consistent with that of high turbulence diffusivity in the profile, from which we can infer that this place may be the periphery of mesoscale eddy; its sub-mesoscale instability can strengthen the destabilization and then lead to strong mixing. The high internal wave diffusivity appears near the Luzon Arc, which indicates that the internal wave breaking occurs when it meets the island arc, leading to strong internal wave mixing.

Key words： Kuroshio; seismic reflection; slope spectra; mixing

吕宋海峡是连接南海与西太平洋的唯一深水通道, 对南海与西太平洋之间的淡水输送、热量交换等起到了至关重要的作用。大量的观测研究表明, 在垂直结构上, 吕宋海峡呈“类三明治”的结构, 即上层水和底层水向西流向南海, 中层水向东流向西太平洋(Tian et al, 2006; Sheu et al, 2009)。
黑潮作为西北太平洋最强的暖流, 携带着暖而咸的海水由赤道向北输送, 对吕宋海峡的水文特征和环流状况有很大影响, 尤其是各种形式的入侵(Hu et al, 2000; Caruso et al, 2006; Nan et al, 2015), 加速了海洋的混合过程(Yang et al, 2014; Mao et al, 2016)。
混合过程是海水的一种普遍存在形式, 即海水的各种特性(如热量、浓度、动量等)逐渐趋于均匀的过程(冯士筰 等, 1999)。海水混合有以下3种形式: 分子混合、湍流混合和对流混合, 其中湍流混合由海水微团随机运动造成, 交换强度比分子混合大许多量级, 是海洋混合的主要形式, 对流混合也通常表现为湍流混合状态(范植松, 2002)。因此, 湍流混合的研究是海洋混合研究中的关键问题。
海洋混合可以改变热量的分布和物质的输送, 影响海洋环流结构, 引起气候变化(Rahmstorf, 2003)。例如混合为大西洋南北方向上的经向翻转流提供了部分动力, 维持其运转所需的扩散率至少要达到10-4.0m2·s-1, 但实际大洋的观测结果仅为10-5.0m2·s-1, 这种“混合赤字”的现象在其他大洋中也存在。学者们希望通过小尺度的物理过程研究来弥补这一缺口。同时, 混合过程的研究对海洋参数化有重要作用。
众所周知, 海洋中存在不同尺度的运动, 能量由大尺度串级到小尺度, 经由湍流混合的形式耗散。所以对海洋混合的研究, 不仅可以解释海洋中能量的收支平衡问题, 还对海洋模式中耗散项的参数化有推动作用。除此之外, 海洋混合的研究对ENSO动力学、混合层动力学、海洋沉积动力学、河口海岸动力学、遥感海洋学、海洋声学及海洋生态学也有积极的推动作用(范植松, 2002)。在传统的海洋混合观测中, 海洋细微结构测量仪器如剪切流传感器、垂向微结构湍流剖面仪 (vertical microstructure profiler, VMP)、湍流微结构剖面仪(turbulence ocean microstructure acquisition profiler, TurboMAP)、锚定式湍流观测仪(moored turbulence measuring instrument, MTMI)等发挥了重要的作用。其他研究发现, 采用斜率谱分析方法也可以定量计算得到海洋中的耗散率和扩散率(Klymak et al, 2007a, b)。
地震海洋学是一门新兴的交叉学科, 是用反射地震的方法来研究海洋水体的结构特征, 自Holbrook等(2003)的工作以来, 用地震海洋学方法来研究物理海洋学问题得到了广泛的应用, 包括海洋锋面、水团边界、地转流、内波、涡旋、温盐阶梯、混合等。与传统的物理海洋学观测相比, 地震海洋资料具有高横向分辨率、短时间内对海水剖面快速成像的特点(宋海斌 等, 2008), 这些特点决定了它对物理海洋学中小尺度的研究具有极大优势。
本文利用反射地震的资料, 通过拾取反射同相轴信息用以代替水体波动, 计算得到水平斜率谱, 进而求出内波段和湍流段的耗散率和扩散率, 并对高值区产生原因做了推测分析。国内外在这方面也展开了广泛研究。Sheen等(2009)首次利用该方法计算了福克兰群岛东侧海域的一条地震剖面的海水耗散率和扩散率。Holbrook等(2013)利用反射地震资料计算了湍流耗散率, 给出了计算流程, 并详细讨论了影响耗散率计算的因素。拜阳 等(2015)分别以南海内波和地中海涡旋为研究对象, 计算了其中的混合参数, 并在斜率谱拟合方法上采用效果更好的集合经验模态分解方法。Falder等(2016)指出存在一个从内波到湍流段的动力学转化机制, 但内波段与湍流段之间没有重叠区域。Sallares等(2016)通过对阿尔沃兰海的地震资料作谱分析, 发现内波段与湍流段之间存在K-H不稳定, 对能量串级和次中尺度过程有了进一步认识。

1 数据资料

本文收集并使用了位于吕宋海峡以东黑潮区的4种数据资料: 反射地震数据, 用来获取海水等密线的位移; 卫星遥感海表面高度异常 (sea level anomaly, SLA)数据, 用来得到航次测线周围的海流情况; Argo温、盐数据, 用来计算不同深度海水的浮力频率; 垂向混合坐标海洋模式 (hybrid coordinate ocean model, HYCOM)模式数据, 用来得到航次测线在不同深度上的流速。数据详细介绍如下。

1.1 反射地震数据

2001年8月4日至6日, 由国家973“中国边缘海形成演化及其重大资源的关键问题”项目委托广州海洋地质调查局, 利用“探宝号”调查船在南海东北部进行了地震数据采集工作。整个测线横穿中国大陆架、大陆坡、南海东北部、恒春海脊、吕宋岛弧、花东海盆和加瓜海脊, 如图1a中黑色线段所示。其中红色部分为本文研究所用的测线段, 总长180km。地震数据使用道间距为12.5m的240道水听器记录, 采样间隔为2ms, 记录长度为10s, 震源位于水深8m处, 容量为3000in3(约0.049m3)的枪阵, 炮间距为50m, 近炮检距为250m。
图1 2001年8月6日航次测线周围海域水深地形(a)和SLA(海表面高度异常)(b)

黑色线段表示完整的航次测线, 红色线段表示本文研究使用的部分测线, 箭头表示由SLA计算得到的地转流

Fig. 1 Bathymetry of the region around the survey line (a) and SLA (b) during the cruise on 6th August 2001. The black line indicates the whole seismic transect during the cruise, the red line indicates the part of the transect used in this study, and the black arrows indicate the geostrophic current calculated by SLA

地震数据的处理流程如下: 1)观测系统定义; 2)去除坏道; 3)去除直达波; 4)带通滤波处理(15~90Hz); 5)结合速度分析的叠前深度偏移; 6)偏移数据叠加; 7)谐波噪声压制(Tang et al, 2011; Holbrook et al, 2013)。最终得到了150~1500m深度范围的反射地震图像(图2a)。
进一步利用希尔伯特-黄(Hilbert-Huang)变换对图2a中的反射地震图像进行同相轴的自动追踪拾取(图2b)。由于观测系统最小偏移距的影响, 在反射地震剖面150m以上的表层看不到有效反射, 也拾取不到同相轴, 为保守起见, 在数据处理分析中剔除0~200m的上层同相轴。同时发现, 在左下侧吕宋岛弧地形处有错误拾取到的同相轴; 800m以下的同相轴由于噪声的干扰, 倾角过大且不符合物理规律, 同样属于错误拾取到的同相轴。将上述错误的同相轴剔除之后, 得到本文研究所用水体的同相轴, 深度范围200~800m, 水平跨度180km(图2c)。
图2 水体反射地震剖面图像(a)、追踪拾取到的同相轴(b)、本文研究使用的同相轴(c)

对应图1中红色测线段

Fig. 2 Seismic reflection profile (a), automatically tracked reflections (b), and tracked reflections used in this study. The profile corresponds to the red line in Fig. 1

1.2 卫星遥感SLA数据

本文使用的卫星高度计资料由法国国家空间研究中心AVISO (Archiving Validation and Interpolation of Satellite Oceanography Data; http://marine. copernicus.eu/)提供。数据内容为网格化的海表面高度异常(SLA), 时间分辨率为日平均, 空间分辨率为0.25°×0.25°。文中选取的是2001年8月6日的数据, 与研究区域测线走航时间一致。按照地转理论中压强梯度力与科氏力相平衡(Stewart, 2008), 可以将上述SLA数据计算得到地转流, 计算公式为:
$\left\{ \begin{matrix} u=-\frac{g}{f}\frac{\partial h}{\partial y} \\ v=\frac{g}{f}\frac{\partial h}{\partial x} \\\end{matrix} \right.$
其中, u、v分别是地转流沿经向和纬向的速度; g是重力加速度, 取值为9.8; f是科氏参数, h是海表面高度异常, $\frac{\partial h}{\partial x}$和$\frac{\partial h}{\partial y}$分别表示h在经向和纬向上的偏导。SLA填色及地转流如图1b所示。

1.3 Argo温、盐数据

本文用来计算浮力频率所使用的温、盐数据来自中国Argo 实时资料中心(http://www.argo.org.cn/), 选择的数据产品为《Argo网格化温盐产品(G-Argo)》, 时间范围为2000—2008年, 时间分辨率为每周一次, 空间范围为120°E—90°W、90°S—60°N, 空间分辨率为3°×3°, 垂向0~2000m共36层。
本文选择测线周围4个节点: (121°30′E, 19°30′N)、(124°30′E, 19°30′N)、(121°30′E, 22°30′N)、(124°30′E, 22°30′N), 通过这4个节点将2000年1月5日至2008年12月31日共计470周的数据计算得到浮力频率值, 先后作时间平均和4个节点的平均, 将结果作为测线不同深度的浮力频率值(图3), 用于后面水平斜率谱的校正。
图3 浮力频率随深度的变化曲线

竖直红线表示参考浮力频率值(3.0cph), 竖直绿线表示研究区域浮力频率均值(3.9cph), cph为周频率cycles per hour的缩写,1cph=(1/3600)s

Fig. 3 Buoyancy frequency curve varying with depth. The vertical red line indicates the reference buoyancy frequency (3.0cph), and the vertical green line indicates the average buoyancy frequency (3.9cph) in the study region, cph is the abbreviation of cycles per hour, 1cph=(1/3600)s

1.4 HYCOM模式数据

本文使用的模式由HYCOM提供, 下载网址为https://hycom.org/, 选择下载再分析的流场数据, 数据全称为Global Ocean Forecasting System (GOFS) 3.1 output on the GLBv0.08 grid, 时间范围为1994年1月1日至2015年12月31日, 时间分辨率为3h; 在空间上: 40°S—40°N间的空间分辨率为0.08°, 其余区域为0.04°。最后选取与走航测线同期的流场数据。

2 计算方法

2.1 内波段与湍流段区间划分

本文选取数据中水平长度大于1.25km的同相轴, 共计916条, 将每一条同相轴做距平处理(即减去其所在深度的平均值), 得到同相轴的垂直位移。由于同相轴位移近似于等密线位移(Holbrook et al, 2013), 因此用同相轴位移代替等密线位移来作斜率谱分析。
采用Welch方法并选取128个点的傅里叶变换长度和半重叠的Hanning窗, 对上述每一条同相轴位移计算其功率谱密度, 得到水平波数谱。每条水平波数谱均乘以(2πk)2(其中k为波数), 得到水平斜率谱。这样做能够更好地把内波段和湍流段区分开来, 因为在双对数坐标下, 斜率谱内波段斜率为-1/2 (Garrett et al, 1975), 湍流段斜率为1/3 (Batchelor, 1959), 正负值交替更容易区分。除此之外, 由于内波的垂向位移依赖于局地的海水层结, 因此斜率谱需要通过浮力频率校正来消除层结的影响(Sheen et al, 2009), 即在水平斜率谱上乘以(N/N0), 其中N为该同相轴所在深度的浮力频率, N0=3cph为参考浮力频率。以下分析均使用校正后的水平斜率谱。
上述916条水平斜率谱的平均如图4所示, 左侧内波段(斜率为-1/2)和右侧湍流段(斜率为1/3)分界明显。由于内波段和湍流段之间转换部分的动力过程同时受内波和湍流影响(Sheen et al, 2009; Falder et al, 2016), 斜率谱最低波数段受带宽的限制, 最高波数段受噪声的影响, 因此本文剔除上述区间, 从而划定内波段范围为: 0.0013~0.0044cpm (cpm为周波数cycles per meter的缩写), 湍流段范围为: 0.0088~0.03cpm(图4), 并应用于所有水平斜率谱(图5)。

2.2 混合参数估算

完整的斜率谱模型包括3部分的运动形式: 内波、惯性对流和惯性扩散(Klymak et al, 2007a, b):
$\begin{align} {{\phi }_{\zeta }}\left( {{k}_{x}} \right)={{\phi }_{0}}k_{x}^{-1/2}+ 4\text{ }\!\!\pi\!\!\text{ }\frac{\Gamma \varepsilon }{N_{0}^{2}}\left[ {{C}_{T}}{{\varepsilon }^{-1/3}}{{\left( 2\text{ }\!\!\pi\!\!\text{ }{{k}_{x}} \right)}^{1/3}}+q{{\nu }^{1/2}}{{\varepsilon }^{-1/2}}\left( 2\text{ }\!\!\pi\!\!\text{ }{{k}_{x}} \right) \right] \\ \end{align}$
其中, $4\text{ }\!\!\pi\!\!\text{ }\frac{\Gamma \varepsilon }{N_{0}^{2}}$$\left[ {{C}_{T}}{{\varepsilon }^{-1/3}}{{\left( 2\text{ }\!\!\pi\!\!\text{ }{{k}_{x}} \right)}^{1/3}} \right]$是惯性对流子区斜率谱, 存在于湍流低波数段; $4\text{ }\!\!\pi\!\!\text{ }\frac{\Gamma \varepsilon }{N_{0}^{2}}\left[ q{{\nu }^{1/2}}{{\varepsilon }^{-1/2}}\left( 2\text{ }\!\!\pi\!\!\text{ }{{k}_{x}} \right) \right]$是惯性扩散子区斜率谱, 存在于湍流高波数段; ${{\phi }_{0}}k_{x}^{-1/2}$是内波段斜率谱;。因为惯性扩散子区波数范围高于反射地震资料的空间分辨率, 本文地震数据的斜率谱仅分为内波段和惯性对流子区两段(图4)。
图4 水平斜率谱的平均

共使用916条水平斜率谱计算谱的平均, 误差棒表示95%置信区间; 竖直虚线内为内波段区间, 竖直实线内为湍流段区间; 两条倾斜线段斜率为1/3和-1/2, 分别用来拟合湍流段和内波段的斜率谱

Fig. 4 The average horizontal slope spectra. A total of 916 slope spectra are used to calculate average slope spectra, and error bars represent 95% bootstrap confidence interval. The internal wave interval is within the two vertical dashed lines, and the turbulence interval is within the two vertical solid lines. The slopes of two tilted lines are 1/3 and -1/2, respectively, which are used for fitting the slope spectra of turbulence and internal wave intervals

在高波数的湍流段, 波数谱的斜率与类Kolmogorov律(-5/3)相一致, 采用Batchelor模型(Batchelor, 1959)拟合, 即湍流惯性对流子区:
$\psi _{\zeta }^{T}=\frac{4\text{ }\!\!\pi\!\!\text{ }\Gamma }{{{N}^{2}}}{{C}_{T}}{{\varepsilon }^{2/3}}{{\left( 2\text{ }\!\!\pi\!\!\text{ }k \right)}^{-5/3}}$
其中, $\psi _{\zeta }^{T}$是湍流段波数谱; Γ是混合效率, 根据经验取值为0.2; N是斜率谱所在深度的浮力频率; ${{C}_{T}}$是常数, 约为0.4(Sreenivasan, 1996); $\varepsilon $是湍流耗散率; k是水平波数。
在低波数的内波段, 功率谱与内波GM75模型(Garrett et al, 1975)相一致, 波数谱斜率为-2.5, 采用Gregg-Henyey参数化模型拟合。有充足证据表明, 在开阔大洋中, 耗散率与内波能量的二次方成正比, 根据内波的波-波相互作用理论, ε等于其在高波数段通过内波谱的净能通量, 据此得到与观测结果符合最好的平均湍动能耗散率计算公式(Henyey et al, 1986, Gregg, 1989):
$\left\langle \varepsilon \right\rangle =7\times {{10}^{-10}}\frac{{{N}^{2}}}{N_{0}^{2}}{{\left\langle \frac{\psi _{\zeta }^{\text{IW}}}{\psi _{\zeta }^{\text{GM}}} \right\rangle }^{2}}$
其中$\left\langle \varepsilon \right\rangle $是待求的内波耗散率; ${{N}_{0}}=3\text{cph}$是大洋参考浮力频率; 7×10-10W·kg-1是GM75谱$\psi _{\zeta }^{\text{GM}}$在纬度为30°的背景耗散率; $\psi _{\zeta }^{\text{IW}}$是内波段斜率谱; $\psi _{\zeta }^{\text{GM}}$是GM75内波斜率谱; $\left\langle \ldots \right\rangle $表示在内波段内求平均。
在斜率谱中, 湍流段给定斜率为1/3, 其耗散率根据最小二乘法原理和公式(3)拟合得到; 内波段耗散率由公式(4)的参数化方法计算得到(图5)。由于反射地震数据采集方法具有快速横穿断面的特点, 波场随时间变化缓慢, 因此计算得到的耗散率反映了能量从大尺度向小尺度的转化速率(Sheen et al, 2009)。
图5 两条典型的水平斜率谱

误差棒表示95%置信区间, 竖直虚线和竖直实线的含义同图4, 内波段和湍流段斜率分别为-1/2、1/3(红色线段)

Fig. 5 Two typical slope spectra. Error bars represent 95% bootstrap confidence intervals. The meaning of vertical dashed lines and vertical solid lines are explained in Fig.4. The slopes of internal wave interval and turbulence interval are -1/2 and 1/3, respectively (as shown by red lines)

在假设研究区域处于稳定剪切状态的前提下, Osborn经验公式(Osborn, 1980)可将耗散率转化成跨密度面扩散率, 且适用于湍流段和内波段:
${{K}_{\rho }}=\frac{\Gamma \varepsilon }{{{N}^{2}}}$
其中, ${{K}_{\rho }}$是跨等密度面扩散率; Γ是混合效率, 根据经验取值为0.2; $\varepsilon $是耗散率。因此, 湍流段和内波段的扩散率由公式(5)进一步计算得到。

3 结果与讨论

3.1 混合参数空间分布

将每一条同相轴反演得到的结果填图到其所在的线段上, 得到该区域湍流段和内波段的耗散率及扩散率, 结果如图6和图7所示。研究区域的湍流段和内波段混合参数统计均值如下: 湍流段的耗散率和扩散率均值分别为8.7×10-8W·kg-1、4.7×10-4m2·s-1; 内波段的耗散率和扩散率均值分别为1.9×10-8W·kg-1、1.1×10-4m2·s-1。两段的扩散率与开阔大洋背景场相比均高一个量级以上。
图6 湍流段计算得到的耗散率(a)和扩散率(b)

图中虚线方框内为混合强烈区域

Fig. 6 Calculated dissipation rate (a) and diffusivity (b) by turbulence interval. The rectangle marks the strong mixing area

图7 内波段计算得到的耗散率(a)和扩散率(b)

图中虚线方框内为混合强烈区域

Fig. 7 Calculated dissipation rate (a) and diffusivity (b) by internal wave interval. The rectangle marks the strong mixing area

本文对研究区域利用水体反射地震资料, 拾取其同相轴代替等密线, 通过斜率谱的方法计算得到该区域湍流段和内波的耗散率及扩散率。此前, 也有物理海洋学家用其他方法对南海东北部和吕宋海峡的湍流混合进行研究, 如Tian等(2009)在IWME (Internal Wave and Mixing Experiment)实验中得到该海域两个剖面的湍流扩散率。Tian等使用低频声学多普勒流速剖面仪(lowered acoustic Doppler current profiler, LADCP)、温盐深剖面测量仪(conductance temperature depth, CTD)等仪器获得剖面的流场和温盐数据, 通过参数化模型计算得到湍流扩散率, 并与垂向微结构湍流剖面仪vertical microstructure profiler, VMP)实测数据对比分析, 验证了参数化结果的可靠性。从结果上来看, Tian等在与本文地震测线相近位置上得到的湍流扩散率量级约为10-4m2·s-1, 与本文的结果相近, 表明本文计算得到的混合参数值是可靠的。

3.2 混合机制分析

南海与太平洋之间的水交换主要由吕宋口上300m的水体控制, 夏季8月在整个吕宋海峡的表层流是流向太平洋的, 受南海气旋涡影响, 黑潮的入侵被抑制, 表现为“跳跃(leaping)”状路径, 在吕宋口北部300m以下有向外流出的弱流(Shaw et al, 1994; Nan et al, 2011)。卫星高度计资料（图1b）表明, 本文的研究区域恰好在黑潮入侵处。测线段处于黑潮主体活动范围, 且西侧连接吕宋岛弧。
对比图6b和图7b, 湍流段扩散率和内波段扩散率的空间分布存在显著差异, 表明剖面不同区域的海水混合受不同尺度动力过程的影响。湍流扩散率(图6b)显示, 剖面水平90~130km和163~180km处(黑色虚线框所示)的湍流扩散率明显增强, 分别为7.8×10-4m2·s-1和8.3×10-4m2·s-1, 远高于湍流段平均值4.7×10-4m2·s-1表明该区域存在明显的湍流混合现象。而该区域湍流耗散率同样存在高值, 表明有明显的能量正向串级过程发生(图6a)。卫星高度计资料(图1b)和HYCOM模式再分析流场数据(图8)显示, 湍流扩散率高值区分别与测线周围的反气旋涡和气旋涡边缘区域一致。而该区域无法探测到的有效反射(黑色虚线框中的空白)(Holbrook et al, 2009), 可能是由海水充分混合并趋于均匀引起。已有的观测和模式结果指出, 涡旋边缘锋面区域存在次中尺度过程, 能显著增强水体垂向交换, 是能量从中尺度向小尺度串级的有效途径之一(D’Asaro et al, 2011; Thomas et al, 2013; Mackinnon et al, 2013; Yang et al, 2017; 郑瑞玺 等, 2018)。因此, 湍流混合的增强可能与中尺度涡边缘的次中尺度过程有关。
图8 HYCOM模式2001年8月6日6:00时刻不同深度测线流速(a)和流速剖面(b)

Fig. 8 Current velocity with different depths (a) and current profile (b) of the line at 6:00 am on 6th August 2001, both from the HYCOM

内波扩散率分布(图7b)显示, 水平0~20km间(黑色虚线框所示)的内波扩散率明显高于其他区域, 达到6.4×10-4m2·s-1。该区域与吕宋岛弧相接壤, 处于吕宋海峡。而吕宋海峡受潮汐、地形和海水层结的共同影响, 是全球海洋内波运动最为活跃的区域之一(Simmons et al, 2011; Zhao, 2014)。内波在生成和传播过程中会发生破碎, 并引起跨等密面扩散与混合, 且在粗糙地形处的强烈混合能够使内波能量串级至耗散尺度(Müller et al, 2000; Thorpe, 2010; van Haren et al, 2012; Pinkel et al, 2012)。因此吕宋岛弧附近内波段扩散率的显著增强与其所处的地形密切相关, 其附近活跃的内波活动及内波破碎过程是该区域海水混合和能量正向串级的重要原因之一。

4 结论

本文利用高分辨率的地震反射资料, 对吕宋海峡以东黑潮入侵南海区域的混合情况进行了初步研究, 通过斜率谱的方法反演得到该区域内波段和湍流段的耗散率及扩散率。分析结果表明, 吕宋海峡以东黑潮区域的混合现象明显, 且在斜率谱上呈现出内波段和湍流段两段区间, 湍流段扩散率平均值为10-3.3m2·s-1, 内波段扩散率平均值为10-4m2·s-1, 比大洋平均扩散率10-5m2·s-1高1~2个量级。从计算结果中看到, 吕宋岛弧的存在有利于内波破碎并诱发强烈的内波混合。结合卫星高度计和HYCOM流场资料分析, 发现观测期间该海域存在多个中尺度涡旋, 在其边缘的次中尺度不稳定加强了观测区域的湍流扰动, 并导致湍流混合的增强。本文仅仅应用了反射地震资料对研究区域的混合现象进行了分析, 希望今后能在地震数据采集的同时进行同步水文观测, 进一步利用多学科观测资料分析探讨海洋中小尺度过程及其相互作用。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
拜阳, 宋海斌, 董崇志 , 等, 2015. 地震海洋学方法在海洋混合参数提取中的研究与应用——以南海内波和地中海涡旋为例[J]. 地球物理学报, 58(7):2473-2485.

BAI YANG, SONG HAIBIN, DONG CHONGZHI , et al, 2015. Extraction of mixing parameters by seismic oceanography and applications: Case study of the internal waves in South China Sea and Mediterranean eddy[J]. Chinese Journal of Geophysics, 58(7):2473-2485 (in Chinese with English abstract).

[2]
范植松 , 2002. 海洋内部混合研究基础[M]. 北京: 海洋出版社.

FAN ZHISONG. Research fundamentals of ocean interior mixing[M]. Beijing: China Ocean Press (in Chinese).

[3]
冯士筰, 李凤岐, 李少菁 , 1999. 海洋科学导论[M]. 北京: 高等教育出版社.

FENG SHIZUO, LI FENGQI, LI SHAOJING. An introduction to marine science[M]. Beijing: Higher Education Press (in Chinese).

[4]
宋海斌, 董崇志, 陈林 , 等, 2008. 用反射地震方法研究物理海洋—地震海洋学简介[J]. 地球物理学进展, 23(4):1156-1164.

SONG HAIBIN, DONG CHONGZHI, CHEN LIN , et al, 2008. Reflection seismic methods for studying physical oceanography: Introduction of seismic oceanography[J]. Progress in Geophysics, 23(4):1156-1164 (in Chinese with English abstract).

[5]
郑瑞玺, 经志友, 罗士浩 , 2018. 南海北部反气旋涡旋边缘的次中尺度动力过程分析[J]. 热带海洋学报, 37(3):19-25.

ZHENG RUIXI, JING ZHIYOU, LUO SHIHAO , 2018. Analysis of sub-mesoscale dynamic processes in the periphery of anticyclonic eddy in the northern South China Sea[J]. Journal of Tropical Oceanography, 37(3):19-25 (in Chinese with English abstract).

[6]
BATCHELOR G K , 1959. Small-scale variation of convected quantities like temperature in turbulent fluid Part 1. General discussion and the case of small conductivity[J]. Journal of Fluid Mechanics, 5(1):113-133.

DOI

[7]
CARUSO M J, GAWARKIEWICZ G G, BEARDSLEY R C , 2006. Interannual variability of the Kuroshio intrusion in the South China Sea[J]. Journal of Oceanography, 62(4):559-575.

DOI

[8]
D'ASARO E, LEE C, RAINVILLE L , et al, 2011. Enhanced turbulence and energy dissipation at ocean fronts[J]. Science, 332(6027):318-322.

DOI

[9]
FALDER M, WHITE N J, CAULFIELD C P , 2016. Seismic imaging of rapid onset of stratified turbulence in the south Atlantic Ocean[J]. Journal of Physical Oceanography, 46(4):1023-1044.

DOI

[10]
GARRETT C, MUNK W , 1975. Space-time scales of internal waves: a progress report[J]. Journal of Geophysical Research, 80(3):291-297.

DOI

[11]
GREGG M C , 1989. Scaling turbulent dissipation in the thermocline[J]. Journal of Geophysical Research: Oceans, 94(C7):9686-9698.

DOI

[12]
HENYEY F S, WRIGHT J, FLATTÉ S M , 1986. Energy and action flow through the internal wave field: an eikonal approach[J]. Journal of Geophysical Research: Oceans, 91(C7):8487-8495.

DOI

[13]
HOLBROOK W S, PÁRAMO P, PEARSE S , et al, 2003. Thermohaline fine structure in an oceanographic front from seismic reflection profiling[J]. Science, 301(5634):821-824.

DOI

[14]
HOLBROOK W S, FER I, SCHMITT R W , 2009. Images of internal tides near the Norwegian continental slope[J]. Geophysical Research Letters, 36(24): L00D10.

[15]
HOLBROOK W S, FER I, SCHMITT R W , et al, 2013. Estimating oceanic turbulence dissipation from seismic images[J]. Journal of Atmospheric and Oceanic Technology, 30(8):1767-1788.

DOI

[16]
HU JIANYU, KAWAMURA H, HONG HUASHENG , et al, 2000. A review on the currents in the South China Sea: seasonal circulation, South China Sea warm current and Kuroshio intrusion[J]. Journal of Oceanography, 56(6):607-624.

DOI

[17]
KLYMAK J M, MOUM J N , 2007a . Oceanic isopycnal slope spectra. Part I: internal waves[J]. Journal of Physical Oceanography, 37(5):1215-1231.

DOI

[18]
KLYMAK J M, MOUM J N , 2007b . Oceanic isopycnal slope spectra. Part II: turbulence[J]. Journal of Physical Oceanography, 37(5):1232-1245.

DOI

[19]
MACKINNON J, ST LAURENT L, GARABATO A C N , 2013. Diapycnal mixing processes in the ocean interior[J]. International Geophysics, 103:159-183.

DOI

[20]
MAO H, SUN Z, LI X, 2016. Observations of water mass mixing in the Luzon strait [C]//Proceedings of the American Geophysical Union, Ocean Sciences Meeting, 2016, abstract# A34B-2648.

[21]
MÜLLER P, BRISCOE M , 2000. Diapycnal mixing and internal waves[J]. Oceanography, 13(2):98-103.

DOI

[22]
NAN F, XUE H J, CHAI F , et al, 2011. Identification of different types of Kuroshio intrusion into the South China Sea[J]. Ocean Dynamics, 61(9):1291-1304.

DOI

[23]
NAN FENG, XUE HUIJIE, YU FEI , 2015. Kuroshio intrusion into the South China Sea: a review[J]. Progress in Oceanography, 137:314-333.

DOI

[24]
OSBORN T R , 1980. Estimates of the local rate of vertical diffusion from dissipation measurements[J]. Journal of Physical Oceanography, 10(1):83-89.

DOI

[25]
PINKEL R, BUIJSMAN M, KLYMAK J M , 2012. Breaking topographic lee waves in a tidal channel in Luzon Strait[J]. Oceanography, 25(2):160-165.

DOI

[26]
RAHMSTORF S , 2003. Thermohaline circulation: the current climate[J]. Nature, 421(6924):699.

DOI

[27]
SALLARES V, MOJICA J F, BIESCAS B , et al, 2016. Characterization of the submesoscale energy cascade in the Alboran Sea thermocline from spectral analysis of high- resolution MCS data[J]. Geophysical Research Letters, 43(12):6461-6468.

DOI

[28]
SHAW P T, CHAO S Y , 1994. Surface circulation in the South China sea[J]. Deep Sea Research Part I: Oceanographic Research Papers, 41(11-12):1663-1683.

DOI

[29]
SHEEN K L, WHITE N J, HOBBS R W , 2009. Estimating mixing rates from seismic images of oceanic structure[J]. Geophysical Research Letters, 36(24): L00D04.

[30]
SHEU D D, CHOU W C, CHEN C T A , et al, 2009. Riding over the Kuroshio from the South to the East China Sea: mixing and transport of DIC[J]. Geophysical Research Letters, 36(7):L07603.

[31]
SIMMONS H, CHANG M H, CHANG Y T , et al, 2011. Modeling and prediction of internal waves in the South China Sea[J]. Oceanography, 24(4):88-99.

DOI

[32]
SREENIVASAN K R , 1996. The passive scalar spectrum and the Obukhov-Corrsin constant[J]. Physics of Fluids, 8(1):189-196.

DOI

[33]
STEWART R H , 2008. Introduction to physical oceanography[M]. Texas: Texas A & M University.

[34]
TANG QUNSHU, ZHENG CHAN , 2011. Thermohaline structures across the Luzon Strait from seismic reflection data[J]. Dynamics of Atmospheres and Oceans, 51(3):94-108.

DOI

[35]
THOMAS L N, TAYLOR J R, FERRARI R , et al, 2013. Symmetric instability in the gulf stream[J]. Deep Sea Research Part II: Topical Studies in Oceanography, 91:96-110.

DOI

[36]
THORPE S A , 2010. Breaking internal waves and turbulent dissipation[J]. Journal of Marine Research, 68(6):851-880.

DOI

[37]
TIAN JIWEI, YANG QINGXUAN, LIANG XINFENG , et al, 2006. Observation of Luzon Strait transport[J]. Geophysical Research Letters, 33(19):L19607.

DOI

[38]
TIAN JIWEI, YANG QINGXUAN, ZHAO WEI , 2009. Enhanced Diapycnal mixing in the South China Sea[J]. Journal of Physical Oceanography, 39(12):3191-3203.

DOI

[39]
VAN HAREN H, GOSTIAUX L , 2012. Energy release through internal wave breaking[J]. Oceanography, 25(2):124-131.

DOI

[40]
YANG QINGXUAN, ZHOU LEI, TIAN JIWEI , et al, 2014. The roles of Kuroshio intrusion and mesoscale eddy in upper mixing in the northern South China Sea[J]. Journal of Coastal Research, 30(1):192-198.

[41]
YANG QINGXUAN, ZHAO WEI, LIANG XINFENG , et al, 2017. Elevated mixing in the periphery of mesoscale eddies in the South China Sea[J]. Journal of Physical Oceanography, 47(4):895-907.

DOI

[42]
ZHAO ZHONGXIANG , 2014. Internal tide radiation from the Luzon Strait[J]. Journal of Geophysical Research: Oceans, 119(8):5434-5448.

DOI

Options
文章导航

/

〈 〉