We show that, in the Hardy space on the unit disc, any Bessel system of pairwise separated model spaces is a Riesz system. In case such model spaces are generated by finite Blaschke products, we will see how this gives a characterization for interpolating sequences of square matrices.