Melt migration modeling in partially molten upper mantle

The objective of this thesis is to investigate the importance of melt migration in shaping major characteristics of geological features associated with the partial melting of the upper mantle, such as sea-floor spreading, continental flood basalts and rifting. The partial melting produces permeable partially molten rocks and a buoyant low viscosity melt. Melt migrates through the partially molten rocks, and transfers mass and heat. Due to its much faster velocity and appreciable buoyancy, melt migration has the potential to modify dynamics of the upwelling partially molten plumes. I develop a 2-D, two-phase flow model and apply it to investigate effects of melt migration on the dynamics and melt generation of upwelling mantle plumes and focusing of melt migration beneath mid-ocean ridges. Melt migration changes distribution of the melt-retention buoyancy force and therefore affects the dynamics of the upwelling plume. This is investigated by modeling a plume with a constant initial melt of 10% where no further melting is considered. Melt migration polarizes melt-retention buoyancy force into high and low melt fraction regions at the top and bottom portions of the plume and therefore results in formation of a more slender and faster upwelling plume. Allowing the plume to melt as it ascends through the upper mantle also produces a slender and faster plume. It is shown that melt produced by decompressional melting of the plume migrates to the upper horizons of the plume, increases the upwelling velocity and thus, the volume of melt generated by the plume. Melt migration produces a plume which lacks the mushroom shape observed for the plume models without melt migration.;Melt migration forms a high melt fraction layer beneath the sloping base of the impermeable oceanic lithosphere. Using realistic conditions of melting, freezing and melt extraction, I examine whether the high melt fraction layer is able to focus melt from a wide partial melting zone to a narrow region beneath the observed neo-volcanic zone. My models consist of three parts; lithosphere, asthenosphere and a melt extraction region. It is shown that melt migrates vertically within the asthenosphere, and forms a high melt fraction layer beneath the sloping base of the impermeable lithosphere. Within the sloping high melt fraction layer, melt migrates laterally towards the ridge. In order to simulate melt migration via crustal fractures and cracks, melt is extracted from a melt extraction region extending to the base of the crust. Performance of the melt focusing mechanism is not significantly sensitive to the size of melt extraction region, melt extraction threshold and spreading rate. In all of the models, about half of the total melt production freezes beneath the cooling base of the lithosphere, and the rest is effectively focused towards the ridge and forms the crust.;To meet the computational demand for a precise tracing of the deforming upwelling plume and including the chemical buoyancy of the partially molten zone in my models, a new numerical method is developed to solve the related pure advection equations. The numerical method is based on Second Moment numerical method of Egan and Mahoney [1972] which is improved to maintain a high numerical accuracy in shear and rotational flow fields. In comparison with previous numerical methods, my numerical method is a cost-effective, non-diffusive and shape preserving method, and it can also be used to trace a deforming body in compressible flow fields.