Three different ratios between surface sites and added Mg2+ were used in the sorption experiments viz. 0.5, 1, and 2 Mg2+ site?1. Surface complexation modeling of the Mg2+ sorption onto maghemite and magnetite was restricted to pH conditions where the interference from Mg(OH)2(s) precipitation could be ruled out. The model calculations showed that Mg2+ sorb onto the magnetite and maghemite surfaces as a mixture of mono- or bidentate surface complexes at 0.5 Mg2+ site?1 and as monodentate complexes at 1 and 2 Mg2+ site?1 conditions. Mg2+ was also found to adsorb more readily at the maghemite surfaces in comparison with magnetite surfaces. For experiments with excess Mg2+ relative to the number of surface sites, the calculations suggested the formation of polynuclear surface complexes on maghemite.