Computer modeling of granite magma diapirism in the Earth’s crust

A new point of view describing processes of partial melting and development of gravitational instability in a thickening crust with increased thickness of the granite layer is suggested. Numeral experiments support the following main conclusions. The critical volume of partially melted material should be formed for the beginning of flotation in a gravitational field. Due to model estimations, the height of the melting area in the granite crust should be not less than 6–7 km. A mushroom-shaped form of the floating body was observed in all models regardless of the thermal source size (fixed or variable width): the high temperature channel (magma leader) and head body of the diapir are formed. The height of diapir floating depends on rheological features of the surrounding crust: 10 times increase in the yield strength (from 1 to 10 MPa) while temperature decrease confines the possible level of rising to a depth of 15–16 km. An elevation of about 750 m is formed in the day surface relief above the axis part of the diapir.