Introduction
Paleogeography describes the distribution of highlands, lowlands, shallow seas, and deep ocean basins on the Earth's surface during different geological periods. As such Paleogeography combines datasets and interpretations from widely different fields of Earth, Environment and Life sciences that are integrated in paleogeographic reconstruction maps produced since the early days of geology until today (e.g., Hunt, 1837; Scotese, 2021). Paleogeographic reconstructions provide fundamental constraints for climate and geodynamic models (e.g. Gurnis et al., 2012; Chandra et al., 2021), for the exploration of resources (e.g. Markwick, 2019; Dutkiewicz et al., 2020; Wrobel-Daveau et al., 2022) and are widely used to convey our understanding of Earth evolution in outreach and teaching at all levels (e.g. Seton et al., 2023). Paleogeographic reconstructions sharply developed in accuracy following the revolution of the plate tectonic theory and the growing demand in resources related to the industrial revolution, particularly oil and gas (Wegener, 1912, Du Toit, 1937; Wilson 1966; Dewey and Bird, 1970; Ziegler et al., 1985). In the current context of climate change, the interest in climate models aiming to predict future conditions by making hindcasts of the present day and pre-industrial climate is booming (e.g., Lunt et al., 2021). Paleoclimatic and palaeoceanography studies using numerical model simulations particularly depend on the geographic boundary conditions that are implemented in the form of paleo-Digital Elevation Models (paleo-DEMs), describing the land-sea and elevation distribution on the Earth's surface at a geological time under investigation (e.g., Tardif et al., 2021; Straume et al., 2022). The accuracy of these models has increased again tremendously in the last decade in relation to numerical power and the inclusion of complex physical processes, yet their reliability is still strongly dependent on the accuracy of paleogeographic reconstructions because land-sea distributions, bathymetric and topographic patterns considerably determine climatic and biotic evolution and their interactions (Fluteau and Sepulchre, 2021). Despite their importance for Earth System Sciences, paleogeographic reconstructions integrating a wide range of geophysical and geological data from various fields remain complex to do and therefore can hardly keep up with the growing amount of data and models available to update them.
The ongoing development of Geographic Information System (GIS) software that provides advanced utilities for storing and working with geospatial data has allowed scientists to not only visualize their discoveries and interpretations in past coordinates but also build precise models of the motion of tectonic plates, hence reconstructing their position in the geological past (Müller et al., 2016; Scotese, 2016; Stampfli and Borel, 2002; Vérard, 2021; Wrobel-Daveau et al., 2022). Today, the most advanced open-access software package for tectonic reconstruction is GPlates (Williams et al., 2012), which is fully featured to model tectonic plate motions and reconstruct the past position of continents. It also allows to model the deformation at plate boundaries and provide estimates on shortening of the Earth crust and uplift of mountain ranges and plateaux (Müller et al., 2019). GPlates has advanced input - output features that allow rapid integration of vast arrays of geological and geophysical data into a single tectonic model that enables users to explore past plate tectonic evolution that have resulted in complex present-day tectonic configuration. However, additional tools beyond reconstructing the past positions of tectonic blocks, are required to make paleographic reconstructions. The past topography on continents, the past shorelines, the gaps and overlaps arising between the rotated tectonic blocks or the parts of the Earth crust that once existed but were destroyed during geological processes (i.e., subduction, crustal shortening etc.) remain to be reconstructed in the final paleogeographic map. Different developing approaches have been recently described to reconstruct complete paleogeographic maps using the rotated tectonic blocks as an input (e.g., Baatsen et al., 2016; Golonka et al., 2006; Poblete et al., 2021; Scotese, 2014a, Scotese, 2016, Scotese, 2021; Seton et al., 2020; Stampfli and Borel, 2002; Torsvik et al., 2021; Vérard et al., 2011, Vérard et al., 2015, 2021; Wrobel-Daveau et al., 2022; Zahirovic et al., 2019). These excellent advanced approaches remain often difficult to access because they can require private licenses and rely on complex geologic and geophysical expertise in combination with different GIS, graphics software and scripting languages.
To make paleogeographic reconstructions more accessible to a growing community in various disciplines we have developed Terra Antiqua, a free user-friendly software for producing simple regional and global paleo-DEMs. Terra Antiqua takes as input the tectonic model reconstructed using GPlates to make paleogeographic reconstructions using a set of tools. These include combining various raster and vector files (e.g. continents, modern topography, bathymetry), defining shorelines or modifying and creating topo/bathymetric features. Terra Antiqua is a plugin for QGIS - a GIS software, that provides all utilities for input, output and manipulation of geospatial data. Within QGIS, it is possible to interactively add data, modify symbologies, arrange data into layers, and create and modify interpretations. Terra Antiqua is also geared with an intuitive Graphical User Interface (GUI) and guiding tips that make reconstruction easy and fast. The outcome of Terra Antiqua can then be saved in any GIS supported format, including NetCDF that is widely used in climate models to introduce boundary conditions. Terra Antiqua comes with guidelines, user manuals, downloads, video tutorials and a forum on www.paleoenvironment.eu/terra-antiqua and https://jaminzoda.github.io/terra-antiqua-documentation/.
Here we review paleogeographic methods, concepts and software that provide a framework for the user-driven approach used in Terra Antiqua. We then describe the workflow of the Terra Antiqua software for paleogeographic reconstructions and map creation. This is then illustrated by an example application to the global paleogeographic reconstructions at 30 and 50 Ma. Finally, we review the potential developments of paleogeographic reconstructions with GIS in general and for Terra Antiqua in particular.
Section snippets
A review of paleogeographic reconstruction approaches leading to Terra Antiqua
Paleogeographic reconstruction methods have evolved to integrate an increasing range of tools and knowledge from data and models providing more accurate constraints on the temporal and spatial evolution of the shape of the Earth surface. We divided the presentation of the integration of these constraints in two main parts following the general procedure of most paleogeographic reconstruction approaches including Terra Antiqua. In the first part, we present the definition of a plate tectonic
Plate modeling software
Scotese (1976) made one of the first reconstruction programs for plates that used Euler's theorem and was named CALCOMP. It was written in FORTRAN for the IBM 370 and served as a basis for many other modeling programs that were made in the years after, including Terra Mobilis (1985) for Macintosh, Plate Trekker (1994) for Windows, PaleoGIS (2000) for ArcGIS and PointTracker (2001–2010) for iMac (Vérard, 2019). Unfortunately, there is not much information available on these tools because some of
Compilation of data sources
Paleogeographic reconstruction will require to merge various sets of raster data into a single DEM. Terra Antiqua enables doing this using the “Compile Topo/Bathymetry” utility. In Terra Antiqua’s workflow the first step to start a reconstruction is to compile the available topography and bathymetry data sources in a single model. This utility is designed in a way that it takes any number of raster files as an input so it can be used to combine all kinds of raster data. The data is inputted in
Application to global reconstructions at 30 and 50 Ma
As an illustration, we use work in progress on global paleo-reconstructions being produced with Terra Antiqua in the frame of a MSc’s thesis (Ruiz, 2020). The process used, as a starting point, modern topography cut along continental block boundaries that were rotated to 50 Ma and 30 Ma according to a tectonic model by Seton et al. (2012) that was modified by Baatsen et al. (2016) and Poblete et al. (2021). These are combined with the bathymetry at 30 and 50 Ma obtained from Müller et al.
Potential developments paleogeographic reconstructions with Terra Antiqua
Terra Antiqua is designed for a broader scientific community to enable easy access to making paleogeographic maps in formats depending on the user's need and application. It is still currently rudimentary, including only some of the many refined tools that have been developed for paleogeographic reconstructions. However, it provides a basic open access framework on which existing and new tools can be implemented easily for user-friendly applications. As such it has the potential to evolve to
Conclusions
With the QGIS plugin Terra Antiqua we aimed to build a software for paleogeographic reconstructions with a user-friendly graphical user interface available to a broader scientific community. To facilitate the paleogeographic reconstructions, it is an intuitive software that utilizes existing GIS APIs and makes the reconstruction process easy and fast. It enables to make high-resolution paleogeographic maps that are becoming an increasingly important mean for the visualization of geological data
Uncited references
Chamberlain et al., 2021
Deng and Ding, 2015
Dupont-Nivet et al., 2010
Eagles and Hoang, 2014
Li et al., 2008
Mulch et al., 2018
NOAA National Geophysical Data Center, 2009
Passey and Ji, 2019
Peters and McClennen, 2015
Spicer et al., 2018
The Paleobiology Database [WWW Document], 2020
Westerweel et al., 2019
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgements
This project was supported by the ERC grant MAGIC [grant number 649081 to Guillaume Dupont-Nivet], a PIFI postdoctoral fellowship from Chinese Academy of Sciences [grant number 2019PM0039 to Jovid Aminov] and a DAAD research grant (to Jovid Aminov). We are grateful to Dr. Douwe van Hinsbergen for his editorial work on this manuscript. We also appreciate the constructive feedback from two anonymous reviewers that helped us to enhance the quality of this manuscript.
© 2023 Elsevier B.V. All rights reserved.