Sabin Zahirovic

Insights on the kinematics of the India‐Eurasia collision from global geodynamic models

The Eocene India-Eurasia collision is a first order tectonic event whose nature and chronology remains controversial. We test two end-member collision scenarios using coupled global plate motion-subduction models. The first, conventional model, invokes a continental collision soon after ∼60 Ma between a maximum extent Greater India and an Andean-style Eurasian margin. The alternative scenario involves a collision between a minimum extent Greater India and a NeoTethyan back-arc at ∼60 Ma that is subsequently subducted along southern Lhasa at an Andean-style margin, culminating with continent-continent contact at ∼40 Ma. Our numerical models suggest the conventional scenario does not adequately reproduce mantle structure related to Tethyan convergence. The alternative scenario better reproduces the discrete slab volumes and their lateral and vertical distribution in the mantle, and is also supported by the distribution of ophiolites indicative of Tethyan intraoceanic subduction, magmatic gaps along southern Lhasa and a two-stage slowdown of India. Our models show a strong component of southward mantle return flow for the Tethyan region, suggesting that the common assumption of near-vertical slab sinking is an oversimplification with significant consequences for interpretations of seismic tomography in the context of subduction reference frames.

Mantle-induced subsidence and compression in SE Asia since the early Miocene

Rift basins developed extensively across Sundaland, the continental core of Southeast Asia, since the Eocene. Beginning in the early Miocene, basins in southern Sundaland experienced widespread synchronous compression (inversion) and marine inundation, despite a large drop in long-term global sea level. The mechanism for this large-scale synchronous regional sea level rise, basin inversion, and subsidence is not well understood and contrary to expectations from traditional basin models and eustatic sea level trends. We present geodynamic models of mantle convection with both deformable and rigid plate reconstructions to investigate this enigma. Models suggest that a slab stagnates within the transition zone beneath Southeast Asia before the Miocene. The stagnant slab penetrated through the 660 km mantle discontinuity during the early Miocene and formed a slab avalanche event, due to continuous subduction and accumulation of negatively buoyant slabs. This avalanche may have induced large-scale marine inundation, regional compression, and basin inversion across southern Sundaland. We argue mantle convection induced large-scale basin compression, in contrast to conventional plate margin-induced compression; this suggests mantle convection may exert a much stronger control on surface processes than previously recognized.

Obituary: Sandy Island (1876–2012)

In October 2012, scientists investigating the tectonic evolution of the eastern Coral Sea aboard the R/V Southern Surveyor uncovered a quirky discrepancy in maps of seafloor topography during their 25-day voyage. While on a transit leg between dredge sites, the ship passed near a purported island between the Chesterfield Islands and Nereus Reef that appeared in numerous scientific data sets and in Google Earth™ with the label “Sandy Island.” However, this 25-kilometer-long and 5-kilometer-wide feature was absent from the hydrographic charts used by the master onboard the ship for navigation.

GPlates: Building a Virtual Earth Through Deep Time

GPlates is an open‐source, cross‐platform plate tectonic geographic information system, enabling the interactive manipulation of plate‐tectonic reconstructions and the visualization of geodata through geological time. GPlates allows the building of topological plate models representing the mosaic of evolving plate boundary networks through time, useful for computing plate velocity fields as surface boundary conditions for mantle convection models and for investigating physical and chemical exchanges of material between the surface and the deep Earth along tectonic plate boundaries. The ability of GPlates to visualize subsurface 3‐D scalar fields together with traditional geological surface data enables researchers to analyze their relationships through geological time in a common plate tectonic reference frame. To achieve this, a hierarchical cube map framework is used for rendering reconstructed surface raster data to support the rendering of subsurface 3‐D scalar fields using graphics‐hardware‐accelerated ray‐tracing techniques. GPlates enables the construction of plate deformation zones—regions combining extension, compression, and shearing that accommodate the relative motion between rigid blocks. Users can explore how strain rates, stretching/shortening factors, and crustal thickness evolve through space and time and interactively update the kinematics associated with deformation. Where data sets described by geometries (points, lines, or polygons) fall within deformation regions, the deformation can be applied to these geometries. Together, these tools allow users to build virtual Earth models that quantitatively describe continental assembly, fragmentation and dispersal and are interoperable with many other mapping and modeling tools, enabling applications in tectonics, geodynamics, basin evolution, orogenesis, deep Earth resource exploration, paleobiology, paleoceanography, and paleoclimate.