As part III of this series, the model is extended to iron oxide-containing melts. All available experimental data in the FeO-Fe2O3-Na2O-K2O-MgO-CaO-MnO-Al2O3-SiO2 system were critically evaluated based on the experimental condition. The variations of FeO and Fe2O3 in the melts were taken into account by using FactSage to calculate the Fe2+/Fe3+ distribution. The molar volume model with unary and binary model parameters can be used to predict the molar volume of the molten oxide of the Li2O-Na2O-K2O-MgO-CaO-MnO-PbO-FeO-Fe2O3-Al2O3-SiO2 system in the entire range of compositions, temperatures, and oxygen partial pressures from Fe saturation to 1 atm pressure.