During three low-altitude (99, 66, 66 km) flybys through the Enceladus plume in 2010 and 2011, Cassini's ion neutral mass spectrometer (INMS) made its first high spatial resolution measurements of the plume's gas density and distribution, detecting in situ the individual gas jets within the broad plume. Since those flybys, more detailed Imaging Science Subsystem (ISS) imaging observations of the plume's icy component have been reported, which constrain the locations and orientations of the numerous gas/grain jets. In the present study, we used these ISS imaging results, together with ultraviolet imaging spectrograph stellar and solar occultation measurements and modeling of the three-dimensional structure of the vapor cloud, to constrain the magnitudes, velocities, and time variability of the plume gas sources from the INMS data. Our results confirm a mixture of both low and high Mach gas emission from Enceladus' surface tiger stripes, with gas accelerated as fast as Mach 10 before escaping the surface. The vapor source fluxes and jet intensities/densities vary dramatically and stochastically, up to a factor 10, both spatially along the tiger stripes and over time between flyby observations. This complex spatial variability and dynamics may result from time-variable tidal stress fields interacting with subsurface fissure geometry and tortuosity beyond detectability, including changing gas pathways to the surface, and fluid flow and boiling in response evolving lithostatic stress conditions. The total plume gas source has 30% uncertainty depending on the contributions assumed for adiabatic and nonadiabatic gas expansion/acceleration to the high Mach emission. The overall vapor plume source rate exhibits stochastic time variability up to a factor ∼5 between observations, reflecting that found in the individual gas sources/jets. Key Words: Cassini at Saturn-Geysers-Enceladus-Gas dynamics-Icy satellites. Astrobiology 17, 926-940.