loop species richness model over all sites