Application of a fish habitat model considering mesoscale oceanographic features in evaluating climatic impact on distribution and abundance of Pacific saury (Cololabis saira)

•A new predictor describing MOFs is firstly introduced to improve fish habitat models.•The key role of MOFs in fish habitat is illustrated in an example of Pacific saury.•Fish habitat models incorporating MOFs can explain the abundance variations of saury.•Both oceanic isotherms and the Oyashio extension regulate habitat’s northward shift.•Recent suitable habitat changes may induce the rapid decline of saury catch. The continuing development of appropriate environmental predictors is important to improving the performance of fish habitat models. Mesoscale oceanographic features (MOFs) such as fronts, eddies, or upwelling zones, sometimes appearing as “oases in a fluid desert” occupying >50 km ocean surface, have customarily been regarded as key factors driving fish distributions. However, previous algorithms quantifying oceanographic features have not provided greatly enhanced predictive power to the fish habitat models which are characterized by a low variable importance due to inadequate attention having been paid to its scale-specific aspects. Here, a new predictor representing MOFs is introduced and a reappraisal by common fish habitat models of the role of the MOFs in fish distribution within a case study on Pacific saury (Cololabis saira) in the northwestern Pacific is reported. When the new predictor with MOFs was introduced, the performances of three commonly used fish habitat models were all improved significantly (26%–30%), together with the paramount importance of MOFs among all predictive variables, suggesting that MOFs may exert significant effects on Pacific saury distribution. Using the optimal model (Random Forest) selected among the three, the habitat distribution of Pacific saury in Oyashio water was then reconstructed for the period of 1993–2020 to explore its spatio-temporal variations and its relationship with the abundance variation. The “suitable habitats” (areas of reconstructed catch > 12 tons) during major fishing seasons (August–November) estimated by the optimal habitat model for Pacific saury showed a clear northward shift of 0.045°/year for the period of 1993–2020, associated with the poleward movement of both oceanic isotherms and the Oyashio extension due to northward movement of the whole wind field in concert with global warming. The estimated annual suitable habitat area (SHA) during the early fishing seasons (June–September) showed large inter-annual variations with a peak around 2010 and a valley in 2015, and correlated