“天河二号”大力相助 科学家看清古生物多样性演化精彩“瞬间”

在地球历史中，曾经发生过重大的生态系统和环境的突变，导致了多次生物大灭绝事件。利用古生物数据库重建地质历史全球生物多样性模式，揭示地球生命的演化历史，以及与环境变化之间的关系，可以为了解当前的地球生物多样性危机提供重要启示。

最近，南京大学、中国科学院南京地质古生物所的樊隽轩教授、沈树忠院士团队，借助“天河二号”超级计算机，运用大数据、超算、遗传算法等全新的方法和手段，获得了全新的寒武纪—早三叠世海洋无脊椎动物的复合多样性曲线，将其统计时间分辨率从1000万年精细到约2.6万年，更加准确地重现了地质历史中三次生物大灭绝事件和两次重大生物辐射事件的精细过程。近日，该成果在《科学》杂志在线发表。

此前时间分辨率只能精确到800万—1000万年

生命演化中有3个基本问题，即地球上复杂的生态系统是如何形成的?生命如何从一种形式演化到另一种形式?人类出现之前的地球生物圈是否存在灾变和生物多样性事件?

如今，地球上曾经生活过的生物中99%以上已经灭绝，但只有很少一部分能保存为化石。如何通过不完整的化石记录重建地球历史生物多样性的变化规律，是一个重大科学难题。

“上世纪八九十年代，科级、属级的显生宙生物多样性曲线，已经可以将地层剖面中保存的海洋生物化石所反映的不同时代信息，精确到800万—1000万年。”南京大学地球科学与工程学院副教授史宇坤介绍。

但是前人使用的低分辨率且不均一的时间标尺，难以准确重现生命演化的精细过程，会直接影响对古生物多样性的估算，导致无法准确评估生物多样性的变化速率和模式，并可能掩盖突发性的重大事件以及短时间的剧烈波动。

“如果想建立更加精确的时间标尺，看到更小时间段内地球发生了哪些生命演化，就需要完成大量的化石生物多样性的计算。”史宇坤说，当统计的物种数量上升到1万个时，用传统的计算机需要计算17.3年。

超算等技术让分辨率精度提高400倍左右

为了建立约5.4亿年—2.4亿年前即寒武纪至三叠纪早期海洋生物多样性曲线，南京大学教授樊隽轩等收集了大量的地层剖面和化石记录，他们从2013年起，将1.1万个中国化石剖面、2.7万个化石名称、30.5万个化石记录录入数据库系统，并最终从中遴选了3112个地层剖面、11268个海洋化石物种的26万个化石数据进行分析。

如何利用这些海量数据精确、迅速地分析地球生命演化信息?史宇坤介绍，研究团队结合模拟退火算法和遗传算法，自主开发了基于并行计算的约束最优化方法。

“我们利用‘天河二号’超级计算机，反复计算和验证。2017年的高性能超算版本，将原来需要17年才能完成的一次计算，压缩到2—3天，计算速率提高了6000倍。”史宇坤说，这是首次在“天河二号”上运行地层学程序，也是古生物研究领域率先引入大数据和超算算法的成功尝试，并获得巨大突破。

最终，科研团队获得了全新的寒武纪—早三叠世海洋无脊椎动物的复合多样性曲线，将生物多样性的时间分辨率从1000多万年精细到2.6万年，精度提高了400倍左右。

“以往有关化石生物多样性的统计分析的数据很粗，每个时间段大约跨越1000万年，但我们通过大数据、超算、定量化技术，可以精细地看到几万年内发生的生物多样性的变化以及环境的变化，这给当今生物多样性如何受环境和气候的影响，提供了重要的启示。”中国科学院院士、南京大学教授沈树忠说。

准确重现三次生物灭绝、两次生物辐射过程

新建立的多样性变化曲线，更准确地重现了地质历史中最大的三次生物灭绝事件和两次重大生物辐射事件的精细过程。

例如，2.52亿年前，地球发生了人类迄今为止识别出的最大规模的生物灭绝事件，导致约80%的海洋生物在数万年内迅速灭亡。这一事件的发生，与当时全球气候的快速升温密切相关。

此外，在此次研究中，研究人员发现，4.9亿—4.7亿年前的中晚奥陶世和3.4亿—3亿年前的石炭—二叠纪，全球气候变冷，海洋平均温度下降5摄氏度左右，但生物多样性却大幅增加。“这对于我们认识人类活动导致的第六次生物大灭绝以及与全球气候变化之间的关系具有重要科学价值。”沈树忠说。

地质历史中生物多样性的重大变化，通常也伴随着环境的剧烈波动。论文选取了6种与气候变化密切相关的环境指标，包括碳、氧、锶同位素、沉积物质总量、大气二氧化碳含量等。

“我们希望能寻找到在最关键的生物多样性变化节点上，对应的环境变化是多少。”史宇坤说，虽然目前这些环境指标还缺少更精细的时间尺度，但初步的分析表明，大气二氧化碳含量与生物多样性密切相关。“例如，生物多样性的变化与二氧化碳的浓度在早古生代的对应关系是0.7，但在晚古生代是0.52，二氧化碳浓度变化与生物多样性变化，可能由相同的因素导致，但到底是什么，还需要进一步研究。”