A hybrid ensemble data assimilation scheme (HYBRID), combining a flow-dependent with a static background covariance, was developed and implemented for assimilating satellite (SeaWiFS) Chl-a data into a marine ecosystem model of the Mediterranean. The performance of HYBRID was assessed against a model free-run, the ensemble-based singular evolutive interpolated Kalman (SEIK) and its variant with static covariance (SFEK), with regard to the assimilated variable (Chl-a) and non-assimilated variables (dissolved inorganic nutrients). HYBRID was found more efficient than both SEIK and SFEK, reducing the Chl-a error by more than 40% in most areas, as compared to the free-run. Data assimilation had a positive overall impact on nutrients, except for a deterioration of nitrates simulation by SEIK in the most productive area (Adriatic). This was related to SEIK pronounced update in this area and the phytoplankton limitation on phosphate that lead to a built up of excess nitrates. SEIK was found more efficient in productive and variable areas, where its ensemble exhibited important spread. SFEK had an effect mostly on Chl-a, performing better than SEIK in less dynamic areas, adequately described by the dominant modes of its static covariance. HYBRID performed well in all areas, due to its “blended” covariance. Its flow-dependent component appears to track changes in the system dynamics, while its static covariance helps maintaining sufficient spread in the forecast. HYBRID sensitivity experiments showed that an increased contribution from the flow-dependent covariance results in a deterioration of nitrates, similar to SEIK, while the improvement of HYBRID with increasing flow-dependent ensemble size quickly levels off.