Upper mantle seismic wave velocity: Effects of realistic partial melt geometries

We investigate seismic wave velocity reduction resulting from the presence of partial melt in the upper mantle. The amount; of shear and bulk modulus reduction produced by the presence of a connected network of realistically shaped and naturally organized melt inclusions is found using finite element calculations. The geometries of the inclusions are taken directly from laboratory experiments of mantle melting, with finite element meshes constructed to conform to these shapes. The shear and bulk moduli of the composite material are found for both the unrelaxed (isolated inclusions) and relaxed (pressure equalized inclusions) cases by assigning appropriate material properties to the fluid, Modulus reduction from deformation simulations of a solid containing realistically shaped and ellipse-shaped melt inclusions quantify the effect of melt pocket cuspateness and melt pocket organization on seismic velocity reduction. The three-dimensional response is estimated from two-dimensional distributions of the melt phase by determining the mode II and mode III components of elastic modulus reduction separately and summing their effects. In general, cuspate and naturally organized melt inclusions cause greater velocity reduction. It is shown that V-P and V-S reduction per percent partial melt are at least 3.6% and 7.9%, respectively. Even higher values for velocity reduction are possible above 1% melt fraction if melt exists only in tubules below 1% melt fraction. The lower, more conservative values of velocity reduction are similar to 70% greater for V-P and 84% greater for V-S than the analytically determined values for ellipsoidal inclusions. Somewhat greater effects are possible if nonrandom organization of melt occurs on scales greater than our model.