Three different models were developed to simulate the effect of contact metamorphism and fluid-rock interaction on the prograde mineralogical and O isotopic evolution of calcareous argillites from the Notch Peak aureole, Utah. All models assume local mineral-fluid equilibrium, a steady-state temperature profile corresponding to peak metamorphic values, and the thermodynamic data for minerals and fluid of Berman (1988). The preferred model, metamorphism with flow of a time-integrated fluid flux of 2 +/- 0.5.10(4) mol/cm2 in the direction of increasing temperature, successfully reproduces the principal petrologic and isotopic features of the aureole: (1) occurrence and positions (in map view) of diopside-in, tremolite-out, grossular-in, wollastonite-in, and quartz-out isograds; (2) stable coexistence of tremolite + calcite + quartz + diopside over an almost-equal-to 1 km distance between the diopside-in and tremolite-out isograds; (3) variable whole-rock O-18 depletions of almost-equal-to 6-9% adjacent to the contact; and (4) a gradual and irregular increase in delta O-18 with increasing distance from the pluton. One unsuccessful model considered interaction of rock with a tiny, stagnant, fluid-filled porosity with fluid pressure equal to either lithostatic or hydrostatic pressure. Regardless of fluid pressure, the model fails to predict the position of the tremolite-out isograd; the occurrence of grossular, wollastonite, and quartz-out isograds at any position; and O-18 depletions > 0.3%. The other unsuccessful model, flow of magmatic fluid in the direction of decreasing temperature through the aureole, fails to predict the observed range of whole-rock isotopic compositions within 0.5 km of the contact, the observed position of the diopside isograd, and the observed spatial distribution of tremolite + calcite + quartz + diopside. Results demonstrate how isotopic and petrologic data for contact aureoles can be integrated to provide quantitative constraints on the magnitude and geometry of metamorphic fluid flow.
