Many ocean acoustic propagation models assume an idealized pressure-release boundary at the ocean surface. This is easily accomplished numerically using finite element, finite difference, and split-step Fourier techniques. For models based on the split-step Fourier algorithm, rough surfaces can also be treated, but require additional complexity in the definitions of the field and propagator functions through a field transformation technique. For this reason, it may be advantageous to model a rough water/air interface in a manner analogous to the bottom treatment by extending the calculation into the air medium, thereby simplifying the definitions of the propagator functions. However, standard approaches to treat density discontinuities in split-step Fourier algorithms invoke smoothing functions, which have been shown to introduce phase errors in range. In this work, the hybrid split-step/finite-difference approach introduced by Yevick and Thomson (1996) is implemented in the Monterey-Miami Parabolic Equation (MMPE) model for both the water/sediment and water/air interfaces. Particular attention is paid to comparisons between the rough surface scattering results from the field transformation technique (applied to a pressure release surface) and the hybrid split-step/finite-difference approach (applied to the water/air interface).Many ocean acoustic propagation models assume an idealized pressure-release boundary at the ocean surface. This is easily accomplished numerically using finite element, finite difference, and split-step Fourier techniques. For models based on the split-step Fourier algorithm, rough surfaces can also be treated, but require additional complexity in the definitions of the field and propagator functions through a field transformation technique. For this reason, it may be advantageous to model a rough water/air interface in a manner analogous to the bottom treatment by extending the calculation into the air medium, thereby simplifying the definitions of the propagator functions. However, standard approaches to treat density discontinuities in split-step Fourier algorithms invoke smoothing functions, which have been shown to introduce phase errors in range. In this work, the hybrid split-step/finite-difference approach introduced by Yevick and Thomson (1996) is implemented in the Monterey-Miami Parabolic Equat...