On the Role of Slab Pull in the Cenozoic Motion of the Pacific Plate Claudio Faccenna, Thorsten W
Total Page:16
File Type:pdf, Size:1020Kb
On the role of slab pull in the Cenozoic motion of the Pacific plate Claudio Faccenna, Thorsten W. Becker, Serge Lallemand, Bernhard Steinberger To cite this version: Claudio Faccenna, Thorsten W. Becker, Serge Lallemand, Bernhard Steinberger. On the role of slab pull in the Cenozoic motion of the Pacific plate. Geophysical Research Letters, American Geophysical Union, 2012, 39, pp.L03305. 10.1029/2011GL050155. hal-00682150 HAL Id: hal-00682150 https://hal.archives-ouvertes.fr/hal-00682150 Submitted on 25 Jan 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. GEOPHYSICAL RESEARCH LETTERS, VOL. 39, L03305, doi:10.1029/2011GL050155, 2012 On the role of slab pull in the Cenozoic motion of the Pacific plate Claudio Faccenna,1,2 Thorsten W. Becker,3 Serge Lallemand,4 and Bernhard Steinberger5,6 Received 27 October 2011; revised 29 December 2011; accepted 3 January 2012; published 8 February 2012. [1] We analyze the role of slab pull acting on the Pacific plate using plate reconstructions [e.g., Lithgow-Bertelloni and plate during its early Tertiary change in motion. Slab pull Richards, 1998; Conrad and Lithgow-Bertelloni, 2002, 2004]. forces are estimated by integrating the negative buoyancy [3] The Pacific plate represents a key site to analyze of a 700 km long slab along a revised subduction boundary the different force contributions as it is the fastest plate on model adopting the Müller et al. (2008) seafloor age Earth and is half bounded by subduction zones. In addition, reconstructions. Our results indicate that torques predicted hot spot tracks record a sharp bend that is likely related to from a simple slab pull model match the Pacific plate force re-organization. Along the Hawaiian-Emperor chain, Euler vectors during the Tertiary fairly well. The change this bend in the track from NNW to WNW direction (the of the Pacific motion at 50–40 Ma appears to be driven Hawaiian-Emperor Bend, HEB) occurred between 50 and by the onset of the Izu-Bonin-Mariana system and, soon 43 Ma [Steinberger et al., 2004; Sharp and Clague, 2006; afterwards, by the Tonga-Kermadec subduction zones. Whittaker et al., 2007]. Different hypotheses have been Citation: Faccenna, C., T. W. Becker, S. Lallemand, and proposed to explain the Pacific plate reorganization and in B. Steinberger (2012), On the role of slab pull in the Cenozoic turn related to different processes occurring in or around motion of the Pacific plate, Geophys. Res. Lett., 39, L03305, the Pacific during Cenozoic times [Lithgow-Bertelloni and doi:10.1029/2011GL050155. Richards, 1998; Conrad and Lithgow-Bertelloni, 2004]. [4] Here, we revise the plate geometry and the tectonic 1. Introduction event history in/around the Pacific plate, drawing exten- sively on the Müller et al. [2008] plate and seafloor age [2] The motion of tectonic plates is generally considered reconstructions and the compilation by Gurnis et al. [2012]. to be driven by mantle convection, and the plates themselves We compute simplified estimates of the expected slab pull can be viewed as the top thermo-chemical boundary layer force acting on the Pacific plate and show that the directions of the mantle. However, it is nonetheless useful to evaluate of Pacific plate motion are directly related to the slab pull the different force contributions separately [Forsyth and acting on its boundary. Those forces increase rapidly during Uyeda, 1975] because this facilitates regional analysis and the inception of the Izu-Bonin-Mariana subduction, moving allows inferences on important parameters such as slab the Pacific to a more westerly direction. strength. Plates can be: (1) pulled directly into the mantle by the negative buoyancy of the attached, subducting oceanic 2. Methods lithosphere (slab pull [e.g., Forsyth and Uyeda, 1975; Chapple and Tullis, 1977]), (2) dragged at their base by [5] Several regional tectonic events took place in and mantle flow excited by density (thermal) anomalies (slab around the Pacific during the early Tertiary: the onset of the suction [e.g., McKenzie, 1969; Turcotte and Oxburgh, 1967; volcanic arc in the Western Pacific along the Izu Bonin Forsyth and Uyeda, 1975; Hager and O’Connell, 1981]), Mariana (IBM) trench at 55 Ma, the accretion of the (3) propelled by potential energy variations over the oceanic Okhotsk and Aleutian terranes at the junction between plate (ridge push [Forsyth and Uyeda, 1975]) or moved by Kamchatka and the Aleutians (65 to 55 Ma), the subduc- a combination of some or all these forces [e.g., Wortel et al., tion of the Kula-Pacific Ridge beneath Kamchatka-Kuril- 1991; Becker and O’Connell, 2001; Conrad and Lithgow- Japan and Ryuku (60 to 50 Ma), and the collision of India Bertelloni, 2004]. Determining the mode of force transmis- with Asia at 55–50 Ma [Patriat and Achache, 1984]. sion and discerning between the different contributions At least some of these events could have profound con- from resisting forces within the mantle is complicated by sequences on the force equilibrium on the plate. uncertainties such as those about lithospheric and mantle [6] Here, we test the idea that the change in the plate- rheology [e.g., Conrad and Lithgow-Bertelloni,2004;Becker, mantle force distribution that led to the plate reorganization 2006]. One possible solution is to analyze the motion of the is reflected in inferred seafloor ages and plate boundaries. The most direct consequence of changing plate configura- tions is represented by the variation of the slab pull torque 1 acting on the Pacific. We compute the slab pull torque by Dipartimento di Scienze Geologiche, Università degli Studi Roma Tre, Rome, Italy. integrating the inferred slab pull along the entire subduction 2IGAG, CNR, Rome, Italy. zone. Following Turcotte and Schubert [1982], slab pull can 3Department of Earth Sciences, University of Southern California, Los be estimated as: Angeles, California, USA. 4 = Laboratoire Géosciences Montpellier, Université Montpellier 2, F ¼ 0:25Drgl Á 2:32ðÞ kA 1 2 ð1Þ CNRS, Montpellier, France. sp 5GeoForschungsZentrum Potsdam, Potsdam, Germany. 6PGP, University of Oslo, Oslo, Norway. where l is slab length, assumed constant at 700 km consid- ering the present mean slab dip of 55° of the western Pacific Copyright 2012 by the American Geophysical Union. slab [Lallemand et al., 2005], g is gravitational acceleration 0094-8276/12/2011GL050155 L03305 1of6 L03305 FACCENNA ET AL.: SLAB PULL IN THE PACIFIC L03305 (≈10 msÀ2), Dr is the density contrast between the slab soon after the subsidence episode, the inactive arc was and the mantle, held constant at 70 kg mÀ3. Slab thickness obducted onto the continental margin (53–50 Ma). Arc is estimated from half-space cooling using the square volcanism along the Kuril arc resumed then at 46 Ma root of age of subducting lithosphere at trench (A), taken [Konstantinovskaya, 2011]. from Müller et al. [2008] for thermal diffusivity (k) equal to [10] The rest of the Kula-Pacific ridge southwest and 10À6 m2 sÀ1. The coefficient 0.25 accounts for the decreas- northeast of the proto-Kuril trench reached the trenches also ing of the excess mass due to conductive heating of the shortly before 60 Ma. In the Müller et al. [2008] recon- slab with depth [Carlson et al., 1983]. We integrate the force struction, the Kula-Pacific ridges system in fact was almost Fsp along the convergent plate boundary segments to deter- parallel to trenches [Whittaker et al., 2007]. Hence, an mine the total torque due to the subduction zone edge forces, almost contemporaneous surge of volcanism is expected all and also express the torque as force per area evaluated at along the Pacific margin. Previous models, conversely, sug- the Pacific plate centroid. Our reconstruction is based on the gested that a ridge subduction episode occurred 10–20 Myrs revised subduction zone geometry and on the geometry of earlier [Maruyama and Seno, 1986; Kinoshita, 2002]. One Gurnis et al. [2012], and computed slab pull is compared to other possible consequence of the arrival of the ridge into the velocity reconstruction from Müller et al. [2008, hereinafter subduction zone, northeast of the paleo-Kuriles at 55 Ma, M08]. We considered the direction of subduction and is the accretion of the Aleutia arc terrane. This terrane detailed shape of the plate boundary based on a new updated belonged to the former Kula plate and was accreted at reconstruction and computed forces as discussed below. 50 Ma [Scholl, 2007; Minyuk and Stone, 2009]. A new [7] For comparison, we also compute slab-pull estimated subduction zone thus developed possibly along a Kula-Pacific using subduction zones as compiled by the model of Gurnis ridge segment, becoming the proto-Aleutian arc, where the et al. [2012], with ages of subducted plates as compiled by oldest volcanic rocks were dated at 46 Ma [Scholl,2007]. Steinberger and Torsvik [2010] (see auxiliary material) and [11] At the same time 55–50 Ma, several thousands compare to absolute plate motions from Steinberger and of kilometers to the south near the paleo-equator, the proto- Gaina [2007, hereinafter SG07].1 The Müller et al. [2008] Izu-Bonin-Mariana (IBM) trench was initiating, probably model uses a global plate circuit to relate Pacific to African along a former transform boundary [e.g., Uyeda and Ben- plate motion and, as a consequence, the direction of Pacific Avraham, 1972; Stern and Bloomer, 1992].