Seismic tomography of Shatsky Rise by adaptive importance sampling

[1] We present results from a wide-angle seismic refraction survey over Shatsky Rise, a large oceanic plateau in the northwestern Pacific. A new Monte Carlo sampling scheme is developed to explore comprehensively the model space of crustal velocity structure by joint refraction and reflection seismic tomography. The new scheme, which is founded on the notion of adaptive importance sampling, is made possible by the combination of several independent developments, including the introduction of effective model parameters, the implementation of automated model regularization, and a polynomial expansion of the probability density function. From 21 ocean-bottom instruments deployed along a 420-km-long refraction transect, we collect 16446 Pg and 4735 PmP travel times, which are inverted to construct a two-dimensional compressional velocity model across the major axis of Shatsky Rise. The convergence of Monte Carlo sampling is tested by running parallel sampling chains, and the effective dimensions of the model space are estimated to be ∼10 by principal component analysis. The thickness of the rise crust varies from ∼9 km to ∼30 km, and the average velocity of the lower crust ranges from ∼7.0 km s−1 to ∼7.2 km s−1. One standard deviation uncertainty for whole-crustal thickness and average lower-crustal velocity is <1 km and ∼0.05 km s−1, respectively, for most of the model domain. Shatsky Rise crustal structure exhibits a negative correlation between crustal thickness and velocity, which implies that a chemically anomalous mantle may have been responsible for the formation of the rise.

[1]  W. Sliter,et al.  1. SHATSKY RISE: SEISMIC STRATIGRAPHY AND SEDIMENTARY RECORD OF PACIFIC PALEOCEANOGRAPHY SINCE THE EARLY CRETACEOUS1 , 1993 .

[2]  Douglas R. Toomey,et al.  Tomographic inversion of local earthquake data from the Hengill‐Grensdalur Central Volcano Complex, Iceland , 1989 .

[3]  R. White,et al.  Crustal structure of the Hatton and the conjugate east Greenland rifted volcanic continental margins, NE Atlantic , 2009 .

[4]  Jinho Kim,et al.  Bathymetry of Shatsky Rise, northwest Pacific Ocean: Implications for ocean plateau development at a triple junction , 1999 .

[5]  T. Vincenty.  DIRECT AND INVERSE SOLUTIONS OF GEODESICS ON THE ELLIPSOID WITH APPLICATION OF NESTED EQUATIONS , 1975 .

[6]  Kazuo Kobayashi,et al.  Mesozoic magnetic anomaly lineations and seafloor spreading history of the northwestern Pacific , 1989 .

[7]  J. Tarduno,et al.  Rapid Formation of Ontong Java Plateau by Aptian Mantle Plume Volcanism , 1991, Science.

[8]  M. Toksöz,et al.  Nonlinear refraction and reflection travel time tomography , 1998 .

[9]  H. Samuel,et al.  Beyond the thermal plume paradigm , 2005 .

[10]  W. Ludwig,et al.  Seismic refraction measurements on the Ontong Java Plateau northeast of New Ireland , 1973 .

[11]  P. Kelemen,et al.  Methods for resolving the origin of large igneous provinces from crustal seismology , 2002 .

[12]  W. Wilcock,et al.  Segment-scale variations in the crustal structure of 150–300 kyr old fast spreading oceanic crust (East Pacific Rise, 8°15′N–10°5′N) from wide-angle seismic refraction profiles , 2003 .

[13]  L. Kroenke,et al.  Crustal structure of the Shatsky Rise from seismic refraction measurements , 1980 .

[14]  Paul D. Asimow,et al.  Temperatures in ambient mantle and plumes: Constraints from basalts, picrites, and komatiites , 2007 .

[15]  Jun Korenaga,et al.  Thermal cracking and the deep hydration of oceanic lithosphere: A key to the generation of plate tectonics? , 2007 .

[16]  M. Tolstoy,et al.  Crustal thickness variations along the Southeast Indian Ridge (100°–116°E) from 2‐D body wave tomography , 2008 .

[17]  B. Mason.  Composition of the Earth , 1966, Nature.

[18]  H. C. Larsen,et al.  Crustal structure of the southeast Greenland margin from joint refraction and reflection seismic tomography , 2000 .

[19]  Eiichi Takahahshi,et al.  Origin of the Columbia River basalts: melting model of a heterogeneous plume head , 1998 .

[20]  P. Charvis,et al.  Seismic structure of Cocos and Malpelo Volcanic Ridges and implications for hot spot‐ridge interaction , 2003 .

[21]  J. Korenaga.  Velocity–depth ambiguity and the seismic structure of large igneous provinces: a case study from the Ontong Java Plateau , 2011 .

[22]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[23]  Shu-Chuan Lin,et al.  Dynamics of thermochemical plumes: 2. Complexity of plume structures and its implications for mapping mantle plumes , 2006 .

[24]  J. Korenaga.  Mantle mixing and continental breakup magmatism , 2004 .

[25]  A. Sobolev,et al.  The Amount of Recycled Crust in Sources of Mantle-Derived Melts , 2007, Science.

[26]  R. White,et al.  Oceanic crustal thickness from seismic measurements and rare earth element inversions , 1992 .

[27]  Neil Gershenfeld,et al.  The nature of mathematical modeling , 1998 .

[28]  D. Hussong,et al.  The crustal structure of the Ontong Java and Manihiki oceanic plateaus , 1979 .

[29]  P. Kelemen,et al.  Origin of thick, high‐velocity igneous crust along the U.S. East Coast Margin , 1995 .

[30]  Don L. Anderson,et al.  The thermal state of the upper mantle; No role for mantle plumes , 2000 .

[31]  J. Korenaga,et al.  Chemical composition of Earth's primitive mantle and its variance: 1. Method and results , 2007 .

[32]  M. Nafi Toksöz,et al.  Nonlinear refraction traveltime tomography , 1998 .

[33]  J. Korenaga.  Why did not the Ontong Java Plateau form subaerially , 2005 .

[34]  M. Nafi Toksöz,et al.  Nonlinear refraction traveltime tomography , 1996 .

[35]  P. Kelemen,et al.  Major element heterogeneity in the mantle source of the North Atlantic igneous province , 2000 .

[36]  R. White,et al.  Magmatism at rift zones: The generation of volcanic continental margins and flood basalts , 1989 .

[37]  A. Klaus,et al.  Magnetic lineations within Shatsky Rise, northwest Pacific Ocean: Implications for hot spot‐triple junction interaction and oceanic plateau formation , 1999 .

[38]  M. Odegard,et al.  Seismic studies on the Ontong Java Plateau, 1970☆ , 1976 .

[39]  R. Detrick,et al.  Crustal evolution over the last 2 m.y. at the Mid‐Atlantic Ridge OH‐1 segment, 35°N , 2001 .

[40]  S. McNutt,et al.  Isopach Map of Sediments in the Pacific Ocean Basin and Marginal Sea Basins , 1980 .

[41]  S. Operto,et al.  Deep structure of the southern Kerguelen Plateau (southern Indian Ocean) from ocean bottom seismometer wide‐angle seismic data , 1996 .

[42]  R. White,et al.  Influence of the Iceland mantle plume on oceanic crust generation in the North Atlantic , 2008 .

[43]  M. Richards,et al.  Flood Basalts and Hot-Spot Tracks: Plume Heads and Tails , 1989, Science.

[44]  H. Shimamura,et al.  Asymmetric plume‐ridge interaction around Iceland: The Kolbeinsey Ridge Iceland Seismic Experiment , 2006 .

[45]  J. Mahoney,et al.  Geochemistry and Age of the Ontong Java Plateau , 2013 .

[46]  Paul J. Wallace,et al.  Proceedings of the Ocean Drilling Program, 192 Initial Reports , 2001 .

[47]  E. Flueh,et al.  Crustal intrusion beneath the Louisville hotspot track , 2010 .

[48]  M. Ewing,et al.  Lower Cretaceous Sediments From The Northwest Pacific , 1966, Science.

[49]  W. Sager,et al.  How Do Oceanic Plateaus Form? Clues From Drilling at Shatsky Rise , 2011 .

[50]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[51]  H. Kopp,et al.  Structural architecture of oceanic plateau subduction offshore Eastern Java and the potential implications for geohazards , 2011 .

[52]  Paul J. Wallace,et al.  Proceedings of the Ocean Drilling Program, 183 Initial Reports , 2000 .

[53]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[54]  Robert B. Smith,et al.  Seismic traveltime inversion for 2-D crustal velocity structure , 1992 .

[55]  Barbara Romanowicz,et al.  GLOBAL MANTLE TOMOGRAPHY: Progress Status in the Past 10 Years , 2003 .

[56]  A. Taira,et al.  Seismological structure and implications of collision between the Ontong Java Plateau and Solomon Island Arc from ocean bottom seismometer-airgun data , 2004 .

[57]  P. Charvis,et al.  Seismic structure of the Carnegie ridge and the nature of the Galápagos hotspot , 2005 .

[58]  J. Ewing,et al.  Seismic‐refraction measurements in the Northwest Pacific Basin , 1969 .

[59]  T. Hilde,et al.  Tectonic evolution of the northern Pacific plate and Pacific-Farallon Izanagi triple junction in the Late Jurassic and Early Cretaceous (M21-M10) , 1988 .

[60]  Ian H. Campbe.  Implications of mantle plume structure for the evolution of flood basalts , 2002 .