The analysis scheme was cycled for 90 days, with new observations assimilated every 12 h. Generally, it was found that the analysis performs best when background error covariances are estimated almost fully from the ensemble, especially when the ensemble size was large. When small-sized ensembles are used, some lessened weighting of ensemble-based covariances is desirable to reduce the impact of spurious analysis corrections over long distances. The relative improvement over 3D-Var analyses was dependent upon the observational data density; generally, there is less improvement for data-rich networks than for data poor networks, with the largest improvement for the network with the data void. As expected, errors depend on the size of the ensemble, with errors decreasing as more ensemble members are added. The sets of initial conditions generated from the hybrid are generally well-calibrated and provide an improved set of initial conditions for ensemble forecasts.