羅曉凡，魏 皓

（天津科技大學(xué) 海洋科學(xué)與工程學(xué)院，天津 300457）

碳元素（主要是CO2）在大氣圈、水圈、巖石圈、生物圈之間轉(zhuǎn)移和交換的過程稱為碳循環(huán)。由于人類活動的影響，全球海面上大氣CO2的年均濃度已由1980年的338.80μatm升高至2012年的392.53μatm，并保持持續(xù)攀升的狀態(tài)（http：／／www.esrl.noaa.gov／gmd／ccgg／trends／global.html）。在1958—2008年間，平均約43%的人為CO2停留在大氣中，其余被陸地和海洋吸收［1］。據(jù)最新觀測數(shù)據(jù)估算全球海洋每年吸收人為CO22.0±1.0PgC［2］，一定程度上緩解了大氣CO2濃度的上升速度以及人為驅(qū)動的氣候變化幅度［3］，但這同時也擾動了海洋碳酸鹽體系，導(dǎo)致了海水pH的下降趨勢，由此引發(fā)的海洋環(huán)境問題對海洋生態(tài)系統(tǒng)產(chǎn)生了潛在的影響［4－5］。海洋作為全球碳循環(huán)重要的組成部分，在保持其自身環(huán)境健康的前提下，能夠在吸收大氣CO2、減緩氣候變暖中起多大的作用，未來又將扮演什么樣的角色，這已然成為當(dāng)前研究的熱點問題。

陸架海是陸地與大洋的連接帶，受人類活動影響最為顯著，并且具有復(fù)雜的潮波系統(tǒng)、環(huán)流結(jié)構(gòu)、季節(jié)性河流輸入、較高的生物生產(chǎn)力等特點，各種過程變化劇烈，空間梯度大。陸架海的面積雖僅占全球海洋總面積的8%，其初級生產(chǎn)卻是全球總初級生產(chǎn)的19%～28%［6］，因此陸架海生態(tài)系統(tǒng)的固碳作用對全球碳收支平衡有不可忽略的貢獻(xiàn)，尤其是“碳失匯”的問題被提出后，陸架海對大氣CO2的源／匯作用成為平衡全球CO2收支不確定性因素之一，引起熱議［7－11］。Takahashi［2］給出的全球氣候態(tài)年均凈海－氣 CO2交換的數(shù)據(jù)中，變化迅速且復(fù)雜的近岸水體數(shù)據(jù)相對匱乏。IPCC－AR4所使用的氣候－碳循環(huán)耦合模式比較計劃的多模式數(shù)據(jù)（Coupled Climate－Carbon Cycle Model Intercomparison Project，C4MIP）也沒有包含時空變化劇烈的陸架海過程，給氣候預(yù)測帶來了不確定性［12］。而為了更準(zhǔn)確地理解海洋在大氣CO2濃度不斷上升情景下的作用與響應(yīng)，陸架海及近岸水體碳通量的量化研究得到越來越多的關(guān)注［7－11］。

本研究綜述了近年來對陸架海CO2源／匯格局的爭論，指出通過建立3維－物理－生物－碳化學(xué)耦合數(shù)值模式量化陸架海碳循環(huán)關(guān)鍵過程的科學(xué)發(fā)展趨勢?；陉懠芎Ｌ佳h(huán)模式的研究現(xiàn)狀，從物理場、生物泵、碳化學(xué)系統(tǒng)等方面總結(jié)了建模中應(yīng)考慮的關(guān)鍵控制過程，提出了構(gòu)建東中國陸架海碳循環(huán)模式的基本思路以及擬解決的科學(xué)問題。

1 陸架海對大氣CO2的源／匯作用

陸架海對大氣CO2的源／匯作用是碳循環(huán)研究中最基礎(chǔ)也是最重要的部分，它的強弱反映陸架海對氣候變化的調(diào)節(jié)能力，目前還處于熱烈的爭論之中。部分學(xué)者認(rèn)為全球陸架海是向大氣釋放CO2的源［13－15］，其根據(jù)是河流輸入的有機碳通量大于海洋有機碳埋藏速率，剩余部分很有可能被分解為CO2返回大氣。Tsunogai等［16］在1999年基于東海觀測結(jié)果提出了“陸架泵”的概念，認(rèn)為如果以東中國海吸收大氣CO2的效率去估算，全球陸架海每年可吸收1PgC。Yool等［17］采用GCM模型量化“陸架泵”作用探究其在全球陸架海的潛在強度并指出“陸架泵”效率在不同海區(qū)變化不一，外推估算全球陸架海每年可吸收0.6PgC，但該模型粗糙的分辨率對研究陸架海問題存在很大局限性與不確定性。Thomas等［18］與Bozec等［19］使用觀測數(shù)據(jù)對北海的碳通量及碳輸出的研究證實并支持了“陸架泵”假說。Bates等［20］指出陸架碳泵作用使得楚科奇海陸架成為大氣CO2較強的碳匯海區(qū)。然而，陸架海是多過程共同作用的海區(qū)，其物理、生物、化學(xué)機制不盡相同，因此Cai等［7］根據(jù)陸架海的流場、形態(tài)、緯度等進(jìn)行分類研究并指出中高緯陸架海區(qū)是大氣CO2匯，吸收CO2約0.33PgC yr－1；低緯陸架海區(qū)是大氣CO2源，釋放CO2約0.11PgC yr－1。Chen等［10］則搜集了大量前人關(guān)于陸架海、河口、鹽沼濕地及紅樹林等海－氣CO2通量的論文，給出一個認(rèn)識：陸架海是大氣CO2匯，近岸生態(tài)系統(tǒng)是大氣CO2源。但Xue等［21］指出地處中緯度的開闊陸架?！宵S海中部缺少“陸架泵”需要的物理環(huán)境進(jìn)行有效的碳輸出，其年均尺度上體現(xiàn)為大氣CO2凈源，這與早期Cai等［7］和Chen等［10］的結(jié)論有出入。這些爭論主要源于陸架海表層海水CO2分壓（pCO2）強烈的時間變化與巨大的空間分布不均勻性，同時，一次調(diào)查數(shù)據(jù)不能絕對代表一個月份或是一個季節(jié)的特征，這便是實測數(shù)據(jù)的局限性所在，而“陸架泵”假說是否成立，還需結(jié)合數(shù)值模式進(jìn)行驗證。此外，Chen等［11］總結(jié)了全球106個河口的海－氣CO2交換，指出大多數(shù)河流與河口都是大氣CO2源，但是大河輸送的營養(yǎng)鹽在河流羽流區(qū)被用來進(jìn)行光合作用，一定程度上使水體呈現(xiàn)CO2不飽和狀態(tài)，因此，大河的羽流區(qū)，甚至是距離河口幾百公里遠(yuǎn)的海區(qū)都可能成為CO2匯。Tseng等［22］認(rèn)為長江沖淡水流量影響著CO2的吸收，三峽大壩建成后，長江沖淡水流量減小，可能導(dǎo)致CO2吸收能力的下降，甚至?xí)箹|海陸架由大氣CO2匯轉(zhuǎn)為源，這個結(jié)論有待驗證。Dai等［23］根據(jù)對南海的研究提出了大洋主控型邊緣海碳循環(huán)概念，認(rèn)為來自大洋深部的外源溶解無機碳通過上升流和躍密混合被輸送至真光層后，如果不被同樣來源的營養(yǎng)鹽所支持的有機碳生產(chǎn)完全消耗，其過剩部分就會溢出至大氣中，而在河流主控型邊緣海，陸源有機物與營養(yǎng)鹽輸入含量的競爭調(diào)控該海區(qū)生物的新陳代謝作用，從而影響CO2通量?？梢?，不同陸架海區(qū)水文環(huán)境、動力過程存在差別，不同生態(tài)系統(tǒng)生物泵運作不同，陸架海源／匯格局的分布各具特色。碳酸鹽系統(tǒng)實測資料的匱乏及其時空局限性，限制了人們對海洋碳循環(huán)整體的進(jìn)一步認(rèn)識，這使得海洋碳循環(huán)模式快速發(fā)展，同時也促進(jìn)了人們對海洋碳循環(huán)關(guān)鍵過程進(jìn)行量化的研究。

2 陸架海碳循環(huán)模式的發(fā)展

相較于大洋，陸架海系統(tǒng)的復(fù)雜性使得陸架海碳循環(huán)模式的研究更為困難，目前，國際上對陸架海碳循環(huán)模式的報道不多，我國對東中國陸架海3維－物理－生物－碳循環(huán)模式的研究尚在起步階段，陸架海碳循環(huán)模式的建立已成為各國氣候預(yù)測模型急需攻克的難題。

自20世紀(jì)60年代起，大洋碳循環(huán)模式經(jīng)歷了由簡單到復(fù)雜的發(fā)展。從早期的箱式收支模式到三維海洋無機碳循環(huán)模式，再到基于海洋動力學(xué)耦合生物地球化學(xué)過程的三維碳循環(huán)模式，直至今天的氣候－碳循環(huán)耦合模式。然而直至20世紀(jì)90年代，陸架海碳循環(huán)模式才開始發(fā)展。CO2在海洋中的轉(zhuǎn)移和歸宿，即海洋吸收、轉(zhuǎn)移大氣CO2的能力及CO2在海洋中的循環(huán)機制均集中體現(xiàn)為其物理、化學(xué)和生物過程［24］。這也是陸架海碳循環(huán)模式的三大模塊，三者的耦合完成了溶解無機碳（DIC）、溶解有機碳（DOC）、顆粒有機碳（POC）等在陸架海的遷移轉(zhuǎn)化。

本研究在Hofmann等［25］對陸架海碳循環(huán)模式綜述的基礎(chǔ)上，列出近年來有關(guān)模擬陸架海碳收支的文章（表1），雖然不同海區(qū)模式所含變量以及邊界條件不同，但模式基本包括了陸源輸入的淡水、營養(yǎng)鹽、碳，海－氣交換的CO2及大氣沉降的營養(yǎng)鹽，沉積物－水界面交換通量，生物作用下營養(yǎng)鹽及碳的遷移轉(zhuǎn)化，跨陸架輸運及陸架海與外海的交換。同時，基于地球物理流體力學(xué)原理所建立的海洋環(huán)流模式提供了更為真實，結(jié)構(gòu)更復(fù)雜的海洋物理場，使得對陸架海碳循環(huán)的模擬也更為準(zhǔn)確。

表1 部分研究陸架邊緣海海－氣交換通量及碳輸出的模式Table 1 Sea－Air CO2flux（SACF）and carbon transportation models developed for the continental shelf seas

3 陸架海碳循環(huán)模式中的關(guān)鍵過程

3.1 關(guān)鍵物理過程

現(xiàn)有的數(shù)值模式對海洋碳循環(huán)的模擬都存在缺點和不足之處，Doney等［38］對海洋碳循環(huán)模式比較計劃——階段2（the phase 2of the Ocean Carbon－cycle Model Intercomparison Project，OCMIP－2）的13個全球碳循環(huán)模式的水動力場以及子午圈翻轉(zhuǎn)流進(jìn)行比較指出，這些缺點與不足主要與物理環(huán)流有關(guān)，碳循環(huán)模式結(jié)果在生物地球化學(xué)過程的差異有40±25%都是來自海洋環(huán)流模式［38］。生源要素必然隨海水運動而遷移，除了由海洋各個界面輸入、輸出外，海洋環(huán)流是決定各種物質(zhì)時空分布的主要因素［39］。因此，陸架海碳循環(huán)模式對海洋動力模式有較高的要求，同時也促進(jìn)了海洋模式向高分辨率、多過程耦合發(fā)展。陸架海多上升、下降流系統(tǒng)，其垂直方向輸送的營養(yǎng)鹽以及碳通量也影響著表層pCO2的分布［27］。在南海，上升流是導(dǎo)致越南沿岸與呂宋島西北海區(qū)pCO2空間分布差異的因素之一［30］，南海北部陸架海區(qū)碳酸鹽系統(tǒng)的動態(tài)格局也主要由河流羽流和近岸上升流影響下的物理和生物強迫控制［40］。此外，在有利于下降流形成的風(fēng)驅(qū)動下，沿岸流和底層離岸輸運有機碳也是陸架海碳循環(huán)模型中需要考慮的物理過程。

海洋的暖化降低了海水CO2溶解度，可能使海洋碳存儲減小，從而對大氣CO2濃度起正反饋作用［41］，加劇全球氣候變化。因此除環(huán)流輸運外，對海水溫度長期變化的模擬、抓住區(qū)域海溫度對氣候變化的響應(yīng)特征是非常重要的。Chai等［30］模式結(jié)果表明水溫對海水pCO2季節(jié)變化的影響比生物活動更強，是控制南海pCO2的主導(dǎo)因素。已有研究表明東中國陸架海增溫更加顯著［42－43］，由此而帶來的對碳循環(huán)的影響還需科學(xué)評估。

在歐洲北海，顯著不同的南北兩個生物地球化學(xué)省的海－氣CO2交換是溫度與生物作用影響平衡的結(jié)果，而這很大程度上取決于水體的層化以及在層化影響下生物作用的差異［31］。夏末，北海南部水體混合均勻，呈非自養(yǎng)型，向大氣凈釋放CO2，而在出現(xiàn)季節(jié)性分層的北海北部，表層水體呈自養(yǎng)型，該區(qū)是大氣CO2的匯［19］。可見，陸架海碳循環(huán)模式對海洋模式中混合的參數(shù)化方案提出了更精細(xì)的要求，以能夠反演時空快速變化的層化／混合、鋒面、河流羽流結(jié)構(gòu)等。

3.2 生物泵

生物泵作用是陸架海碳循環(huán)中重要的環(huán)節(jié)，碳從真光層至深層的傳遞要依靠生物的作用。在陸架海，大約一半的有機碳通過生物泵作用轉(zhuǎn)移到深海［25］。生物過程驅(qū)動了海洋對CO2的吸收，Kühn等［32］和Lorkowski等［36］對北海碳收支的模擬中指出，如果模式不包含生物過程，北海北部將會和北海南部一樣表現(xiàn)為大氣CO2的源，這與實際調(diào)查結(jié)果不符。在羅斯海，海冰消融的海區(qū)因較高的生物初級生產(chǎn)而成為大氣CO2的強匯區(qū)［29］。因此陸架海碳循環(huán)模式耦合低營養(yǎng)級生態(tài)模式是當(dāng)前碳循環(huán)模式發(fā)展主流，包含詳細(xì)的生物過程與營養(yǎng)鹽循環(huán)的碳循環(huán)模式可以更準(zhǔn)確的對碳化學(xué)系統(tǒng)收支進(jìn)行估算。Borges等［44］指出在具有高初級生產(chǎn)力的近岸環(huán)境，碳化學(xué)系統(tǒng)對富營養(yǎng)化的響應(yīng)要強于對酸化的響應(yīng)。此外，底棲生物過程影響著底層營養(yǎng)鹽的循環(huán)以及堿度的變化，模式中對碳化學(xué)系統(tǒng)估算的不確定性有一部分是來自底棲作用［4］。在經(jīng)典生物泵中，微生物通常被認(rèn)為是分解者，而研究表明，細(xì)菌還可以把活性DOC高效地轉(zhuǎn)化為惰性DOC，形成巨大的惰性DOC碳庫，進(jìn)入的碳循環(huán)的慢性循環(huán)（千年尺度）中，起到儲碳作用，這一過程就是“微生物碳泵”［45］。目前，雖然對微生物碳泵有了理論的認(rèn)識，但將這一過程進(jìn)行模式化，還是面臨很大的挑戰(zhàn)。復(fù)雜的生物過程不僅增加了模式的難度，同時也增加了誤差在模式中的傳遞。此外，碳循環(huán)與營養(yǎng)鹽循環(huán)都有其獨特的控制過程，模式中應(yīng)考慮將兩者解耦，這樣也可以使浮游植物的生長不受碳而只受營養(yǎng)鹽的限制［4，36］。

3.3 碳化學(xué)系統(tǒng)

3.3.1 界面碳通量

在現(xiàn)有日趨完善的生態(tài)動力學(xué)模型的發(fā)展基礎(chǔ)之上耦合碳化學(xué)系統(tǒng)是當(dāng)前陸架海碳循環(huán)模式的發(fā)展方向，因此，對碳化學(xué)系統(tǒng)的整體認(rèn)知是模式研究者建模中重要的環(huán)節(jié)。陸架海碳化學(xué)系統(tǒng)參與了海－氣界面、陸－海界面、沉積物－水界面以及水體內(nèi)部的化學(xué)反應(yīng)，調(diào)控著水體pH，并影響其他化學(xué)反應(yīng)，對形成和維持海洋生命活動有重要作用。在海洋中，碳的存在形式主要為DIC、DOC、POC，且約98%的碳是以溶解無機碳的形式存在［12］。

評估陸架海在全球氣候變化中的作用，海－氣CO2交換是最基礎(chǔ)的一環(huán)，也是對碳酸鹽系統(tǒng)的平衡起重要作用的界面過程。海水碳酸鹽系統(tǒng)是由DIC、總堿度（Talk）、pH值、pCO2等通過基本的熱力學(xué)關(guān)系相聯(lián)系，包括了碳酸、硼酸、磷酸、硅酸、水等弱酸的電離平衡［46］。在模式中，為了減小其復(fù)雜度，通常只考慮碳酸－硼酸－水體系。目前，模式估算海－氣CO2通量（SACF）多使用公式（1）：

其中，fice為海冰覆蓋率；為本文章上述所指的pCO2，即海水表層CO2分壓；為海面上大氣CO2分壓；α為CO2的溶解度，模式中通常使用Weiss［47］給出的與溫鹽相關(guān)的經(jīng)驗公式；k為液相傳輸速率，模式中通常使用Wanninkhof［48］給出的與海面10m風(fēng)速以及溫度相關(guān)的經(jīng)驗公式。盡管使用上述方法存在不確定性，主要體現(xiàn)在氣體交換系數(shù)（k﹒α），但使用該方法估算SACF是被普遍接受的［46］。海－氣界面是碳循環(huán)與人類活動影響直接相關(guān)的關(guān)鍵界面，Omstedt等［49］對波羅的海pCO2長期變化的模擬表明，工業(yè)革命前波羅的海表現(xiàn)為凈源作用，而工業(yè)革命后，pCO2季節(jié)變化增強，年際變異上源／匯作用都有出現(xiàn)；同時對pH的長期模擬指出，在工業(yè)化之前pH維持一個穩(wěn)定狀態(tài)，而工業(yè)化后，隨著大氣CO2濃度的增加，pH有小幅下降。在地中海，1 980s其表現(xiàn)為大氣CO2的凈源，每年釋放1.5×1012mol C，而到了2 000s，大氣CO2濃度較1 980s增加了30μatm，使得大氣向海洋的CO2通量增加，地中海表現(xiàn)出弱匯作用，平均每年吸收0.2×1012mol C［50］。Lu等［51］對南海SEATS站位建立1維模型的模擬結(jié)果表明，海－氣CO2交換調(diào)控pCO2的季節(jié)變化，甚至超過SST的影響，且因其為寡營養(yǎng)鹽海區(qū)，生物過程的貢獻(xiàn)較??；在日變化上，pCO2的調(diào)控因子與局地物理和生物環(huán)境有關(guān)。

與大洋不同，陸架海直接受人類擾動的另一個界面就是陸－海界面，主要表現(xiàn)為河流輸入帶來的影響。河流碳通量在特定時間和地區(qū)是單向輸送的，主要包括DIC、DOC、POC、PIC（顆粒無機碳）的輸送［52］。2004—2008年，長江平均每年向東海輸送DIC為1.56×1013g［53］。Wang等［54］研究表明2009年，黃海和長江向海輸送的DOC分別為3.20×1010g和1.58×1012g，POC分別為3.89×1011g和1.52×1012g，可能對河口及其近岸的生物化學(xué)循環(huán)和生態(tài)系統(tǒng)產(chǎn)生重要影響。Raymond等［55］指出當(dāng)考慮近岸區(qū)域DIC和Talk的季節(jié)、年際、長期變化趨勢時，河流輸入的準(zhǔn)確性是重要的。但在模式中，河流碳輸入資料的匱乏，時間變化的不確定性，模式分辨率不足等都可能使得河口低鹽區(qū)的模擬不甚理想［4，25］，另外河口過程、濕地等增加了問題的復(fù)雜性，陸－海界面過程是陸架海碳循環(huán)模式研究的一個難點。由于河流帶來流域植被代謝的碳、營養(yǎng)物質(zhì)以及農(nóng)耕區(qū)的化肥等，因此，將陸架海碳循環(huán)模型與陸地碳循環(huán)模型、水文模型銜接，獲得連續(xù)的河流邊界條件是河流主控型陸架海碳循環(huán)模型的一個發(fā)展方向。此外，海底地下水（Submarine Groundwater Discharge－SGD）也是陸地向海洋輸送水量和物質(zhì)的重要通道，被認(rèn)為是近海營養(yǎng)鹽、碳等溶解物質(zhì)的源［56］。全球尺度上，SGD碳通量輸送粗略估算約為0.23PgC yr－1，相當(dāng)于全球河流碳輸送的30%［57］。在南海北部，SGD雖占珠江徑流量的12%～21%，但其攜帶的DIC可能高達(dá)（153～347）×109mol yr－1，相當(dāng)于珠江DIC輸入量的23%～53%［58］?？梢姡琒GD對近岸碳收支平衡的潛在影響。然而現(xiàn)今對SGD碳通量的研究還不全面，在陸架海碳循環(huán)模式中對SGD的定量化是比較困難的，SGD數(shù)值模式的發(fā)展雖然為研究SGD開拓了視野，但目前可以用來對模式進(jìn)行有效驗證的實測數(shù)據(jù)及經(jīng)驗數(shù)據(jù)非常少，由此而帶來的不確定性限制了模式的廣泛應(yīng)用與發(fā)展［59］。

沉積物－水界面交換體現(xiàn)在有機碳沉積與再懸浮以及有機碳在微生物作用下的再礦化作用，在動力較強、海底有機污染嚴(yán)重的陸架海，這個界面尤其重要。大洋沉積通過對有機碳和CaCO3的埋藏可以作為永久碳庫，然而在東中國陸架海，海底是非河流影響區(qū)域重要的營養(yǎng)鹽來源，強烈的潮流、風(fēng)場的擾動使得顆粒物再懸浮嚴(yán)重，加之活躍的底棲生物作用，有機碳的保存效率大幅降低［60－62］。由于該界面精細(xì)過程觀測難以實施，僅有零星底界面交換實驗數(shù)據(jù)，而交換又與動力過程密切相關(guān)，目前還沒有一個完善的底界面交換模型，這個過程的參數(shù)化是陸架碳循環(huán)模型中另一難點，需要突破。

陸架海－外海界面碳通量和營養(yǎng)鹽通量受開邊界條件的影響。Kühn等［32］指出使用來自二氧化碳信息處理中心（Carbon Dioxide Information Analysis Center，CDIAC）的氣候態(tài)DIC、Talk數(shù)據(jù)作為北海的開邊界條件會高估北海DIC濃度。Artioli等［4］指出完善初始條件與邊界條件，尤其是河流輸入以及與波羅的海通量交換，可以顯著減小對北海營養(yǎng)鹽模擬的偏差。在中國東海，高溫高鹽的黑潮水入侵為陸架海輸送了可觀的營養(yǎng)鹽，其中硝酸鹽輸送入東海的平均值為170.8±41.6kmol s－1，約為磷酸鹽輸送量的13.6倍［63］，這對于陸架海營養(yǎng)鹽循環(huán)、生物生產(chǎn)以及整個碳循環(huán)系統(tǒng)都有重要影響。

3.3.2 碳的遷移轉(zhuǎn)化

在碳的遷移轉(zhuǎn)化過程中（圖1），生物作用主導(dǎo)有機碳與無機碳的相互轉(zhuǎn)化。而DOM的光降解主要表現(xiàn)為有色溶解有機物（CDOM）的光降解，陸源CDOM入海后通過光降解作用在近海海域發(fā)生轉(zhuǎn)化與消除，潛在地影響海洋生態(tài)系統(tǒng)及海洋碳循環(huán)［64－66］。賀志鵬等［67］指出黃海表層至20m層DOC濃度呈上升趨勢可能就是表層DOC發(fā)生光降解作用的直接證據(jù)。Yang等［65］估算得黃東海每年約5.4～10.5TgC的DOC發(fā)生光礦化，約占初級生產(chǎn)的2.9%～5.7%。DIC以碳酸鹽形式移出水體主要表現(xiàn)為CaCO3的產(chǎn)出，包括珊瑚、顆石藻外殼、某些浮游動物的碳酸鹽外殼和骨架等。全球海洋碳酸鹽的生產(chǎn)和輸出以顆粒無機碳的形式創(chuàng)造了較高的碳輸出效率（0.4～1.8GtC yr－1）［68］。Cao等［69］指出CaCO3的形成與溶解調(diào)節(jié)著海洋 DIC系統(tǒng)并最終調(diào)控海洋對大氣CO2的緩解能力，在南海，近岸淺水海區(qū)CaCO3的溶解對外海上層水體中累積的Ca2＋和總堿度有貢獻(xiàn)。近岸水體酸化很可能會減少碳酸鹽的生產(chǎn)，加速碳酸鹽基底的溶解，從而打破碳化學(xué)系統(tǒng)的平衡，不利于鈣化生物的生存。因此，對于碳酸鹽沉積底質(zhì)的陸架海以及珊瑚生態(tài)系統(tǒng)，在模式中應(yīng)酌情考慮碳酸鹽的生成與溶解過程。對于POC的源匯項，除受陸源與生物作用影響外，還受到沉降及再懸浮過程的影響。海洋對大氣CO2的調(diào)節(jié)作用重點在于碳的垂直遷移并最終被埋藏于海底的過程。郭志剛等［70］等指出有機碳從水體向海底沉積物轉(zhuǎn)移并保存的條件要滿足4點：1）有充足的POC向海底轉(zhuǎn)移；2）有一定的沉積速率；3）適于細(xì)顆粒沉積的弱水環(huán)境；4）相對還原的海底環(huán)境。然而陸架海動力調(diào)整變化快，沉降并保存在沉積物中的POC在受到動力擾動時，便會發(fā)生再懸浮，使POC沉降受限。據(jù)統(tǒng)計，大陸架及大陸坡的沉積物中，有40%～85%的有機碳要發(fā)生再懸?。?1］。長江口及鄰近海域海－氣交換的CO2僅有13%以顆粒有機碳的形式垂直轉(zhuǎn)移而形成表層沉積物的有機碳，其余87%沒有被埋藏的碳通過生物泵或動力過程進(jìn)入水中再循環(huán)［72］。

圖1 陸架海碳的遷移轉(zhuǎn)化示意圖Fig.1 A sketch map showing the migration and transformation of carbon in the continental shelf seas

4 東中國陸架海碳循環(huán)模式的構(gòu)建思想

東中國海地處東亞季風(fēng)區(qū)，渤海、黃海的全部以及東海的2／3都屬于寬闊的陸架海，受到長江、黃河等陸地徑流的強烈影響以及臺灣暖流、黑潮入侵等帶來的強烈作用。目前，基于水動力模型的營養(yǎng)鹽－低營養(yǎng)級生態(tài)動力學(xué)模型在東中國海已取得較好的應(yīng)用［73］，在此基礎(chǔ)上耦合碳化學(xué)系統(tǒng)是碳循環(huán)模型發(fā)展的基本方向。東中國陸架海碳循環(huán)模型的基本構(gòu)想是在生態(tài)動力學(xué)模型中增加DIC作為主要的狀態(tài)變量，同時增加Talk、pCO2、SACF作為診斷變量，其中為減小模式復(fù)雜度，可以考慮將Talk與鹽度建立關(guān)系，pCO2通過碳酸鹽體系平衡計算［74］，SACF的計算前文已給出（詳見3.3.1節(jié)）。相較于大洋模型，陸架海碳循環(huán)模型顯著的特點在于陸－海界面和沉積物－水界面碳通量的參數(shù)化。目前對河流碳通量動態(tài)變化規(guī)律的研究并不細(xì)致，但河流徑流量主導(dǎo)其無機碳和有機碳入海通量變化［75－76］，由此可根據(jù)徑流量以及DIC、DOC、POC在不同年代入海通量的年際特征值推算碳通量長時間序列作為河流邊界條件。沉積物－水界面所表現(xiàn)的固碳及去碳作用，在年或年際尺度上處于動態(tài)平衡狀態(tài)，季節(jié)尺度上有變化，這在現(xiàn)代陸架環(huán)境下主要與生物過程相關(guān)［70］。這一界面過程目前還沒有精細(xì)的數(shù)值模型可供采用，結(jié)合水動力過程的沉積物－水界面物質(zhì)交換、再懸浮過程的模型量化，應(yīng)該是東中國陸架海碳循環(huán)模型的一個主要特色。

東中國陸架海3維－物理－生物－碳循環(huán)模型的建立，可進(jìn)一步認(rèn)識東中國海表層pCO2、海－氣碳通量、向陸架外碳輸出效率的季節(jié)與年際變化；定量評估東中國海源／匯動態(tài)格局以及生物泵、“陸架泵”的貢獻(xiàn)；河流徑流量變化所致的物理、生物效應(yīng)對河口及其近岸碳通量的影響；沉積－水界面交換通量動力學(xué)參數(shù)化在碳循環(huán)中的作用；并預(yù)測未來陸架海碳循環(huán)系統(tǒng)的演變規(guī)律，為全球碳循環(huán)模型提供可嵌入的東中國陸架海碳循環(huán)模塊，這對于改善全球碳循環(huán)模型，減少其不確定性，準(zhǔn)確評估東亞海區(qū)在氣候變化中的作用都是非常必要的。

（References）：

［1］ LEQUéRéC，RAUPACH M R，CANADELL J G，et al.Trends in the sources and sinks of carbon dioxide［J］.Nature Geoscience，2009，2（12）：831－836.

［2］ TAKAHASHI T，SUTHERLAND S C，WANNINKHOF R，et al.Climatological mean and decadal change in surface ocean pCO2，and net sea－air CO2flux over the global oceans［J］.Deep－Sea Research II，2009，56＊8－10）：554－577.

［3］ FRIEDLINGSTEIN P，COX P，BETTS R，et al.Climate－carbon cycle feedback analysis：results from the C4MIP model intercomparison［J］.Journal of Climate，2006，19（14）：3337－3353.

［4］ ARTIOLI Y，BLACKFORD J C，BUTENSCH?N M，et al.The carbonate system in the North Sea：Sensitivity and model validation［J］.Journal of Marine Systems，2012，102－104：1－13.

［5］ MCGRATH T，KIVIMāE C，TANHUA T，et al.Inorganic carbon and pH levels in the Rockall Trough 1991－2010［J］.Deep－Sea Research I，2012，68：79－91.

［6］ LONGHURST A.Seasonal cycles of pelagic production and consumption［J］.Progress in Oceanography，1995，36：77－167.

［7］ CAI W J，DAI M H，WANG Y C.Air－sea exchange of carbon dioxide in ocean margins：A province based synthesis［J］.Geophysical Research Letters，2006，33（12）：L12603.

［8］ ZHAI W D，DAI M H.On the seasonal variation of air－sea CO2fluxes in the outer Changjiang（Yangtze River）Estuary，East China Sea［J］.Marine Chemistry，2009，117：2－10.

［9］ GUO X H，DAI M H，ZHAI W D，et al.CO2flux and seasonal variability in a large subtropical estuarine system，the Pearl River Estuary，China［J］.Journal of Geophysical Research，2009，114（G3）：G03013.

［10］ CHEN C T A，BORGES A V.Reconciling opposing views on carbon cycling in the coastal ocean：Continental shelves as sinks and nearshore ecosystems as sources of atmospheric CO2［J］.Deep－Sea Research II，2009，56（8－10）：578－590.

［11］ CHEN C T A，HUANG T H，F(xiàn)U Y H，et al.Strong sources of CO2in upper estuaries become sinks of CO2in large river plumes［J］.Current Opinion in Environmental Sustainability，2012，4：179－185.

［12］ IPCC.Summary for Policymakers of the Synthesis Report of the IPCC Fourth Assessment Report［M］.Cambridge，UK：Cambridge University Press，2007.

［13］ VER M B，MACKENZIE F T，LERMAN A.Biogeochemical responses of the carbon cycle to natural and human perturbations：Past，present，and future［J］.American Journal of Science，1999，299：762－801.

［14］ MACKENZIE F T，VER L M，LERMAN A.Coastal－zone biogeochemical dynamics under global warming［J］.International Geology Review，2000，42（3）：193－206.

［15］ FASHAM M J R，BALINO B M ，BOWLES C，et al.A new vision of ocean biogeochemistry after a decade of the Joint Global Ocean Flux Study（JGOFS）［J］.Ambio（Special Issue），2001，10：4－30.

［16］ TSUNOGAI S，WATANABE S，SATO T.Is there a"continental shelf pump"for the absorption of atmospheric CO2？［J］.Tellus B，1999，51（3）：701－712.

［17］ YOOL A，F(xiàn)ASHAM M J R.An examination of the"continental shelf pump"in an open ocean general circulation model［J］.Global Biogeochemical cycles，2001，15（4）：831－844.

［18］ THOMAS H，BOZEC Y，ELKALAY K，et al.Enhanced open ocean storage of CO2from shelf sea pumping［J］.Science，2004，304（5673）：1005－1008.

［19］ BOZEC Y，THOMAS H，ELKALAY K，et al.The continental shelf pump for CO2in the North Sea－evidence from summer observation［J］.Marine Chemistry，2005，93：131－147.

［20］ BATES N R.Air－sea CO2fluxes and the continental shelf pump of carbon in the Chukchi Sea adjacent to the Arctic Ocean［J］.Journal of Geophysical Research，2006，111（C10）：C10013.

［21］ XUE L，ZHANG L J，CAI W J，et al.Air－sea CO2fluxes in the southern Yellow Sea：An examination of the continental shelf pump hypothesis［J］.Continental Shelf Research，2011，31（18）：1904－1914.

［22］ TSENG C M，LIU K K，GONG G C，et al.CO2uptake in the East China Sea relying on Changjiang runoff is prone to change［J］.Geophysical Research Letters，2011，38（24）：L24609.

［23］ DAI M H，CAO Z M，GUO X H，et al.Why are some marginal seas sources of atmospheric CO2？［J］.Geophysical Research Letters，2013，40（10）：2154－2158.

［24］ XU Y F，PU Y F，ZHAO L.Advances in the studies of ocean carbon－cycle model［J］.Advances in Earth Science，2005，20（10）：1106－1115.徐永福，浦一芬，趙亮.海洋碳循環(huán)模式的研究進(jìn)展［J］.地球科學(xué)進(jìn)展，2005，20（10）：1106－1115.

［25］ HOFMANN E E，CAHILL B，F(xiàn)ENNEL K，et al.Modeling the Dynamics of Continental Shelf Carbon［J］.Annual Review of Marine Science，2011，3：93－122.

［26］ WALSH J J，DIETERLE D A.CO2cycling in the coastal ocean.I－A numerical analysis of the southeastern Bering Sea with applications to the Chukchi Sea and the northern Gulf of Mexico［J］.Progress in Oceanography，1994，34（4）：335－392.

［27］ IANSON D，ALLEN S E.A two－dimensional nitrogen and carbon flux model in a coastal upwelling region［J］.Global Biogeochemical Cycles，2002，16（1）：1－16.

［28］ FENNEL K，WILKIN J，PREVIDI M，et al.Denitrification effects on air－sea CO2flux in the coastal ocean：Simulations for the northwest North Atlantic［J］.Geophysical Research Letters，2008，35（24）：L24608.

［29］ ARRIGO K R，DIJKEN G V，LONG M.Coastal Southern Ocean：A strong anthropogenic CO2sink［J］.Geophysical Research Letters，2008，35（1）：L21602.

［30］ CHAI F，LIU G M，XUE H J.Seasonal and Interannual Variability of Carbon Cycle in South China Sea：A Three－Dimensional Physical－Biogeochemical Modeling Study［J］.Journal of Oceanography，2009，65：703－720.

［31］ PROWE A E F，THOMAS H，PāTSCH J，et al.Mechanisms controlling the air－sea CO2flux in the North Sea［J］.Continental Shelf Research，2009，29：1801－1808.

［32］ KüHN W，PāTSCH J，THOMAS H，et al.Nitrogen and carbon cycling in the North Sea and exchange with the North Atlantic－A model study，Part II：Carbon budget and fluxes［J］.Continental Shelf Research，2010，30：1701－1716.

［33］ DRUON J N，MANNINO A，SIGNORINI S，et al.Modeling the dynamics and export of dissolved organic matter in the Northeastern U.S.continental shelf［J］.Coastal and Shelf Science，2010，88（4）：488－507.

［34］ KUZNETSOV I，NEUMANN T.Simulation of carbon dynamics in the Baltic Sea with a 3Dmodel［J］.Journal of Marine Systems，2013，111－112：167－174.

［35］ W?HLSTR?M I，OMSTEDT A，BJ?RK G，et al.Modelling the CO2dynamics in the Laptev Sea，Arctic Ocean：Part I［J］.Journal of Marine Systems，2012，102－104：29－38.

［36］ LORKOWSKI I，PāTSCH J，MOLL A，et al.Interannual variability of carbon fluxes in the North Sea from 1970to 2006－Competing effects of abiotic and biotic drivers on the gas－exchange of CO2［J］.Estuarine，Coastal and Shelf Science，2012，100：38－57.

［37］ ?RTHUN M，BELLERBY R G J，OMAR A M，et al.Spatiotemporal variability of air－sea CO2fluxes in the Barents Sea，as determined from empirical relationships and modeled hydrography［J］.Journal of Marine Systems，2012，98－99：40－50.

［38］ DONEY S C，LINDSAY K，CALDEIRA K，et al.Evaluating global ocean carbon models：The importance of realistic physics［J］.Global Biogeochemical Cycles，2004，18（3）：GB3017.

［39］ SUN W X，JIANG W S，LI L.Numerical Models of Hydrodynamics in the Coastal Environment［M］.Beijing：Science Press，2004.孫文心，江文勝，李磊.近海環(huán)境流體動力學(xué)數(shù)值模型［M］.北京：科學(xué)出版社，2004.

［40］ CAO Z M，DAI M H，ZHENG N，et al.Dynamics of the carbonate system in a large continental shelf system under the influence of both a river plume and coastal upwelling［J］.Journal of Geophysical Research，2011，116（G2）：G02010.

［41］ DONEY S C，TILBROOK B，ROY S，et al.Surface－ocean CO2variability and vulnerability［J］.Deep－Sea Research II，2009，56（8－10）：504－511.

［42］ WEI H，SHI J，LU Y Y，et al.Interannual and long－term hydrographic changes in the Yellow Sea during 1977—1998［J］.Deep－Sea Research II，2010，57（11－12）：1025－1034.

［43］ LI J X，WEI H，LUO X F.Spatial and temporal inter－annual variations of sea surface temperature in the continental shelf of the East China Sea［J］.Oceanologia et Limnologia Sinica，2012，43（3）：624－634.李家星，魏皓，羅曉凡.東海陸架表層水溫年際變化時空特征分析［J］.海洋與湖沼，2012，43（3）：624－634.

［44］ BORGES A V，GYPENS N.Carbonate chemistry in the coastal zone responds more strongly to eutrophication than to ocean acidification［J］.Limnology ＆ Oceanography，2010，55（1）：346－353.

［45］ JIAO N Z.Carbon fixation and sequestration in the ocean，with special reference to the microbial carbon pump［J］.Scientia Sinica Terrae，2012，42（10）：1473－1486.焦念志.海洋固碳與儲碳－并論微生物在其中的重要作用［J］.中國科學(xué)：地球科學(xué)，2012，42（10）：1473－1486.

［46］ XU Y F，ZHAO L，PU Y F，et al.Uncertainties in the estimate of the air－sea exchange flux of carbon dioxide［J］.Earth Science Frontiers，2004，11（2）：565－571.徐永福，趙亮，浦一芬，等.二氧化碳海氣交換通量估計的不確定性［J］.地學(xué)前緣，2004，11（2）：565－571.

［47］ WEISS R F.Carbon dioxide in water and seawater：the solubility of a non－ideal gas［J］.Marine Chemistry，1974，2：203－215.

［48］ WANNINKHOF R.Relationship between wind speed and gas exchange over the ocean［J］.Journal of Geophysical Research，1992，97（C5）：7373－7382.

［49］ OMSTEDT A，GUSTAFSSON E，WESSLANDER K.Modelling the uptake and release of carbon dioxide in the Baltic Sea surface water［J］.Continental Shelf Research，2009，29（7）：870－885.

［50］ TAILLANDIER V，ORTENZIO F D，ANTOINE D.Carbon fluxes in the mixed layer of the Mediterranean Sea in the1980sand the2000s［J］.Deep－Sea Research I，2012，65：73－84.

［51］ LU Z M，GAN J P，DAI M H.Modeling seasonal and diurnal pCO2variations in the northern South China Sea［J］.Journal of Marine Systems，2012，92（1）：30－41.

［52］ GAO Q Z，SHEN C D.Riverine carbon flux and continental erosion［J］.Advances in Earth Science，1998，13（4）：369－375.高全洲，沈承德.河流碳通量與陸地侵蝕研究［J］.地球科學(xué)進(jìn)展，1998，13（4）：369－375.

［53］ LI D.The Study on the Hydro－Chemical Characteristics and the Flux to the sea about the Rivers in the East of China［D］.Shanghai：East China Normal University，2009.李丹.中國東部若干入海河流水化學(xué)特征與入海通量研究［D］.上海：華東師范大學(xué)，2009.

［54］ WANG X C，MA H Q，LI R H，et al.Seasonal fluxes and source variation of organic carbon transported by two major Chinese Rivers：The Yellow River and Changjiang（Yangtze）River［J］.Global Biogeochemical Cycles，2012，26：GB2025.

［55］ RAYMOND P A，OH N H，TUMER R E，et al.Anthropogenically enhanced fluxes of water and carbon from the Mississippi River［J］.Nature，2008，451：449－452.

［56］ MOORE W S.The role of submarine groundwater discharge in coastal biogeochemistry［J］.Journal of Geochemical Exploration，2006，88（1－3）：389－393.

［57］ ZHANG J，MANDAL A K.Linkages between submarine groundwater systems and the environment［J］.Current Opinion in Environmental Sustainability，2012，4（2）：219－226.

［58］ LIU Q，DAI M H，CHEN W，et al.How significant is submarine groundwater discharge and its associated dissolved inorganic carbon in a river－dominated shelf system？［J］.Biogeosciences，2012，9：1777－1795.

［59］ SANTOS I R，EYRE B D，HUETTEL M.The driving forces of porewater and groundwater flow in permeable coastal sediments：A review［J］.Estuarine，Coastal and Shelf Science，2012，98：1－15.

［60］ KIM G，RYU J W，YANG H S，et al.Submarine groundwater discharge（SGD）into the Yellow Sea revealed by 228Ra and 226Ra isotopes：implications for global silicate fluxes［J］.Earth and Planetary Science Letters，2005，237（1－2）：156－166.

［61］ YUAN Y，WEI H，ZHAO L，et al.Observations of sediment resuspension and settling off the mouth of JiaozhouBay，Yellow Sea［J］.Continental Shelf Research，2008，28（19）：2630－2643.

［62］ HUNG J J，LIN P L，LIU K K.Dissolved and particulate organic carbon in the southern East China Sea［J］.Continental shelf Research，2000，20（4－5）：545－569.

［63］ GUO X Y，ZHU X H，WU Q S，et al.The Kuroshio nutrient stream and its temporal variation in the East China Sea［J］.Journal of Geophysical Research，2012，117（C1）：C01026.

［64］ GUO W D，CHENG Y Y，YU X X，et al.An overview of the Photochemistry of marine chromophoric dissolved organic matter［J］.Marine Science Bulletin，2008，27（3）：107－114.郭衛(wèi)東，程遠(yuǎn)月，余翔翔，等.海洋有色溶解有機物的光化學(xué)研究進(jìn)展［J］.海洋通報，2008，27（3）：107－114.

［65］ YANG G P，REN C Y，LU X L，et al.Distribution，flux，and photoproduction of carbon monoxide in the East China Sea and Yellow Sea in spring［J］.Journal of geophysical research，2011，116（C2）：C02001.

［66］ XIE H X，AUBRY C，BéLANGER S，et al.The dynamics of absorption coefficients of CDOM and particles in the St.Lawrence estuarine system：Biogeochemical and physical implications［J］.Marine Chemistry，2012，128－129：44－56.

［67］ HE Z P，SONG J M，ZHANG N X.Analysis of biogeochemical characteristics of dissolved organic carbon in the southern Yellow Sea［J］.Advances in Marine Science，2006，24（4）：477－488.賀志鵬，宋金明，張乃星.南黃海溶解有機物的生物地球化學(xué)特征分析［J］.海洋科學(xué)進(jìn)展，2006，24（4）：477－488.

［68］ BERELSON，W M，BALCH W M，NAJJAR R，et al.Relating estimates of CaCO3production，export，and dissolution in the water column to measurements of CaCO3rain into sediment traps and dissolution on the sea floor：A revised global carbonate budget［J］.Global Biogeochemical Cycles，2007，21：GB1024.

［69］ CAO Z M，DAI M H.Shallow－depth CaCO3dissolution：Evidence from excess calcium in the South China Sea and its export to the Pacific Ocean［J］.Global Biogeochemical Cycles，2011，25（2）：GB2019.

［70］ GUO Z G，YANG Z S.Vertical transfer processes for the particulate carbon in the East China Sea［J］.Oceanologia et Limnologia Sinica，1997，28（6）：659－664.郭志剛，楊作升.東海顆粒碳的垂直轉(zhuǎn)移過程［J］.海洋與湖沼，1997，28（6）：659－664.

［71］ ZHANG N X，SONG J M，HE Z P.Biogeochemical mechanism of particulate organic carbon（POC）variations in sea waters［J］.Acta Ecologica Sinica，2006，26（7）：2328－2339.張乃星，宋金明，賀志鵬.海水顆粒有機碳（POC）變化的生物地球化學(xué)機制［J］.生態(tài)學(xué)報，2006，26（7）：2328－2339.

［72］ SONG J M.Biogeochemistry of China Seas［M］.Jinan：Shandong Science ＆ Technology Press，2004.宋金明.中國近海生物地球化學(xué)［M］.濟南：山東科技出版社，2004.

［73］ ZHAO L.GUO X Y.Influence of cross－shelf water transport on nutrients and phytoplankton in the East China Sea：a model study［J］.Ocean Science，2011，7（1）：27－43.

［74］ KANTHA L H.A general ecosystem model for applications to primary productivity and carbon cycle studies in the global oceans［J］.O－cean Modelling，2004，6（3－4）：285－334.

［75］ DAI M H，YIN Z Q，MENG F F，et al.Spatial distribution of riverine DOC inputs to the ocean：an updated global synthesis［J］.Current Opinion in Environmental Sustainability，2012，4：170－178.

［76］ LIU X C，SHEN H T，HUANG Q H.Concentration variation and flux estimation of dissolved inorganic nutrient from the changjiang river into its estuary［J］.Oceanologia et Limnologia Sinica，2002，33（3）：332－340.劉新成，沈煥庭，黃清輝.長江入河口區(qū)生源要素的濃度變化及通量估算［J］.海洋與湖沼，2002，33（3）：332－340.