# Role of Kuroshio Current in fish resource variability off southwest Japan

## Abstract

Western boundary currents in the subtropics play a pivotal role in transporting warm water from the tropics that contribute to development of highly diverse marine ecosystem in the coastal regions. As one of the western boundary currents in the North Pacific, the Kuroshio Current (hereafter the Kuroshio) exerts great influences on biological resource variability off southwest Japan, but few studies have examined physical processes that attribute the coastal fish resource variability to the basin-scale Kuroshio variability. Using the high-quality fish catch data and high-resolution ocean reanalysis results, this study identifies statistical links of interannual fish resource variability off Sukumo Bay, Shikoku island of Japan, to subsurface ocean temperature variability in the Kuroshio. The subsurface ocean temperature variability off the south of Sukumo Bay exhibits vertically coherent structure with sea-surface height variability, which originates from the westward-propagating oceanic Rossby waves generated through surface wind anomalies in the Northwest Pacific. Although potential sources of the atmospheric variability remain unclarified, the remotely-induced oceanic Rossby waves contribute to fish resource variability off Sukumo Bay. These findings have potential applications to other coastal regions along the western boundary currents in the subtropics where the westward-propagating oceanic Rossby waves may contribute to coastal ocean temperature variability.

## Introduction

Western boundary currents in the subtropics contribute to the establishment of rich marine ecosystem along the coastal regions through poleward transport of warm water from the tropics1. The warm water along the western boundary currents also plays a crucial role in the formation of a relatively warm and moist climate in the western side of ocean basin2,3. One of the western boundary currents in the subtropical North Pacific, called the Kuroshio Current (hereafter the Kuroshio), provides favorable environment for highly diverse marine ecosystems off Philippines, Taiwan, and the south of Japan4,5,6. Since these countries are traditionally dependent on fish resources from the Kuroshio as sources of protein foods, understanding physical processes in the Kuroshio and its potential link with fish resource variability is greatly important.

The Kuroshio variability has been extensively studied across different temporal and spatial scales using observation data and modelling techniques7,8,9. In particular, off the south of Japan, the Kuroshio frequently undergoes small and large meanderings in a meridional orientation10,11,12,13,14,15, which have great influences on spatial distribution of fish species4,16. Although there are several physical processes proposed to explain the meridional fluctuations of the Kuroshio11,17,18,19,20,21,22, westward-propagating oceanic Rossby waves from the Northwest Pacific associated with ocean density variations are considered to play an important role in triggering small meanderings off Kyushu island in southwest Japan19,21,22,23. This variation may further cause small meanders off the southwest of Shikoku island24 and large meanders off Kii Peninsula13,17 on interannual timescales. However, most of the previous studies limit their discussions to the underlying physical processes behind the Kuroshio variability, and potential impacts on the amount of fish resources off the south of Japan have been poorly understood.

One of the coastal bays located off southwest Japan, Sukumo Bay, is known to be greatly influenced by the Kuroshio variability25,26. Climatologically, the Kuroshio off Kyushu island first approaches south of Sukumo Bay, then turns eastward along the south coast of Japan. The northward approach of the Kuroshio brings warm water into Sukumo Bay, and the interaction with nutrient-rich cold water from Bungo Channel provides favorable conditions for biological activity off Sukumo Bay. Physical processes and predictability of warm water intrusion in Sukumo Bay are well-examined in the previous studies25,26,27,28,29, but the potential influence on the fish resource variability has yet to be investigated in the context of the links with basin-scale Kuroshio variability.

To bridge a gap in our current understanding of the relationship between coastal fish resources and basin-scale ocean current variability, this study aims to (i) establish statistical links of the Kuroshio variability to the fish resources off Sukumo Bay and (ii) identify potential sources of the Kuroshio variability. For these purposes, we utilize high-resolution ocean reanalysis results in the Northwest Pacific with the ability to resolve the Kuroshio variability involving mesoscale eddies. Also, high-quality fish catch data with daily catch efforts off Sukumo Bay are used to estimate monthly fish resources and their variability. The identified relationship would be beneficial for establishing fish resource prediction and management off Sukumo Bay, and have potential applications to other coastal regions along the western boundary currents in the subtropics.

## Results

### Impact of local processes on fish resource variability

Fish resources off Sukumo Bay undergo pronounced seasonal to interannual variability. Monthly fish catch per unit effort (CPUE) during 2006–2018 exhibits remarkable year-to-year variations with high values in 2007 and 2015 and low values in 2009 and 2010 (Fig. 1a). The interannual variability of the CPUE off Sukumo Bay is found to be due mostly to the wintertime variability. This can be clearly seen in Fig. 1b showing a strong seasonality of the CPUE with the large standard deviation during November-January. Further analysis of fish species data also shows that during boreal winter, small fishes such as sardine and horse mackerel become the dominant species off Sukumo Bay due probably to the combined effects of the nutrient-rich southward current from Bungo Channel and the northward warm water intrusion of the Kuroshio. Given that most of the annual CPUE is explained by the wintertime catches, this study focuses on identifying physical processes that control interannual variability in the wintertime CPUE off Sukumo Bay.

Spatial patterns of November-January mean sea-surface height (SSH) and subsurface ocean temperature at 150 m depth (T150) are presented in Fig. 2a,b. The Kuroshio, the core of which is estimated by the strongest gradient of the SSH, flows northeastward off Shikoku island. The spatial distribution of T150 mostly follows a similar pathway of the Kuroshio, although some coastal region (132–132.5°E and 32–32.5°N) on the western flank of the Kuroshio shows relatively warm temperature, representing northward intrusion of warm water from the Kuroshio to Bungo Channel. This region well corresponds to the south of major fishing areas during boreal winter.

The wintertime CPUE off Sukumo Bay shows moderately high correlation (0.4–0.5) with the SSH and T150 off the south of Sukumo Bay (Fig. 2c,d). Here we calculated the correlation between the November-January mean SSH/T150 data from 2006 to 2018 (i.e. 13 values) as two-dimensional (i.e. latitude/longitude) oceanic variables from the JCOPE reanalysis and the November-January mean one-dimensional CPUE data after aggregating all the CPUE data reported in different regions off the Sukumo Bay. Assuming that the one-dimensional CPUE data represents the area-averaged fish resource off the Sukumo Bay, we constructed the two-dimensional correlation coefficients between the CPUE and the SSH/T150 over the grid cells of the JCOPE model, so the estimated correlation coefficients have 11 degrees of freedom. It should also be noted that lag autocorrelation of November-January CPUE shows a very weak relation with the CPUE beyond two months (Fig. S1), so the wintertime CPUE variability is not much affected by other seasons’ CPUE variability. Although some region near the southwestern tip of Shikoku island shows relatively high correlation, both the SSH and the T150 off the south of Sukumo Bay (132–132.5°E and 32–32.5°N) show a coherent structure with moderately high correlations along the Kuroshio. In particular, the correlation coefficient with the T150 is statistically significant at 90% confidence level, although that with the SSH is not significant. Even if the outlier year of 2016 associated with extremely high CPUE (Fig. 1a) is removed from the analysis, the correlation values with the SSH and T150 remain relatively high above 0.4 off the south of Sukumo Bay (Fig. S2a,b), although the correlation coefficient with the T150 is statistically significant at 80% due to a decrease in the number of available data. This significant relationship with the T150 suggests that the above-normal subsurface ocean temperature associated with the northward approach of the Kuroshio may play an important role in providing favorable conditions for increase in the wintertime fish resources off Sukumo Bay. Since there is no clear relationship between the wintertime CPUE and the area-averaged Chlorophyll-a off the south of Sukumo Bay (Fig. S3), subsurface ocean temperature variations may directly contribute to changes in spatial distribution of small fishes off Sukumo Bay.

### Remote influences on Kuroshio variability

Subsurface ocean temperature variability off Sukumo Bay shows strong association to the overlying SSH variability. This is evident in Fig. 3 showing time series of the T150 and SSH anomalies during boreal winter (November-January). It should be noted that we extend the time series back to 1993 when the satellite SSH data become available and incorporated into the ocean reanalysis results. The T150 anomalies remarkably fluctuate year by year, and exhibit a distinct in-phase relationship with the SSH anomalies. To investigate a physical link with the SSH variability and its potential sources, we define positive and negative events as years when the T150 anomalies during boreal winter exceed above one standard deviation and below negative one standard deviation, respectively. This leads to seven positive and six negative events, including higher than normal CPUE years (2007, 2015, and 2016) and lower than normal CPUE years (2013) in Fig. 1a. Since the negative events include the limited number of lower-than-normal CPUE year, we focus on the analysis of physical processes for the positive events.

The in-phase relationship with the SSH variability is found for subsurface ocean temperature at different depths. During positive events, the ocean temperature in the upper 700 m shows a warmer-than-average condition compared to that in all the analysis years (Fig. S4a). Since the ocean density becomes lower than that in all the analysis years but does not show much difference below 700 m (Fig. S4b), the thermocline located at around 200–300 m, defined as the maximum vertical gradient of ocean density below the surface mixed-layer, appears to become deeper. The associated downwelling induces above-normal subsurface ocean temperature through moving the isotherm layers downward. A similar but opposite relationship can be seen for the negative events. The strong coherence between the SSH and subsurface ocean temperature involving the thermocline variability is also reported in the observational studies over the Kuroshio recirculation region30.

To obtain useful insights into the generation of the SSH anomalies, spatial patterns of the SSH anomalies composited with every three-month lag during the positive events are presented in Fig. 4. The positive SSH anomalies off the south of Sukumo Bay during boreal winter (November-January; NDJ (0)) exhibit strong association with positive SSH anomalies along the Kuroshio pathway. In particular, the positive SSH anomalies off the south of Sukumo Bay appear to be a part of positive SSH anomalies east of Kyushu island. These anomalies seem to originate from the positive SSH anomalies to the southeast at three-month lag (132–136°E and 29–31°N in Fig. 4b) and further east at six-month lag (138–142°E and 29–31°N in Fig. 4c). On the other hand, there is very weak contribution of the SSH anomalies off the southwest Kyushu island at these lags. This lag-composite analysis suggests that most of the positive SSH anomalies off the south of Sukumo Bay may be related to westward migration of the positive SSH anomalies from the Northwest Pacific.

To highlight the migration of the SSH anomalies from the Northwest Pacific, we plot two Hovmöller diagrams for north-south SSH anomalies along 132°E and west-east SSH anomalies along 30°N as a function of time lag (Fig. 5). The positive SSH anomalies at zero-month lag (i.e. Dec (0)) off the south of Sukumo Bay are associated with clear northward propagation of the anomalies with two-month lag (left panel in Fig. 5). This indicates that the positive SSH anomalies migrate northward along the east coast of Kyushu island under the influence of the northward Kuroshio (Fig. 4a). The positive SSH anomalies further originate from the westward-propagating SSH anomalies with 12-month lag in the region of 150–160°E at 30°N (right panel in Fig. 5). The westward propagation speed is estimated to be around 10 cm s−1. This is higher, by a factor of two, than the theoretical phase speed (around 5 cm s−1) for the first baroclinic Rossby waves at 30°N, probably due to the acceleration of the oceanic Rossby waves interacting with the mean current field31. As such, the oceanic Rossby waves generated in the Northwest Pacific may contribute to the generation of the positive SSH anomalies off Sukumo Bay.

To explore the potential sources of the oceanic Rossby waves from the perspective of atmospheric forcing, we calculated composite anomalies with a one-year lag prior to the positive events for sea-level pressure (SLP) in Fig. 6a and surface wind stress curl in Fig. S5, respectively. The SLP anomalies exhibit significant positive values in the region of the Northwest Pacific (150–160°E at 30°N; Fig. 6a). This leads to anomalous anticyclonic wind stress curl that tends to induce downwelling oceanic Rossby waves, although the negative wind stress curl anomalies along 30°N show significant values in the very limited areas of 140–150°E (Fig. S5). The positive SLP anomalies in the Northwest Pacific are associated with negative sea-surface temperature (SST) anomalies to the northeast and positive SST anomalies to the southwest (Fig. 6b). These SST anomalies appear not to force the atmospheric variability, rather to be driven by the atmospheric forcing in such a way that anomalously dry northwesterly wind on the northeastern flank of the positive SLP anomalies in the Northwest Pacific increases evaporation and deepens the wintertime mixed layer, which enhances entrainment of cold water from deeper ocean. A similar but opposite process seems to operate for the positive SST anomalies on the southwestern flank of the positive SLP anomalies. Therefore, the SLP anomalies in the Northwest Pacific may be driven by remote forcing outside the Northwest Pacific, for example, atmospheric teleconnection from the tropical Pacific associated with El Niño-like condition (Fig. 6b), but potential sources of SLP anomalies require further detailed analysis and modelling studies.

## Discussions

This study has identified the potential roles of the Kuroshio variability on the wintertime fish resource variability estimated off Sukumo Bay and the remote influence from the atmospheric variability in the Northwest Pacific. Previous studies have attributed fish resource variability to regional variations in the physical and bio-geochemical conditions4,5,6,16, but little attention has been paid to its remote link with basin-scale ocean and atmosphere variability. The Kuroshio undergoes interannual variability under the influence of westward-propagating oceanic Rossby waves from the Northwest Pacific19,21,22,23. Since the Kuroshio variability is suggested to cause interannual modulation of small meanders off the southwest of Shikoku island24, the remotely induced oceanic Rossby waves have the potential to modulate the ocean circulation near the southwest coast of Japan and the amount of fish resource through changes in the subsurface ocean temperature. The identified relationship can be applied to other coastal regions along the western boundary currents in the subtropics where the westward-propagating oceanic Rossby waves may contribute to coastal ocean temperature variability32.

However, some issues behind the fish resource variability remain to be addressed. First, due to the limited availability of high-quality fish catch data (2006–2018), the established relationship between the fish resource and the ocean variables may be influenced by the outlier year, for example, 2016 with extremely high fish catch (Fig. 1a), but more prolonged datasets in the near future will help to verify the robustness of the statistical relationship. Second, the relationship between fish resource increase and anomalous warm water intrusion due to the northward approach of the Kuroshio may be straightforward, but due to the lack of coastal ocean observations and reanalysis results for bio-geochemical components, it is difficult to examine how the warm water affects the phyto- and zoo-plankton activities and the amount of fish species. Since the area off Sukumo Bay is largely influenced by nutrient-rich southward current from Bungo Channel during boreal winter, anomalous warm water advection from the Kuroshio into the channel may provide favorable temperature conditions for spawning grounds of small fishes such as sardines33 and enhance biomass and food availability for small larvae34. The warm water intrusion into Bungo Channel may be associated with more frequent occurrence of the Kyucho25, a coastal phenomenon with a sudden increase in coastal ocean current, but the Kyucho in Bungo Channel rarely occurs during boreal winter. Along this line, further observational and modelling studies involving the interaction of dynamical and bio-geochemical processes would advance our understanding of the fish resource variability off Sukumo Bay.

Previous studies have mainly focused on westward propagation of cyclonic oceanic Rossby waves and eddies along 30°N accompanied with the Kuroshio pathway variations southeast of Kyushu island23,24. On the other hand, the present study has also identified potential roles of westward-propagating anticyclonic eddies in the Kuroshio variations. The oceanic Rossby waves that influence the Kuroshio variability are not solely generated by atmospheric forcing but through internal ocean processes such as eddy-mean current interaction31. Since the subtropical Northwest Pacific at 30°N also receives influence of southward current associated with re-circulation from the Kuroshio extension current system (Fig. S6), the anticyclonic eddies traveling from the north may propagate westward via interaction with the re-circulation current. However, the relative contributions from the internal oceanic processes and the atmospheric forcing remain unclear. This needs to await further ocean modelling studies in which the atmospheric forcing such as wind stress curl is prescribed with or without interannual variations.

This study provides further implication for the development of fish resource prediction based on the ocean current information. Since there is a one-year lag relationship between anticyclonic eddies in the Northwest Pacific and anomalous increase in fish resource off Sukumo Bay, monitoring the SSH variability in the Northwest Pacific is imperative for predicting fish resource variability off Sukumo Bay one year ahead. Given that the SSH variability in the Northwest Pacific is driven mostly by the atmospheric variations, seasonal climate prediction over the Northwest Pacific using a global ocean-atmosphere coupled model may help extend the prediction lead time beyond one year. The long-term prediction information for fish resource would benefit fishery people to efficiently establish fishing plan as well as sustainably use and manage fish resources.

## Methodology

We analyzed daily fish catch data based on mid-size surrounding nets off Sukumo Bay, Shikoku island of Japan from 2006 to 2018. The data covers the period since 2006 when all the local fishery cooperatives near Sukumo Bay were integrated into the Sukumo Bay fishery cooperative. To estimate fish resource variability, we calculated the monthly fish catch per unit effort (CPUE) defined as the fish catch divided by the number of fishing days per person. For oceanic and atmospheric data over the global domain, we used the sea-surface temperature (SST) from the Optimum Interpolation SST version 2 (OISST V2)35 and the basic atmospheric variables from the ERA-Interim reanalysis results36. To estimate biological activity in the upper ocean, we used the monthly Chlorophyll-a data with 4-km horizontal resolution obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite37. Here we analyzed all the above datasets with the same horizontal resolution of 1° × 1° over the domain.

To examine coastal ocean variability, we utilized monthly reanalysis results from the Japan Coastal Ocean Predictability Experiment 2 (JCOPE2)38 with a high horizontal resolution of 1/12° in the Northwest Pacific (108–180°E and 10.5–62°N). The JCOPE2 system is based on the Princeton Ocean Model (POM)39 with 46 vertical levels of sigma coordinate. The boundary condition of the JCOPE2 system is provided by relatively low-resolution model with a 1/4° horizontal resolution and 21 vertical levels in the entire Pacific. The low-resolution model was spun-up for 15 years using the monthly mean surface forcing from an initial condition with no motion, annual mean ocean temperature and salinity40. The spin-up results of the low-resolution model over the last five years are used for the lateral boundary condition of the high-resolution model, then both the low and high resolution models are integrated from Oct 1992 using surface wind stress, heat and salt fluxes of the six-hourly NCEP/NCAR reanalysis data41 via bulk formulae42. Besides atmospheric forcing, the ocean model benefited from data assimilation via three-dimensional variational assimilation (3DVAR)43 method using ocean observation data of sea-surface height (SSH) anomaly derived from several satellites, the SST from the Advanced Very High Resolution Radiometer/Multi-Channel SST (AVHRR/MCSST), and subsurface ocean temperature/salinity from the NOAA Global Temperature-Salinity Profile Program (GTSPP). Here we analyzed JCOPE2 reanalysis results during 1993–2018. To calculate monthly anomalies, we subtracted monthly climatology and removed a linear trend using a least squares method.

For the analysis of correlation between the fish resource (x) and the oceanic variables (y), we used Pearson product-moment correlation in which the least-squares regression line was calculated and the degree of the line fitting was evaluated using the least-squares method. Here, the correlation coefficient, R, was defined using the following equation:

$${\rm{R}}=\frac{\sum x{\text{'}}_{i}\,y{\text{'}}_{i}}{\sqrt{\sum x{\text{'}}_{i}^{2}}\sqrt{\sum y{\text{'}}_{i}^{2}}}$$
(1)

where the superscript’ of each variable means monthly detrended anomalies and the subscript i indicates the monthly time series.

## References

1. 1.

Tittensor, D. P. et al. Global patterns and predictors of marine biodiversity across taxa. Nature. 466, 1098–1101 (2010).

2. 2.

Minobe, S., Kuwano-Yoshida, A., Komori, N., Xie, S. P. & Small, R. J. Influence of the Gulf Stream on the troposphere. Nature. 452, 206–209 (2008).

3. 3.

Beal, L. M. et al. On the role of the Agulhas system in ocean circulation and climate. Nature. 472, 429–436 (2011).

4. 4.

Nakata, K., Hada, A. & Matsukawa, Y. Variations in food abundance for Japanese sardine larvae related to the Kuroshio meander. Fish. Oceanogr. 3, 39–49 (1994).

5. 5.

Zenimoto, K. et al. The effects of seasonal and interannual variability of oceanic structure in the western Pacific North Equatorial Current on larval transport of the Japanese eel Anguilla japonica. J. Fish. Biol. 74, 1878–1890 (2009).

6. 6.

Hsiao, S. H., Fang, T. H., Shih, C. T. & Hwang, J. S. Effects of the Kuroshio Current on copepod assemblages in Taiwan. Zool. Stud. 50, 475–490 (2011).

7. 7.

Mizuno, K. & White, W. B. Annual and interannual variability in the Kuroshio current system. J. Phys. Oceanogr. 13, 1847–1867 (1983).

8. 8.

Qiu, B. & Lukas, R. Seasonal and interannual variability of the North Equatorial Current, the Mindanao Current, and the Kuroshio along the Pacific western boundary. J. Geophys. Res. Oce. 101, 12315–12330 (1996).

9. 9.

Feng, M., Mitsudera, H. & Yoshikawa, Y. Structure and variability of the Kuroshio Current in Tokara Strait. J. Phys. Oceanogr. 30, 2257–2276 (2000).

10. 10.

Nitani, H. Variation of the Kuroshio south of Japan. J. Oceanogr. Soc. Japan. 31, 154–173 (1975).

11. 11.

White, W. B. & McCreary, J. P. On the formation of the Kuroshio meander and its relationship to the large-scale ocean circulation. Deep-Sea Res. 23, 33–47 (1976).

12. 12.

Solomon, H. Occurrence of small “trigger” meanders in the Kuroshio off southern Kyushu. J. Oceanogr. 34, 81–84 (1978).

13. 13.

Kawabe, M. Sea level variations along the south coast of Japan and the large meander in the Kuroshio. J. Oceanogr. Soc. Japan. 36, 97–104 (1980).

14. 14.

Sekine, Y. & Toba, Y. Velocity variation of the Kuroshio during formation of the small meander south of Kyushu. J. Oceanogr. Soc. Japan. 37, 87–93 (1981).

15. 15.

Kawabe, M. Variations of current path, velocity, and volume transport of the Kuroshio in relation with the large meander. J. Phys. Oceanogr. 25, 3103–3117 (1995).

16. 16.

Nakata, H., Funakoshi, S. & Nakamura, M. Alternating dominance of postlarval sardine and anchovy caught by coastal fishery in relation to the Kuroshio meander in the Enshu‐nada Sea. Fish. Oceanogr. 9, 248–258 (2000).

17. 17.

Sekine, Y. A numerical experiment on the path dynamics of the Kuroshio with reference to the formation of the large meander path south of Japan. Deep-Sea Res. 37, 359–380 (1990).

18. 18.

Qiu, B. & Miao, W. Kuroshio path variations south of Japan: Bimodality as a self-sustained internal oscillation. J. Phys. Oceanogr. 30, 2124–2137 (2000).

19. 19.

Mitsudera, H., Waseda, T., Yoshikawa, Y. & Taguchi, B. Anticyclonic eddies and Kuroshio meander formation. Geophys. Res. Lett. 28, 2025–2028 (2001).

20. 20.

Ebuchi, N. & Hanawa, K. Influence of mesoscale eddies on variations of the Kuroshio path south of Japan. J. Oceanogr. 59, 25–36 (2003).

21. 21.

Usui, N., Tsujino, H., Fujii, Y. & Kamachi, M. Generation of a trigger meander for the 2004 Kuroshio large meander. J. Geophys. Res. Oce. 113, C01012, https://doi.org/10.1029/2007JC004266 (2008).

22. 22.

Usui, N., Tsujino, H., Nakano, H. & Fujii, Y. Formation process of the Kuroshio large meander in 2004. J. Geophys. Res. Oce 113, C08047, https://doi.org/10.1029/2007JC004675 (2008).

23. 23.

Usui, N., Tsujino, H., Nakano, H. & Matsumoto, S. Long-term variability of the Kuroshio path south of Japan. J. Oceanogr. 69, 647–670 (2013).

24. 24.

Kashima, M. et al. Quasiperiodic small meanders of the Kuroshio off Cape Ashizuri and their inter-annual modulation caused by quasiperiodic arrivals of mesoscale eddies. J. Oceanogr. 65, 73–80 (2009).

25. 25.

Takeoka, H., Akiyama, H. & Kikuchi, T. The Kyucho in the Bungo Channel, Japan—Periodic intrusion of oceanic warm water. J. Oceanogr. 49, 369–382 (1993).

26. 26.

Akiyama, H. & Saitoh, S. I. TheKyucho in Sukumo Bay induced by Kuroshio warm filament intrusion. J. Oceanogr. 49, 667–682 (1993).

27. 27.

Isobe, A., Guo, X. & Takeoka, H. Hindcast and predictability of sporadic Kuroshio‐water intrusion (kyucho in the Bungo Channel) into the shelf and coastal waters. J. Geophys. Res. Oce. 115, C04023, https://doi.org/10.1029/2009JC005818 (2010).

28. 28.

Nagai, T. & Hibiya, T. Numerical simulation of tidally induced eddies in the Bungo Channel: A possible role for sporadic Kuroshio-water intrusion (kyucho). J. Oceanogr. 68, 797–806 (2012).

29. 29.

Nagai, T. & Hibiya, T. Effects of tidally induced eddies on sporadic Kuroshio-water intrusion (kyucho). J. Oceanogr. 69, 369–377 (2013).

30. 30.

Ebuchi, N. & Hanawa, K. Mesoscale eddies observed by TOLEX-ADCP and TOPEX/POSEIDON altimeter in the Kuroshio recirculation region south of Japan. J. Oceanogr. 56, 43–57 (2000).

31. 31.

Dewar, W. K. On “too fast” baroclinic planetary waves in the general circulation. J. Phys. Oceanogr. 28, 1739–1758 (1998).

32. 32.

Hill, K. L., Robinson, I. S. & Cipollini, P. Propagation characteristics of extratropical planetary waves observed in the ATSR global sea surface temperature record. J. Geophys. Res. Oce. 105, 21927–21945 (2000).

33. 33.

Kimura, R., Watanabe, Y. & Zenitani, H. Nutritional condition of first-feeding larvae of Japanese sardine in the coastal and oceanic waters along the Kuroshio Current. ICES J. Mar. Sci. 57, 240–248 (2000).

34. 34.

Nakata, K., Zenitani, H. & Inagake, D. Differences in food availability for Japanese sardine larvae between the frontal region and the waters on the offshore side of Kuroshio. Fish. Oceanogr. 4, 68–79 (1995).

35. 35.

Reynolds, R. W., Rayner, N. A., Smith, T. M., Stokes, D. C. & Wang, W. An improved in situ and satellite SST analysis for climate. J. Climate. 15, 1609–1625 (2002).

36. 36.

Dee, D. P. et al. The ERA‐Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteorol. Soc. 137, 553–597 (2011).

37. 37.

Hu, C., Lee, Z. & Franz, B. Chlorophyll a algorithms for oligotrophic oceans: A novel approach based on three‐band reflectance difference. J. Geophys. Res. Oce. 117, C01011, https://doi.org/10.1029/2011JC007395 (2012).

38. 38.

Miyazawa, Y. et al. Water mass variability in the western North Pacific detected in a 15-year eddy resolving ocean reanalysis. J. Oceanogr. 65, 737–756 (2009).

39. 39.

Mellor, G. L. Häkkinen, S. M. Ezer, T. & Patchen, R. C. A generalization of a sigma coordinate ocean model and an intercomparison of model vertical grids. Oce. Forecast. 55–72 (2002).

40. 40.

Conkright, M. E. et al. World Ocean Database, 2001. Volume 1: Introduction. NOAA Atlas NESDIS 42 (ed. Levitus, S.) 1–167 (U.S. Government Printing Office, 2002).

41. 41.

Kalnay, E. et al. The NCEP/NCAR 40-year reanalysis project. Bull. Amer. Meteorol. Soc. 77, 437–472 (1996).

42. 42.

Kagimoto, T. Miyazawa, Y. Guo, X. & Kawajiri, H. High resolution Kuroshio forecast system: Description and its applications. High Resolution Numerical Modelling of the Atmosphere and Ocean (ed. Ohfuchi, W. & Hamilton, K.) 209–239 (Springer, 2008).

43. 43.

Fujii, Y. & Kamachi, M. A reconstruction of observed profiles in the sea east of Japan using vertical coupled temperature-salinity EOF modes. J. Oceanogr. 59, 173–186 (2003).

## Acknowledgements

The JCOPE2 was run on the Earth Simulator at the Japan Agency for Marine Earth Science and Technology (JAMSTEC). All figures were generated using Grid Analysis and Display System (GrADS) Version 2.1.a3 (http://cola.gmu.edu/grads/downloads.php). We appreciate the Sukumo Bay fishery cooperative and the Sukumo Bay fishery guidance office for providing high-quality fish catch data. We are also grateful to Drs. Toshio Yamagata and Toru Miyama, and two anonymous reviewers for providing insightful comments to improve the current research. The present research is supported by the Ocean Policy Research Institute, the Sasakawa Peace Foundation, and in part by the Environment Research and Technology Development Fund S-15 (Predicting and Assessing Natural Capital and Ecosystem Services; PANCES) of the Ministry of the Environment, Japan.

## Author information

Y.M. conducted research and wrote main manuscript. S.V. and Y.M. guided the whole research. All the authors reviewed the manuscript.

Correspondence to Yushi Morioka.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions