<p>&#160; &#160; &#160; &#160;Studying magma reservoir processes is one of the keys to understand the evolution of igneous systems. One of the main processes, magma differentiation, governs the thermal evolution and chemical composition of the melt-crystal assemblage (magma or mush depending on the relative proportions), and therefore exerts a first order control over its physical properties (density, viscosity), and thus on reservoir dynamics. Various approaches have been implemented to model differentiation in an attempt to benchmark all the involved variables like initial and phase compositions, temperature, pressure, and oxygen fugacity (C<sub>0</sub>, X, T, P, <em>f</em>O<sub>2</sub>). Those approaches are among others mass balance calculations considering partition coefficients (D) values, experimental studies, thermodynamic models or a combination of those. In any of those cases, the evolution of trace elements is governed by the value of the D that is known to be dependent on (P, T, X, <em>f</em>O<sub>2</sub>). However, most of the present-day studies still use fixed values of D to provide first order estimates.</p><p>&#160; &#160; &#160; Here, we present an approach combining thermodynamic modeling (relying on Rhyolite-MELTS, Gualda et al., 2012), that integrates X-T-P-<em>f</em>O<sub>2</sub>-dependent D for Rare Earth Elements (REE). We applied this new approach to a MORB system, with olivine, clinopyroxene and plagioclase as main mineral phases, and compared results to more classical approaches. D are derived from the models of Sun & Liang (2012, 2013, 2014) and Sun et al. (2017). The resulting model highlights that T & X effects on the D values can add or counterbalance each other depending on the mineral considered. In any cases our results emphasize the gain of using thermodynamic models along with both T- & X-dependent D values to properly model the evolution of igneous systems. Relying on our results, and on the corresponding thermodynamics constraints, we were also able to provide D for any mineral composition crystallized from this MORB system. Results bring to light that an error of ~1 order of magnitude of the D<sup>mineral-melt</sup> value could be introduced when considering a fixed value of D.</p><p>&#160;</p><p>Gualda et al. (2012) <em>Journal of Petrology</em>, 53-5, 875-890; Sun & Liang (2012) <em>Contrib Mineral Petrol</em> 163-5: 807-823; (2013) <em>Chem Geol</em> 358: 23-36; (2014) <em>Chem Geol</em> 372: 80-91; Sun et al. (2017) <em>Geochim Cosmochim Acta</em> 206: 273-295.</p><p>&#160;</p>