Developing accurate and computationally efficient models for ocean acoustics is inherently challenging due to several factors including the complex physical processes and the need to provide results on a large range of scales. Furthermore, the ocean itself is an inherently dynamic environment within the multiple scales. Even if we could measure the exact properties at a specific instant, the ocean will continue to change in the smallest temporal scales, ever increasing the uncertainty in the ocean prediction. In this work, we explore ocean acoustic prediction from the basics of the wave equation and its derivation. We then explain the deterministic implementations of the Parabolic Equation, Ray Theory, and Level Sets methods for ocean acoustic computation. We investigate methods for evolving stochastic fields using direct Monte Carlo, Empirical Orthogonal Functions, and adaptive Dynamically Orthogonal (DO) differential equations. As we evaluate the potential of Reduced-Order Models for stochastic ocean acoustics prediction, for the first time, we derive and implement the stochastic DO differential equations for Ray Tracing (DO-Ray), starting from the differential equations of Ray theory. With a stochastic DO-Ray implementation, we can start from non-Gaussian environmental uncertainties and compute the stochastic acoustic ray fields in a reduced order fashion, all while preserving the complex statistics of the ocean environment and the nonlinear relations with stochastic ray tracing. We outline a deterministic Ray-Tracing model, validate our implementation, and perform Monte Carlo stochastic computation as a basis for comparison. We then present the stochastic DO-Ray methodology with detailed derivations. We develop varied algorithms and discuss implementation challenges and solutions, using again direct Monte Carlo for comparison. We apply the stochastic DO-Ray methodology to three idealized cases of stochastic sound-speed profiles (SSPs): constant-gradients, uncertain deep-sound channel, and a varied sonic layer depth. Through this implementation with non-Gaussian examples, we observe the ability to represent the stochastic ray trace field in a reduced order fashion.