Mixing relations for the “high-pressure” parameters of the equation of state of Ghiorso (2004a) are developed and compositional coefficients are optimized to permit calculation of melt density in portions of the system SiO<sub>2</sub>-TiO<sub>2</sub>-Al<sub>2</sub>O<sub>3</sub>-FeO-MgO-CaO-Na<sub>2</sub>O-K<sub>2</sub>O to pressures in excess of 40 GPa and temperatures up to 2500°C. Four data sets are analyzed and fitted to yield an internally consistent model: (1) density estimates made from measurements of the sinking/floating of reference mineral markers in silicate liquids at known temperatures and pressures, (2) density estimates obtained from shock compression studies on molten liquids, (3) liquid densities inferred from the temperature and pressure dependence of the slopes of mineral fusion curves, and (4) estimates of densities of molten silicate liquids obtained by molecular dynamics simulations. Calibration compositions include chemically complex liquids (komatiite, peridotite and MORB bulk compostions) as well as simple liquids with mineral-like stoichiometry. The model recovers density with an average error of ∼2 percent. The model is limited by <i>not</i> including the effects of volatiles or oxidized iron at elevated pressure.