Our method produces an ensemble of cloud index fields from which probabilities are calculated. The cloud index fields are derived from geostationary satellite images that are then advected using an ensemble of cloud motion fields that are the result of assimilating optical flow vectors and NWP winds. Probabilities of cloud cover are then estimated based on the proportion of ensemble members that contain clouds over a specified region. We explore different methods to calculate probabilities for a given point or region from the ensemble and the performance of forecasts over large and small areas. We analyze the performance of this forecasting method over Tucson, AZ using case studies, Brier scores, and reliability diagrams.