楊順, 宋海斌, 范文豪, 吳迪
海洋地質國家重點實驗室, 同濟大學海洋與地球科學學院, 上海 200092
隨著觀測技術的發展,特別是衛星遙感技術在20世紀70年代以來的進展,人們逐漸認識到中尺度渦旋運動才是海洋表層運動的主體.單個渦旋可以攜帶其核心水,并可以長距離輸送核心水的性質.因此,海洋渦旋對于全球海洋中熱量、CO2、營養鹽等物質與能量的輸送與再分配起著至關重要的作用,進而能調節全球氣候變化和影響海洋資源分布(Zhang et al., 2014).近年來隨著觀測與模擬手段的進步,關于亞中尺度渦的研究也越來越多(Su et al., 2018;Onken et al., 2020;Gula et al., 2016).與中尺度渦相比,亞中尺度渦的空間尺度和時間尺度都更小,對觀測和模擬都提出了更高的要求.亞中尺度過程具有較強的垂向流速,對上層海洋的營養物質做了重要補給(Mcgillicuddy et al., 2007).與大中尺度洋流的地轉平衡假設不同,亞中尺度流通常超越了準地轉這一理論框架,同時仍然受到地球自轉和海洋中普遍存在的穩定密度層結的強烈影響(McWilliams, 2016).
地震海洋學是一門利用傳統反射地震方法研究物理海洋現象及其過程的交叉學科(宋海斌 等, 2008;Holbrook et al., 2003).它利用不同深度海水層的溫度和鹽度差異引起海水聲速和密度差異產生的波阻抗差,從而對海水層進行成像.海水層地震圖像可以近似看做是平滑后的海水垂向溫度梯度(Ruddick et al., 2009).相比于傳統物理海洋學方法,地震海洋學具有采集速度快、橫向分辨率高、海水全深度成像等優點.目前,地震海洋學已經在各種定性和定量問題上做出了大量研究,比如海洋鋒面(Tsuji et al., 2005;Nakamura et al., 2006)、內波(Holbrook and Fer, 2005;Tang et al., 2014b)、溫鹽參數反演(Blacic et al., 2016;Papenberg et al., 2010)、混合參數提取(Holbrook et al., 2013;Sheen et al., 2009)、地轉流速估算(Sheen et al., 2011;黃興輝等, 2013)等.
已有很多學者利用地震海洋學方法對渦旋進行了研究,早期主要為地中海渦旋(Meddy).Biescas等(2008)首次對地中海渦旋的精細結構進行地震成像,揭示了地中海渦旋的透鏡狀形態特征以及核心弱、邊界強的反射特征,打開了利用反射地震方法研究海洋渦旋的大門.Pinheiro等(2010)利用地震剖面結合聲浮子、SLA(海表高度異常)、SST(海表溫度)資料對地中海渦旋以及地中海潛流進行了研究,揭示了傳統物理海洋學測量不容易觀察到的熱鹽結構橫向變化的復雜性與細節.Song等(2011)首次發現地中海渦旋旋臂,并對渦旋以及其混合動力機制進行了研究,結果顯示渦旋核心的混合作用缺失,而渦旋頂底及周圍鋒區的混合異常高.后來又不斷有人利用反射地震數據反演出地中海渦旋的溫度、鹽度和位密斷面(Papenberg et al., 2010; 黃興輝 等, 2011; Biescas et al., 2014),這加強了物理海洋學與反射地震學之間的聯系,使得雙方對于對方的研究有了更深的認識和理解.近年來,隨著研究的深入,有學者開始對其他海域的渦旋進行地震成像.黃興輝等(2013)基于地轉平衡假設從反射同相軸估算了南海一個中尺度渦的地轉流速.傳統測量海水流速利用ADCP(聲學多普勒海流剖面儀)獲得垂向剖面的絕對海水流速,或者利用衛星遙感獲得的海表高度數據基于地轉平衡假設獲得海表地轉流速數據.ADCP測量深度較淺,衛星遙感只能獲得海表面的地轉流速,二者都存在較大局限性.而從地震反射同相軸估算的海水流速數據則可以貫穿整個海水深度,極大彌補了物理海洋學對流速測量的缺陷.Tang等(2014a)結合地震數據、衛星高度數據和原位水文數據對阿拉斯加灣的渦旋進行了綜合分析,三者表現出較好的一致性.Gunn等(2018)對南大洋Bellingshausen海的地震成像表明這里存在大量中尺度到亞中尺度暖渦,揭示了該海域的溫鹽環流特征,并且表明以前對該海域渦旋的出現頻率可能被低估了1個數量級,這對計算南極半島西部陸架冰質量損失具有重要意義.Gula等(2019)利用地震數據對灣流附近的一個透鏡狀亞中尺度相干渦進行了成像,并結合Glider剖面、同步采集的XBT/CTD以及亞中尺度分辨率的流體動力學數值模擬分析了其產生機制,模擬結果和地震剖面特征較一致,渦旋是由灣流和海底山脊的摩擦作用形成的.
大洋邊界廣泛存在邊界流,例如典型的太平洋西邊界的黑潮和大西洋西邊界的灣流,邊界流和地形相互作用形成了大量各種尺度的渦旋.前人對地中海西部、南海、阿拉斯加灣、南大洋和灣流區的渦旋進行了地震成像,由于水團性質存在較大差異,不同海域渦旋反射特征各不相同.雖然地震海洋學對渦旋的研究已較為廣泛,但更多大洋邊界處的渦旋還未被對其進行高分辨率反射地震成像,因此對研究海域進行拓展十分必要.
渦旋可能從海面延伸至海底,也可能集中在表層,或者只出現在海水次表層(Talley et al., 2011).這里,我們利用多道反射地震數據在中美洲鸚鵡灣捕捉到一個次表層透鏡狀反射結構,結合再分析水文數據,認為這是一個氣旋渦,并對其反射特征進行了分析.最后結合同航次采集的其他地震剖面,分析了該氣旋渦不同部位的反射特征.
研究區域位于中美洲尼加拉瓜岸外的鸚鵡灣(圖1).赤道逆流(North Equatorial Counter Current,NECC)自西向東到大洋東邊界遇到陡峭的陸坡后向北轉向,并在此形成上升流區,即哥斯達黎加穹窿區(Costa Rica Dome,CRD)(Reyes-Hernández et al., 2016;Brenes et al., 2008).NECC的一部分在中美洲海岸,匯入了哥斯達黎加沿岸流(Costa Rica Coastal Current,CRCC).CRCC位于CRD和海岸之間,并沿海岸向西北流去,經過Tehuantepec灣后離開海岸(Fiedler, 2002).由于局地表層風驅動作用或者海底地形的摩擦作用,該海域發育有大量各種尺度的氣旋渦和反氣旋渦(Reyes-Hernández et al., 2016).這里的主要水團有:熱帶表層水(TSW)、副熱帶底層水(STUW)和南極中層水(AAIW).水團各自溫鹽范圍大致為:TSW,T>25 ℃,S<34;STUW,13 ℃ 圖1 (a)尼加拉瓜岸外地形圖(彩色底圖為高程圖,黑色實線為L115地震測線,黃色五角星表示地震采集起點,紅色圓點表示歷史CTD測站位置);(b)L115地震測線采集期間(2000年6月23日)的日平均海表高程(彩色底圖)和日平均海水流速(灰色箭頭),來自CMEMSFig.1 (a) Topography offshore Nicaragua. Base map is the elevation, black line is seismic survey line L115, yellow pentagram is the starting point of seismic acquisition, and red dots are the location of historical CTD stations. (b) Daily average sea surface height (color map) and daily average sea water velocity (gray arrows) during the acquisition of L115 seismic survey line (June 23, 2000), data from CMEMS 圖2 研究區域歷史CTD數據,來自NCEI(a) 溫度剖面; (b) 鹽度剖面; (c) 聲速剖面, 藍線為擬合聲速; (d) T-S圖(灰色等值線為對應的位密),TSW、STUW、AAIW分別表示熱帶表層水、副熱帶下層水和南極中層水.顏色表示緯度.Fig.2 Historical CTD data of the study area, from NCEI(a) Temperature profile; (b) Salinity profile; (c) Sound velocity profile, the blue line is the fitted sound velocity, and (d) T-S diagram (gray contours are potential density). TSW, STUW, AAIW represent Tropical Surface Water, Subtropical Underwater and Antarctic Intermediate Water, respectively. The color indicates the latitude. 地震數據來自于海洋地球科學數據系統(Marine Geoscience Data System,MGDS),選用航次為EW0005.該航次地震數據原為用作勘查尼加拉瓜岸外俯沖帶構造的多道地震數據.這里主要選用大致平行海岸線的測線L115.此外,在分析渦旋的結構時,還使用了7條垂直海岸線的地震測線.8條地震測線的采集日期與時間如表1所示.地震采集參數如下:震源BOLT氣槍(體積6817 in3),炮間距50 m,最小偏移距188 m,每炮240道接收,道間距25 m(故CMP間距為12.5 m),采樣間隔4 ms,記錄長度4 s. 本文所使用的地震數據處理流程(Yilmaz, 2001)包括:觀測系統定義,去除低頻干擾,直達波壓制,CMP道集分選,速度分析,動校正,CMP水平疊加,疊后去噪.由于研究對象是海水,其聲波速度一般在1540~1480 m·s-1范圍內,相對變化較小,故這里使用1500 m·s-1聲速進行常速動校正.此外,海水層反射同相軸傾角一般小于1°,偏移效果不明顯,故這里直接使用疊加剖面而不進行任何偏移處理.一般疊加剖面上存在低頻噪聲,此外還有頻域濾波難以去除的涌浪噪聲(地震剖面上表現為垂向強干擾),適當使用F-K域傾角濾波可以有效壓制這些噪聲.具體流程可以參考地震海洋學數據處理流程(Ruddick et al., 2009). 表1 地震測線的采集日期與采集時間(UTC,2000年)Table 1 Acquisition date and time of seismic survey lines (UTC, 2000) 除了地震數據,還使用了水文數據.水文數據分為再分析數據和歷史CTD數據.再分析數據由歐盟哥白尼海洋環境監測中心(Copernicus Marine Environment Monitoring Service,CMEMS)提供.該數據為對各種來源(CTD、XBT、ADCP、Argo、Glider、Mooring、遙感等)的觀測資料進行質量控制和同化處理、結合相關物理海洋理論模型整合而成的再分析數據.它不僅包含的參數多,空間范圍廣,而且延伸的時段長,是一個綜合的數據集.這里使用的參數均為日平均數據,包括:海水不同空間位置(經緯度及深度)對應的位溫、緯向流速、經向流速和海表高度等.數據水平網格分辨率為1/12°,垂向深度最大可達5000 m.歷史CTD數據來自美國海洋與大氣管理局國家環境信息中心(NOAA-NCEI).CTD測站位置如圖1a所示,一共50個.由CTD數據獲得區域的溫鹽剖面和T-S圖,并且由溫度和鹽度數據并利用海水熱力學狀態方程(TEOS-10)計算出對應的聲速.擬合聲速數據用于地震剖面的時深轉換. 利用CMEMS的再分析數據的水平流速即可計算出垂直地震剖面的水平流速、地震剖面處的Rossby數以及OW參數. 圖3a為L115地震剖面,剖面上部明顯存在一個凸透鏡狀反射.透鏡狀反射的寬度約為160 km,對應剖面上30~190 km處厚度約500 m,深度從100 m延伸到約600 m,核心深度約200 m,橫向位置120 km.反射體大致左右對稱,左上方同相軸向西傾斜,右上方同相軸向東傾斜.此外,中心表現為強反射,中心的東西兩側表現為弱反射,弱反射外圍又顯示為強反射.上邊界反射強度和坡度均較大,而下邊界為較弱的近水平反射帶.結合前人研究(Biescas et al., 2008;Pinheiro et al., 2010;Song et al., 2011),初步判斷這個透鏡狀反射是一個渦旋.渦旋上部兩側的傾斜反射即是渦旋和周圍水團由于橫向溫度梯度而形成溫度鋒.鋒面處以熱鹽侵入為主,橫向熱鹽侵入將橫向溫度梯度轉化為垂向溫度梯度. 下面結合來自CMEMS再分析數據計算的垂直地震剖面的水平流速、沿地震剖面的Rossby數和OW參數進行分析.再分析數據為2000年6月23日的日平均數據,比地震采集時間相比略晚.圖3b為垂直地震剖面的水平流速.200 m以上深度存在較強的水平流速,且東西兩側的流速方向相反:西側流出,東側流入.流速中心分別在約40 km和180 km處.流速方向相反說明剖面存在水平流速剪切.圖3c為被科里奧利參數f標準化后的相對渦度ζ,即Rossby數Ro.北半球f為正值,若存在渦旋,則正、負Ro分別對應氣旋渦和反氣旋渦.剖面200 m以上深度的偏東側存在Ro正值中心,位置與透鏡狀反射的核心相比略偏東,且最大值Romax>0.5.圖3d為剖面對應的OW參數W,與圖3c的Ro正值中心的相同位置表現為負值中心.W<0,表明渦旋主導該區域動力特征.這些剖面特征初步驗證了上述推斷. 為進一步分析,需考慮剖面周圍的海洋學背景.這里選擇了海水溫度、Rossby數、OW參數以及水平流速矢量在不同深度(100、200、300、500 m)的水平切片(圖4).當深度為100 m時,地震測線剛好穿過一個低溫中心,在相同位置存在Ro正值中心和OW參數負值中心,核心處Ro接近1,表明這里亞中尺度過程已相當發育.并且渦旋周圍的正值OW參數說明這里存在水平剪切,和地震剖面上鋒面處(圖3a)位置吻合,這是渦旋邊緣存在強的拉伸和剪切而導致的亞中尺度鋒面,對于渦旋來說也是渦絲存在的位置.流速矢量也表現出氣旋(逆時針旋轉)的特征.氣旋的中心沒有被地震測線穿過,而是位于后者正西南方約25 km處.當深度為200 m時,四種參數的特征與深度為100 m時基本一致,只是強度都有所減弱,但是低溫中心略微向東南移動.當深度為300 m時,各種參數的強度進一步減弱,且氣旋中心向西北移動,氣旋空間規模明顯變小.當深度為500 m時,各種參數的強度繼續減弱,且氣旋中心繼續沿測線向西北移動,氣旋規模已經相當小. 圖3 (a) L115地震剖面; (b) 垂直剖面的水平流速,紅色和藍色分別表示流入和流出; (c)和(d) 分別為地震剖面處的Rossby數和Okubo-Weiss參數Fig.3 (a) Seismic transect L115; (b) Horizontal velocity perpendicular to the transect, red and blue indicate inflow and outflow respectively; (c) and (d) are Rossby number and Okubo-Weiss parameter at the seismic transect, respectively 上述剖面上的流速特征、正Rossby數,以及不同深度切片上的低溫中心、正Rossby數和流速矢量,綜合表明L115地震測線穿過了一個氣旋渦.不同深度的流速矢量表明垂向上渦中心從上向下逐漸向西北傾斜,并且渦旋的強度隨著深度增加而逐漸減弱. 一般普遍使用Rossby變形半徑或者Rossby數來判別亞中尺度渦.例如,McWilliams (1985)按照渦旋水平半徑小于Rossby變形半徑這一標準,認為地中海渦旋為亞中尺度相干渦.但目前還沒有一個公認的標準將海洋現象定義為亞中尺度過程(Zheng et al., 2020).這里定義Rossby數為0.5~1或者水平半徑小于Rossby變形半徑的渦旋為亞中尺度渦.Rossby變形半徑與緯度、水體高度的特征尺度、密度分布有關.次表層亞中尺度渦又叫做亞中尺度相干渦(McWilliams, 1985).尼加拉瓜岸外的Rossby變形半徑約100 km(Chelton et al., 1998),這里L115剖面顯示的渦旋橫向半徑約80 km.一般海表高度的低值中心和高值中心分別對應氣旋渦和反氣旋渦(Pinheiro et al., 2010; Pingree and Le Cann, 1992).但在海表高度圖上(圖1b),渦旋上方并沒有顯示對應于氣旋渦的低異常,即該渦旋屬于海洋次表層氣旋.此外,該氣旋渦中心對應的Rossby數大于0.5(圖4b),所以該渦旋應該屬于亞中尺度相干渦. 一般認為渦旋中心和渦旋邊緣動力特征不同.渦旋邊緣由于與周圍水團的交換,會存在劇烈的熱鹽侵入,在地震剖面上會表現為強反射特征.而渦旋中心由于水團性質較為均一,故而在地震剖面上表現為弱反射.渦旋的這種中心弱、邊界強的反射特征在以往的研究中普遍存在(Biescas et al., 2008;Pinheiro et al., 2010;Song et al., 2011;Gunn et al., 2018;Gula et al., 2019;Tang et al., 2014a).然而圖3a地震剖面中卻顯示在渦旋核心處顯示出強反射特征,這與以往研究的渦旋特征不同.Biescas等(2008)也曾在地震剖面上發現了一個核心不是弱反射的渦旋,而且從渦旋邊緣到核心有多個相間的同心圓狀強反射帶.他們認為有兩種可能性:地震測線可能沒有穿過均勻混合的渦旋核心,或者渦旋已經存在較長時間,從而周圍水團對核心的橫向熱鹽侵入較強烈.但是L115剖面渦旋反射特征與Biescas等人的卻不同.這里橫向上從渦旋核心到渦旋邊緣依次是強反射、弱反射、強反射,且大致關于核心左右對稱;而在垂向上核心下方存在近水平中等強度反射帶.所以可以排除第二種可能,因為若是渦核已被周圍水團橫向侵入,渦核和邊緣之間不應該存在弱反射.但是渦核下方強反射表明可能是渦旋核心和頂底水團發生了劇烈的雙擴散作用形成雙擴散層,從而引起混合. 圖4 L115地震測線采集期間(2000年6月23日)不同深度(100、200、300、500 m)的(a)海水溫度、(b)Rossby數、(c)Okubo-Weiss參數和(d)水平流速矢量品紅色實線表示地震測線,紅色圓點表示采集起點.需要說明的是不同深度的溫度色標范圍不一樣.Fig.4 (a) Sea water temperature, (b) Rossby number, (c) Okubo-Weiss parameter and (d) horizontal velocity at different depths (100, 200, 300, 500 m) during the acquisition of seismic survey line L115 (June 23, 2000)The magenta line is the seismic survey line, and the red dot is the acquisition starting point.Note that the colorbar ranges of temperature at different depths are distinct. 第一種解釋可能存在.因為圖4d平面流速矢量顯示的渦旋中心和地震測線的確有偏移,即使二者對應時間并不完全吻合.也就是說,地震剖面并未捕獲到渦旋核心.如圖5a所示,渦核1在地震剖面的另一側且未被后者穿過,而2—4層則在地震剖面上反映出來.其中2層即為地震剖面上的強反射核心,可能為頂底的雙擴散作用引起的.3層為2層外側的弱反射層,仍保持有來自核心的水團,但是被2層隔開.4層則為最外側的強反射層,橫向熱鹽侵入和垂向雙擴散都存在.另外可以發現,地震剖面上的強反射核心并不是同心圓狀,而是呈上下相鄰的兩個核心,即為雙核結構(圖5b).兩個核心的寬度皆約40 km,上、下方核心厚度分別約80 m和100 m.這里推斷,由于渦旋內部水團和頂底水團的雙擴散作用,使得內部水團性質不再均一,從而形成了上下相鄰的強反射雙核結構. 再分析水文數據可以用來判定渦旋位置,但由于橫向分辨率(1/12°網格分辨率)限制,無法揭示渦旋的精細結構,高分辨地震數據則彌補了這一缺陷,二者結合則可以分析渦旋位置和內部精細結構. 除了大致平行海岸線的L115測線,EW0005航次還有7條與L115正交的地震測線穿過了該氣旋渦.如圖6所示,根據200 m深度的海水水平流速圈出氣旋渦的中心.按照地震測線采集時間順序(表1)和測線相對位置(圖6)可以看出,前7條垂直海岸線的測線按照從東南到西北順序逐漸往返采集.由地震測線和渦旋中心的位置關系(圖6b—i)可以看出,8條測線中有4條穿過了渦旋核心,分別為L22、L26、L24和L28,剩余4條則只穿過了渦旋邊緣.需要說明的是,除了L115測線采集時間較晚(6月23日),其他7條測線依次在6月1—4日共4天內采集.對于上百公里尺度的海水運動來說,可以認為渦旋4天內的變化極小(Tsuji et al., 2005),故可以利用這些不同的地震剖面分析渦旋不同部位(邊緣、核心)的地震反射特征. 圖7為7條穿過渦旋的地震剖面.根據一般渦旋的反射特征,圈定出地震剖面中的渦旋邊界.其中上邊界均表現為明顯的強反射帶,且渦旋右上邊界為傾斜反射,這是由于渦旋水團和周圍水團發生橫向熱鹽侵入而形成的水平溫度梯度帶,即亞中尺度鋒面,鋒面處海水垂向運動較為強烈,是海洋營養物質垂向輸送的重要通道.而渦旋下邊界反射則不明顯,除了L115剖面和L10剖面的渦旋下邊界分別為中等強度反射和傾斜反射外,其他剖面下邊界均表現為近水平弱反射.這種強反射上邊界和弱反射下邊界的特征可能是由于氣旋冷渦等溫線的上凸特性導致的(Zhang et al., 2014).等溫線上凸會使得渦旋上邊界的垂向溫度梯度大,下邊界的垂向溫度梯度小,從而造成這種反射強度的差異. 前7個剖面上渦旋中心均表現為中等強度反射,局部有強反射,且除L10渦旋中心同相軸較為傾斜外,其余反射同相軸均近水平.根據前人研究,渦旋中心一般表現為弱反射或者同心圓狀中等強度反射.這里地震剖面上的氣旋渦中心卻普遍為一系列近水平的中等強度反射,甚至出現局部強反射,與L115的核心強反射特征也不相同.因為海水地震圖像主要是海水垂向溫度梯度的反映(Ruddick et al., 2009),這些近水平反射說明渦旋內部水團不太均勻發生了部分混合,并且表現出良好的分層性.這表明該氣旋渦可能已經存在了較長時間. 圖5 (a)三維渦旋反射層示意圖.1—4分別表示渦核、垂向為主的強反射雙擴散層、弱反射層、強反射熱鹽侵入層.四邊形表示地震剖面,地震剖面未穿過渦核.(b)L115地震剖面上渦旋的強反射雙核結構.藍色區域C1和C2為強反射所圈定的兩個核心Fig.5 (a) The cartoon of 3D eddy reflection layers.Number 1—4 represent eddy core, strong reflection double diffusion layer, weak reflection layer and strong reflection themohaline intrusion layer, respectively.The quadrilateral represents the seismic transect, which does not pass through the eddy core. (b) Dual-core strong reflection structure of the eddy on seismic transect L115. The blue areas C1 and C2 are the two cores delineated by strong reflection 圖7 7個地震剖面中的渦旋圖像(黑色實線大致圈出渦旋邊界,藍色虛方框為水團邊界,藍色虛橢圓為小渦旋,海底反射被切除)Fig.7 The cyclonic eddy images on 7 seismic transects.The black curves roughly delineate the eddy boundaries, the blue dotted rectangles delineate the water boundaries, and the blue dotted circle is small eddy. Seafloor reflections are muted 從渦旋的視厚度可以看出,前7個剖面中,穿過渦旋核心比穿過渦旋邊緣的渦旋視厚度更大,這是由渦旋整體的凸透鏡形態決定的.但是L115的渦旋視厚度明顯較大,這是因為和前7條測線相比,其采集時間較晚,在向西北移動的過程中渦旋的強度逐漸增強、規模逐漸變大.此外,由其他水文參數也可以看出(未顯示),地震數據采集期間,該氣旋渦緩慢向西北方向移動,其空間尺度逐漸增大、強度(流速大小、Rossby數)逐漸增強. 此外,深度為800~900 m時,部分剖面(L22、L26、L24和L28)存在局部較連續的中等強度近水平反射帶,這可能是渦旋下方的部分小尺度的水團邊界(圖7中的藍色虛線方框).根據深度范圍,結合圖2中歷史CTD數據,推斷可能是上方的副熱帶下層水(STUW)和下方的南極中層水(AAIW)形成的水團邊界.L30剖面的渦旋下方存在一個內部弱反射的小渦旋,結合位置分析,推斷該小渦旋可能是從上方氣旋渦中脫落形成.7個剖面渦旋邊界右上方均表現為近水平中等—強反射,總體一致性較好.Brenes等(2008)分析,沿岸流CRCC存在于CRD和海岸之間,且由于大量渦旋的存在,CRCC向北的流路相當曲折.所以這里認為,渦旋右上方的近水平反射即為氣旋渦西北側向極流去的CRCC.這里CRCC將低緯度的熱帶表層水(TSW)沿著海岸向北輸送,對區域氣候影響較大. 總體來說,與L115剖面正交的7個剖面和前者的渦旋反射特征表現出了較大差異.這一方面體現了該氣旋渦空間結構的不對稱性,也說明了其時間演化特征. 一個在中美洲鸚鵡灣采集的地震剖面L115捕獲到了海水層的一個透鏡狀反射結構.透鏡狀反射的寬度約150 km,厚度約500 m,深度從100 m延伸到約600 m,核心深約200 m.結合和地震采集時間幾乎同期的CMEMS再分析數據中的流速和海水溫度數據以及計算的Rossby數和Okubo-Weiss參數,根據剖面上和不同深度切片上的數據,將這一透鏡狀反射綜合解釋為氣旋渦,且可能是亞中尺度相干渦(SCV).和一般渦旋中心的弱反射特征不同,L115地震剖面上該氣旋渦的中心為強反射,并且進一步發現強反射核心為上下相鄰的雙核結構,推測這是由于地震剖面并未穿過渦旋核心水、剖面上的核心強反射可能只是由于渦旋頂底與核心水外層的水團雙擴散作用造成的. 表2 8個地震剖面上的渦旋反射特征Table 2 Reflection characteristics of the eddy on the eight seismic transects 7條和L115大致正交的地震剖面也捕捉到了這個氣旋渦.這些剖面整體反射特征一致性較好,但和L115剖面差異較大.渦旋內部普遍表現為近水平的中等強度反射,這與以往研究的渦旋特征不同.此外,渦旋上邊界為強反射,下邊界為弱反射,這可能是由于氣旋渦的等溫線上凸,使得渦旋上、下邊界垂向溫度梯度分別較大和較小.前7個剖面和L115剖面的渦旋反射特征的差異表明了該氣旋渦空間結構的不對稱性,也體現了氣旋渦的時間演化特征. 過去地震海洋學對海洋渦旋的研究分布在地中海周邊、南海、南大洋、阿拉斯加灣和墨西哥灣等海域.這里我們又對中美洲西海岸的渦旋進行了地震成像,進一步拓寬了海洋渦旋的地震海洋學研究的海域范圍.不同海域渦旋的反射特征有較大差別,未來可以利用地震海洋學方法對其他海域的各種性質的渦旋進行研究,充分利用地震數據高分辨率全海水深度快速成像的優勢,加深對海洋渦旋結構和性質的認識. 致謝感謝海洋地球科學數據系統(Marine Geoscience Data System,http:∥www.marine-geo.org)提供的反射地震數據,歐盟哥白尼海洋環境監測中心(Copernicus Marine Environment Monitoring Service,http:∥marine.copernicus.eu/services-portfolio/access-to-products/)提供的再分析水文數據,以及美國海洋與大氣管理局國家環境信息中心(NOAA-NCEI,https:∥www.nodc.noaa.gov/OC5/indprod.html)提供的歷史CTD數據.感謝外審專家給出的建設性修改意見.

2 數據與方法
2.1 數據

2.2 方法

3 結果與討論
3.1 亞中尺度相干渦

3.2 核心強反射

3.3 正交剖面反射特征


4 結論與展望
