A Computational Fluid Dynamics (CFD) model was developed to describe the behavior of a subsea gas release and the subsequent rising gas plume. Four numerical approaches were assessed for their suitability to capture the characteristic behaviors in a rising gas plume by comparing the CFD results with experimental data obtained from an underwater gas release experiment carried out in a towing tank basin with 10 m depth. The k-epsilon turbulence model was found to be unsatisfactory in capturing random wandering behavior of the subsea gas plume due to the inherent Reynolds-Averaged Navier-Stokes (RANS) nature of the approach. The result is an over-prediction of the plume central line velocity and an under -prediction of the plume width as there was no mechanism to distribute and dissipate the high momentum gained during the initial gas release phase. The results obtained using the Large Eddy Simulation (LES) approach show the inherently random wandering behavior of the plume is successfully captured and both the centerline velocity and the velocity profile are in much better agreement with the experimental data. The study was carried out using air as released gas. This was done to compare with the available experimental data where air was used as the source. However the CFD model developed is applicable for actual subsea gases. (C) 2017 Elsevier Ltd. All rights reserved.