基于系統動力學的雛菊世界模型氣候控制敏感性分析

陳海濱，唐海萍

(1.地表過程與資源生態國家重點實驗室，北京 100875;2.西北農林科技大學經濟管理學院，楊凌 712100;3.北京師范大學資源學院，北京 100875)

英國大氣物理學家J.E.Lovelock于20世紀60年代受NASA邀請，幫助尋找火星上的生命痕跡。在任務完成之后，他繼續思考生命維持所必需的條件:如地球表面的溫度經歷了36億年后仍然保持相對不變，而同一時期的太陽溫度卻升高了25%;地球是如何能夠使大氣中氧氣濃度保持在21%左右，這個值是他認為的生命安全上限。Lovelock認為解決這樣的問題需要采用一個系統的觀點，而不能由生物化學或生物地球化學等學科來解答，也不能由新達爾文主義生物學來解答，答案來自于生理學或控制論［1］。

Lovelock對于Walter Cannon描述的有關產生體內平衡的協同生理過程印象深刻。Cannon認為，生命體努力維持其內部環境恒定不變［2］。Lovelock相信地球也遵循著Cannon的原理:巖石、草、鳥類、海洋和大氣都齊力協力，像一個巨型的有機體一樣運行來調控生存條件［3］。地球這個“超級生命有機體”需要一個名字。Lovelock的鄰居，諾貝爾文學獎獲得者William Golding，推薦以古希臘大地女神命名，稱之為蓋婭［1］。

提出伊始，蓋婭假說的“超級生命有機體”這樣略帶神秘的概念遭到了科學保守主義者，尤其是新達爾文主義者的強烈批評，認為它是偽科學，是目的論。例如，Richard Dawkins認為，“蓋婭不能產生于達爾文式的生命進化，地球作為物理實體不是自然選擇的單元”［4］。為了反擊對蓋婭假說的批評，Lovelock于1983年提出了雛菊世界模型［5-6］。Lovelock用簡單的微分方程就揭示了，假想的世界不需要一個清晰的、公認的、可以描述其內部相互作用的數學方程就可以擁有自我平衡的性質。

在雛菊世界模型中，Lovelock論證了光照強度不斷增強時星球溫度的變化情況。本文在此基礎上，嘗試用系統動力學方法進一步分析雛菊世界的氣候條件對于不同初始面積分配、不同反射率組合條件的敏感性。一方面可以驗證蓋婭假說，另一方面通過探討模型存在的不足，為下一步的模型改善以及蓋婭假說的發展提供一些參考。

1 方法原理

1.1 蓋婭假說

Lenton和Wilkinson認為，蓋婭假說的發展經歷了經典的3個階段:一是“它明顯是錯誤的”;二是“它可能有些是正確的”;三是“我們都知道它是正確的”［7］。

第一階段 蓋婭假說認為，地球有機體能調節地球系統的物理、化學成分，以使星球成為生命最優的棲息地［1，3］。這個階段的蓋婭假說是目的論的，因為他沒有提出一個機制來解釋，個體生物之間似乎需要某種秘密的共識來達成行星水平的環境自我調節［8］。Lovelock在其后來的文獻中也承認:“20世紀70年代引入的蓋婭假說認為，地球的大氣、海洋、氣候和地殼是由生物調節處于穩定狀態，并適合于生物生存。特別地，蓋婭假說認為巖石和水的溫度、氧化狀態、酸性和其他某些方面在任何時候都保持恒定，這個自我平衡狀態是由地球表面的有機體維持的。應該認識到，蓋婭假說如此闡述是錯誤的”［9］。

第二階段 為了應對批評，Lovelock引入雛菊世界模型。模型的模擬結果表明，雖然太陽輻射增強了近兩倍，但是行星表面溫度卻能夠維持在雛菊生長最優溫度的幾度范圍之內。行星水平的溫度自我調節可以產生于假定的生物個體水平的溫度-增長率關系，蓋婭假說不是目的論［5-6］。但是，這仍然無法令進化生物學家滿意。

第三階段 這階段的蓋婭假說在批評與證明中不斷完善，新版本的蓋婭假說定義為“地球系統的物理、化學和生物組成部分可以調節星球，以維持其作為生命的棲息地”［9］。生命有機體對氣候系統有重要影響已經不再有爭議［10-11］。Lenton在Nature上發文指出:“物理系統是穩定的，生物系統可不斷自我增殖，因此存在一個生物生長的物理最優狀態;行星自我調節的反饋機制產生于生物體的自然選擇”［12］。而進化生物學家也承認蓋婭假說值得嚴肅思考，因為它可以解釋很多重大問題［13-16］。

當然，爭議依然存在。主要原因在于復雜多變的地球系統充滿了不確定性，有很多案例并不支持生物自我調節地球環境的觀點［17］。Lovelock在《蓋婭的復仇》一書中引用William Hamilton的話說:“正如哥白尼的發現需要牛頓來解釋一樣，我們也需要另一個牛頓來解釋為什么達爾文的自然選擇會導致一個宜居的星球”［18］。無論如何，蓋婭假說改變了科學界對自然、地球的觀點，并推動了地球系統科學的發展。通過對蓋婭假說的爭論，促進了地質學、大氣物理學、氣象學、生物學等相關學科的發展，同時也促進了多學科的交叉、融合。

1.2 雛菊世界模型

為了證明蓋婭假說，Lovelock提出了一個雛菊世界的寓言。雛菊世界是一個假想的世界，它圍繞著一顆輻射能量緩慢增大的恒星做軌道運動，它意圖模擬太陽-地球系統。在1983年版的雛菊世界中，只種植著兩種雛菊以代表生命形式——白雛菊和黑雛菊。白雛菊和黑雛菊通過反射較多和較少太陽光，改變著各自的局地溫度，后者反過來又決定了雛菊的增長率。該模型經過模擬，追蹤了當太陽輻射緩慢增強時兩種雛菊的種群動態以及雛菊世界的溫度變化。結果證明，在太陽輻射變化的相當大范圍內(從當前輻射水平的0.6倍到1.1倍)，雛菊世界的溫度保持幾乎恒定［5-6］。在后來的雛菊世界模擬中，加入兔子、狐貍和其他物種。結果表明，物種數量越多，整個行星的穩定性越強;也就是說，溫度調節能力越強;甚至當行星受到干擾時，系統也可以保持強健和穩定［19］。這些發現都支持了生物多樣性與地球系統穩定性存在正向關系的觀點。

雛菊世界模型的優點在于它對紛繁復雜的地球系統的簡化，使其成為回答“生物與環境相互作用會怎樣”這類問題的杰出模型［20］。雖然它無法解釋地球內部是如何自組織的，但是它卻提供了地球系統中生物與環境緊密耦合的一般經驗。在雛菊世界模型出現以后的近30年時間里，很多文獻努力通過模型的演變來證明或反對蓋婭假說［21］。進化生物學家通過引入特殊的案例來證明在真實的世界里面有些物種在干擾或破壞自我調節作用;而支持者通過引入其他案例證明地球自我調節作用的強健與穩定;當然也有很多理論研究關注的是自我調節實現與否的關鍵條件是什么。

當前，雛菊世界模型已經擴展到眾多學科領域，包括非線性動態、生態系統與食物網、進化論、生理學、最大熵產生以及人造生命等。雖然仍然與蓋婭理論緊密相關，但雛菊世界模型已經發展成為一個獨立的模型［21］。因此，在討論雛菊世界模型與蓋婭假說的關系之外，該模型也可以用于發展整合人類活動的地球系統模型的研究。

1.3 系統動力學與Stella軟件概述

系統動力學是20世紀50年代末由美國麻省理工學院教授福雷斯特(J.W.Forrester)提出［22-23］。它是一種分析研究復雜反饋系統動態行為的科學方法，是系統科學的一個分支，也是一門溝通自然科學和社會科學領域的橫向學科，實質上就是分析研究復雜反饋大系統的計算仿真方法。目前，它已經廣泛應用于工業、經濟、生態、環境等諸多學科領域［24］。

生態系統也是一個復雜的動態系統，因此，許多生態學者認為應用系統動力學進行生態學研究將有極大潛力，并積極開展了種群數量動態、大氣組成變化、土地利用與土地覆蓋變化等方面的推廣應用工作，各類生態系統動力學模型應運而生。

Stella軟件是美國Isee公司于是1985年開發的系統模擬工具。因其圖形界面十分友好，在國外已成為一個構造系統模型和模擬復雜系統動態相互關系的良好工具，被廣泛應用于科研、教學、管理等諸多研究領域［25］。

2 模型假設與參數設定

2.1 模型假設

關于生物與環境條件的相互作用，雛菊世界模型假設:

(1)生物影響環境條件 如上所述，由黑、白雛菊、裸地分布的星球影響了光的吸收率，因此影響了黑、白雛菊的局地溫度和星球的溫度。

(2)環境條件反作用于生物 局地溫度影響雛菊生長的速率。當局地溫度接近于最適溫度時，雛菊增長率最高;但是，如果局地溫度太高或太低，雛菊覆蓋的區域就會隨時間減少。黑、白雛菊的年增長率與局地溫度的關系表示為(1)式(圖1)。

式中，R為年增長率，T為雛菊局地溫度。

圖1 黑、白雛菊年增長率-溫度關系圖Fig.1 Annual growth rate-temperature relationship of black and white daisy

2.2 參數設定

(1)星球總表面積為1000英畝;

(2)雛菊衰敗速率為30%/a;

(3)為了簡化模擬，假設到達地面的太陽輻射值恒定，為太陽常數值的70%，約等于947.1W/m2;

(4)星球的反射率為:

式中，A為星球反射率;αb為黑雛菊所占面積比例，Ab為黑雛菊反射率;αg為裸地所占面積比例，Ag為裸地反射率;αw為白雛菊所占面積比例，Aw為白雛菊反射率。

(5)星球整體的太陽光吸收率:

(6)太陽光吸收率與星球平均溫度之間的關系用下式表達:

式中，σ =5.67 ×10－8W m－2K－4，Te為星球平均溫度;SL為到達地面的太陽輻射;A為星球平均反射率。其關系如圖2所示。

圖2 星球平均溫度-太陽光吸收比率關系圖Fig.2 Planet's average temperature-sunlight absorption ratio relationship

(7)由于太陽光吸收率的不同，導致黑雛菊局部升溫，而白雛菊局部降溫:

式中，Tb，w表示黑雛菊或白雛菊的局地溫度;A為星球平均反射率;Ab，w表示黑雛菊或白雛菊的反射率;Te表示星球的平均溫度。

根據Lovelock的模型，q′的大小表示較高溫度區域向較低溫度區域的熱能傳導能力，在這里q′取20℃。

(8)星球的裸地面積對于雛菊的現實增長率有限制作用，以增長率衰減乘數來表示這個限制作用:

增長率衰減乘數=裸地面積/星球總面積

2.3 雛菊世界因果回路圖

根據上述模型假設與參數關系，繪出雛菊世界因果回路圖。圖3中共有8個反饋回路，白雛菊由其中4個反饋回路組合控制，黑雛菊由另外4個反饋回路組合控制。

白雛菊的4個反饋回路分別是與分布面積有關的2個正、負反饋回路;另外白雛菊的增長率受剩余空地面積負反饋回路的影響，同時也受星球平均反射率負反饋回路的影響。星球平均反射率是由黑、白雛菊的反射率與覆蓋面積決定的。

黑雛菊的4個反饋回路中，2個也是與分布面積有關的正、負反饋回路;另外2個是與增長率相關的負反饋回路。

圖3 雛菊世界因果反饋回路圖(改自安德魯·福特［26］)Fig.3 Casuality feedback loop diagram of Daisyworld(revised from Fort，2011［26］)

3 敏感性分析

3.1 初始面積的敏感性分析

假設不同的地面覆蓋物的反射率分別為:白雛菊0.75;裸地0.5;黑雛菊0.25。

在Stella軟件環境下，模擬不同初始面積分配條件下的雛菊世界溫度變化情況，得到結果如圖4所示。結果顯示，在A、B、C三種初始面積組合條件下，雖然起始星球平均溫度不同，但達到穩定狀態后的星球平均溫度均為21.55℃，表明雛菊世界穩定時的星球平均溫度對于初始面積分配不敏感。

3.2 反射率的敏感性分析

圖4 雛菊世界溫度調節初始面積敏感性分析Fig.4 Initial area sensibility analysis of Daisyworld's temperature regulation

在Stella軟件環境下，模擬3組不同地面覆蓋物反射率條件下的雛菊世界溫度變化情況，得到結果如圖5:在A、B、C3種反射率組合條件下，雖然起始星球平均溫度不同，但達到穩定狀態后的星球平均溫度均為21.55℃，表明雛菊世界穩定時的星球溫度對于地面覆蓋物反射率條件不敏感。

值得一提的是，在C反射率組合下，星球平均溫度經過約500次模擬(即約500a)才達到穩定狀態，而其他兩種組合只需要幾十次模擬就可以達到穩定狀態。究其原因，是因為C組合中，黑、白雛菊的反射率與裸地差別不大，導致黑、白雛菊面積變化相對較慢。表明，雖然雛菊世界穩定狀態時的星球溫度對于地面覆蓋物反射率條件不敏感，但是覆蓋物反射率組合卻影響雛菊世界達到穩定狀態的時間。

4 討論

本文基于雛菊世界模型，運用系統動力學的方法，通過stella軟件求解不同初始面積分配與反射率組合條件下星球溫度的穩定值，用簡單且易于操作的方法論證了蓋婭假說。結果證明，雛菊世界的星球溫度對于初始面積分配、光照反射率條件均不敏感，在不同的初始面積分配以及不同的反射率條件下，雛菊世界穩定狀態時最終都達到21.55℃。同時，通過模擬發現，星球覆蓋物反射率組合影響星球達到穩定狀態所需要的時間。

雛菊世界模型的優點在于通過簡單的反饋關系揭示了在沒有人為控制的情況下，地球系統水平的環境條件可以通過生物個體水平的自然選擇與競爭實現自我調節，保持穩定，從而使環境條件有利于生物的群體生長。這對于那些沒有受過系統分析訓練，或不習慣系統思考的人，確實是個很好的啟示，能夠改變人們的思維習慣。因此，雛菊世界模型現在被廣泛地用做教學案例。

然而，相比復雜的真實地球系統，雛菊世界模型還是過于簡單。綜合體現在以下5個方面:

(1)在雛菊世界原型中，太陽光照強度是不斷增強的，但在本研究中，為了簡化模擬，假定太陽輻射值為恒定。因此，本研究沒有模擬太陽輻射強度變化對星球溫度的長期(幾十億年)影響。

(2)只有溫度一個環境變量。現實世界中，對生物有影響的環境變量還有水文、地形、地貌、養分、大氣組成等，這些變量模型都沒有包含進去。

(3)調節環境變量的生物只考慮非常簡單的黑、白雛菊，不包含其他植物、動物、微生物。在現實世界中，黑雛菊很少見，而且單黑、白兩種雛菊很難代表所有植物對光的反映，更不用說用雛菊來代表動物和微生物。而從調節大氣溫度的角度考慮，微生物所起的作用更關鍵。既使修正后的模型把兔子、狐貍等考慮進去，但比起真實地球系統的復雜性，顯然遠遠不夠。

圖5 雛菊世界溫度調節反射率敏感性分析Fig.5 Albedo sensibility analysis of Daisyworld's temperature regulation

(4)生物與環境反饋關系過于簡單。雛菊模型只包含一個溫度-生長的反饋關系，現實世界中則存在著多種環境變量與不同的生物之間的高度復雜的反饋關系。生物的某些行為，例如呼吸、新陳代謝、應激反應、生殖等，顯然與環境變化也有重要關系。另外，生物多樣性與環境調節功能的反饋關系也是未來需要進一步研究的焦點問題。

(5)雛菊模型只考慮生物-環境關系，沒有考慮人為活動的影響。當前，人類活動顯著改變著地球陸地、海洋、大氣、巖石等圈層，并影響生物組成、結構與分布。在這樣的背景下，蓋婭假說如果不考慮人為因素的影響，其得出的結論往往與現實相悖。例如，蓋婭假說預測大氣組成將受生物過程的調節，并保持適合生物生存的水平。但自從工業革命以來，地球大氣CO2濃度已經上升了35%，生物圈吸收CO2的速度卻只提高了2%［27］。上述現象，Lovelock在《蓋婭的復仇》和《蓋婭消失的面孔:最后的警告》中也提到了，他認為“地球即將進入一個長達10萬年的病態發熱期”［18，28］。

因此，如前演示，這樣一個極度簡化的雛菊世界模型必然呈現簡單的系統行為。由于蓋婭假說主要考慮生物-環境的關系，沒有將人類活動考慮進去，最近全球變化的許多事實都在蓋婭假說的討論范圍之外。為了使模型與現實世界更相符，迫切需要設計更為精細的模型，包含更多的環境、生物變量以及反饋關系和人類活動因素，最終能夠將蓋婭系統融入到自然-人類耦合的大系統中。

［1］ Lovelock J E.Gaia，A New Look at Life on Earth.Oxford:Oxford University Press，1979.

［2］ Cannon W B.The Wisdom of the Body.New York:W.W.Norton＆ Company，Inc，1932.

［3］ Lovelock J E，Margulis L.Atmospheric homeostasis by and for the biosphere:the gaia hypothesis.Tellus，1973，26(1/2):2-10.

［4］ Dawkins R.The Extended Phenotype:The Long Reach of the Gene(Popular Science).Revised Edition.Oxford:Oxford University Press，1999:236-236.

［5］ Lovelock J E.Daisy World-a cybernetic proof of the Gaia hypothesis.The Co-evolution Quarterly Summer，1983:66-72.

［6］ Watson A J，Lovelock J E.Biological homeostasis of the global environment:the parable of Daisyworld.Tellus B，1983，35B(4):284-289.

［7］ Lenton T M，Wilkinson D M.Developing the Gaia Theory:a response to the criticisms of Kirchner and Volk.Climatic Change，2003，58(1/2):1-12.

［8］ Doolittle W F.Is Nature Really Motherly?Spring:The CoEvolution Quarterly，1981:58-63.

［9］ Lovelock J E.The Ages of Gaia:A Biography of Our Living Earth，2nded.Oxford:Oxford University Press，1995:213-216.

［10］ Kirchner J W.The Gaia Hypothesis:fact，theory，and wishful thinking.Climatic Change，2002，52(4):391-408.

［11］ Lawton J.Earth system science.Science，2001，292(5524):1965-1965.

［12］ Lenton T M.Gaia and natural selection.Nature，1998，394(6692):439-447.

［13］ Hamilton W D.Ecology in the large:Gaia and Genghis Khan.Journal of Applied Ecology，1995，32(3):451-453.

［14］ Wilkinson D M.Is Gaia really conventional ecology?Oikos，1999，84(3):533-536.

［15］ Cohen J E，Rich A D.Interspecific competition affects temperature stability in daisyworld.Tellus，2000，52(3):980-984.

［16］ Dagg J L.Unconventional bed mates:Gaia and the selfish gene.Oikos，2002，96(1):182-186.

［17］ Kirchner J W.The Gaia Hypothesis:conjectures and refutations.Climatic Change，2003，58(1/2):21-45.

［18］ Lovelock J E.The Revenge of Gaia.New York:Basic Books，2006.

［19］ Lovelock J E.A numerical model for biodiversity.Philosophical Transactions:Biological Sciences，1992，338(1286):383-391.

［20］ Lenton T M，Lovelock J E.Daisyworld revisited:quantifying biological effects on planetary self-regulation.Tellus b，2001，53(3):288-305.

［21］ World A J，Ackland G J，Dyke J G，Williams W T P，Lenton T M.Daisyworld:a review.Reviews of Geophysics，2008，46，RG1001，doi:10.1029/2006RG000217.

［22］ Forrester J W.Industrial dynamics:a major breakthrough for decision makers.Harvard Business Review，1958，36(4):37-66.

［23］ Forrester J W.Industrial Dynamics.Cambridge Mass:MIT Press，1961.

［24］ Wang Q F.System Dynamics，Revised Edition.Beijing:Tsinghua University Press，1994.

［25］ Richmond B.An Introduction to Systems Thinking.Isee Systems，ISBN 0-9704921-1-1.

［26］ Andrew F.Modeling the Environment:An Introduction To System Dynamics Modeling Of Environmental Systems，2nd edition.Translated by Tang H P，Shi P J.Beijing:Science Press，2011.

［27］ Schimel D S，House J I，Hibbard K A，Bousquet P，Ciais P，Peylin P，Braswell B H，Apps M J，Baker D，Bondeau A，Canadell J，Churkina G，Cramer W，Denning A S，Field C B，Friedlingstein P，Goodale C，Heimann M，Houghton R A，Melillo J M，Moore B III，Murdiyarso D，Noble I，Pacala S W，Prentice I C，Raupach M R，Rayner P J，Scholes R J，Steffen W L，Wirth C.Recent patterns and mechanisms of carbon exchange by terrestrial ecosystems.Nature，2001，414(6860):169-172.

［28］ Lovelock J E.The Vanishing Face of Gaia:A Final Warning.New York:Basic Books，2009.

參考文獻:

［24］ 王其潘.系統動力學(修訂版).北京:清華大學出版社，1994.

［26］ 安德魯·福特著.環境模擬-環境系統的系統動力學模型導論(第二版).唐海萍，史培軍譯.北京:科學出版社，2011:265-277.