Titanomagnetite–melt partitioning of Mg, Mn, Al, Ti, Sc, V, Co, Ni, Cu, Zn, Ga, Zr, Nb, Mo, Hf and Ta was investigated experimentally as a function of oxygen fugacity (fO2) and temperature (T) in an andesitic–dacitic bulk-chemical compositional range. In these bulk systems, at constant T, there are strong increases in the titanomagnetite–melt partitioning of the divalent cations (Mg2+, Mn2+, Co2+, Ni2+, Zn2+) and Cu2+/Cu+ with increasing fO2 between 0.2 and 3.7 log units above the fayalite–magnetite–quartz buffer. This is attributed to a coupling between magnetite crystallisation and melt composition. Although melt structure has been invoked to explain the patterns of mineral–melt partitioning of divalent cations, a more rigorous justification of magnetite–melt partitioning can be derived from thermodynamic principles, which accounts for much of the supposed influence ascribed to melt structure. The presence of magnetite-rich spinel in equilibrium with melt over a range of fO2 implies a reciprocal relationship between a(Fe2+O) and a(Fe3+O1.5) in the melt. We show that this relationship accounts for the observed dependence of titanomagnetite–melt partitioning of divalent cations with fO2 in magnetite-rich spinel. As a result of this, titanomagnetite–melt partitioning of divalent cations is indirectly sensitive to changes in fO2 in silicic, but less so in mafic bulk systems

Lattice diffusion coefficients and partition coefficients have been determined for Li, Mg, Al, Sc, Ti, Cr, V, Mn, Co, Ni, Cu, Zn, Ga, Y, Zr, Nb, Mo, In, Lu, Hf, Ta and U in single crystals of natural magnetite as a function of oxygen fugacity (fO 2 ) at 1150 °C and 1 bar by equilibration with a synthetic silicate melt reservoir. Most experiments were run for twelve hours, which was sufficient to generate diffusion profiles from 25 to > 1000 µm in length. The results were checked at one condition with two additional experiments at 66.9 and 161 h. The profiles were analysed using scanning laser-ablation inductivelycoupled-plasma mass-spectrometry. Diffusion coefficients (D) were calculated by fitting data from individual element diffusion profiles to the conventional diffusion equation for one-dimensional diffusion into a semi−infinite slab with constant composition maintained in the melt at the interface. Equilibrium magnetite/melt partition coefficients are given by the ratio of the interface concentrations to those in the melt. Plots of log D as a function of log fO 2 produce V-shaped trends for all the investigated elements, representing two different mechanisms of diffusion that depend on (fO 2 ) −2/3 and (fO 2 ) 2/3 . Diffusion coefficients at a given fO 2 generally increase in the order: Cr < Mo ≈ Ta < V < Ti < Al < Hf ≈ Nb < Sc ≈ Zr ≈ Ga < In < Lu ≈ Y < Ni < U ≈ Zn < Mn ≈ Mg < Co < Li < Cu. Thus, Cu contents of magnetites are most susceptible to diffusive reequilibration, whereas the original content of Cr should be best preserved.

Protracted magma storage in the deep crust is a key stage in the formation of evolved, hydrous arc magmas that can result in explosive volcanism and the formation of economically valuable magmatic-hydrothermal ore deposits. High magmatic water content in the deep crust results in extensive amphibole ± garnet fractionation and the suppression of plagioclase crystallization as recorded by elevated Sr/Y ratios and high Eu (high Eu/Eu*) in the melt. Here, we use a novel approach to track the petrogenesis of arc magmas using apatite trace element chemistry in volcanic formations from the Cenozoic arc of central Chile. These rocks formed in a magmatic cycle that culminated in high-Sr/Y magmatism and porphyry ore deposit formation in the Miocene. We use Sr/Y, Eu/Eu*, and Mg in apatite to track discrete stages of arc magma evolution. We apply fractional crystallization modeling to show that early-crystallizing apatite can inherit a high-Sr/Y and high-Eu/Eu* melt chemistry signature that is predetermined by amphibole-dominated fractional crystallization in the lower crust. Our modeling shows that crystallization of the in situ host-rock mineral assemblage in the shallow crust causes competition for trace elements in the melt that leads to apatite compositions diverging from bulk-magma chemistry. Understanding this decoupling behavior is important for the use of apatite as an indicator of metallogenic fertility in arcs and for interpretation of provenance in detrital studies.