Operational flood forecasting requires that accurate estimates of the uncertainty associated with modelgenerated streamflow forecasts be provided along with the probable flow levels. This paper demonstrates a stochastic ensemble implementation of the Sacramento model used routinely by the National Weather Service for deterministic streamflow forecasting. The approach, the simultaneous optimization and data assimilation method (SODA), uses an ensemble Kalman filter (EnKF) for recursive state estimation allowing for treatment of streamflow data error, model structural error, and parameter uncertainty, while enabling implementation of the Sacramento model without major modification to its current structural form. Model parameters are estimated in batch using the shuffled complex evolution metropolis stochasticensemble optimization approach (SCEM-UA). The SODA approach was implemented using parallel computing to handle the increased computational requirements. Studies using data from the Leaf River, Mississippi, indicate that forecast performance improvements on the order of 30% to 50% can be realized even with a suboptimal implementation of the filter. Further, the SODA parameter estimates appear to be less biased, which may increase the prospects for finding useful regionalization relationships.