Potential impact of ocean circulation on the declining Japanese eel catches
Recruitment of Japanese eels, Anguilla japonica, has declined in recent decades possibly due to both anthropogenic and ocean-atmosphere factors. The potential impact of ocean circulation on the decreasing Japanese eel catches in the western North Pacific was examined based on a three-dimensional particle-tracking method, in which virtual larvae (v-larvae) were programmed to swim horizontally and vertically, in addition to being transported by ocean currents after being released in their North Equatorial Current (NEC) spawning area. Transport patterns varied among years between 1993 and 2013, and dispersion of v-larvae towards East Asia decreased in the last two decades, especially for the western Taiwan and Japan regions. In recent years, instead of entering the Kuroshio and moving towards East Asia as in the 1990s’, more v-larvae tended to enter the southern areas due to the weakening of the NEC and strengthening of subsurface southward flow near the spawning area. Changes in ocean circulation in the western Pacific appear to be caused by the weakening of subtropical and tropical wind stress curl in the past two decades. This suggests that decadal changes in ocean circulation have occurred that affect the larval migration success of the Japanese eel to their recruitment areas.
Anguillid eels are widely distributed in the Indo-Pacific and North Atlantic regions where the juveniles live in freshwater and estuarine habitats1,2. Their adults migrate offshore to spawn over deep water and their larvae, called leptocephali have a long larval dispersal stage that varies among tropical and temperate species3,4,5. World eel recruitment has significantly declined in the past few decades6, raising concerns about what are the major causes. Overfishing and habitat loss due to human activities were considered to be one of the primary reasons causing the eel recruitment decline7. Changes in ocean circulation or oceanic conditions may also significantly affect eel recruitment8,9,10,11,12,13,14,15, or both types of factors have likely contributed16.
The Japanese eel, Anguilla japonica, is one of the most important eel species for fisheries and aquaculture, and it has been listed as endangered on the IUCN red list17. Among the Northern Hemisphere species of eels that have all experienced declines, the Japanese eel seems to have shown the earliest recruitment declines that started in the 1970s18. After about 2010, the annual recruitment had decreased by as much as 90% compared to eel catches in the 1960s17.
The Japanese eel is distributed across East Asia in areas that are adjacent to the western Pacific Ocean. Their spawning area is located along the West Mariana Ridge within the westward flowing North Equatorial Current (NEC)19,20. The NEC forms a zone of westward flow from about 10–20°N, which then bifurcates into northward flow (the Kuroshio) and southward flow (the Mindanao Current)21,22. The eel larvae are carried primarily by the NEC and Kuroshio toward their growth habitats in East Asia, although it is possible that some active swimming occurs5. The maturing silver eels migrate thousands of kilometers to return to their spawning area for reproduction3. The migration process of Japanese eels, however, is not well understood due to limited direct observational evidence on adult migration23 or larval dispersal24. Numerical modelling methods, therefore, have been used to simulate the potential migration paths of Japanese eels, and to evaluate their linkages to ocean-atmosphere changes12,14,15,25,26.
Previous numerical modelling studies on Japanese eel leptocephali examined the effect of ocean circulation and ocean-atmosphere conditions on larval transport. For example, some studies suggested that the transport of larvae carried by currents from the NEC to the Kuroshio is lowest in the El Niño years compared to the non-El Niño periods11,12. Other recent modelling studies evaluated the interannual migration success of Japanese eel larvae from their spawning area to East Asia and suggested that more complex ocean circulation factors related to the NEC, Kuroshio, and eddies can impact the migration, which are connected to the Philippine-Taiwan Oscillations (PTO)14,15. The strength and position of the NEC appear to be important factors14,15,27, and tropical rainfall affecting the position of the salinity front that can form in the spawning area may also affect recruitment in some years27. Therefore, ocean-current and climate-related factors need further investigation using long-term simulations.
The western Pacific Ocean is one of the most variable ocean regions, having the fastest sea level rise in the world28. The western North Pacific Ocean is characterized by major oceanic features, such as the NEC, the Kuroshio, the Subtropical Countercurrent (STCC), and eddies, which will likely vary during a changing climate. For example, the Kuroshio has warmed two times faster than the global warming rate in the past century29. The NEC bifurcation latitude is shifting southward in the past 60 years30. The Kuroshio along the East China Sea has gradually moved onshore31, and the Luzon Strait intrusion into the South China Sea is weakening in the last two decades32.
The recruitment of the Japanese eel has been declining during these years even after the some of the major anthropogenic impacts have been reduced in at least some areas, so the effects of the changing ocean need further exploration. The present study investigates the potential impact of ocean circulation on the decreasing Japanese eel catches in the western North Pacific. The 22-year (1993–2013) common period for model reanalysis and annual glass eel catches is used. Note that the present work does not strictly obey the definition of the time period required to study climate effects due to the limitation of data. The long-term (1993–2013) dispersal of virtual larvae (v-larvae) in the western Pacific is simulated based on a three-dimensional (3D) particle tracking method, in which swimming behavior is included in addition to transport by ocean currents. Biological factors affecting mortality (i.e., feeding success or predation) are not considered in the present work that focuses on the effects of ocean circulation. In addition to the yearly variation that was reported by the earlier studies, a long-term trend of decreasing dispersal towards parts of East Asia is observed in our simulations. We also tentatively examine the possible causes of recent ocean circulation changes that could be connected to climate changes.
Data and Methods
Annual Japanese glass eel (recruitment-stage early juveniles) catch data from 1993 to 2013 for Taiwan and Japan were obtained from the Taiwan Fisheries Yearbook (Fisheries Agency, Council of Agriculture, Taiwan) and Japan Aquaculture Information News (The Nihon Yoshoku Shimbun, Tokyo, Japan, which collected the official data from Ministry of Agriculture, Forestry, and Fisheries in Japan), respectively25. Ocean drifter data were obtained from the Global Drifter Program (http://www.aoml.noaa.gov/phod/dac/index.php). Surface drifters released in the NEC region from May to September during two periods (1993–1997 and 2009–2013) were selected. There were 17 and 29 drifters in the early and later periods, respectively.
The data-assimilative ocean circulation model known as the Japan Coastal Ocean Predictability Experiment 2 (JCOPE2) provides the three-dimensional currents and hydrological fields that were used in particle tracking in the present study. JCOPE2 was constructed from the Princeton Ocean Model with a generalized coordinate system33. The model domain of JCOPE2 encompasses the western North Pacific (10.5–62°N and 108–180°E), with a horizontal resolution of 1/12° (8–9 km) and 46 vertical layers. The external forcing to drive JCOPE2 includes wind stresses and net heat/freshwater fluxes at the sea surface converted from the six-hourly atmospheric reanalysis produced by the National Centers for Environmental Prediction/National Center for Atmospheric Research. Satellite and in situ temperature and salinity data were assimilated into the model based on a three-dimension variational method33. The daily JCOPE2 reanalysis fields cover the period from January 1993 to the present. Comparison of simulated trajectories of passive particles carried by JCOPE2 and observed trajectories was performed in a previous study14, which showed the satisfactory performance of JCOPE2 in simulating the 3D circulation over the western North Pacific Ocean.
The wind data used for calculating wind stress curl and Sverdrup transport was based on the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis products (https://www.esrl.noaa.gov/psd/data/reanalysis/reanalysis.shtml).
A three-dimensional (3D) particle-tracking method was used to simulate the movement of virtual eel larvae (v-larvae). The particles were carried by ocean currents in addition to having their own swimming behavior. The 3D particle-tracking scheme developed by Ohashi and Sheng34 was used in this study. The tracking scheme was based on the fourth-order Runge–Kutta method35. The tracking time step was three hours. The same tracking scheme was used previously by Chang et al.14,26,36 for investigating the migration of Japanese eel larvae and adults in the western Pacific Ocean and was also applied to simulations of the long-distance migration of adult American eels in the Atlantic Ocean37,38.
Two important swimming behaviors were considered in this study: diel vertical migration (DVM) and horizontal swimming. Eel larvae likely show DVM behavior39; i.e., they remain in upper surface waters at night and swim to deeper waters probably to avoid predators during the daytime. The age-dependent DVM and horizontal swimming followed a previous study14,40. V-larvae at night stayed at a fixed depth of 50 m. The swimming speed and daytime diving depth linearly increased with age. Swimming speed increased by 0.0006 m s−1 per day, and daytime DVM deepened by 0.75 m per day. (i.e., on day 0, swimming speed was zero, and v-larvae did not perform DVM; On day 100, swimming speed was 0.06 m s−1, and DVM was between 50 to 125 m). The maximum swimming speed was 0.14 m s−1, and the maximum daytime diving depth was 230 m to the end of the simulation. A random walk displacement was included to represent unresolved sub-grid turbulent flow and other local processes34. The estimated maximum horizontal and vertical displacements due to the random walk were 600 m and 20 m, respectively. The duration of day and night were determined by the time of sunrise (6 am) and sunset (6 pm) throughout the simulation. The swimming direction of the v-larvae was set to be the same as the local flow in the open ocean with water depths greater than 100 m as in the previous study to allow intercomparison of the results14. When the v-larvae approached the coastal and shelf waters with water depths shallower than 100 m, they were set to search for lower salinity and swim toward coastal fresher waters41. As actual eel larvae likely do not swim continuously throughout their early life history and may use a directional swimming orientation5, an additional simulation using different larval behaviors was conducted as a control to establish that larval behavior was not affecting the overall trend of decreasing transport to East Asia over time (see Discussion section).
Numerical experiments were conducted to examine the interannual variation of v-larvae transport during a simulation period from 1993 to 2013. The release region and time of v-larvae were chosen on the basis of the observed spawning area and season based on the collections of Japanese eel eggs, newly hatched preleptocephali, and spawning-condition adults along of the southern West Mariana Ridge3. The particles (v-larvae) were released at locations spread over the region of 140 to 143°E and 12 to 15°N with a separation distance of 10 km in both zonal and meridional directions, and one being released at each location. Leptocephali estimated to be 10–40 days old were found in June and July42. The releases time were set from May 1 to July 31 and each release was staggered by a time interval of 5 days during the three-month period. About 18,000 v-larvae were released each year. The Japanese eels born in summer likely reach the East China Sea and south of Japan in winter and early spring24. The migration period estimated from observations is about six to eight months. The tracking duration was set to be eight months to track v-larvae during their migration14.
The number of v-larvae reaching each of 8 designated regions (Fig. 1) were calculated as visitation frequency. The regions were (1) Subtropical Countercurrent (2) Upstream Kuroshio (3) the South China Sea and Taiwan Strait (4) the Philippine Sea (5) East China Sea and Okinawa Trough (6) South of Japan (7) Northern North Equatorial Current (8) Southern North Equatorial Current. Visitation frequency was defined as the number of v-larvae that visited each individual grid, in which each of the v-larvae would only be counted once on the same grid in each simulation.
The model sub-regions (labeled 1–8, separated by white lines) shown over the bathymetry (color) in the western North Pacific Ocean. Gray vectors show the average ocean currents in upper 200 m from year 1993 to 2013. The yellow box represents the Japanese eel spawning area. The green dots mark the rivers included in tracking scheme. Figure is created using MATLAB R2011b (http://www.mathworks.com/).
Full size image
The linear trend (percentage per year) for observed glass eel recruitment and simulated v-larvae dispersion referred to the 22-year (1993–2014) mean. The significance of linear trends was examined based on Mann-Kendall tests. Significance was defined as p < 0.05.
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
The observed annual glass eel catches in Japan and Taiwan showed inter-annual fluctuations (Fig. 2). The eel catch variations between Taiwan and Japan were significantly correlated (r = 0.69, p < 0.01). Apart from the biological factors, such as changes in demographic dynamics could potentially influence the annual recruitment, the high correlation also suggests the eel catches in Taiwan and Japan are likely influenced by the same environmental factors25, such as ocean-atmosphere conditions. In addition to year-to-year variations, a long-term decreasing trend can be seen (p < 0.01). The annual eel catches in Taiwan and Japan showed overall decreases of about 5% per year in the past two decades.
Observed recruitment-stage glass eel catches of the Japanese eel, Anguilla japonica, in Japan (red) and Taiwan (blue). Dashed lines show the linear trends (p < 0.01). Unit is in t. Figure is created using MATLAB R2011b (http://www.mathworks.com/).
Full size image
Simulated v-larvae transport paths were widely distributed in the western North Pacific during each of the 21 simulation years with some clear variations among years (Fig. 3). The greatest visitation frequencies of v-larvae occurred in the NEC where larvae were transported away from their spawning area, but in some years many larvae entered areas north of 15°N. Some v-larvae entered the Kuroshio and arrived to East Asia, but others were advected southward or remained in the NEC region until the end of the 8-month simulations. The amount of v-larvae arriving in the region south of Japan was greatest in 1993, 1997, 2002, 2004, and 2008 and was lowest or none in 2007, 2009, 2010, 2011 and 2013 as well as in some earlier years (Figs 3 and 4). For those v-larvae that reached south of Japan, the main transport routes were sometimes more near-shore (i.e., 2002) or farther away from the coast (i.e., 2004).
V-larvae visitation frequency from the year 1993 to 2013 after being released in the spawning area. Unit is number of v-larvae. Figure is created using MATLAB R2011b (http://www.mathworks.com/).
Full size image
Visitation frequency v-larvae in each sub-region. Dashed lines show the linear trends after being released in the spawning area. 1 normalized unit is 20,000 v-larvae. Percentage numbers indicate the overall rate of change per year. Figure is created using MATLAB R2011b (http://www.mathworks.com/).
Full size image
Classifications of the v-larvae dispersion by region (Figs 1 and 4) showed that there were decreasing trends of dispersion into many of the sub-regions north of the spawning area (region 1–7, Fig. 4). Due to the alternations between higher and lower numbers of v-larvae entering the regions, the regressions were not significant (p = 0.10–0.78) except for in the South China Sea and Taiwan Strait region (p < 0.01) and the South of Japan region (p = 0.04), both of which have a direct impact on the eel catches in Taiwan and Japan, respectively. In contrast, the number of v-larvae were increasing to the south of the spawning area (region 8) in the past 21 years. As the release number of v-larvae per year was fixed, the changes of v-larvae distributions indicated that v-larvae tended to be advected more southward in recent years in comparison to the earlier decades, resulting in reductions of v-larvae recruitment along the Kuroshio. Cross-correlations of v-larvae visitation frequencies were observed between most sub-regions (Table 1) except the STCC (region 1). V-larvae variation in the STCC was not related to the Kuroshio system but was linked to the NEC. The amount of v-larvae remaining in the STCC region depends on the bifurcation of the NEC as described in a recent study14. V-larvae in the Kuroshio associated regions (region 2–6) were all correlated with each other. The v-larvae entry into the areas south of Japan (region 6) was significantly associated with v-larvae variation near Taiwan (region 2 and 3), similar to the relationship observed in glass eel catches (Fig. 2). V-larvae variation in the southern NEC was negatively-correlated with all other sub-regions, suggesting that increases of v-larvae to the south of the spawning area corresponds to decreases of v-larvae to East Asia.
Table 1 Cross-correlation between the subregions shown in Fig. 1.
Full size table
V-larvae dispersion towards East Asia has been reduced in recent years as illustrated by a comparison between 1993–1997 and 2009–2013 (Fig. 5). In the 1990’s, v-larvae were able to reach the East China Sea in 7 months, whereas v-larvae over the past few years could only approach the east coast of the Philippines in that amount of time. The intrusion through Luzon Strait to the South China Sea, as well as the final arrival to areas south of Japan both declined in recent years.
Mean ages of all v-larvae crossing a given grid (days) for the (a) years 1993–1997, and (b) years 2009–2013 of v-larvae released in the spawning area. Figure is created using MATLAB R2011b (http://www.mathworks.com/).
Full size image
The decreases of v-larvae being transported to the margins of East Asia appear to be caused by changes of ocean current patterns. In the past two decades, the westward flow of the northern NEC has been getting weaker (Fig. 6a–c, positive eastward trend between 12–15°N suggested that the westward NEC was weakening), while the subsurface (below 50 m) southward flow south of 14°N was getting stronger (Fig. 6d–f, negative trend south of 14°N indicated more southward flow). The weakening of the NEC would make v-larvae stay in the NEC longer, and the enhanced southward flow, although weak, could effectively bring more v-larvae southward. We traced the origin of the v-larvae entering the southern NEC (Fig. 6g), and many originated from 12–14°N where the enhanced southward current was observed (Fig. 6d). In addition to v-larvae trajectories, the limited numbers of observed drifters showed generally consistent patterns with the changes of v-larvae transport patterns in the past two decades (Fig. 6h). Drifters in the early decade (1993–1997) mostly joined the Kuroshio after departing from the northern NEC, but drifters launched in recent years (2009–2013) partly moved southward without entering the Kuroshio.
Linear trends from 1993 to 2013 (a,d) and 5-year averaged vertical profiles (b,c,e,f) in the NEC (130–143°E mean) region of zonal (a–c) and meridional velocities (d–f) for the years 1993–1997 (b,e) and years 2009–2013 (c,f). Positive and negative values for zonal (meridional) velocity represent eastward (northward) and westward (southward) directions, respectively. The trend for those v-larvae entering southern NEC is shown in (g). Positive (negative) trends suggested increase (decrease) of v-larvae entering southern NEC. Observed ocean drifter trajectories in the years 1993–1997 (blue) and years 2009–2013 (red) is shown in (h). Figure is created using MATLAB R2011b (http://www.mathworks.com/).
Full size image
The changing ocean circulation in the western North Pacific is likely caused by the forcing of the large-scale winds. The classic theory of Sverdrup dynamics explains the relationship between the wind stress curl (WSC) and the western boundary current (i.e., strengthening of WSC leads to the enhancement of western boundary current)43. The dipole of WSCs exist in the subtropical and tropical Pacific, with anti-cyclonic WSC to the north of 15°N and cyclonic WSC to the south of 15°N (Fig. 7a). The corresponding southward Sverdrup transports in the subtropical Pacific, and northward Sverdrup transports in the tropical Pacific then merge and form the western boundary currents that are the Kuroshio and the Mindanao Current. Decadal variations appeared in both subtropical and tropical WSC and the corresponding Sverdrup transports in the past 60 years (Fig. 7b,c). The subtropical and tropical WSC weakened in the previous two decades, the corresponding weakening of Sverdrup transport lead to the reduction of ocean interior currents, which contributed to the weakening of the NEC, and the weakening of subtropical WSC also led to the weakening of the Kuroshio. The weakening of northward Sverdrup transport in the tropical Pacific could serve as a candidate, explaining the stronger southward flow in the subsurface water south of 14°N if the vertical shear was fixed (i.e., the weakening of vertically integrated northward transport would result in weaker northward current near the surface and stronger southward current in the subsurface, Fig. 6e,f). The strengthening of subsurface southward flow could also occur with the change of vertical shear, in order to balance with the change of the horizontal density field43. Overall, the weakening of the NEC together with the strengthening of southward subsurface current could have resulted in greater transport of v-larvae to the areas south of the spawning area in recent years. Furthermore, the weakening of the Kuroshio in the recent decade also contributed to reduction of v-larvae northward dispersion. The Kuroshio along the East China Sea also appears to be weakening in association with the changes of wind stress curl at 20–30°N44, which also explains why the decreasing rate of v-larvae in the downstream Kuroshio region (Fig. 4, South of Japan region 6) is greater than that in the upstream Kuroshio (Fig. 4, region 2). As a result, v-larvae dispersions in recent years were confined further to the south in comparison to the earlier decade. Apart from the Kuroshio, glass eel catches in Taiwan are also affected by the intrusion of water through the Luzon Strait and into the Taiwan Strait. While the subtropical and tropical WSCs weakened, the Luzon Strait intrusion also decreased due to the weakening of westward current and reduced eddy activity45. The weakening of Luzon Strait intrusion together with the weakening of the Kuroshio likely leads to decreases in eel catches in Taiwan.
1953–2014 mean (a) wind stress curl (color) and Sverdrup transport (black vectors) based on the NCEP reanalysis data. Annual mean integrated Sverdrup transport (Sv, 1 Sv = 106 m3s) at (b) 24°N and (c) 12°N. Red and blue dashed lines are 10-yr and 20-yr linear trends, respectively. Percentages indicate the overall rate of change per year.
Full size image
The present study explored the potential impact of western North Pacific Ocean circulation on the declining and fluctuating Japanese eel catches in East Asia based on a three-dimensional particle-tracking method, in which v-larvae were programmed to swim (horizontal swimming and vertical DVM), in addition to movements resulting from transport by ocean currents. The simulated v-larvae dispersion to different regions of the western North Pacific showed clear variations among years, and in recent decades there was a decrease in the number of v-larvae reaching areas near Japan and Taiwan where the species recruits. Previous larval transport modelling studies have found similar interannual variability and that ocean-atmosphere factors such as El Niño or the Philippines–Taiwan Oscillation (PTO) appear to influence the number of larvae that may successfully reach their recruitment areas12,14,15. Those studies and the present study examined the importance of ocean circulation to Japanese eel larvae dispersion, however we focused on the possible effect of changing winds. Our findings further suggest that changes in ocean circulation could be playing a role in the decreasing trend of glass eel catches in Japan and Taiwan using long-term simulation results. An important factor in the decreases in recruitment could be that instead of entering the Kuroshio and moving northward, more larvae may be moved southward due to weakening of the westward flow of the NEC or a southward shifting of its position and strengthening of subsurface southward flow near the spawning area in recent years. These changes in ocean circulation may be in response to changes in winds resulting from a weakening of subtropical and tropical wind stress curl in the past two decades. Other factors, such as changes in horizontal temperature or salinity balanced by the changes of vertical shear may also lead to the strengthening of subsurface southward current. A detailed dynamical analysis is necessary for understanding the actual cause for change of ocean circulation.
The weakening of wind stress curl was proposed to be connected to global warming46, and climate projections have also suggested that weakening of the NEC and the upstream Kuroshio (south of 25°N) could be associated with the increase of greenhouse gases47. However, the actual dynamical forcing and responses remain unclear, and uncertainty also exists within different models47. Moreover, the global temperature warmed up rapidly starting in the 1980′s, yet the wind stress curl showed decadal oscillations prior to the 1990’s. The 22-years of ocean reanalysis available for our study could be too short to fully represent all aspects of the long-term climate variation, therefore, the role of the warming climate in changing atmospheric conditions that influence ocean circulation requires more exploration.
Although the present study and previous transport modelling studies illustrate the possible role of ocean circulation in causing distinct variability in the transport patterns that could affect recruitment levels, there are other factors that could be involved in the actual levels of recruitment each year. For example, the general decreasing trends of the time-series of glass eel catches in Taiwan and Japan were similar to the number of larvae reaching those areas in the simulations, but their year by year variations were not significantly correlated (p > 0.1). Other factors that could be involved include the number of adult eels migrating to the spawning area each year, which has the potential to be influenced by climatic factors affecting the freshwater habitats where the juveniles live48. That could affect the number of larvae being transported out of the spawning area. The fishing effort could also vary year by year, which could potentially influence the amount of annual recruitment. Another factor that appears to be unique for the Japanese eel is that the latitude of their spawning locations appear to be influenced by the latitude of a shallow salinity front that often forms in the NEC as a result of tropical rainfall42,49,50. When a distinct front is present, spawning seems to occur south of the front, but when there is no front, the spawning occurs over a wider range of latitude51. It is also possible that changes in productivity might affect early larval survival through reductions in their marine snow food source16. Changes in productivity-related indicators have changed in the NEC region52,53, and chlorophyll-a has steadily decreased and surface temperatures have increased in the North Pacific subtropical gyre since 199854. Distinguishing the possible contributions of these types of biologically linked factors and the effects of changing currents will require further detailed research efforts.
Another biological factor that could be more important than previously considered is that of active swimming by the larvae, whose bodies are filled with energy storage material that could be used for swimming5. In order to confirm that the general findings of dispersion decreases were not sensitive to the particular swimming strategy used in the simulations, different swimming behaviors were examined. Considering that the 24-hour swimming used in the present study and previous ones14,26 could be energy consuming, and therefore not biologically accurate, and the newborn larvae may just drift with ocean currents while they are feeding and growing, an additional experiment (Exp. A) set v-larvae to start swimming when they reached 30-days old, and they only swam at night when they probably stop feeding. The linearly increasing swimming speed, DVM and swimming direction were the same as the standard simulation run. Another experiment (Exp. B) used the same strategy as Exp. A, but set a fixed swimming direction of northwestward for v-larvae to facilitate movement towards and across the Kuroshio. Both experiments showed decreasing trends of reaching areas near East Asia (Fig. 8), especially with northwest swimming (p < 0.01), which also showed less interannual variability, less southward transport, and higher numbers of v-larvae crossing out of the Kuroshio system and approaching their recruitment areas. This suggests the decreasing dispersion caused by the ocean circulation changes is robust and that larval behavior is important. The result also suggests that new research on the influence of directional swimming in biologically realistic contexts will be useful in future studies.
V-larvae visitation frequency for the years (top) 1993–1997 and (bottom) 2009–2013, and the visitation frequency in the regions adjacent to East Asia (sub-regions 2–6) for the (left) standard experiment, (middle) Exp. A and (right) Exp. B. 1 normalized unit is 20,000 v-larvae. Figure is created using MATLAB R2011b (http://www.mathworks.com/).
Full size image
Further research on the effects of ocean current changes and the influence of larval behaviors are needed not only for the Japanese eel, but also on other species, such as A. marmorata, A. bicolor pacifica, and A. luzonensis in the western Pacific as well as the eel species in the other subtropical gyres where equatorial currents and western boundary currents are important for transporting larvae5. How changes in the ocean circulation or in biological factors are influencing larval survival will be essential to help efforts to conserve these interesting species in a changing ocean-atmosphere system in the future.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Aida, K., Tsukamoto, K. & Yamauchi, K. Eel Biology. Springer Japan, 497, https://doi.org/10.1007/978-4-431-65907-5 (2003).
Tesch, F. The eel. Blackwell publishing, 408 (2003).
Tsukamoto, K. Oceanic migration and spawning of anguillid eels. Journal of fish biology 74, 1833–1852, https://doi.org/10.1111/j.1095-8649.2009.02242.x (2009).
CASArticlePubMed Google Scholar4.
Kuroki, M., Miller, M. J. & Tsukamoto, K. Diversity of early life history traits in freshwater eels and the evolution of their oceanic migrations. Can. J. Zool. 92, 749–770, https://doi.org/10.1139/cjz-2013-0303 (2014).
Article Google Scholar5.
Miller, M. J. & Tsukamoto, K. The ecology of oceanic dispersal and survival of anguillid leptocephali. Canadian Journal of Fisheries and Aquatic Sciences 74, 958–971, https://doi.org/10.1139/cjfas-2016-0281 (2017).
CASArticle Google Scholar6.
ICES. Reports of the Eifac/ICES working group on eels (WGEEL), 18–22 March 2013 in Sukarietta, Spain, 4–10 September 2013 in Copenhagen, Denmark ICES CM2013/ACOM:18. 253 pp (2013).
Dekker, W. Did lack of spawners cause the collapse of the European eel, Anguilla anguilla? Fisheries Management and Ecology 10, 365–376, https://doi.org/10.1111/j.1365-2400.2003.00352.x (2003).
Article Google Scholar8.
Baltazar-Soares, M. et al. Recruitment collapse and population structure of the European eel shaped by local ocean current dynamics. Current biology: CB 24, 104–108, https://doi.org/10.1016/j.cub.2013.11.031 (2014).
CASArticlePubMed Google Scholar9.
Bonhommeau, S. et al. Impact of climate on eel populations of the Northern Hemisphere. Marine Ecology Progress Series 373, 71–80, https://doi.org/10.3354/meps07696 (2008).
ADSArticle Google Scholar10.
Knights, B. A review of the possible impacts of long-term oceanic and climate changes and fishing mortality on recruitment of anguillid eels of the Northern Hemisphere. Science of The Total Environment 310, 237–244, https://doi.org/10.1016/s0048-9697(02)00644-7 (2003).
ADSCASArticlePubMed Google Scholar11.
Kim, H. et al. Effect of El Niño on migration and larval transport of the Japanese eel (Anguilla japonica). ICES Journal of Marine Science: Journal du Conseil 64, 1387–1395, https://doi.org/10.1093/icesjms/fsm091 (2007).
Article Google Scholar12.
Zenimoto, K., Sasai, Y., Sasaki, H. & Kimura, S. Estimation of larval duration in Anguilla spp., based on cohort analysis, otolith microstructure, and Lagrangian simulations. Marine Ecology Progress Series 438, 219–228, https://doi.org/10.3354/meps09255 (2011).
ADSArticle Google Scholar13.
Friedland, K. D., Miller, M. J. & Knights, B. Oceanic changes in the Sargasso Sea and declines in recruitment of the European eel. ICES Journal of Marine Science 64, 519–530, https://doi.org/10.1093/icesjms/fsm022 (2007).
Article Google Scholar14.
Chang, Y.-L., Sheng, J., Ohashi, K., Béguer-Pon, M. & Miyazawa, Y. Impacts of Interannual Ocean Circulation Variability on Japanese Eel Larval Migration in the Western North Pacific Ocean. PLoS ONE 10, e0144423, https://doi.org/10.1371/journal.pone.0144423 (2015).
ArticlePubMedPubMed Central Google Scholar15.
Hsu, A. C., Xue, H., Chai, F., Xiu, P. & Han, Y.-S. Variability of the Pacific North Equatorial Current and its implications on Japanese eel (Anguilla japonica) larval migration. Fisheries Oceanography 26, 251–267, https://doi.org/10.1111/fog.12189 (2017).
Article Google Scholar16.
Miller, M. J., Feunteun, E. & Tsukamoto, K. Did a “perfect storm” of oceanic changes and continental anthropogenic impacts cause northern hemisphere anguillid recruitment reductions? ICES Journal of Marine Science 73, 43–56, https://doi.org/10.1093/icesjms/fsv063 (2016).
Article Google Scholar17.
Jacoby, D. & Gollock, M. Anguilla japonica. The IUCN Red List of Threatened Species Version 2015–3 (2014).
Dekker, W. Slipping through our hands: population dynamics of the European eel, s.n. (2004).
Tsukamoto, K. Oceanic biology: Spawning of eels near a seamount. Nature 439, 929–929, http://www.nature.com/nature/journal/v439/n7079/suppinfo/439929a_S1.html (2006).
Tsukamoto, K. et al. Oceanic spawning ecology of freshwater eels in the western North Pacific. Nature communications 2, 179, https://doi.org/10.1038/ncomms1174 (2011).
ArticlePubMedPubMed Central Google Scholar21.
Qiu, B. & Chen, S. M. S. Level and Gyre Circulation Variability in the Northwestern Tropical Pacific Ocean. Journal of Physical Oceanography 42, 193–206, https://doi.org/10.1175/jpo-d-11-061.1 (2012).
ADSArticle Google Scholar22.
Qiu, B. et al. The Pacific North EquatorialCurrent: New insights from the origins of the Kuroshio and Mindanao Currents (OKMC) Project. Oceanography 28(4), 24–33, https://doi.org/10.5670/oceanog.2015.78 (2015).
Article Google Scholar23.
Manabe, R. et al. First observations of the oceanic migration of Japanese eel, from pop-up archival transmitting tags. Marine Ecology Progress Series 437, 229–240, https://doi.org/10.3354/meps09266 (2011).
ADSArticle Google Scholar24.
Shinoda, A. et al. Evaluation of the larval distribution and migration of the Japanese eel in the western North Pacific. Rev Fish Biol Fisheries 21, 591–611, https://doi.org/10.1007/s11160-010-9195-1 (2011).
Article Google Scholar25.
Han, Y.-S., Zhang, H., Tseng, Y.-H. & Shen, M.-L. Larval Japanese eel (Anguilla japonica) as sub-surface current bio-tracers on the East Asia continental shelf. Fisheries Oceanography 21, 281–290, https://doi.org/10.1111/j.1365-2419.2012.00624.x (2012).
Article Google Scholar26.
Chang, Y.-L., Miyazawa, Y. & Béguer-Pon, M. Simulating the Oceanic Migration of Silver Japanese Eels. PLoS ONE 11, e0150187, https://doi.org/10.1371/journal.pone.0150187 (2016).
ArticlePubMedPubMed Central Google Scholar27.
Tzeng, W. N. et al. Evaluation of multi-scale climate effects on annual recruitment levels of the Japanese eel, Anguilla japonica, to Taiwan. PLoS One 7, e30805, https://doi.org/10.1371/journal.pone.0030805 (2012).
ADSCASArticlePubMedPubMed Central Google Scholar28.
Peyser, C. E., Yin, J., Landerer, F. W. & Cole, J. E. Pacific sea level rise patterns and global surface temperature variability. Geophysical Research Letters 43, 8662–8669, https://doi.org/10.1002/2016GL069401 (2016).
ADSArticle Google Scholar29.
Wu, L. et al. Enhanced warming over the global subtropical western boundary currents. Nature Climate Change 2, 161–166, https://doi.org/10.1038/nclimate1353 (2012).
ADSCASArticle Google Scholar30.
Chen, Z. & Wu, L. Long-term change of the Pacific North Equatorial Current bifurcation in SODA. Journal of Geophysical Research: Oceans 117, n/a-n/a, https://doi.org/10.1029/2011JC007814 (2012).
Wang, J. & Oey, L.-Y. Inter-annual and decadal fluctuations of the Kuroshio in East China Sea and connection with surface fluxes of momentum and heat. Geophysical Research Letters 41, 8538–8546, https://doi.org/10.1002/2014GL062118 (2014).
ADSArticle Google Scholar32.
Nan, F. et al. Weakening of the Kuroshio Intrusion into the South China Sea over the Past Two Decades. Journal of Climate 26, 8097–8110, https://doi.org/10.1175/jcli-d-12-00315.1 (2013).
ADSArticle Google Scholar33.
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, https://doi.org/10.1007/s10872-009-0063-3 (2009).
Article Google Scholar34.
Ohashi, K. & Sheng, J. Investigating the Effect of Oceanographic Conditions and Swimming Behaviours on the Movement of Particles in the Gulf of St. Lawrence Using an Individual-Based Numerical Model. Atmosphere-Ocean, 1–21, https://doi.org/10.1080/07055900.2015.1090390 (2015).
Press, W., Teukolsky, S., Vetterling, W. & Flannery, B. Numerical recipes in Fortran 77: the art of scientific computing. Cambridge University Press (1992).
Chang, Y.-L., Miyazawa, Y. & Béguer-Pon, M. The dynamical impact of mesoscale eddies on migration of Japanese eel larvae. PLOS ONE 12, e0172501, https://doi.org/10.1371/journal.pone.0172501 (2017).
ArticlePubMedPubMed Central Google Scholar37.
Béguer-Pon, M. et al. Exploring the role of the physical marine environment in silver eel migrations using a biophysical particle tracking model. ICES Journal of Marine Science: Journal du Conseil, https://doi.org/10.1093/icesjms/fsv169 (2015).
Béguer-Pon, M., Ohashi, K., Sheng, J., Castonguay, M. & Dodson, J. J. Modeling the migration of the American eel in the Gulf of St. Lawrence. Marine Ecology Progress Series 549, 183–198, https://doi.org/10.3354/meps11706 (2016).
ADSArticle Google Scholar39.
Castonguay, M. & McCleave, J. Vertical distributions, diel and ontogenetic vertical migration and net avoidance of leptocephali of Anguilla and of other common species in the Sargasso Sea. Journal of Plankton Research 9, 195–214 (1987).
Article Google Scholar40.
Melià, P. et al. Integrating field data into individual-based models of the migration of European eel larvae. Marine Ecology Progress Series 487, 135–149 (2013).
ADSArticle Google Scholar41.
Tosi, L., Spampanato, A., Sola, C. & Tongiorgi, P. Relation of water odour, salinity and temperature to ascent of glass-eels, Anguilla anguilla (L.): a laboratory study. Journal of fish biology 36, 327–340, https://doi.org/10.1111/j.1095-8649.1990.tb05613.x (1990).
Article Google Scholar42.
Tsukamoto, K. Discovery of the spawning area for Japanese eel. Nature 356, 789–791 (1992).
ADSArticle Google Scholar43.
Marshall, J. P. A. Atmosphere, Ocean and Climate Dynamics. Academic Press 344 (2007).
Wang, Y.-L., Wu, C.-R. & Chao, S.-Y. Warming and weakening trends of the Kuroshio during 1993-2013. Geophysical Research Letters 43, 9200–9207, https://doi.org/10.1002/2016gl069432 (2016).
ADSArticle Google Scholar45.
Chang, Y.-L. & Oey, L.-Y. The Philippines–Taiwan Oscillation: Monsoonlike Interannual Oscillation of the Subtropical–Tropical Western North Pacific Wind System and Its Impact on the Ocean. Journal of Climate 25, 1597–1618, https://doi.org/10.1175/JCLI-D-11-00158.1 (2012).
ADSArticle Google Scholar46.
Saenko, O. A., Fyfe, J. C. & England, M. H. On the response of the oceanic wind-driven circulation to atmospheric CO2 increase. Climate Dynamics 25, 415–426, https://doi.org/10.1007/s00382-005-0032-5 (2005).
ADSArticle Google Scholar47.
Hu, D. et al. Pacific western boundary currents and their roles in climate. Nature 522, 299–308, https://doi.org/10.1038/nature14504 (2015).
ADSCASArticlePubMed Google Scholar48.
Kettle, A. J., Asbjørn Vøllestad, L. & Wibig, J. Where once the eel and the elephant were together: decline of the European eel because of changing hydrology in southwest Europe and northwest Africa? Fish and Fisheries 12, 380–411, https://doi.org/10.1111/j.1467-2979.2010.00400.x (2011).
Article Google Scholar49.
Kimura, S., Inoue, T. & Sugimoto, T. Fluctuation in the distribution of low-salinity water in the North Equatorial Current and its effect on the larval transport of the Japanese eel. Fisheries Oceanography 10, 51–60, https://doi.org/10.1046/j.1365-2419.2001.00159.x (2001).
Article Google Scholar50.
Kimura, S. & Tsukamoto, K. The salinity front in the North EquatorialCurrent: A landmark for the spawning migration of the Japanese eel (Anguilla japonica) related to the stock recruitment. Deep Sea Research Part II: Topical Studies in Oceanography 53, 315–325, https://doi.org/10.1016/j.dsr2.2006.01.009 (2006).
ADSArticle Google Scholar51.
Aoyama, J. et al. Spawning Sites of the Japanese Eel in Relation to Oceanographic Structure and the West Mariana Ridge. PLoS ONE 9, e88759, https://doi.org/10.1371/journal.pone.0088759 (2014).
ADSArticlePubMedPubMed Central Google Scholar52.
Sugimoto, T. & Tadokoro, K. Interdecadal variations of plankton biomass and physical environment in the North Pacific. Fisheries Oceanography 7, 289–299, https://doi.org/10.1046/j.1365-2419.1998.00070.x (1998).
Article Google Scholar53.
Watanabe, Y. W., Ishida, H., Nakano, T. & Nagai, N. Spatiotemporal Decreases of Nutrients and Chlorophyll-a in the Surface Mixed Layer of the Western North Pacific from 1971 to 2000. J Oceanogr 61, 1011–1016, https://doi.org/10.1007/s10872-006-0017-y (2005).
CASArticle Google Scholar54.
Signorini, S. R., Franz, B. A. & McClain, C. R. Chlorophyll variability in the oligotrophic gyres: mechanisms, seasonality and trends. Frontiers in Marine Science 2, https://doi.org/10.3389/fmars.2015.00001 (2015).
Y.M. was supported by JSPS KAKENHI (26287116). The annual eel catch data provided by Dr. Y.-S. Han and the particle tracking scheme from Drs. K. Ohashi and J. Sheng are gratefully acknowledged.
AffiliationsApplication Laboratory, Japan Agency for Marine-Earth Science and Technology, Yokohama, 236-0001, JapanYu-Lin K. Chang & Yasumasa MiyazawaDepartment of Marine Science and Resources, College of Bioresource Sciences, Nihon University, Fujisawa, Kanagawa, 252-0880, JapanMichael J. Miller & Katsumi Tsukamoto
AuthorsSearch for Yu-Lin K. Chang in:Nature Research journals • PubMed • Google Scholar Search for Yasumasa Miyazawa in:Nature Research journals • PubMed • Google Scholar Search for Michael J. Miller in:Nature Research journals • PubMed • Google Scholar Search for Katsumi Tsukamoto in:Nature Research journals • PubMed • Google Scholar
Y.L.C. conceived, designed the experiments, and analyzed the data. Y.L.C. and Y.M. performed the numerical experiments. Y.L.C., Y.M., M.M. and K.T. wrote the paper.
The authors declare no competing interests.
Correspondence to Yu-Lin K. Chang.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
08 November 2017
21 March 2018
03 April 2018