Mineral equilibria calculations in the system K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (KFMASHTO) using thermocalc and its internally consistent thermodynamic dataset constrain the effect of TiO₂ and Fe₂O₃ on greenschist and amphibolite facies mineral equilibria in metapelites. The end‐member data and activity–composition relationships for biotite and chloritoid, calibrated with natural rock data, and activity–composition data for garnet, calibrated using experimental data, provide new constraints on the effects of TiO₂ and Fe₂O₃ on the stability of these minerals. Thermodynamic models for ilmenite–hematite and magnetite–ulvospinel solid solutions accounting for order–disorder in these phases allow the distribution of TiO₂ and Fe₂O₃ between oxide minerals and silicate minerals to be calculated. The calculations indicate that small to moderate amounts of TiO₂ and Fe₂O₃ in typical metapelitic bulk compositions have little effect on silicate mineral equilibria in metapelites at greenschist to amphibolite facies, compared with those calculated in KFMASH. The addition of large amounts of TiO₂ to typical pelitic bulk compositions has little effect on the stability of silicate assemblages; in contrast, rocks rich in Fe₂O₃ develop a markedly different metamorphic succession from that of common Barrovian sequences. In particular, Fe₂O₃‐rich metapelites show a marked reduction in the stability fields of staurolite and garnet to higher pressures, in comparison to those predicted by KFMASH grids.