Climate change often leads to shifts in the distribution of small pelagic fish, likely by changing the match-mismatch dynamics between these sensitive species within their environmental optima. Using present-day habitat suitability, we projected how different scenarios of climate change (IPCC Representative Concentration Pathways 2.6, 4.5 and 8.5) may alter the large scale distribution of European sardine Sardina pilchardus (a model species) by 2050 and 2100. We evaluated the variability of species-specific environmental optima allowing a comparison between present-day and future scenarios. Regardless of the scenario, sea surface temperature and salinity and the interaction between current velocity and distance to the nearest coast were the main descriptors responsible for the main effects on sardine's distribution. Present-day and future potential "hotspots" for sardine were neritic zones (<250 km) with water currents <0.4 m s, where SST was between 10 and 22 °C and SSS > 20 (PSU), on average. Most variability in projected shifts among climatic scenarios was in habitats with moderate to low suitability. By the end of this century, habitat suitability was projected to increase in the Canary Islands, Iberian Peninsula, central North Sea, northern Mediterranean, and eastern Black Sea and to decrease in the Atlantic African coast, southwest Mediterranean, English Channel, northern North Sea and Western U.K. A gradual poleward-eastward shift in sardine distribution was also projected among scenarios. This shift was most pronounced in 2100 under RCP 8.5. In that scenario, sardines had a 9.6% range expansion which included waters along the entire coast of Norway up and into the White Sea. As habitat suitability is mediated by the synergic effects of climate variability and change on species fitness, it is critical to apply models with robust underlying species-habitat data that integrate knowledge on the full range of processes shaping species productivity and distribution.