Two-dimensional patterns of metamorphic fluid flow and isotopic resetting in layered and fractured rocks

One-dimensional advection-dispersion models predict that characteristic delta(18)O vs. distance and delta(18)O vs. delta(13)C profiles should be produced during isothermal metamorphic fluid flow under equilibrium conditions. However, the patterns of isotopic resetting in rocks that have experienced fluid flow are often different from the predictions. Two-dimensional advection-dispersion simulations in systems with simple geometries suggest that such differences may be as a result of fluid channelling and need not indicate disequilibrium, high dispersivities, or polythermal flow. The patterns of isotopic resetting are a function of: (1) the permeability contrast between more permeable layers ('channels') and less permeable layers ('matrix'); (2) the width and spacing of the channels; (3) the width and spacing of discrete fractures; and (4) the orientation of the pressure gradient with respect to layering. In fractured systems, the efficiency of isotopic transport depends on the fracture aperture and the permeability of the surrounding rock. Resetting initially occurs along and immediately adjacent to the fractures, but with time isotopic resetting because of flow through the rock as a whole increases in importance. Application of the one-dimensional advection-dispersion equations to metamorphic fluid flow systems may yield incorrect estimates of fluid fluxes, intrinsic permeabilities, dispersivities, and permeability contrasts unless fluid flow occurred through zones of high permeability that were separated by relatively impermeable layers.