A data-driven framework for paleomagnetic Euler pole analysis

Owing to the inherent axial symmetry of the Earth’s magnetic field, paleomagnetic data only directly record the latitudinal and azimuthal positions of crustal blocks in the past, but paleolongitude cannot be constrained. An ability to overcome this obstacle is fundamental to paleogeographic reconstruction. The paleomagnetic Euler pole (PEP) analysis presents a unique means to recover such information in deep-time. However, prior applications of the PEP method have invariably incorporated subjective decisions into its execution, undercutting its fidelity and rigor. Here we present a data-driven approach to PEP analysis that addresses some of these deficiencies---namely the objective identification of change-points and small-circle arcs that together approximate an apparent polar wander path. We elaborate on our novel methodology and conduct some experiments with synthetic data to demonstrate its performance. We furthermore present implementations of our methods both as adaptable, stand-alone scripts and as a streamlined interactive workflow that can be operated through a web browser.

U-Pb speleothem geochronology reveals a major 6 Ma uplift phase along the western margin of Dead Sea Transform

The timing of vertical motions adjacent to the Dead Sea Transform plate boundary is not yet firmly established. We utilize laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U-Pb geochronology of carbonate cave deposits (speleothems) to constrain paleo-groundwater levels along the western margin of the Dead Sea Transform and provide a proxy for the timing of large-scale incision and tectonic uplift. Phreatic speleothems can form in caves that are located slightly below the groundwater level. Tectonic uplift and/or base level subsidence can trigger incision of canyons and induce a drop in the groundwater table. This can cause dewatering of the caves, cessation of the deposition of phreatic speleothems, and initiation of growth of vadose speleothems. The transition between deposition of phreatic and vadose speleothems can therefore reflect tectonic or erosive events. We obtained 102 U-Pb ages from 32 speleothems collected from three cave complexes across a 150-km-long, north-to-south transect. These ages indicate that phreatic deposition began between 14.68 ± 1.33 and 11.34 ± 1.62and ended by 6.21 ± 0.59 Ma. Later, vadose speleothems grew intermittently until the Quaternary. These results suggest an abrupt drop in the water table starting at ca. 6 Ma with no re-submergence of the caves. We interpret this to indicate river incision of ∼150−200 m that was driven by uplift and folding of the western margin of the Dead Sea Transform and by inland morpho-tectonic, base-level subsidence in the Dead Sea area. The observed timing corresponds with a change in the Euler pole of the plates motion along the Dead Sea Transform. The growth period of phreatic speleothems suggests groundwater level stability and limited vertical tectonic motions between 14 Ma and 6 Ma.

Kinematics of Deformable Blocks: Application to the Opening of the Tyrrhenian Basin and the Formation of the Apennine Chain

We describe the opening of back-arc basins and the associated formation of accretionary wedges through the application of techniques of deformable plate kinematics. These methods have proven to be suitable to describe complex tectonic processes, such as those that are observed along the Africa–Europe collision belt. In the central Mediterranean area, these processes result from the passive subduction of the lithosphere belonging to the Alpine Tethys and Ionian Ocean. In particular, we focus on the opening of the Tyrrhenian basin and the contemporary formation of the Apennine chain. We divide the area of the Apennine Chain and the Tyrrhenian basin into deformable polygons that are identified on the basis of sets of extensional structures that are coherent with unique Euler pole grids. The boundaries between these polygons coincide with large tectonic lineaments that characterize the Tyrrhenian–Apennine area. The tectonic style along these structures reflects the variability of relative velocity vectors between two adjacent blocks. The deformation of tectonic elements is accomplished, allowing different rotation velocities of lines that compose these blocks about the same stable stage poles. The angular velocities of extension are determined on the basis of the stratigraphic records of syn-rift sequences, while the rotation angles are obtained by crustal balancing.

Tectonic escape of Sicily Microplate in the context of Africa-Europe collision

<p>Sicily is in the centre of an area where complex geodynamic processes work together, these are: the Tyrrhenian-Apennine System evolution, the African-Ionian slab subduction and Africa-Europe collision.</p><p>During the last 5 Ma it was involved in a process of escape towards east-southeast: while on one side Africa acted as an intender pushing toward north, on the other side the fragmentation and retreat of the African-Ionian slab created space to the east.</p><p>The aim of this study is to reconstruct the kinematic evolution of Sicily, here considered as an independent plate starting from 5 Ma ago, and its role in the context of the Tyrrhenian-Apennine system.</p><p>The plates and microplate involved in the evolution are Europe, Africa and Calabria. The boundaries between these and Sicily are the margin of the Sicily microplate and are lithospheric structures known from the literature and identifiable from high resolution bathymetric maps, seismic sections, geodetic data, focal mechanism of recent earthquakes, gravimetric maps, lithosphere thickness maps and so on.</p><p>Briefly the margin between Sicily and Europe is along the Elimi chain, a E-W trending morpho-structure with transpressive kinematics, the margin with Calabria microplate is along the right-lateral Taormina line and the margin with Africa is expressed along the Malta Escarpment, south of Etna Mount, with transpressive kinematics and along the Sicily Channel, where a series of troughs (Pantelleria, Linosa and Malta) were interpreted in literature as pull-apart basins related to a dextral trascurrent zone.</p><p>The Euler pole of rotation between Sicily and Africa was found starting from the structures in the Sicily Channel and using the GPlates software, then we were able to find also Sicily-Europe and Sicily-Calabria poles and the respective velocity vectors and to compare these with the geological data and better refine the model.</p>

Comparing onshore and offshore volumes of large igneous provinces associated with passive continental margins

<p>The rifting of continents can lead to the initiation of seafloor spreading and the formation of passive margins. Magma-rich passive margins, which are defined as being underlain by enormous volumes of igneous rocks, are often associated with large igneous provinces (LIPs). However, the relationship between the igneous units found along these magma-rich passive margins, rifting processes, and LIPs is poorly understood.</p><p>We have developed the VOLMIR (VOLcanic passive Margin Igneous Rocks) dataset to investigate these relationships. VOLMIR is based on seismic reflection profiles in which the volumes and geometries of both shallow seaward dipping reflector (SDR) and deeper high velocity lower crustal (HVLC) units can be measured. We find a relatively consistent ratio of SDR to HVLC volumes, with SDR volumes about one third that of HVLC. This consistency suggests that the units are related during formation and may be used to provide insight into how such units form during continental rifting and breakup. Presumably, as magmas rise and erupt to the surface to form SDRs, the remaining high-density residuum or cumulate becomes the HVLC. The volumes of SDR units are moderately positively correlated with distance from the Euler pole, and weakly negatively correlated with distance from the nearest hotspot, suggesting that lithospheric processes play more of a role in late-stage continental rifting and breakup than hotspot/mantle plume processes.</p><p>The Mid- and South Atlantic Ocean breakups, and the resulting offshore volcanic passive margins, are temporally and spatially associated with the Central Atlantic Magmatic Province (CAMP) and Paraná-Etendeka LIP. Using VOLMIR, we estimate the amount of igneous material in the offshore volcanic passive and compare it to estimates for the adjacent on-land LIPs. The results indicate that a significant volume of volcanics exist in the offshore passive margins, increasing the inferred amount of volcanic output associated with the LIPs. Further studies will provide insight into the relationship between offshore passive margins and on-land LIPs, and perhaps provide further understanding on the role of magmatism in rifting processes.</p>

MORGEN – The Mid Ocean Ridge GENerating algorithm

<p>Paleo-digital elevation models (paleoDEM) based on plate tectonic and paleogeographic reconstructions use age grids of ocean floor to determine ocean bathymetry. In recent years, such age grids have also been developed for now-subducted oceans from the far geological past, as far back as the Neoproterozoic, using geology and paleomagnetism-based estimates of ocean opening. In such reconstructions, mid ocean ridges are drawn based on estimated Euler poles and rotations, and conceptual knowledge on the geometry consisting of spreading ridges and transform faults.</p><p>Current procedures to draw mid ocean ridges in plate tectonic reconstructions are laborious, as new ridges are drawn every time the Euler pole location changes. Fortunately this is also a task that can be automated. We have written an algorithm using pyGPlates that takes as input a smooth curve at the approximate position of the reconstructed mid ocean ridge at the moment of its formation, and then calculates spreading and transform segments according to their typical geometries in modern oceans, assuming symmetric spreading. The algorithm allows gradual readjustment of ridge orientations upon Euler pole changes comparable to documented cases in the modern oceans (e.g., in the Weddell Sea). The algorithm also contains modules that can convert the calculated mid ocean ridges with other plate boundaries to boundary topologies – which can be used as input for the recently published TracerTectonics algorithm, produce isochrons which can be converted to age grids, check for subduction of isochrons and subsequently create bathymetry grids. We illustrate the use of the MORGEN algorithm with recently published reconstructions of subducted, as well as future oceans.</p>

Some evidence for a wide fan-shaped extension of the East Antarctic plate at the Mesozoic-Cenozoic transition

<p>The Transantarctic Mountains (TAM) separate the Mesozoic to recent West Antarctic rift system (WARS) from a wide and depressed triangular sector of East Antarctica spanning from 100° E to 160° E in longitude and from the Oates, George V and Adelie coastlines to 85° S in latitude. The sub-ice bedrock of this sector shows a basin and range style topography comprising two major basins of continental proportions -the Wilkes Basin and the Aurora Basin complex- and many smaller basins such as the Adventure, Concordia, Aurora and Vostok trenches. Most of these basins and trenches exhibit a triangular shape with the acutest angle pointing approximatively to a single pole towards the South, giving a fan shaped pattern of significant dimensions. We name here this region as the East Antarctic Fan shaped Basin Province (EAFBP). To the West, this province is limited by the intraplate Gamburtsev Mountains (GM).</p><p>Origins and inter-relationships between these four fundamental Antarctic tectonic units (WARS, TAM, EAFBP, GM) are still poorly understood and strongly debated. In the EAFBP, very little is known about the mechanism generating the basins, their formation time, whether they are all coeval and if and how they relate to Australia basins before Antarctica-Australia rifting. Present genetic hypotheses for some of the basins span from continental rifting to a purely flexural origin or a combination of the two. Also, post-tectonic erosional and depositional processes may have had a significant impact on the present-day topographic configuration.</p><p>Here we investigate the possibility that the EAFBP is the result of a single genetic mechanism: a wide fan-shaped intra-continental extension around a pivot point at about 135° E, 85° S that occurred at the Mesozoic-Cenozoic transition. We discuss evidence from the sub-ice topography and potential field airborne and satellite data.</p><p>We have used international community-based Antarctic compilations in public domain, including BedMachine (Morlighem et al., 2020), AntGG (Scheinert et al., 2016) and ADMAP 2.0 (Golynsky et al., 2018). We have applied image segmentation techniques to the rebounded sub-ice topography to automatically trace the first order shape of the sub-ice basins. Then we have fitted the edges of the basins by maximum circles and we have estimated the best Euler pole identified by their intersection. Potential field anomalies have been taken into account in order to enlighten major discontinuities not revealed by the sub-ice topography.</p><p>Software simulations of the EAFBP opening in the frame of global plate tectonics reconstructions indicate that it may be inserted in the frame of the later phase of the Antarctica-Australia rifting, giving constraints on timing that allow us to date the EAFBP opening at the Mesozoic-Cenozoic transition.</p><p>The reconnaissance of the EAFBP as the result of a continental-scale fan-shaped extension may have deep implications on global and regional tectonics plate reconstructions, plate deformation assumptions and new tectonic evolutionary models of WARS, TAM and GM.</p>

GNSS-A seafloor geodetic observation capability in 2021 and its applicability to global geodesy

<p>Recently, the GNSS-A (Global Navigation Satellite System – Acoustic combination technique) seafloor geodetic observation technology, developed in the Hydrographic and Oceanographic department in Japan Coast Guard (JCG), was upgraded to be able to monitor a subseafloor interplate coupling condition of about 1 cm/year and an interplate shallow slow slip event of about 5 cm (e.g., Yokota et al., 2018, Scientific Data; Ishikawa et al., 2020, Front. Ear. Sci.). By observing such small-scale seafloor crustal movements, GNSS-A technology makes a decisive contribution to subduction seismology and disaster prevention sciences (e.g., Yokota et al., 2016, Nature; Yokota and Ishikawa, 2020, Sci. Adv.). This technology was achieved by connecting high-precision underwater acoustic ranging technology and high-rate GNSS on a vessel at sea surface.</p><p>The GNSS-A, which is carried out all over the world, especially in the Pacific Rim, has been constructed for observation of plate boundary subduction processes and fault movement processes. Unlike the GNSS network, GNSS-A has never contributed to global geodesy within the framework of the Global Geodetic Observing System (GGOS). However, it can be a unique observation method for the construction of the International Terrestrial Reference Frame (ITRF). It can make an important contribution in determining the movement and Euler pole on an oceanic plate that have few land area.</p><p>In the future, if an extensive seafloor geodetic observation network as shown by Kato et al. (2018, JDR) will be established, there is a possibility of constructing a next-generation reference frame that completely explains the plate motion on the earth's surface. This presentation will present the current state of the GNSS-A ability and cost and future prospects for the contribution to global geodesy.</p>

An analytical model concerning the possible initiation of subduction between the India and Africa plates, caused by the Morondova plume

<p>Identifying the geodynamic processes that trigger the formation of new subduction zones is key to understand what keeps the plate tectonic cycle going, and how plate tectonics once started. Here we discuss the possibility of plume-induced subduction initiation. Previously, our numerical modeling revealed that mantle upwelling and radial push induced by plume rise may trigger plate motion change, and plate divergence as much as 15-20 My prior to LIP eruption. Here we show that, depending on the geometry of plates, the distribution of cratonic keels and where the plume rises, it may also cause a plate rotation around a pole that is located close to the same plate boundary where the plume head impinges: If that occurs near one end of the plate boundary, an Euler pole of the rotation may form along that plate boundary, with extension on one side, and convergence on the other.  This concept is applied to the India-Africa plate boundary and the Morondova plume, which erupted around 90 Ma, but may have influenced plate motions as early as 105-110 Ma. If there is negligible friction, i.e. there is a pre-existing weak plate boundary, we estimate that the total amount of convergence generated in the northern part of the India-Africa plate boundary can exceed 100 km, which is widely thought to be sufficient to initiate forced, self-sustaining subduction. This may especially occur if the India continental craton acts like an “anchor” causing a comparatively southern location of the rotation pole of the India plate. Geology and paleomagnetism-based reconstructions of subduction initiation below ophiolites from Pakistan, through Oman, to the eastern Mediterranean reveal that E-W convergence around 105 Ma caused forced subduction initiation, and we tentatively postulate that this is triggered by Morondova plume head rise. Whether the timing of this convergence is appropriate to match observations on subduction initiation as early as 105 Ma depends on the timing of plume head arrival, which may predate eruption of the earliest volcanics. It also depends on whether a plume head already can exert substantial torque on the plate while it is still rising – for example, if the plate is coupled to the induced mantle flow by a thick craton.</p>