Efficient numerical simulation of stochastic internal-wave-induced sound-speed perturbation fields

An efficient method is presented to numerically simulate stochastic internal-wave-induced sound-speed perturbation fields in deep ocean environments. The sound-speed perturbation field is represented as an internal-wave eigenfunction expansion in which WKB amplitude scaling and stretching of the depth coordinate are exploited. Individual realizations of the sound-speed perturbation field are constructed by evaluating a multidimensional fast Fourier transform of a complex-valued function whose modulus has a known simple form and whose phase is random. Approximations made are shown to be consistent with approximations built into the Garrett–Munk internal-wave spectrum, which is the starting point of this analysis. Both time-varying internal-wave fields in three space dimensions and frozen fields in a vertical plane are considered.