A Monte Carlo algorithm has been developed to investigate the effects of multiple scattering on the volume scattering function measured by the LISST-VSF instrument. The developed algorithm is compared to experimental results obtained from bench-top measurements using 508nm spherical polystyrene beads and Arizona test dust as scattering agents. The Monte Carlo simulation predicts measured volume scattering functions at all concentrations. We demonstrate that multiple scattered light can be a major contributor to the detected signal, resulting in errors in the measured volume scattering function and its derived inherent optical properties. We find a relative error of 10% in the scattering coefficient for optical depths ∼0.4, and it can reach 100% at optical depths ∼2.