skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Contact metamorphism surrounding the Alta stock: Finite element model simulation of heat- and {sup 18}O/{sup 16}O mass-transport during prograde metamorphism

Abstract

Results of calcite-dolomite geothermometry and oxygen isotope studies of marbles in the southern portion of the contact aureole surrounding the Alta stock (Utah) provide evidence for extensive hydrothermal metamorphism in this part of the aureole. Simulation of these two independent data sets with two-dimensional, finite element fluid flow and heat transport models constrains the pattern of fluid flow, minimum permeability, and the permeability structure in this part of the aureole. Model results demonstrate that intrusion of the stock into a homogeneous, isotropic permeability medium yields peak metamorphic temperatures significantly lower than those measured in the marbles and significant {sup 18}O depletions both above and below the Alta-Grizzly thrust system. The latter contradicts the observations in the south aureole that {sup 18}O depletions in the marbles are restricted to marbles below the Alta-Grizzly thrust; dolomitic marbles above the thrust retain original sedimentary values up to the intrusive contact. Models with horizontal permeability barriers above the Alta-Grizzly thrust and extending over the top of the Alta stock are capable of reproducing the observed thermal and {delta}{sup 18}O profiles in the southern aureole. The presence of such horizontal barriers reduces the predominantly vertical fluid flow and heat transfer that would occur in amore » homogeneous and isotropic permeability medium, forcing fluid flow and heat transfer laterally away from the upper flanks of the stock. Such horizontal flow patterns are necessary to produce significant {sup 18}O depletion above the thrust, and to provide the necessary lateral heat transfer to duplicate the observed temperature profile. Best fit model results to the observed thermal and {delta}{sup 18}O profiles provide several new insights into the dynamics of fluid circulation and hydrogeologic characteristics of the southern Alta aureole during prograde metamorphism.« less

Authors:
; ;  [1]
  1. Univ. of Utah, Salt Lake City, UT (United States)
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
530654
Resource Type:
Journal Article
Resource Relation:
Journal Name: American Journal of Science; Journal Volume: 297; Journal Issue: 1; Other Information: PBD: Jan 1997
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; METAMORPHISM; MATHEMATICAL MODELS; MARBLE; UTAH; GEOTHERMOMETRY

Citation Formats

Cook, S.J., Bowman, J.R., and Forster, C.B.. Contact metamorphism surrounding the Alta stock: Finite element model simulation of heat- and {sup 18}O/{sup 16}O mass-transport during prograde metamorphism. United States: N. p., 1997. Web. doi:10.2475/ajs.297.1.1.
Cook, S.J., Bowman, J.R., & Forster, C.B.. Contact metamorphism surrounding the Alta stock: Finite element model simulation of heat- and {sup 18}O/{sup 16}O mass-transport during prograde metamorphism. United States. doi:10.2475/ajs.297.1.1.
Cook, S.J., Bowman, J.R., and Forster, C.B.. 1997. "Contact metamorphism surrounding the Alta stock: Finite element model simulation of heat- and {sup 18}O/{sup 16}O mass-transport during prograde metamorphism". United States. doi:10.2475/ajs.297.1.1.
@article{osti_530654,
title = {Contact metamorphism surrounding the Alta stock: Finite element model simulation of heat- and {sup 18}O/{sup 16}O mass-transport during prograde metamorphism},
author = {Cook, S.J. and Bowman, J.R. and Forster, C.B.},
abstractNote = {Results of calcite-dolomite geothermometry and oxygen isotope studies of marbles in the southern portion of the contact aureole surrounding the Alta stock (Utah) provide evidence for extensive hydrothermal metamorphism in this part of the aureole. Simulation of these two independent data sets with two-dimensional, finite element fluid flow and heat transport models constrains the pattern of fluid flow, minimum permeability, and the permeability structure in this part of the aureole. Model results demonstrate that intrusion of the stock into a homogeneous, isotropic permeability medium yields peak metamorphic temperatures significantly lower than those measured in the marbles and significant {sup 18}O depletions both above and below the Alta-Grizzly thrust system. The latter contradicts the observations in the south aureole that {sup 18}O depletions in the marbles are restricted to marbles below the Alta-Grizzly thrust; dolomitic marbles above the thrust retain original sedimentary values up to the intrusive contact. Models with horizontal permeability barriers above the Alta-Grizzly thrust and extending over the top of the Alta stock are capable of reproducing the observed thermal and {delta}{sup 18}O profiles in the southern aureole. The presence of such horizontal barriers reduces the predominantly vertical fluid flow and heat transfer that would occur in a homogeneous and isotropic permeability medium, forcing fluid flow and heat transfer laterally away from the upper flanks of the stock. Such horizontal flow patterns are necessary to produce significant {sup 18}O depletion above the thrust, and to provide the necessary lateral heat transfer to duplicate the observed temperature profile. Best fit model results to the observed thermal and {delta}{sup 18}O profiles provide several new insights into the dynamics of fluid circulation and hydrogeologic characteristics of the southern Alta aureole during prograde metamorphism.},
doi = {10.2475/ajs.297.1.1},
journal = {American Journal of Science},
number = 1,
volume = 297,
place = {United States},
year = 1997,
month = 1
}
  • The study area extends for 80 km through regionally metamorphosed rocks in the metamorphic high of southern New England, in which Acadian (Devonian) metamorphism ranged from sillimanite-muscovite grade through sillimanite-garnet-cordierite-orthoclase grade. Over this metamorphic interval, amphibolites containing the assemblages plagioclase + hornblende {plus minus} quartz {plus minus} cummingtonite {plus minus} orthoamphibole {plus minus} garnet were transformed by a series of reactions into pyroxene granulites with assemblages of plagioclase + orthopyroxene + augite {plus minus} quartz {plus minus} hornblende {plus minus} garnet. Three particularly important reactions were examined in detail. These reactions, along with others, define an assemblage series in quartz-bearingmore » rocks in the transition from amphibolite facies to granulite facies. This assemblage series is based on the coexistence of diagnostic minerals in common mafic rocks with increasing grade. This transition is estimated to have occurred over a temperature interval of 580-730C.« less
  • An accelerated Finite Element Contact Block Reduction (FECBR) approach is presented for computational analysis of ballistic transport in nanoscale electronic devices with arbitrary geometry and unstructured mesh. Finite element formulation is developed for the theoretical CBR/Poisson model. The FECBR approach is accelerated through eigen-pair reduction, lead mode space projection, and component mode synthesis techniques. The accelerated FECBR is applied to perform quantum mechanical ballistic transport analysis of a DG-MOSFET with taper-shaped extensions and a DG-MOSFET with Si/SiO{sub 2} interface roughness. The computed electrical transport properties of the devices obtained from the accelerated FECBR approach and associated computational cost as amore » function of system degrees of freedom are compared with those obtained from the original CBR and direct inversion methods. The performance of the accelerated FECBR in both its accuracy and efficiency is demonstrated.« less
  • The finite element code FEHMN, developed by scientists at Los Alamos National Laboratory (LANL), is a three-dimensional finite element heat and mass transport simulator that can handle complex stratigraphy and nonlinear processes such as vadose zone flow, heat flow and solute transport. Scientists at LANL have been developing hydrologic flow and transport models of the Yucca Mountain site using FEHMN. Previous FEHMN simulations have used an equivalent Kd model to model solute transport. In this thesis, FEHMN is modified making it possible to simulate the transport of a species with a rigorous chemical model. Including the rigorous chemical equations intomore » FEHMN simulations should provide for more representative transport models for highly reactive chemical species. A fully kinetic formulation is chosen for the FEHMN reactive transport model. Several methods are available to computationally implement a fully kinetic formulation. Different numerical algorithms are investigated in order to optimize computational efficiency and memory requirements of the reactive transport model. The best algorithm of those investigated is then incorporated into FEHMN. The algorithm chosen requires for the user to place strongly coupled species into groups which are then solved for simultaneously using FEHMN. The complete reactive transport model is verified over a wide variety of problems and is shown to be working properly. The new chemical capabilities of FEHMN are illustrated by using Los Alamos National Laboratory`s site scale model of Yucca Mountain to model two-dimensional, vadose zone {sup 14}C transport. The simulations demonstrate that gas flow and carbonate chemistry can significantly affect {sup 14}C transport at Yucca Mountain. The simulations also prove that the new capabilities of FEHMN can be used to refine and buttress already existing Yucca Mountain radionuclide transport studies.« less
  • The finite element code FEHMN is a three-dimensional finite element heat and mass transport simulator that can handle complex stratigraphy and nonlinear processes such as vadose zone flow, heat flow and solute transport. Scientists at LANL have been developed hydrologic flow and transport models of the Yucca Mountain site using FEHMN. Previous FEHMN simulations have used an equivalent K{sub d} model to model solute transport. In this thesis, FEHMN is modified making it possible to simulate the transport of a species with a rigorous chemical model. Including the rigorous chemical equations into FEHMN simulations should provide for more representative transportmore » models for highly reactive chemical species. A fully kinetic formulation is chosen for the FEHMN reactive transport model. Several methods are available to computationally implement a fully kinetic formulation. Different numerical algorithms are investigated in order to optimize computational efficiency and memory requirements of the reactive transport model. The best algorithm of those investigated is then incorporated into FEHMN. The algorithm chosen requires for the user to place strongly coupled species into groups which are then solved for simultaneously using FEHMN. The complete reactive transport model is verified over a wide variety of problems and is shown to be working properly. The simulations demonstrate that gas flow and carbonate chemistry can significantly affect {sup 14}C transport at Yucca Mountain. The simulations also provide that the new capabilities of FEHMN can be used to refine and buttress already existing Yucca Mountain radionuclide transport studies.« less
  • The aim of this research was to study the laser alloying process of iron with chromium. In the paper, a multiphase model of mass and heat transfer for the laser alloying is presented. Laser melting of the chromium layer and the substrate was performed using a continuous laser source operated with a TEM{sub 10} mode, with constant beam diameter ({phi}), scanning velocity (V) and varied output beam power. The partial differential equations of the conservation of mass, momentum and energy in the melted pool for multiphase system were solved. The distribution of chromium in iron after laser alloying was obtainedmore » by including the Volume of Fluid algorithm in the model. The results of the computations were compared with the experimental evaluation of the microstructure and the chromium concentration, which were based on scanning electron microscopy and x-ray microanalysis (Energy Dispersive Spectroscopy) of the laser alloyed layers. The comparison of computational calculations and experimental results is presented and a good accuracy of the proposed model is shown.« less