We present new 3-D P wave and S wave velocity models of the upper 20 km of the Mount St. Helens (MSH) region. These were obtained using local-source arrival time tomography from earthquakes and explosions recorded at 70 broadband stations deployed as part of the imaging Magma Under St. Helens (iMUSH) project and augmented by several data sets. Principal features of our models include (1) low P wave and S wave velocities along the St. Helens seismic zone to depths of at least 20 km corresponding to high conductivity imaged by iMUSH magnetotelluric studies. This delineates a zone of weakness that magma can exploit at the location of MSH; (2) a 5- to 7-km diameter, 6–15 km deep, 3–6% negative P wave and S wave velocity anomaly beneath MSH, consistent with previous estimates of the source region for recent eruptions. We interpret this as a magma storage region containing up to 15–20 km³ of partial melt, which is about 5 times more than the largest documented eruption at MSH; (3) a broad region of low P wave velocity below 10-km depth extending between Mount Adams and Mount Rainier along and to the east of the main Cascade arc, which is likely due to high-temperature arc crust and possible presence of fluids or melt; (4) several anomalies associated with surface-mapped features, including high-velocity igneous units such as the Spud Mountain and Spirit Lake plutons and low velocities in the Chehalis sedimentary basin and the Indian Heaven volcanic field. Our results place further constraints on the geometry of these features at depth.