Lava flows comprise three-phase mixtures of melt, crystals, and bubbles. While existing one-phase treatments allow melt phase viscosity to be assessed on the basis of composition, water content, and/or temperature, two-phase treatments constrain the effects of crystallinity or vesicularity on mixture viscosity. However, three-phase treatments, allowing for the effects of coexisting crystallinity and vesicularity, are not well understood. We investigate existing one- and two-phase treatments using lava flow case studies from Mauna Loa (Hawaii) and Mount Etna (Italy) and compare these with a three-phase treatment that has not been applied previously to basaltic mixtures. At Etna, melt viscosities of 425 ± 30 Pa s are expected for well-degassed (0.1 w. % H(2)O), and 135 ± 10 Pa s for less well-degassed (0.4 wt % H(2)O), melt at 1080°C. Application of a three-phase model yields mixture viscosities (45% crystals, 25-35% vesicles) in the range 5600-12,500 Pa s. This compares with a measured value for Etnean lava of 9400 ± 1500 Pa s. At Mauna Loa, the three-phase treatment provides a fit with the full range of field measured viscosities, giving three-phase mixture viscosities, upon eruption, of 110-140 Pa s (5% crystals, no bubble effect due to sheared vesicles) to 850-1400 Pa s (25-30% crystals, 40-60% spherical vesicles). The ability of the three-phase treatment to characterize the full range of melt-crystal-bubble mixture viscosities in both settings indicates the potential of this method in characterizing basaltic lava mixture viscosity.