Anchovy Engraulis encrasicolus distribution in European waters spans from the Mediterranean Sea to the North Sea, and is expected to expand further north with global warming. Observations from the eastern Mediterranean (North Aegean Sea), the Bay of Biscay and the North Sea reveal latitudinal differences in growth, maximum size, fecundity and timing of reproduction. We set up a mechanistic framework combining a bioenergetics model with regional physical-biogeochemical models providing temperature and zooplankton biomass to investigate the underlying mechanisms of variation in these traits. The bioenergetics model, based on the Dynamic Energy Budget theory and initially calibrated in the Bay of Biscay, was used to simulate growth and reproduction patterns. Environment partly explained the increased growth rate and larger body size towards the north. However, regional calibration of the maximum assimilation rate was necessary to obtain the best model fit. This suggests a genetic adaptation, with a pattern of cogradient variation with increasing resource towards the north, in addition to a countergradient thermal adaptation. Overall, the seasonal energy dynamics supports the pattern of body-size scaling with latitude, i.e. food-limited growth but low maintenance costs in the warm Aegean Sea, and larger size in the North Sea allowing sufficient storage capacity for overwintering. Further, the model suggests a synchronisation of reproductive timing with environmental seasonality as a trade-off between thresholds of temperature and reserves for spawning and overwintering, respectively. Finally, low temperature, short productive and spawning seasons, and insufficient reserves for overwintering appear to be current limitations for an expansion of anchovy to the Norwegian Sea.