Estimating the local paleo-fluid flow velocity: New textural method and application to metasomatism

Crystal growth driven by a flowing solution is modeled for a flow with low Reynolds number using a computational dynamic software. Considering equivalent crystallographic faces, the chemical flux is calculated along upstream and downstream faces. Upstream flux is higher compared to downstream and leads to a symmetry breakdown of the crystal shape and develops mirror symmetry parallel to the flow velocity. Moreover the ratio of these two fluxes (upstream/downstream) gives a quantitative relation between the relative crystal growth rate and the flow velocity. Thus, using an inverse method, the flow direction and velocity can be deduced by the study of the variation of the growth band thicknesses of equivalent crystallographic faces.This new method was applied to the formation of metasomatic tourmalinite associated with a leucogranite sill. The approach is complemented by a study of the chemistry of the tourmaline. In the studied case, the application of the new method gives the high fluid velocities in pores (10− 3–10− 4 m/s) during metasomatism. Equivalent Darcy velocities are estimated and discussed accounting for the major role played by the regional deformation. Finally, a two-stage tectono-hydrodynamic model is proposed for the metasomatism. The first stage is genetically linked to the sill injection, and the second is characterized by a wider event with hydrothermal flow passing along the leucogranite sills.