A heterogeneous model was developed to describe interactions between ozone and hydrophobic organic compounds, exemplified by pentachlorophenol, in highly gas-saturated vadose zones where water moisture was limited to a thin film on soil particle surfaces. The soil was assumed to be free of soil organic matter. The model included a set of transient equations considering diffusion with simultaneous chemical reaction and hydrophobic partitioning. From dimensionless analysis, it was found that the film concentrations of ozone and the hydrophobic organic component were dependent on the Damköhler numbers. Effects of Damköhler numbers on the film profiles of components were examined. With the interfacial flux of ozone calculated from film profiles, dimensionless governing equations of ozone transport and contaminant removal across an experimental column were established. These equations were dependent on the Stanton number. One-dimensional column experiments were conducted to test the model. The optimal time for flow rate adjustment during the process was approximated. Finally, effects of ozone velocity and ozone gas concentration on the Stanton number were evaluated.