Research Article

The Emperor Seamounts: Southward Motion of the Hawaiian Hotspot Plume in Earth's Mantle

See allHide authors and affiliations

Science  22 Aug 2003:
Vol. 301, Issue 5636, pp. 1064-1069
DOI: 10.1126/science.1086442

Abstract

The Hawaiian-Emperor hotspot track has a prominent bend, which has served as the basis for the theory that the Hawaiian hotspot, fixed in the deep mantle, traced a change in plate motion. However, paleomagnetic and radiometric age data from samples recovered by ocean drilling define an age-progressive paleolatitude history, indicating that the Emperor Seamount trend was principally formed by the rapid motion (over 40 millimeters per year) of the Hawaiian hotspot plume during Late Cretaceous to early-Tertiary times (81 to 47 million years ago). Evidence for motion of the Hawaiian plume affects models of mantle convection and plate tectonics, changing our understanding of terrestrial dynamics.

The concept of an age-progressive set of volcanic islands, atolls, and seamounts produced by a hotspot plume fixed in the deep mantle was first developed to explain the Hawaiian Islands (1). The bend separating the westward-trending Hawaiian island chain from the northward-trending Emper- or Seamounts has most often been interpreted as an example of a change in plate motion recorded in a fixed-hotspot frame of reference (2).

However, global plate circuits suggest large relative motions between Hawaii and hotspots in the Atlantic and Indian Oceans (36). Improved mapping of marine magnetic anomalies in the Pacific has failed to define the directional change at 43 million years ago (Ma) (7) that would be expected if such a large change in plate motion had occurred. There was also a general lack of circum-Pacific tectonic events (8) documented for this time. Recent age data suggest a slightly older age for the bend [∼47 Ma (9)], but this revised timing still does not correspond to an episode of profound plate motion change recorded within the Pacific basin or on its margins.

One approach to examine hotspot fixity is to determine the age and paleolatitude of volcanoes that form a given hotspot track. For the Hawaiian hotspot, the paleolatitudes of extinct volcanic edifices of the Emperor chain should match the presentday latitude of Hawaii (∼19°N) if the hotspot has remained fixed with respect to Earth's spin axis. The most reliable indicators of paleolatitude are basaltic rocks, but their reliability depends on each section spanning enough time to sample geomagnetic secular variation. Recovery of such samples requires ocean-drilling technology, and only a few seamounts have been sampled to date.

Paleomagnetic analyses of 81-million-year-old basalt recovered from Detroit Seamount (Site 884) yielded a paleolatitude of ∼36°N (10), which is discordant with Hawaii. Data from ∼61-million-year-old basalt (9) from Suiko Seamount define a paleolatitude of 27°N (11). These data suggest that the Emperor Seamounts record southward motion of the hotspot plume in the mantle (10).

A paleomagnetic test. The Ocean Drilling Program (ODP) Leg 197 (12) sought to test the hypothesis of southward motion of the Hawaiian hotspot by drilling additional basement sites in the Emperor chain (Fig. 1). We collected detailed stepwise alternating field (AF) demagnetization data aboard the drilling ship JOIDES Resolution. Although these shipboard data are of high resolution, they alone are insufficient to define paleolatitudes. Magnetic minerals with intermediate to high coercivities, carrying magnetizations resistant to AF demagnetization, are commonly formed during subaerial or seafloor weathering. The magnetizations of these mineral phases are easily resolvable in thermal demagnetization data, which we also discuss here (13).

Fig. 1.

Hawaiian-Emperor chain shown with ODP Leg 197 sites (12) and marine magnetic-anomaly identifications (40).

The geomagnetic field at a radius r, colatitude θ, and longitude ϕ can be described by the gradient of the scalar potential (Φ): Math Math where Math are partially normalized Schmidt functions, re is the radius of Earth, and the Gauss coefficients Math and Math describe the size of spatially varying fields. At least several millennia must be sampled such that the axial dipole term (Math) becomes dominant, allowing an estimate of the paleolatitude.

We used the angular dispersion of inclination averages (14) from independent lava flows (inclination groups) and compared this with global lava data of the same age (15) to examine whether secular variation has been adequately sampled (13, 16). Observations of the physical aspects of the lava flows, as well as petrologic and geochemical data, were used to group cooling units into lava units (12, 13).

Marine sediments can also provide useful paleolatitude information, but they generally provide only minimum values of paleolatitude because of potential compaction-induced inclination flattening (17). Sediment magnetizations, however, can record long time intervals. Similarly, chemical remanent magnetizations (CRMs), carried by minerals formed during weathering, can preserve stable magnetizations that provide insight into the time-averaged field.

Koko Seamount (Site 1206). Site 1206 (Fig. 1) was positioned on the southeastern side of the lower summit terrace on Koko Seamount with the use of crossing underway seismic profiles (fig. S1). The base of the thin sediment cover contains nannofossils of Zones NP14 and NP15, which provide an early- to middle-Eocene minimum age for the volcanic basement [43.5 to 49.7 Ma (18)]. Fifteen volcanic formations (including pahoehoe flows, flow foot breccias, and subaerial a,a units) were recovered in 278 m of basement penetration (12). Thin intercalations of limestone, volcaniclastic sandstone, and a deeply weathered flow top were also recovered, providing geological evidence of time between lava flow units. The lavas are mainly of tholeiitic composition, although two alkalic flows were noted. Plateaus in 40Ar/39Ar incremental heating spectra from six whole-rock samples yield a mean age of 49.15 ± 0.21 Ma (2σ uncertainty quoted; table S1).

Reflected-light microscopy revealed relatively fresh titanomagnetite grains in the lavas (fig. S2A). AF demagnetization of basalt samples (n = 74) showed the removal of a small low-coercivity component, followed by the definition of a characteristic remanent magnetization (12). Two normal-polarity intervals were defined, separated by a thin reversed-polarity zone; the radiometric age data are most compatible with the chron 21n-21r-22n sequence. This assignment is confirmed by thermal demagnetization (n = 113), which revealed univectorial decay after the removal of a small overprint (fig. S2, B and C).

Seventeen inclination groups were identified in the thermal demagnetization data (table S2), with a mean thermal inclination (Math; hereafter, all uncertainty regions are the 95% confidence interval unless otherwise noted) that was nearly identical to that isolated by AF treatment (Math, based on the 14 inclination groups sampled; (fig. S2D). The angular dispersion (13) of the thermal data (Math) is within the error of that predicted by global lava flows of 45 to 80 Ma (15). A comparison of the inclination units based on thermal demagnetization with a synthetic Fisher distribution (19) (Fig. 2) suggests that the basalt sequence well represents the time-averaged geomagnetic field. Furthermore, a stable magnetization carried by hematite [likely a CRM with unblocking temperatures (TUB) > 580°C] in samples of the deeply weathered basalt yields a mean inclination (Math, n = 10) that is indistinguishable from that of the lava flows.

Fig. 2.

Paleomagnetic inclination groups for Sites 1206 (A), 1205 (B), 1204 (C), and 1203 (D), based on thermal demagnetization data (blue) shown compared with a synthetic Fisher distribution (magenta) that has the same mean as the experimental data and a dispersion predicted by global lava data (15). Latitude of Hawaii, mean basalt inclination value (thermal demagnetization) for Hole A, Site 1204 [(C), green], and mean sediment inclination value (thermal demagnetization) for Site 1203 (D) are also shown.

Nintoku Seamount (Site 1205). Site 1205 on Nintoku Seamount (Fig. 1) was selected with the use of two crossing seismic profiles that defined a ∼43-m-thick sedimentary sequence above a flat igneous basement on the northwestern edge of the summit region (fig. S1). Nannofossils of Zone NP10 in the sediment immediately overlying the basement provide a minimum age [53.6 to 54.7 Ma (18)] similar to that obtained by radiometric dating of basalt from the nearby Deep Sea Drilling Project (DSDP) Site 432 [56.2 ± 0.6 Ma; 1σ error (20)].

A sequence of 25 subaerially erupted a'a and pahoehoe lavas and interbedded sediment and soil horizons were recovered in 283 m of basement penetration (12). Two flows of tholeiitic basalt are intercalated within the lower part of the sequence, which is otherwise composed of alkalic basalt. Plateaus in 40Ar/39Ar incremental heating spectra from six whole-rock basalt samples (table S1) spanning the section give a mean age of 55.59 ± 0.25 Ma.

The recovered basalt contains fresh to slightly altered titanomagnetite grains (fig. S2E). AF demagnetization (n = 144) defined a reversed polarity magnetization after the removal of a normal polarity overprint. Biostratigraphic and radiometric age constraints indicate that the basalt was erupted during chron 24r (18). Thermal demagnetization of basalt samples (n = 85) in general showed univectorial decay after the removal of a small overprint (fig. S2F). Unblocking temperatures sometimes extended to 625°C, indicating the presence of some hematite. Low unblocking temperatures (<325°C) were seen in other samples, indicating the dominance of high-Ti titanomagnetite carriers (fig. S2G) (13).

The thermal demagnetization data (22 inclination groups) yield a mean (Math) (table S3) that is similar to that of the AF data (Math) (fig. S2H). The estimated angular dispersion of the thermal data (Math) is higher than that predicted by models (15) based on the average data of lavas from 45 to 80 Ma, and may point to higher frequency changes in secular variation with time. Nevertheless, the geologic evidence for time, together with a comparison of the data versus a Fisherian distribution (Fig. 2), suggests that the mean value represents the time-averaged field.

Detroit Seamount (Site 1204). Site 1204 (Fig. 1) was chosen along two crossing seismic profiles that revealed a flat basement surface beneath ∼850 m of sediments (fig. S1). Two holes were drilled (A and B), penetrating 60 and 138.5 m of basalt, respectively. Diamictite and volcanic ash–rich sediment directly overlying the basement contain Campanian nannofossils (CC22/23) of ∼73 to 76 Ma (19); nannofossils of this zone were also recovered in a sediment interbed in the last core from Hole B, suggesting that the entire basement section drilled is of this age (12).

The basement units are pillow lavas consisting of multiple lobes (which are sometimes separated by calcareous sediment) intercalated with a 46-m-thick interval of massive lava (diabase) and a 10-m-thick hyaloclastite lapilli breccia interval (12). The basalt has undergone low-temperature alteration, consistent with late-stage seafloor weathering. 40Ar/39Ar incremental heating spectra from four whole-rock basalt samples and a feldspar separate failed to show plateaus in the step-heating ages. The apparent ages for total fusion range from 33 to 58 Ma, indicating variable loss of radiogenic 40Ar (table S1).

Reflected-light microscopy showed a dominance of titanomaghemite in basalt samples (fig. S2I). Linear decay during AF demagnetization was observed in only ∼70% of the samples examined (Hole A: n = 9; Hole B: n = 39). In these cases, a normal polarity magnetization was defined (12). Thermal demagnetization, however, revealed a clear two-component structure in many samples (13). A reversed-polarity component was defined at low unblocking temperatures (TUB ranging from 100°C to between 275° and 350°C), whereas a normal-polarity component was defined at higher temperatures (fig. S2, J and K). We interpret the high TUB component as the primary magnetization because it is the only magnetization seen in the least altered rocks after the removal of viscous components; we believe the reversed-polarity component is a later magnetization carried by titanomaghemite.

We identified only a single inclination unit in Hole A. Thermal demagnetization data (n = 10) yield a mean inclination (Math) that differs from that observed by AF treatment (Math), presumably because the latter fails to separate the two components of magnetization. On the basis of a lithofacies succession in Hole B (12), we identified five inclination groups. Our thermal results (n = 36) indicate a mean inclination of Math, based on AF data) for Hole B (table S4 and fig. S2L). Angular-dispersion estimates are low (Math) and at face value suggest that the sequence does not adequately sample secular variation (Fig. 2). However, rather than random directions, AF and thermal demagnetization of breccia samples often revealed a consistent normal-polarity direction, similar to that of the basalt flows. This observation, together with the presence of titanomaghemite, suggests that CRMs might be important at this site.

Detroit Seamount (Site 1203). Site 1203 (Fig. 1) on Detroit Seamount was located in an area of flat basement imaged in three crossing seismic lines (fig. S1). A thick (462-m) pelagic sedimentary sequence overlies the basement. Eighteen compound lava-flow units and 14 volcaniclastic sedimentary interbeds were drilled in 453 m of basement penetration (12). Nannofossils of Zone CC22 were identified in sedimentary beds within the basement sequence, indicating an age of 75 to 76 Ma (12). Plateaus in 40Ar/39Ar incremental heating spectra from three wholerock basalt samples and two feldspar separates (table S1) yield a mean age of 75.82 ± 0.62 Ma.

The upper part of the basement sequence consists of nonvesicular pillow lavas and thick, sparsely vesicular pahoehoe lava flows. These lavas were deposited far from eruptive centers at relatively shallow depths (fig. S1). The lower part of the sequence is dominated by highly vesicular compound pahoehoe lavas (up to 65 m thick) but also contains pillow lavas and hyaloclastite, indicative of shallow marine to subaerial environments. Tholeiitic to transitional basalt dominates the section, but some alkali basalt is intercalated in the lower sequence. This occurrence of alkalic lavas is unexpected; it may indicate the interfingering of volcanoes in different stages of evolution, as can be seen today between Mauna Loa and Mauna Kea on Hawaii.

The lavas have a range of magnetic mineralogies, from relatively fresh titanomagnetite to grains showing titanomaghemitization (fig. S2M). AF demagnetization of most of the basalt samples (n = 199) and sediment samples (n = 34) that we examined showed a stable single component of normal polarity after the removal of a low-coercivity (<10-mT) overprint (12). In the two uppermost flows, AF demagnetization patterns were more complex, indicating the presence of additional, unresolved components. Titanomaghemite, detected by reflected light microscopy, is common in the upper flows.

Thermal demagnetization data of lava samples (n = 87) (fig. S2N) and sediment samples (n = 28) (fig. S2O) mostly showed univectorial decay after the removal of a small overprint at low unblocking temperatures. In the two uppermost flows of the sequence, however, a distinct reversed component of magnetization was isolated (TUB ranging from ∼100°C to between 250° and 325°C) before definition of a normal polarity magnetization at higher temperatures. As in the Site 1204 basalts, we interpret the higher TUB magnetization as the primary remanence.

Thermal demagnetization results from lava samples (16 inclination groups) yield a mean (Math; table S5) similar to that calculated from the AF data (Math, for 14 inclination units, excluding the uppermost lava flows; fig. S2P). The estimated angular dispersion of the thermal demagnetization data (Math) is slightly higher than that expected from global lava-flow data (15) (Fig. 2). The inclination average from AF demagnetization of sediment samples (Math) (12, 13) is confirmed by thermal demagnetization (Math) and is somewhat steeper than the basalt mean.

Paleolatitude history. The inclination groups, averaged by site, form a progressive sequence of decreasing paleolatitudes with time (Fig. 3) that is inconsistent with the fixed-hotspot hypothesis. We did not recover coral reef material north of Koko Guyot (nor did we find evidence of such material as sedimentary debris). This is consistent with the idea that the hotspot was once located farther north, beyond the latitudinal zone supporting reef growth (21).

Fig. 3.

(A) Paleolatitude data from ODP Leg 197 sites (1206, Koko Seamount; 1205, Nintoku Seamount; 1204B and 1203, Detroit Seamount), ODP Site 884 (Detroit Seamount) (10), and DSDP Site 433 (Suiko Seamount) (11). Orange, results of thermal demagnetization; blue, results of alternating field demagnetization. Result from 433 is based on AF and thermal data. Magenta, magnetization carried by hematite from weathered basalt from Site 1206. (B) Average paleolatitude value for Detroit Seamount (square), based on inclination groups derived from basalts of Sites 884, 1203, and 1204B (Model B, see text), plotted with select values from other seamounts (A). Also shown is a least-squares fit to the data (orange) and several paleolatitude trajectories representing combinations of plate and hotspot motion.

Four data sets are now available from Detroit Seamount based on thermal demagnetization (Fig. 3). Lava emplacement may have been less frequent at Site 884 relative to the other sites because of its flank position. The Site 1204 (Hole B) lavas record a low angular dispersion, but these rocks might carry a CRM, explaining the agreement of their mean inclination with that of the Site 884 basalt section and the Site 1203 sediments. Because of potential inclination shallowing, the mean inclination from the Site 1203 sediments should be a minimum. This suggests further that the mean derived from the basalts at the same site is shallower because the available lavas underrepresent higher inclination values.

We consider two scenarios: one in which the paleomagnetic results from the Site 884 basalts, 1204 (Hole B) basalts, and Site 1203 sediments best represent the field (paleolatitude model A; Math, N = 3) and another in which we combine all the individual basalt inclination units from Detroit Seamount into a mean (Model B; Math, N = 32). With either model, the paleolatitude and age data yield average rates (Model A: 57.7 ± 19.2 mm year1; Model B: 43.1 ± 22.6 mm year1) that are consistent with the hypothesis that the Hawaiian hotspot moved rapidly southward from 81 to 47 Ma (10). The values are consistent with updated estimates of hotspot motion based on independent relative plate motions (5). Both paleolatitude models suggest that most of the motion occurred before the time of the Hawaiian-Emperor bend (Model A, >44 Ma; Model B, >43 Ma). This is further supported by the mean paleolatitude value from Koko Seamount (based on both thermoremanent magnetizations and CRMs), which is only 2.5° north of the fixed-hotspot prediction.

Modeling of hotspot motion. Crust ages available from marine magnetic anomalies, radiometric age data from drill sites, and geochemical data (22) indicate that the Hawaiian hotspot was close to a spreading ridge during the formation of Detroit Seamount (23). Hence, asthenospheric channeling of the plume (24) from a position to the south toward a more northerly ridge could have played some role in the difference between the paleomagnetic data and the prediction of the fixed-hotspot model. The monotonic age progression of lavas recovered from Detroit to Koko Seamounts, however, leads us to believe that this potential channeling of plume material was limited to the region at or north of Detroit Seamount. Furthermore, the similarity of the Hawaiian-Emperor chain with the Louisville chain of the South Pacific suggests that asthenospheric channeling was not the sole cause of the paleolatitude progression.

We examined whether the observed paleolatitude motion can be explained by a geodynamic model of the interaction of a plume with large-scale mantle flow. The flow calculation (25) requires a mantle density and viscosity model and a surface-velocity boundary condition (13). A tomography model was used to infer mantle density variations (26), and a viscosity structure based on an optimized fit to the geoid (with additional constraints from heat flow) (27) was applied. Both moving- and fixed-plume sources (28, 29) that originate at the top of the low-viscosity layer at the base of the mantle were considered (13).

Fast motion occurs when a conduit is sheared and tilted in the large-scale flow and the tilted conduit rises to the surface aided by a large-scale upwelling. Thus, in these models, fast hotspot motion corresponds to slower mantle flow rates (∼10 to 20 mm year1). Most computations (25) yield a hotspot motion of 5° to 10° toward the south to southeast during the past 100 million years (My). It is possible to achieve a good fit to the paleomagnetic data, because the age of the initiation of the Hawaii hotspot is unknown (and can hence be used as a free parameter). For the moving-source model, southward motion tends to be faster if an earlier plume origin is assumed. However, plume initiation ages from 180 to 120 Ma (which imply the subduction of volcanic edifices older than the oldest extant seamount, Meiji Guyot) generally yield the best fits. For the fixed-source model, the computed hotspot motion consists of two distinct phases. During the first phase, which lasts 100 to 150 My, southward motion can be rapid. The second phase begins when the first conduit elements that arise from the fixed source reach the surface. The computed hotspot motion is slow during the second phase and, in the example shown (Fig. 4), somewhat toward the north (13).

Fig. 4.

(A) Computed Hawaiian hotspot motion for fixed-source model (colored line), and tracks for fixed-source models (continuous line; plume initiation at 160 Ma) and moving-source models (dashed line; plume initiation at 170 Ma). Tickmark interval is 10 Ma for both. (B) Computed changes of hotspot latitude for fixed-source plume model (13) (continuous red lines) for plume-initiation ages of 150, 160, and 170 Ma (upper to lower). Moving-source model (13) results (dashed purple lines) are shown for plume initiation at 180, 170, and 160 Ma (upper to lower). Paleolatitude means for Koko, Nintoku, Suiko, and Detroit Seamounts (Fig. 3) are also shown.

Overall, the results of the large-scale flow modeling approach described above are consistent with the Leg 197 paleomagnetic data. Potentially important differences lie in the total motion predicted since ∼80 Ma (13) and in the need to incorporate in the modeling results a change in plate motion at or near the time of the bend. The paleomagnetic data do not require a change in plate motion, although a small change is not excluded.

Implications of hotspot motion. The hotspot motion defined by the new paleomagnetic and radiometric age data has implications for a wide variety of issues, including true polar wander (TPW) (30), the morphology of the past geomagnetic field, and the history of plate motions. Some investigators (31) have proposed that as much as 30° of TPW (rotation of the entire solid Earth) has accumulated during the past 200 My. However, a fixed-hotspot reference frame is used to define TPW in these studies. The data presented here, together with other tests (32, 33), indicate that TPW has been overestimated; Earth has been relatively stable with respect to the spin axis since the Early Cretaceous Epoch. Similarly, some changes in the morphology of the geomagnetic field with time (34) that have relied on fixed hotspots to anchor data from global sites are probably artificial. One recent analysis that has not relied on the fixed-hotspot reference frame has called for a significant axial octopole contribution (Math) to the time-averaged field (35). This conclusion is controversial, but if correct, it would imply that our paleolatitude calculations underestimate the true hotspot motion.

Backtracking the position of early-Tertiary and older Pacific basin sites, an essential aspect of some paleoclimate and tectonic studies, requires rethought, given that previous efforts have also relied on fixed hotspots. The northerly position of the Late Cretaceous Hawaiian hotspot (23) casts doubt on the southern option for the Kula-Farallon ridge [a plate configuration that is typically called on to create high rates of northward transport for tectonostratigraphic terranes in Alaska and British Columbia (36, 37)].

The fixed-hotspot interpretation of the Hawaiian-Emperor bend implies that huge plates can undergo large changes in direction rapidly. But such changes cannot be associated with internal buoyancy forces (such as subducting slabs) because these require many millions of years to develop. This has led to the suggestion that plate-boundary forces might be responsible (38). The new paleolatitude and radiometric age data (9) suggest that changes of plate motion at the time of the Hawaiian-Emperor bend were much smaller and more gradual than previously thought. Given the central role the Hawaiian-Emperor bend has played as an example of plate motion change, these observations now raise the question of whether major plates can undergo large changes in direction rapidly, and whether plate boundary forces alone can play a dominant role in controlling plate motion.

The similarity of the Hawaiian and Louisville hotspot tracks implies that the motion we are tracking by the new paleomagnetic data is of large scale. This Late Cretaceous to early-Tertiary episode of hotspot motion was not isolated; motion of the Atlantic hotspots relative to those in the Pacific occurred at similar rates during mid-Cretaceous times (39). These data sets indicate a much more active role of mantle convection in controlling the distribution of volcanic islands. At times, it is this large-scale mantle convection that is the principal signal recorded by hotspot tracks.

Supporting Online Material

www.sciencemag.org/cgi/content/full/1086442/DC1

Materials and Methods

Figs. S1 and S2

Tables S1 to S5

References and Notes

View Abstract

Navigate This Article