Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
to make a “hockey puck” style geometry (figure 4.3). The gabbro-diorite unit is<br />
added to the base and is given the same shape and thickness <strong>of</strong> 1km. Model<br />
boundaries must be placed at a substantial distance from the region <strong>of</strong> interest<br />
so as to not interfere with the heat transfer model. Therefore, the outer model<br />
boundaries lie 100 km away from the pluton center and are thermally insulated.<br />
The interface between the intrusion and the country rock is open and allows for<br />
free conductive heat transfer. A tetragonal mesh is used, with a fine boundary<br />
mesh located at the intrusive contact for increased resolution <strong>of</strong> the thermal<br />
evolution adjacent to the heat source.<br />
Initial wall rock temperature is 150°C at an emplacement depth <strong>of</strong> 5km<br />
based on a typical upper crustal geotherm <strong>of</strong> 30°C/km. I am concerned with the<br />
lateral evolution <strong>of</strong> the thermal gradient, so there is no need to calculate the local<br />
geothermal gradient for this model. The intrusion temperature was approximately<br />
900°C, higher than average for a granitic magma (Wiebe, 1997). The gabbro<br />
diorite sheet is given an initial temperature <strong>of</strong> 1200°C. The thermal diffusivity ( )<br />
<strong>of</strong> Bar Harbor Formation is set at 1.30E-6 m 2 /s calculated from<br />
(4.1)<br />
Where is thermal conductivity, is density, and is specific heat. Thermal<br />
diffusivities for Cadillac Mountain Granite and gabbro are given values <strong>of</strong>1.34E-6<br />
m 2 /s and 1.01E-6 m 2 /s, respectively. A series <strong>of</strong> 300 time steps are taken from<br />
time = 0s to time = 1E14s, producing solutions for the evolving thermal gradient<br />
nearly up to a point <strong>of</strong> thermal equilibrium (Figure 4.4).<br />
36