Numerical simulations reproduce field observations showing transient weakening during shear zone formation by diffusional hydrogen influx and H2O inflow

Abstract of the corresponding paper:

Exposures on Holsnøy island (Bergen Arcs, Norway) indicate fluid infiltration through fractures into a dry, metastable granulite, which triggered a kinetically delayed eclogitization, a transient weakening during fluid-rock interaction, and formation of shear zones that widened during shearing. It remains unclear whether the effects of grain boundary-assisted aqueous fluid inflow on the duration of granulite hydration were influenced by a diffusional hydrogen influx accompanying the fluid inflow. To better estimate the fluid infiltration efficiencies and the parameter interdependencies, a 1D numerical model of a viscous shear zone is utilized and validated using measured mineral phase abundance distributions and H2O-contents in nominally anhydrous minerals (NAMs) of the original granulite assemblage to constrain the hydration by aqueous fluid inflow and diffusional hydrogen influx, respectively. Both hydrations are described with a diffusion equation and affect the effective viscosity. Shear zone kinematics are constrained by the observed shear strain and thickness. The model fits the phase abundance and H2O-content profiles if the effective hydrogen diffusivity is approximately one order of magnitude higher than the diffusivity for aqueous fluid inflow. The observed shear zone thickness is reproduced if the viscosity ratio between dry granulite and deforming, re-equilibrating eclogite is ~104 and that between dry granulite and hydrated granulite is ~102. The results suggest shear velocities <10-2 cm/a, hydrogen diffusivities of ~10-13±1 m2/s, and a shearing duration of <10 years. This study successfully links and validates field data to a shear zone model and highlights the importance of hydrogen diffusion for shear zone widening and eclogitization.

The files uploaded here represents the numerical codes and dataset utilized to produce the results presented in the associated paper.
Find a description of the single code files ([name].m) in the ST1.doxc file.