Predictions of the distribution of groundfish species are needed to support ongoing marine spatial planning initiatives in Canadian Pacific waters. Data to inform species distribution models are available from several fishery-independent surveys. However, no single survey covers the entire region and different gear types are required to survey the range of relevant habitat. Here, we demonstrated a method for integrating presence–absence data across surveys and gear types that allows us to predict the coastwide distributions of 65 groundfish species in British Columbia. Our model leverages data from multiple surveys to estimate how species respond to environmental gradients while accounting for differences in survey catchability. We find that this method has two main benefits: (1) it increases the accuracy of predictions in data-limited surveys and regions while having negligible impacts on accuracy when data are already sufficient, and (2) it reduces uncertainty, resulting in tighter confidence intervals on predicted occurrences. These benefits are particularly relevant in areas of our coast where our understanding of habitat suitability is limited due to a lack of spatially comprehensive long-term groundfish surveys.