Positive magnetic anomalies seen in MAGSAT average scalar anomaly data overlying some subduction zones can be explained in terms of the magnetization contrast between the cold subducted oceanic slab and the surrounding hotter, nonmagnetic mantle. Three-dimensional modeling studies show that peak anomaly amplitude and location depend on slab length and dip. A model for the Aleutian Arc anomaly matches the general trend of the observed MAGSAT anomaly if a slab thickness of 7 km and a relatively high (induced plus viscous) magnetization contrast of 4 A/m are used. A second source body along the present day continental margin is required to match the observed anomaly in detail, and may be modeled as a relic slab from subduction prior to 60 m.y. ago.