White Rose University Consortium logo
University of Leeds logo University of Sheffield logo York University logo

Modeling metamorphic fluid flow with reaction-compaction permeability feedbacks

Yardley, B. and Balashov, V.N. (1998) Modeling metamorphic fluid flow with reaction-compaction permeability feedbacks. American Journal of Science, 298. pp. 441-470. ISSN 0022-9599

Available under licence : See the attached licence file.

Download (1355Kb)


Existing models of metasomatic flow do not allow for the effect that reaction has on the flow patterns. Instead, it is assumed that the volatiles produced are negligible in volume compared to those infiltrated and that reaction does not modify permeability. This is clearly unlikely to be true for infiltration-driven decarbonation reactions. The rates of porosity creation by reaction and porosity loss by creep have been calculated for a representative volume of calcite –quartz-wollastonite marble, and it is found that, even for a weak calcite matrix, the rate of porosity generation by reaction is likely to outstrip the collapse of porosity, as long as the system is out of equilibrium. We have applied a self-consistent 1D finite-difference model to the reaction of calcite + quartz to wollastonite in a 10m thick marble, in response to influx of H2O rich fluid, with fixed boundary conditions. The model allows us to evaluate the effect of reaction on the porosity structure and fluid pressure variation across the layer, from which local Darcy fluxes can be evaluated. The progress of reaction that we model is constrained by hydrological considerations, with the requisite parameters recalculated as reaction progresses, assuming creep compaction of rock under the stress difference between lithostatic and fluid pressures. We fnd that the volume of fluid realised by decarbonation, driven by influx of H20, is sufficient to create a back-flow, so that further advancement of the reaction front is only possible as a result of diffusion of water against the Darcy flux. The effect of creep driven by differences between fluid pressure and lithostatic pressure is to reduce the permeability of the layer and especially reduce the secondary porosity developed in the zone at and behind the advancing reaction front. We predict that in a 3D situation, the porous zone of reacted marble becomes a conduit for layer-parallel flow, and the secondary porosity is infilled by calc-silicate minerals due to silica metasomatism.

Item Type: Article
Institution: The University of Leeds
Academic Units: The University of Leeds > Faculty of Environment (Leeds) > School of Earth and Environment (Leeds)
Depositing User: Sherpa Assistant
Date Deposited: 08 Apr 2005
Last Modified: 07 Jun 2014 11:33
Published Version: http://www.geology.yale.edu/journals/Ajs.html/
Status: Published
Refereed: Yes
URI: http://eprints.whiterose.ac.uk/id/eprint/382

Actions (repository staff only: login required)