Spatio-temporal multivariate occupancy models can handle multiple species in occupancy models. This method for fitting such models is described in Hepler and Erhardt (2021) "A spatiotemporal model for multivariate occupancy data".
Please use the canonical form https://CRAN.R-project.org/package=multiocc to link to this page.