Numerical Modelling of Melting Processes and Induced Diapirism In the Lower Crust

Abstract

SUMMARY The processes of partial melting and magmatic diapirism within the lower crust are evaluated using a numerical underplating model. Fully molten basalt (T= 1200°C) is emplaced at the Moho beneath a solid granite (T= 750OC) in order that a melt front grows into the granite. If diapirism does not occur, this melt front in the granite reaches a minimal depth in the crust before (like in the molten basalt) crystallization takes place. The density contrast between the partially molten granite layer and the overlying solid granite can lead to a Rayleigh-Taylor instability (RTI) which results in diapiric rise of the partially molten granite. Assuming a binary eutectic system for both the granite and the underplating basalt and a temperature- and stress-dependent rheology for the granite, we numerically solve the governing equations and find (a) that diapirism occurs only within a certain but possibly realistic range of parameters, and (b) that if diapirs occur, they do not rise to levels shallower than 15 or perhaps 12km. The growth rate depends on the degree of melting and the thickness of the partially molten layer, as well as the viscosity of the solid and the partially molten granite. From a comparison of the growth rate with the velocity of a Stefan front it is possible to predict whether a melt front will become unstable and result in diapiric ascent or whether a partially molten layer is created, which remains at depth. We carry out such a comparison using our thermodynamically and thermomechanically consistent model of melting and diapirism.