Three-dimensional dynamic gas-liquid flow simulations that accurately track the phase interface are numerically challenging. This article presents a numerical study of the performance of the level-set phase-interface tracking method when combined with extremely high order (7th to 11th) weighted essentially non-oscillatory (WENO) advection schemes for gas-liquid free surface flows. Comparisons between simulation results and prior benchmark results suggest that such a combination of methods can be satisfactorily applied to the level-set and Navier-Stokes equations for free surface flow simulations when volume conservation is enforced at every time step, and minor numerical oscillations are suppressed through use of an artificial viscosity term. In particular, simulations of solid body rotation, the unsteady flow following an ideal dam break, tank sloshing, and the rise of a single bubble all agree with analytical or experimental results to within + 3.12% when the level-set method is combined with an 11th order WENO scheme. Furthermore, use of an 11th order WENO advection scheme actually has a computational cost advantage because, for the same accuracy, it can be used on a coarser grid when compared with a more-common second-order advection scheme; computational savings of up to 87% are possible.