Spatiotemporal variations and influential factors affecting reproductive dynamics of Pacific saury in the Northwest Pacific during fishing season.

This study examined the spatial and monthly variations in reproductive dynamics of Pacific saury collected in the high sea (40°N-49°N, 149°E-168°E) of the northwestern Pacific during the fishing season of the Taiwanese stick-held dip-net fisheries (July-November) from 2018 to 2019. Generalized additive mixed-effects models (GAMMs) were applied to explore the relationship between the probability of maturing and spawning occurrence (PMOS ) and the explanatory variables (month, body length, sea-surface temperature [SST] as fixed effects, and fishing set as a random effect) for the age 0 and age 1 fish, respectively. In addition, the impact of the geographical difference in growth rates of age 0 fish, quantified as the radius of the otolith annual ring, on the PMOS of age 1 fish was explored in the GAMMs. Results showed that the mean values of the condition factor (CF) varied by months, and higher mean CFs were observed north of 45°N. However, the mean values of the gonado-somatic index did not show a clear spatiotemporal pattern. The spatial distribution of the ovarian maturation states revealed that the spawning ground of Pacific saury during July-November was broader than that previously known and had partly overlapped with the high-sea fishing ground. The best GAMM showed that the predicted PMOS of the age 0 fish were generally low (mean = 19%, standard deviation [SD] = 10%) and tended to increase with increasing body length. In contrast, the predicted PMOS were relatively higher (mean = 37%, SD = 10%) for the age 1 fish and increased with increasing SST and decreasing latitude. Two candidate GAMMs, which are equally supported as the best model (ΔAIC < 2), provide evidence that a higher growth rate of age 0 fish may result in higher reproductive activity in age 1 fish. This is one of few studies focusing on the reproductive dynamics of Pacific saury during the fishing season; current limitations, future directions, and conservation implications were discussed.

Evaluation of cooking effects on otolith stable carbon and oxygen isotope values of teleostean fish Pomadasys kaakan (Cuvier, 1830).

RATIONALE: For years, archaeologists, climatologists, and ecologists have used stable oxygen isotope values (δ13 C, δ18 O) in fish otoliths from archaeological sites to reconstruct the habitats, paleo-temperature, and seasonality of the fish captures. Otoliths from archaeological sites might have been heated when ancient people cooked the fish for food. Therefore, there are debates as to whether the cooking behaviors would cause further isotopic fractionations of the carbonate in the otoliths. METHODS: In this study, we have evaluated the effects of the cooking methods on the otolith δ13 C and δ18 O values by comparing the otoliths of the javelin grunter (Pomadasys kaakan) from the same individuals, with the left otoliths taken out before the different cooking processes. Otolith sections of the fish were then made and several subsamples were milled along the microstructures visible in the otolith pairs, mostly annual check rings, followed by the stable isotope analyses. RESULTS: There were no morphological changes between the cooked and uncooked otoliths. The δ13 C and δ18 O values were highly consistent for the otolith subsamples between the cooked and uncooked pairs, suggesting none or trivial effects of the cooking processes on the isotopic values of the otoliths. In addition, some javelin grunters showed lower δ13 C (-5‰ to -6‰) and δ18 O (-4‰ to -5‰) values deposited in the wide translucent zone of the otoliths, suggesting seasonal migration of the fish to estuaries during the summer. CONCLUSIONS: This study suggests that cooking processes do not change otolith stable isotopic compositions and will therefore allow for future research to use the otolith δ13 C and δ18 O values to better understand the life history and used habitats of the preserved fish remains in the middens.

Mercury bioaccumulation and trophic transfer in Pacific Saury from the North Pacific Ocean.

This study investigates mercury (Hg) dynamics in Pacific Saury (Cololabis saira) across the North Pacific Ocean, specifically off East Japan in 2018. Saury traits vary with total mercury (THg) concentrations in muscle tissues ranging from 0.017 to 0.082 µg g-1 w. w., averaging of 0.042 (n = 46). A positive correlation between THg and saury length (Knob length, 270-319 mm) indicates increased Hg concentration with size. Stable isotopic tracers suggest Pacific Euphausiids (Krill) are significant contributors to the saury diet (>70% of total). Significant correlations between logarithm THg concentration (Log THg) and δ15N (‰) (R2 = 0.70) demonstrate Hg trophic biomagnification, with regional variations. Comparative analysis between the eastern (ENPO) and western North Pacific Ocean (WNPO) indicates differences, with WNPO saury exhibiting lower δ15N values and higher THg levels than ENPO saury. This suggests that the WNPO, located near East Asia, the world's largest Hg emitter, experiences elevated Hg levels in seawater due to anthropogenic release. Overall, this study advances understanding of Pacific Saury's ecological interactions and Hg bioaccumulations, emphasizing the importance of species-specific behaviors and regional influences in ecological studies.

Inferring the ecology of north-Pacific albacore tuna from catch-and-effort data.

Catch-and-effort data are among the primary sources of information for assessing the status of terrestrial wildlife and fish. In fishery science, elaborate stock-assessment models are fitted to such data in order to estimate fish-population sizes and guide management decisions. Given the importance of catch-and-effort data, we scoured a comprehensive dataset pertaining to albacore tuna (Thunnus alalunga) in the north Pacific Ocean for novel ecological information content about this commercially valuable species. Specifically, we used unsupervised learning based on finite mixture modelling to reveal that the north Pacific albacore-tuna stock can be divided into four pseudo-cohorts. We discovered that smaller body mass pseudo-cohorts inhabit relatively high-subtropical to temperate-latitudes, with hotspots off the coast of Japan. Larger body mass pseudo-cohorts inhabit lower-tropical to subtropical-latitudes, with hotspots in the western and central north Pacific. These results offer evidence that albacore tuna prefer different habitats depending on their body mass, and point to long-term migratory routes for the species that the current tagging technology is unlikely to capture in full. We discuss the implications of the results for data-driven modelling of albacore tuna in the north Pacific, as well as the management of the north Pacific albacore-tuna fishery.

Effects of decadal climate variability on spatiotemporal distribution of Indo-Pacific yellowfin tuna population.

Spatial variations in tuna population and abundance are strongly linked to large-scale climate fluctuations, such as the Pacific decadal oscillation (PDO) and Atlantic multidecadal oscillation (AMO). However, the mechanisms underlying the association of climate indices with yellowfin tuna (YFT) abundance and habitat preference remain unclear. We analysed long-term longline fishery data for YFT and oceanic climate variability index data for 1971-2018. The standardized catch per unit effort (CPUE) of Indo-Pacific Ocean YFT was higher during negative AMO and positive PDO phases. In tropical Pacific Ocean, the trend of YFT habitat preference exhibited seesaw patterns because of the distinct environmental factors influenced by the PDO phase. The PDO changed the environmental parameters throughout the tropical Indian Ocean such that the habitat preference of YFT remained consistent throughout. However, the variations in habitat suitability did not correspond to the distribution or standardized CPUE of YFT throughout the Pacific Ocean during AMO events. Moreover, the changes in habitat suitability had a positive periodicity of 8-16 years with AMO in the Indian Ocean, but revealed opposite trends with the distribution or standardized CPUE of YFT. Our results provide sufficient information to distinguish the variations between PDO phase changing and YFT standardized CPUE/ habitat preference. Furthermore, the AMO phase shift period 60-100 years longer than that of the PDO (20-30 years), and models employing time series of fishery and environmental data must be extended the time period of our study to make the AMO match the fishery data more complete.