A parameterization of numerical weather prediction uncertainty is presented for use by atmospheric transport and dispersion models. The theoretical development applies Taylor dispersion concepts to diagnose dispersion metrics from numerical wind field ensembles, where the ensemble variability approximates the wind field uncertainty. This analysis identifies persistent wind direction differences in the wind field ensemble as a leading source of enhanced ''virtual'' dispersion, and thus enhanced uncertainty for the ensemble-mean contaminant plume. This dispersion is characterized by the Lagrangian integral time scale for the gridresolved, large-scale, ''outer'' flow that is imposed through the initial and boundary conditions and by the ensemble deviation-velocity variance. Excellent agreement is demonstrated between an explicit ensemblemean contaminant plume generated from a Gaussian plume model applied to the individual wind field ensemble members and the modeled ensemble-mean plume formed from the one Gaussian plume simulation enhanced with the new ensemble dispersion metrics.