Article Geochemistry 3 Volume 2 Geophysics April XX, 2001 Geosystems Paper number 2000GC000106 G ISSN: 1525-2027 AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society

A three-plate kinematic model for Lau Basin opening

Kirsten E. Zellmer and Brian Taylor Department of Geology and Geophysics, School of Ocean Earth Science and Technology, University of Hawaii, Honolulu, Hawaii 96822 ([email protected]; [email protected])

[1] Abstract: We present revised compilations of bathymetry, magnetization, acoustic imagery, and seismicity data for the Lau Basin and surrounds. We interpret these data to more precisely locate the plate boundaries in the region and to derive a three-plate kinematic model for the opening of the Lau Basin during the Brunhes Chron. Our tectonic model includes the Niuafo'ou microplate, which is separated from the along the Peggy Ridge-Lau Extensional Transform Zone-Central Lau Spreading Center and from the Tongan Plate along the newly discovered Fonualei Rift and Spreading Center. Our model shows that Australia-Tonga plate motion between 15.58S and 198S is partitioned across the microplate, and it resolves the former apparent conflict between geodetic versus velocities. Regionally, this implies that the angular rate of opening of the Lau Basin has been faster than that of the Havre Trough to the south and that the rapid current rates of Pacific at the (240 mm/yr at 168S) have been sustained for at least 0.78 m.y.

Keywords: Lau Basin; plate kinematics; microplate; seafloor spreading; geodetic vectors. Index terms: Marine geology and geophysics; ; plate motions±recent; Pacific Ocean. Received September 1, 2000; Revised January 25, 2001; Accepted March 1, 2001; Published April XX, 2001.

Zellmer, K. E., and B. Taylor, 2001. A three-plate kinematic model for Lau Basin opening, Geochem. Geophys. Geosyst., vol. 2, Paper number 2000GC000106 [8651 words, 13 figures, 7 tables]. Published April XX, 2001.

1. Introduction to be recognized as a young ocean basin formed by splitting a volcanic island arc [2] The Tonga-Fiji region has been key to the [Karig, 1970]. Although subsequent models understanding of Earth dynamics since the of back arc basin evolution were often two- formulation of the plate tectonics paradigm dimensional and thus assumed only two plates [e.g., Isacks et al., 1968]. The 11-km-deep above the subducting lithosphere (see review Tonga Trench is the locus of fastest subduction by Taylor and Karner [1983]; also see Parson [Bevis et al., 1995], and the subducted Pacific and Wright [1996]), Chase [1971] recognized lithosphere generates a majority of all deep the existence of multiple plates and triple earthquakes [Isacks et al.,1968].TheLau junctions in the opening of the (Figure 1) was the first back arc basin Basin. After two decades of controversy

Copyright 2001 by the American Geophysical Union Geochemistry Geophysics 3 zellmer and taylor: three-plate kinematic model 2000GC000106 Geosystems G

175ûE 180û 175ûW 170ûW

Samoa

15ûS 15ûS North Fiji Fiji Basin Lau Basin 20ûS 20ûS

Tonga Trench Lau Ridge Tonga Ridge

25ûS 25ûS

Australian Plate

30ûS 30ûS

Havre Trough 35ûS 35ûS

Pacific Plate

Rift Taupo

40ûS 40ûS 175ûE 180û 175ûW 170ûW km -10 -4 -3 -2 -1 0 Figure 1. Regional bathymetric map showing the location of the Lau Basin and neighboring geologic features. The image is constructed from the predicted bathymetry of Smith and Sandwell [1997]; seafloor relief is shaded from the east, and islands are shown in black. Geochemistry Geophysics 3 zellmer and taylor: three-plate kinematic model 2000GC000106 Geosystems G

[Lawver et al., 1976; Lawver and Hawkins, were collected using a Hydrosweep system on 1978; Weissel, 1977, 1981; Hamburger and MMAJ9507, the HAWAII MR1 (HMR1) side Isacks, 1988; Parson and Hawkins, 1994], scan system on MW9603, and a SeaBeam 2000 Taylor et al. [1996] used a magnetization system on BMRG08MV. Surveys were navi- inversion of bathymetry and magnetic data to gated using the Global Positioning System definitively show that the opening of the Lau (GPS). Detailed descriptions of these three Basin was by seafloor spreading processes surveys and their results will be published fundamentally similar to those occurring at elsewhere (for BMRG08MV, see Wright et al. mid-ocean ridges. Nevertheless, the tectonics [2000]). of the Lau Basin region north of 188S remained poorly understood, and Brunhes Chron spread- ing rates were reported to be only 75% of the 2.1. Bathymetry opening rates measured geodetically [Bevis et [5] Using GMT [Wessel and Smith, 1991] and al., 1995]. Revised geodetic results [Bevis, MB-system software [Caress and Chayes, 1997] partially closed the gap between the 1995], we created a 0.0028 (200 m) grid two data sets, but a quantitative kinematic of bathymetry (Figure 2). The compilation model of Lau Basin opening has not been incorporates the MMAJ9507, MW9603, and attempted previously. BMRG08MV swath bathymetry with (1) the predicted bathymetry of Smith and Sandwell [3] Weuseasetofrecompiledbathymetry, [1997], (2) data obtained from the U. S. magnetization, acoustic imagery, and seismicity National Geophysical Data Center (NGDC), data to define a multiple plate kinematic model (3) data obtained from the South Pacific that (1) explains a majority of the tectonic Applied Geoscience Commission (SOPAC), features throughout the Lau Basin, (2) is con- (4) digitized contours of published bathymetry sistent with both spreading rates determined maps [Wiedicke and Collier, 1993; Wiedicke from the marine magnetic data and the revised and Habler, 1993; Defense Mapping Agency, geodetic velocities, and (3) satisfies rigid plate 1990a, 1990b], and (5) additional unpublished criteria. We address the recent tectonics (0.78 swath bathymetry data (R/V Ewing 9512 and Ma to present) and implications of a three-plate 9914; R/V Yokosuka 9614, refer to Fujiwara system for Lau Basin opening. et al. [2001]). The bathymetry compilation process was completed in stages similar to 2. Data Compilation those described by Martinez et al. [1995] in order to obtain a grid that best merges all the [4] We include new results from three system- data of varying quality and resolution. To atic surveys to update the Lau Basin data achieve this end, all bathymetric data sets compilation and magnetization inversion of were first combined and median filtered in 1 Taylor et al. [1996]. R/V Moana Wave cruise arc min (1.85 km) blocks to generate a MW9603 surveyed the basin from edge to edge regional base data set that smoothly merges between 178S and 188350S. A R/V Hakerei with the more detailed data. This regional data Maru 2 cruise (hereafter referred to as set was then combined with the original data, MMAJ9507) surveyed the spreading axis and median filtered in 0.0028 blocks, and gridded nearby ridge flanks between 198250S and 228S. at 0.0028 node spacing. High-resolution swath R/V Melville cruise BMRG08MV principally data that were considered to be the most surveyed the forearc region [Wright et al., reliable and best navigated were given greatest 2000]. Swath bathymetry and acoustic imagery priority within the grid. Geochemistry Geophysics 3 zellmer and taylor: three-plate kinematic model 2000GC000106 Geosystems G

180û 179ûW 178ûW 177ûW 176ûW 175ûW 174ûW 173ûW 172ûW 14ûS 14ûS

15ûS 15ûS

16ûS 16ûS

17ûS 17ûS

18ûS 18ûS

19ûS 19ûS

20ûS 20ûS

21ûS 21ûS

22ûS 22ûS

23ûS 23ûS 180û 179ûW 178ûW 177ûW 176ûW 175ûW 174ûW 173ûW 172ûW

km

-10 -4 -3 -2 -1 0 Figure 2. Bathymetry of the Lau Basin region, with