郭雙喜,魯遠征,岑顯榮,屈玲,Sharon L.Walker ,周生啟*

東北太平洋Explorer Ridge熱液羽狀流位溫濁度異常和物質(zhì)能量通量估算

郭雙喜1,魯遠征1,岑顯榮1,屈玲1,Sharon L.Walker2,周生啟1*

(1.中國科學(xué)院南海海洋研究所 熱帶海洋環(huán)境國家重點實驗室,廣東 廣州510301;2.美國國家海洋和大氣局 太平洋海洋環(huán)境實驗室,美國 西雅圖98115)

深海熱液流體與周圍海水之間存在明顯的物理和化學(xué)差異,通過檢測海水的位溫濁度異常是探測深海熱液活動的重要手段之一。本文采用“海底火山帶項目(Sub marine Ring of Fire 2002)”拖曳式溫鹽深測量儀數(shù)據(jù)資料,研究了東北太平洋Explorer Ridge熱液場的水文特征及物質(zhì)能量通量的釋放。結(jié)果表明Explorer Ridge熱液場熱液羽狀流中性浮力層所在深度范圍約為1 600～1 900 m,距離海底的高度約為200 m,最大位溫、鹽度和濁度異常分別為0.04℃、0.004和0.18 NTU;中性浮力層熱液羽狀流帽呈橢圓結(jié)構(gòu),其長軸與洋中脊線重合,羽狀流帽總面積約為27 k m2;熱液羽狀流在中性層范圍內(nèi)存在明顯的分層現(xiàn)象,通過經(jīng)驗公式計算得到Explorer Ridge熱液場觀測范圍內(nèi)熱液噴口的總的浮力通量為6.19×10-2m4/s3,平均值為2.063×10-2m4/s3;總的體積通量為9.884×10-2m3/s,平均值為3.295×10-2m3/s;總的熱通量為194.9 MW,平均值為64.967 MW。

東北太平洋;Explorer Ridge熱液場;位溫異常;濁度異常;熱通量

1 引言

海底熱液活動(Seafloor Hydrother mal Activity)普遍發(fā)育于全球大洋中脊、弧后盆地的板塊增生帶、板內(nèi)火山和轉(zhuǎn)換斷層等地質(zhì)構(gòu)造環(huán)境,是巖石圈和海洋之間進行能量和物質(zhì)交換的主要通道。海底熱液活動對海洋地質(zhì)、生物和化學(xué)環(huán)境、大洋環(huán)流以及全球氣候都有著重要的影響?；谏詈ｃ@探計劃(DSDP)和大洋鉆探計劃(ODP),科學(xué)家在大洋中脊和弧后盆地擴張脊的板塊邊界等處發(fā)現(xiàn)大量金屬硫化物礦體,這種多金屬礦體的形成被認為與熱液活動密切相關(guān)[1]。同時,觀測發(fā)現(xiàn)深海熱液噴口聚集了眾多特殊的生物群落和生態(tài)系統(tǒng),由于熱液系統(tǒng)獨特的物理環(huán)境,這些生物系統(tǒng)被認為可能和生命的起源有關(guān)[2]。另外,海底熱液活動為海底提供了巨大的能量,據(jù)估計全球海洋熱液活動釋放的總熱量為1012～1013W,這足以推動大洋中層水的循環(huán),甚至在更深層次上影響全球海洋環(huán)境和全球氣候變化[3]。因此,研究海底熱液活動,對于進一步研究海底熱液成礦、熱液噴口生物群落、熱液活動對海洋和氣候環(huán)境的影響以及推進深海調(diào)查技術(shù)的發(fā)展等都有重要的科學(xué)意義。

海底熱液活動的主要特征表現(xiàn)為攜帶大量礦物質(zhì)的高溫?zé)嵋毫鲝暮５谉嵋簣龊跓焽鑷姵?在浮力作用下向上運動數(shù)百米,并大量卷挾周圍海水,當(dāng)熱液密度達到與周圍海水密度相當(dāng)時,開始側(cè)向擴散,最終形成巨大的蘑菇狀熱液羽狀流結(jié)構(gòu)。盡管在上升過程中不斷被稀釋,熱液羽狀流與周圍海水之間仍然存在明顯的物理和化學(xué)差異。到目前為止,人們有多種手段來探測深海熱液活動,如海底采樣[4]、聲學(xué)成像[5-6]、光學(xué)反演[7]等,而通過檢測海水的物理和化學(xué)異常仍是當(dāng)前探測熱液活動最主要和高效的手段之一[8-9]。熱液流體從海底流向海洋,向海洋注入大量的物質(zhì)和能量,據(jù)估計,全球熱液活動向海洋傳送的熱量占地球總散失熱量的25%,其量值與全球火山爆發(fā)釋放的能量相當(dāng)[10-11],因此,合理估算熱液場物質(zhì)和能量通量是理解海底熱液活動對深海物理化學(xué)環(huán)境的乃至對全球氣候影響的基礎(chǔ)。東北太平洋是全球海底熱液活動的高發(fā)區(qū)域,對該海域海底熱液活動的探測一直是海洋學(xué)家研究的熱點。本文主要通過東北太平洋Explorer Ridge熱液場的觀測資料,定量分析其熱液羽狀流物理異常的空間分布,并通過經(jīng)驗公式估算其物質(zhì)和能量通量的釋放。

2 數(shù)據(jù)來源

東太平洋海隆(East Pacific Rise,EPR)位于太平洋東部張裂性板塊的邊界,是全球洋中脊的一部分,并將西面的太平洋板塊與東面的北美洲板塊和科克斯板塊等分隔開。東太平洋海隆北部由于板塊擴展速率的差別而形成Explorer Ridge、Juan de Fuca Ridge和Gorda Rigde 3個海脊,其中Explorer Ridge是東太平洋海隆的開端,位于北美溫哥華島西部200余千米處(圖1a)。

圖1 東北太平洋地形圖(a)和Explorer Ridge海域地形及各拖曳CTD測線示意圖(b),星號表示觀測到的熱液羽狀流核心位置Fig.1 Topographic map of Northeast Pacific Ocean(a)and measuring line of towed CTD(b).The red star is the core position of the measured hydrother mal plume

1984年6月加拿大海洋漁業(yè)局通過載人深潛器 首次發(fā)現(xiàn)Explorer Ridge存在熱液場的證據(jù)[12]。2002年6-8月,由美國和加拿大的海洋學(xué)家聯(lián)合執(zhí)行的海底火山帶項目通過海底聲吶成像系統(tǒng)、自動水下深潛器以及拖曳式溫鹽深測量儀(CTD)等觀測手段再次對Explorer Ridge熱液區(qū)進行了深入觀測。拖曳式CTD總共執(zhí)行了7條測線(圖1b),其中5條測線是沿著(或平行于)洋中脊,另兩條測線橫跨洋中脊。圖2為T02B-01測線拖曳CTD的運行軌跡示意圖,為充分獲取熱液區(qū)水文數(shù)據(jù),拖曳式CTD主要集中在海底以上約500 m范圍進行觀測。CTD主要攜帶溫度、鹽度、濁度和壓強等傳感器,其溫度、鹽度和濁度測量精度分別為0.001℃、0.001和0.001 NTU,空間分辨率為1 m,經(jīng)過光滑處理后為5 m。表1為7條拖曳式CTD測線的具體信息。

圖2 T02B-01測線拖曳CTD在緯度-水深平面的運行軌跡示意圖Fig.2 Moving trajectory of towed CTDin latitude-depth plane for T02B-01 measuring line

表1 拖曳CTD各測線信息Tab.1 Measuring line infor mation of towed CTD

3 水文特征分析

3.1 背景水文特征

拖曳式CTD在下放和回收位置可獲得接近全水深的觀測數(shù)據(jù)。圖3a和b為7條測線起始和終止位置的位溫和鹽度廓線。從圖中可以看出,東北太平洋水體從上往下可分為3個主要部分:混合層(0～30 m),躍層(30～200 m)和下層(200 m以深),總體上呈現(xiàn)為表層高溫低鹽、底層低溫高鹽特征?；旌蠈雍蛙S層易受海表風(fēng)能輸入、太陽輻射、蒸發(fā)降水和淡水輸入等影響,在時間和空間上變化明顯;而底層受海表變化的影響幾乎可以忽略,其變化主要是受底層平流的影響[13]。值得注意的是,在躍層區(qū),位溫隨著深度的增加先急劇減小后增加,在180 m深度處出現(xiàn)局部位溫極大值,這是北太平洋亞北極區(qū)特有的次表層溫度反轉(zhuǎn)現(xiàn)象。北太平洋溫度反轉(zhuǎn)一般被認為是由于冬季海表冷卻、降水盈余以及存在強的鹽躍層而形成的,其隨時間和空間的變化也得到了充分的研究[14-17]。盡管存在溫度反轉(zhuǎn),但由于強鹽躍層的存在,水體仍然表現(xiàn)為穩(wěn)定層結(jié)。圖3c為相應(yīng)的位溫-鹽度曲線,圖中等值線為背景位密分布。從圖中可以看出,表層和次表層位溫和鹽度數(shù)據(jù)較為分散,而下層位溫和鹽度數(shù)據(jù)較為集中,這說明表層和次表層易受海表影響,在時間和空間上變化明顯,而下層受到海表的影響較弱。

圖3 各測線起始和終止位置的位溫-水深(a)和鹽度-水深廓線(b)及相應(yīng)的鹽度-位溫分布(c),其中等值線為背景位密(單位:kg/m3)分布Fig.3 Potential temperature-depth(a),salinity-depth(b)profiles at the starting and end points of all measuring lines,and corresponding salinity-potential temperature(c),the contour lines denote potential density(unit:kg/m3)distribution

3.2 位溫和濁度異常

熱液羽狀流在中性浮力層形成巨大的蘑菇狀結(jié)構(gòu),與周圍背景流體相比具有相對較高的溫度、鹽度和濁度,這為深海熱液場的探尋提供了非常有利的條件。拖曳式CTD在經(jīng)過熱液羽狀流中性浮力層的時候,可探測到明顯的溫度、鹽度和濁度異常。

熱液場位溫(鹽度)異常值一般定義為熱液場位溫(鹽度)觀測值與相同位密層背景海水的位溫(鹽度)的差值,而濁度異常值定義為熱液場濁度觀測值與背景海水濁度的差值[18-20]。圖4為T02B-01測線1 500 m以深位密-位溫、位密-鹽度和位密-濁度數(shù)據(jù)分布。從圖中可以看出,背景海水的位溫、鹽度與位密之間存在明顯的線性關(guān)系,這與前人在深海底層所觀測到的結(jié)果一致[18-20];而背景海水的濁度因不含異常濃度的固體顆粒物而表現(xiàn)為常值。在位密27.602～27.635 kg/m3之間(對應(yīng)的深度1 630～1 840 m)出現(xiàn)了數(shù)據(jù)異常。通過將觀測到的數(shù)據(jù)減去線性擬合所得到的背景溫、鹽和濁度表達式,可計算出相應(yīng)的異常值,如圖5所示,考慮到背景海水本身的擾動和測量噪聲,為了較為準(zhǔn)確地確定熱液羽狀流存在的范圍,以位密小于27.6 kg/m3情況下的位溫、鹽度和濁度的標(biāo)準(zhǔn)差的3倍值(3σ)作為背景噪聲閾值,其值分別為0.008℃、0.001和0.015 NTU(如圖5中的垂直虛線所示),然后以大于3σ的數(shù)據(jù)的范圍來定為存在熱液羽狀流的范圍(如圖5中的水平虛線所示)?？梢钥闯鲈谖幻?7.602 kg/m3與27.637 kg/m3之間存在明顯的溫、鹽和濁度異常值,T02B-01測線最大位溫異常、鹽度異常和濁度異常分別為0.04℃、0.004和0.18 NTU。假設(shè)熱液口噴出的熱液流體溫度為300～400℃,根據(jù)最大位溫異常值0.04℃可估算出熱液羽狀流總共所卷挾的周圍海水與熱液口噴出的熱液流體體積比約為104,這也與Lupton等[21]根據(jù)東北太平洋Explorer Ridge附近的Endeavor Ridge熱液噴口數(shù)據(jù)估算的結(jié)果是一致的。

圖4 T02B-01測線1 500 m以深位密—位溫(a)、位密—鹽度(b)和位密—濁度(c)分布,灰色直線為背景數(shù)據(jù)的線性擬合,兩條虛線之間為出現(xiàn)數(shù)據(jù)異常的區(qū)間Fig.4 Dataset of potential density-potential temperature(a),potential density-salinity(b)and potential density-turbidity(c)below 1 500 mfor T02B-01.The gray solid lines are thelinear fitting and the anomalies appear bet ween the t wo gray dashed lines

圖5 T02B-01測線1 500 m以深位溫異常(a)、鹽度異常(b)和濁度異常(c),垂向虛線為各自的背景噪聲閾值,兩條水平虛線為異常值區(qū)間的邊界Fig.5 The ano malies of potential temperature(a),salinity(b)and turbidity(c)below 1 500 mfor T02B-01.The vertical dashed lines denote corresponding backgroud disturbance thresholds.The anomalies appear bet ween the t wo horizontal dashed lines

3.3 觀測斷面的位溫和濁度異常

為了進一步確定熱液羽狀流的空間位置,可將拖曳CTD的觀測數(shù)據(jù)進行克里金(Kriging)插值,從而得到測線斷面的觀測數(shù)據(jù)。克里金插值法是基于自相關(guān)統(tǒng)計模型對周圍觀測值進行加權(quán),以得出未測量位置的預(yù)測值,是廣泛應(yīng)用于地學(xué)領(lǐng)域的一種最優(yōu)內(nèi)插法[22]。由于熱液噴口所噴出的高溫流體攜帶大量固體顆粒物,因此容易在海底以上200～300 m的熱液羽狀流中性浮力層通過CTD檢測到溫度和濁度異常;但另一方面,美國Alvin深潛器在Explorer Ridge附近的Endeavor Ridge熱液噴口采樣發(fā)現(xiàn)噴出流體的鹽度與周圍海水并沒有明顯的差異[21,23],Lupton等[21]將中性浮力層所觀測到的鹽度異常歸結(jié)為熱液流體在上升過程中卷挾周圍的高鹽度水并攜帶至中性浮力層所致;同時,通過圖5可以看出鹽度異常的量值并不如溫度和濁度更為明顯。因此,我們主要關(guān)注位溫和濁度異常值的斷面結(jié)果。由圖5a,5c已知背景海水的位溫和濁度波動約為0.008℃和0.015 NTU,因此我們將位溫異常大于0.008℃或者濁度異常大于0.015 NTU的區(qū)域認為是熱液羽狀流所在的位置。

圖6至圖8為各測線所在斷面的位溫和濁度異常等值線圖。由于T02B-01、T02B-02和T02B-03均沿著Explorer Ridge的中脊線且在同一直線上,因此將其顯示在同一個斷面上,如圖6所示,可以看到T02B-01測線斷面出現(xiàn)了明顯的位溫和濁度異常,即在此斷面上存在著熱液羽狀流;而T02B-02和T02B-03測線斷面沒有出現(xiàn)明顯的位溫和濁度異常。T02B-01測線斷面羽狀流寬度從49.73°N延伸至49.8°N(約10.1 k m距離),深度范圍從1 630 m至1 840 m。T02B-06和T02B-09測線橫跨Explorer Ridge中脊線,從圖7中可以看出,T02B-09測線斷面沒有檢測到明顯的位溫和濁度異常,而T02B-06測線檢測到明顯的位溫和濁度異常,羽狀流寬度從130.23°W延伸至130.27°W(約3.8 k m距離),深度范圍從1 650 m至1 900 m,其羽狀流寬度明顯小于T02B-01測線斷面測得的羽狀流寬度,而深度范圍相當(dāng)。T02B-10和T02B-11測線為位于T02B-01測線東西兩側(cè)且與T02B-01平行的兩條斷面,如圖8所示,T02B-11測線斷面沒有檢測到明顯的位溫和濁度異常,而T02B-10測線檢測到明顯的位溫和濁度異常,羽狀流寬度從49.73°N延伸至49.77°N(49.77°N以北仍存在羽狀流,但超出了CTD測線范圍),深度范圍從1 620 m至1 900 m,同樣與T02B-01測線斷面測得的羽狀流深度范圍相當(dāng)。從圖6至圖8中可以看到熱液羽狀流所在位置的形狀不對稱,這有可能是因為存在多個熱液噴口或者深層平流的緣故。在深海弱層結(jié)條件下,忽略背景流的影響,熱液羽狀流在中性浮力層主要通過分子擴散形式向四周水平擴散,而熱擴散系數(shù)大于羽狀流攜帶懸浮顆粒物的擴散系數(shù),這可能是造成圖6至圖8通過位溫異常得到的羽狀流面積相對較大、而濁度異常得到的羽狀流面積更為集中且異常值相對更大的原因之一。

圖6 T02B-01、T02B-02和T02B-03測線緯度-水深斷面位溫異常(a)和濁度異常(b),虛線為熱液羽狀流的核心位置Fig.6 Potential temperature anomaly(a)and turbidity ano maly(b)for T02B-01,T02B-02 and T02B-03 sections in latitude-depth plane.The dashed line denotes the core position of hydrother mal plu me

圖7 T02B-06(a,b)和T02B-09(c,d)測線經(jīng)度-水深斷面位溫異常 (a,c)和濁度異常(b,d)Fig.7 Potential temperature(a,c)and turbidity(b,d)anomalies for T02B-06(a,b)and T02B-09(c,d)sections in longitude-depth plane

圖8 T02B-10(a,b)T02B-11(c,d)測線緯度-水深斷面位溫異常(a,c)和濁度異常(b,d)Fig.8 Potential temperature(a,c)and turbidity(b,d)anomalies for T02B-10(a,b)and T02B-11(c,d)sections in latitude-depth plane

4 熱液羽狀流帽與通量估算

4.1 熱液羽狀流帽

熱液柱在上升過程中卷挾周圍海水使得自身溫度迅速降低,當(dāng)上升到中性浮力層后向四周水平擴散,形成具有巨大的羽狀流帽(plu me hat)的蘑菇狀羽狀流結(jié)構(gòu),羽狀流帽中的溫度、鹽度和濁度等物理特性相對于周圍背景海水有明顯異常,因而容易被CTD等儀器所檢測到。根據(jù)T02B-01、T02B-06和T02B-10 3條測線及其位溫異常所檢測到的羽狀流位置(圖9中黑色實線),假設(shè)羽狀流帽為橢圓結(jié)構(gòu),3條羽狀流位置線恰好可構(gòu)成一個規(guī)則的橢圓,其長軸沿著T02B-01測線(即洋中脊線),長度約為10.1 k m,短軸垂直于洋中脊線,長度約為3.4 k m。海底熱液活動區(qū)一般會沿洋中脊線存在多個熱液噴口,而噴口的分布、背景海流的大小和方向等都可能是影響羽狀流帽形狀的重要因素。通過平面幾何計算,容易得到Explorer Ridge熱液場本航次觀測范圍內(nèi)中性浮力層熱液羽狀流帽面積約為27 k m2。

圖9 T02B-1、T02B-6和T02B-10測線及其檢測到的熱液羽狀流位置Fig.9 Detected plume area for T02B-1,T02B-6 and T02B-10灰色實線為3條測線位置,黑色實線為各測線上檢測到羽狀流的位置,3條黑色線構(gòu)成規(guī)則的橢圓形羽狀流帽(由于T02B-10測線起始位置所限,故用虛線將其羽狀流線補齊)The gray solid lines denote t he measuring lines,the black solid lines denote the detected plu me sections,and the ellipse denotes the plu me hat(t he plu me line of T02B-10 is complemented with dashed one due to the li mitation of t he initial observation position)

4.2 通量估算

在不考慮深海平流的情況下,熱液羽狀流能夠上升的高度取決于熱液噴口流體的浮力通量B0和背景海水的浮力頻率N,前者體現(xiàn)了熱液口流體自身所具有的能量,后者表現(xiàn)為層結(jié)海水對羽狀流上升的阻礙作用。Turner[24]提出了靜止海水中熱液羽狀流上升的最大高度Zmax與初始浮力通量B0以及浮力頻率N之間的標(biāo)度關(guān)系為:

式中,浮力頻率N=[(-g/ρ)(dρ/d z)]1/2,浮力通量B0的表達式為:

式中,g、ρ、α和β分別為流體重力加速度、密度、熱膨脹系數(shù)和鹽度收縮系數(shù);Δθ0、ΔS0以及V0分別為熱液噴口處熱液流體與周圍海水之間的溫度差、鹽度差以及噴口流體的體積通量。

由式(1)和式(2)可以得到浮力通量B0和體積通量V0為:

熱液噴口流體熱通量H0表達式為:

式中,ρ0為噴口流體密度;Cp為流體比熱。熱液區(qū)流體物性參數(shù)為g=9.8 m/s2,α=-2.13×10-4℃-1,Cp=6 400 J/(kg·℃)[25]。根據(jù)東北太平洋熱液場熱液噴口的溫度和鹽度特性,取Δθ0=300℃,ΔS0=0[23]。通過羽狀流所在高度范圍的背景海水的密度數(shù)據(jù)可計算得到浮力頻率N=1.171×10-3s-1。因此,只需要知道羽狀流的最大高度Zmax,即可通過式(3)至式(5)得到噴口流體的浮力通量B0、體積通量V0以及熱通量H0。

取圖6中虛線(T02B-01測線)所示羽狀流核心位置的位溫和濁度異常廓線,其位置為49.785 1°N,130.263 3°W,深度為1 887 m,如圖10所示。從圖中可以看出位溫異常和濁度異常廓線在深度為1 663 m、1 681 m和1 719 m位置存在明顯不同的位溫和濁度異常峰值,從而可將中性層羽狀流厚度范圍分成3個不同的羽狀流層Layer1、Layer2和Layer3。Rudnicki和Elderfield[25]在大西洋中脊的TAG熱液場同樣發(fā)現(xiàn)了熱液羽狀流的分層現(xiàn)象,不同層內(nèi)流體的水文性質(zhì)存在明顯差別。Rudnicki和Elderfield[25]將其歸因于3種可能的解釋:(1)羽狀流內(nèi)部存在對流現(xiàn)象;(2)在不同高度處羽狀流的水平擴散率不同;(3)不同的羽狀流層來自不同的熱液噴口。Rudnicki和Elderfield[25]認為第三種解釋可能更為合理,因為在深海弱混合條件下,不同噴口噴出的熱液流體類似于獨立的水團,在向上發(fā)展過程中能夠保持相對獨立的物理特性。

假設(shè)羽狀流層Layer1、Layer2和Layer3來自于3個不同的熱液噴口,由圖5已知背景海水的位溫和濁度擾動分別為0.008℃和0.015 NTU(圖10中垂向虛線),以濁度異常為例,將背景濁度擾動與真實濁度擾動廓線的上下交點視為Layer3的上邊界和Layer1的下邊界,Layer1的上邊界(即Layer2的下邊界)和Layer2的上邊界(即Layer3的下邊界)由廓線的局部最小值確定,從而得到各羽狀流層的上下邊界。由于無法確定各個噴口的具體位置,因此將各羽狀流層上邊界到該廓線所在位置的海底的垂向距離近似視為各羽狀流層的最大高度,從而得到各羽狀流層的最大高度Zmax分別為193 m、223 m和249 m。再通過式(3)至式(5)可以計算3個不同噴口流體的浮力通量B0、體積通量V0以及熱通量H0,其結(jié)果如表2所示,浮力通量B0值分別為1.115×10-2m4/s3、1.987×10-2m4/s3和3.088×10-2m4/s3,總的浮力通量為6.19×10-2m4/s3,平均值為2.063×10-2m4/s3;體積通量V0值分別為1.779×10-2m3/s、3.174×10-2m3/s和4.931×10-2m3/s,總的體積通量為9.884×10-2m3/s,平均值為3.295×10-2m3/s;熱通量H0值分別為35.098 MW、62.559 MW和97.243 MW,總的熱通量為194.9 MW,平均值為64.967 MW。

圖10 羽狀流核心位置(圖6中虛線)的位溫異常(a)和濁度異常(b)廓線Fig.10 Profiles of potential temperature(a)and turbidity(b)anomalies at the detected plu me core(dashed line in Fig.6)垂向虛線(0.008℃和0.015 NTU)分別為背景海水的位溫和濁度擾動,橫向虛線為各羽狀流層的邊界線The vertical dashed lines respectively denote background disturbances(0.008℃and 0.015 NTU),the horizontal dashed lines denote the boundaries for plumelayers

表2 通量計算結(jié)果Tab.2 The estimated fluxes of hydrother mal vents for layers

5 結(jié)論和討論

通過對東北太平洋Explorer Ridge熱液場拖曳CTD數(shù)據(jù)進行位溫和濁度異常分析,并對噴口物質(zhì)能量通量進行估算,得到如下結(jié)論:(1)Explorer Ridge熱液場熱液羽狀流中性浮力層所檢測到的最大位溫、鹽度和濁度異常分別為0.04℃、0.004和0.18 NTU,熱液羽狀流總共所卷挾的周圍海水與熱液口噴出的熱液流體體積比約為104。熱液羽狀流中性浮力層所在深度范圍約為1 600～1 900 m,其核心位置距離海底的高度約為200 m。(2)對檢測到的熱液羽狀流位置進行幾何分析發(fā)現(xiàn)中性浮力層羽狀流帽呈橢圓結(jié)構(gòu),其長軸與洋中脊線重合,長度約為10.1 k m,短軸垂直于洋中脊線,長度約為3.4 k m,從而得到Explorer Ridge熱液場本航次觀測范圍內(nèi)中性浮力層熱液羽狀流帽面積約為27 k m2。(3)發(fā)現(xiàn)熱液羽狀流在中性層范圍內(nèi)存在明顯的分層現(xiàn)象并分析其可能的形成機制。通過經(jīng)驗公式結(jié)合觀測到的最大羽狀流高度計算得到Explorer Ridge熱液場觀測范圍內(nèi)熱液噴口的總的浮力通量為6.19×10-2m4/s3,平均值為2.063×10-2m4/s3;總的體積通量為9.884×10-2m3/s,平均值為3.295×10-2m3/s;總的熱通量為194.9 MW,平均值為64.967 MW。

盡管到目前為止Explorer Ridge熱液場觀測資料仍然較為缺乏,但Explorer Ridge附近同屬東太平洋海隆開端的Juan de Fuca Ridge(JDFR,位置見圖1)熱液場積累了相對充分的觀測數(shù)據(jù)。已有的觀測資料表明JDFR熱液中性浮力層深度約為1 900 m,溫度鹽度異常分別為0.02～0.05℃和0.005[21,26-27];39個獨立羽狀流噴口的總熱通量約為534 MW[28],和本文得到的Explorer Ridge熱液場數(shù)據(jù)相當(dāng),這說明在這兩個相鄰的熱液區(qū)熱液羽狀流具有相似的結(jié)構(gòu)和物理特征。另外,熱液羽狀流到達中性浮力層后其水平擴散過程可能會受到底層流的影響,Baker[29]通過長時間觀測發(fā)現(xiàn)JDFR熱液羽狀流在跨洋脊方向近似為對稱結(jié)構(gòu),并指出熱液場底層流主要沿著洋脊方向;而Thomson等[27]觀測發(fā)現(xiàn)JDFR熱液場沿洋脊方向的典型流速為v=0.35 c m/s?；诖宋覀兛梢越乒浪鉋xplorer Ridge熱液場本航次觀測區(qū)域內(nèi)中性浮力層與背景海水之間的熱通量Hn=r Anv Cpd(式中,密度r=103kg/m3,跨洋脊方向羽狀流縱向截面積An=3.4 k m×200 m,比熱Cp=3 900 J/(kg·℃),溫度異常d=0.04℃),得到Hn為371.3 MW,與經(jīng)驗公式估算得到的噴口的總熱通量也較為接近,說明本航次觀測區(qū)域的熱液活動處于穩(wěn)定地活躍狀態(tài)。

致謝:感謝美國國家海洋和大氣局(NOAA)提供的拖曳CTD觀測數(shù)據(jù)。

[1] 孫樞.大洋鉆探與中國地球科學(xué)[J].地球科學(xué)進展,1995,10(3):213-214.Sun Shu.Ocean drilling and earth science in China[J].Advance in Earth Sciences,1995,10(3):213-214.

[2] Reysenbach A L,Cady S L.Microbiology of ancient and modern hydrot her mal systems[J].Trends in Microbiology,2001,9(2):79-86.

[3] Elderfield H,Schultz A.Mid-ocean ridge hydrother mal fluxes and the chemical co mposition of the ocean[J].Annual Review of Earth and Planetary Sciences,1996,24(1):191-224.

[4] Hah m D,Baker E T,Siek Rhee T,et al.First hydrother mal discoveries on the Australian-Antarctic Ridge:Discharge sites,plu me chemistry,and vent organisms[J].Geochemistry,Geophysics,Geosystems,2015,16(9):3061-3075.

[5] Xu G,Jackson D R,Bemis K G,et al.Ti me-series measurement of hydrother mal heat flux at the Grotto mound,Endeavour Segment,Juan de Fuca Ridge[J].Earth and Planetary Science Letters,2014,404:220-231.

[6] Bemis K G,Silver D,Xu G,et al.The pat h to COVIS:A review of acoustic i maging of hydrother mal flow regi mes[J].Deep-Sea Research PartⅡ:Topical Studiesin Oceanography,2015,121:159-176.

[7] Mittelstaedt E,Escartín J,Gracias N,et al.Quantifying diffuse and discrete venting at the Tour Eiffel vent site,Lucky Strike hydrother mal field[J].Geochemistry,Geophysics,Geosystems,2012,13(4):Q04008.

[8] 欒錫武,秦蘊珊.現(xiàn)代海底熱液活動的調(diào)查研究方法[J].地球物理學(xué)進展,2002,17(4):592-597.Luan Xiwu,Qin Yunshan.Survey methods of modern hydrother mal activity[J].Progressin Geophysics,2002,17(4):592-597.

[9] Baker E T.Exploring the ocean for hydrother mal venting:New techniques,new discoveries,newinsights[J].Ore Geology Reviews,2017,86:55-69.

[10] Lowell R P,Rona P A,Von Herzen R P.Seafloor hydrother mal systems[J].Jour nal of Geophysical Research:Solid Earth,1995,100(B1):327-352.

[11] 翟世奎,王興濤,于增慧,等.現(xiàn)代海底熱液活動的熱和物質(zhì)通量估算[J].海洋學(xué)報,2005,27(2):115-121.Zhai Shikui,Wang Xingtao,Yu Zenghui,et al.Heat and mass flux esti mation of modern seafloor hydrother mal activity[J].Haiyang Xuebao,2005,27(2):115-121.

[12] Tunnicliffe V,Botros M,De Burgh M E,et al.Hydrother mal vents of Explorer ridge,northeast Pacific[J].Deep-Sea Research Part A:Oceanographic Research Papers,1986,33(3):401-412.

[13] Uda M.Oceanography of the subarctic Pacific Ocean[J].Journal of the Fisheries Board of Canada,1963,20(1):119-179.

[14] Ueno H,Yasuda I.Distribution and for mation of the mesother mal structure(temperature inversions)in the North Pacific subarctic region[J].Journal of Geophysical Research:Oceans,2000,105(C7):16885-16897.

[15] Ueno H,Oka E,Suga T,et al.Seasonal and interannual variability of temperatureinversions in the subarctic North Pacific[J].Geophysical Research Letters,2005,32(20):2444.

[16] Ueno H,Yasuda I.Temperature inversions in the subarctic North Pacific[J].Journal of Physical Oceanography,2005,35(12):2444-2456.

[17] Masuda S,Awaji T,Sugiura N,et al.Interannual variability of temperature inversions in the subarctic Nort h Pacific[J].Geophysical Research Letters,2006,33(24):194-199.

[18] Baker E T,Ger man CR,Elderfield H.Hydrother mal plumes over spreading-center axes:Global distributions and geological inferences[M]//Sea-floor Hydrother mal Systems:Physical,Chemical,Biological,and Geological Interactions.Washington,D.C.:American Geophysical Union,1995:47-71.

[19] 王曉媛,武力,曾志剛,等.海底熱液柱溫度異常自動化計算方法探討[J].海洋學(xué)報,2012,33(2):185-191.Wang Xiaoyuan,Wu Li,Zeng Zhigang,et al.Automatic calculation on t hetemperature anomaly of marine hydrother mal plu me[J].Haiyang Xuebao,2012,33(2):185-191.

[20] 陳小丹,梁楚進,董昌明.西南印度洋龍旂熱液區(qū)羽流信號的檢測與通量估算[J].海洋學(xué)研究,2015,33(4):43-52.Chen Xiaodan,Liang Chujin,Dong Changming.Detection and flux esti mation of hydrother mal plu mes in the Longqi hydrother mal field in the Southwest Indian Ocean[J].Journal of Marine Sciences,2015,33(4):43-52.

[21] Lupton J E,Delaney J R,Johnson H P,et al.Entrain ment and vertical transport of deep-ocean water by buoyant hydrother mal plu mes[J].Nature,1985,316(6029):621-623.

[22] Oliver M A,Webster R.Kriging:a method of interpolation for geographical infor mation systems[J].International Journal of Geographical Information System,1990,4(3):313-332.

[23] Mc Duff R E,Lupton J E,Kadko D,et al.Chemistry of hydrother mal fluids,Endeavor Ridge,northeast Pacific[J].Eos Trans.AGU,1984,65:975.

[24] Turner J S.Buoyancy Effects in Fluids[M].Cambridge:Cambridge University Press,1979.

[25] Rudnicki M D,Elderfield H.Theory applied to t he Mid-Atlantic Ridge hydrother mal plu mes:The finite-difference approach[J].Journal of Volcanology and Geot her mal Research,1992,50(1/2):161-172.

[26] Baker E T,Massoth GJ.Characteristics of hydrother mal plumes fromt wo vent fields on the Juan de Fuca Ridge,northeast Pacific Ocean[J].Earth and Planetary Science Letters,1987,85(1/3):59-73.

[27] Thomson R E,Delaney J R,Mc Duff R E,et al.Physical characteristics of the Endeavour Ridge hydrother mal plu me during July 1988[J].Eart h and Planetary Science Letters,1992,111(1):141-154.

[28] Baker E T,Massoth GJ,Walker S L,et al.A method for quantitatively esti mating diffuse and discrete hydrother mal discharge[J].Earth and Planetary Science Letters,1993,118(1/4):235-249.

[29] Baker E T.A 6-year ti me series of hydrother mal plu mes over the Cleft segment of the Juan de Fuca Ridge[J].Journal of Geophysical Research:Solid Earth,1994,99(B3):4889-4904.

Temperature and turbidity anomalies and flux estimation of hydrother mal plume in Explorer Ridge in the Northeast Pacific Ocean

Guo Shuangxi1,Lu Yuanzheng1,Cen Xianrong1,Qu Ling1,Sharon L.Walker2,Zhou Shengqi1

(1.State Key Labor ator y of Tropical Oceanography,South China Sea Institute of Oceanology,Chinese Academy of Sciences,Guangzhou 510301,China;2.Pacif ic Marine Environ mental Labor ator y,NOAA,Seattle,WA 98115,USA)

There are significant physical and chemical differences bet ween deep-sea hydrother mal fluids and ambient seawater.Detecting thetemperature and turbidity anomaly is one of the main methods to explore the hydrother mal fields.Based on the towed CTDdataset of the project of Sub marine Ring of Fire 2002,the hydrological characteristics and flux discharge of hydrother mal vents in Explorer Ridge were analyzed.The neutral buoyancy layer in Explorer Ridge hydrother mal filed is 1 600 to 1 900 m,about 200 m above the seafloor.The maxi mu ms of temperature,salinity and turbidity anomaly are 0.04℃,0.004 and 0.18 NTU,respectively.The shape of the neutral buoyancy layer is approxi mately elliptical,the major axis overlaps with the ocean ridge,and the total area is about 27 k m2.It is found that the neutral buoyancy layer is distinctly layered.It is esti mated with empirical for mulas that the total buoyancy flux,volu me flux and heat flux of the measured hydrother mal vents are 6.19×10-2m4/s3,9.884×10-2m3/s and 194.9 MW,respectively,and their average values are 2.063×10-2m4/s3,3.295×10-2m3/s and 64.967 MW,respectively.

the Northeast Pacific Ocean;Explorer Ridge;potential temperature anomaly;potential turbidity anomaly;heat flux

P738.6

A

0253-4193(2017)12-0001-11

郭雙喜,魯遠征,岑顯榮,等.東北太平洋Explorer Ridge熱液羽狀流位溫濁度異常和物質(zhì)能量通量估算[J].海洋學(xué)報,2017,39(12):1-11,

10.3969/j.issn.0253-4193.2017.12.001

Guo Shuangxi,Lu Yuanzheng,Cen Xianrong,et al.Temperature and turbidity ano malies and flux esti mation of hydrother mal plu me in Explorer Ridge in the Northeast Pacific Ocean[J].Haiyang Xuebao,2017,39(12):1-11,doi:10.3969/j.issn.0253-4193.2017.12.001

2017-03-03;

2017-05-15。

國家自然科學(xué)基金(41406035,41476167,41606010);中科院戰(zhàn)略性先導(dǎo)專項(XDA11030301);廣東省自然科學(xué)基金(2016 A030311042,2016 A030313155)。

郭雙喜(1983—),男,湖北省武漢市人,博士,從事海洋中小尺度動力過程的研究。E-mail:sxguo@scsio.ac.cn

*通信作者:周生啟,男,研究員,主要從事深海動力過程,海洋小尺度過程及海洋混合和輸運研究。E-mail:sqzhou@scsio.ac.cn