The Role of Environmental Factors on the Fishery Catch of the Squid Uroteuthis chinensis in the Pearl River Estuary, China

The Pearl River Estuary (PRE) is one of the major fishing grounds for the squid Uroteuthis chinensis. Taking that into consideration, this study analyzes the environmental effects on the spatiotemporal variability of U. chinensis in the PRE, on the basis of the Generalized Additive Model (GAM) and Clustering Fishing Tactics (CFT), using satellite and in situ observations. Results show that 63.1% of the total variation in U. chinensis Catch Per Unit Effort (CPUE) in the PRE could be explained by looking into outside factors. The most important one was the interaction of sea surface temperature (SST) and month, with a contribution of 26.7%, followed by the interaction effect of depth and month, fishermen’s fishing tactics, sea surface salinity (SSS), chlorophyll a concentration (Chl a), and year, with contributions of 12.8%, 8.5%, 7.7%, 4.0%, and 3.1%, respectively. In summary, U. chinensis in the PRE was mainly distributed over areas with an SST of 22–29 °C, SSS of 32.5–34‰, Chl a of 0–0.3 mg × m−3, and water depth of 40–140 m. The distribution of U. chinensis in the PRE was affected by the western Guangdong coastal current, distribution of marine primary productivity, and variation of habitat conditions. Lower stock of U. chinensis in the PRE was connected with La Niña in 2008.


Introduction
Uroteuthis chinensis (Gray, 1849) (Cephalopoda: Loginidae), a species of squid, lives in warm continental shelf waters and is widely distributed in the South China Sea (SCS), East China Sea, and Japan [1]. It is a fast-growing and highly reproductive species with a short life span (less than 7 months) and high yield (accounting for about 3/5 of the total production of the family Loliginidae) [2,3]. Because of these characteristics, the U. chinensis is considered to be an ecological opportunist that can increase the population rapidly under a suitable environment [4]. For this reason, environmental factors play a critical role in the life cycle of U. chinensis. Studies showed that cephalopods are sensitive to water temperature [5][6][7][8][9], marine primary productivity [10], and food supply [11], with temperature being the key factor affecting the population biomass and the species' distribution [12,13]. Such conditions could impact the population dynamics by acting on the spawning activity and recruitment [7,8,14]. In the east of the Ionian Sea, the population structure and distribution of Illex coindetii depend on temperature, salinity, and circulation [15]. Any change in water temperature, chlorophyll a concentration, and salinity would affect the catch of Ommastrephes bartramii in the northwestern Pacific to a great extent [16,17]. The available studies on the U. chinensis in the SCS are all focused on biological characteristics [18][19][20], migration the catch of Ommastrephes bartramii in the northwestern Pacific to a great extent [16,17]. The available studies on the U. chinensis in the SCS are all focused on biological characteristics [18][19][20], migration characteristics [21], feeding behavior [22], and resource status [23]. Studies showed that U. chinensis did not migrate on a large scale, but moved in short distances according to local water temperature. It moved northward from the SCS to the Taiwan bank and other places, with the increase of water temperature every spring, and moved southward to the SCS, looking for suitable conditions during winter [21,24]. These studies helped to understand more about the migration characters of the species, but its spatiotemporal distribution and quantitative relationship with marine environment in the Pearl River Estuary (PRE) remains unclear. It is still necessary to understand the impact of environmental variability on U. chinensis abundance in the PRE.
Due to the monsoon, the PRE has two major currents: the Guangdong Coastal Current (GCC) and the South China Sea Warm Current (SCSWC), and two coastal upwellings (in the eastern Hainan Island and the western Guangdong waters) [25][26][27]. In addition, a large amount of fresh water is discharged from the PRE every year, mixing with seawater to form the Pearl River plume [28], which radiates to the coastal waters [29]. Because of these dynamic characteristics, the PRE boasts a high primary productivity [22,30], making it one of the major fishing grounds for U. chinensis in the world [31]. Therefore, the status of U. chinensis stock in the PRE is of great research value. This study is based on six years of survey and satellite remote sensing data and looks into the relationship between the spatiotemporal distribution of the U. chinensis and some environmental factors by conducting a quantitative analysis. The possible mechanism driving the spatiotemporal distribution of U. chinensis in the PRE is also discussed. Results of this research are helpful to understand the migration characteristics, to predict the center of the fishing ground, and to protect the key habitat of the species.

Fishery Data
The research area considered in this research is located between 19.15-22.15° N and 111.12-115.37° E (as shown in Figure 1). The U. chinensis data in the PRE was obtained from the monitoring records of light falling-net fishing vessels from August 2005 to May 2010, and the fishery data was collected at a spatial resolution of 0.25° × 0.25° and summarized in days. The dataset is composed by 359 catch data in total and 20 zero catch observations. The statistics covered longitude and latitude, operating date, voyage, catch species and catches (Table 1). Since 1999, China has suspended fishing activities in the SCS from 1 June to 31 July. In 2009, the closed fishing season was extended from May 16 to August 1. During these periods, ships are prohibited from fishing in the SCS [32,33] and this is why data from June to July are missing in this study.    2005  8  A1, A2, A3, A4, A5  25  10  6  9  A5, A6, A7, A8  14  6  3   2006  8  B1, B2, B3, B4, B5  25  13  6  9  B5, B6  8  4  4  1  E15  16  6  5  2  E17, E18, E19  14  7  4  3  E20  6  2  5 Note: Number of catch species did not include unidentified species.

Satellite Remote Sensing Dataremote Sensing Data
Satellite remote sensing data includes sea surface temperature (SST), sea surface salinity (SSS), sea surface chlorophyll a concentration (Chl a), current, and net primary productivity (NPP). Among them, SST was obtained from the MODIS Aqua Level 3 data products (https://oceancolor.gsfc. nasa.gov/), with temporal resolution set at 1 d and spatial resolution set at 4 km. The SSS and current data could be collected with the help of Global Ocean Physical Reanalysis Product of the Copernicus Marine Environment Management Service (CMEMS, http://marine.copernicus.eu/), with the same temporal resolution and spatial resolution of 0.083 • × 0.083 • . The NPP and Chl a were obtained from the CMEMS (http://marine.copernicus.eu/), with a spatial resolution of 0.25 • × 0.25 • and also under the same temporal conditions. The software R v.4.0.0 (R Development Core Team, 2020) was used to spatially fuse and match SST, SSS, Chl a, and fishery data, and ArcGIS 10.3 (Esri, Redlands, CA, USA) was used to plot the distribution of the flow field and NPP. Data fusion for the remote sensing data with different resolutions could be derived through the following algorithm [34]: where, Ave j is the average value of each environmental factor in the research area after data fusion, with a resolution of 0.25 • × 0.25 • ; m is the number of pixels of each environmental factor in the area with a resolution of 0.25 • × 0.25 • ; the value(i) is the unit pixel value in the study area, and j represents the fishing area, all with the same aforementioned spatial resolution of 0.25 • × 0.25 • .

Catch Per Unit Effort
The fish stock is expressed by catch per unit effort (CPUE), calculated with: where, ∑ Catch (units: tons [t]) is the sum of catches within a fishing grid of 0.25 • × 0.25 • . ∑ Fishing_days is the sum of fishing days for the fishing vessel in the fishing grid, and the unit of CPUE is marked as t day −1 . Day was chosen as the time step in grouping CPUE for each fishing grid.

Clustering Fishing Tactics
Clustering Fishing Tactics (CFT) is a method widely used to identify fishing tactics from commercial data by capturing component records through clustering techniques [35][36][37]. A data matrix containing CPUE records of each species was constructed in the research area, and then the data was normalized to the relative proportion according to weight to eliminate the impact of catch rate. Next, the square root transformation of the standardized matrix was carried out so that the species with lower dominance could make similar contributions to the catch composition. After the initial steps, a principal component analysis (PCA) was applied to the multidimensional data matrix. It is worth noting that most of the changes in the data were explained by the first few axes of the PCA-converted data. Finally, all PCA axes were reserved for clustering analysis [36,38,39]. In this study, CFT identified the fishing tactics adopted by fishermen, so as to evaluate the impact of their behavior and other human factors in fishing activities on the stock of U. chinensis. R 4.0.0 (R Development Core Team, 2020) was used to conduct PCA and cluster analysis on fishery data, and the fishing tactic (FT) result taken as one of the explanation factors of GAM [38]. The fishing strategy was converted into a numerical variable for processing in order to calculate its impact on the change of CPUE.

GAM Analysis
A Generalized Additive Model (GAM) is a nonparametric extension of the generalized linear model [40]. In fishery, GAM is widely used in quantitative analysis of the relationship between fishery resources and environmental factors, which is a non-linear regression [13,41,42]. The distribution of Tweedie was determined in 1984 [43], and it is suitable to describe nonnegative and biased random variables. The process includes several common important distributions (p = 0 normal distribution, p = 1 Poisson distribution, 1 < p < 2 Poisson gamma composite distribution, p = 2 gamma distribution, p = 3 inverse Gaussian distribution) [44]. If 1 < p < 2, it is considered that Tweedie distribution is suitable for analyzing CPUE [35,45]. This specific method is a special probability distribution part of the exponential distribution family, generally expressed by Tw P (θ, ϕ) and completely determined by the variance function V(µ) = µ P , where θ is a gauge parameter and ϕ is a decentralized parameter [43,44]. Since Tweedie distribution still belongs to the exponential distribution family, the statistical relationship model between the corresponding variables and the influencing factors can be established within the framework of the generalized linear model [46]. To avoid the zero-catch problem of CPUE, the GAM was applied together with the Tweedie distribution to analyze the environmental influence on the U. chinensis stock [45]. The GAM basic model constructed in this study is: where, µ CPUE is the mean of Tweedie distribution of CPUE; s is the function of smoothing splines; s(Month, SST) represents the interaction effect of month and SST; s(Month, Depth) is the interaction effect of month and water depth. Furthermore, s(Chl a) stands for the effect of Chl a; s(SSS) is the effect of SSS; s(FT) is the effect of fishing tactic; s(year) is the effect of year. Finally, ε represents the model error. The GAM was constructed and tested by the 'mgcv' package of software R v.4.0.0 [47]. A forward stepwise method was used to select variables that have significant influence on the model, so the specific expression of GAM could be determined.

Analysis of Fishing Tactics
On the basis of the cluster analysis, the catch matrix was identified in five groups, assigned to five different fishing tactics. Among them, FT1 and FT3 mainly included offshore pelagic economic fishes. FT1 was represented by Decapterus (77.08%) and FT3 represented by U. chinensis (62.02%). FT2 mainly included middle and demersal fishes, represented by Trichiurus haumela (75.57%). Furthermore, FT4 mainly considers unidentified species (68.44%), and had a high catch of U. chinensis (30.98%), without a bycatch of Decapterus and Scombridae. Last, the majority of FT5 included oceanodromous fishes, represented by Scombridae (49.04%), mainly composed of Katsuwonus pelamis ( Table 2).

GAM Analysis
The cumulative bias explanatory rate of spatiotemporal and environmental factors on U. chinensis CPUE was obtained by using the GAM to fit and predict the effect of adding variables to the model ( Table 2). The explanatory variables selected for the model included the interaction effect of month and sea surface temperature (Month, SST), the interaction effect of month and water depth (Month, Depth), chlorophyll a concentration (Chl a), sea surface salinity (SSS), fishing tactic (FT), and year (Year). The cumulative bias explanatory rate of these factors for U. chinensis CPUE was 63.1%, and the correlation coefficient R 2 was 0.51 (Table 3).
In GAMs, the contribution of selected factors to CPUE represents the influencing degree of each factor on U. chinensis CPUE (Table 4). Among them, the interaction effect of month and SST was the most influencing factor, with a contribution of 26.7%; followed by the interaction effects of month and depth, FT, SSS, Chl a, and year, with contribution of 12.8%, 8.5%, 7.7%, 4.0%, and 3.1%, respectively. The Chi-square and F test showed that the selected explanatory variables were significantly correlated with CPUE (p < 0.05; Table 4).
The interaction effect of month and SST based on GAM showed that from December to February in the waters with SST of 21-22 • C and from July to August in the waters with SST of 28-30 • C, month and SST had a positive effect on U. chinensis CPUE, which showed an upward trend in this range (Figure 2a). The result of the interaction effect of month and depth indicates that the positive effect of month and water depth on U. chinensis CPUE was the highest in the area of 40-60 m from January to June, when U. chinensis CPUE increased (Figure 2b). In other months (from July to December), no obvious tendency towards depth was shown. The species was evenly distributed in the water depth of 40-130 m. the selected explanatory variables were significantly correlated with CPUE (p < 0.05; Table  4). The interaction effect of month and SST based on GAM showed that from December to February in the waters with SST of 21-22 °C and from July to August in the waters with SST of 28-30 °C, month and SST had a positive effect on U. chinensis CPUE, which showed an upward trend in this range (Figure 2a). The result of the interaction effect of month and depth indicates that the positive effect of month and water depth on U. chinensis CPUE was the highest in the area of 40-60 m from January to June, when U. chinensis CPUE increased ( Figure 2b). In other months (from July to December), no obvious tendency towards depth was shown. The species was evenly distributed in the water depth of 40-130 m. The influence of such factors as FT, SSS, Chl a, and year on the U. chinensis CPUE was analyzed on the basis of GAMs (Figure 3). Results showed that in FT1-FT3, the U. chinensis CPUE increased with the increase of FT; in FT3-FT5, the U. chinensis CPUE decreased with the increase of FT. FT3 had a positive effect on U. chinensis CPUE (Figure 3a). When the SSS was in 30.0-33.0‰, the U. chinensis CPUE increased together with the, SSS but when the SSS was in 33-34‰, it got diminished with the increase of SSS. It was also observed The influence of such factors as FT, SSS, Chl a, and year on the U. chinensis CPUE was analyzed on the basis of GAMs (Figure 3). Results showed that in FT1-FT3, the U. chinensis CPUE increased with the increase of FT; in FT3-FT5, the U. chinensis CPUE decreased with the increase of FT. FT3 had a positive effect on U. chinensis CPUE (Figure 3a). When the SSS was in 30.0-33.0‰, the U. chinensis CPUE increased together with the, SSS but when the SSS was in 33-34‰, it got diminished with the increase of SSS. It was also observed that SSS in 33.0‰ had a positive effect on U. chinensis CPUE (Figure 3b). When Chl a was 0-0.5 mg × m −3 , the U. chinensis CPUE decreased as the levels of Chl a enhanced. Differently, when Chl a was 0.5-0.8 × mg m −3 , the U. chinensis CPUE increased in accordance with Chl a. Another possibility that was observed involved Chl a in the range of 0.8-2 mg × m −3 , and the U. chinensis CPUE decreasing with it. Furthermore, Chl a of 0-0.2 mg × m −3 had a positive effect on U. chinensis CPUE (as shown in Figure 3c). with Chl a. Another possibility that was observed involved Chl a in the range of 0 × m −3 , and the U. chinensis CPUE decreasing with it. Furthermore, Chl a of 0-0.2 m had a positive effect on U. chinensis CPUE (as shown in Figure 3c). From 2005 to 2 U. chinensis CPUE decreased as the years passed. Oppositely, from 2008 to 2010 chinensis CPUE increased with the year. Being more specific, the year of 2008 had tive effect on U. chinensis CPUE (Figure 3d).

Seasonal Variation of Uroteuthis ChinensisU. chinensis
In summer (from August to September in that zone), along the outside of the S the SST was higher and Chl a was lower. Furthermore, SSS in the area were distrib a ladder shape along the PRE to the SCS (Figure 4e,g,i). The fishing grounds of U. c in the PRE were mainly distributed in the sea area at a longitude of 111-114° E and of 19.  (Figures 4a,c,e,g,i).
In winter (from December to February), SSS gradually decreased from east along the GCC, together with Chl a that also got slowly reduced along the coast of dong to the open sea. Major U. chinensis concentrated in the area outside the GC higher SST and lower Chl a (Figure 4f,h,j). The fishing grounds of U. chinensis in were mainly distributed in the sea area at a longitude of 111-115. 5 (Figures 4b,d,f,h,j).

Seasonal Variation of U. chinensis
In summer (from August to September in that zone), along the outside of the SCSWC, the SST was higher and Chl a was lower. Furthermore, SSS in the area were distributed in a ladder shape along the PRE to the SCS (Figure 4e (Figure 4a,c,e,g,i).
In winter (from December to February), SSS gradually decreased from east to west along the GCC, together with Chl a that also got slowly reduced along the coast of Guangdong to the open sea. Major U. chinensis concentrated in the area outside the GCC with higher SST and lower Chl a (Figure 4f,h,j). The fishing grounds of U. chinensis in the PRE were mainly distributed in the sea area at a longitude of 111-115.

Relationship between the U. chinensis and Marine Environment
The fact that water temperature affected each stage of the life history of squid from embryos, juveniles to adults [7,14,48], was one of the important factors indicating its spatiotemporal distribution [12,41]. It was verified that temperature variation could significantly change the growth and development of the reproductive systems of squid [49,50].

Relationship between the U. chinensis and Marine Environment
The fact that water temperature affected each stage of the life history of squid from embryos, juveniles to adults [7,14,48], was one of the important factors indicating its spatiotemporal distribution [12,41]. It was verified that temperature variation could significantly change the growth and development of the reproductive systems of squid [49,50].
Warmer temperatures may lead to the rapid development of embryos and juveniles, shorten their hatching period [48], and thus affect the recruitment of the squid population [7]. In addition, the water temperature affected fish stock in some specific months [50,51]. Therefore, the interaction effect of month and temperature was considered for GAM fitting. Results showed that the interaction effect of SST and month had the greatest impact on the U. chinensis CPUE (with a contribution of 26.7%; Table 4). SST played a positive role on the U. chinensis CPUE in summer (from August to September; Figure 3a). This happened partly because warm water speeded up the embryonic development and the growth process of juvenile squid [14]. In addition, due to the short life span of the squid, its abundance directly depended on the recruitment capacity of the current year [52]. The CPUE increased if the number of individuals hatched in spring and summer-under the appropriate water temperature from August to September. In the East China Sea, a higher SST in summer, especially at 29 • C, could improve the growth rate and reproductive capacity of U. edulis [50]. In winter (from December to February), SST showed a positive effect on the U. chinensis CPUE (Figure 2a). This might be related to the thermophilic-feeding migration of the species, hatched in the autumn, considered to be the spawning season. In the fishing grounds in the south of the Taiwan Strait, the U. chinensis migrated to warmer areas for feeding every winter [53]. According to its age analysis caught by fishing vessels [2,54], generally, U. chinensis will become the target size of the fishing vessel and be captured when they grow to be 2-5 months old. Therefore, from December to February, SST 21-23 • C can be one of the indicators for the distribution of U. chinensis in the PRE.
It was also observed that the levels of salinity affected the physiological activities of marine organisms by changing their osmotic pressure [55]. The GAM analysis showed that SSS contributes 7.7% to the U. chinensis stock ( Table 4). The species was mainly distributed in the waters with a SSS of 32.5-34‰ (Figure 3b). The narrow suitable salinity range for U. chinensis is mainly a result of the poor osmolality regulation mechanism of most cephalopods, since it is a stenohaline species and tends to prefer medium or high salinity environment [31]. Moreover, the salt levels also affected the survival rate of squid embryos, consequently affecting the squid stock as well [56]. When the salinity was lower than 32‰ or higher than 38‰, the survival rate of squid embryos was greatly reduced. It is also noteworthy that large fluctuations in salinity may cause the eggs to die [57].
When it comes to the chlorophyll a concentration, it reflected the phytoplankton in stock in the sea area. As the major feeding source of zooplankton and some marine organisms, phytoplankton is an important part of marine primary productivity and reflects the primary productivity level of the sea area [42,58]. The GAM analysis indicated that the contribution of Chl a to the U. chinensis stock was 4.0% (Table 4). The trophic level of the U. chinensis was 2.7-3.6 [18], as it mainly fed on fish and cephalopods [59]. In the PRE, U. chinensis tended to live in waters with low Chl a (Figure 3c), and the U. chinensis CPUE decreased, while the Chl a did the opposite and kept increasing. This might be related to the hysteresis in the response of the species to chlorophyll a [12]. In addition, the U. chinensis vision was affected by the transparency of sea water. In the southeast waters of Brazil, the abundance of U. chinensis was negatively correlated with the chlorophyll a concentration, as the sea area with a low concentration of the pigment had a high transparency, causing U. chinensis to be easier to be attracted by baits. The consequence of all of this was an increase in the catch [8].

Relationship between U. chinensis and CFT
The GAM analysis also showed that the contribution of fishing tactics to the U. chinensis CPUE was 8.5% (Table 4). Considering that the northern SCS is a typical sea area with multi-species fishes, and with a variety of resource species and complex composition [60,61], fishermen might target species on the basis of economic benefits and the condition of the sea they are acting on. The major species caught in the PRE were U. chinensis (FT3), Euthynnus alletteratus (FT5), Decapterus (FT1), and hairtail (FT2). The examination of these five clustering fishing tactics showed that four clustering fishing tactics (FT1, FT2, FT3, FT5) were aimed at a single species (Table 2). This might be due to the spatial segregation of habitats of major commercial species in the PRE. The spatial segregation of habitat usually led to a negative correlation between the CPUE of the bycatch species and target species [62]. Studies on habitat utilization and feeding behavior of four different species of squid (U. duvaucelii, U. edulis, U. chinensis and Loliolus uyii) in the northern SCS showed that although the feeding strategy was similar, their habitats did not overlap. This spatial segregation could reduce their competition for resources and buffer their trophic interactions, thus increasing the possibility of the species coexisting in the region [22].

Spatiotemporal Variation of U. chinensis
Water depth also affected the spawning activity of squid by changing the level of dissolved oxygen in water [63]. Here, it is important to consider that there is an appropriate season when squid migrate to shallow waters for spawning [31,41,64]. The GAM analysis showed that the interaction effect of month and depth contributes 12.8% to the U. chinensis CPUE (Table 4) and that the species was widely distributed in the water depth of 40-140 m. From January to March, in waters that are 40-60 m deep, water depth played a positive role on the U. chinensis CPUE (Figure 2b). This phenomenon might be related to the mass spawning behavior of the U. chinensis. In the coastal waters of Guangdong Province, the period between February to May was the peak spawning period of the U. chinensis as they migrate to the shallow sea to spawn [31]. This seasonal spawning activity led to a regional fishing season [53]. U. chinensis can spawn in any season, but the peak happens when there is a suitable and appropriate environment for the process to occur [2,65]. In the PRE, the SST during January and March ranged from 20-24 • C, meeting the spawning requirements of the species [31]. This also resulted in an increase of its CPUE in waters from 40-60 m depth and within the period from January to March (Figure 2a,b). From August to November, the water in the PRE became warmer, boosting the metabolism of the U. chinensis, and driving their demand for food, which led to increased feeding activity. Moreover, warmer temperatures provided a more suitable condition for the early lives of U. chinensis [14,48], partly leading to a more even distribution of the animals along the study area (Figure 2b).
The GAM analysis also provided that the contribution of interannual variation to the fluctuation of U. chinensis CPUE was 3.1%. The CPUE decreased significantly from 2007 to 2008 (Figure 3d). It was also observed that large-scale climate variability, such as the North Atlantic Oscillation, affects the abundance of a myriad squid species (e.g., Illex illecebrosus, Loligo pealeii) [66]. For instance, the yield of Loligo Opalesens decreased after the El Nino-Southern Oscillation event occurred in the Southern California Bay [67]. La Niña also affected the recruitment of squid by changing the environment of spawning grounds, see also [16]. In Antarctica, squid got more abundant in the sea area, due to a change in the water temperature caused by La Niña [68]. In the offshore area of the sea area, this resulted in less squid abundance [69]. From August 2007 to April 2008, a La Niña event [70] was formed, which led to the decrease of SST in the SCS. Furthermore, the center of the Symlectoteuthis oualaniensis fishing ground in the Xisha-Zhongsha waters shifted to the south by about 2 • N [13]. Furthermore, La Niña also made rain happen more often and also affected the wind speed in the region [70], as well as the wind speed that strengthened the wave field in the area in terms of cycle and frequency. This new characteristic of the wave field disturbed the near shore spawning grounds, deviating mature fish swarm from these areas, imposing a negative impact on the fish stock [71]. Therefore, it can be said that the decrease of U. chinensis stock in 2008 may be connected with La Niña.

Analysis of Seasonal Variation of U. chinensis
Current and wind fields are also considered important factors in the life span of cephalopods [72] that affect the variation of water temperature and salinity [73,74], and may further change the composition, structure, and thermophilic characteristics of fish communities [75]. The fronts and vortices formed by the flow field served as key spawning habitats for pelagic fishes [76]. Fishing grounds with abundant resources can be formed around the large-scale current, such as the Illex argentinus in the southwest Atlantic Ocean, the Ommastrephes bartramii in the North Pacific Ocean, and the Todarodes pacificus in the waters around Japan, all of them distributed in the western boundary current [1,77]. The current velocity was an important factor for the success of fish migration and feeding [78]. Studies showed that in summer, driven by the southwest monsoon [79], SCSWC was stronger in the PRE. In the southern boundary of the current where the velocity was about 0.3 m × s −1 , a number of areas with high U. chinensis CPUE were formed (Figure 4a). In winter, under the control of the northeast monsoon [27], the southwesterly Guangdong Coast Current (GCC) started to form. At the southern boundary of the current where the velocity is 0.2 m × s −1 , several areas with high U. chinensis CPUE were identified. In the PRE, the current velocity in different seasons can be considered one of the indicators to predict the distribution of U. chinensis fishing grounds.
As one of most important indicators of the nutritional potential of the basic links in the marine food chain, primary productivity plays an important role in marine ecology [80]. Its size determines the potential yield of marine fishery [81] and can be used to indicate the spatiotemporal variation of fish stock [10]. In the spawning season, high NPP suggests high plankton biomass, an important feature for the growth and recruitment of juvenile fish. In the catching season, the NPP acts like the comprehensive index of biological concentration of squid bait, and the interaction between them can directly affect the annual stock level of squid [82]. A significant correlation between the interannual stock variation of Ommastraphes bartrami and the average NPP was observed in spawning and catching months [82]. Hence, high primary productivity areas caused by strong upwelling were generally considered as potential spawning grounds [83]. However, the relationship between primary productivity and squid was also affected by the cascade effect. The increase of primary productivity did not act on squid directly, and high levels of primary productivity can eventually lead to the increase of both prey and predator. On the western coast of California, the increase of primary productivity anomaly led to the decline of squid stock after three months [84]. In the PRE, the suitable NPP of U. chinensis differed depending on the month. In summer, the high resource was mainly distributed in the narrow zone with a NPP of about 3-5 mg × m −3 day −1 , while in winter it was mainly distributed in the water area with a NPP of 16-19 mg m −3 day −1 (Figure 4c,d).
In summer, the PRE was majorly affected by the SCSWC and the Pearl River plume (PPP) (Figure 4) and the environmental characteristics on both sides of the SCSWC were significantly different. During this period, the PRE was in the wet season and a great amount of fresh water from land discharged to the PRE [85]. Therefore, the distribution of SSS was mainly affected by the PPP. The U. chinensis was mainly distributed in the waters outside the SCSWC, located in the area where the diluted water from PRE and the high salinity water from SCS were mixed. In winter, the area was mainly influenced by GCC [85]. In this case, the species could be found in the area where high salinity and low salinity seawater mixed and its CPUE was lower than in summer ( Figure 4). Therefore, the relationship between the U. chinensis and the environment was similar in winter and summer, and it was adjusted by the changes of currents and marine environments. In the SCS, the high Sthenoteuthis oualaniensis CPUE was concentrated in the area with a high SST. With the seasonal warming, the same species shifted to lower latitudes [13], similar to the migration characteristics of the U. chinensis presented throughout this study. In addition, the U. chinensis CPUE in summer was higher than in winter, an observation that is consistent with the peak fishing season of squid in the Beibu Gulf [86]. However, in the Gulf of Cadiz, according to the bottom trawl survey, the abundance of Loligo vulgaris was the highest in autumn [52]. This difference may be related to the different fishing methods, latitude, and spawning periods, which need to be further explored in related future research.

Conclusions
This paper studied the effect of marine environmental factors on the spatiotemporal distribution of U. chinensis in the PRE on the basis of the long-term satellite remote sensing and survey data. It was observed that the interaction effect of SST and month was the most important environmental factor affecting the U. chinensis stock (accounting for 26.7% of the U. chinensis CPUE), followed by the interaction effect of depth and month (accounting for 12.8% of the U. chinensis CPUE). In the PRE, U. chinensis was mainly distributed in the sea area with an SST of 22-29 • C, SSS of 32.5-34‰, Chl a of 0-0.3 mg × m −3 , and water depth of 40-140 m. It is important to clarify that, in this study, only data related to SST, SSS, Chl a, depth, NPP, and current available from satellite remote sensing were analyzed. Morphological characteristics such as body length, age, and parameters such as dissolved oxygen and transparency were not considered and represent possibilities for follow-up studies to improve the accuracy of the model, thus providing a scientific basis for protecting the habitat of the U. chinensis.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy policy.

Conflicts of Interest:
The authors declare no conflict of interest.