Geodetically modeled reservoir volume changes during volcanic eruptions are commonly much smaller than the observed eruptive volumes. This discrepancy is thought to be partially due to the compressibility of magma, which is largely controlled by the presence of exsolved volatiles. The 2006 eruption of Augustine Volcano, Alaska, produced an eruptive volume that was ∼3 times larger than the geodetically estimated syn-eruptive subsurface volume change. In this study, we use a multistep methodology that combines constraints from geodetic, volcanic gas, geologic, and petrologic data together with equations relating physical processes to observable parameters. We apply a Monte Carlo approach to quantify uncertainties. Ultimately, we solve for the exsolved volatile volume fraction and the magma compressibility. We estimate Augustine's 2006 pre-eruptive exsolved volatile phase to be ∼5.5 vol% of the magma at storage depths, yielding a bulk magma compressibility of ∼3.8 × 10⁻¹⁰ Pa⁻¹. We develop a novel approach to estimate the H₂O/CO₂ ratio of the syn-eruptive gas emissions in the absence of direct H₂O emission measurements which are hard to obtain due to the high background levels in ambient air. We find a best-fit H₂O/CO₂ molar ratio of 29. We also investigate the effects of applying different equations of state to our model. We find that the Ideal Gas Law might be used as a first approximation due to its simplicity; however, it overestimates volatile density and compressibility significantly at storage depths. This project capitalizes on the insights that can be gained by integrating multidisciplinary data with models of physical processes.