HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Full Text Full Text (PDF) Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Wing, B. A. Articles by Ferry, J. M. Search for Related Content Articles by Wing, B. A. Articles by Ferry, J. M. GeoRef GeoRef Citation

Magnitude and geometry of reactive fluid flow from direct inversion of spatial patterns of geochemical alteration

^* Department of Earth and Planetary Sciences, McGill University, Montreal, Quebec, Canada H3A 2A7; wing{at}eps.mcgill.ca
^** Department of Earth and Planetary Sciences, Johns Hopkins University, Baltimore, Maryland 21218 USA

Time-integrated equations governing tracer mass-balance provide a framework for interpreting spatial patterns of mineralogical, elemental, and isotopic alteration in rocks caused by reactive fluid flow, and linear inverse theory allows quantitative estimates of the magnitude and three-dimensional (3D) geometry of fluid flow. We demonstrate the inverse technique with a simple model system that involves fluid advection and fluid-rock oxygen isotope exchange only. Analytic forward models of reactive transport enable patterns of isotope alteration to be calculated exactly, and inversions of the resulting spatial patterns of synthetic isotope composition demonstrate the effectiveness of our approach.

For a field application of the inverse method, we consider regional metamorphic rocks from the Waits River Formation, southeastern Vermont, that expose a mineralogical and stable isotopic record of reactive fluid flow over an area of 120 km². Spatial patterns of prograde changes in whole-rock CO₂, ¹⁸O, and ¹³C allow a direct inversion for the magnitude and 3D geometry of reactive flow at the time of final prograde mineral reaction. Consistency tests suggest that the inverse estimate is primarily controlled by patterns of CO₂ changes, and resolution tests indicate that the smallest features in the inverse flux field that can be reliably detected have an east-west dimension of 5 km and a north-south dimension of 14 km. The regionally-averaged time-integrated fluid flux vector returned by the inversion has a magnitude of 2.2x10⁴ mol fluid/cm² rock, trends 255° in the horizontal plane, and is directed upward at 32° from the horizontal. The geometry of terrain-scale reactive flow was largely controlled by regional structure. Reactive fluids that entered the study area appear to have been tectonically-driven horizontally from the east at a depth 0.5 to 2.0 km beneath the present level of exposure.

This article has been cited by other articles:

				S. C. Penniston-Dorland and J. M. Ferry Element mobility and scale of mass transport in the formation of quartz veins during regional metamorphism of the Waits River Formation, east-central Vermont American Mineralogist, January 1, 2008; 93(1): 7 - 21. [Abstract] [Full Text] [PDF]