Seasonal Movement Patterns of the Bigfin Reef Squid Sepioteuthis lessoniana Predicted Using Statolith δ18O Values

Sepioteuthis lessoniana is a widely distributed neritic squid in the Indo-Pacific Ocean. It is an important species in fisheries, but species management is difficult because of inadequate information regarding its life history. The daily growth and δ18O values from the core to the edge of statoliths of S. lessoniana collected in northern Taiwan were analyzed to predict the experienced temperature history along with ontogenetic stages. The probability of occurrence in a given area at each life stage in three seasonal groups was determined using salinity values, deduced and measured temperatures, and the known ecology of S. lessoniana. The results showed that ontogenetic variation in the statolith δ18O values in S. lessoniana reflected the seasonal temperature fluctuation observed in Taiwanese waters, which indicated the reliability of the prediction method. Complex and diverse distribution and movement patterns were observed in the three seasonal groups. The results also indicated the importance of the waters near the coast of northeastern Taiwan as a spawning ground. Based on a model prediction, the distribution of S. lessoniana is likely associated with water temperature and current. A potential migration route from the Penghu Islands to northeastern Taiwan suggests a high level of population connectivity in S. lessoniana in Taiwan. This study provides information on the spatial and vertical distribution of S. lessoniana at various ontogenetic stages, which is essential for resource management and conservation of this commercial species.


INTRODUCTION
Planktonic paralarvae of loliginid squids are distributed widely after hatching by coastal trends and ocean currents (Young and Harman, 1988;Zeidberg and Hamner, 2002). By contrast, juveniles and adults live in coastal and inshore waters and are commonly aggregated in schools for spawning (Jereb and Roper, 2005). The distribution range of paralarvae is a critical parameter for recruitment success and gene flow in cephalopod populations (O'dor, 1998;Moreno et al., 2008), and understanding movement patterns is crucial for studies on population structure, fishery management, and marine protected area evaluation (O'dor, 1992;Semmens et al., 2007). Temperature is one of the major factors affecting the movements of squids because squids are highly vulnerable to drastic changes in temperature. Direct data of the thermal movement patterns of squids have been collected primarily through external tagging (e.g., Ueta and Jo, 1990;Jackson et al., 2005;Gilly et al., 2006;Kanamaru et al., 2007;Bazzino et al., 2010). Low recapture rates and the size of the electronic device usually limit the success of squid tagging experiments (Semmens et al., 2007). Alternatively, elemental signatures (e.g., Sr:Ca ratio) of the hard structures of squid provide a potential source of information regarding the ambient temperatures that they have been exposed to (Ikeda et al., 2003;Yamaguchi et al., 2015Yamaguchi et al., , 2018Liu et al., 2016). However, studies on the relationship between the Sr:Ca ratio and temperature have provided inconsistent results in different species Gillanders, 2002, 2003;Gillanders et al., 2013). A technique that is more commonly applied to teleosts should be considered as a supplementary method for improving the accuracy and reliability of predictions.
In recent years, stable oxygen isotopic ratios (expressed as δ 18 O) in biogenic carbonates have been used to examine the spatial distribution and movement patterns of teleosts (e.g., Steer et al., 2010;Currey et al., 2014;Shiao et al., 2017). The δ 18 O value in fish otoliths is in equilibrium with that in the ambient water, and δ 18 O value uniformly increases with a decrease in seawater temperature in various species of fish (Elsdon and Gillanders, 2002;Høie et al., 2004a). Squid statoliths are composed of calcium carbonate, and the relationships between statolith and ambient water chemistry are similar to that of fish otoliths (Arkhipkin et al., 2004;Arkhipkin, 2005;Gillanders et al., 2013). Therefore, statolith δ 18 O value can be analogously used as a temperature proxy for reconstructing the ontogenetic preferences of squids (Radtke, 1983;Landman et al., 2004;Trasviña-Carrillo et al., 2018). For example, Trasviña-Carrillo et al. (2018) analyzed the δ 18 O value in the entire statolith of the jumbo squid Dosidicus gigas and found that temperature preference did not differ between sexes but did among ontogenetic stages because of differences in living depth. Landman et al. (2004) also predicated the experienced temperature of the giant squid Architeuthis sanctipauli to be in the range of 10.5-12.9 • C and average living depths to be 125-250 m by analyzing the δ 18 O values in one entire statolith, which are considered to provide lifelong average temperature data. However, micromill sampling for higher temporal resolution of δ 18 O information (Høie et al., 2004b), which was used on fish otoliths, has never been used on squid statoliths to increase the temporal resolution of their ontogenetic thermal migration history.
The bigfin reef squid, S. lessoniana, is widely distributed in the neritic waters of the Indo-Pacific Ocean, which include the waters surrounding Taiwan (Roper et al., 1984;Okutani, 2015). Field observations have revealed that adults migrate from offshore areas to shallow inshore areas for spawning (Segawa, 1987). In the waters off the northern coast of Taiwan, S. lessoniana hatching occurs almost throughout the entire year; hatching peaks in May and August-September (Chen et al., 2015). Hatching can be divided into at least two seasonal groups (spring and autumn) on the basis of hatching date and life history traits of individuals (Ching et al., 2017). This species of squid is important to fisheries because of its high economic value, however, empirical evidence concerning the movement and distribution during the ontogenetic stages is scarce. Ueta and Jo (1990) studied the migration of subadultadult individuals of S. lessoniana around Tokushima Prefecture, Japan, by using the tag-recapture method and fishery data. Their results suggested that this species stays in inshore waters and migrates to offshore waters for overwintering. Kanamaru et al. (2007) evaluated past tagging studies for S. lessoniana migration in Japan and reported that individuals released in autumn moved farther than those released in spring. However, additional empirical evidence to support these hypotheses is required because conventional tagging studies did not provide complete information on movement patterns throughout the life history of squids. In addition, high-resolution data of seasonal and geological movement patterns of S. lessoniana populations are required for species management.
The present study investigated the ontogenetic movement patterns of different seasonal groups of S. lessoniana in northern Taiwan through statolith δ 18 O analyses. This study is the first attempt to generate a 3D habitat preference model, including spatial and depth preferences, and it is expected to provide new comprehensive spatial and temporal insights into migration patterns of cephalopods.

Squid Collection and Age Estimation
Twenty-two adult S. lessoniana individuals were collected by jigging from the inshore waters of northeastern Taiwan (Figure 1) between November 2017 and March 2018. The mantle length (ML; in mm) of each individual was measured. Statoliths were extracted, cleaned ultrasonically using 70% hydrogen peroxide, rinsed, oven-dried, and embedded in Epofix resin (Struers, Denmark). The left statolith of each squid was ground and polished along the posterior side to approximately 50-100 µm above the core by using a metallographic grinding and polishing machine (P20FR-HA; Top Tech Machines co. Ltd., Taiwan). This thickness ensured that a sufficient volume of milling powder was available for isotopic analysis (>40 µg per sample). Alternative formation of translucent and opaque growth zones has been previously validated as occurring on a daily basis (Jackson, 1990); thus, the growth increment at the lateral dome region of statolith was examined from the photographs recorded under a compound microscope (400×, DM-2500, Leica Microsystems GmbH, Germany) by using a digital camera (DFC-450, Leica Microsystems Lt., Switzerland) to estimate age. The growth increments were counted twice, and an average value was used as daily age. If the difference between two counts was >5%, a third count was recorded to minimize measurement error (Arkhipkin and Shcherbich, 2012). The hatching date of each individual was back-calculated from the deduced age (in days) and the date of collection. Each individual was further categorized into the nearest seasonal group according to its hatching month, namely spring (March-May), summer (June-August) and autumn (September-November). Based on the morphological and ecological changes in accordance with growth and age (Segawa, 1987), we defined four stages in S. lessoniana corresponding to age, namely the embryonic-paralarval (age: 0-20 days), juvenile (age: 20-60 days), juvenile-subadult (age: 60-110 days) and subadult-adult (age: >110 days) stages.

Isotopic Analysis
After age counting, statolith slides were drilled using an ESI New Wave Research Micromill and carbonate powders were collected sequentially from the edge to the core at the lateral dome region at intervals of 150-190 µm (Figure 2). The tip of the drill was approximately 200 µm in diameter (H23RS, Comet, Germany), and milling depth was set at approximately 150 µm. Most statoliths involved four drilling paths (except a smaller statolith, K180313003, which had three paths only), and each path was recorded to estimate covered daily growth.
The powder samples of each drilling path were transferred to glass vials and reacted with 100% orthophosphoric acid at 70 • C in an automated online system (Kiel Carbonate IV, Thermo Electron Corporation, Germany) to produce CO 2 . The values of δ 18 O were determined by analyzing the released CO 2 gas by using a mass spectrometer (Finnigan MAT 253, Thermo Electron Corporation, Germany) at National Taiwan University. The long-term reproducibility of the Finnigan MAT 253 is higher than ±0.08% (one standard deviation [SD]) for δ 18 O, based on repeat samples of international reference standards . The values of δ 18 O (%) were reported in standard notation relative to standards Vienna Pee Dee Belemnite (VPDB) after calibration against the NBS-19 standard:  (Chang, 2000): To evaluate the feasibility of the equation, the relationship of measured and deduced temperature derived from the edge of statoliths in 22 crossed-season captured individuals were analyzed to understand if it followed the 1:1 correspondence. The duration (days) which statolith δ 18 O values represented was considered, and the measured temperature was averaged based on the duration at the depth of 50 m.

Prediction of Movement Patterns
To accurately predict the experienced temperature and corresponding living areas of S. lessoniana, individual differences in living period and temporal and spatial variations in seawater temperature were considered in the evaluation. The individual living period and season corresponding to each data point of statolith δ 18 O values were established through the examination of microstructure. During the defined period, water temperature and salinity were obtained from the HYbrid Coordinate Ocean Model website (HYCOM) 1 and are presented in a spatial resolution of 0.08 • × 0.08 • at the depths of 0, 30, 50, and 100 m, according to the living depth of this species (Roper et al., 1984;Tomano et al., 2016;Ammar and Maaroof, 2019). Next, we set the unit of the spatial grid at 0.4 • × 0.4 • to establish the movement pattern because this spatial grid covers the minimum range of squids captured in the adult stage in the present study (Figure 1). For setting units, the average values of temperature and salinity in each grid at each depths (0, 30, 50, and 100 m) was used to produce δ 18 O water,VSMOW values based on the equation established by Chang (2000).
The deduced temperature depended on the estimated statolith δ 18 O value and the variation in water δ 18 O values among grids at different depths. Therefore, each data point of statolith δ 18 O values was used to deduce the experienced temperature in every grid at various depths, each of which had their specific water δ 18 O values. When the deduced temperature matched with the measured temperature (from the HYCOM), the grid thus obtained was considered the possible living area of the squid during a specific period; this period was estimated using the growth rings in the milling area for the statolith δ 18 O measurement. The extent of match between the deduced and measured temperature was based on the comparison of probability distributions between these two values determined using Student's t test (Figure 3). The probability distribution of deduced temperature was modeled using a known living period (days) as well as an average and the uncertainty (SD) associated with the deduced temperature. The uncertainties were estimated by running Monte Carlo simulations 1,000 times and included variations in salinity, instrumental measurements, and parameters of the temperature-dependent equation. Similarly, the probability distribution of measured temperature was modeled using known living periods (days) and average values of measured temperature and associated SDs. Once the t test resulted in a p value larger than 0.05, we accepted that the two temperatures did not differ significantly and inferred a possibility of the occurrence of a living area in the grid (Figure 3). We repeated modeling to determine the matched area and depth of each statolith δ 18 O value from individuals at each life stage.
Furthermore, we considered the collection location, spawning site, and movement ability of S. lessoniana to determine the possible living areas. The δ 18 O value in the outermost portion of statoliths reflected the occurrence of S. lessoniana near the collection location. If this value represented the 10-day average of the signal, the maximum movement of S. lessoniana was calculated to be approximately 50 km away from the collection location because the mean swimming speed of adult S. lessoniana individuals is 5 km per day, according to a study by Kanamaru et al. (2007). Consequently, grids located farther than 50 km from the collection site were eliminated. In addition to this "backward" evaluation, based on the collection location, a "forward" method was also used to evaluate the living area at the embryonicparalarval stage. We excluded the matched areas that were not adjacent to the coast because the egg capsules of S. lessoniana are always found in inshore waters (Segawa, 1987). Combining the results of forward and backward calculations and mobility, the possible occurrence of four life stages of individuals were further selected from the matched areas found based on the statolith δ 18 O value.
After examining the living area at an individual level, we calculated the probability of occurrence in each setting area for FIGURE 3 | Determination of the extent of matching between deduced and measured temperature by using the two-sample t test. The deduced temperature was based on the statolith δ 18 O values, and measured temperature was the data extracted from the HYbrid Coordinate Ocean Model website. The test was conducted in each spatial grid (0.4 • × 0.4 • ) at depths corresponding to each specific life stage of S. lessoniana. The criteria of matching were set at the p value of 0.05.
Frontiers in Marine Science | www.frontiersin.org each seasonal group. The probability was calculated using the following equation: where C kj represents the occurrence of kth squid in area j, and k j is the total number of individuals in grid j. If squid sample number 1 exists in grid 1 at a specific life stage, the C kj value is 1; otherwise, the C kj value is 0. Furthermore, the P j value indicates the probability of occurrence in grid j. If the P j value is 1, all squid samples at the same ontogenetic stage occur in grid j. If an area with a probability (p value) of >0.5 is identified, it would be considered a major residential area for S. lessoniana seasonal groups around Taiwan.

Statistical Analysis
To

RESULTS
After statolith age determination, seven individuals were found to belong to the spring group with a mean (±SD) ML of 264 ± 49 mm and an age ranging from 150 to 167 days. Six individuals were found to belong to the summer group with a mean ML of 260 ± 24 mm and an age ranging from 146 to 192 days. Nine individuals were found to belong to the autumn group with a mean ML of 298 ± 45 mm and an age ranging from 140 to 177 days (  Figure 4A). Larger variations were observed in the summer and autumn groups  Table 2). The values decreased slightly from embryonic-paralarval stage to the juvenile stage and then increased until the subadult-adult stage ( Figure 4B). The values of δ 18 O statolith in embryonic-paralarval stage were significantly lower than those in the subadult-adult stage (Dunn's tests, Z = −4.078, p < 0.001). In addition, the values in the juvenile stage were significantly different from those in the juvenile-subadult and subadult-adult stages. No interaction was observed between the seasonal group and ontogenetic stages in δ 18 O statolith values ( Table 2).
Individual ontogenetic trends showed that the δ 18 O statolith values in the spring group ranged between −2.87 and -1.61 and were relatively stable until the statolith edge (Figures 5A,B). In the summer group, the δ 18 O statolith values in the statolith core varied from −2.54 to −1.87 ; they decreased to the lowest levels, between −2.32 and −2.83 , at the juvenile stage and then increased to the highest levels, between −1.06 and −0.65 , in the statolith edge (Figures 5C,D). The δ 18 O statolith values in the autumn group showed a pattern similar to that of the summer group (Figures 5E-G), but the δ 18 O statolith values were higher at the juvenile-subadult stage (−2.23 to −0.97 ) and subadultadult stage (−0.50 to −0.12 ) than those in the summer group.

Predicted Occurrence Area and Movement Pattern
Individual values of deduced and measured temperature were close and the linear regression showed a well correspondent (Figure 6, p < 0.001). The slope was not significantly different from the line of 1:1 (F 1,34 = 0.83, p = 0.37) indicating that the deduced temperature could reflect reasonably the experienced temperature of S. lessoniana. Figures 7-9 present the distribution of occurrence probability of seasonal groups at different depths by comparing deduced and measured temperatures while considering collection location, spawning site, and movement ability. The probability of occurrence results indicated that the three seasonal groups of S. lessoniana in northeastern Taiwan have diverse distributions and movement patterns. The individuals in the spring group had the highest possibility of hatching in neritic waters (approximate depth: 0-50 m) near the coast, extending from northeastern to eastern Taiwan and the Ryukyu Islands (Figure 7). After hatching and reaching the paralarval stage, individuals were widely distributed in the inshore waters of northeastern Taiwan, possibly from the sea surface to an approximate depth of 50 m; these individuals may have subsequently migrated to relatively deeper waters near the coast of northeastern Taiwan or to the offshore waters of eastern Taiwan. In addition, the waters at a depth of 30-50 m along southern Taiwan were also potential hatchling grounds, and the hatched individuals moved northward as the growth proceeded. Finally, they mainly remained in the northeastern waters during the subadult-adult stages (Figure 7).
The summer group was most likely to hatch in the areas near northeastern Taiwan, extending southward to the Penghu Islands (Figure 8). However, the predicted distribution at the juvenile stage covered a wide area, owing to constant seawater temperatures in summer. Similarly, a wide and deep living area (approximate depth: 100 m), including the inshore waters of northern Taiwan and the waters near the Ryukyu Islands (approximate depth: 100 m), was observed in the last two stages, and it differed from the pattern of the spring group.
In the autumn group, the predicted hatching sites were similar to those in the summer group (Figure 9). However, juvenile individuals might have been distributed in the southern waters (Figure 9), but they were not detected in the waters, with the probability >0.5, because the deduced temperature at this stage did not completely match the water temperature in the study area. The autumn group then spent their subadult and adult stage (approximately 3-4 months) in the inshore waters of northern Taiwan at depths between 0 and 100 m before capture.

DISCUSSION
Stable oxygen isotope ratios recorded in statoliths provide information regarding diverse movement patterns of S. lessoniana; differences in ontogenetic distributions were observed among the three seasonal groups. Although the approach for investigating the movement of animals by using δ 18 O values in biogenic carbonates has been widely applied in studies on fish (Trueman et al., 2012;Currey et al., 2014;Shiao et al., 2017;Darnaude and Hunter, 2018), thus far, studies have not adapted this method to determine the movement history of cephalopods. We demonstrated the potential of using this approach in studies on cephalopod ecology. Because of the significant effects of larval dispersal and demographic population connectivity on cephalopod resources (O'dor, 1992;Boyle and Boletzky, 1996;Semmens et al., 2007), the results of the present study may improve the fishery management and conservation for the bigfin reef squid.

Variation in Statolith δ 18 O Values in Relation to Experienced Temperature
The variation in statolith δ 18 O values among individuals has been observed in several squid species and is associated with differences in experienced temperature (Radtke, 1983;Landman et al., 2004;Trasviña-Carrillo et al., 2018). The general trend of statolith δ 18 O values is to decrease with an increase in temperature; this trend follows the theoretical predictions and observations from other biogenic carbonates (Rexfort and Mutterlose, 2006;Trueman et al., 2012;Kitagawa et al., 2013;Linzmeier et al., 2016). However, assessing the reliability of using statolith δ 18 O values to reconstruct experienced temperature in cephalopods is challenging because previous studies have analyzed δ 18 O value from an entire statolith, which indicates the mean experienced temperature throughout its lifespan. Thus, in our study, we evaluated the deduced temperature compared with the seawater temperature and temperature preference of S. lessoniana at different life stages.
The deduced temperature of the embryonic stage and an approximate 20-day paralarval stage from the statolith core ranged from 22 to 28 • C, regardless of seasonal group. This finding was consistent with those of other studies that found that S. lessoniana hatches in a warm environment of approximately 20 to 30 • C (Segawa, 1987;Walsh et al., 2002;Ikeda et al., 2009). Each squid species has its optimum living temperature, which supports its growth and survival, particularly in its early life stages (Jackson and Choat, 1992;Forsythe et al., 2001). This suggests that the population of S. lessoniana in northern Taiwan shows a preference for this temperature at hatch.
The statolith δ 18 O values exhibited variation after the paralarval stage among seasonal groups because of seawater temperature varying with the seasons. For example, the spring group experienced warmer (summer) temperatures at the juvenile and adult stages than at the hatching stage in spring. By contrast, the summer and autumn groups reached their juvenile and adult stages in autumn and winter, respectively hence, they experienced lower temperatures in their juvenile and adult stages than in their hatching stage. The summer and autumn groups exhibited more positive δ 18 O values than the spring group at the juvenile and adult stages because of lower temperatures. Thus, the ontogenetic variation in statolith δ 18 O values in the spring group was less obvious than in the summer and autumn groups (Figure 4). The statolith δ 18 O values of S. lessoniana ontogenetic variation mirrored the seasonal temperature fluctuation off Taiwan waters, thus indicating the reliability of the prediction method used in the present study.

Prediction of Ontogenetic Movement and Geographical Distribution
The seasonal and vertical migration of squid species has been described in the literatures (Bazzino et al., 2010;Argüelles et al., 2012;Yamaguchi et al., 2019); this behavior was also observed in our study. Interpreting time-series data in squids is generally more difficult than interpreting data obtained from bivalves (e.g., Nakashima et al., 2004;Owen et al., 2008;Nishida et al., 2015) because the squid species move freely, unlike bivalves. We converted the measured statolith δ 18 O values to the deduced temperature and cross-matched them with the seasonal and depth changes in the seawater temperatures to estimate the probability of the geographical distribution of S. lessoniana by using a procedure derived from a widely used method to study fish migration (involving the use of otolith δ 18 O values; Thorrold et al., 1997;Weidel et al., 2007;Shiao et al., 2014). In our case, the movement patterns of S. lessoniana in four life stages exhibited diversity among the seasonal groups (Figures 7-9). These findings are comparable to the known ecology of S. lessoniana, including spawning (Segawa, 1987), migration (Ueta and Jo, 1990;Kanamaru et al., 2007) and depth distribution (Roper et al., 1984;Tomano et al., 2016).
According to statolith δ 18 O values, most probably the individuals of the spring and summer groups hatched near the coasts of northeastern Taiwan, with an occurrence probability of 1; these findings support empirical evidence that the waters near the coast of northeastern Taiwan are one of main spawning grounds for S. lessoniana (Chen et al., 2015;Ching et al., 2017). The dominant topography in northern Taiwan is that of an eroded coastline with complicated topographical features and structures (Song et al., 1997), forming macroalgae-rich and coral-rich environments for spawning. In addition, the quantity of nutrients supplied by the year-round upwelling of the Kuroshio Current off northeastern Taiwan supports high primary production and sustains sequential consumers in the waters (Liu et al., 1992;Chen, 1997;Gong et al., 2003). Therefore, this biomass could be a crucial factor sustaining the abundance of S. lessoniana near Taiwan (e.g., Otero et al., 2008;Rodhouse et al., 2014). However, the individuals of the autumn group had a high probability (>0.8) of hatching in the waters of southern and southwestern Taiwan (i.e., the Penghu Islands). The paralarvae were transported northward by the Kuroshio Current or the South China Warm Current into the Taiwan Strait during September and October (Tang et al., 2000;Jan et al., 2002Jan et al., , 2006. Although the results suggest that waters at a depth of 100 m near the coast of southern Taiwan are also potential hatching grounds, adult females are unlikely to have laid eggs at a depth of down to 100 m because the environment is unsuitable for planktonic paralarvae and no egg capsules have been observed on the seabed at the aforementioned depth. Distribution of juvenile individuals appeared to be widespread around Taiwan, both horizontally and vertically. This distribution pattern has two possible explanations. First, the juvenile squids were passively shifted by ocean currents, thus reflecting the seawater temperatures of a larger region. Second, the consistent seawater temperatures around Taiwan in summer and autumn reduced the precision of distribution prediction. In particular, the estimated probabilities of juvenile distribution in autumn group were all less than 50%. Compared with the juvenile stage, the subadult-adult stages showed a narrower area of predicted distribution. In the spring group, the subadult individuals remained at a depth of approximately 50 m near the coast of northeastern Taiwan or at a depth of approximately 100 m in the inshore waters of eastern Taiwan. Subsequently, adult individuals migrated to the coasts of northeastern for mating and spawning. To our knowledge, larger individuals are rarely captured in the eastern waters, which suggests that the eastern waters are not a principal habitat for S. lessoniana individuals in the spring group. Nevertheless, the occurrence of S. lessoniana in the eastern water needs to be understood through additional surveys using systematic fishery records. By contrast, individuals in the juvenile-subadult stage in the summer group were mostly found in the areas between the coastal waters of China to northeastern Taiwan; furthermore, individuals in this stage in the autumn group might remain in the waters in northeastern Taiwan or around the Penghu Islands. Both the summer and autumn groups used the northern waters as their main habitat during the adult stage in winter. In addition, adult individuals in the autumn group appeared to move to further offshore waters (about 100 km). During winter months, drastic reductions in temperature occur in northern Taiwan primarily during strong northeasterly monsoon and cold surge events (Chen and Huang, 1999;Chen et al., 2002). In such a turbulent state, the individuals move offshore for overwintering and return to inshore areas for feeding and spawning when the environment becomes relatively stable. This explanation is consistent with the findings of Ueta and Jo (1990), who studied the migration of S. lessoniana subadult-adult individuals around Tokushima Prefecture.
The individuals in the seasonal groups in the present study exhibited a possible migration route from the Penghu Islands to northeastern Taiwan and considerable habitat overlap. In general, loliginid squids spawn throughout the year and consequently exhibit multicohort formation along with highly diverse dispersion patterns; hence, high levels of genetic diversity are achieved in the populations of these cephalopods (O'dor, 1998). A study reported that the elemental signatures in the entire statolith exhibited less variation between S. lessoniana samples from northern Taiwan and the Penghu Islands for the same season (Ching et al., 2017). This finding suggests a high level of population connectivity in S. lessoniana in Taiwan. In recent years, the coexistence of three cryptic lineages of S. lessoniana has been reported in the Indo-Pacific Ocean (Cheng et al., 2014;Tomano et al., 2016). These cryptic lineages of S. lessoniana exhibit similar morphology but are genetically distinct (Akasaki et al., 2006;Hsiao et al., 2016;Shen et al., 2016). Although the extent of cryptic diversity within the S. lessoniana species complex in Taiwan remains unclear, a single cryptic lineage may be predominant in northern Taiwan and the Penghu Islands, based on the migration pattern prediction. Additional studies with larger areas of geographic sampling and using a combination of molecular methods are needed to provide more knowledge regarding the population structure of the bigfin reef squid over its distribution range.

Method Improvement and Applications in Future
As the first study to determine the ontogenetic movement of squid by using statolith δ 18 O values, we provided information on the life history of S. lessoniana, but we acknowledge that our method can be improved considerably in the future.
First, the δ 18 O values in the juvenile statoliths of the autumn group suggested high experienced temperatures (approximately 25-29 • C), and this does not satisfactorily accord with the observed water temperature range from the sea surface to a depth of 100 m depth in autumn. We emphasize that the mechanisms of isotopic fractionation in statoliths [such as the results of fish otolith by Thorrold et al. (1997) and Høie et al. (2004a)] and other potential sources of variability (e.g., Høie et al., 2004b;Darnaude and Hunter, 2018;Linzmeier, 2019) should be carefully considered. Second, the seawater δ 18 O values vary by <1% across the surface in the present study area and by approximately 1% with depth (LeGrande and Schmidt, 2006), slightly biasing the prediction of experienced temperature. We used salinity to predict seawater δ 18 O values, based on the equation established from the waters in Taiwan Strait (Chang, 2000). The development of location-specific relationships between salinity and seawater δ 18 O values can minimize the bias of reconstructed temperature. Third, reducing the amount of sample powder from statolith can increase the temporal sampling resolution and enhance the precision of reconstructed temperature (Leder et al., 1996;Høie et al., 2004b). As reported by Sakamoto et al. (2019) in their otolith study, the weight of drilling powder can be as low as 0.3-11.4 µg, representing a temporal resolution of 10-30 days, which is considerably higher than that in our study (>40 µg and approximately 30-50 days). A statolith analysis with higher temporal resolution can significantly benefit the studies on stock discrimination and individual migration. Fourth, defining minimum and maximum probabilities that can indicate the existence of S. lessoniana populations requires additional statistical support. In addition, the size of the spatial grid can be reduced to match the movement behavior of S. lessoniana for a higher precision of habitat determination.
Seawater temperature directly and strongly affects cephalopod ecology and fisheries (Jackson and Moltschaniwskyj, 2002;Forsythe, 2004). Statolith δ 18 O is a crucial parameter that provides evidence of general patterns of distributional extent and movement. Statolith, which are involved in orientation and balance, are found in all cephalopod species; hence, they can be widely used for ecological research on cephalopods (Clarke, 1978;Arkhipkin, 2005). For example, the predicted geographical distributions based on statolith δ 18 O signatures are comparable to the estimated stock boundaries determined using fishery data or tagging methods. Precise geographic boundary and habitat use (e.g., spawning ground) allow managers to implement suitable management measures to conserve targeted species (Gislason et al., 2000;Hobday et al., 2010). In recent years, in response to varaitions in the behavior of oceans and difficulties in managing resources, scientists have highlighted the importance of technological improvements through use of finer spatial and temporal scales for near real-time animal tracking (Maxwell et al., 2015;Dunn et al., 2016). This study combined the deduced temperatures of the individuals over ontogenetic stages to describe continuous movement patterns of S. lessoniana lifespan. The intra-and inter-annual movement patterns also support the decisions pertaining to the establishment of fishing grounds, forecasting of catches, and dynamic fishery management for cephalopods (e.g., Yamaguchi et al., 2019). However, the extent of geographic distribution resolution using statolith δ 18 O may be species specific. Our results showed that accurately predicting residence waters during the periods when variations in water temperature of the region are not obvious (from summer to early autumn in Taiwan) is difficult. Thus, statolith δ 18 O is unlikely to serve as the main indicator of distribution for tropical cephalopod species. Species with a large diel vertical migration (e.g., Young, 1978;Hunt and Seibel, 2000) also probably biases the prediction of experienced temperature, thus increasing the risk of misinterpretation when this method is applied to a species. Establishing a temperature-dependent relationship of δ 18 O for specific cephalopod species is also necessary for improving the approach in the future.

CONCLUSION
All wild squids possess natural isotopic signatures that are incorporated into their statoliths and reflect the environmental temperature or isotopes composition. Migration between waters may produce shifts in isotopic composition. This study applied sequential isotopic sampling to statoliths, and the findings suggest the ecological features of seasonal movement strategies and population connectivity in S. lessoniana in Taiwan. As a result, the waters off northeastern Taiwan have proven to be an important spawning ground for S. lessoniana. Implementing fishery management by restrictions or closed fishing zone during spawning seasons will ensure the recruitment of S. lessoniana population in Taiwan sustainable. Flexible life history traits and a large distribution range are critical for cephalopods and support a high level of genetic diversity and ensure population abundance. These findings extend the limited knowledge about the life history of S. lessoniana within a year. Future developments can reduce the uncertainty associated with this approach and provide more accurate species-specific interpretations of the variations of statolith δ 18 O values within individuals and stocks of freemoving cephalopods.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

ETHICS STATEMENT
Ethical review and approval was not required for the animal study because we collected the squid samples from the commercial fishery and did not conduct a control or treatment experiment.

AUTHOR CONTRIBUTIONS
C-HW led this project. C-HW, T-YC, and AY coordinated the experiment. T-YC and AY assisted in wild squid collections. C-IC wrote the manuscript and prepared the samples. C-IC and M-TC were responsible for data analysis. J-CS and P-LW contributed the isotopic analyses. All authors revised, reviewed, and finally approved the manuscript.

FUNDING
This study was funded by a project from the Ministry of Science and Technology, Taiwan (MOST 106-2611-M-019-004) to C-HW.