Abstract
This study investigated the relationship of the catch rates (CRs) of Spanish mackerel (Scomberomorus commerson) with oceanographic factors in the waters around Taiwan by using high-resolution fishery and environmental data for the period 2011–2016. The investigation results revealed that trammel nets accounted for 69.79% of the total catch of S. commerson and were operated mostly in the Taiwan Strait (TS). We noted seasonal variations in the distribution of high CRs. These CRs were observed in the southwestern TS, including the waters along the southwestern coast of Taiwan and around the Penghu Islands, and extended to the Taiwan Bank during autumn; they increased in winter. To predict the spatial and temporal patterns of Spanish mackerel density and their relationship with oceanographic and spatiotemporal variables, generalized additive models were used. These models explained 48.4% of the total deviance, which was consistent with the assumed Gaussian distribution. Moreover, all variables examined were significant CR predictors (p < 0.05). Latitude and longitude were the key factors influencing the spatiotemporal distribution of S. commerson, and sea surface chlorophyll a concentration was a key oceanographic factor. Observing projected changes in El Niño/Southern Oscillation events for S. commerson revealed that CRs were higher and distributed further southward during La Niña events than during other events. We inferred that the S. commerson distribution gradually moved toward the southwest with the northeast monsoon, which was enhanced during La Niña in winter.
Introduction
The waters around Taiwan (hereafter abbreviated as WT) are affected by several currents, including the China Coastal Current, Kuroshio Current, Kuroshio Branch Current, and South China Sea Surface Current (, ). Moreover, the fluctuation of these currents is strongly influenced by the seasonal monsoons and varied topography; consequently, the currents change the oceanographic conditions in the WT in different seasons (; ). Biogeographically, the Taiwan Strait (TS) and waters off eastern Taiwan are migration routes for numerous pelagic fish species; therefore, their ecosystem dynamics and biogeochemical and physical processes vary substantially in space and time (; ; ). The WT constitute excellent spawning, feeding, and overwintering habitats (; ; ). Overall, these unique features contribute to the abundant and multiple living fishery resources in the coastal waters off Taiwan, rendering these waters a valuable fishing ground (; ).
Scomberomorus species, which include the narrow-barred Spanish mackerel (Scomberomorus commerson), Japanese Spanish mackerel (Scomberomorus niphonus), and spotted Spanish mackerel (Scomberomorus koreanus), are crucial commercial species in the WT (; ). According to the life-history features of Scomberomorus spp., the swimming dispersal of larvae and adults of these species is extensive, ensuring connectivity among stocks and causing genetic homogenization across distant populations (; ; ). Among these species, S. commerson is the most commonly caught, which accounts for 80% of the total production of Scomberomorus species in Taiwan (Figure 1; ). S. commerson is a high-migration species and epipelagic predator that is extensively distributed in the Indo-Pacific region, from shallow coastal waters to the edge of the continental shelf at depths of 10–70 m (; ; ).
FIGURE 1
Trammel nets, longlines, and trolling lines are primarily used for catching S. commerson in the WT (). In recent years, numerous challenges have been countered in the production of Scomberomorus species, and the catch of S. commerson decreased from 6600 metric tons in 2002 to 508 metric tons in 2018 (Figure 1). This dramatic decline might be due to the high demand, unregulated fishing practices, or overexploitation of these species. evaluated the stock status of S. commerson under fishing pressure in the WT and found that three stocks had collapsed. S. commerson species were classified as highly susceptible to the high levels of exploitation in the waters off eastern Taiwan (). Furthermore, climate-driven environmental variability has long been recognized as another factor affecting fishery resources and continues to be critical in the WT (; ; ).
The catch of S. commerson is significantly related to sea surface temperature (SST) changes; the recommended SST and salinity for achieving optimal S. commerson catch rates (CRs) are 14–31°C and 23–35 psu, respectively (; ). Several studies have suggested that chlorophyll-a concentrations and sea surface height (SSH) influence the distribution patterns of pelagic fish species and the variability of pelagic fish species abundance (; ; ). Furthermore, changes in marine environments affect fishery resources at different spatial and temporal scales (; ; ). The most extensively studied climatic events affecting fish include those on interannual and decadal scales, such as the El Niño/Southern Oscillation (ENSO) and Pacific Decadal Oscillation, which can cause SST changes on different time scales in the western Pacific Ocean (; , ; Wu et al., 2020). ENSO events also exert a major influence on the wind stress of monsoons in the WT, and they engender changes in SST and upwelling intensity on an interannual scale (; ). Moreover, habitat changes may exert multiple economic effects on coastal communities through reduced availability of ecosystem services, such as fishery landings and ecotourism ().
Species distribution models (SDMs) have been extensively used to predict the distribution of a species under climate change; they have also been used to assess the past distribution of a species by analyzing evolutionary relationships and the management of the species’ habitat (; ). In addition, computer algorithms based on such models have been applied to predict the spatiotemporal habitat distribution of a species by using environmental data and to explore spatiotemporal variations in species and environments. Generalized additive models (GAMs) are among the most commonly used SDMs (), which are nonparametric extensions of generalized linear models (). GAMs enable the incorporation of smoothing functions to model the nonlinear effect of continuous explanatory variables. They can be used to effectively determine the multiple nonlinear relationships between covariates and response variables in a semiparametric manner (). Moreover, GAMs can determine highly nonlinear and nonmonotonic relationships between responses and sets of explanatory variables (Wood, 2006); accordingly, such models are ideal for expressing underlying relationships in ecological systems. Numerous studies have applied GAMs to predict the spatial distribution of fishing grounds for a single species or the distribution of fisheries’ resources (; ; ; ).
In the literature, data on the relationship between the distribution of S. commerson and environmental variations are limited because high-resolution fishery data are insufficient. To address this gap in the literature, the present study collected high-spatial-resolution catch and effort data to investigate the temporal and spatial variations in fishing grounds and CRs of S. commerson in relation to oceanographic conditions. The study also explored the mechanisms underlying the ENSO events that influence the interannual CRs and spatiotemporal distribution of this species. Estimating the preferred environmental conditions of crucial commercial species is critical for determining their migratory patterns, and it is an essential step toward the ecosystem-based management of fisheries.
Materials and Methods
Spanish Mackerel Fishery Data
Daily high-resolution fishery data for S. commerson were collected from the voyage data recorders and logbooks of fishing vessels operating in the WT from 2011 to 2016. The fishery data comprised daily operating positions for 0.1° spatial grids, including latitude and longitude information, fishing dates, fishing methods, working hours, and total catch (kg) of S. commerson.
To determine the fishery dynamics and spatiotemporal variation of S. commerson in the WT, this study calculated the annual catch percentage of each fishing method. According to the data, S. commerson are mainly caught using trammel nets; therefore, the monthly observed CRs of trammel nets were calculated for individual 0.1° spatial grids across the study region by using the following equation:
where i represents the latitude and j represents the longitude of each 0.1° spatial grid.
The spatial distributions of catch percentages associated with the various fishing methods and seasonal CRs were mapped using Quantum GIS version 3.6 (). Moreover, the longitudinal and latitudinal gravitational centers of the observed CRs (G) were estimated using monthly longitudinal and latitudinal locations of the fishing vessels (L) and monthly observed CRs () as follows:
where i and j denote the latitude and longitude, respectively.
Environmental Data
Remote sensing data of SST, sea surface chlorophyll-a (CHL), sea surface salinity (SSS), and SSH in the WT were collected. Satellite-derived SST data were extracted from National Oceanic and Atmospheric Administration (NOAA) Advanced Very-High-Resolution Radiometer (AVHRR) SST images with a spatial resolution of 1.1 km. The NOAA High-Resolution Picture Transmission data, including the AVHRR scenes, were received at the ground station of National Taiwan Ocean University. Daily CHL data were downloaded with a spatial resolution of 1.1 km from the National Aeronautics and Space Administration’s Ocean Color Web (). The SSS and SSH data were downloaded with a 0.08° spatial resolution from the Asia Pacific Data Research Center.1 Both types of data were obtained from the HYbrid Coordinate Ocean Model (HYCOM). To compensate for missing data due to cloud coverage, this study constructed monthly averaged composite maps of the remotely sensed environmental data to fit the fishery data by using Interactive Data Language (version 7.0).
Statistical Models for Spatial and Temporal Predictions of Catch Rates
To analyze the relationships between environmental variation and CR, statistical models were used to predict the spatiotemporal pattern of Spanish mackerel density. SDMs were developed, and on the basis of these models, GAMs were developed to examine potential seasonal fishing grounds. The GAMs were constructed through R (version 3.6; ) by using the GAM function of the “mgcv” package. The observed CRs constituted the response variable, and spatiotemporal factors (year, month, longitude, and latitude) and environmental factors (SST, SSH, SSS, and CHL) constituted the predictor variables. The GAMs can be expressed as follows:
where a0 is a constant and s(xi) is a spline smoothing function for each model covariate xn or the interaction between two covariates. All covariates were considered to be continuous, and the effective degrees of freedom were estimated for each main factor. Because the log-link function cannot handle zeros, we added a constant value of 0.1 (c) to all CRs; this value is commonly used in CR standardization processes (e.g., ; , ). Time and location were treated as interaction terms to account for the possible interannual variability engendered by environmental spatial distribution variations.
The model with the optimal conformation was selected using a stepwise procedure that was based on the lowest value of Akaike’s information criterion (AIC), and the p-value for the final set of variables was lower than 0.05. The selected GAMs were used to predict the relative abundance of S. commerson in the WT from 2011 to 2016. Furthermore, this study applied the Oceanic Niño Index (ONI) as the primary indicator for monitoring El Niño and La Niña alternating phases and examined how the climate variability in ENSO events affects the variation of S. commerson. The ONI was downloaded from the Climate Prediction Center of the National Weather Service (National Oceanic and Atmosphere Administration).
Results
Spatiotemporal Distribution of Scomberomorus commerson in the Waters Around Taiwan
The spatial distribution of catch percentages of S. commerson caught using different fishing methods revealed that trammel nets and longlines accounted for a wider distribution among all fisheries (Figure 2). Longlines were mostly operated off eastern Taiwan, whereas trammel nets were mostly operated in the TS. Of the annual total catch of S. commerson, trammel nets accounted for 69.79 ± 11.98% and longlines accounted for only 12.71 ± 4.48% (Figure 3).
FIGURE 2
FIGURE 3
Accordingly, this study analyzed the influence of trammel nets on the distributions of and relationship between CR and environmental factors. Concerning the seasonal mean spatial distributions of CRs associated with trammel nets, the results indicated that higher CRs were mainly distributed around the southwestern waters of the TS during autumn; these CRs increased during winter but decreased gradually during spring and summer (Figure 4).
FIGURE 4
Time-Series Variation of Catch Rates and Fishing Locations of Trammel Nets
The time series of the monthly catches and trammel net CRs for S. commerson during 2011–2016 revealed that the most catches and highest CRs occurred in autumn and winter (Figure 5). The highest catch and CR occurred in January 2013 (nearly 130 metric tons and 28.64 kg/h, respectively). The lowest catch and CR occurred in July 2011 (0.43 metric tons) and May 2011 (0.91 kg/h), respectively. These results indicate that the major fishing season for S. commerson in the WT begins in early autumn and ends in winter. Trammel fisheries are distributed around the WT and are used year round, with some seasonal variation.
FIGURE 5
The monthly latitudinal and longitudinal G of the CRs revealed seasonal variations, and the latitudinal and longitudinal G gradually moved from northern and western areas (23.7°N, 119.81°E) in summer to southern and eastern areas (23.35°N, 119.66°E) in winter (Figure 6). By contrast, the latitudinal and longitudinal G of the CRs exhibited no obvious variations from 2014 to 2015.
FIGURE 6
Environmental Effect on the Catch Rates of Scomberomorus commerson
Our analyses included a total of 8,859 data points regarding S. commerson catch. Concerning the environmental variables affecting the CR of S. commerson, the selected GAMs explained 48.4% of the deviance. Normal quantile–quantile plots (Figure 7A) indicated that the distributions of the residuals for the selected GAMs adequately conformed to the assumed Gaussian distribution.
FIGURE 7
The addition of predictor variables at different levels increased the explained deviance, which was attributed to the decreased AIC value (Table 1). The results revealed that the oceanographic and spatiotemporal variables significantly affected the CRs of S. commerson (p < 0.01). The interaction of year and month with latitude (24.2 and 26.2%, respectively) and that of longitude with latitude (29.8%) explained a large proportion of the variance in CR. High CRs were observed in the southwest part of the TS (Figure 8C; 22°–25°N, 118–120°E) during the major fishing season: early autumn to winter (Figure 8B). Analyzing annual variation data indicated that the high CRs in the south of the TS (Figure 8A; 22°N–25°N) decreased from 2013 to 2015.
TABLE 1
| One variable model explained (%) | AIC | p-Value | ||
| Log (CR+c) = | s(SST) | 5.82 | 16,944.38 | <0.01 |
| +s(SSH) | 2.58 | 16,729.14 | <0.01 | |
| +s(SSS) | 2.32 | 16,541.91 | <0.01 | |
| +s(CHL) | 6.01 | 16,112.57 | <0.01 | |
| +s(YY, LAT) | 24.2 | 14,211.19 | <0.01 | |
| +s(MM, LAT) | 26.2 | 13,505.94 | <0.01 | |
| +s(LAT, LON) | 29.8 | 12,609.22 | <0.01 | |
| Total deviance explained | 43.6 | |||
Individual and total deviance explained and Akaike’s information criterion value in the model with the optimal conformation, as determined using a stepwise procedure.
FIGURE 8
Regarding the oceanographic variables, CHL explained the largest proportion of the variance in CR, followed by SST, SSH, and SSS (Table 1). The results also revealed that high CRs were positively associated with SSTs of approximately 20–25°C, CHL concentrations of 1–2 mg/m3, SSHs of approximately 0.6–0.7 m, and SSSs of >34.0 psu (Figures 7B–E). SST and CHL were the main variables determining the distribution of S. commerson; specifically, the abundance of S. commerson decreased as the SST and CHL concentration increased.
Projected Changes in El Niño/Southern Oscillation Events on the Fishing Grounds
The GAM results regarding annual variations revealed that the high CRs in the south of the TS (Figure 8A) decreased from 2013 to 2015 during El Niño events (Figure 9B). A box plot of the CRs in winter during various ENSO events indicated the same pattern (Figure 9A). The CR was observed to be 27.01 ± 37.6 kg/h during La Niña events (2011 and 2012), and it decreased to 17.61 ± 30.08 kg/h during El Niño events (2015 and 2016). Therefore, the CRs differed significantly between the various ENSO events (p < 0.05). To examine the effect of ENSO events on the CRs and distribution of S. commerson, this study used the optimal model to predict CRs during winter for the 2011–2016 period. This study also compared the spatial distributions of the predicted CRs for S. commerson during the El Niño events in January and February 2016 with those during the La Niña events in January and February 2012, and the results are presented in Figure 10. During the La Niña events, high CRs (>7 kg/h) were distributed southward, which extended south of 22°N. By contrast, during the El Niño events, high CRs were distributed mostly toward the north of 22°N.
FIGURE 9
FIGURE 10
Discussion
Using daily catch data for S. commerson along with vessel data recorder data, this study investigated the spatiotemporal variations of fishing grounds in relation to oceanographic conditions in the WT. Trammel nets are operated mostly in the TS, whereas longlines are operated off eastern Taiwan. Moreover, revealed that S. commerson are caught using drift gill nets and trolling lines in the Taiwan Bank. Trammel nets are composed of single, double, or triple layers of netting kept vertical by floats on the headrope and weights on the groundrope (). According to the results of a bathymetric survey in the WT, the average water depths in eastern areas exceed 2,000 m, but those in western areas are less than 50 m that is, the water depth off eastern Taiwan is obviously more than that of western Taiwan (). Therefore, the use of trammel nets in the waters off western Taiwan is highly efficient. Through the analysis of vessel operational data and bathymetric data, this study demonstrated that various fishing methods are used to effectively catch S. commerson in the WT. Although the data on trammel nets were only used to determine the relationship between the CRs of S. commerson and the oceanographic factors, the commercial catch data for S. commerson included data on various fishing methods (gears). The differences in effort among the considered fishing methods may explain the variations of the calculated CRs; however, it may incorrectly attributed varying catchability to differences in the abundance of fishery resources (; ). described an estimation model that accounts for variations in catchability and can thus aid in identifying unidimensional indices to represent variables that reflect aspects of an underlying domain. The effort data (working hours) used in this study were collected from voyage data recorders. Thus, subsequent studies should eliminate bias due to confounding among the various methods (gears).
This study revealed seasonal variations in the spatial distribution of S. commerson, which was distributed from northern and western areas in summer to southern and eastern areas in winter. speculated that S. commerson may move southeast from the north in the TS in winter. The present study also indicated that the distribution patterns of S. commerson changed with seasonal progression. The WT are strongly influenced by the East Asian monsoon; therefore, wind stress explains a majority of the transport reversals (; ). The mixed China coastal water, with its low temperature and low salinity, is conveyed by the northeast monsoon–driven flow to the northern TS in winter (; ). In this study, the GAM results reveal relatively high CRs at SSTs of 20–25°C during autumn and winter (the major fishing season). Thus, we suggest that S. commerson migrate abundantly into the TS from early winter to late spring by following the mixed China coastal water; subsequently, they migrate northward and southward as the SST gradually increases.
During autumn, high CRs were distributed widely across the southwestern TS, including in the waters along the southwestern coast of Taiwan and between the Penghu Islands and the Taiwan Bank, and concentrated gradually in winter. Previous studies on upwelling in the TS have reported that upwelling is a year-round phenomenon, with the degree of upwelling varying in the water around the Penghu Islands and Taiwan Bank (; ; ). Fronts and mesoscale eddies are formed in the Taiwan Bank through the interaction of different water masses and monsoons, which all influence larval species due to the dispersion and convergence of nutrients (Zhang et al., 2014; ). revealed that the fishery resource structures in the frontal habitat of the Taiwan Bank are supported by both warm- and cold-water species, including S. commerson. reported the presence of hydrated and postovulatory oocytes for S. commerson in the central TS, implying that this area could be a spawning ground for this species. Thus, the aggregative response of biological interactions links primary productivity (PP) and top predators could be feeding and spawning grounds for S. commerson among different seasons in the upwelling of southwestern TS.
The GAM results revealed that CHL and SST were the dominant environmental factors affecting the CRs. The spawning season of S. commerson was estimated to be from March through August in the TS, and the SST might influence the species’ gametogenesis, gonad atresia, and spawning behavior (; ). The area between the Taiwan Bank and the Penghu Islands is a major upwelling region and key fishing ground for S. commerson. The thermal fronts in the Taiwan Bank are characterized by two frontal belts on the northern and southern edges. A previous study indicated that the formation of the thermal front along the southern edge is linked to the interaction between topography and upwelling phenomena and that the thermal front along the northern edge may be tidally generated (). In the present study, we determined that physical oceanographic elements such as SSH can influence the degree of productivity of S. commerson habitats. SSH reflects oceanic features such as current dynamics, fronts, eddies, and convergences (). Our study demonstrated that the highest CRs corresponded to the areas with SSTs of 20–25°C, CHL concentrations of 1–2 mg/m3, SSHs of 0.6–0.7 m, and SSSs of 34.0–34.5 psu. These areas are potential habitats for S. commerson and their prey because of their physical oceanographic structure. Nutrients are supplied by the year-round upwelling of the southwestern TS, and the CHL pattern in the frontal region varies with the supply of nutrient-rich water from subsurface upwelling (). Our findings are similar to those of a study on S. commerson in the northern Persian Gulf and Java Sea, which observed that S. commerson is sensitive to SSTs ranging from 13 to 29°C, SSSs ranging from 23 to 35 psu, and low CHL concentrations (; ; ).
Numerous studies have identified a positive correlation between the abundance of marine species and CHL. However, we found that the abundance of S. commerson decreased as the CHL concentration increased. Previous studies have also identified a negative correlation between abundance and CHL. revealed that the catch per unit effort of yellowfin tuna (Thunnus albacares) increased as the CHL concentration decreased. reported that the likelihood of a nonzero blue marlin (Makaira nigricans) catch increased with the CHL concentration but that the CRs decreased. found that Bull sharks (Carcharhinus leucas) were present when CHL with low level. CHL concentrations in early stage prior to the present abundance index possibly influenced the CRs, and this time-lag phenomenon could be related to differences in trophic level species transformation (). In addition, PP is an essential measure of an ocean’s capacity to transform carbon dioxide into particulate organic carbon at the base of the food web and is an effective predictor of the potential yield of the world’s oceans. Studies (, ) have suggested that PP is a key factor in the variation in biological resources at each trophic level. We used CHL as a phytoplankton biomass index rather than PP. Although CHL exhibited similar features to PP in the study area, CHL data derived from satellites should be analyzed by reparametrizing the original relationships on the basis of in situ data (). Therefore, PP could be incorporated into models to explore the relationship between the CRs of S. commerson and oceanographic factors.
Generalized additive models could aid in exploring the possible responses of a highly migratory species such as S. commerson to oceanographic variables. Spatial and temporal variables are essential to GAMs because they can be used to determine whether changes in the CRs of a species are related to environmental variables, spatiotemporal variables, or interactive variables. Our results suggest that CR assessment models that incorporate spatiotemporal and environmental variables can be used to identify the possible effects of such variables on the distribution of S. commerson. In addition, climate variability may lead to variations in the distribution of S. commerson. Climate-driven environmental variabilities of ENSO are pivotal factors affecting the intensity of upwelling regions in the southwestern TS (Wu et al., 2020; ). Variations in fishery resource structures and the spatial distribution of pelagic species are also influenced by ENSO events in the Taiwan Bank and change the upwelling size during summer (). Moreover, the northeast monsoons in the TS are stronger during La Niña winters than during other seasons, increasing the intensity of variations in hydrographic features relative to those observed during El Niño events (; Wu et al., 2020). The monthly spatial and temporal predictions for S. commerson during winter revealed that the CRs were higher and distributed further south during La Niña events. This study compared the 20°C isotherms obtained for the winters of 2012 (La Niña) and 2016 (El Niño). As displayed in Figure 11, the 20°C isotherms were located further north during the El Niño winter and extended south of the TS during the La Niña winter. We inferred that S. commerson was concentrated in the southwest when the northeast monsoon intensified during the La Niña winter and that the S. commerson population dispersed toward the northeast when the northeast monsoon weakened.
FIGURE 11
The Gulf Cooperation Council, a regional intergovernmental union, developed several regulations for S. commerson fisheries, such as license requirements, bans on fishing during certain periods, mesh size limitations, landing size requirements, and restrictions on time at sea. Numerous studies have investigated S. commerson in this region and proposed suggestions for the management of resources (; ; ). However, information regarding the population dynamics, age, growth, and reproduction of the fish in our study area is required to formulate regulations for individual or regional management. In addition, approaches to managing fishery resources must account for factors that influence the spatial distribution of the fish. Oceanographic variations associated with regional and global phenomena, such as ENSO events, result in large-scale shifts in the dynamics and distribution of highly migratory species. Therefore, comprehensively reviewing the results of related studies and investigating unexplored topics are crucial for the formulation of effective management strategies.
Conclusion
This study combined daily high-resolution fishery data with environmental data by using GAMs to preliminarily analyze the relationship between the CRs of S. commerson and environmental factors. Using these data, we applied the models to effectively predict the distribution of S. commerson in the WT under the effects of ENSO events. Analyzing fishery dynamics indicated that trammel nets accounted for most of the catch of S. commerson in the WT. Trammel nets are operated mostly in the TS, whereas longlines are operated off eastern Taiwan. The overall CRs of S. commerson in the WT increased from early autumn to the end of winter relative to those in the other seasons. Relatively high CRs were concentrated in the waters off southwestern Taiwan, including the waters along the southwestern coast of Taiwan and between the Penghu Islands and the Taiwan Bank. All environmental factors significantly influenced the CRs of S. commerson and exhibited seasonal variations. CHL and SST were the most influential environmental factors in the WT and could be related to the spawning behavior of this species.
Examining the projected changes in ENSO events for S. commerson in our study revealed that the fishing conditions and fishing grounds were affected during different phases. However, comprehensively determining the fishery dynamics would require long-term and systematically recorded catch data, which can further reveal the degree of influence of variations in environmental characteristics under climate change on specific species. Moreover, interspecific relationships must be considered because such relationships may affect the fishing conditions of single species. The movements and distributions of top predators are governed by bottom-up processes and direct environmental effects (e.g., physiological tolerance to anoxia and thermal preferences; ). Top marine predators can have high phenotypic plasticity and adaptive capabilities that mitigate the effects of climate change on them; however, climate change may still affect these predators through their prey (). Hence, understanding how low-trophic marine ecosystems are related to each life stage of S. commerson can help researchers clearly understand the mechanism through which climate change influences the distribution and abundance of S. commerson. Future research should thus apply additional long-term oceanographic variables and prey species in different prediction models, considering the influence of interspecific relationships.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this article are not readily available because due to legal restrictions imposed by the government of Taiwan in relation to the “Personal Information Protection Act,” data cannot be made publicly available. Requests for data can be sent as a formal proposal to the FA, COA (http://www.fa.gov.tw/). Requests to access the datasets should be directed to http://www.fa.gov.tw/.
Author contributions
L-CC: formal analysis, data curation, and writing – original draft preparation. K-WL: conceptualization, methodology, writing – review and editing, and supervision. P-YH and MN: software and visualization. J-SW and C-TT: data curation and supervision. C-CC: formal analysis. All authors have read and approved the published version of the manuscript.
Funding
This study was financially supported by Taiwan’s Council of Agriculture (109AS-1.1.5-ST-aC) and the National Science Council (MOST 108-2611-M-019-007 and MOST 109-2611-M-019-007).
Acknowledgments
We give special thanks to Director General J.-R. Chen and Chief H.-J. Hsieh of the Fisheries Research Institute for their support and guidance. We are also grateful to P.-H. Huang and Y.-L. Wu of National Taiwan Ocean University for their help with this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
References
1
AschR. G.CheungW. W. L.ReygondeauG. (2018). Future marine ecosystem drivers, biodiversity, and fisheries maximum catch potential in Pacific Island countries and territories under climate change.Mar. Policy88285–294. 10.1016/j.marpol.2017.08.015
2
BellJ. D.GanachaudA.GehrkeP. C.GriffithsS. P.HobdayA. J.Hoegh-GuldbergO.et al (2013). Mixed responses of tropical Pacific fisheries and aquaculture to climate change.Nat. Climate Change3591–599. 10.1038/NCLIMATE1838
3
BishopJ. (2006). Standardizing fishery-dependent catch and effort data in complex fisheries with technology change.Rev. Fish Biol. Fish.1621–38. 10.1007/s11160-006-0004-9
4
ChangY.ShimadaT.LeeM. A.LuH. J.SakaidaF.KawamuraH. (2006). Wintertime sea surface temperature fronts in the Taiwan Strait. Geophys. Res. Lett.33:L23603. 10.1029/2006GL027415
5
ChassotE.BonhommeauS.DulvyN. K.MélinF.WatsonR.GascuelD.et al (2010). Global marine primary production constrains fisheries catches.Ecol. Lett.13195–505. 10.1111/j.1461-0248.2010.01443.x
6
ChassotE.MélinF.WatsonLe PapeO.GascuelD. (2007). Bottom-up control regulates fisheries production at the scale of eco-regions in European seas.Mar. Ecol. Prog. Series34345–505. 10.3354/meps06919
7
ColletteB. B. (2001). “Scombridae. tunas (also, albacore, bonitos, mackerels, seerfishes, and wahoo),” in FAO Species Identification Guide for Fishery Purposes. The Living Marine Resources of the Western Central Pacific. Bony Fishes Part 4 (Labridae to Latimeriidae), Estuarine Crocodiles, edsCarpenterK. E.NiemV. (Rome: FAO), 3721–3756.
8
ColletteB. B.RussoJ. L. (1979). “An introduction to the Spanish mackerels Genus scomberomorus,” in Proceedings of the Colloquium on the Spanish and king mackerel resources of the Gulf of Mexico, (Ocean Springs, MS: Gulf States Marine Fisheries Commission),
9
DudleyR.AghanashinikarA. P.BrothersE. B. (1992). Management of the Indo-Pacific Spanish mackerel (Scomberomorus commerson) in Oman.Fish. Res.1517–43. 10.1016/0165-7836(92)90003-C
10
Fisheries Agency (2019). Taiwan Fisheries Yearbook. Available online at: https://www.fa.gov.tw/cht/PublicationsFishYear/index.aspx
11
HarahapS. A.SyamsuddinM. L.PurbaN. P. (2020). Range of sea surface temperature and chlorophyll-α values based on mackerel catches in the northern waters of West Java, Indonesia.AACL Bioflux132265–2272.
12
HastieT. J.TibshiraniR. J. (1990). Generalized Additive Models.London: Chapman and Hall.
13
HazenE. L.JorgensenS.RykaczewskiR. R.BogradS. J.FoleyD. G.JonsenI. D.et al (2013). Predicted habitat shifts of Pacific top predators in a changing climate.Nat. Climate Change3234–238. 10.1038/nclimate1686
14
HoC. H.LurH. S.YaoM. H.LiaoF. C.LinY. T.YagiN.et al (2018). The impact on food security and future adaptation under climate variation: a case study of Taiwan’s agriculture and fisheries.Mitigation Adaptation Strategies Global Change23311–347. 10.1007/s11027-017-9742-3
15
HongH. S.ChaiF.ZhangC. Y.HuangB. Q.JiangY. W.HuJ. Y. (2011). An overview of physical and biogeochemical processes and ecosystem dynamics in the Taiwan Strait.Continental Shelf Res.313–12. 10.1016/j.csr.2011.02.002
16
HsiaoP. Y.ShimadaT.LanK. W.LeeM. A.LiaoC. H. (2021). Assessing summertime primary production required in changed marine environments in upwelling ecosystems around the Taiwan Bank.Remote Sensing13:765. 10.3390/rs13040765
17
HsuC. C.HanY. S.TzengW. N. (2007). Evidence of flathead mullet Mugil cephalus L. spawning in Waters Northeast of Taiwan.Zool. Stud.46717–725.
18
JanS.SheuD. D.KuoH. M. (2006). Water mass and throughflow transport variability in the Taiwan Strait.J. Geophys. Research-Oceans111:C12012. 10.1029/2006JC003656
19
JanS.WangJ.ChernC. S.ChaoS. Y. (2002). Seasonal variation of the circulation in the Taiwan Strait.J. Mar. Systems35249–268.
20
JuP. L.TianY. J.ChenM. R.YangS. Y.LiuY.XingQ. W.et al (2020). Evaluating stock status of 16 commercial fish species in the coastal and offshore waters of Taiwan using the CMSY and BSM methods.Front. Mar. Sci.7:618. 10.3389/fmars.2020.00618
21
KarakulakF. S.ErkH. (2008). Gill net and trammel net selectivity in the northern Aegean Sea. Turkey.Sci. Mar.72527–540. 10.3989/scimar.2008.72n3527
22
KuoN. J.HoC. R. (2004). ENSO effect on the sea surface wind and sea surface temperature in the Taiwan Strait.Geophys. Res. Lett.31:L13309. 10.1029/2004GL020303
23
LamT. J. (1983). “Environmental influences on gonadal activity in fish,” in Fish Physiology, edsHoarW. S.RandallD. J.DonaldsonE. M. (New York, NY: Academic Press, Inc).
24
LanK. W.ChouC. P.LeeM. A.VayghanA. H. (2018). Association between the interannual variation in the oceanic environment and catch rates of bigeye tuna (Thunnus obesus) in the Atlantic ocean.Fish. Oceanogr.27, 1–13. 10.1111/fog.12259
25
LanK. W.KawamuraH.LeeM. A.ChangY.ChanJ. W.LiaoC. H. (2009). Summertime sea surface temperature fronts associated with upwelling around the Taiwan Bank.Continental Shelf Res.29903–910. 10.1016/j.csr.2009.01.015
26
LanK. W.NishidaT.LeeM. A.LuH. J.HuangH. W.ChangS. K.et al (2012). Influence of the marine environment variability on the yellowfin tuna (Thunnus albacares) catch rate by the Taiwanese longline fishery in the Arabian Sea, with special reference to the high catch in 2004.J. Mar. Sci. Technol.20514–524. 10.6119/JMST-011-0506-1
27
LanK. W.LeeM. A.ZhangC. I.WangP. Y.WuL. J.LeeK. T. (2014). Effects of climate variability and climate change on the fishing conditions for grey mullet (Mugil cephalus L.) in the Taiwan Strait.Climatic Change126189–202. 10.1007/s10584-014-1208-y
28
LanK. W.LianL. J.LiC. H.HsiaoP. Y.ChengS. Y. (2020b). Validation of a primary production algorithm of vertically generalized production model derived from multi-satellite data around the waters of Taiwan.Remote Sensing12:1627. 10.3390/rs12101627
29
LanK. W.ChangY. J.WuY. L. (2020a). Influence of oceanographic and climatic variability on the catch rate of yellowfin tuna (Thunnus albacares) cohorts in the Indian Ocean.Deep-Sea Res. Part II: Top. Stud. Oceanography175:104681. 10.1016/j.dsr2.2019.104681
30
LanK. W.WuY. L.ChenL. C.NaimullahM.LinT. H. (2021). Effects of climate change in marine ecosystems based on the spatiotemporal age structure of top predators: a case study of bigeye tuna in the Pacific Ocean.Front. Mar. Sci.8:614594. 10.3389/fmars.2021.614594
31
LanK. W.ZhangC. I.KangH. J.WuL. J.LianL. J. (2017). Impact of fishing exploitation and climate change on the grey mullet (Mugil cephalus L.) stock in the Taiwan Strait.Mar. Coastal Fish.9271–281. 10.1080/19425120.2017.1317680
32
LeeK. A.SmootheyA. F.HarcourtR. G.RoughanM.ButcherP. A.PeddemorsV. M. (2019). Environmental drivers of abundance and residency of a large migratory shark, Carcharhinus leucas, inshore of a dynamic western boundary current.Mar. Ecol. Prog. Series622121–137. 10.3354/meps13052
33
LehodeyP.BertignacM.HamptonJ.LewisA.PicautJ. (1997). El Niño Southern Oscillation and tuna in the western Pacific.Nature389715–718. 10.1038/39575
34
LiaoC. H.LanK. W.HoH. Y.WangK. Y.WuY. L. (2018). Variation in the catch rate and distribution of swordtip squid (Uroteuthis edulis) associated with factors of the oceanic environment in the southern East China.Mar. Coastal Fish.10452–464. 10.1002/mcf2.10039
35
LinC. Y.WangS. P.ChiangW. C.GriffithsS.YehH. M. (2020). Ecological risk assessment of species impacted by fisheries in waters off eastern Taiwan.Fish. Manag. Ecol.27345–356. 10.1111/fme.12417
36
LiuC. S.LiuS. Y.LallemandS. E.LundbergN.ReedD. L. (1998). Digital elevation model offshores Taiwan and its tectonic implications.Terrestrial Atmospheric Oceanic Sci.9705–738.
37
MadiganD. J.ChiangW. C.WallsgroveN. J.PoppB. N.KitagawaT.ChoyC. A.et al (2016). Intrinsic tracers reveal recent foraging ecology of giant Pacific bluefin tuna at their primary spawning grounds.Mar. Ecol. Prog. Series553253–266. 10.3354/meps11782
38
ManuderM. N.PuntA. E. (2004). Standardizing catch and effort data: a review of recent approaches.Fish. Res.70141–159. 10.1016/j.fishres.2004.08.002
39
McPhersonG. R. (1985). Development of the northern Queensland mackerel fishery.Australian Fish.4415–17.
40
MeriemS. B.Al-MarzouqiA.Al-MamryJ. (2006). Fisheries exploitation pattern of narrow-barred Spanish mackerel, Scomberomorus commerson, in Oman and potential management options.J. Appl. Lchthyol.22218–224. 10.1111/j.1439-0426.2006.00739.x
41
National Aeronautics and Space Administration [NASA] (2017). Moderate-Resolution Imaging Spectroradiometer (MODIS) Aqua Data.Maryland, MD: NASA.
42
NelderJ. A.WedderburnR. W. N. (1972). Generalized linear models.J. R. Statist. Soc. A135370–384.
43
NguyenK. Q.NguyenaV. Y. (2017). Changing of sea surface temperature affects catch of spanish mackerel Scomberomorus Commerson in the set-net fishery.Fish. Aquaculture J.8:1000231. 10.4172/2150-3508.1000231
44
NiamaimandiN.KaymaramF.HoolihanJ. P.MohammadiG. H.FatemiS. M. R. (2015). Population dynamics parameters of narrow-barred Spanish mackerel, Scomberomorus commerson (Lacèpéde, 1800), from commercial catch in the northern Persian Gulf.Global Ecol. Conserv.4666–672. 10.1016/j.gecco.2015.10.012
45
PuntA. E.WalkerT. I.TaylorB. L.PribacF. (2000). Standardization of catch and effort data in a spatially-structured shark fishery.Fish. Res.45129–145.
46
QGIS Development Team (2019). QGIS 3.6 Noosa. 2019. Available online at: https://qgis.org/en/site(accessed 30 September, 2019).
47
R Core Team (2018). R: A Language and Environment for Statistical Computing.Vienna: R Core Team.
48
RandallJ. E. (1995). Coastal fish of Oman.Honolulu: University of Hawaii Press.
49
Roa-UretaR. H.LinY. J.RabaouiL.Al-AbdulkaderK.QurbanM. A. (2019). Life history traits of the narrow-barred Spanish mackerel (Scomberemorus commerson) across jurisdictions of the southeast Arabian Peninsula: implications for regional management policies.Regional Stud. Mar. Sci.31:100797. 10.1016/j.rsma.2019.100797
50
ShojiJ.TanakaM. (2005). Distribution, feeding condition, and growth of Japanese Spanish mackerel (Scomberomorus niphonius) larvae in the Seto Inland Sea.Fish. Bull. Natl. Oceanic Atmospheric Administration103371–379.
51
ShuiB. N.HanZ. Q.GaoT. X.MiaoZ. Q.TakashiY. (2009). Mitochondrial DNA variation in the East China Sea and Yellow Sea populations of Japanese Spanish mackerel Scomberomorus niphonius.Fish. Sci.75593–600. 10.1007/s12562-009-0083-3
52
SolankiH. U.BhatpuriaD.ChauhanP. (2017). Applications of generalized additive model (GAM) to satellite-derived variables and fishery data for prediction of fishery resources distributions in the Arabian sea. Geocarto Int.32, 30–43. 10.1080/10106049.2015.1120357
53
SuN. J.SunC. L.PuntA. E.YehS. Z.DiNardoG. (2011). Modelling the impacts of environmental variation on the distribution of blue marlin, Makaira nigricans, in the Pacific Ocean.ICES J. Mar. Sci.681072–1080. 10.1093/icesjms/fsr028
54
TsengH. C.YouW. L.HuangW.ChungC. C.TsaiA. Y.ChenT. Y.et al (2020). Seasonal variations of marine environment and primary production in the Taiwan Strait.Front. Mar. Sci.7:38. 10.3389/fmars.2020.00038
55
TsengW. Y.ChenC. H.HuS. H.ChenG. S.PaoW. H. (1971). Preliminary study on Spanish mackerel of Taiwan.Bull. Taiwan Fish. Res. Institute1889–114.
56
WangT. W.ChanT. Y.ChanB. K. K. (2013). Diversity and community structure of decapod crustaceans at hydrothermal vents and nearby deep-water fishing grounds Off Kueishan Island, Taiwan: a high biodiversity deep-sea area in the NW pacific.Bull. Mar. Sci.89505–528. 10.5343/bms.2012.1036
57
WangY. F.YaoL. J.ChenP. M.YuJingWuQ. (2020). Environmental influence on the spatiotemporal variability of fishing grounds in the Beibu Gulf, South China Sea.J. Mar. Sci. Eng.8:957. 10.3390/jmse8120957
58
WengJ. S.WuS. F.LoY. S.ShiaoJ. C.LeeM. A.LiuK. M.et al (2020). Reproductive biology of the narrow-barred Spanish mackerel (Scomberomorus commerson) in the central Taiwan Strait, western Pacific.Deep-Sea Res. Part II: Top. Stud. Oceanography175:104755. 10.1016/j.dsr2.2020.104755
59
WoodS. M. (2006). Generalized Additive Models, an Introduction with R.London: Chapman and Hall.
60
WuY. L.LeeM. A.ChenL. C.ChanJ. W.LanK. W. (2020). Evaluating a Suitable Aquaculture site selection model for cobia (Rachycentron canadum) during extreme events in the inner Bay of the Penghu Islands. Taiwan.Remote Sensing12:2689. 10.3390/rs12172689
61
ZhangF.LiX.HuJ.SunZ.ZhuJ.ChenZ. (2014). Summertime sea surface temperature and salinity fronts in the southern Taiwan Strait.Int. J. Remote Sensing354452–4466. 10.1080/01431161.2014.916454
Summary
Keywords
catch rate, Spanish mackerel, generalized additive model, fishing dynamic, environmental factor, climatic variability, waters around Taiwan
Citation
Chen L-C, Weng J-S, Naimullah M, Hsiao P-Y, Tseng C-T, Lan K-W and Chuang C-C (2021) Distribution and Catch Rate Characteristics of Narrow-Barred Spanish Mackerel (Scomberomorus commerson) in Relation to Oceanographic Factors in the Waters Around Taiwan. Front. Mar. Sci. 8:770722. doi: 10.3389/fmars.2021.770722
Received
04 September 2021
Accepted
29 November 2021
Published
15 December 2021
Volume
8 - 2021
Edited by
Inna Senina, Collecte Localisation Satellites (CLS), France
Reviewed by
Guillaume Briand, Collecte Localisation Satellites (CLS), France; Anindya Wirasatriya, Diponegoro University, Indonesia
Updates
Copyright
© 2021 Chen, Weng, Naimullah, Hsiao, Tseng, Lan and Chuang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kuo-Wei Lan, kwlan@mail.ntou.edu.tw
This article was submitted to Marine Fisheries, Aquaculture and Living Resources, a section of the journal Frontiers in Marine Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.