GEN3SIS: An engine for simulating eco-evolutionary processes in the context of plate tectonics and deep-time climate variations

Explaining the origin of large-scale biodiversity gradients has been a key aspiration of early naturalists such as Wegener, Darwin and Humboldt; who looked at natural processes in an integrated way. Early on, these naturalists acknowledged the role of plate tectonics and climate variations in shaping modern day biodiversity patterns. 

As science advanced, the complexity of ecological, evolutionary, geological and climatological processes became evident while research became increasingly fragmented across different disciplines. Nevertheless, recent development in mechanistic modeling approaches now enable bringing disciplines back together, opening a new interdisciplinary scientific pathway.

Here, we present GEN3SIS, the GENeral Engine for Eco-Evolutionary SImulationS. It is the first spatially explicit eco-evolutionary model that incorporates deep-time Earth history, including plate tectonics, as well as climate variations in a modular way. The modular design allows exploring the consequences of user-defined biological processes that act across “real world” spatio-temporal landscapes. Emerging from the model are specie’s ranges, alpha and beta diversity patterns, ecological traits as well as phylogenies. Subsequently, these patterns can be compared to empirical data. Furthermore, GEN3SIS allows assessing paleoclimatic and paleogeographic hypotheses by using different Earth history scenarios and comparing simulation outputs with empirical biological data.

As a case study, we explore the cold-adapted plant biodiversity dynamics throughout the Earth’s Cenozoic history, based on a deep-time tectonic and climate reconstruction. The Cenozoic India-Asia collision formed the Himalayan mountain range. In this highly elevated region, the first cold niches of the Cenozoic appeared, demanding adaptation from the local living flora. We hindcast diversification of cold-adapted species with GEN3SIS, for which we use a topo-climatic reconstruction for the last 55 Myr. The model predicts the emergence of current cold-species richness patterns. Moreover, simulations indicate that cold-adapted flora emerged in the Oligocene, first in the Himalayas, followed by a spread to the Arctic. This agrees with observed low species richness and high nestedness of Arctic assemblages compared to those of the Himalayan mountain ranges. Under ongoing climate change a major loss of cold-adapted plant diversity is expected by the end of the century, particularly in lower latitude mountain ranges. Hindcasting and forecasting dynamics of cold-adapted lineages highlights the transient fate of cold organisms in a warming world.

GEN3SIS is made available as an R package, which allows customizing (i) the simulated landscape including environmental variables and (ii) all the processes interacting under different spatial and temporal scales. Consequently, GEN3SIS fosters collaborations between different natural disciplines and therefore contributes to an interdisciplinary understanding of the processes that shaped Earth’s history.