任湘湘，李海，吳輝碇

（國家海洋環(huán)境預(yù)報中心，北京 100081）

海洋生態(tài)系統(tǒng)動力學(xué)模型研究進(jìn)展

任湘湘，李海，吳輝碇

（國家海洋環(huán)境預(yù)報中心，北京 100081）

海洋生態(tài)系統(tǒng)動力學(xué)模型作為定量地認(rèn)識和分析海洋生態(tài)系統(tǒng)現(xiàn)象的有力工具，近年來得到了長足發(fā)展。本文首先回顧了海洋生態(tài)動力學(xué)模型的發(fā)展歷史，著重介紹了21世紀(jì)以來生態(tài)系統(tǒng)動力學(xué)模型的三大發(fā)展趨勢：一是進(jìn)一步探索海洋生態(tài)系統(tǒng)復(fù)雜性，二是全球氣候變化與海洋生態(tài)系統(tǒng)的相互作用；三是不再局限于理論研究，而進(jìn)入于災(zāi)害預(yù)報與評估、公共決策等應(yīng)用領(lǐng)域。其次介紹了海洋生態(tài)動力學(xué)模型的分類及典型海洋生態(tài)動力學(xué)數(shù)值模型COHERENS的特點、功能和最新的應(yīng)用情況。最后總結(jié)歸納了目前海洋生態(tài)動力學(xué)模型研究領(lǐng)域的幾大問題與挑戰(zhàn)，展望了該研究領(lǐng)域未來的發(fā)展趨勢和方向。

海洋生態(tài)系統(tǒng)動力學(xué)；模型；研究進(jìn)展

1 海洋生態(tài)系統(tǒng)動力學(xué)模型發(fā)展簡史

海洋生態(tài)系統(tǒng)動力學(xué)的發(fā)展過程自始至終與現(xiàn)代海洋學(xué)的發(fā)展過程交織在一起，密不可分，這是由海洋學(xué)的特性所決定的。海洋學(xué)從建立之初便是一門物理、生物、化學(xué)和地質(zhì)等學(xué)科相互交叉的綜合性應(yīng)用科學(xué)，正如MOLL所指出的，海洋環(huán)境要素及其相互之間的關(guān)系錯綜復(fù)雜，一方面隨著我們對海洋環(huán)境認(rèn)識的不斷深入，解決這些難題的必要性也日益提升，但另外一方面隨著我們所要解決的環(huán)境問題日趨復(fù)雜，我們通過現(xiàn)場實驗等手段調(diào)查海洋環(huán)境的能力卻相對下降，在這種形勢下發(fā)展生態(tài)系統(tǒng)動力學(xué)模型成為對復(fù)雜海洋環(huán)境進(jìn)行研究的一種重要手段[1]。

上世紀(jì)40年代末，隨著競爭和捕食—被捕食關(guān)系等新生態(tài)理論的誕生，物理海洋學(xué)家和海洋生態(tài)學(xué)家嘗試將數(shù)學(xué)模型引入海洋生態(tài)學(xué)領(lǐng)域，從此開始了第一代海洋生態(tài)動力學(xué)模型的研究工作。1949年Riley與Stommel等人建立一個垂直一維數(shù)學(xué)模型來描述歐洲北海浮游生物的季節(jié)變化[2]，標(biāo)志著海洋生態(tài)學(xué)從定性描述進(jìn)入了定量模擬的時代。此后直到上世紀(jì)80年代前半期，海洋生態(tài)動力學(xué)模型雖然發(fā)展迅速，但在理論和方法上并無質(zhì)的飛躍。這一時期的生態(tài)動力學(xué)模型，均以垂直一維水柱生態(tài)模型為基礎(chǔ)，可被視為Riley模型的延續(xù)[3]。

1986年前后，美國和日本的科學(xué)家開始了三維海洋生態(tài)動力學(xué)模型的研究工作。Kishi（1986）建立了一個瀨戶內(nèi)海三維有限差分海洋生態(tài)動力學(xué)模型，Walsh（1988）建立了北大西洋海灣三維海洋生態(tài)動力學(xué)模型[1]。

進(jìn)入20世紀(jì)90年代，一系列國際海洋研究計劃陸續(xù)啟動。熱帶海洋與全球大氣計劃（TOGA2-Tropical Ocean Global Atmosphere 2）、海岸帶陸海相互作用（LOICZ-Land-Ocean Interaction in the Coastal Zone）和全球海洋觀測系統(tǒng)（GOOS-Global Ocean Observing System）等計劃的實施大大地推動了海洋生態(tài)系統(tǒng)動力學(xué)的發(fā)展[4]。1989年1月美國推出了全球變化研究計劃(USGCRP-United States Global Change Research Program)，其中全球海洋生態(tài)系統(tǒng)動力學(xué)研究計劃(GLOBEC-GLOBal Ocean ECosystems Dynamics)于1995年被納入國際地圈生物圈計劃(IGBP-International Geosphere-Biosphere Program)的核心計劃，海洋生態(tài)系統(tǒng)動力學(xué)研究成為海洋科學(xué)的前沿和熱點。在這一背景下，圍繞歐洲北海誕生了一系列各具特點的模型。1993年Skogen等將三維斜壓水動力學(xué)模型POM(Princeton Ocean Model)與生物—化學(xué)模型耦合建立了NORWECOM(NORWegian ECOlogical Model system)。1994年Delhez等人將一個簡單的生態(tài)模型與一個中尺度環(huán)流模型耦合建立了GHER(GeoHydrodynamics and Environmental Research)模型。1995年MOLL等人在原有一維水柱生態(tài)模型的基礎(chǔ)上耦合三維斜壓水動力學(xué)模型HAMSOM(Hamburg Shelf Ocean Model)，建立了三維生態(tài)動力學(xué)模型ECOHAM(Ecological North Sea Model Hamburg)。1999年Luyten等人將物理模型、生物模型、沉積物模型、污染物輸運模型耦合在一起建立了多功能的COHERENS(COupled HydrodynamicalEcologicalmodelforREgioNal Shelf seas)模型[1]。我們可以看到這一時期的海洋生態(tài)動力學(xué)模型往往由物理—化學(xué)—生物模型耦合而成，其對海水溫度、鹽度、流場等物理要素的刻畫日臻成熟和細(xì)致，但對生態(tài)系統(tǒng)本身的發(fā)展演變做了很大程度的簡化，可以籠統(tǒng)地概括為營養(yǎng)鹽-浮游植物-浮游動物-碎屑（NPZD-Nutrient-Phytoplankton-Zooplankton-DetritusModel）模型。營養(yǎng)鹽方面考慮氮、磷的濃度及浮游植物生長所受的營養(yǎng)鹽限制，生物方面則只考慮浮游植物、浮游動物、細(xì)菌三類[5]。但無論如何，海洋生態(tài)模型進(jìn)入了強調(diào)“動力學(xué)”的時代，可以被稱為是真正“系統(tǒng)”意義上的海洋生態(tài)動力學(xué)模型。

20世紀(jì)90年代后期至今，海洋生態(tài)動力學(xué)模型的研究出現(xiàn)了三大新趨勢。

第一個新趨勢是對海洋生態(tài)系統(tǒng)的復(fù)雜性有了更多考慮。一些模型對引入了Fe，Si等營養(yǎng)鹽的限制（Fe是微量元素），將浮游植物分為不同的粒級，浮游動物也分為大型和小型，從而引入不同的攝食方式[5]，有的模型甚至涵蓋了魚類等更高營養(yǎng)級，更完整地再現(xiàn)了海洋食物鏈。例如Armstrong建立的模型就側(cè)重于評價Fe在浮游植物生長過程中的限制作用[6]。Chai等建立的一維模式除了考慮N、Si循環(huán)外還將Fe的限制作用隱含在硅藻生長速率中[7]。Chai等還將CSM（Climate System Models）模型的海洋學(xué)部分與一個一維碳、硅、氮生態(tài)模型CoSINE（Carbon,Silicon,Nitrogen Ecosystem）耦合，模擬試驗將赤道太平洋東部海區(qū)鐵濃度提高后，硅藻生產(chǎn)力的變化情況[8]。Fennel近年來嘗試將簡單的浮游生物動力學(xué)模型和魚類生產(chǎn)模型耦合在一起，經(jīng)驗證，該模型可以對魚類捕獲量的年際變化進(jìn)行模擬[9—10]。

第二個新趨勢是從全球氣候變化的角度，關(guān)注CO2濃度變化對海洋生態(tài)系統(tǒng)的影響，以及海洋生態(tài)系統(tǒng)對全球氣候變化的響應(yīng)[11]。Leonard等用一個含有Fe限制的模型很好地再現(xiàn)了1990—1994年Enso事件期間浮游植物豐度下降與群落結(jié)構(gòu)的變化[12]。Skogen和MOLL分別使用NORWECOM和ECOHAM對1984—1994年間歐洲北海年平均初級生產(chǎn)力進(jìn)行了長期模擬，模擬結(jié)果顯示歐洲北海十年間年平均初級生產(chǎn)力在120—140 gCm-2y-1間波動,沒有顯著變化[1]。

第三個新趨勢是突破理論研究的局限，應(yīng)用于災(zāi)害預(yù)報與評估、公共決策制定等領(lǐng)域。近來Stow和Lynch等人則進(jìn)一步提出了對物理-生物耦合模型的精確度進(jìn)行評估的方法，使得海洋生態(tài)動力學(xué)模型可應(yīng)用于預(yù)報，并最終服務(wù)于環(huán)境評估和公共決策[13—14]。Dyble等人則用物理-生物耦合模型來模擬感染病原體和有害藻類地生長、輸運和毒性，從而預(yù)測其對食品安全和人類健康的影響[15]。

我國的海洋生態(tài)學(xué)研究始于1958年的全國海洋普查，但直到上世紀(jì)80年代，該領(lǐng)域的研究工作仍是以現(xiàn)場調(diào)查為主，如1981—1985年的“渤海水域漁業(yè)資源、生態(tài)環(huán)境及其增殖潛力的調(diào)查研究”，1985—1989年的“黃海大海洋生態(tài)系統(tǒng)調(diào)查”，這一系列的調(diào)查工作為我國海洋生態(tài)動力學(xué)的發(fā)展積累了寶貴的物理、化學(xué)、生物資料[4]。

我國真正意義上的海洋生態(tài)動力學(xué)模型研究起步于上世紀(jì)90年代后期。隨著國家自然科學(xué)基金重大項目“渤海生態(tài)系統(tǒng)動力學(xué)與生物資源可持續(xù)利用研究”以及國家重點基礎(chǔ)研究發(fā)展規(guī)劃（973）項目“東、黃海生態(tài)系統(tǒng)動力學(xué)與生物資源可持續(xù)利用”等一系列重大科研項目的帶動，我國學(xué)者在海洋生態(tài)學(xué)現(xiàn)場調(diào)查、過程和模型研究中做了大量的工作[16]。

近年來，俞光耀（1999）[17—18]、高會旺（1999）[19]、崔茂常（2000）[20—21]、張書文（2002）[22]、魏皓（2003）[23—24]等人在海洋生態(tài)動力學(xué)建模領(lǐng)域進(jìn)行了一系列嘗試，對渤海、黃海區(qū)域的營養(yǎng)鹽分布、浮游植物生物量、初級生產(chǎn)力等生態(tài)系統(tǒng)要素進(jìn)行了較成功的模擬。夏潔等（2006）將POM模式與ERSEM（European Regional Seas Ecosystem Model）耦合建立垂直一維生態(tài)動力學(xué)模型對南黃海東部營養(yǎng)鹽、葉綠素濃度的季節(jié)變化進(jìn)行了模擬[25]；劉浩（2006）等建立了一個與POM耦合同時考慮了太陽輻射和河流輸入的NPZD生態(tài)模型并將其應(yīng)用于渤海，驗證了渤海f-ratio特征及硝酸鹽、銨鹽對渤海浮游植物生長的相對貢獻(xiàn)率[26]；劉浩等（2008）運用物理-生物耦合模型對渤海的初級生產(chǎn)力、浮游植物生物量和N、P的季節(jié)變化特征進(jìn)行了模擬[27]。

特別值得一提的是，曾慶存院士等在2008年發(fā)表的《地球系統(tǒng)動力學(xué)模式及模擬研究》一文中將海洋生物地球化學(xué)模式作為地球系統(tǒng)動力學(xué)模式的重要組成部分進(jìn)行了介紹[28]。這種將地球上的大氣、海洋、巖石和生態(tài)各圈視作一個有機整體進(jìn)行研究的觀點正是當(dāng)前地球科學(xué)發(fā)展的主要趨勢。我國海洋生態(tài)動力學(xué)模型的研究已經(jīng)開始從全球氣候變化的高度關(guān)注海洋生態(tài)系統(tǒng)在地球動力系統(tǒng)碳循環(huán)中的作用[29—30]，這標(biāo)志著我國的海洋生態(tài)動力學(xué)模型研究進(jìn)入了一個新的階段。

2 海洋生態(tài)系統(tǒng)動力學(xué)數(shù)值模型分類

根據(jù)科學(xué)目的、研究對象、研究區(qū)域、數(shù)值求解方法等不同標(biāo)準(zhǔn)，海洋生態(tài)動力學(xué)模型有許多分類方法。

2.1 根據(jù)研究范圍不同，海洋生態(tài)動力學(xué)模型可分為：

（1）全球或海盆尺度海洋生態(tài)系統(tǒng)動力學(xué)模型目前主流的全球或海盆尺度海洋生態(tài)系統(tǒng)動力學(xué)模型大多是生物—物理耦合的三維模型，建立在海洋動力學(xué)環(huán)流模型的基礎(chǔ)之上，并力圖全面地納入各種生物—化學(xué)變化過程。但必須承認(rèn)目前的全球或海盆尺度海洋生態(tài)動力學(xué)模型仍停留在定性描述階段。近年來，人們逐漸認(rèn)識到大尺度的氣候變化對全球生態(tài)系統(tǒng)有著深刻影響，而生態(tài)系統(tǒng)也以其特有的方式反作用于全球氣候系統(tǒng)。但目前地球科學(xué)界尚未揭示出生態(tài)系統(tǒng)與氣候變化之間相互作用的根本機制。因此從全球氣候變化的意義上探索海洋生態(tài)系統(tǒng)和全球氣候變化之間的關(guān)系是許多全球或大尺度海洋生態(tài)系統(tǒng)動力學(xué)數(shù)值模型重點關(guān)注的問題，特別是人為產(chǎn)生的CO2對氣候變化和海洋生態(tài)系統(tǒng)帶來的影響是研究的熱點。

（2）近岸區(qū)域海洋、海灣生態(tài)系統(tǒng)動力學(xué)模型迄今為止，歐洲北海及其某些海灣是海洋生態(tài)動力學(xué)模型研究最頻繁的區(qū)域，并由此誕生了NORWECOM、 GHER、 ECOHAM、 ERSEM、COHERENS、 POL3dERSEM（Proudman Oceanographic Laboratory 3d European Regional Seas Ecosystem Model）等一系列各具特色的模型。上世紀(jì)70—80年代，近岸海洋生態(tài)動力學(xué)模型的研究重點是富營養(yǎng)化條件下的近岸生態(tài)系統(tǒng)。但近年來，隨著研究的不斷深入和環(huán)境、資源管理要求的不斷升級，模型除了關(guān)注浮游生態(tài)系統(tǒng)、初級生產(chǎn)力模擬之外，近岸海洋生態(tài)動力學(xué)模型也開始應(yīng)用于漁業(yè)產(chǎn)量測算、近岸環(huán)境變化預(yù)測等領(lǐng)域[3]。Franks和陳長勝等人將海灣近岸海洋模型ECOM-si與NPZ生態(tài)模型耦合對喬治淺灘初級生產(chǎn)力和春季浮游植物水華等一系列問題進(jìn)行了研究，并取得了一些符合實際情況的研究結(jié)論[31—33]。

（3）河口生態(tài)系統(tǒng)動力學(xué)模型

河口地區(qū)由于自身環(huán)境的特殊性，其流場結(jié)構(gòu)、水質(zhì)、生物種群結(jié)構(gòu)等均具有與其他海區(qū)明顯不同的特點，故此形成其特有的河口地區(qū)生態(tài)系統(tǒng)。GEMBASE是具有代表性的河口生態(tài)動力學(xué)模型之一。Chen等人利用生物—物理耦合模型模擬了河流淡水輸入對于美國路易斯安娜與得克薩斯近岸生態(tài)系統(tǒng)的影響[34]。

2.2 根據(jù)研究對象不同，海洋生態(tài)動力學(xué)模型可分為

2.2.1 過程模型

過程模型一般針對特定的生物過程（如攝食、種群競爭等），描述生物個體的生理參數(shù)與生物功能之間的關(guān)系（如攝食率與餌料濃度的線性關(guān)系，代謝率與溫度的關(guān)系等），并以線性或非線性回歸方程表示。它能夠較好地反映數(shù)據(jù)的統(tǒng)計特征，但難以刻畫生物過程的動力學(xué)規(guī)律[4]。過程模型是建立生態(tài)系統(tǒng)模型的基礎(chǔ)。

2.2.2 個體模型、種群模型和種間模型

個體模型通常用來描述生物個體與環(huán)境之間的物質(zhì)能量交換和生物個體本身的生長、發(fā)育并最終死亡的整個過程。種群模型是通過描述種群密度、出生率、死亡率、年齡分布、生長形式等特征屬性來刻畫某一種群的形態(tài)、生長、發(fā)展和變化。種間模型則是建立在種群模型的基礎(chǔ)之上，除了刻畫幾個不同種群的自身特征屬性之外，同時描述了不同種群之間的競爭、捕食等錯綜復(fù)雜的關(guān)系。Chen等利用海灣近岸海洋模型ECOM-si（Estuarine,Coastal and Ocean Model semiimplicit）與個體模型IBMs的耦合模型模擬了喬治淺灘潮流混合對橈足類的生長繁殖過程的影響[35]。He等人將區(qū)域海洋模型ROMS與一個種群動力學(xué)模型耦合對2005年緬因灣西部亞歷山大藻暴發(fā)進(jìn)行了后報模擬[36]。

2.2.3 系統(tǒng)模型

系統(tǒng)模型是針對整個系統(tǒng)，在一定時空尺度內(nèi)所建立的模型，著眼于模擬整個生態(tài)系統(tǒng)的行為[37]。系統(tǒng)模型往往包含物理、化學(xué)、生物等多種變量，或著眼于刻畫海洋生態(tài)系統(tǒng)中的物質(zhì)和能量流動。系統(tǒng)模型有的重在模擬大洋或局地海洋生態(tài)系統(tǒng)的總結(jié)構(gòu)，刻畫整個系統(tǒng)內(nèi)的物質(zhì)循環(huán)、能量傳遞及穩(wěn)態(tài)調(diào)節(jié)機制；有的重在研究生態(tài)系統(tǒng)中物理和生物要素之間的相互作用關(guān)系[38]。近年來,系統(tǒng)模型已經(jīng)成為海洋生態(tài)動力學(xué)模型的主流，NORWECOM、ECOHAM、COHERENS、ERSEM等各有特點的模型均為系統(tǒng)模型。

2.3 根據(jù)模型空間結(jié)構(gòu)不同，海洋生態(tài)動力學(xué)模型可分為：

2.3.1 箱式模型

將某特定區(qū)域內(nèi)的水體視為若干個相互連接的箱體，箱內(nèi)所有生態(tài)變量都是均勻的，僅為時間的函數(shù)，不同的箱之間以及箱與外界環(huán)境之間存在物質(zhì)交換，根據(jù)這種思路建立起來的生態(tài)動力學(xué)模型為箱式模型[39]。其中最有代表性的是歐洲區(qū)域海生態(tài)模型ERSEM，它將歐洲北海表層分為85個箱，底層分為45個箱。歐洲各國科學(xué)家利用該模型在歐洲北海生態(tài)學(xué)研究方面做了大量卓有成效的工作。Varela等（1995）用ERSEM計算了北海的總初級生產(chǎn)力[1]； Bryant等人(1995)利用ERSEM模擬了北海北部魚類的繁殖、生長和種群變遷[1]； Patsch等（1997）模擬再現(xiàn)了歐洲北海1955—1993年間的海富營養(yǎng)化現(xiàn)象[1]。Lenhart（1997）用改進(jìn)的ERSEMⅡ箱式模型模擬了河流輸入的營養(yǎng)鹽對歐洲北海營養(yǎng)鹽和初級生產(chǎn)力空間分布的影響[1]。Allen等人（2001）將ERSEM箱式生態(tài)動力學(xué)模型嵌入到一個三維環(huán)流模型，建立了三維生態(tài)動力學(xué)模型POL3dERSEM[1]。夏潔將ERSEM與一維POM模式進(jìn)行耦合，模擬了南黃海東部浮游生態(tài)系統(tǒng)的季節(jié)變化規(guī)律[20]。如今將箱式模型嵌入到物理模型中已成為一種常見的做法。

2.3.2 一維模型

垂直一維生態(tài)動力學(xué)模型不考慮水平輸運對生態(tài)系統(tǒng)的影響，因此一般適用于生態(tài)系統(tǒng)水平方向變化不顯著的大洋區(qū)域，對描述生態(tài)變量年際變化非常有效。1949年Riley開發(fā)的模型就是垂直一維模型。1985年Radach等人則利用垂直一維模型模擬了北海中部25年間（1962—1986年）生產(chǎn)力的年際變化[1]。

2.3.3 三維模型

三維海洋生態(tài)系統(tǒng)動力學(xué)模型對大氣、海洋等物理背景的模擬更加精確，能夠反映出生態(tài)系統(tǒng)的三維空間分布特征和隨時間演變的特征，NORWECOM，ECOHAM，COHERENS,FVCOM都是三維海洋生態(tài)動力學(xué)模型的經(jīng)典代表。

3 典型的三維海洋生態(tài)系統(tǒng)動力學(xué)模型——COHERENS

COHERENS（Coupled Hydrodynamical Ecological model for REgionNal Shelf seas）是1990—1998年間由 NOMADS（National Model Archive and Distribution System）等一系列歐盟計劃共同資助建立的三維多功能沿岸和陸架海模型。它主要由物理模型、生物模型、沉積物模型和污染物輸運模型耦合而成。

COHERENS的物理模型以POM模式為基礎(chǔ),具有以下基本特點：

垂直方向采用σ坐標(biāo)；

水平采用曲線正交坐標(biāo)和Arakawa C網(wǎng)格差分方案；

水平時間差分是顯式格式而垂直差分是隱式格式；

模式包含二階湍流閉合子模式提供垂直混合系數(shù)；

模式具有自由表面，時間積分采用分裂解法。模式外模是二維的，根據(jù)CFL條件和外波速采用短時間步長；內(nèi)模是三維的，根據(jù)CFL條件和內(nèi)波速采用長時間步長；

包含完備的熱力學(xué)過程；

采用了Boussinesq近似和準(zhǔn)靜力近似。

COHERENS的生物模型是以Tett（1990）提出的“微生物艙”為基礎(chǔ)建立的，共包含8個狀態(tài)變量：微型浮游植物碳、微型浮游植物氮、碎屑碳、碎屑氮、硝酸鹽、銨鹽、溶解氧和浮游動物氮。生物模型主要考慮了微型浮游生物部分和碎屑部分之間的有機碳和有機氮的循環(huán)過程，此外還考慮了溶解硝酸鹽、銨鹽和溶解氧濃度的變化。通過求解這些狀態(tài)變量的輸送擴散方程得到它們隨時間變化的三維分布，在方程中將這些生物過程作為輸運方程的源匯項，同時考慮了垂直沉降以及平流、擴散產(chǎn)生的物理輸送。生物模型結(jié)構(gòu)見圖1。

COHERENS自誕生以來在海洋環(huán)境模擬和海洋生態(tài)系統(tǒng)模擬方面得到了廣泛的應(yīng)用。Lee等（2002）將其應(yīng)用在歐洲北海北部，并與之前用ERSEM、ECOHAM1模擬的結(jié)果進(jìn)行了比較。COHERENS模型只考慮了生物圈氮循環(huán)，而ECOHAM1模型則只考慮了磷循環(huán)，ERSEM不但包含氮、磷、硅循環(huán)還納入了底棲動物。ERSEM在生物部分自有其優(yōu)勢，但COHERENS模型的物理模塊更加成熟。鑒于COHERENS自帶的生態(tài)模塊未將硅藻考慮在內(nèi)，Geneviève Lacroix等（2005）將COHERENS與生態(tài)模型MIRO[40]耦合建立三維海洋生態(tài)動力學(xué)模型MIRO&CO-3D對1991—2003年間北海南部海灣的無機碳、有機碳、營養(yǎng)鹽、浮游植物、浮游動物、細(xì)菌的年變化進(jìn)行了模擬。在時間序列上，模型能夠反映營養(yǎng)鹽和浮游植物的季節(jié)變化規(guī)律，但模擬的春季葉綠素a平均濃度偏低，硅循環(huán)的模擬則失敗了，因為模型模擬的春季硅藻暴發(fā)比實際延遲了；在空間分布上，冬季海表面營養(yǎng)鹽濃度特別是PO4模擬得比較成功，但萊茵河口的模擬結(jié)果偏低，法國海岸的模擬結(jié)果則偏高。陳秀華等（2007）以COHERENS為基礎(chǔ)建立了σ坐標(biāo)下的三維渤黃東海夏季環(huán)流模型并在此基礎(chǔ)上進(jìn)行了拉格朗日顆粒追蹤和歐拉輸運數(shù)值模擬來推測長江口及其臨近海域赤潮高發(fā)區(qū)的可能赤潮源[41—42]。Liang B.C.等（2006）將COHERENS與SWAN耦合后應(yīng)用于渤海，模擬了海浪對表面風(fēng)壓和表面混合高度的影響[43]。李艷蕓等（2006）運用COHERENS建立了渤海天文潮預(yù)報模型,模擬了9712號熱帶風(fēng)暴作用下渤海的增水過程[44]。

圖1 COHERENS模式生態(tài)模型結(jié)構(gòu)圖[34]

4 總結(jié)與展望

綜上所述，上世紀(jì)90年代特別是進(jìn)入新世紀(jì)以來，海洋生態(tài)動力學(xué)模型取得了長足的發(fā)展。對海洋水動力環(huán)境的描述日益細(xì)化，例如對垂直方向的湍流描述多采用湍流閉合方案[45]。生態(tài)系統(tǒng)基本上分為營養(yǎng)鹽、溶解有機物、浮游植物、浮游動物、細(xì)菌等幾大單元，過程上則注重描述營養(yǎng)鹽和碳的循環(huán)過程，近年來更是將鐵、硅的營養(yǎng)鹽限制作用還有不同粒級浮游生物的營養(yǎng)和攝食考慮進(jìn)來[5]。

但必須承認(rèn)，海洋生態(tài)動力學(xué)模型經(jīng)過幾十年的發(fā)展仍有一些問題未能得到解決。

（1）觀測資料的匱乏一直是困擾世界海洋生態(tài)學(xué)界的一個難題。魏皓曾指出在過去的生態(tài)模型中，缺乏針對特定海域的生物化學(xué)過程參數(shù)[38]，張素香則認(rèn)為參數(shù)的缺乏使得模擬結(jié)果的準(zhǔn)確性難以提高[46]。

（2）模型對海洋生態(tài)系統(tǒng)特別是更高營養(yǎng)級的描述仍然不夠全面也不夠準(zhǔn)確。目前的海洋生態(tài)動力學(xué)模型有兩個分支，一是重點描述生物化學(xué)過程的浮游生物生態(tài)動力學(xué)模型，魚類一般作為“捕食壓力”隱含其中；二是魚類生產(chǎn)模型，以浮游生物為代表的低營養(yǎng)級一般作為“食物”隱含其中[9—10]。浮游生物生態(tài)動力學(xué)模型的發(fā)展雖相對成熟，但仍不能全面準(zhǔn)確地模擬許多化學(xué)、生物變化過程。而同時涵蓋浮游生物和魚類的海洋生態(tài)動力學(xué)模型仍停留在起步階段，劉桂梅認(rèn)為目前在多營養(yǎng)級模型建立方面仍然存在較大不足[4]。

（3）模型的時間、空間尺度難以選擇。傳統(tǒng)的浮游生物動力學(xué)模型其物理變化過程的時間和空間尺度往往比化學(xué)、生物變化過程的尺度大得多，這成為制約該類型模型發(fā)展的一個因素。商少凌指出，尺度不同的物理過程和生物過程的模擬不能匹配是造成模型結(jié)果與觀測不符的重要原因之一[5]。而隨著涵蓋魚類等更高營養(yǎng)級的海洋生態(tài)動力學(xué)模型不斷發(fā)展，必須將魚類的生命周期長（往往超過一年），遷徙范圍廣（覆蓋整個大洋）等因素考慮在內(nèi)，因此物理模型的空間尺度則應(yīng)當(dāng)選擇海盆級，而時間尺度則應(yīng)當(dāng)選擇跨年或十年[47]。因此考慮到整個海洋生態(tài)系統(tǒng)中物理、化學(xué)、生物等各因素的特性和相互作用，海洋生態(tài)動力學(xué)模型的時空尺度選擇確實是一項挑戰(zhàn)。

（4）全球尺度或者海盆尺度的模型只能進(jìn)行趨勢描述而無法進(jìn)行精確模擬。目前的區(qū)域性海洋生態(tài)動力學(xué)模型無法拓展到全球范圍，而全球尺度模型的精度往往很難提高，這些都與觀測資料的匱乏和物理模型的限制是分不開的。

針對這些問題和挑戰(zhàn)，人們對海洋生態(tài)動力學(xué)模型的未來發(fā)展有以下幾點考慮：

（1）在發(fā)展模型的同時注重實測資料的獲取和處理

資料匱乏是貫穿海洋生態(tài)動力學(xué)模型發(fā)展過程的老問題。王輝曾指出在我國目前的生態(tài)動力學(xué)模型發(fā)展水平下，應(yīng)針對生態(tài)系統(tǒng)的特點設(shè)計現(xiàn)場調(diào)查方案，采取先進(jìn)的取樣和分析手段[37]。同時還應(yīng)該發(fā)展參數(shù)優(yōu)化技術(shù)等資料處理方法，使得現(xiàn)有資料獲得最大限度的利用，為模型的發(fā)展提供基礎(chǔ)。

（2）更加全面、準(zhǔn)確地描述化學(xué)、生物變化過程

海洋生態(tài)系統(tǒng)中物理、化學(xué)、生物變化過程及其相互作用的錯綜復(fù)雜是公認(rèn)的。高會旺曾指出，上世紀(jì)90年代以后的海洋生態(tài)動力學(xué)模型狀態(tài)變量已經(jīng)從2—4個發(fā)展到3—7個[45]?？梢灶A(yù)見，如果狀態(tài)變量太少，則無法全面反映海洋生態(tài)系統(tǒng)的復(fù)雜性，而動輒十幾個，甚至幾十個狀態(tài)變量也必將嚴(yán)重影響模型的計算速度。因此在全面真實反映海洋生態(tài)系統(tǒng)狀況和計算效率之間，研究者必須取得平衡。

（3）涵蓋海洋食物網(wǎng)中的更高營養(yǎng)級

多營養(yǎng)級的海洋生態(tài)動力學(xué)模型研究雖然剛剛起步，但已是大勢所趨，浮游生物生態(tài)動力學(xué)模型和魚類生產(chǎn)模型都不能完整地模擬整個海洋食物網(wǎng)，故未來海洋生態(tài)動力學(xué)模型的發(fā)展將瞄準(zhǔn)多營養(yǎng)級的食物網(wǎng)動力學(xué)并體現(xiàn)生物多樣性[9]。

（4）全球變化意義上的海洋生態(tài)動力學(xué)模型應(yīng)當(dāng)對中尺度時空變化過程給予充分重視[5][37]。

海洋生態(tài)學(xué)界非常重視全球氣候變化和海洋生態(tài)系統(tǒng)之間的相互作用，全球變化意義上的海洋生態(tài)動力學(xué)模型發(fā)展是大勢所趨，但目前全球或海盆尺度的海洋生態(tài)動力學(xué)模型仍停留在定性描述的階段，對中尺度時空變化過程的忽略可能是一個重要原因。中尺度渦旋，熱帶非穩(wěn)定波等中尺度問題的解決將有可能使得全球尺度的海洋生態(tài)動力學(xué)模型進(jìn)入定量模擬的時代。

[1]Moll A,Radach G.Review of three-dimensional ecological modeling related to the North Sea Shelf system Part1:models and their results[J].Progress in Oceanography,2003,57(2003):175-217.

[2]Riley G A,Stommel H,Bumpus D F.Quantitative ecology of plankton of Western North atlantic[J].Bull Bingham Oceanogr Coll,12:1,169.

[3]王海黎,洪華生.海洋生態(tài)動力學(xué)模式[J].海洋科學(xué),1996,2:16-18.

[4]劉桂梅,孫松,王輝.海洋生態(tài)系統(tǒng)動力學(xué)模型及其研究進(jìn)展[J].地球科學(xué)進(jìn)展,2003,18(3):427-432.

[5]商少凌,柴扉,洪華生.海洋生物地球化學(xué)模式研究進(jìn)展[J].地球科學(xué)進(jìn)展,2004,19(4),621-629.

[6]Armstrong R A.An optimization-based model of iron-light-ammonium colimitation of nitrate uptake and phytoplankton growth[J].Limnology Oceanography,1999,44:1436-1446.

[7]Chai F,Dugdale R C,Peng T H,et al.One dimensional ecosystem model of the equatorial Pacific Upwelling system.PartⅠ:Model development and silicon and nitrogen cycle[J].Deep Sea ResearchⅡ,2002,49:2713-2734.

[8]Chai F,Jiang M S,Chao Y,et al.Modeling responses of diatom productivity and biogenic silica export to iron environment in the equatorial Pacific Ocean[J].GLOBAL BIOGEOCHEMICAL CYCLES, 2007, VOL.21, GB3S90, doi: 10.1029/2006GB002804.

[9]Fennel W.Parameterization of truncated food web models from the perspective of an end to end model approach[J].Journal of Marine Systems,2009,76(1-2):171-185.

[10]Fennel W.Towards bridging biogeochemical and fish-production models[J].Journal of Marine Systems,2008,71(1-2):171-194.

[11]Harley C D G,Hughes A R,Hultgren K M,et al.The impacts of climate change in coastal marine systems[J].Ecology Letters,2006,9:228-241.

[12]Leonard C L,McClain C R,Murtugudde R,et al.An iron based ecosystem model of central equatorial Pacific[J].Journal of Geophysical Research,1999,104:1325-1341.

[13]Stow C L,Jolliff J,McGillicuddy Jr.D J,et al.Skill assessment for coupled biological/physical models of marine systems.[J].Journal of Marine Systems,2009,76:4-15.

[14]Lynch D R,McGillicuddy Jr.D J,Wenner F E.Preface Skill assessment for coupled biological/physical models of marine systems[J].Journal of Marine Systems,2009,76:1-3.

[15]Dyble J,Bienfang P,Dusek E,et al.Environmental controls,oceanography and population dynamics of pathogens and harmful algal blooms:connecting sources to human exposure[J].Environmental Health,2008,7(S2):S5.

[16]蘇紀(jì)蘭,唐啟升.我國海洋生態(tài)系統(tǒng)基礎(chǔ)研究的發(fā)展——國際趨勢和國內(nèi)需求[J].地球科學(xué)進(jìn)展,2005,20(2):139-143.

[17]余光耀,吳增茂,張志南,等.膠州灣北部水層生態(tài)動力學(xué)模型與模擬Ⅰ.膠州灣北部水層生態(tài)動力學(xué)模型.青島海洋大學(xué)學(xué)報,1999,29(3):421-428.

[18]吳增茂,余光耀,張志南,等.膠州灣北部水層生態(tài)動力學(xué)模型與模擬Ⅱ.膠州灣北部水層生態(tài)動力學(xué)的模擬研究[J].青島海洋大學(xué)學(xué)報,1999,29(3):429-435.

[19]高會旺,王強.1999年渤海浮游植物生物量的數(shù)值模擬[J].中國海洋大學(xué)學(xué)報,2004,34(5):867-873.

[20]Zhu Hai,Cui Maochang.Coupled physical-ecological modeling of the central part of Jiaozhou bayⅠ:Physical modeling[J].Chinese Journal of Oceanography Liminology,2000,18(4):309-314.

[21]Cui Maochang,Zhu Hai.Coupled physical-ecological modeling of the central part of Jiaozhou bayⅡ:Coupled with an ecological model[J].Chinese Journal of Oceanography Liminology,2001,19(1):21-28.

[22]張書文,夏長水,袁業(yè)立.黃海冷水團水域物理－生態(tài)耦合數(shù)值模式研究[J].自然科學(xué)進(jìn)展,2002,12(3):315-319.

[23]魏皓,趙亮,馮士筰.渤海碳循環(huán)與浮游植物動力學(xué)過程研究[J].海洋學(xué)報,2003,25(S2):151-156.

[24] 魏皓,趙亮,馮士筰.渤海浮游植物生物量與初級生產(chǎn)力變化的三維模擬[J].海洋學(xué)報,2003,25(S2):66-72.

[25]夏潔,高會旺.南黃海東部海域浮游生態(tài)系統(tǒng)要素季節(jié)變化的模擬研究[J].安全與環(huán)境學(xué)報，2006,6(4):59-65.

[26]劉浩,尹寶樹.渤海生態(tài)動力學(xué)過程的模型研究Ⅰ.模型描述[J].海洋學(xué)報,2006,28(6):21-31.

[27]劉浩,潘偉然.營養(yǎng)鹽負(fù)荷對浮游植物水華影響的模型研究[J].水科學(xué)進(jìn)展,2008,19(3):345-351.

[28]曾慶存,周廣慶,浦一芬,等.地球系統(tǒng)動力學(xué)模式及模擬研究[J].大氣科學(xué),2008,32(4):653-690.

[29]徐永福,浦一芬,趙亮.海洋碳循環(huán)模式的研究進(jìn)展[J].地球科學(xué)進(jìn)展,2005,20(10):1106-1115.

[30]徐永福,趙亮,李陽春.海洋碳循環(huán)與海洋生態(tài)系統(tǒng)動力學(xué)[J].海洋環(huán)境科學(xué),2007,26(5),495-500.

[31]Franks P J S,Chen C.A 3-D prognostic numerical model study of the Georges Bank in Summer[J].Journal of Marine Research,2001,54:631-651.

[32]Chen C,Beardsley R C,Franks P J S.A 3-D prognostic model study of the ecosystem over Georges Bank and adjacent coastal regions[J].Part 1:Physical model.Deep-Sea Research,2001,48:419-456.

[33]Ji R,x Chen C.Spring phytoplankton bloom and associated lower trophic level food web dynamics on Georges Bank:1-D and 2-D model stydies[J].Deep-Sea ResearchⅡ,2006,53(23-24):2656-2683.

[34]Chen C,Denis A,Wiesenburg,et al.Influences of river discharge on biological production in the inner shelf:A coupled biological and physical model of the Louisiana-Texas Shelf[J].Journal of Marine Resesrch,1997,55(2):293-320.

[35]陳長勝.海洋生態(tài)系統(tǒng)動力學(xué)與模型[M].北京:高等教育出版社,2003.

[36]Ruoying He,Dennis J,McGillicuddy J,et al.Historic 2005 toxic bloom of Alexandrium fundyense in the western Gulf of Maine:2 Coupled biophysical numerical modeling[J].Journal of Geophysical Research,2008,113(C7):

[37]王輝.海洋生態(tài)系統(tǒng)模型研究的幾個基本問題[J].海洋與湖沼,1998,29(4):341-346.

[38]魏皓,趙亮,武建平.浮游植物動力學(xué)模型及其在海域富營養(yǎng)化研究中的應(yīng)用[J].地球科學(xué)進(jìn)展,2001,16(2):220-225.

[39]Luyten P J,Jones J E,Proctor R,et al.1999.COHERENS-Acoupled hydrodynamical-ecological model for regional and shelf seas:User Documentation[M].Belgium:MUMM Report,Management Unit of the Mathematical Models of the north Sea,1999.

[40]Lacroix G,Ruddick K,Park Y,et al.Validation of the 3D biogeochemical model MIRO&CO with field nutrient and phytoplankton data and MERIS-derived surface chlorophyll a images[J].Journal of Marine Systems,2007,64:66-88.

[41]Chen X H.Numerical simulation of summer circulation in the East China Sea and its application in estimating the sources of red tides in the Yangtze River estuary and adjacent sea Areas[J].Journal of Hydrodynamics,2007,19(3):272-281.

[42] 陳秀華,朱良生.COHERENS模式在長江口赤潮源推測中的應(yīng)用[J].海洋通報,2007,26(3):3-11.

[43]Liang B C,Li H J,Lee D Y.Numerical study of wave effects on surface wind stress and surface mixing length by three-dimensional circulation modeling[J].Journal of Hydrodynamics,2006,18:97-404.

[44]李艷蕓,李紹武,風(fēng)暴潮預(yù)報模式在渤海海域中的應(yīng)用研究[J],海洋技術(shù),2006,25(1):101-106.

[45]高會旺,馮士筰,管玉平.海洋浮游生態(tài)系統(tǒng)動力學(xué)模式的研究[J].海洋與湖沼,2000,31(3):341-347.

[46]張素香,李瑞杰,羅鋒,等.海洋生態(tài)動力學(xué)模型的研究進(jìn)展[J].海洋湖沼通報,2006,4:121-127.

[47]Young B,Heath M,Werner F.et al.Chanllenges of Modeling Ocean Basin Ecosystems[J].SCIENCE,2004,304:1463-1466.

Areview of marine ecosystem dynamics model

REN Xiang-xiang，LI Hai，WU Hui-ding

(National Marine Environmental Forecasting Center,Beijing 100081 China)

With increasing pressure for a profound understanding of marine ecosystems,numerical modeling becomes a powerful tool for the research.The development of marine ecosystem dynamics model in the last decades is reviewed.In general,marine ecosystem dynamics model could be classified into several different categories according to various features.Typical model COHERENS（COupled Hydrodynamical Ecological model for REgioNal Shelf seas）is introduced in this article.The latest study focussed on marine ecosystem dynamics model are summarized,for example,interaction between marine ecosystem and global climate change,biological approach models including higher trophic levels,application of ecosystem models in forecasting and public policy.Finally,issues and challenges in the marine ecosystem model in the near future are also discussed.

marine ecosystem dynamics，numerical model，research progress

Q178

A

1003-0239（2012）01-0065-08

2011-01-27

國家自然科學(xué)基金“南海初級生產(chǎn)對季風(fēng)變動響應(yīng)的數(shù)值模擬”（40531006）;科技部"赤潮預(yù)警預(yù)報、應(yīng)急及損害評估技術(shù)"（2007AA092003-01）;國家科技支撐項目“重大海洋災(zāi)害預(yù)警及應(yīng)急技術(shù)研究”（2006BAC03B-01-05）；國家自然科學(xué)基金“南海海洋生態(tài)系統(tǒng)碳循環(huán)的數(shù)值模擬研究”（41076011）

任湘湘(1981-)，女，助理研究員，主要從事物理海洋和海洋生態(tài)學(xué)模型研究。E-mail:rxx@nmefc.gov.cn