Rapid Microplate Rotations and Backarc Rifting at the Transition Between Collision and Subduction
Total Page:16
File Type:pdf, Size:1020Kb
Rapid microplate rotations and backarc rifting at the transition between collision and subduction Laura M. Wallace Institute of Geological and Nuclear Science, Lower Hutt 6315, New Zealand Robert McCaffrey Department of Earth and Environmental Science, Rensselaer Polytechnic Institute, 110 8th Street, Troy, New York 12180, USA John Beavan Institute of Geological and Nuclear Science, Lower Hutt 6315, New Zealand Susan Ellis ABSTRACT In all regions studied, we de®ne the geometry Using global positioning system velocities from convergent plate boundaries in Papua of block boundaries (faults) based on earth- New Guinea, New Zealand, Tonga, Vanuatu, and the Marianas, we note a spatial corre- quakes and geology. To place the regional lation between rapid tectonic block rotations and the transition from subduction to col- GPS velocities into larger plate kinematic con- lision. We present a mechanism for the block rotations, in which the change from collision texts, we use GPS velocities from sites on the of a buoyant indentor to normal subduction exerts a torque on the upper-plate microplate. bounding plates (e.g., Australian, Paci®c, and This work improves our understanding of the causes of rapid vertical axis rotations, often Philippine Sea plates) given in Beavan et al. observed in paleomagnetic studies. We also show how collision-induced rotations may lead (2002) and Sella et al. (2002). to backarc rifting. Papua New Guinea Keywords: GPS, global positioning system, collision, subduction, microplate rotation, backarc Papua New Guinea forms the boundary rifting, tectonics. zone between the obliquely and rapidly (;110 mm/yr) converging Australian and Paci®c INTRODUCTION the interseismic period are hampered by the plates (Fig. 1). The southern boundary of the The fastest modern-day tectonic block ro- effects of elastic strain due to locking on faults South Bismarck microplate changes (from tations, as much as 98/m.y., occur in the fore- (particularly the subduction zone). To deal west to east) from collision between the Fin- arcs of convergent plate margins (e.g., Cal- with this problem, we simultaneously invert isterre Range and the underthrusting New mant et al., 2003; Wallace et al., 2004a). It GPS velocities, earthquake slip vector azi- Guinea Highlands to subduction of oceanic has been suggested that the bending and ro- muths, transform fault orientations, and long- crust at the New Britain trench in the Solomon tation of subduction-zone forearcs can be term geological fault slip rates (including sea- Sea. The northern boundary of the South Bis- caused by collision of buoyant features with ¯oor spreading rates) for angular velocities of marck microplate is marked by the Manus and the subduction zone (e.g., Vogt et al., 1976; tectonic blocks and the degree of coupling on New Guinea backarc rifts, which are connect- McCabe, 1984). Earlier studies of forearc faults bounding the blocks (e.g., McCaffrey, ed by left-lateral transform faults along the block rotations were based on paleomagnetic 1995, 2002). Throughout this paper we refer Bismarck Sea seismic lineation (e.g., Taylor, declination anomalies that can provide the to the point where the angular velocity inter- 1979). GPS velocities (Wallace et al., 2004a) long-term vertical axis rotation rate, but can- sects Earth's surface as the ``pole of rotation.'' show that the South Bismarck microplate ro- not constrain the full angular velocity, so that We refer to the tectonic blocks or microplates tates clockwise relative to the New Guinea the orientation of the spin axis is unknown. (often coinciding with the forearc and/or arc) Highlands microplate (the lower plate of the Using published global positioning system as ``convergent plate boundary microblocks.'' collision) at ;98/m.y. about a pole where the (GPS) velocities from ®ve convergent margins around the world, we estimate the angular ve- locities of forearc blocks to reveal an apparent Figure 1. Tectonic setting spatial relationship between the rotation of up- of Papua New Guinea and per-plate microplates and the collisions of vectors showing rotation buoyant crustal masses (e.g., oceanic plateaus, of South Bismarck plate seamount chains, or continental fragments) (SBS) relative to lower plate (see Wallace et al., with the subduction zones. We suggest that the 2004a, for details of mod- change from collision to subduction exerts a eling results). Pole of ro- torque on the upper-plate microplate, causing tation for South Bismarck the microplate to rotate rapidly (relative to the microplate relative to New Guinea Highlands subducting plate) about a point near where the (NGH) is shown by semi- collision occurs. These rotations often result circular arrow, clockwise. in backarc rifting. Our observations have im- Rates of rifting in Manus plications for the knowledge of what drives Basin (white arrows) are rapid crustal block rotations, and the causes of given in mm/yr. BSSLÐ Bismarck Sea seismic lin- backarc rifting. eation. In best-®tting x2 5 model, we obtain n KINEMATIC MODELING OF BLOCK 1.99, using 333 earth- ROTATIONS AT CONVERGENT quake slip vectors and 66 PLATE BOUNDARIES global positioning sys- tem velocities to estimate Estimates of tectonic block rotations at sub- 44 parameters. duction margins from GPS velocities during q 2005 Geological Society of America. For permission to copy, contact Copyright Permissions, GSA, or [email protected]. Geology; November 2005; v. 33; no. 11; p. 857±860; doi: 10.1130/G21834.1; 5 ®gures. 857 Figure 3. Tectonic setting of Tonga (Tong) and Va- nuatu (Van), global posi- tioning system (GPS) ve- locities (observed, black; model, red), and rotation poles of Vanuatu and Tonga tectonic blocks rel- ative to subducting plate. GPS velocity vectors shown relative to Austra- lian (Aus) plate. White vectors with rates next to them show relative mo- tion in backarc (rates are in mm/yr, predicted from block model). For best x2 5 model, we obtain n 2.13, using 952 earth- quake slip vectors, 4 spreading rates, and 83 GPS velocities to esti- mate 42 free parameters. NHTÐNew Hebrides Trench; EFTÐErromango-Futuna Trough; WTPÐWest Torres Plateau; DERÐD'Entrecasteaux Ridge; VANÐVanuatu; WNFBÐwestern North Fiji Basin; HHRÐHazel Holme Ridge; CSRÐcentral spreading ridge. velocities show that several crustal blocks in wise rotation rate of Tonga estimated from pa- the eastern North Island rotate clockwise at leomagnetic studies (Sager et al., 1994). 3.88±5.08/m.y. relative to the subducting Pa- ci®c plate (Wallace et al., 2004b). Paleomag- Vanuatu Figure 2. Rotation of New Zealand tectonic netic declinations and geological strain esti- Vanuatu straddles the Australia±Paci®c blocks relative to Paci®c plate and tectonic mates also suggest clockwise rotation of the plate boundary zone between the Tonga-Fiji setting (see Wallace et al., 2004b, for details eastern North Island (e.g., Walcott, 1984; region and New Guinea (Fig. 3). The major of modeling results). Tectonic blocks and Beanland and Haines, 1998). Backarc rifting tectonic features in this region are the New their respective poles of rotation are labeled (block names: WAIRÐWairarapa, RAUKÐ in the Taupo volcanic zone (central North Is- Hebrides Trench, where Australia subducts be- Raukumara, CHKGÐCentral Hikurangi, land) is related to eastern North Island rota- neath Vanuatu, and a complex series of rifts AXIRÐAxial Ranges). White vectors are la- tion. The poles of rotation of the eastern North and transforms in the North Fiji Basin. Sub- beled with extension rates in Taupo volcanic Island tectonic blocks relative to the subduct- duction of various buoyant features at the New zone (TVZ) and convergence rates at Hiku- ing Paci®c plate are located adjacent to where Hebrides Trench (e.g., D'Entrecasteaux Ridge; rangi Trough (rates in mm/yr). Note: arrows in TVZ are not to scale. In our best-®tting the buoyant Chatham Rise impinges upon the West Torres Plateau) has in¯uenced the recent x2 5 model, we obtain n 1.32 from 367 GPS margin (Fig. 2). deformation of the Vanuatu island chain (e.g., velocities, 10 earthquake slip vectors, and Pelletier et al., 1998; Calmant et al., 2003). 138 free parameters. Tonga Schellart et al. (2002) suggested that the North The Paci®c plate subducts westward be- Fiji Basin is produced by asymmetric opening Finisterre collision is most advanced (Fig. 1). neath the Australian plate (Fig. 3) at the Tonga about a hinge point located near 118S, 1658E, South Bismarck microplate rotation results in Trench east of the Tonga Islands. West of Ton- and that this asymmetric opening is produced rapid rifting in the Manus and New Guinea ga, some of the highest documented rates of by variation in the degree of rollback of the Basins. Paleomagnetic evidence indicates that backarc rifting in the world are in the Lau Ba- subducting Australian plate. the South Bismarck microplate has been ro- sin (e.g., Bevis et al., 1995; Taylor et al., Following Calmant et al. (2003), we divide tating clockwise at an average rate of 88/m.y. 1996). The southern limit of rifting in the Lau Vanuatu into central and southern tectonic since the Finisterre collision initiated (Weiler Basin is near 268S, adjacent to where the Lou- blocks and include the western North Fiji Ba- and Coe, 2000), consistent with the geodetic isville Seamount chain subducts, and rifting sin as a separate block (Fig. 3). The eastern rates. propagates south at ;100 mm/yr (Parson and boundaries of the southern and central blocks Wright, 1996; Taylor et al., 1996) (Fig. 3). coincide with the zone of backarc thrusting New Zealand Much of the rifting history of the eastern Lau and extension just to the east of the island The plate boundary in the North Island, Basin is linked to the 128 mm/yr southward chain. We invert GPS velocities (Calmant et New Zealand, accommodates subduction of migration of the Louisville Seamount chain± al., 2003), spreading rates (Pelletier et al., the Paci®c plate beneath the North Island at Tonga Trench intersection (Pelletier et al., 1998), and earthquake slip vectors for angular the Hikurangi Trough (Fig.