28.12.2013 Views

sgr ms thesis - University of Maine

sgr ms thesis - University of Maine

sgr ms thesis - University of Maine

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!