Modelling Subduction Zone Magmatism Due to Hydraulic Fracture
Abstract:
To take this further we have developed a thermal model of a subduction zone. The model uses a finite difference, marker-in-cell method to solve the heat equation (Gerya, 2010). The velocity field was prescribed using the analytical expression of cornerflow (Batchelor, 1967). The markers contained within the fixed grid are used to track the different compositions and their properties. The subduction zone thermal model was benchmarked (Van Keken, 2008). We used the hydrous melting parameterization of Katz et.al., (2003) to calculate the degree of melting caused by the addition of water to the wedge. We investigate models where the hydrofractures, with properties constrained by estimated water fluxes, have random end points. The model predicts degree of melting, magma productivity, temperature of the melt and water content in the melt for different initial water fluxes. Future models will also include the buoyancy effect of the melt and residue.
Batchelor, Cambridge UP, 1967.
Davies, Nature, 398: 142-145, 1999.
Gerya, Cambridge UP, 2010.
Katz, Geochem. Geophys. Geosy, 4(9), 2003
Van Keken et.al. Phys. Earth. Planet. In., 171:187-197, 2008.
