A practical method is presented to estimate rare earth element (REE) concentrations in a magmatic fluid phase in equilibrium with water-saturated granitic melts based on empirical fluid–melt partition coefficients of REE ( k P R E E ). The values of k P R E E can be calculated from a set of new polynomial equations linking to the Cl molality ( m C l v ) of the magmatic fluid phase associated with granitic melts, which are established via a statistical analysis of the existing experimental dataset. These equations may be applied to the entire pressure range (0.1 to 10.0 kb) within the continental crust. Also, the results indicate that light REEs (LREE) behave differently in magmatic fluids, i.e., either being fluid compatible with higher m C l v or fluid incompatible with lower m C l v values. In contrast, heavy REEs (HREE) are exclusively fluid incompatible, and partition favorably into granitic melts. Consequently, magmatic fluids tend to be rich in LREE relative to HREE, leading to REE fractionation during the evolution of magmatic hydrothermal systems. The maximum k P R E E value for each element is predicted and presented in a REE distribution diagram constrained by the threshold value of m C l v . The REE contents of the granitic melt are approximated by whole-rock analysis, so that REE concentrations in the associated magmatic fluid phase would be estimated from the value of k P R E E given chemical equilibrium. Two examples are provided, which show the use of this method as a REE tracer to fingerprint the source of ore fluids responsible for the Lake George intrusion-related Au–Sb deposit in New Brunswick (Canada), and the Bakircay Cu–Au (–Mo) porphyry systems in northern Turkey.