The reliability of eight Fe-Mg exchange geothermobarometers for garnet-bearing peridotites, pyroxenites and eclogites has been examined using a database comprised of more than 300 published peridotite, pyroxenite and eclogite experiments conducted from 10 to 70 kbar and 850 to > 1650 degrees C. We have tested Fe-Mg exchange geothermometers suitable for a range of mantle lithologies, including websterite, harzburgite, wehrlite and eclogite. All geothermometers maintained an average difference in experimental and calculated temperature (T) (Delta T = Texp-Tcalc) of less than +/- 50 degrees C with a standard deviation of Delta T between +/- 50 to 150 degrees C. Most geothermometers performed well across a narrow range in ln Kd(Fe-Mg)(A-B) (where Kd(Fe-Mg)(A-B) =(Fe(A)xMg(B))/(Fe(B)xMg(A))), however, systematic overestimation and underestimation of T were observed outside of the optimal range of lnKd(Fe-Mg)(A-B). Increases in experimental pressure (P) adversely affected several geothermometers, particularly those calibrated empirically using natural samples. All previously published calibrations of the garnet-clinopyroxene geothermometer were unable to reliably reproduce the experimental T for both peridotite and eclogite experimental compositions, which hinders their confident application to natural datasets. To improve the state of mantle geothermobarometry we have used our experimental database to recalibrate the (1) garnet-clinopyroxene Fe-Mg exchange geothermometer, and (2) garnet-orthopyroxene Fe-Mg exchange geothermometer. Each geothermometer has been recalibrated across an extended P, T, and compositional range. The inclusion of eclogitic experiments in the calibration for the garnet-clinopyroxene geothermometer permits application to both eclogitic and peridotitic/pyroxenitic assemblages equilibrated under a wide range of PT conditions in the upper mantle. Using multiple linear regression to solve for lnKd, we found the following expressions best reproduced the experimental T (degrees C) of our dataset: T-Fe-Mg(grt-cpx)(degrees C) = 3356.34/((-0.008 x P (kbar)) + (0.259 x X-Ca(grt)) + (0.914 x X-Mg(grt)) + (-0.159 x Jd(cpx)) + (ln(Kd(Fe-Mg)(grt-cpx)) + 1.265)) - 273 T-Fe-Mg(grt-opx)(degrees C) = 1851.85/((-0.007 x P (kbar)) + (-1.83 x X-Ca(grt)) + (ln(Kd(Fe-Mg)(grt-cpx)) + 1.08)) - 273 where, X-Ca(grt) = Ca/(Ca+Fe+Mg), Kd(Fe-Mg)(grt-opx) = (Fe-grt xMg(opx))/(Fe(opx)xMg(grt)), X-Mg(grt) = Mg/(Ca+Fe+Mg), Jd(cpx) = Na - Cr - 2 x Ti Kd(Fe-Mg)(grt-cpx) = (Fe-grt xMg(cpx))/(Fe(cpx)xMg(grt)), with all elements calculated on the basis of 12 oxygen anions in garnet and 6 oxygen anions in clino- and orthopyroxene. Fe2+ = total Fe. Our updated calibrations resolve several issues with earlier calibrations, including a poor performance at elevated P and compositional limitations. An improvement in precision and accuracy has been demonstrated through application to the experimental calibration dataset, a second independent set of published experimental data, and to natural peridotites, pyroxenites and eclogites from on and off craton settings. Iterative PT estimates on natural datasets calculated using our updated calibrations compare well with estimates from widely used calibrations such as the Taylor (1998) two-pyroxene solvus geothermometer. We anticipate that this contribution will provide an important reference for the reliability of mantle geothermometers and that our updated calibrations will be used in future studies on peridotite, pyroxenite and eclogite inclusions in diamond and mantle-derived xenoliths.
