鄭靜靜，劉桂梅，高　姍（.國家海洋環(huán)境預報中心，北京0008；2.廈門大學近海海洋環(huán)境科學國家重點實驗室，福建廈門36005）

海洋缺氧現(xiàn)象的研究進展

鄭靜靜1，2，劉桂梅1，高姍1
（1.國家海洋環(huán)境預報中心，北京100081；2.廈門大學近海海洋環(huán)境科學國家重點實驗室，福建廈門361005）

綜述海洋缺氧研究的國內(nèi)外現(xiàn)狀，缺氧形成生物化學機理，以及從物理方面分別分析徑流沖淡水、上升流、風、氣候變暖等因素對缺氧影響；并分別從箱式模型、一維模型、二維模型、三維模式這4個角度，探討了缺氧數(shù)值模型的發(fā)展和現(xiàn)狀。基于目前取得的主要研究進展，提出未來我國缺氧研究發(fā)展重點和方向。

缺氧；機理；模型

1　引言

海洋中的溶解氧是重要的生源要素參數(shù)，也是海洋生態(tài)系統(tǒng)得以維持發(fā)展的關鍵因子。近幾十年來，受人類活動的影響，大量的污染物排入近岸海域，造成水體富營養(yǎng)化逐年加劇，致使的近岸底層水體缺氧現(xiàn)象也呈不斷上升趨勢。通常定義水體中的溶解氧（DO）濃度＜2.0 mg/L為缺氧狀態(tài)［1］，當水體中溶解氧的濃度＜2.0 mg/L時，海洋中大部分水生生物將面臨死亡，海底拖曳無法捕捉到魚蝦種群。

水體缺氧嚴重危害河口、近海環(huán)境，表現(xiàn)在多個方面。首先，海洋中大多數(shù)生物都需要溶解氧來維持。如魚類生長需要6 mg/L的溶解氧，蝦、蟹生長所需的溶解氧為2—3.5 mg/L［2］。缺氧對海洋中底棲生物的影響最為顯著，當溶解氧濃度小于2 mg/L，且持續(xù)時間較長時，海洋中的大多數(shù)魚類、浮游動物，特別是運動能力較弱的底棲生物群落將面臨大規(guī)模死亡［3］。缺氧事件發(fā)生還會降低海洋物種多樣性，改變海洋生物的群落結(jié)構(gòu)，減少魚類和底棲動物的豐富度，從而影響漁業(yè)生產(chǎn)，帶來直接或間接的經(jīng)濟損失［4］。如波羅的海的“死亡區(qū)”因為長時間持續(xù)缺氧，每年損失的碳有2.64×105t，占整個波羅的?？偝跫壣a(chǎn)力的30%，并造成漁業(yè)減產(chǎn)總量為1.06×105t［5］。

隨著全球缺氧區(qū)域的擴張，釋放越來越多的溫室氣體，如NO2、H2S，進而對全球氣候產(chǎn)生潛在影響［6］。Naqvi等對全球缺氧區(qū)每年向大氣釋放的N2O和CH4總量進行估算，結(jié)果顯示N2O和CH4分別為1.48—3.11 Tg和0.36—0.72 Tg［7］。

海洋缺氧事件不但影響了海洋生物化學過程，而且也會改變?nèi)虻奶嫉h(huán)，從而對整個生態(tài)系統(tǒng)產(chǎn)生嚴重的危害［8］。海洋的缺氧現(xiàn)象已經(jīng)在全球范圍內(nèi)愈演愈烈，引起科學家們的高度重視，缺氧形成機制及成因成為當今海洋領域研究的熱點之一。文章綜述了海洋缺氧的現(xiàn)狀和形成機理，以及利用模型來研究海洋缺氧形成過程，并提出現(xiàn)階段缺氧研究的側(cè)重點。

2　缺氧的現(xiàn)狀

近幾十年來，由于人為排污加劇，發(fā)生缺氧的海域數(shù)量以每年5.54%的指數(shù)速率迅速增長［9］。另外缺氧發(fā)生的頻率、范圍和持續(xù)時間均有顯著增加。如西大西洋最大的缺氧區(qū)出現(xiàn)在墨西哥灣北部，密西西比河入?？谔帲?002年的缺氧面積達到2.2×104km2，而1993到2009年缺氧區(qū)域的平均面積為1.6×104km2［10］。Conley等發(fā)現(xiàn)波羅的海最早出現(xiàn)缺氧區(qū)是20世紀60年代，其1991—2000年缺氧區(qū)年均面積為4.9×104km2［11］。地中海海灣從1987年開始出現(xiàn)缺氧［12］，近幾年缺氧現(xiàn)象愈發(fā)嚴重，導致大量海洋生物死亡［13］。也有學者研究發(fā)現(xiàn)秘魯近海［14］、美國新澤西洲海岸［15］、弗吉尼亞州的切薩皮克灣［16］、卡羅萊納州沿岸河流（大西洋西海岸）［17］、韓國南部海灣［18］均出現(xiàn)不同程度的缺氧現(xiàn)象。這些海域的缺氧現(xiàn)象已經(jīng)不是海洋生態(tài)系統(tǒng)循環(huán)自然發(fā)生的生態(tài)現(xiàn)象，人類活動造成的環(huán)境污染是導致缺氧區(qū)形成的一個重要原因。

我國珠江口外和長江口外海域也存在底層水體明顯的缺氧現(xiàn)象。1999夏季的調(diào)查顯示長江口外存在一處面積高達1.37×104km2的缺氧區(qū)域，氧虧損總量高達1.59×106t［19］。20世紀50年代，長江口外發(fā)生夏季缺氧事件的頻率為60%，而到了1990年后缺氧事件發(fā)生頻率達到90%，并且缺氧面積大于5.0×103km2的基本上都發(fā)生在20世紀90年代末，這說明近幾十年來長江口缺氧區(qū)范圍不斷擴大，嚴重程度也是逐年加劇。夏季，在我國珠江口外海域亦發(fā)現(xiàn)有底層水體缺氧現(xiàn)象，1985年的調(diào)查發(fā)現(xiàn)珠江口缺氧區(qū)主要位于橫琴島和高欄島附近水域，底層溶解氧含量最低為1.76 mg/L，此后，珠江口缺氧區(qū)呈現(xiàn)不斷擴大化、嚴重化趨勢［20］。王麗芳等于2005年1月和8月的兩次調(diào)查發(fā)現(xiàn)，珠江口上游廣州河段的底層溶解氧濃度平均為3.09 mg/L和1.21 mg/L，在廣州附近水域底層溶解氧最低值分別為0.48 mg/L和0.21 mg/L［21］。Helly等研究了近岸海域缺氧的全球分布，結(jié)果表明近岸缺氧區(qū)的面積高達1.0×106km2［22］。這些缺氧區(qū)主要分布在西歐、北歐沿海，美國東、西部海岸，中國東部沿海，以及日本沿海等人口密集、經(jīng)濟發(fā)達的海岸帶區(qū)域，這也進一步印證了人類活動對缺氧區(qū)的形成有著不可推卸的責任。

3　缺氧形成機理

缺氧實際上就是溶解氧消耗大于補給過程的持續(xù)。要弄清海洋中缺氧的問題，首先需要知道海洋中溶解氧的源和匯。

海洋中溶解氧的源，主要有兩方面：大氣復氧和浮游植物光合作用。大氣復氧是大氣中的氧通過海氣界面進入到海洋中的過程，是海水溶解氧重要來源之一。當海水中的溶解氧不飽和時，氧氣就會從大氣進入到海洋，當海水中溶解氧為過飽和狀態(tài)時，水中的氧氣就會釋放到空氣中。目前廣泛使用的水體大氣復氧方程見式（1）：

式中：k2為復氧系數(shù)，Os為溶解氧飽和濃度，O為實際溶解氧濃度。這個公式表明大氣對水體的復氧速率與溶解氧的虧損（Os-O）成正比［23］。此外，海水溫度、紊動（如風浪、環(huán)流）等對海氣交換也有一定影響。水體表面?zhèn)髻|(zhì)系數(shù)與溫度的關系滿足式（2）：

式中：KT為溫度T時的傳質(zhì)系數(shù)，K20為20℃時的表面?zhèn)髻|(zhì)系數(shù)，θ為常數(shù)，T為溫度［24］。由此可見海水溫度越高，越有利于氧氣與海水之間的氣體交換。然而在海氣交換過程中，溫度越高，氧氣的溶解度越低，這樣會部分抵消溫度對復氧速率的影響。另外溫度升高，溶解氧飽和度會降低，將不斷從海洋中釋放氧。紊動（風浪、環(huán)流）不僅影響海氣界面氧氣傳輸，而且對水體內(nèi)部溶解氧的擴散、輸運也有著重要作用。海洋中溶解氧的另一種來源是浮游植物光合作用，其光合作用強弱和光照強度有關，隨著深度的增加，透光層以下的光合作用急劇減小，因此光合作用主要發(fā)生在上面的透光層區(qū)。海洋中浮游植物的光合作用過程可由式（3）的正過程表示：

海洋中溶解氧的匯，主要有生物呼吸作用，有機物的分解耗氧，無機物的氧化作用，底泥耗氧等過程。浮游植物的呼吸作用是光合作用的逆過程，在補償深度（呼吸耗氧量等于光合產(chǎn)氧量）以下，特別是在透光層以下，由于光線減弱，浮游植物的呼吸作用占主導地位，消耗溶解氧。有機物的降解過程也消耗大量的溶解氧，氧的消耗量主要取決于有機物的含量。河口地區(qū)有機物來源豐富，含量高，其分解耗氧量多，所以容易出現(xiàn)缺氧現(xiàn)象。研究表明珠江口外海域底層水體缺氧主要是由有機物分解消耗了底層溶解氧導致的，尤其是密度躍層以上水域［25］。此外海洋中的一些還原態(tài)無機物如Fe2+、Mn2+在氧氣充足的情況下會發(fā)生氧化反應，轉(zhuǎn)化為更穩(wěn)定的高價態(tài)，此過程氧消耗量很小。另外，NH4+通過硝化作用耗氧，珠江口缺氧區(qū)形成很大一部分是由硝化反應耗氧貢獻的。據(jù)估計，珠江口上游缺氧區(qū)硝化作用的耗氧量占總耗氧量的20%—30%［26］。海洋中的底泥是指沉積到海底的顆粒態(tài)有機物，其經(jīng)礦化作用消耗溶解氧。底泥耗氧在水體耗氧中也有著非常重要的貢獻，在珠江口密度躍層以下水域，底泥耗氧占據(jù)主導地位［25］。

3.1生物化學機制

缺氧現(xiàn)象的形成是個復雜的過程，是物理和生物化學共同作用的結(jié)果，受溫度、鹽度、水體層化、浮游植物的生物量以及有機物等諸多因素的影響。河口近岸區(qū)域的缺氧研究均指出河口水域，夏季溫度適宜，河流攜帶的大量氮、磷營養(yǎng)鹽促使的浮游植物爆發(fā)性繁殖，初級生產(chǎn)力增加，浮游植物通過光合作用產(chǎn)生的氧氣也隨之增加，與此同時，大量繁殖的浮游植物通過呼吸作用的耗氧量也對應增加。除了呼吸作用耗氧外，有機物的氧化分解是最為主要的生物化學耗氧過程。有機物的降解耗氧可以用下面的經(jīng)驗方程（4）表達：

河口中有機物的來源主要有浮游植物死亡的貢獻和河流輸入。夏季，一方面河流徑流量增加，河流輸入的有機物也隨之增加；另一方面浮游植物初級生產(chǎn)力提高，浮游植物死亡貢獻的有機物也相應增加。以上兩種有機物進入河口后，溶解態(tài)有機物會直接在水體中分解消耗溶解氧。顆粒態(tài)有機物中有一部分在細菌等微生物的作用下轉(zhuǎn)化為溶解態(tài)有機物消耗溶解氧，另一難分解的部分在輸運過程中沉降到底層，發(fā)生復雜的礦化反應，大量消耗溶解氧。雖然目前已經(jīng)明確有機物降解耗氧是水體缺氧的一個重要因素，但是關于導致河口缺氧的有機物來源一直備受爭論。Jia等通過使用放射性同位素δ13C對珠江口伶仃洋中的沉積物柱狀樣進行研究，發(fā)現(xiàn)伶仃洋內(nèi)河流輸入的陸源有機物對底泥中有機物的貢獻遠大于浮游植物死亡后對底泥有機物的貢獻［27］。但是也有研究認為河流輸入的有機物對低氧的貢獻微小。Dagg等認為河流輸入的陸源有機物中，只有極少一部分的有機物是能被細菌利用分解，而來自海源的有機物（即浮游植物殘體），在短時間內(nèi)細菌能對其大部分進行分解。實驗表明，河流輸入的有機物，細菌在經(jīng)過620 d的分解后，僅有34%被利用；而相對于海源有機物來講，細菌僅分解3 d就利用了58%的有機物［28］。越來越多研究表明缺氧的形成與浮游植物的初級生產(chǎn)力有著直接的關系。Kemp等對切薩皮克灣的總有機碳收支研究發(fā)現(xiàn)，陸源輸入的有機物遠小于海源輸入的有機物，缺氧形成主要是由海源有機物分解消耗溶解氧導致［29］。亞得里亞海北部，波羅的海和黑海，其缺氧程度的增加與河流輸入的N通量、初級生產(chǎn)力的增加直接相關［30］。

近幾十年來，隨著人為排污加劇，尤其是營養(yǎng)鹽負荷增加的近岸海域，缺氧事件頻發(fā)，很多研究認為水體富營養(yǎng)化是造成水體缺氧重要的“原始驅(qū)動力”。從直接影響來看，富營養(yǎng)化的近海為浮游植物的爆發(fā)性繁殖提供了極其有利的物質(zhì)基礎，這些生物在生長過程中通過呼吸作用會消耗水體中的溶解氧，而其死亡后，殘體會沉降到水體底層被細菌分解，也會消耗大量的溶解氧，加劇低氧狀態(tài)；從間接影響來看，營養(yǎng)鹽比率的變化，使浮游植物的種群結(jié)構(gòu)發(fā)生改變，新演替出的優(yōu)勢種由于一般不受捕食者喜愛而不被捕食，最終下沉到底層被細菌分解耗氧。另外水體底層缺氧又會誘發(fā)底泥磷的釋放，釋放的磷上升至透光層促進浮游植物生長，形成一種正反饋效應，進一步加劇低氧趨勢。水體富營養(yǎng)化導致浮游植物大量繁殖為缺氧的形成提供了物質(zhì)基礎，如果物理條件成熟，使水體出現(xiàn)層化，阻止溶解氧的垂向交換，底層溶解氧得不到表層溶解氧的補充，水體底層就會出現(xiàn)缺氧現(xiàn)象。同時如果水動力條件能夠維持較長的水體停留時間，那么水體缺氧就會進一步惡化，導致生態(tài)系統(tǒng)崩潰。我國長江口缺氧區(qū)形成的生化過程，主要是由于長江徑流帶來的大量陸源營養(yǎng)鹽，促進表層浮游植物生長，產(chǎn)生的大量顆粒態(tài)有機碳向底層輸送，底層溶解氧的消耗主要包括陸源有機物，浮游植物死亡貢獻以及底泥再懸?。?1］。切薩皮克灣浮游植物生長主要受N限制，Testa等研究發(fā)現(xiàn)，當?shù)斎朐黾訒r，浮游植物初級生產(chǎn)力增加，植物殘體死亡降解以及呼吸作用消耗大量氧氣，發(fā)生缺氧的區(qū)域擴張，而當1—5月氮輸入減少時，7—8月缺氧區(qū)域也對應減?。?2］。Justi?等人利用數(shù)值模型研究墨西哥灣北部水體底層缺氧對密西西比河營養(yǎng)鹽輸入變化的響應，結(jié)果表明如果密西西比河營養(yǎng)鹽通量減少30%，缺氧發(fā)生頻率也隨之減少37%［33］。

3.2物理機制

缺氧現(xiàn)象除了受生物和化學過程影響之外，物理因素在缺氧現(xiàn)象的形成、維持和破壞中也起著重要作用。徑流沖淡水、上升流、風、氣候變暖對海洋缺氧現(xiàn)象有著重要的貢獻。

3.2.1徑流沖淡水

徑流沖淡水不但可以攜帶大量陸源有機物，而且可以在河口區(qū)形成溫鹽躍層，限制表底層的水體交換。夏季大量徑流沖淡水向外擴散，高溫低鹽的淡水浮在低溫高鹽的海水上面，咸淡水交界處有著明顯的密度差異，一方面形成溫鹽躍層，阻礙了表層溶解氧的補充，另一方面形成鋒面，如同一堵墻，阻礙營養(yǎng)物質(zhì)和顆粒物質(zhì)的運輸，使其在鋒面所處位置附近大量沉降。如夏季長江口外海的高溫、低鹽的表層長江沖淡水和底層北上的臺灣暖流高鹽、低溫水形成較強的溫度和鹽度躍層成為垂向溶解氧交換的屏障［34］。近幾十年來美國墨西哥北部灣缺氧現(xiàn)象的越來越嚴重，這與密西西比河、阿查法拉亞河近50a來的營養(yǎng)鹽通量成數(shù)倍增長有著密切的關系［35］。Djakovac等研究發(fā)現(xiàn)亞得里亞海北部自1990年后發(fā)生大規(guī)模缺氧事件的頻率明顯減少，這與波河輸入營養(yǎng)鹽減少有關［36］。

3.2.2上升流

陸架淺海上升流，一方面攜帶外源性低溶解氧入侵，使溶解氧背景值相對偏低，更易發(fā)生缺氧現(xiàn)象；另一方面使底層富含營養(yǎng)鹽的水體上升，加劇富營養(yǎng)化并促使藻華爆發(fā)，植物殘體沉降從而增加底層有機碎屑，其分解需要消耗更多溶解氧。由美國新澤西海岸的西南風和適宜的海底地形共同作用而引起上升流，能夠持續(xù)幾個星期，而由上升流帶來的顆粒態(tài)有機物降解會消耗大量溶解氧，從而致使氧虧損［15］。春夏季，在強烈的南風以及西南季風作用下，秘魯-智利陸架處的高NO3、低氧的赤道次表層水涌升，使得真光層富營養(yǎng)化，浮游植物生長旺盛，植物呼吸作用增強消耗溶解氧，隨后出現(xiàn)缺氧現(xiàn)象［37］。在墨西哥灣北部陸架區(qū)，風引起的上升流加劇了缺氧現(xiàn)象，F(xiàn)eng等利用模型研究發(fā)現(xiàn)由于上升流導致的高初級生產(chǎn)力在底層水體消耗大量的溶解氧［38］。

3.2.3風

風場強弱也影響缺氧現(xiàn)象的發(fā)生和發(fā)展，風速增加有利于增強水體的垂向混合過程，打破水體層化，促進溶解氧的垂向交換，破壞水體底層缺氧的形成。風向影響河口環(huán)流和水體滯留時間，進而影響缺氧水體發(fā)生的位置。美國北卡羅萊納州的Pamlico河口在弱風速條件下，水體混合弱，層化穩(wěn)定持續(xù)時間長，阻止溶解氧的垂向交換，因此更易出現(xiàn)缺氧現(xiàn)象［39］。臺風可以增強水體垂向混合過程，使得表層的溶解氧能有效向水體底層補充，從而使河口水體底層缺氧現(xiàn)象消失；臺風帶來強降雨使得臺風過后的河口與近海水體出現(xiàn)鹽度躍層，同時也使得陸源營養(yǎng)鹽和顆粒有機物增加，浮游植物爆發(fā)，從而加劇低氧事件發(fā)生。墨西哥灣北部密西西比河徑流沖淡水通量直接影響營養(yǎng)鹽通量大小和水體層化程度，臺風過境的頻率與強度也會影響墨西哥灣北部缺氧現(xiàn)象［40］。在路易斯安娜州西北部的德克薩斯陸架處，熱帶風暴使得底層的沉積物再懸浮，降解消耗大量的溶解氧，Justi?等通過數(shù)值模擬發(fā)現(xiàn)，熱帶風暴（Fay）過后，德克薩斯陸架底層的溶解氧消耗速度達到0.5 mg/d，一周后開始出現(xiàn)缺氧現(xiàn)象［41］。

3.2.4氣候變暖

全球變暖，導致海水溫度升高，溫度與氧的溶解度呈顯著負相關，溫度越高，氧的溶解度越低，從而使得通過海氣交換進入到缺氧區(qū)表層的溶解氧降低，限制了缺氧區(qū)溶解氧的有效補充。此外溫度升高，生物的呼吸作用增強，消耗氧氣。Bendtsen等對波羅的海缺氧進行數(shù)值模擬發(fā)現(xiàn)當溫度升高3℃，氧氣的溶解度降低，其對低氧的貢獻為25%，同時生物的呼吸增強，導致發(fā)生缺氧的區(qū)域擴大并且缺氧發(fā)作的時間延長［42］。另一方面海水溫度升高，使得溫躍層更加穩(wěn)定，水體垂直混合減弱，不利于氧的交換，底層海水消耗的溶解氧得不到上層及時的補充，從而導致缺氧區(qū)面積不斷擴大。此外，Shaffer等研究表明隨著大洋溫度的不斷升高，誘發(fā)更多甲烷從底泥中釋放，這樣會形成一種正反饋效應，推進全球變暖趨勢，進一步加劇缺氧惡化趨勢［43］。同時氣候變暖導致冰川融化，海平面上升，進而導致底層水鹽度的上升［44］，從而增加水體的層化程度，降低了底部水體的垂直混合。未來全球變暖將進一步加速缺氧區(qū)域的擴張，尤其是溫帶至副極地區(qū)域。Keeling等通過海洋模型研究預測表明，到下個世紀，全球溶解氧的含量將下降1%—7%，并且溶解氧的含量在未來的一千年里將持續(xù)下降［45］。

總之，生化過程和物理過程共同控制著海洋中溶解氧的變化，缺氧現(xiàn)象是這兩大過程相互作用的結(jié)果。生化過程控制著溶解氧的產(chǎn)生和消耗，主要包括光合作用產(chǎn)氧、大氣復氧，有機物分解耗氧、硝化作用耗氧、底泥耗氧。物理過程控制著溶解氧的水平和垂直輸運，水平輸運主要受沿岸流、鋒面、潮流等過程影響，垂直輸運主要包括上升流，水體層化等過程?？v觀這些缺氧事件，盡管具體發(fā)生時間、地點、生物地理環(huán)境不同，但相互之間卻有共同的特征，海洋水體底層缺氧是由水體層化和有機物的生物氧化過程中大量耗氧共同作用形成的。

4　缺氧模型研究

目前對于缺氧的研究多是基于水文和生化要素調(diào)查資料分析的方法，但是此類方法存在一定局限性，如對缺氧空間范圍和時間變化認識的全面性受限于調(diào)查資料，同時難以獲取缺氧發(fā)生的物理機制。為了對水體缺氧的形成機理及其影響因素有更深入的認識，在進行長期、廣泛的現(xiàn)場監(jiān)測和實驗研究的基礎上，還要使用各種海洋數(shù)值模式來模擬和預測水體缺氧的機理、發(fā)展、影響等相關問題。目前，科學家已開發(fā)出了許多數(shù)值模式來研究海洋缺氧問題。伴隨著海洋生物，物理海洋和化學海洋學的快速發(fā)展，并結(jié)合利用計算機技術的新成果，海洋生態(tài)動力學模型從零維模型逐漸發(fā)展到三維模型。

4.1箱式模型

圖1　箱式模型示意圖

箱式模型，又稱零維模型。主要以研究各變量隨時間的變化為目的，把所研究區(qū)域在空間上劃分成一個或幾個均勻的箱子。箱式模型的優(yōu)點是簡便易行，但也存在諸多不足之處，如不適用于動力機制研究，另外空間分辨率也較低等。Das等利用箱式模型來研究墨西哥灣中北部的巴拉塔里亞河口的碳預算以及缺氧現(xiàn)象［46］。他將研究的區(qū)域劃分為6個箱子，見圖1。每個箱子的質(zhì)量平衡方程見式（5）：

式中：Vi指每個箱子的體積，F(xiàn)i指由于水位變化引起的水體積的輸入及輸出，Pi指降水，Ri指鄰近濕地的河流輸入，Ei指蒸發(fā)，Qi指密西西比河的河流輸入。Das的研究結(jié)果表明巴拉塔里亞河口的碳輸出對缺氧現(xiàn)象的影響很小。Officer等利用溶解氧箱式收支模型，定量分析了各個因素對缺氧形成的貢獻及其相互之間的關系［47］。Officer等的工作無疑是具有開創(chuàng)性的，但他的箱式模型只考慮了時間的變化，而忽略了空間的變化，考慮的生態(tài)因子，物理過程及生物求化學過程也較簡單。

4.2一維模型

一維模型是將生態(tài)模型與一維水動力模型耦合，但只考慮垂直方向的輸運，而忽略水平輸運引起的生態(tài)變化。這種模型適用于生態(tài)變量水平變化不明顯的海區(qū)或是開闊的大洋區(qū)域，常用于研究生態(tài)變量的年際變化。Vanderborght等利用一維CONTRASTE（Coupled，Networked，Transport-Reaction Algorithm Strong Tidal Estuaries）模式對斯凱爾特河口O2、CO2、N2O等生源氣體的通量進行了計算。CONTRASTE模式包含了潮流，淡水輸入，物理化學及生態(tài)的輸運過程。模式中含有鹽度、懸浮物、無機碳、無機氮、有機碳、有機氮、浮游植物等因子。生態(tài)過程主要描述了異養(yǎng)細菌的呼吸作用，光合作用，硝化和反硝化作用［48］。CONTRASTE模式中變量的平衡方程見式（6）：

式中：Cj表示生態(tài)變量的濃度，A表示斷面的面積，Q表示通過斷面的流量，K表示擴散系數(shù)，Rj是指由生態(tài)過程引起的變量Cj的變化。Hofmann等通過建立一維模型來研究斯凱爾特河口的氮和碳收支，并估算了溶解氧的幾個主要源匯過程，結(jié)果表明硝化作用消耗大量的溶解氧，是造成斯凱爾特河口底層水體缺氧的主要生化過程［49］。

4.3二維模型

二維模型是指生態(tài)模型與二維水動力模型耦合，只考慮垂直方向和一個水平方向上水動力或兩個水平方向上水動力輸運引起的狀態(tài)變量的變化。二維模式方面，Valle-Levinson等用了一個二維混合模型來研究美國西部長島海峽缺氧的機理，模式中溶解氧的平衡方程見式（7）：

式中：Avo表示溶解氧的垂直擴散率，R表示溶解氧的產(chǎn)生量與消耗量的差值，溶解氧的瞬時變化受物理輸運和生態(tài)過程共同影響。模擬結(jié)果揭示了生化反應，層化作用和被削弱的重力環(huán)流以及垂直混合作用的減弱致使水體底層出現(xiàn)缺氧現(xiàn)象，該模型，詳盡的分析了缺氧形成的物理機理，然而對于生化耗氧只做簡單的數(shù)值估計［50］。Ishikawa等也利用一個二維模型模擬了日本利根川（Tone River）河口水體停留時間和鹽楔入侵等因素對水體缺氧形成的影響［51］。

4.4三維模型

三維模型更注重物理過程和生化過程相互作用下生態(tài)系統(tǒng)的變化，用來模擬生態(tài)系統(tǒng)在三維空間的分布特征，更加符合生態(tài)系統(tǒng)的實際情況。耦合了物理與生態(tài)過程的海洋生態(tài)動力學模型中，生態(tài)相關狀態(tài)變量的濃度變化除了受到物理輸運過程的控制，還受到生物、化學作用的影響，見式（8）：

式中：v→代表三維的流速，kT代表混合系數(shù)，C代表生態(tài)相關的狀態(tài)變量。第一項為對流項，表示由于水體流動造成的狀態(tài)變量的水平和垂向的輸運。第二項為混合項，表示由于湍流作用造成的水平與垂向擴散。第三項Bc代表生態(tài)過程對C濃度的變化率。Fennel等在區(qū)域海洋模式系統(tǒng)（Regional Ocean Modeling System，ROMS）基礎上，耦合生態(tài)模式來研究不同的模型精度，物理邊界條件，以及底泥耗氧過程對墨西哥灣北部缺氧區(qū)域的影響［52］。Fennel生態(tài)模式中影響溶解氧的過程有海氣交換，浮游植物光合作用產(chǎn)氧，硝化作用耗氧，浮游動物新代謝耗氧，碎屑礦化分解耗氧以及海底沉積物再礦化耗氧，如圖2所示。研究結(jié)果發(fā)現(xiàn)墨西哥灣北部缺氧區(qū)域的模擬對不同底泥耗氧過程很敏感，另外改變模型物理邊界條件可以通過影響水體層化進而影響缺氧區(qū)域的模擬。Sohma等把河口生態(tài)模型應用到日本的Mikawa海灣中，在Sohma的模型中應用了3個模型，分別是新海洋生態(tài)模型，水力模型，潮汐水質(zhì)模型，模擬系統(tǒng)能夠較好揭示夏季Mikawa海灣的水體動力學以及水體中的碳、氮、磷和氧的生態(tài)系統(tǒng)動力學，從理論上說明缺氧現(xiàn)象形成的物理原因和生化原因［53］。羅琳等使用三維斜壓水動力—生態(tài)耦合模式模擬分析了珠江口底層水體缺氧現(xiàn)象的形成及其原因，研究結(jié)果表明過量的營養(yǎng)鹽、有機質(zhì)的輸入以及徑流沖淡水引起的強化是珠江口底層水體缺氧的主要原因，而由潮汐和鋒面提供的垂向混合控制著低氧水團的范圍和強度［54］。Justi?利用非結(jié)構(gòu)化網(wǎng)格的三維水動力模型（Finite Volume Coastal Ocean Model，F(xiàn)VCOM）耦合經(jīng)修正水質(zhì)分析模擬程序模擬了德克薩斯州陸架處的缺氧的時空變化及其缺氧形成的機制［41］。生態(tài)動力模型能夠用于表達缺氧發(fā)生的物理及生化過程，因此缺氧區(qū)的生態(tài)動力學模型開發(fā)是未來發(fā)展的重點。

圖2　Fennel模式的缺氧形成過程

5　缺氧現(xiàn)象的研究展望

海洋缺氧區(qū)的形成是一系列自然和人為因素共同作用的結(jié)果，缺氧的存在和缺氧程度的惡化，必然會嚴重威脅生態(tài)系統(tǒng)的健康，對沿海地區(qū)的水產(chǎn)養(yǎng)殖產(chǎn)業(yè)造成重大的經(jīng)濟損失。缺氧的研究涉及的學科眾多，覆蓋了生物、化學和物理等諸多學科，它的發(fā)展應強調(diào)學科間的交叉、滲透與綜合。目前，國外在缺氧形成機制，生態(tài)效應，影響因素以及防治等方面的研究更加成熟，國內(nèi)對缺氧的研究比較落后，各個方面的技術有待提高?？v觀我國缺氧的發(fā)展歷史和研究現(xiàn)狀，其發(fā)展趨勢應在以下幾個方面有所側(cè)重：

（1）定量分析物理及生物化學過程對缺氧現(xiàn)象的影響。河口水體底層的缺氧現(xiàn)象是物理過程和生物化學過程綜合作用的結(jié)果，物理過程調(diào)控著生物化學過程，這兩者都是缺氧形成的必要條件，缺一不可。但在不同河口及近岸海域，二者的貢獻有所不同。定量估算各物理過程、生化過程對溶解氧收支的貢獻，找出導致缺氧的主導因素，有助于加深對缺氧形成機制的理解，為未來制定有效的治理方案提供科學依據(jù)；

（2）建立缺氧的生態(tài)系統(tǒng)動力學模型，開展缺氧的數(shù)值預報系統(tǒng)，提高缺氧區(qū)災害的預警能力。雖然近幾年來我國海洋生態(tài)動力數(shù)值模型研究取得了較大的進展，但由于受海洋觀測調(diào)查資料、海洋環(huán)境認知不足等因素的限制，至今尚未開展業(yè)務化的缺氧生態(tài)動力學數(shù)值預報系統(tǒng)。而美國在數(shù)值預報方面已取得一定進展，美國通過政府機構(gòu)與科研院所的廣泛合作，針對墨西哥灣和切薩皮克灣兩個缺氧最嚴重的區(qū)域開展監(jiān)測，預報和研究，并于2002年開始了墨西哥灣的缺氧區(qū)的數(shù)值預報，利用硝酸鹽的排放負荷來預測當年的缺氧區(qū)規(guī)模，預報結(jié)果較為準確。我國需要加強缺氧區(qū)水動力和生化要素的現(xiàn)場觀測，通過對現(xiàn)場監(jiān)測和實驗研究的綜合分析來確定模型所需生態(tài)參數(shù)、初始場、邊界條件等，反過來再通過生態(tài)動力學模型來驗證參數(shù)的可行性與適用性，從而完善和改進生態(tài)動力學模式，為缺氧預警提供科學依據(jù)；

（3）全球變暖與缺氧關系。我國沿海海洋環(huán)境要素變化，如海表溫度升高、海平面變化等，與全球氣候變化有著密不可分的聯(lián)系。全球變暖引起的海水溫度升高，導致水體層化加強、氧的溶解度降低、生物呼吸作用增強等進而促進近岸海域缺氧現(xiàn)象呈惡化趨勢。目前全球氣候變化與缺氧的研究還處于比較初級階段，缺氧對全球氣候變化響應的內(nèi)在機理尚不清晰，相關的研究內(nèi)容主要局限于氣候變化通過營養(yǎng)鹽和徑流量的改變影響缺氧過程。未來隨著全球氣溫上升，缺氧區(qū)域?qū)⑦M一步擴張，因此需要加強物理生態(tài)耦合過程模型研發(fā)，深入研究氣候變化對缺氧的影響，提高預測未來全球變暖下缺氧發(fā)展趨勢的能力。

近年來我國近海針對缺氧區(qū)的調(diào)查和監(jiān)測，往往只是針對個別的研究課題，缺乏系統(tǒng)性和連續(xù)性的研究，對缺氧的情況尚未有一個全面深刻的認識，不同部門之間的數(shù)據(jù)共享機制尚未建立。美國的海洋缺氧區(qū)研究工作，有一部《赤潮和缺氧研究控制法案》作為總的原則，事實上不止這項工作，在很多領域我們可以看到美國都是法律和戰(zhàn)略規(guī)劃先行。我國可以借鑒這種思路，在海洋缺氧區(qū)研究領域上，站在戰(zhàn)略高度上制定全盤規(guī)劃，制定恰當?shù)臓款^部門，然后引入多部門多機構(gòu)的合作，共同解決我國近海缺氧問題。

［1］Wei H，He Y C，Li Q J，et al.Summer hypoxia adjacent to the Changjiang Estuary［J］.Journal of Marine Systems，2007，67（3-4）: 292-303.

［2］Gray J S，Wu R S S，Or Y Y.Effects of hypoxia and organic enrichment on the coastal marine environment［J］.Marine EcologyProgress Series，2002，238:249-279.

［3］Karlson K，Rosenberg R，Bonsdorff E，et al.Temporal and spatial large-scale effects of eutrophication and oxygen deficiency on benthic fauna in Scandinavian and Baltic waters-a review［J］. Oceanography and Marine Biology，2002，40:427-489.

［4］Yin K D，Lin Z F，Ke Z Y.Temporal and spatial distribution of dissolved oxygen in the Pearl River Estuary and adjacent coastal waters［J］.Continental Shelf Research，2004，24（16）:1935-1948.

［5］Diaz R J，Rosenberg R.Spreading dead zones and consequences for marine ecosystems［J］.Science，2008，321（5891）:926-929.

［6］Wang F F，Liu J，Qiu J D，et al.Historical evolution of hypoxia in the East China Sea off the Changjiang（Yangtze River）estuary for the last～13，000 years:evidence from the benthic foraminiferal community［J］.Continental Shelf Research，2014，90:151-162.

［7］Naqvi S W A，Bange H W，F(xiàn)arías L，et al.Marine hypoxia/anoxia as a source of CH4and N2O［J］.Biogeosciences，2010，7（7）: 2159-2190.

［8］Lam P，Lavik G，Jensen M M，et al.Revising the nitrogen cycle in the Peruvian oxygen minimum zone［J］.Proceedings of the National Academy of Sciences of the United States of America，2009，106（12）:4752-4757.

［9］Vaquer-Sunyer R，Duarte C M.Thresholds of hypoxia for marine biodiversity［J］.Proceedings of the National Academy of Sciences of the United States ofAmerica，2008，105（40）:15452-15457.

［10］Obenour D R，Scavia D，Rabalais N N，et al.Retrospective analysis of midsummer hypoxic area and volume in the northern GulfofMexico，1985-2011［J］.EnvironmentalScience& Technology，2013，47（17）:9808-9815.

［11］Conley D，Bj?rck S，Bonsdorff E，et al.Hypoxia-related processes in the baltic sea［J］.Environmental Science&Technology，2009，43（10）:3412-3420.

［12］Friligos N，Psilidou R，Xatzigewrgiou E，et al.Seasonal variations on nutrients and dissolved oxygen［M］//Tsiavos C.Oceanographic Study of the Amvrakikos Gulf.Chemical Oceanography Final Report 3.Athens，Greece:Hellenic Centre of Marine Research（HCMR），1989（in Greek）.

［13］Kountoura K，Zacharias I.Temporal and spatial distribution of hypoxic/seasonal anoxic zone in Amvrakikos Gulf，Western Greece［J］.Estuarine，Coastal and Shelf Science，2011，94（2）: 123-128.

［14］Neira C，Sellanes J，Levin L A，et al.Meiofaunal distributions on the Peru margin:relationship to oxygen and organic matter availability［J］.DeepSeaResearchPartI:Oceanographic Research Papers，2001，48（11）:2453-2472.

［15］Glenn S，Arnone R，Bergmann T，et al.Biogeochemical impact of summertime coastal upwelling on the New Jersey Shelf［J］. Journal of Geophysical Research:Oceans，2004，109（C12）: C12S02.

［16］Murphy R R，Kemp W M，Ball W P.Long-term trends in chesapeake bay seasonal hypoxia，stratification，and nutrient loading［J］.Estuaries and Coasts，2011，34（6）:1293-1309.

［17］Mallin M A，Johnson V L，Ensign S H，et al.Factors contributing to hypoxia in rivers，lakes，and streams［J］.Limnology and Oceanography，2006，51（1）:690-701.

［18］Lim H S，Diaz R J，Hong J S，et al.Hypoxia and benthic communityrecoveryinKoreancoastalwaters［J］.Marine Pollution Bulletin，2006，52（11）:1517-1526.

［19］李道季，張經(jīng)，黃大吉，等.長江口外氧的虧損［J］.中國科學D輯，2002，32（8）:686-694.

［20］林洪瑛，劉勝，韓舞鷹.珠江口底層海水季節(jié)性缺氧現(xiàn)象及其引發(fā)CTB的潛在威脅［J］.湛江海洋大學學報，2001，21（S1）: 25-29.

［21］王麗芳，戴民漢，翟惟東.近岸、河口缺氧區(qū)域的主要生物地球化學耗氧過程［J］.廈門大學學報（自然科學版），2007，46（S1）: 33-37.

［22］Helly J J，Levin L A.Global distribution of naturally occurring marine hypoxia on continental margins［J］.Deep Sea Research Part I:Oceanographic Research Papers，2004，51（9）:1159-1168.

［23］雒文生，李莉紅，賀濤.水體大氣復氧理論和復氧系數(shù)研究進展與展望［J］.水利學報，2003，34（11）:64-70.

［24］李然，李嘉，李克鋒.紊動水體大氣復氧系數(shù)研究進展［J］.上海環(huán)境科學，2000，19（4）:176-179.

［25］張恒，李適宇.生化過程對夏季珠江口底層缺氧影響的模擬研究［C］//中國環(huán)境科學學會2009年學術年會論文集.武漢:中國環(huán)境科學學會，2009:376-382.

［26］王麗芳.珠江口水域硝化作用與無機氮分布特征［D］.廈門:廈門大學，2007.

［27］Jia G D，Peng P A.Temporal and spatial variations in signatures of sedimented organic matter in Lingding Bay（Pearl estuary），southern China［J］.Marine Chemistry，2003，82（1-2）:47-54.

［28］Dagg M J，Ammerman J W，Amon R M W，et al.A review of water column processes influencing hypoxia in the northern Gulf of Mexico［J］.Estuaries and Coasts，2007，30（5）:735-752.

［29］Kemp W M，Smith E M，Marvin-DiPasquale M，et al.Organic carbon balance and net ecosystem metabolism in Chesapeake Bay［J］.Marine Ecology Progress Series，1997，150:229-248.

［30］Cloern J E.Our evolving conceptual model of the coastal eutrophication problem［J］.Marine Ecology Progress Series，2001，210:223-253.

［31］劉海霞，李道季，高磊，等.長江口夏季低氧區(qū)形成及加劇的成因分析［J］.海洋科學進展，2012，30（2）:186-197.

［32］Testa J M，Li Y，Lee Y J，et al.Quantifying the effects of nutrient loading on dissolved O2cycling and hypoxia in Chesapeake Bay using a coupled hydrodynamic-biogeochemical model［J］.Journal of Marine Systems，2014，139:139-158.

［33］Justi? D，Rabalais N N，Turner R E.Simulated responses of the GulfofMexicohypoxiatovariationsinclimateand anthropogenic nutrient loading［J］.Journal of Marine Systems，2003，42（3-4）:115-126.

［34］Li D J，Zhang J，Huang D J，et al.Oxygen depletion off the Changjiang（Yangtze River）Estuary［J］.Science in China Series D:Earth Sciences，2002，45（12）:1137-1146.

［35］Rabalais N N，Turner R E，Wiseman Jr W J.Gulf of Mexico hypoxia，A.K.A."The dead zone"［J］.Annual Review of Ecology and Systematics，2002，33:235-263.

［36］Djakovac T，Supi? N，Aubry F B，et al.Mechanisms of hypoxia frequency changes in the northern Adriatic Sea during the period 1972-2012［J］.Journal of Marine Systems，2015，141:179-189.

［37］Farías L，Cornejo M.Effect of seasonal changes in bottom water oxygenation on sediment N oxides and N2O cycling in the coastal upwelling regime off central Chile（36.5°S）［J］.Progress in Oceanography，2007，75（3）:561-575.

［38］Feng Y，F(xiàn)ennel K，Jackson G A，et al.A model study of the response of hypoxia to upwelling-favorable wind on the northern Gulf of Mexico shelf［J］.Journal of Marine Systems，2014，131: 63-73.

［39］Lin J，Xu H Z，Cudaback C，et al.Inter-annual variability of hypoxic conditions in a shallow estuary［J］.Journal of Marine Systems，2008，73（1-2）:169-184.

［40］Go?i M，Gordon E S，Monacci N M，et al.The effect of Hurricane Lili on the distribution of organic matter along the inner Louisiana shelf（Gulf of Mexico，USA）［J］.Continental Shelf Research，2006，26（17-18）:2260-2280.

［41］Justi? D，Wang L X.Assessing temporal and spatial variability of hypoxia over the inner Louisiana-upper Texas shelf:application of an unstructured-grid three-dimensional coupled hydrodynamicwater quality model［J］.Continental Shelf Research，2014，72: 163-179.

［42］Bendtsen J，Hansen J L S.Effects of global warming on hypoxia in the Baltic Sea-North Sea transition zone［J］.Ecological Modelling，2013，264:17-26.

［43］Shaffer G，Olsen S M，Pedersen J O P.Long-term ocean oxygen depletion in response to carbon dioxide emissions from fossil fuels［J］.Nature Geoscience，2009，2（2）:105-109.

［44］Hilton T W，Najjar R G，Zhong L，et al.Is there a signal of sea-levelriseinChesapeakeBaysalinity？［J］.Journalof Geophysical Research:Oceans，2008，113（C9）:C09002.

［45］Keeling R F，K?rtzinger A，Gruber N.Ocean deoxygenation in a warming world［J］.Annual Review of Marine Science，2010，2（1）: 199-229.

［46］Das A，Justi? D，Swenson E.Modeling estuarine-shelf exchanges in a deltaic estuary:implications for coastal carbon budgets and hypoxia［J］.Ecological Modelling，2010，221（7）:978-985.

［47］Officer C B，Biggs R B，Taft J L，et al.Chesapeake bay anoxia: origin，development，and significance［J］.Science，1984，223（4631）:22-27.

［48］Vanderborght J P，Wollast R，Loijens M，et al.Application of a transport-reaction model to the estimation of biogas fluxes in the Scheldt estuary［J］.Biogeochemistry，2002，59（1-2）:207-237.

［49］Hofmann A F，Soetaert K，Middelburg J J.Present nitrogen and carbon dynamics in the Scheldt estuary using a novel 1-D model［J］.Biogeosciences，2008，5（4）:981-1006.

［50］Valle-Levinson A，Wilson R E，Swanson R L.Physical mechanisms leading to hypoxia and anoxia in western long island sound［J］.Environment International，1995，21（5）:657-666.

［51］Ishikawa T，Suzuki T，Qian X.Hydraulic study of the onset of hypoxia in the Tone River Estuary［J］.Journal of Environmental Engineering，2004，130（5）:551-561.

［52］Fennel K，Hu J T，Laurent A，et al.Sensitivity of hypoxia predictions for the northern Gulf of Mexico to sediment oxygen consumption and model nesting［J］.Journal of Geophysical Research:Oceans，2013，118（2）:990-1002.

［53］Sohma A，Sekiguchi Y，Yamada H，et al.Erratum to“a new coastalmarineecosystemmodelstudycoupledwith hydrodynamics and tidal flat ecosystem effect”［Marine Pollution Bulletin 43（2001）187-208］［J］.Marine Pollution Bulletin，2002，44（5）:432-439.

［54］羅琳，李適宇，王東曉.珠江河口夏季缺氧現(xiàn)象的模擬［J］.水科學進展，2008，19（5）:729-735.

Review on the marine hypoxia

ZHENG Jing-jing1，2，LIU Gui-mei1，GAO Shan1
（1.National Marine Environmental Forecasting Center，Beijing 100081 China；2.State Key Laboratory of Marine Environmental Science，Xiamen University，Xiamen 361005 China）

Many semi-enclosed bays，estuaries and coasts of the world have different degrees of hypoxia. Influenced by the natural changes and human activities，the phenomenon of marine hypoxia has become increasingly serious，threatening the marine ecological environment.The problem of marine hypoxia has caused wide attention in recently years.In this paper，the domestic and world-wide situation and biochemical mechanism of marine hypoxia is reviewed.Meanwhile，the physical aspects，the effects of the factors such as runoff，upwelling，wind，climate warming and so on，are also analyzed.The development and situation of the numerical model of hypoxia are discussed respectively，including box model，one dimensional model，two dimensional model and three dimensional model.Based on the main study progress at present，this research puts forward the key points and directions of the future research of hypoxia in China.

hypoxia；mechanism；model

P731

A

1003-0239（2016）04-0088-10

10.11737/j.issn.1003-0239.2016.04.011

2015-06-01

國家自然科學基金（41222038）；國家海洋局海洋公益性行業(yè)科研專項項目（201205018）；中國科學院戰(zhàn)略性先導科技專項（XDA11020104）；國家重點基礎研究發(fā)展計劃“973計劃”（2011CB403606）。

鄭靜靜（1990-），女，碩士研究生在讀，主要從事海洋生態(tài)動力學模型研究。E-mail：jingjing.zheng@foxmail.com