ebook img

Modeling magmatic accumulations in the upper crust: Metamorphic implications for the country rock PDF

54 Pages·2016·4.35 MB·English
Save to my drive
Quick download
Download
Most books are stored in the elastic cloud where traffic is expensive. For this reason, we have a limit on daily download.

Preview Modeling magmatic accumulations in the upper crust: Metamorphic implications for the country rock

(cid:0)(cid:2)(cid:2)(cid:3)(cid:4)(cid:5)(cid:3)(cid:6) (cid:7)(cid:8)(cid:9)(cid:10)(cid:11)(cid:2)(cid:12)(cid:13)(cid:4)(cid:5) Modelingmagmaticaccumulationsintheuppercrust: Metamorphicimplica- tionsforthecountryrock Madison M. Douglas, Adelina Geyer, Antonio M. A´lvarez-Valero, Joan Mart´ı PII: S0377-0273(16)30018-X DOI: doi: 10.1016/j.jvolgeores.2016.03.008 Reference: VOLGEO5788 Toappearin: JournalofVolcanologyandGeothermalResearch Receiveddate: 23October2015 Reviseddate: 29February2016 Accepteddate: 8March2016 Pleasecitethisarticleas: Douglas,MadisonM.,Geyer,Adelina,A´lvarez-Valero,Antonio M., Mart´ı, Joan, Modeling magmatic accumulations in the upper crust: Metamorphic implications for the country rock, JournalofVolcanologyandGeothermalResearch (2016), doi: 10.1016/j.jvolgeores.2016.03.008 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. ACCEPTED MANUSCRIPT T P I R C S U N A M D E Corresponding author: Adelina Geyer, Group of Volcanology, SIMGEO (UB-CSIC), T Institute of Earth Sciences PJaume Almera, ICTJA-CSIC, Lluis Solé i Sabaris s/n, 08028 E Barcelona, Spain. ([email protected]) C C A 1 ACCEPTED MANUSCRIPT T P I R C S U N A M D E T P E C C A 2 ACCEPTED MANUSCRIPT T P I R C S U N A M D E T P E C C A The extent and degree of the generated aureole depend on many parameters, including: temperature, composition and thermal properties (i.e. thermal conductivity, specific heat capacity) of the magma and country rock; the degree of crystallization of the magma; the size and geometry of the reservoir; the latent heat of fusion for the crystallizing magma; and the amount of heat energy consumed by endothermic metamorphic reactions taking place in the aureole (Johnson et 3 ACCEPTED MANUSCRIPT al., 2011;Peacock and Spear, 2013). Additionally, at shallow crustal levels, hydrothermal processes may play an important role in incorporating cations and H O into the system, 2 as well as redistributing the heat released from the cooling magTma reservoir that can P affect contact metamorphism equilibria (e.g. Peacock and Spear, 2013). I R C S U N A ) M D E , 2013 T P E C C A 4 ACCEPTED MANUSCRIPT T P I R C S U N A M D E T P E C C A 5 ACCEPTED MANUSCRIPT T P I R Heat transfer models C S The temperature distribution within the mUagma chamber and the host rock is N calculated using FEM, by solving the pure conductive heat transfer equation, and A assuming the effects of viscous heating and pressure-volume work are negligible (see M Rodríguez et al., 2015 for further details): D [1] E T P where the equation parameters refer to density (), specific heat capacity at constant E pressure (Cp), tempeCrature (T), time (t), thermal conductivity (k) and Q denotes heat C sources other than viscous heating (see Table 1 for details of the thermal and physical A parameters). 6 ACCEPTED MANUSCRIPT T Symbol Value Variable P SI Unit A 109 Dorn parameter Pa s D I Cp Specific heat capacity at cRonstant pressure J/(kg K) Cp Specific heat capacity of crust4 J/(kg K) crust Supp. Mat. Table S8 C Cp Supp. Mat. Tables S4-6 Specific heat capacity of magma (melt phase) J/(kg K) melt S Cp Supp. Mat. Tables S4-6 Specific heat capacity of magma (solid phase) J/(kg K) solid Supp. Mat. Table S7 MagUma chamber depth km d 0.25-2 MaNgma chamber aspect ratio - e E 120,000 Activation energy J/mol A A 1.061-9.141 Magma chamber height km h M magma Thermal diffusivity of magma 4 m2/s crust Supp. Mat. Table S8 D Thermal diffusivity of crust4 m2/s E Thermal conductivity W/(m K) k  · Cp ·  L 400 T Latent heat of crystallization1 kJ/kg  Supp. Mat. Tables S1-3 Melt fraction3 - P Q 0.03 Subcrustal heat flow1 W/m2 mantle Q 0.0E7 Surface heat flow1 W/m2 surface crust C2700 Crust density kg/m3  2344-2694 Magma density3 kg/m3 magma ϕ C 1-  Solid fraction3 - R A 8.314 J/(mol K) Molar gas constant thermal_grad 25-30 Crustal thermal gradient K/km T_liquidus Supp. Mat. Tables S1-3 Magma liquidus temperature3 K T_solidus Supp. Mat. Tables S1-3 Magma solidus temperature3 K T_surface 298.15 Surface temperature K T Initial magma temperature K 0_magma T Initial crustal temperature K 0_crust P 1-3106 Magma chamber internal pressure at h/2 Pa  Eq. 3 Country rock viscosity Pa s 4.243-19.696 Magma chamber width w km 1, Bea (2010); 2, Whittington et al. (2009); 3, use Rhyolite MELTS code by Gualda et al. (2012); 4 Rodríguez et al. (2015) Physical and thermal parameters 7 ACCEPTED MANUSCRIPT We consider the modeled magma chamber to be emplaced immediately before it begins cooling. The completely liquefied magma reservoir (Gelman et al., 2013) is T P assumed to intrude in a single episode without further replenishment (Gutiérrez and I Parada, 2010). The geometric modeling, mesh discretizationR and numerical computations C are carried out with COMSOL Multiphysics v4.4 software package S (http://www.comsol.com). To simulate the solidifying magma, we use the heat transfer U with phase change module, which allows uNs to solve Equation 1 after setting the A properties of a phase-change material (from liquid to solid) according to the Apparent M Heat Capacity formulation (AHC) (see Rodríguez et al., 2015 and references therein). D The AHC formulation is considered to be the best representation of a naturally occurring E wide phase-change temperature interval, as occurs during magma cooling (Rodríguez et T al., 2015). The latent heaPt of crystallization is accounted for by increasing the heat E capacity of the material within the phase change temperature range. For further details C on the methodology, the reader is referred to the work by Rodríguez et al. (2015). A C comprehensive report detailing one of the models, which includes all settings within A COMSOL (e.g. model properties; physics settings; geometry; mesh), as well as an example of one of the models, are available from the corresponding author upon request (COMSOL Multiphysics commercial software V4.4 or greater and heat transfer with phase change module are required). The performed FE models are axisymmetric and constructed over a cylindrical coordinate system with positive z values related to altitudes above sea level (Fig. 1a). The magma chamber geometry is oblate, spherical or prolate in shape with height h and width w (Fig. 1a). The computational domain corresponds to a section of the crust with a 40 km 8 ACCEPTED MANUSCRIPT radius stretching to a depth of 40 km below sea level (Fig. 1a). The FE mesh consists of between 200,000 and 400,000 linear triangular elements up to 250 m in size farther from the magma reservoir, and 2 m in size near the edge of the magmTa chamber (Fig. 1b). P Time steps taken to perform the calculations range from less than a year at the beginning I R of the cooling processes, up to 300 years to speed calculations once reached the C metamorphic peak. We have tested that the chosen time steps do not affect the results S obtained. U N The selected starting magmatic compositions for the numerical simulations A correspond to three representative types that range from alkaline basalt to rhyolite (Table M 2). The former, typical of intraplate magmatism, is taken from Polat (2009). The D intermediate magma is of andesitic composition and from a subduction zone taken from E Yogodzinski et al. (1994). The rhyolitic composition is from an archetypal case of pelitic T P crustal melting taken from Sensarna et al. (2004). Depending on the model, pressure at E the magma chamber center ranges from 1 to 3 kbar (Supplementary Material, Table S1). C C A 9

Description:
implications for the country rock, Journal of Volcanology and Geothermal This is a PDF file of an unedited manuscript that has been accepted for
See more

The list of books you might like

Most books are stored in the elastic cloud where traffic is expensive. For this reason, we have a limit on daily download.