Simple models of coupled fluid infiltration and redox reactions in the crust

Simple one-dimensional numerical models are presented for coupled advection-hydrodynamic dispersion and kinetically controlled oxidation-reduction reactions in graphite-free porous media containing magnetite coexisting with silicate assemblages. Fluid-solid interactions involving either O-H (O-2-H2O-H-2) or C-O-H (O-2-H-2-O-H-2-CO2-CO-CH4) fluids are considered at similar to 500 degrees C and 5 kbar. The major implications of the modeling are as follows: (1) Regional (km scale) reduction of typical magnetite-bearing rocks originally at fo(2) near NNO may be possible during long-term metamorphic fluid flow if the infiltrating fluids have sufficiently ion fo(2) and sufficiently large concentrations of CH4 and/or H-2. Regional oxidation of such rocks by highly oxidized O-H or C-O-H fluids appears to be difficult to achieve, (2) Nearly identical mineral assemblages and modes may be produced by very different kinetic reaction pathways. The model implies that "equilibrium" assemblages preserved in rocks may not always reflect the true kinetic reaction path that evolved during fluid flow, and highlights the need for quantitative measurements of metamorphic reaction rates. (3) Preservation of sharp lithologic contacts between rocks of very different redox states containing accessory amounts of oxides may be unlikely if fluid-rock interaction times exceed 10(3)-10(4) years. Substantial contact disruption over these times scales is predicted even for oxide-rich rocks if redox contrasts between layers are large. Flow across lithologic contacts malt produce asymmetric patterns of metasomatic mineral zonation that may prove useful for mapping flow directions in metamorphic sequences. (4) For fluid flow in typical T gradients through originally homogeneous rock, significant major element metasomatism (e.g.. K, Na. Ca) may be possible without producing large changes in oxide abundances.