面向气候变化的复杂地球系统建模与模拟探索(上).pdfVIP

自然杂志第34卷第3期 特约专稿 0253-9608．2012．03．001 doi：10．3969／j．issn 面向气候变化的复杂地球系统建模与模拟探索(上)* 李琦①查传捷② ①教授，②研究员，北京大学数字地球工作宣，北京100871 *创新方法基础性设计与实施(2008IM05100) 关键词 数字化离散非线性 回递函数 不可逆反应混沌 节能减排 地球是开放系统的集合．在现实世界中。任何取自于地球的模拟对象皆包含着耗散与回递的因子．所以它们亦是 非线性的。非线性方程式基本上无解。这就是建模与数字化面临的难题。因此本文给出系统的严谨数学定义与其参 数．透过自动化模式、线性机械线路、形式语言与其等价的文法规则及混沌系统理论，直接以时空拓扑结构及回递函 数、■子方程式与化学反应式之本质。取代偏微积分方程式，以非线性离散数学直接描述所有取自于地球的模拟对 象，令数字化模式与模拟对象一致，以解决建模问题，更为科学打开通往数学物理化的门扉。进而从根本上消弭数学 与物理间不一致的基础危机。 Loranz)开 1960年混沌理论先锋劳伦斯(Edward 事实上。问题的关键不在于数据的搜罗与累积，而 启了数字化电脑天气模拟风气之先，劳伦斯经过一番苦 是数字的“本质”。答案也应该是直接根据自然律的代 数性建立系统模式，而不是拐弯抹角去寻找方程式。 思，以量子力学的普朗克(MaxPlanck)黑体幅射方程 式、渡与粒子二重性、质能不灭定理、盖·吕萨克定律 Law)与查理定律(Charles’sLaw)阐述 (Gay．Lussac’s 1给出系统的定义 温度、湿度及风速与压力之间的关系“。]。劳伦斯选择 了12个偏微积分“原始”方程式将它们电脑化，科学家 凡是系统皆可由一组七个参数(因子)定义，如下： 们声称这些方程式理论上都没有问题，因为几个世纪下 S21，O，Q，正，，0，，q，L，其中，I={fo，i-， 来证明这些方程式都很管用。只不过时代变化得太快， i：。…，i．，i。，…}是输人信号的集合；O={O。，O，， 一旦资料超过期限，预测就无法准确了。 O：。…，O．，O。，…}是输出的集合，又成为下一个状态 历年来动力学家们普遍相信，只要找出系统的偏微 的输入，此即回递系统的特质，亦是系统呈现出混沌之肇 积分方程式就能解决问题，但这些方程式总会带着一点 因；Q2{鼋o，毋，啦，…tq。，目…，…}是系统状态的集 非线性。非线性系统是一种对起始条件非常敏感的东 合．状态的集合小至量子态，大至整个宇宙，当然也包括 西，解这类公式就像在走迷宫，每走出一步后迷宫的墙 地球以及地球中所有的原子与分子。 壁就会重组。这些方程式经过电脑化无可避免的简化 凡系统皆有自己的输入、输出、状态及三组回递函 后，设计师们努力让模拟的结果与现况不至于太离谱。 数：，i，^，厶。这三组回递函数分别对应于上述I，O， 一旦二者差异太大，譬如：2011年3月11日日本东北太Q三个集合，三组回递函数皆包含着耗散与回递的因 平洋近海发生M9．0级的地震，引发海啸；2010年英格子，因此是非线性、离散函数，它们是自然系统回馈机制 兰发生千年大雪，西北方伍德弗德(Woodford)南方班森下源白物理、化学与生化的反应式。将这三组回递函 (Benson)气温，在一夕间骤降至一18℃，改变了地球环数，化为数学语言，如下： 境与动力学的背景，电脑设计师们就会尽力修正相关的 输入的集合：，=，I(i，q)，由即刻之输入信号i． i。∈IV 资料库及附加一些子程式，直到使模拟的结果与事实吻 及目前系统状态q。决定．也就是：，={V