Conserving oceanic apex predators, such as sharks, is of utmost importance. However, scant abundance and distribution data often challenge understanding the population status of many threatened species. Occurrence records are often scarce and opportunistic, and fieldwork aimed to retrieve additional data is expensive and prone to failure. Integrating various data sources becomes crucial to developing species distribution models for informed sampling and conservation purposes. The white shark, for example, is a rare but persistent inhabitant of the Mediterranean Sea. Here, it is considered Critically Endangered by the IUCN, while population abundance, distribution patterns, and habitat use are still poorly known. This study uses available occurrence records from 1985 to 2021 from diverse sources to construct a spatial log‐Gaussian Cox process, with data‐source specific detection functions and thinning, and accounting for physical barriers. This model estimates white shark presence intensity alongside uncertainty through a Bayesian approach with Integrated Nested Laplace Approximation (INLA) and the inlabru R package. For the first time, we projected species occurrence hot spots and landscapes of relative abundance (continuous measure of animal density in space) throughout the Mediterranean Sea. This approach can be used with other rare species for which presence‐only data from different sources are available.