Current ocean acoustic simulation environments often rely on a single depth-varying sound speed profile. This work introduces a robust and configurable Octave/C + + implementation of a four-dimensional sound speed environment (4D SSP) for use with Bellhop in the MIT/LAMSS software package for acoustic simulations. This tool shows considerable advances in that the 4D SSP environment can be created from MIT GCM model output, CTD data, or manually varied from existing historical profiles. By considering the fluctuations in sound speed in longitude, latitude, depth, and time, output pressure data better approximates experimental data. Two experiments in the Arctic and the Santa Barbara Channel (SBC) are simulated in 1-D and 4-D to compare to experimental data taken. A framework for acoustic data assimilation using this new simulated data environment in tandem with experimental data is shown to derive the true sound speed field. [Work supported by ONR under the Information in Ambient Noise MURI.]