Because all models are a simplification of the phenomenon they aim to represent, it is often more useful to estimate the probability of an event rather than a single "best" model result. Previous air quality ensemble approaches have used computationally expensive simulations of separately developed modeling systems. We present an efficient method to generate ensembles with hundreds of members based on several structural configurations of a single air quality modeling system. We use the Decoupled Direct Method in three dimensions to directly calculate how ozone concentrations change as a result of changes in input parameters. The modeled probability estimate is compared to observations and is shown to have a high level of skill and improved resolution and sharpness. This approach can help resolve the practical limits of incorporating uncertainty estimation into deterministic air quality management modeling applications.