Abstract Fish spawning phenology is a major concern for conservation and fisheries management. New intensive data sources, such as GPS-based tracking data and high-resolution catch declaration data, are becoming increasingly available in the field of marine ecology. These data benefit from high spatiotemporal resolution and open new research avenues for investigating the interannual variability in fish phenology. In this paper, we demonstrate how an integrated species distribution model informed by commercial catch data combined with spatiotemporal dimension reduction methods known as empirical orthogonal functions (EOFs) can be used to synthesize spatiotemporal signals in fish reproduction phenology. Specifically, we address the following questions: (1) Can we identify seasonal spatial patterns that can be interpreted in terms of reproductive phenology and essential habitats? (2) Can we identify changes in reproductive phenology over time? (3) Are these changes related to environmental drivers? The analysis illustrates the reproductive phenology of three key commercial species in the Bay of Biscay (sole, hake, and sea bass). The EOF analysis emphasized strong seasonal spatiotemporal patterns that correspond to reproduction patterns and feeding patterns. Based on this methodology, we identified seasonal variations in the timing of reproduction, and we related these variations to sea surface temperature, a key driver of fish reproduction.