Understanding broad-scale species distribution and finer-scale ecological interactions is essential for conservation. We assessed species richness, distribution, habitat use and interspecific associations of cetacean in the Bohol Sea, Philippines. During 72 days of dedicated survey (2010 - 2013), we encountered 12 species of cetacean in 291 sightings, 16.8% of which involved mixed species. We used maximum entropy (MaxEnt) models to assess species’ habitat suitability and found slope and distance from the coast to be influential contributors to cetacean distribution. To explore habitat use, through foraging ecology and niche segregation of sympatric species, we compared behavioral budgets across species and found significant differences (chi-sq = 21.44; p-value = 0.044). We then used GLMs to determine the foraging likelihood in relation to oceanographic features, group size and presence of associated species. Results from model selection complimented those derived from MaxEnt. However, some inter-specific exclusion behavior might also occur. Overall, our study suggests that the Bohol Sea supports a high cetacean biodiversity while more complex inter-specific dynamics might further shape species’ ecological niches. These results highlight the importance of multi-species ecology and can be used to develop management actions
