Subduction zone plate boundaries extend for hundreds to thousands of kilometers along strike fueling volcanic arcs on the overriding plate. Slab inputs to the mantle that are continuous along strike give rise to discrete volcanoes with variable distance from the plate boundary and along-strike spacing (e.g., Lee & Wada, 2017;O'Hara et al., 2020) as well as compositional heterogeneity within and between different volcanoes (e.g., Wanke et al., 2019). Heterogeneity also occurs at intermediate scales in which groups of adjacent volcanoes with common geochemical or eruptive characteristics define along-strike segments (O'Hara et al., 2020;Schmidt et al., 2008). It is unclear how these aspects of volcanic arc heterogeneity are linked to deep crustal magma reservoirs, which process mantle melt inputs into their eventual volcanic or intrusive products. Magma reservoirs beneath volcanic arcs are thought to span the entire crustal depth range and create long-lived hot zones, although the specific organization of melt accumulations is transient (Cashman Abstract Volcanic arcs consist of many distinct vents that are ultimately fueled by the common melting processes in the subduction zone mantle wedge. Seismic imaging of crustal-scale magmatic systems can provide insight into how melt is organized in the deep crust and eventually focused beneath distinct vents as it ascends and evolves. Here, we investigate the crustal-scale structure beneath a section of the Cascades arc spanning four major stratovolcanoes: Mt. Hood, Mt. St. Helens (MSH), Mt. Adams (MA), and Mt. Rainier, based on ambient noise data from 234 seismographs. Simultaneous inversion of Rayleigh and Love wave dispersion constrains the isotropic shear velocity (Vs) and identifies radially anisotropic structures. Isotropic Vs shows two sub-parallel low-Vs zones (∼3.45-3.55 km/s) at ∼15-30 km depth with one connecting Mt. Rainier to MA, and another connecting MSH to Mt. Hood, which are interpreted as deep crustal magma reservoirs containing up to ∼2.5%-6% melt, assuming near-equilibrium melt geometry. Negative radial anisotropy, from vertical fractures like dikes, is prevalent in this part of the Cascadia, but is interrupted by positive radial anisotropy, from subhorizontal features like sills, extending vertically beneath MA and Mt. Rainier at ∼10-30 km depth and weaker and west-dipping positive anisotropy beneath MSH. The positive anisotropy regions are adjacent to rather than co-located with the isotropic low-Vs anomalies. Ascending melt that stalled and mostly crystallized in sills with possible compositional differences from the country rock may explain the near-average Vs and positive radial anisotropy adjacent to the active deep crustal magma reservoirs.Plain Language Summary Volcanic arcs, a common result of subduction processes, comprise a large proportion of active volcanoes in the world and pose significant hazards. Seismic tomography measures variations of seismic wave speed in the subsurface, which can then be used to infer important properties of the vo...