The state of the atmosphere is of utmost importance for infrasound propagation. In propagation modeling, the true state of the atmosphere is mainly represented by the analysis. The analysis is the best deterministic estimate of the atmosphere using a data assimilation system existing of a General Circulation Model (GCM). However, the analysis excludes error variances of both model and observations. In addition, the coarse resolution of GCM results in averaging of, e.g., clouds or gravity waves, over larger regions known as parameterization. Consequentially, arrivals due to fine-scale structure in wind and temperature can be missing. Therefore, infrasound propagation including the atmospheric best-estimate error variances based on an ensemble model is proposed. The ensemble system exists of model perturbations with an amplitude comparable to the analysis error estimates to obtain a probability density function. The best-estimate analysis error variances are described by a set of perturbations using the European Centre for Medium-range Weather Forecasts (ECMWF) Ensemble Data Assimilation (EDA) system. Probabilistic infrasound propagation is demonstrated by one year of mining activity, e.g., blasting, in Gällivare, northern Sweden, observed at infrasound array IS37 in Norway, part of the International Monitoring System (IMS) for verification of the Comprehensive Nuclear-Test-Ban Treaty (CTBT). Probabilistic infrasound propagation is compared with the standard deterministic result obtained using the analysis.