Postglacial isostatic adjustment in a self-gravitating spherical earth with power-law rheology

Since microphysics cannot say definitively whether the rheology of the mantle is linear or non-linear, the aim of this paper is to constrain mantle rheology from observations related to the glacial isostatic adjustment (GIA) process—namely relative sea-levels (RSLs), land uplift rate from GPS and gravity-rate-of-change from GRACE. We consider three earth model types that can have power-law rheology (n = 3 or 4) in the upper mantle, the lower mantle or throughout the mantle. For each model type, a range of A parameter in the creep law will be explored and the predicted GIA responses will be compared to the observations to see which value of A has the potential to explain all the data simultaneously. The coupled Laplace finite-element (CLFE) method is used to calculate the response of a 3D spherical self-gravitating viscoelastic Earth to forcing by the ICE-4G ice history model with ocean loads in self-gravitating oceans. Results show that ice thickness in Laurentide needs to increase significantly or delayed by 2 ka, otherwise the predicted uplift rate, gravity rate-of-change and the amplitude of the RSL for sites inside the ice margin of Laurentide are too low to be able to explain the observations. However, the ice thickness elsewhere outside Laurentide needs to be slightly modified in order to explain the global RSL data outside Laurentide. If the ice model is modified in this way, then the results of this paper indicate that models with power-law rheology in the lower mantle (with A  10⁻³⁵ Pa⁻³ s⁻¹ for n = 3) have the highest potential to simultaneously explain all the observed RSL, uplift rate and gravity rate-of-change data than the other model types.     

The large volume multi-anvil press as a highP-T deformation apparatus

The rheological properties of mantle materials are being investigated up to pressures of 16 GPa and temperatures of 1600°C for times up to 24 h, using a new sample assembly for the 6–8 multi-anvil apparatus. Al₂O₃ pistons, together with a liquid confining medium, are used to generate deviatoric stress in the specimen. Strain rates are estimated by monitoring the relative displacement of the guide blocks of the multi-anvil apparatus, scaled to the total axial strain of the sample. The applied stress on the sample is estimated using grain size piezometry. Strain rates and flow stresses of approximately 10^–4 to 10^–6 s^–1 and 50 to 250 MPa respectively, are presently attainable.Preliminary results on San Carlos olivine single crystals, partially dynamically recrystallized to a grain size of 10 to 300 m, indicate that the effective viscosity of polycrystalline olivine is consistent with values obtained from olivine single crystal creep laws. Assuming a dislocation creep mechanism (n3.5) with (010)[001] as the dominant slip system, the data are best fit using a creep activation volume of 5 to 10×10^–6 m³ mol^–1.     

Reconciling strong slab pull and weak plate bending: The plate motion constraint on the strength of mantle slabs

Although subducting slabs undergo a bending deformation that resists tectonic plate motions, the magnitude of this resistance is not known because of poor constraints on slab strength. However, because slab bending slows the relatively rapid motions of oceanic plates, observed plate motions constrain the importance of bending. We estimated the slab pull force and the bending resistance globally for 207 subduction zone transects using new measurements of the bending curvature determined from slab seismicity. Predicting plate motions using a global mantle flow model, we constrain the viscosity of the bending slab to be at most ~ 300 times more viscous than the upper mantle; stronger slabs are intolerably slowed by the bending deformation. Weaker slabs, however, cannot transmit a pull force sufficient to explain rapid trenchward plate motions unless slabs stretch faster than seismically observed rates of ~ 10^− 15 s^− 1. The constrained bending viscosity (~ 2 × 10²³ Pa s) is larger than previous estimates that yielded similar or larger bending resistance (here ~ 25% of forces). This apparent discrepancy occurs because slabs bend more gently than previously thought, with an average radius of curvature of 390 km that permits subduction of strong slabs. This gentle bending may ultimately permit plate tectonics on Earth.     

Lateral variation in upper mantle viscosity: role of water

Differences in the viscosity of the earth's upper mantle beneath the western US (∼10¹⁸-10¹⁹ Pa s) and global average values based on glacial isostatic adjustment and other data (∼10²⁰-10²¹ Pa s) are generally ascribed to differences in temperature. We compile geochemical data on the water contents of western US lavas and mantle xenoliths, compare these data to water solubility in olivine, and calculate the corresponding effective viscosity of olivine, the major constituent of the upper mantle, using a power law creep rheological model. These data and calculations suggest that the low viscosities of the western US upper mantle reflect the combined effect of high water concentration and elevated temperature. The high water content of the western US upper mantle may reflect the long history of Farallon plate subduction, including flat slab subduction, which effectively advected water as far inland as the Colorado Plateau, hydrating and weakening the upper mantle.     

Small-scale convection under the back-arc occurring in the low viscosity wedge

Water released from subducting slabs through a dehydration reaction may lower the viscosity of the mantle significantly. Thus, we may expect a low viscosity wedge (LVW) above the subducting slabs. The LVW coupled with a large-scale flow induced by the subducting slabs may allow the existence of roll-like small-scale convection whose axis is normal to the strike of the plate boundary. Such a roll structure may explain the origin of along-arc variations of mantle temperature proposed recently in northeast Japan. We study this possibility using both 2D and 3D models with/without pressure- and temperature-dependent viscosity. 2D models without pressure and temperature dependence of viscosity show that, with a reasonable geometry of the LVW and subduction speed, small-scale convection is likely to occur when the viscosity of the LVW is less than 10¹⁹ Pa s. Corresponding 3D model studies reveal that the wavelength of rolls depends on the depth of the LVW. The inclusion of temperature-dependent viscosity requires the existence of further low viscosity in the LVW, since temperature dependence suppresses the instability of the cold thermal boundary layer. Pressure (i.e. depth) dependence coupled with temperature dependence of the viscosity promotes short wavelength instabilities. The model, which shows a relatively moderate viscosity decrease in the LVW (most of the LVW viscosity is 10¹⁸∼10¹⁹ Pa s) and a wavelength of roll ∼80 km, has a rather small activation energy and volume (∼130 kJ/mol and ∼4 cm³/mol) of the viscosity. This small activation energy and volume may be possible, if we regard them as an effective viscosity of non-linear rheology.     

Rheology of the lower mantle

The rheology of the lower mantle of the Earth is examined from the viewpoint of solid state physics. Recent developments in high-pressure research suggest that the lower mantle contains a considerable amount of (Mg, Fe)O with Fe/Mg + Fe = 0.2–0.3. The pressure and temperature dependences of diffusion in (Mg, Fe)O are estimated by the theory of diffusion in ionic solids. Of the materials composing the lower mantle, (Mg, Fe)O may be the “softest”, and therefore the rheology of the lower mantle may be that of (Mg, Fe)O, unless the framework effect is important.Temperatures in the lower mantle are inferred from the depths of phase transitions and the melting temperatures of the core materials. A thermal boundary layer at the base of the mantle is suggested. The physical mechanisms of creep are examined based on a grain size-stress relation and non-Newtonian flow is shown to be the dominant flow mechanism in the Earth's mantle.The effective viscosity for the temperature models, with and without the thermal boundary layer, is calculated for constant stress and constant strain rate (with depth). For constant strain rate, which may be appropriate for discussing the mechanics of descending slabs, the increase in effective viscosity with depth is smaller than for the constant-stress case, which may be appropriate for discussing the flow induced by the surface motion of plates.The relatively small depth gradient of viscosity, for constant strain rate, suggests that the lower mantle could also participate in convection. The effective viscosity increases with depth, however, by at least 10² to 10³ from the top to the bottom of the lower mantle, for a reasonable range of activation volumes and temperatures. There will be a low-viscosity layer at the base of the mantle, in contrast to the high-viscosity layer at the top of the mantle (plates), if a thermal boundary layer is present. The constant Newtonian viscosity inferred from rebound data may be an apparent feature resulting from the difference in deformation mechanisms between isostatic rebound and large-scale flow.     

Mantle-derived volatiles in continental crust: the Massif Central of France

CO₂-rich gases and groundwaters from springs and boreholes originating within the basement of the Massif Central have variable³He/⁴He ratios with correspondingR/R_a values ranging from 0.8 to 5.5 and 0.3 to 2.8 respectively, indicating the presence of a significant component of mantle helium. Molar concentrations of rare gases in the CO₂-rich gases are approximately 5 orders of magnitude greater than in the waters and suggest that a near-surface Henry's Law fractionation has occurred between exsolving CO₂ and water.δ¹³C values of the CO₂-rich gases are in the range −4.2 to −6.1‰, i.e. in that range normally attributed to mantle carbon, but which could also represent an average crustal composition and therefore do not discriminate between mantle and crustal sources.C/³He ratios show 4 orders of magnitude variation from 1.4 × 10¹² to 5 × 10⁸ and, compared to a mantleC/³He ratio of 10⁹, indicate that either a complex fractionation has occurred between mantle helium and mantle CO₂ or more likely that mantle rare gases have been diluted by large quantities of CO₂ with an average crustal carbon isotope composition. The regional distribution of³He and C does not show any obvious relationship to age or proximity of volcanic centres or major faults, suggesting that mantle-derived C and He components decoupled from their silicate melt sources at some depth.The results from this area of active fluid circulation suggest that C-isotope data derived from metamorphic terrains should be interpreted with great caution, but that input of some mantle-derived carbon is expected to accompany crustal extension.     

Hydration of an active shear zone: interactions between deformation, metasomatism and magmatism—the spinel-lherzolites from the Montferrier (southern France) Oligocene basalts

Geochemical and textural investigations have been simultaneously performed on spinel-lherzolite xenoliths from the Oligo-Miocene alkali basalts of Montferrier (southern France).All the investigated samples have undergone a deformation very particular by intense shearing under high stresses (up to 1.75 kbar), low temperatures ( 900°C) and strain rates of about 10⁻¹⁸ to 10⁻¹⁵ s⁻¹.Mineral chemistry reveals that the Montferrier lherzolites are fragments of an undepleted relatively shallow upper mantle level located at a depth of 50 km (15 kbar). Moreover, Na and Ti enrichment in diopside would reflect a metasomatic event, also emphasized by the common occurrence of pargasite in 50–70% of the investigated samples.Crystallization of this amphibole is attributed to a hydrous infiltration which is related in time and space to the deformation. Indeed, amphibole is preferentially concentrated in strongly deformed zones and in kink-band boundaries of orthopyroxene porphyroclasts. Moreover, the grain boundaries were used by the pervasive agent to percolate into the lherzolite: significant chemical variations (increase in MgO: 15% and decrease in Al₂O₃: 55%) are observed within the range of 7–5 μm adjacent to the grain boundary.Finally, Sr isotopic data (⁸⁷Sr/⁸⁶Sr) demonstrate that the amphibole, i.e. the metasomatic agent, is genetically related to the host lava of the xenoliths. Thus, the hydrous silicate liquid from which the amphibole has crystallized may be an early percolation of the ascending alkali magma.This silicate liquid hydrated the shear zone, located at a depth of 50 km, induced the hydraulic fracturation of the lherzolite and the magmatic conduit opening. Subsequently, the alkali magma sampled some fragments of this strongly deformed and metasomatized undepleted upper mantle level and brought them to the surface.     

One-dimensional modelling of mantle flow

A one-dimensional model of flow between a fixed boundary at the bottom and a moving one on top with no net flow through vertical sections is tested for geophysically interesting mantle viscosity-depth functions. Such a model, although simplistic, may help in answering the question to what depth the return flow extends, at least in the case of moving plates measuring many thousand kilometers across, such as the Pacific plate.It the viscosity in the asthenosphere is less than three orders of magnitude smaller than that of the mantle below, the return flow extends to great depth and the asthenosphere is a zone of concentrated shear. If the viscosity contrast is greater, the return flow is concentrated in the asthenosphere. For a wide range of model parameters typical flow velocities below the asthenosphere are about one-tenth of the plate velocity. The pressure gradient required by the mantle flow may be manifest in gravity trends across moving plates, but no excessive gravity anomalies are required by the model if the absolute viscosity values conform to those inferred from post-glacial rebound data. A thinner and lower-viscosity layer is favored over a thicker and more viscous layer if both fit glacial rebound evidence. The present model may not be applicable if down to the core the viscosity is as low as about 10²¹ N s m^–2 with a free-slip bottom boundary.     

Ra in the Dead Sea

The²²⁸Ra concentrations of the Dead Sea waters range from 0.13 to 1.48 dpm kg⁻¹, two to three orders of magnitude higher than those of ocean waters and lake waters. However, the²²⁸Ra/²²⁶Ra activity ratios, (0.12–1.29) × 10⁻², are in the range reported for the hydrosphere.The surface waters of the Dead Sea are enriched in²²⁸Ra by a factor of about three over the near-bottom waters. There is a factor of about two spatial variability in the mid-depth Ra concentrations at the two profile stations. The near-bottom²²⁸Ra gradients yield vertical eddy diffusivity coefficient (K) of 2.0 and 0.4 cm² s⁻¹ at profile locations 1 and 2 respectively. These values are comparable to those measured in oceans and lakes.     

Thermal Structure of the Ionian Slab

We propose a thermal model of the subducting Ionian microplate. The slab sinks in an isothermal mantle, and for the boundary conditions we take into account the relation between the maximum depth of seismicity and the thermal parameter L_th of the slab, which is a product of the age of the subducted lithosphere and the vertical component of the convergence rate. The surface heat-flux dataset of the Ionian Sea is reviewed, and a convective geotherm is calculated in its undeformed part for a surface heat flux of 42 mW m^–2, an adiabatic gradient of 0.6 mK m^–1, a mantle kinematic viscosity of 10¹⁷ m² s^–1 and an asthenosphere potential temperature of 1300°C. The calculated temperature-depth distribution compared to the mantle melting temperature indicates the decoupling limit between lithosphere and asthenosphere occurs at a depth of 105 km and a temperature of 1260°C. A 70–km thick mechanical boundary layer is found. By considering that the maximum depth of the seismic events within the slab is 600 km, a L_th of 4725 km is inferred. For a subduction rate equal to the spreading rate, the corresponding assimilation and cooling times of the microplate are about 7 and 90 Myr, respectively. The thermal model assumes that the mantle flow above the slab is parallel and equal to the subducting plate velocity of 6 cm yr^–1, and ignores the heat conduction down the slab dip. The critical temperature, above which the subduced lithosphere cannot sustain the stress necessary to produce seismicity, is determined from the thermal conditions governing the rheology of the plate. The minimum potential temperature at the depth of the deepest earthquake in the slab is 730°C.     

Viscosity of a Teide phonolite in the welding interval

The viscosity of a natural phonolitic composition with variable amounts of H₂O has been experimentally determined. The starting materials were crystal-free phonolitic glasses from Montaña Blanca, situated within the Las Cañadas caldera of Teide. Dry phonolitic melt viscosities were determined using concentric cylinder viscometry in the low viscosity range. The glassy quench products of these runs were then hydrated by high pressure synthesis in a piston–cylinder apparatus to generate a suite of samples with water contents ranging from 0.02 to 3.75 wt%. Samples thus hydrated were quenched rapidly and prepared (cut and polished) for the determination of water contents by infrared spectroscopy before and after experimental viscometry. The viscosities of the melts (dry and hydrated) were determined at 1 bar using a micropenetration technique. Samples were stable under the measurement conditions up to 3.75 wt% H₂O. Homogeneity of water content was confirmed by infrared spectroscopy and total water contents were calculated using absorptivity coefficients for compositions extremely close to that investigated here. The variation of viscosity as a function of water content and temperature can be described in the high viscosity interval of relevance to many welding processes by the non-Arrhenian expression:

(1)

log₁₀ η=−5.900−0.286 ln (H₂O)+(10775.4−394.8(H₂O))/(T−148.7+21.65 ln (H₂O))

whereas the high viscosity range alone is adequately described by the Arrhenian expression

(2)

log₁₀ η=−10.622−0.738 ln (H₂O)+(17114.3−590.4(H₂O))×1/T

where η is the viscosity in Pa s, H₂O is the water content in wt% and T is the temperature in K.These results are particularly useful for the scaling of conditions extant during the welding of phonolitic products of Montaña Blanca. The welding of glassy phonolitic rocks is enhanced by the lower viscosity of these melts with respect to calcalkaline rhyolites. The ratio of viscosities of phonolitic to calcalkaline rhyolitic melts is a complex function of temperature and water content and reaches up to 10^4.5 at 0.1 wt% H₂O and 500°C. Abundant evidence of welding and remobilisation of pyroclastic and spatter products of Teide system volcanism are consistent with these experimental observations.     

Helium accumulation in groundwater, III. Limits on helium transfer across the mantle-crust boundary beneath Australia and the magnitude of mantle degassing

The groundwaters of the Great Artesian Basin (Australia) have been previously shown to be accumulating in-situ production helium for groundwaters ages < 50 kyr and an external helium flux equivalent to whole crustal production for groundwater ages > 100 kyr [1,2]. New helium isotope measurements show that the observed in-situ production helium (³He/⁴He 1.6 × 10⁻⁸) is isotopically distinct from the crustal degassing helium flux (³He/⁴He 6.6 × 10⁻⁸). Furthermore, the crustal degassing helium isotope ratio is marginally in excess of the whole crustal production ratio (³He/⁴He= 3.5 × 10⁻⁸) and the production ratio in a variety of continental rock types. This suggests that the upper limit on volatile transport across the mantle-crust boundary beneath the (relatively) stable and “complacent” Australian continent can be characterized by a “conductive-diffusive” helium/heat flux ratio of 2.6 × 10⁶⁴He atoms mW⁻¹ s⁻¹ which is two orders of magnitude less than the “intrusive-volcanic” ratio of 2.9 × 10⁸⁴He atoms mW⁻¹ s⁻¹ measured at the Galapagos [16]. These results constrain the transcrustal mantle degassing fluxes of⁴He and⁴⁰Ar to be much less than the mid-ocean ridge degassing fluxes; which are much less than the degassing of⁴He and⁴⁰Ar from continental crust. Thus, the degassing of the Earth's interior is dominated by magmatic processes but the dominant fluxes of⁴He and⁴⁰Ar to the atmosphere must come from the continental crust.     

Constraints on the dynamics of mantle plumes from uplift of the Hawaiian Islands

The ∼0.2 mm/yr uplift of Hawaiian islands Lanai and Molokai and Hawaiian swell topography pose important constraints on the structure and dynamics of mantle plumes. We have formulated 3-D models of mantle convection to investigate the effects of plume-plate interactions on surface vertical motions and swell topography. In our models, the controlling parameters are plume radius, excess plume temperature, and upper mantle viscosity. We have found that swell height and swell width constraints limit the radius of the Hawaiian plume to be smaller than 70 km. The additional constraint from the uplift at Lanai requires excess plume temperature to be greater than 400 K. If excess plume temperature is 400 K, models with plume radius between 50 and 70 km and upper mantle viscosity between 10²⁰ and 3×10²⁰ Pa s satisfy all the constraints. Our results indicate that mantle plume in the upper mantle may be significantly hotter than previously suggested. This has important implications for mantle convection and mantle melting. In addition to constraining plume dynamics, our models also provide a mechanism to produce the observed uplift at Lanai and Molokai that has never been satisfactorily explained before.     

Hydrothermal vent waters at 13°N on the East Pacific Rise: isotopic composition and gas concentration

Gas concentrations and isotopic compositions of water have been measured in hydrothermal waters from 13°N on the East Pacific Rise. In the most Mg-depleted samples ( 5 × 10⁻³ moles/kg) the gas concentrations are: 3–4.5 × 10⁻⁵ cm³ STP/kg helium, 0.62–1.24 cm³ STP/kg CH₄, 10.80–16.71 × 10⁻³ moles/kg CO₂. The samples contain large quantities (95–126 cm³/kg) of H₂ and some carbon monoxide (0.26–0.36 cm³/kg) which result from reaction with the titanium sampling bottles. δ¹³C in methane and CO₂ (−16.6 to −19.5 and −4.1 to −5.5 respectively) indicate temperatures between 475 and 550°C, whereas δ¹³C_CO is compatible with formation by reduction of CO₂ on Ti at 350°C close to the sampling temperature.³He/⁴He are very homogeneous at (7.5 ± 0.1)R_A(³He/⁴He = 1.0 × 10⁻⁵) and very similar to already published data as well as CH₄/³He ratios between 1.4 and 2.1 × 10⁶.¹⁸O and D in water show enrichments from 0.39 to 0.69‰ and from 0.62 to 1.49‰ respectively. These values correspond to W/R ratios of 0.4–7. The distinct¹⁸O enrichments indicate that the isotopic composition of the oceans is not completely buffered by the hydrothermal circulations. The³He-enthalpy relationship is discussed in terms of both hydrothermal heat flux and³He mantle flux.     

Mercury's magnetic field: a thermoelectric dynamo?

Permanent magnetism and conventional dynamo theory are possible but problematic explanations for the magnitude of the Mercurian magnetic field. A new model is proposed in which thermoelectric currents driven by temperature differences at a bumpy core-mantle boundary are responsible for the (unobserved) toroidal field, and the helicity of convective motions in a thin outer core (thickness 10² km) induces the observed poloidal field from the toroidal field. The observed field of 3 × 10⁻⁷ T can be reproduced provided the electrical conductivity of Mercury's semiconducting mantle approaches 10³ Ω⁻¹ m⁻¹. This model may be testable by future missions to Mercury because it predicts a more complicated field geometry than conventional dynamo theories. However, it is argued that polar wander may cause the core-mantle topography to migrate so that some aspects of the rotational symmetry may be reflected in the observed field.     

A 100 m laser strainmeter system installed in a 1 km deep tunnel at Kamioka, Gifu, Japan

We have installed a laser strainmeter system in a deep tunnel about 1,000 m below the ground surface at Kamioka, Gifu, Japan. The system consists of three types of independent interferometers: (1) an EW linear strainmeter of the Michelson type with unequal arms, (2) an NS-EW differential strainmeter of the Michelson type with equal arms and (3) a NS absolute strainmeter of the Fabry–Perot type. These are configured in L-shaped vacuum pipes, each of which has a length of 100 m. (1) and (2) are highly sensitive (order of 10⁻¹³ strain) and have wide dynamical range (10⁻¹³–10⁻⁶). Observations with strainmeters (1) and (2) started on June 11, 2003. (3) is a new device for absolute-length measurements of the order of 10⁻⁹ of a long-baseline (100 m) Fabry–Perot cavity by the use of phase-modulated light. This third strainmeter will be ready for operation before the end of 2004. The laser source of strainmeters (1) and (2) is a frequency-doubled YAG laser with a wavelength of 532 nm. The laser frequency is locked onto an iodine absorption line and a stability of 2 × 10⁻¹³ is attained. The light paths of the laser strainmeter system are enclosed in SUS304 stainless steel pipes. The inside pressure is kept to be 10⁻⁴ Pa. Consequently, quantitative measurement of crustal strains of the order of 10⁻¹³ can be attained by employing the laser strainmeter system of (1) and (2) at Kamioka. This resolving power corresponds to that of a superconducting gravimeter. Using the laser strainmeter system, we expect to determine parameters related to fluid core resonance, core modes and core undertone as well as other geodynamic signals such as slow strain changes caused by silent earthquakes or slow earthquakes.     

The polar spirals of Mars may be due to glacier surges deflected by Coriolis forces

All previous accounts of the spiral patterns at the Martian poles emphasize that the north polar spiral is centered about the geographic pole, whereas that of the south polar region is off-set by about 4°. This paper demonstrates that the patterns near both poles are centered on topographic highs rather than the spin poles themselves. This is circumstantial evidence in favour of the relatively unexplored mechanism of radial outflow of viscous rock by gravity spreading.The hypothesis developed here is that the spiral patterns are essentially due to crevasse patterns formed perpendicular to flow lines which are perturbed by Coriolis forces. In order to account for a crevasse pattern that has a form concave to the east the angular deflection of an hypothetical ice flow emanating from the topographic high centered about the geographical north pole, must be about 40° or 0.7 radians in a westward direction at 85°N latitude.The polar cap rock has previously been assumed to consist mainly of either frozen carbondioxide or water ice. Corresponding viscosities (at 190 K) allow for the occurrence of radial outflow or gravity driven tectonics at a maximum rate of 1 cm a⁻¹, but the flow pattern remains unaffected by Coriolis forces.The spiral patterns of the Martian poles can be explained if the flowing mass has an occasional effective kinematic viscosity as low as about 7 × 10⁶ m² s⁻¹, because gravity tectonics will then be deflected by Coriolis forces resulting in appropriately curved flowlines. A tensile fracture pattern, resembling an anticlockwise spiral pattern perpendicular to the clockwise deflected flowlines may subsequently form by local brittle failure.The occasional kinematic viscosity 7 × 10⁶ m² s⁻¹ would cause flow rates of 0.2 m s⁻¹ along the slopes of the topographic highs. This velocity and the corresponding viscosity is tentatively thought to be possible when thermal and pressure runaway occurs in the polar layered deposits. This would mean glacier surges on the Martian poles are two orders of magnitude faster than those hitherto observed on Earth.     

Isotope characteristics of submarine lavas from the Philippine Sea: implications for the origin of arc and basin magmas of the Philippine tectonic plate

Igneous rocks from the Philippine tectonic plate recovered on Deep Sea Drilling Project Legs 31, 58 and 59 have been analyzed for Sr, Nd and Pb isotope ratios. Samples include rocks from the West Philippine Basin, Daito Basin and Benham Rise (40–60 m.y.), the Palau-Kyushu Ridge (29–44 m.y.) and the Parece Vela and Shikoku basins (17–30 m.y.). Samples from the West Philippine, Parece Vela and Shikoku basins are MORB (mid-ocean ridge basalt)-like with ⁸⁷Sr/⁸⁶Sr= 0.7026−0.7032, ¹⁴³Nd/¹⁴⁴Nd= 0.51300−0.51315, and ²⁰⁶Pb/²⁰⁴Pb= 17.8−18.1. Samples from the Daito Basin and Benham Rise are OIB (oceanic island basalt)-like with ⁸⁷Sr/⁸⁶Sr= 0.7038−0.7040, ¹⁴³Nd/¹⁴⁴Nd= 0.51285−0.51291 and ²⁰⁶Pb/²⁰⁴Pb= 18.8−19.2. All of these rocks have elevated ²⁰⁷Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb compared to the Northern Hemisphere Regression Line (NHRL) and have δ²⁰⁷Pb values of 0 to +6 and δ²⁰⁸Pb values of +32 to +65. Lavas from the Palau-Kyushu Ridge, a remnant island arc, have ⁸⁷Sr/⁸⁶Sr= 7032−0.7035, ¹⁴³Nd/¹⁴⁴Nd= 0.51308−0.51310 and ²⁰⁶Pb/²⁰⁴Pb= 18.4−18.5. Unlike the basin magmas erupted before and after them, these lavas plot along the NHRL and have Pb-isotope ratios similar to modern Pacific plate MORB's. This characteristic is shared by other Palau-Kyushu Arc volcanic rocks that have been sampled from submerged and subaerial portions of the Mariana fore-arc.At least four geochemically distinct magma sources are required for these Philippine plate magmas. The basin magmas tap Source 1, a MORB-mantle source that was contaminated by EMI (enriched mantle component 1 [31]) and Source 2, an OIB-like mantle source with some characteristics of EMII (enriched mantle component 2 [31]). The arc lavas are derived from Source 3, a MORB-source or residue mantle including Sr and Pb from the subducted oceanic crust, and Source 4, MORB-source or residue mantle including a component with characteristics of HIMU (mantle component with high U/Pb [31]). These same sources can account for many of the isotopic characteristics of recent Philippine plate arc and basin lavas. The enriched components in these sources which are associated with the DUPAL anomaly were probably introduced into the asthenosphere from the deep mantle when the Philippine plate was located in the Southern Hemisphere 60 m.y.b.p.