我注意到中国地质大学(北京)王成善教授的一次谈话[1],其中提到中生代晚期发生过多次的气候波动变化,出现时冷时热的现象。这一点用通常的理论难以解释。不过从最近发表的有关恐龙排放甲烷气体可能会产生温室效应的论文来看。这种气温的波动应该跟恐龙的活动有关。
我尝试着对晚期恐龙生存的微分方程进行了改进[2],就是除了将重金属含量作为灭绝参数A来计算以外,另外再增加一个参数用来反映温室效应对恐龙数量增长的抑制。
当然温室效应是否会对恐龙数量的增长产生负面作用还存在疑问,但温室效应达到了一定程度以后确实会对恐龙产生危害。另外同重金属产生的负面作用相比,温室效应消失的也比较快,在恐龙数量急剧减少以后,就会出现温室气体被植物吸收并最终消失的现象。
考虑到这些因素以后,这里设重金属含量用A1表示,温室气体含量用A2表示,最终可以得到下面的三个方程:
在上述方程中e1比较小,反映出重金属消失得比较慢。而e2则比较大,反应出温室气体消失得比较快。
方程(3)还可以增加一个受制于资源限制,而导致恐龙数量不可能无限增长的一个项。如我在博文“地球自转轴倾角多大才适合恐龙生存”中的方程(6)中所引入的树高受限制的项一样[3]。这就跟马尔萨斯方程一致了。
通过前面的数值解的分析比较,该方程可能存在一个逐渐衰减的振荡解。该振荡解对应了晚中生代地球温度波动的现象。由于今天时间比较紧张,具体的计算结果明天再给出。
参考文献:
[1] 揭秘远古温室地球——访中国地质大学(北京)王成善教授[EB/OL]. [2012-05-11].
http://www.cugb.edu.cn/readRecord.action?ID=1571.
[2] 科学网—晚期恐龙生存的微分方程模型 - 程智的博文[EB/OL]. [2012-05-09].
http://blog.sciencenet.cn/blog-361477-567775.html.
[3] 科学网—地球自转轴倾角多大才适合恐龙生存 - 程智的博文[EB/OL]. [2012-05-11].
http://blog.sciencenet.cn/blog-361477-565218.html.
https://blog.sciencenet.cn/blog-361477-569598.html
上一篇:
民科的完胜下一篇:
生存微分方程的数值解