Phase change in subducted lithosphere , impulse , and quantizing Earth surface deformations

The new paradigm of plate tectonics began in 1960 with Harry H. Hess’s 1960 realization that new ocean floor was being created today and is not everywhere of Precambrian age as previously thought. In the following decades an unprecedented coming together of bathymetric, topographic, magnetic, gravity, seismicity, seismic profiling data occurred, all supporting and building upon the concept of plate tectonics. Most investigators accepted the premise that there was no net torque amongst the plates. Bowin (2010) demonstrated that plates accelerated and decelerated at rates 10 times smaller than plate velocities, and that globally angular momentum is conserved by plate tectonic motions, but few appeared to note its existence. Here we first summarize how we separate where different mass sources may lie within the Earth and how we can estimate their mass. The Earth’s greatest mass anomalies arise from topography of the boundary between the metallic nickel–iron core and the silicate mantle that dominate the Earth’s spherical harmonic degree 2 and 3 potential field coefficients, and overwhelm all other internal mass anomalies. The mass anomalies due to phase changes in olivine and pyroxene in subducted lithosphere are hidden within the spherical harmonic degree 4–10 packet, and are an order of magnitude smaller than those from the core–mantle boundary. Then we explore the geometry of the Emperor and Hawaiian seamount chains and the 60 bend between them that aids in documenting the slow acceleration during both the Pacific Plate’s northward motion that formed the Emperor seamount chain and its westward motion that formed the Hawaiian seamount chain, but it decelerated at the time of the bend (46 Myr). Although the 60 change in direction of the Pacific Plate at of the bend, there appears to have been nary a pause in a passive spreading history for the North Atlantic Plate, for example. This, too, supports phase change being the single driver for plate tectonics and conservation of angular momentum. Since mountain building we now know results from changes in momentum, we have calculated an experimental deformation index value (1–1000) based on a world topographic grid at 5 arcmin spacing and displayed those results for viewing.


Introduction
Gravity, the weakest force in our universe, has also provided us with one of the hidden keys for understanding plate tectonics.We say "hidden" because the deep positive mass anomalies within the Earth that drive the motions of the planet's tectonic plates are an order of magnitude smaller than, and hence obcured by, the Earth's greatest mass anomalies (see Fig. 5).Those greatest mass anomalies are of 0.4−3.5 × 10 23 kg, which we, as Bowin (2000), interpreted as due to topography at the core-mantle boundary (CBM) at 2900 km depth.The density contrast at the core-mantle boundary, with nickel iron of the core below and silicates of the lower mantle above, is approximately 4.3 gm/cc.Thus mass anomalies due to topography of up to 3 km relief there overwhelm all other mass anomaly contributions within the Earth.The power spectrum for the gravitational field of the Earth (see Fig. 1) shows a marked break in slope between degree 3 and degree 4. A global map view of the geoid anomalies (gravitational potential field) contained in the spherical Published by Copernicus Publications on behalf of the European Geosciences Union.
Figure 1.Power spectra for Earth, Venus, and Mars from sets of spherical harmonic coefficients.Earth (solid line) is from GEM-L2 (Lerch et al., 1982) and GEM-t2 (Marsh et al., 1990), Venus (dotted line) is from Nerem et al. (1993), and Mars (dashed line) is from Smith et al. (1993).If the horizontal harmonic degree axis were plotted as a log scale, then the spectra would be nearly a straight line.
harmonic coefficients of the Goddard Earth Model (GEM) t2 set of coefficients is shown in Fig. 2. That view of the gravitational potential field of the Earth has been a puzzle since it was first observed because it shows no correlation with surface continents nor topography, with the single exception of a positive geoid anomaly over the Andes of South America (SA).In Fig. 2 note the prominent negative geoid low in the Indian Ocean, the Sri Lanka (SL) Low, and the positive geoid, the New Guinea High (NG).And then in Fig. 3, note that the large summed contributions of degrees 2 and 3 are restricted to only those two NG and SL geoid anomaly sites.In Fig. 3, one can also see that the individual degree contributions become decreasingly smaller above degree 10.Thus, we now approach the Earth's second-largest mass anomaly as it remains buried within the spherical harmonic degree 4-10 packet.Besides illustrating the spatial distribution of the degree 4-10 geoid anomalies, Fig. 4 further demonstrates how the large contributions from the degree 2-3 dominate the 2-10, 2-30, and 2-250 geoid (seen in Fig. 2) fields.Further, Fig. 4 also visually introduces the use of vertical derivatives of the Earth's potential field into our discussion for resolving Earth structures.In the gravity degree 2-3 map note the paucity of contour lines compared to those in the degree 2-3 geoid map.Geoid anomalies (1 d −1 ) are most diagnostic of large mass anomalies at depths of thousands of kilometers, whereas gravity anomalies (1 d −2 ) are most sensitive to mass anomalies within 100 km of the surface.Consider a site directly over a point mass, the geoid anomaly (N) is where V is the disturbing potential, g 0 is normal gravity (9.8 m s −2 ), G is the gravitational constant, M is the anomalous mass, and z is the depth of the point mass.The vertical component of gravity g due to the same point mass is The ratio of gravity anomaly to the geoid anomaly directly above the point mass at depth z is then Conversely, the depth can be determined by Then, of course, reentering the value of z into either, or both, of the point-mass equations for the geoid or gravity can provide estimates of the causative mass.It was by this means that the magnitudes of the mass anomalies of the Earth shown in Fig. 5 were estimated.Two additional vertical derivatives of the gravity potential can also aid in defining shallower mass anomaly structures: the vertical gravity gradient (1 d −3 ) and the vertical gradient of the vertical gravity gradient (1 d −4 ).
For each mass element within the Earth, its geoid contribution falls off as 1 d −1 with distance, gravity as 1 d −2 , vertical gravity gradient 1 d −3 , and vertical gravity gradient of vertical gravity gradient 1 d −4 .But because all masses within the Earth contribute to all harmonic degrees, it is impossible to invert any one field into its causative masses.However, at selective locations the ratio of gravity anomaly divided by the geoid anomaly value at center locations (Fig. 5) has already provided constraints on the internal mass anomaly structure of the Earth (Bowin, 1983(Bowin, , 1986(Bowin, , 1991(Bowin, , 2000)).At lithospheric subduction sites (positive geoid anomaly band in degree 4-10 map in Fig. 4), the phase changes of olivine and pyroxene grains in the subducted lithosphere are inferred to produce the positive density contrast with the surrounding upper mantle, which forms the second-greatest mass anomalies within the Earth (2.0−8 × 10 22 kg).
For many decades there has been agreement that viscous flow within the mantle could theoretically lead to topographic and mass perturbations at the Earth's surface and at density horizons at depth (Pekeris, 1935;Hager, 1984).As the blog "Retos Terricolas" (5 June 2014), "Dynamic topography vs. isostasy: The importance of definitions", nicely points out, dynamic topography was originally coined by oceanographers Bruce (1968) and Wyrtki (1975) to refer to deviations of the ocean surface relative to the geoid after subtracting effects due to tides, waves, and wind by time filtering.The remaining deviations of "dynamic topography" related to water currents in the oceans are the oceanographers goal.In the 1980s the term dynamic topography was adopted by solid-Earth scientists but included "static" forces, such as the weight of sinking plates.Hager and O'Connell (1981), Hager et al. (1985), and Ricard et al. (1993) inferred a triplet of mass anomalies by adding negative downward deflections of the ocean surface and at the core-mantle boundary below due to mantle flow.In this scenario the positive mass of the sinking slab would become masked by both of those two dynamic flow-induced negative depressions.For a mass at a single depth, to be the source for a geoid anomaly would require all of its vertical derivatives to have the same equivalent point-mass depths.To utilize that fact for a test, Bowin (1994Bowin ( , 2000) ) computed spherical harmonic coefficients from balanced sets at six different depths (100, 500, 1000, 2000, 2900, 4000 km).This was done for balanced sets of two positive and two negative masses and another balanced set of three positive and three negative masses to provide for both even-and odd-order harmonic values.By then calculating and displaying calculated cumulative contribution curves (CCC) (see, e.g., Bowin (2000), Fig. 12) plotted as percentages of their maximum value, estimates of depths to point-mass sources could be estimated for each potential derivative.The result for the SA geoid high shown in Fig. 6 yielded a consistent point-mass depth of nearly 1200 km depth for all harmonic degrees and all derivatives, indicating a single source as its origin and hence a shallower depth for the actual positive mass anomaly.This 1200 km point-mass depth is also consistent with the SA geoid high source being shallower above the zero contour line at 2900 km depth of the core-mantle topography (Fig. 4, degree 2-3 geoid).This strong evidence for a single source depth for the South American geoid high is strong support opposing dynamic topography as being a significant contributor to the Earth's geoid.

Deriving high derivative Earth potential field model
The gravity field coefficients are estimated with high accuracy in the long-wavelength part (corresponding to the coefficients of low degrees).In order to avoid filtering these low degree coefficients, we shift the damping factors above to a certain degree (30) and obtain a new set of filtered coefficients that have yielded stable estimates of global geoid, gravity, vertical gravity gradient, and vertical gradient of vertical gravity gradient anomaly values for degrees 31-264.These coefficients reflect mass anomalies that exist in the outer 200 km of the Earth.Our derived filtered spherical harmonic coefficients for degrees 2-264 are given in the accompanying Supplement.
Various alternative methods exist for the determination of global gravity field models -such as integration on Earth sphere, least-square adjustment or collocation -and various data sources are used, such as satellite-based data or terrestrial (surface) observations.The error behaviors of the resulting gravity field models depend on the method, the quality and properties of the observations.In general, the gravity field is not homogeneously determined.With satellite gravimetry the high-frequency part (corresponding to the gravity field coefficients of higher degrees) is computed less accurately than the low-frequency part.This results in some short-wavelength noise in the computed gravity field models.It is therefore necessary to apply filtering to extract the signal by suppressing the accompanying noise.
The idea of filtering gravity field models in the spectral domain is based on smooth SH coefficients with some depression factors.Various smooth functions have been proposed such as Pellinen (1966) and Jekeli (1981).With a Gaussian mean filter in the spectral domain as an example, as presented by Jekeli (1981), the recursive formula for a Gaussian smoothing coefficients in the spectral domain is given as with r the radius of a spherical cap to be smoothed and R ⊕ the radius of the Earth.
Using the above filter coefficients, the filtered gravity field coefficients can be derived as The gravity field coefficients are estimated with high accuracy in the long-wavelength part (corresponding to the coefficients of low degrees).In order to avoid filtering these low degree coefficients, we shift the damping factors above to a certain degree n 0 and obtain a new set of filter coefficients such that The gravity field coefficients are estimated with high accuracy in the long-wavelength part (corresponding to the coefficients of low degrees).In order to avoid filtering these low degree coefficients, we shift the damping factors above to a certain degree (30) and obtain a new set of filter coefficients that have yielded stable estimates of global geoid, gravity, vertical gravity gradient, and vertical gradient of vertical gravity gradient anomaly values for degrees 31-264.These coefficients reflect mass anomalies that exist in the outer 200 km of the Earth.
Plots of our results are shown in Fig. 7 for degree 31-264 contributions for geoid, gravity, vertical gravity gradient, and vertical gradient of vertical gravity gradient together all demonstrate that the only principle mass anomalies in the outer 200 km of the Earth are those of the negative mass anomalies associated with fore-arc lithosphere depressions, and the positive mass anomalies of the uplifted island arc/mountain ranges.

Motions of the Earth's tectonic plates and earthquakes
In 1960 Harry H. Hess was the first to recognize that the ocean floor was not everywhere ancient (Precambrian) but is, in fact, locally, being created today (Hess, 1960; it was to be published in "The Sea, Ideas and Observations" but became delayed).Robert S. Dietz received a copy of the Hess (1960) preprint and published a paper, Dietz (1961), on seafloor spreading without acknowledging Hess's prior work.Dietz concluded that ocean basins most likely contained rocks no older than the Mesozoic, and he also noted that linear magnetic anomalies aligned normal to the presumed flow direction of the presumed convection cells were being found in the Atlantic.Hess then had his paper published in the next available Geological Society of America (GSA) publication volume available, Hess (1962).Magnetic anomaly stripes were reported in the Atlantic Ocean by Vine and Mathews (1963).Wilson (1963) in his Fig. 4 referenced forming the Hawaiian chain over a "stable core of a convection cell" as a possible origin of the Hawaiian chain of islands, whose source later was identified as being to hot spot plumes.Looking back with hindsight, we see a remarkable coming-together of oceanic magnetic anomaly stripes that allowed linkage to normal and reverse paleomagnetic age dating of the oceanic crust beneath, as well as the global patterns of rifts and transform faults.And, in turn, aided in the dating of folded and deformed rocks of the former or active subduction collision zones, the depressed sections at rifts, and the horizontal displacements along strike-slip faults, like that of the San Andreas, all fit into the new understanding of plate tectonics and its three types of boundaries.Another consequence of plates rubbing against one another while moving is to increase strains in the rocks of the opposing plates, which, when released, cause an earthquake.In this way, the distribution of epicenters reveals the locations of plates and their boundaries, as in Fig. 8.However, that figure only displays those having earthquakes having a magnitude 6 or greater; hence the less seismically active borders do not stand out.Essentially all the analyses of plate motions (Morgan, 1973;Soloman and Sleep, 1974;Forsyth and Uyeda, 1975;Chapple and Tullis, 1977) have assumed that plates are not undergoing acceleration.Accordingly, every plate must be in dynamic equilibrium, such that the sum of the torques about any axis must be zero.Previously, the paucity of clear evidence for systematic plate velocity changes through time has led researchers to view convective motion in the mantle, with traction on the overlying lithosphere, as likely producing plate tectonics, in which plates would move in stages with a terminal velocity.Different stages might have different plate motion direction and/or speed.Harada and Hamano (2000) made a simple but profound assumption: that hot spots on the Pacific (PA) Plate were "fixed" in time and space relative to each other.With  2003) and earthquake locations with magnitudes greater than magnitude 6.0.
that assumption, they formed a spherical triangle amongst the center points of the Hawaiian, Louisville, and Easter Island hotspots.This triangle was then used to identify which digitized location points along the Emperor-Hawaiian seamount chain correlated with which digitized points along the Louisville ridge chain.Then, with those location points correlated in time, they could begin estimating poles of rotation.They then took each location point along the Emperor-Hawaiian chain, and bisected its arc with the Hawaii hotspot location and at that bisect point projected a line normal to it.They then did the same procedure to the equivalent point along the Louisville chain with its bend, and where the two normal lines intersect is an equivalent stage pole for that selected point.From ages of dredged samples, they estimated Euler poles at 4 and 2 Myr intervals.Figure 9, from Bowin and Kuiper (2005) using the Harada and Hamano (2000) Euler poles, had the Pacific Plate speeding up as it moved northward toward the Aleutians, then after the bend slowed down, and then again increased speed as it moved towards the Marianas (away from Hawaii).
Bowin (2010) demonstrated that plate tectonics conserves angular momentum (see Fig. 10) using the 14-majortectonic-plate history of the Earth (except the Juan de Fuca and Philippine plates).Euler stage pole histories of the past 62 million years of Gripp and Gordon (1990) were given.We believe it is worth restating that "mirroring of the bend in the Emperor-Hawaiian seamount chain in the locations of the 4448 filtered stage pole locations for the absolute motions of the Nazca Plate (Bowin, 2010;Plate No. 11, Fig. 12) gives credence to both the quaternion analyses utilized and the role of impulse in plate motions".Of the 14 plates analyzed, the Nazca Plate is the only one that has a stage pole pattern revealing such a pronounced bend in stage pole trend at 46 Myr.
For this present study of active surface deformation we have chosen to use a more detailed 52 tectonic plate model of the Earth's surface by Bird (2000).Bird's published Euler pole values are relative to a fixed Pacific Plate.An important contribution was provided by Rick Rosson (Table 1) of Mathworks.com,who converted Bird's 52 relative Euler poles to absolute poles of rotation (Table 1) Although the rates of plate acceleration/deceleration are very small, they are real and thus indicate that impulse (force times time) due to changes in plate momentums (plate mass times velocity) are what causes deformations within the Earth and on its surface.Thus one of our early tasks was to display the individual plate motion absolute velocity vectors for the 52 plates utilizing the absolute rotation Euler pole data of Table 1 on a global rectangular grid having 5 arcmin spacing (see Fig. 8).
Earthquakes occur as a result from the release of strain that accumulates in the Earth, and those strains, in turn, result from plate motions.Therefore, this dependency proves that the maximum energy for seismicity is similarly bounded.Thus, we contend that this new knowledge of plate acceleration/deceleration opens a new portal through which earthquakes and lithospheric strain accumulation must now be viewed.We can now state that there will never be a magnitude 12 earthquake.We can also confidently state that there will never be a magnitude 11 earthquake.We cannot rule out the possibility that earthquakes of a very low magnitude 10 might occur.This is because the seismicity time constant for strain build-up in the crust is on the order of thousands of years, and we only have barely a century of seismicity observations.Furthermore, at the present GPS receiver sensitivities, nearly 100 years of GPS observations are required to directly assess present-day plate acceleration.Bowin (2010) demonstrated that it is the sinking of the positive phase change mass anomalies of the subducted lithosphere that drive plate tectonics.The present locations of those positive mass anomalies are most clearly revealed in the spherical harmonic degree 4-10 coefficient packet of the Earth's potential field.As the linear bands of subducted litho-sphere gradually shift locations, so, too, will the degree 4-10 and order coefficients change that delineate their location pattern within the degree 4-10 packet.Indeed, many questions of Earth history could be better answered if we knew how the Earth's past potential coefficients differed from today's.But, alas, we do not and cannot.According to classical geodesy we have to assume that the absolute motion of the Pacific Plate at any time results from instantaneous in-  Figure 8 provides a present-day sample of how the Earth is behaving after its motion changed from the Emperor seamount chain northward mantle motion era to the Hawaiian seamount construction era of westward mantle motion.In particular, note that the absolute easterly motion of the direction of the Pacific Plate, from the north end of the Tonga trench to the south end of the Yap trench, is essentially at 90 • to the northerly motion of the adjacent Australian (Au) Plate to the south.An inability to satisfactorily visualize the interface between a northern down-going subduction limb of an Australian thermal convection roller cell and a side-slipping transcurrent portion of a Pacific roller thermal-driven convection cell led Bowin (2010) to reject that scenario as not being credible.And bear in mind that this geometric setting would have persisted for more than 40 Myr.
Whereas Fig. 8 illustrates the absolute vector motions of the present 52 plates, Fig. 11 summarizes the global patterns of oceanic linear magnetic anomalies and transform faults, and thereby records a history of past events.Particularly striking is the apparent consistency shown by the magnetic record for the separation of North and South America from Africa with the opening of the Atlantic Ocean since the Middle Jurassic.This matching of the coastlines was noted in the initial ideas of continental drift by Wegener (1912).The bold-labeled numbers (2, 2a, 2b, 3a, 3b, 4a, 4b, 4c, 5) illustrate the migration progression of a pattern of plate propagation extended south from the South Pacific, separating Australia from Antarctica and continuing westward across the Indian Ocean into the Gulf of Aden and then up the Red Sea to the Gulf of Aqaba (Bowin, 1974).It was not until  nearly 36 years later that Bowin (2010) recognized that that beginning of that extension of the South Pacific spreading into the Indian Ocean realm began at the time of the bend in the Emperor-Hawaiian seamount chain (46 Myr).And it appears to have begun with the wedge-shaped opening (spreading) of the Tasman Sea off eastern Australia (labeled 2).Then the Pacific Plate spreading progressed southward (2a, 2b), to progressing westward, separating Australia from Antarctica (3a, 3b), and continuing through time to progress westward, moving India northward towards Eurasia, and then opening the Gulf of Aden and hence the Red Sea, to now the Gulf of Aqaba.A southwest-trending branch of a closely spaced young pattern of spreading and transform offsets can also be seen emanating from the triple-junction point in the Indian Ocean east of Madagascar.This branch coincides with that classified as ultraslow-spreading ocean ridges (Dick et al., 2003).
The remnant of a former Mesozoic triple junction remains evident in the magnetic anomaly pattern in the western Pacific Ocean basins (Fig. 9).In the eastern portion of the central and North Pacific Ocean basin, note the appearance of several south-north-trending propagations of axes of spreading there, with some dominating over a former trend.Besides the fact that the magnetic age reversal timescale has given us an extremely valuable record into past events, our brief overview of oceanic magnetic anomaly patterns has demon- strated that plate tectonic motions are not chaotic.The motions seem rather slow and steady.Even the 60 • change in azimuth of lithosphere flow that accompanied the "bend" in the Emperor-Hawaiian seamount chain trend near 46 Myr appeared to have not affected the spreading in the Atlantic Ocean, although it did coincide with the beginning of a significant rearrangement of continents in the Southern Hemisphere.The band of positive geoid anomalies in the degree 4-10 packet (Fig. 4, degrees 4-10) shows good agreement with the absolute plate motion vectors for present-day plate motions (Fig. 8), supporting the phase-change conclusion of Bowin (2010).Unfortunately, we cannot satisfactorily estimate what the Earth's spherical harmonic degree 4-10 coefficients were for the time period before the bend (say, 90-50 Kyr).Without such a quantitative understanding of a period before the bend with a period after the bend, the best estimates of rates of Earth's internal mass redistributions will remain those of Fig. 4.
We have now become forced to assume the Earth's spherical harmonic coefficients for degrees 4-10 must change with time to meet the evolving location pattern of dense phasechange subducted lithosphere and periodotic material.Perhaps one should consider that the coefficients for harmonic degrees 2-3 may also evolve with time, and that the locations of the Indian Low (Sri Lanka Low) and New Guinea High locations may also shift with time.But it it appears to be only the smaller mass anomalies of the falling dense masses of the phase change mass anomalies that conserve the angular momentum of the global plate tectonic system.The unanswered question now shifts to the following: how does such a planetary system of conservation of angular momentum become activated?

Quantifying a deformation index value
Since we can now view "mountain building" deformation as resulting from changes of momentum or impulse, how might such deformations be quantized?Topographic perturbations from a global mean value are a first-order measure of a degree of uplift or subsidence that has occurred in a region, but we sought a deformation index value that would be more relevant than simply a local topographic mean or standard deviation or a combination value.Our selection for this initial test of a deformation index value is to use, at each 5 arcmin pixel, the maximum absolute topographic difference with its eight adjacent neighboring pixels.We used this as a way to estimate the absolute maximum local topographic relief value at each grid location as a stand-in for a deformation index value.
Figure 12 presents our first global view of approximately 9 million maximum topographic relief values at 5 arcmin spacing using a color scale where the topographic relief in meters is divided by 5 so that the resulting deformation index value lies between 1 and 1000.Here, in addition to the orange and reds associated with the mountainous belts of the Himalayan and Andean ranges, and island arcs of the world, they also locate many of the undersea rises of the southwest Indian Ocean, western Pacific Ocean, and North and South Atlantic mid-ocean rises.And in the northwestern Pacific Ocean basin, a series of non-red nearly east-west-trending, quasi-equally spaced fracture zone features stand out.We are also struck (and puzzled) by the great abundance of yellow and orange colors from topographic relief in the western central Pacific region.It stands out as a unique part of the world.If deformation index values less than 50 are ignored, and those above 50 are sorted from lowest to highest, then a plot like Fig. 13 results.Figure 13 also helps clarify that over 6 million grid locations in the world map of Fig. 12 have deformation index values of less than 50 and hence are white on that map.
The Supplement related to this article is available online at doi:10.5194/se-6-1075-2015-supplement.

Figure 3 .
Figure 3. Magnitudes of individual harmonic degree contributions for the Earth's 10 major geoid anomaly values computed from the GEM-9 coefficients referenced with reciprocal flattening of 298.257 (actual Earth flattening).Four positive anomalies are New Guinea (NG), Iceland (I), Crozet (C), and South America (SA).Six negative anomalies are Indian Ocean (Sri Lanka) (SL), west of lower California (WLC), central Asia (CA), and south of New Zealand (SNZ).Reprinted from Bowin (1985) with permission.

Figure 5 .
Figure 5.Initial estimates of the magnitudes of Earth mass anomalies.Estimates are calculated from ratios of gravity anomaly divided by geoid anomaly values at anomaly centers.

Figure 6 .
Figure 6.Percentage contribution curve (CCC) and estimated equivalent point-mass depths for the South American (SA) GEM-t2 geoid high.Depths are estimated from the percentage CCC using spline function coefficients for each degree.Spline functions are from model sets of two and three balanced point masses at depths of 100, 500, 1000, 2900, and 4000 km.

Figure 7 .
Figure 7. Residual geoid, gravity, vertical gravity gradient, and vertical gradient of vertical gravity gradient for spherical harmonic degrees 31-464.Color scale weighted so that each color has same number of data samples.

Figure 8 .
Figure 8. Absolute plate velocity vectors (1 per 1000 plotted black arrows) for the 52 plates of Bird (2003) and earthquake locations with magnitudes greater than magnitude 6.0.

Figure 9 .
Figure 9. Stage pole velocities for Pacific Plate at labeled time on Emperor and Hawaiian seamount chains.From Bowin and Kuiper (2005).

Figure 11 .
Figure 11.Mercator map of oceanic magnetic anomalies and transform faults.Scanned copy of the original 1986 illustration.

Figure 12 .
Figure 12.World view of deformation index value (1-1000) based on a global topographic grid at 5 arcmin spacing.At each grid location, the maximum relief in meters of the location with the eight adjacent surrounding locations is divided by 5 to obtain its deformation index value.

Figure 13 .
Figure 13.Deformation index values sorted lowest to greatest value vs. sort point number.Note that more than the first 6 million values have values less than 50, and that accounts for the white areas in Fig. 12.

Table 1 .
Euler pole and angular velocity estimates; based on Pacific Plate absolute rotation fromMorgan and Morgan (2007).