Structural heterogeneity is a common feature of all glasses, however, little is known about the underlying contributions of chemical fluctuations and modulations in free volume in concrete glass forming systems. In this investigation, we relate the dynamics of structural relaxation of (100-x)NaPO 3-xAlF 3 glasses to their heterogeneous structure as determined from multinuclear magnetic resonance spectroscopic analysis. For this, we evaluate differential scanning calorimetry (DSC) data using the integral isoconversional method to determine the variation in activation energy, E a , of the glass transition as a function of temperature and conversion progress. Specific heat measurements from DSC allow for the determination of the effective size of the cooperatively rearranging region (CRR). From 31 P, 19 F, and 27 Al NMR, we observe that the introduction of AlF 3 into the NaPO 3 network increases the average connectivity (i.e., the number of heteronuclear Al-O-P bonds), rationalizing the higher E a determined from the DSC measurements. We find highly constrained regions of Al(OP) 4 F 2 with Al-F-Al cross-linking (high E a) and, simultaneously, more flexible regions of phosphate chains containing P-F••(Na +) n bonds (low E a); this results in a topologically and dynamically heterogeneous structure as evidenced by the increased variability in E a with higher AlF 3 content. The decreasing size of the CRR reflects the increased heterogeneity: at low AlF 3 the CRR is large, while at high AlF 3 (high heterogeneity), the CRR is significantly smaller (by a factor of 10 3). Finally, we relate the heterogeneity to other macroscopic properties, such as T g and mechanical properties.