Résumé
Mantle rock studies provide evidence of interaction with upwelling magmas. In orogenic lherzolites, one of the most conspicuous effects of these interactions is the development of harzburgite and dunite bands. Recent studies have suggested that these bands were formed at the expense of the host lherzolites by melt-rock reactions associated with magma percolation. In order to better understand the geochemical effects associated with percolation-reaction processes, we propose a numerical model of melt infiltration that takes into account modal variations in time and space resulting from melt-rock reactions. Melt volume variations are considered by means of porosity variations, and a local equilibrium approach is used for trace element exchange between melt and minerals. The transport of trace elements by the interstitial melt is described by a mass balance equation while the modal variations are constrained by the mineralogical trends observed in refractory peridotites massifs. The model is applied to REE, Cr and Ni in percolated peridotites affected by an olivine-forming reaction, with the aim of reproducing the evolution of these elements in refractory peridotites from the Ronda massif. Our modelling can explain the negative correlation between the LREE/HREE ratio and the HREE content and between Cr and Ni in the Ronda refractory peridotites. Our results validate the hypothesis that, in the Ronda, the bands of refractory peridotites represent porous-flow channels formed by olivine-forming melt-rock reaction, at increasing melt volume. Because similar geochemical features are observed in ophiolitic peridotites and in mantle xenoliths, it is likely that melt-rock reactions associated with magmatic infiltration are widespread and represent important mantle processes.