Understanding interannual variability in the distribution of, and transport processes affecting, the early life stages of Todarodes pacificus using behavioral-hydrodynamic modeling approaches

Abstract

To understand interannual variability in the distribution of the early life stages of Todarodes pacificus summer spawning population, and to identify the key transport processes influencing this variability, we used a coupled bio-physical model that combines an individual-based model (IBM) incorporating ontogenetic vertical migration for paralarval behavior and temperature-dependent survival process with a ROMS oceanographic model. Using the distribution of paralarvae observed in the northern East China Sea (ECS) during several field cruises as an end point, the spawning ground for the summer-spawning population was estimated to extend from southeast Jeju Island to the central ECS near 29°N by running the model backwards in time. Running the model forward, interannual variability in the distribution of paralarvae predicted by the model was consistent with that observed in several field surveys; surviving individuals in the northern ECS were substantially more abundant in late July 2006 than in 2007, in agreement with observed paralarval distributions. The total number of surviving individuals at 60days after release based on the simulation throughout summer spawning period (June–August) was 20,329 for 2006, compared with 13,816 for 2007. The surviving individuals were mainly distributed in the East/Japan Sea (EJS), corresponding to a pathway following the nearshore branch of the Tsushima Warm Current flowing along the Japanese coast during both years. In contrast, the abundance of surviving individuals was extremely low in 2007 compared to 2006 on the Pacific side of Japan. Interannual variability in transport and survival processes made a substantial impact on not only the abundance of surviving paralarvae, but also on the flux of paralarvae to adjacent waters. Our simulation results for between-year variation in paralarval abundance coincide with recruitment (year n+1) variability of T. pacificus in the field. The agreement between the simulation and field data indicates our model may be useful for predicting the recruitment of T. pacificus.