Exploring the Spatial and Temporal Distribution of Frigate Tuna ( Auxis thazard ) Habitat in the South China Sea in Spring and Summer during 2015–2019 Using Fishery and Remote Sensing Data

: Frigate tuna Auxis thazard ﬁshery is considered a potential marine resource in the open South China Sea (SCS). However, little is known about the spatial and temporal distribution of its habitat, and how this may respond to environmental changes. Using ﬁsh survey and remote sensing data from 2015–2019, we applied generalized additive models to identify relationships between environmental factors and the distribution of A. thazard in the SCS. To examine seasonal patterns in the habitat of A. thazard in the SCS, we generated a habitat suitability index model using environmental factors screened by generalized additive models. Results showed that A. thazard migrates from south to north in the SCS, and its suitable habitat is patchily distributed. Signiﬁcant environmental factors affecting the habitat of A. thazard distribution were different in different seasons; we demonstrated A. thazard to be sensitive to Chl-a in spring (optimum 0.155, optimal range ~0.1252–0.1840), and in summer to be sensitive to SST (optimum 30.405, optimal range ~29.789–31.021) and SSH (optimum 0.741, optimal range ~0.618–0.864). Suitable habitat in spring occurs mainly in northeastern areas, while in summer it occurs mainly around the southeastern Nansha Islands. Compared with that in spring, the area of suitable habitat increases in summer, as does the habitat suitability index of the entire sea area. These results improve our understanding of environmental factors that affect the distribution of A. thazard habitat in the SCS, and provide a scientiﬁc basis for the development and management of A. thazard resources.


Introduction
Fishing grounds in the South China Sea (SCS) are diverse, including waters on the continental shelf, slope, in Beibu Gulf, and around coral reefs northeast of Hainan Island.Because of its unique geographical location, the fauna of the SCS has a tropical marine biota and high species diversity [1].Most fishery resources of economic value in the SCS have been exploited, with resources in many traditional fishing grounds having declined [2].In response to overfishing and coastal ecosystem degradation, the Chinese government introduced a series of measures aimed to rebuild marine fisheries, including input and output controls, technical limits such as gear specifications and restrictions, spatial and temporal closures, and other ecological instruments [3].Auxis thazard (Thunnidae) is a warm-water tuna that is widely distributed in the tropical and subtropical seas, including the East China seas and South China seas [4]; it is harvested using a wide variety of gear, including small-scale artisanal fisheries, light falling-net, and purse seine.In the SCS, the density of A. thazard in spring and summer is higher than that in autumn, and exhibits a patchy spatial distribution.According to the latest research, there are about 850,000 tonnes of small tuna species in the SCS, of which the A. thazard accounts for about 30% [5].The development of a new, offshore tuna fishery in the SCS could relieve some fishing pressure from inshore waters.
Because habitats represent a collection of environmental conditions that are suitable for species' survival, the marine environment factor can be used to speculate habitat distribution [6].Understanding the temporal and spatial dynamics of fish habitat can greatly improve ecosystem-based fisheries management, and an understanding of the distribution of fish resources [7].The habitat suitability index (HSI) model is a framework model that describes the quality of fish and wildlife habitats [8].It provides one method for linking the abundance of wild animals with habitat conditions [9] and can be used in combination with geographic information system (GIS) technology to create maps significant for research on fish habitat [10].Because the HSI performs well and its predictions in ecological research are reliable [11,12], it has been widely used to predict the spatial and temporal distribution of suitable fish habitats, such as for the identification of environmental conditions that affect the distribution of Chilean jack mackerel (Trachurus murphyi), and to quantify suitable habitat area (using SST, SSH and Chl-a) for this species in the southeast Pacific [13], and to detect habitat of yellowfin tuna in the western and central Pacific Ocean based on remote sensing data [9].
Relationships between habitat and the marine environment are complex.With increased fishing intensity, and changes in the ecological environment and global climate, the resource amount and spatial distribution of many marine species have recently changed, profoundly affecting their habitats and ecosystems [14].To develop a sustainable fishery for A. thazard it is important to identify relationships between the marine environment and this species' catch, and to screen the main environmental factors that affect the distribution of resources required by this species.
Generalized additive models (GAM) provide a flexible method to explore nonlinear relationships between fishery resources and the marine environment.They can be used for different species and different sea areas, and when combined with a variety of environmental factors, can better reveal relationships between the environment and the spatial and temporal distribution of fishery resources [15][16][17][18][19][20].Using GAM, SST was better than the other environmental factors (SSHA and chlorophyll-a concentration) as an oceanographic predictor to explore the relationship of environmental parameters of Bigeye Tuna in the eastern Indian Ocean off Java [21].For Thunnus atlanticus in the Gulf of Mexico, the most influential oceanographic conditions were salinity, temperature and sea surface height [22].
While the morphology, biology, spatial and temporal distribution of fishing grounds, population structure, and stock estimates of A. thazard have been investigated [23][24][25][26].In this study, we use GAMs to analyze environmental factors associated with the spatial and temporal distribution of A. thazard in the SCS during spring and summer, and construct an integrated HSI model based on these environmental factors.Our objectives are to explore the relationship between A. thazard and its environment, and to identify the main temporal and spatial distribution characteristics of habitat suitable for it, to provide a scientific basis for rational development and management of regional A. thazard resources.In the SCS, the fishing periods of A. thazard are especially common from spring to autumn, spring is the breeding season for A. thazard, studying the spawning habitat of A. thazard can help to understand the distribution of this species' larval and juvenile and its relationship to the environment [26].Since 1999, China has implemented the summer fishing moratorium from May to August each year, which prohibits fishing to operate except for fishing gear.Therefore, both main fishing methods for A. thazard (light falling-net and purse seine) are banned.During the summer fishing moratorium, decreased fishing intensity reduces human impacts on habitat distribution.So, we select the spring-summer period for the analyses.

Fisheries and Environmental Data
Fisheries data were collected through an independent fisheries survey in the SCS conducted by the the South China Sea Fisheries Research Institute of the Chinese Academy of Fisheries Science from 2015-2019 (Table 1).Catch data, longitude/latitude, gea, speed and course were recorded for each station.Surveys were conducted in offshore SCS waters between 10-20 • N and 110-118 • E (Figure 1).The survey vessel, with weight 421 t and length 43.6 m, was equipped with metal halide fish lamps (1 kW).The net, ~281.60 m × 81.76 m, had maximum and minimum meshes of 35 mm and 20 mm, respectively, and weighed 2816 kg.When surveying, 300 lights were switched on at 19:00 for ~3 h every night.The catch data were logarithmically processed to satisfy a normal distribution when used.As an environmental database, daily sea surface temperature (SST), sea surface temperature gradient (SSTG), Chlorophyll-a (Chl-a), sea surface height (SSH), sea surface salinity (SSS) and mixed layer depth (MLD) were used in this study.The units of SST, SSTG, Chl-a, SSH, SSS, and MLD were • C, • C, mg m −3 , m, PSU, and m, respectively.The SST and Chl-a concentration data are from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board the satellite Aqua platform and were obtained from National Oceanic and Atmosphere Administration (NOAA) (https://coastwatch.pfeg.noaa.gov/data.html,accessed on 26 March 2022).The data have a daily temporal resolution and a spatial resolution of 0.04 • .The SSH, SSS, and MLD are provided by Copernicus Marine Environmental Monitoring Service (https://resources.marine.copernicus.eu/product-detail/GLOBAL_MULTIYEAR_PHY_001_030/INFORMATION, accessed on 27 March 2022).The data have a daily temporal resolution and a spatial resolution of 0.083 • .In the product, the MLD is defined as a depth where the density increase compared to density at 10 m depth corresponds to a temperature decrease of 0.2 • C in local surface conditions (θ10 m, S10 m, P0 = 0 db, surface pressure).Based on maximum temperature gradient characteristic at the ocean temperature front, the SSTG is calculated using the following equation [27],

Model Construction
2∆Y where Grad SST is the sea surface temperature gradient; D X and D Y are the gradients of adjacent points in the X and Y directions, respectively; ∆X and ∆Y are the pixel spacing in the X and Y directions, respectively; and i and j denote the positions of pixel points on the image, respectively.
Since the spatial resolutions of each dataset are different they were re-sampled to fit the fishery data resolution.Daily values of environmental data were extracted from each pixel corresponding to the location of fishing activities by ArcGIS software (ESRI, Redlands, CA, USA).

Model Construction 2.2.1. GAM
GAM is a nonparametric extension of a generalized linear model (GLM) [28] which is more suitable for dealing with nonlinear problems.GAM focuses on exploring data nonparametrically to capture relationships between the marine environment and resource abundance [29].To avoid excessive computation and over-fitting, the linear correlation between each explanatory variable is measured according to a Variance Inflation Factor (VIF) before GAM is performed.When VIF > 10, explanatory variables have a linear relationship with each other and should be removed from the model [30,31].We used the R mgcv package to fit the GAM [32,33].The general form the model can be written as: where g(µ) is the link function, which is the logarithm of Catch; β is the intercept term; f is a smooth function; i is the number of predictor variables, and ε is an error term.The GAM procedure uses smoothing splines s( ).

HSI
The HSI model objectively estimates the suitable environmental range for a species to survive [34].Assuming the sea area with the largest Catch is the sea area with the most abundant resources, it follows that the sea area with the least Catch is the sea area with the Fishes 2022, 7, 218 5 of 17 most impoverished resources; the former suitability index (SI) is set to 1 and the latter to 0 [35].We derive SI using the following formula: where SI a is the SI of the environmental factor a; Catch i is the catch at station i; Catch max is the maximum catch at all stations in the voyage; Catc min is the minimum catch at all stations in the voyage.We fit the relationship between environmental factors and SI by nonlinear regression and solved by R software.We calculated the integrated habitat index model using the Arithmetic Mean Method (AMM) and verified it with survey data [36].The equation of our HSI model is: where SI 1 , SI 2 , • • • SI n , are the predicted SI values on the fishing ground for each environmental factor; and n is the number of environmental factors in the HSI model.If the differences between the theoretical and actual value exceed 0.3, the forecast is considered to be accurate, and vice versa [37].The HSI spatial distribution of each season was plotted, and the habitat area was calculated.The areas with HSI ≥ 0.6, HSI ≥ 0.8, 0.2 < HSI < 0.6, and HSI ≤ 0.2 were regarded as suitable, optimal, normal, and poor habitats, respectively, for A. thazard [38].

Cross Validation
The spring and summer optimal model were cross validated by dividing or partitioning the data into one portion used to calibrate the model (training data), and one portion to validate predictions, called the testing data.Huberty described a rule of thumb for determining the optimum partitioning of training and testing data for models, with the proportion of testing data being: where ρ is the number of predictors.Other data are used as training data for model calculations.

Time Delay Analysis
To determine if time-lagged effects of significant environmental factors existed on A. thazard habitat, spatiotemporal relationships between monthly average HSI values and environmental factors were further evaluated using cross-correlation plots.

Catch Data Distribution
The overall distribution of the logarithmic catch data is presented in Figure 2. Lg (Catch) of A. thazard was high in summer, peaking at 6.13, with an upper quartile of 5.41, lower quartile of 4.48, and a minimum of 2.25 (Figure 2).Catch of A. thazard (lg) was low in spring, peaking at 5.26, with an upper quartile of 3.37, lower quartile of 2.32, and minimum of 1.48.Comparing the summer and spring kernel density distribution plots of A. thazard lg (Catch) the distribution of A. thazard lg (Catch) in summer peaks at ~4.78, and that in spring peaks at ~3.13.All maxima, upper and lower quartiles, and median values for A. thazard lg (Catch) in summer were higher than those in spring.
A normality test for A. thazard lg (Catch) reveals the lg (Catch) data points to basically form a straight line in the normal QQ-plot, indicating that lg (Catch) conforms to the normal distribution and satisfies the basic premise of GAM construction (Figure 3).
A normality test for A. thazard lg (Catch) reveals the lg (Catch) data points to basically form a straight line in the normal QQ-plot, indicating that lg (Catch) conforms to the normal distribution and satisfies the basic premise of GAM construction (Figure 3).

Environmental Factor Screening
Collinearity diagnostics tests of temporal and environmental variables reveal the VIF of all environmental factors to be <10, indicating no collinearity exists amongst them (Table 2).Therefore, all environmental factors can be used in the GAM, which is constructed using complete subset regression, selecting the optimal combination of explanatory variables according to Mallows' Cp and Adjusted R-square (Tables 3 and 4).Four-fold cross-validation was used to evaluate the model; when the RMSE > 1, the model was deemed reasonable.Significant environmental variables (p ≤ 0.05) were different between spring and summer.In spring, the significant variable was Chl-a, while in summer they were SST and SSH.A normality test for A. thazard lg (Catch) reveals the lg (Catch) data points to basically form a straight line in the normal QQ-plot, indicating that lg (Catch) conforms to the normal distribution and satisfies the basic premise of GAM construction (Figure 3).

Environmental Factor Screening
Collinearity diagnostics tests of temporal and environmental variables reveal the VIF of all environmental factors to be <10, indicating no collinearity exists amongst them (Table 2).Therefore, all environmental factors can be used in the GAM, which is constructed using complete subset regression, selecting the optimal combination of explanatory variables according to Mallows' Cp and Adjusted R-square (Tables 3 and 4) Four-fold cross-validation was used to evaluate the model; when the RMSE > 1, the model was deemed reasonable.Significant environmental variables (p ≤ 0.05) were different between spring and summer.In spring, the significant variable was Chl-a, while in summer they were SST and SSH.

Environmental Factor Screening
Collinearity diagnostics tests of temporal and environmental variables reveal the VIF of all environmental factors to be <10, indicating no collinearity exists amongst them (Table 2).Therefore, all environmental factors can be used in the GAM, which is constructed using complete subset regression, selecting the optimal combination of explanatory variables according to Mallows' Cp and Adjusted R-square (Tables 3 and 4).Four-fold cross-validation was used to evaluate the model; when the RMSE > 1, the model was deemed reasonable.Significant environmental variables (p ≤ 0.05) were different between spring and summer.In spring, the significant variable was Chl-a, while in summer they were SST and SSH.

Analysis of HSI Model
We fitted the relationship between lg (Catch) and environmental factors that differed significantly between seasons (Chl-a, SST, SSH) using unary nonlinear regression (Figure 4, Table 5), and calculated SI value for each significant environmental factor.The significance level of each function expression is p < 0.05.The HSI model is obtained using the arithmetic average method, and the model accuracy verified by the 0.3 standard (Table 6), if the difference between the predicted value and the actual value is less than 0.3, the forecast is considered correct, otherwise, the forecast is considered wrong.In spring, the numbers of correct and incorrect forecasts were 47 and 0, respectively, and the model accuracy was 100%.In summer, the numbers of correct and incorrect forecasts were 72 and 28, respectively, and the model accuracy was 72%.The average model accuracy of both seasons was 86%.Therefore, the fit is reliable and precise.The most suitable value of significant environmental factors in each season was obtained when SI was close to 1, and the greatest range of values occurred when SI ≥ 0.6.The optimum Chl-a value in spring was 0.154604, with an approximate optimum range of 0.1252-0.1840; the optimum SST value in summer was 30.405269, with an approximate optimum range of 29.7892-31.0214;and the optimum SSH value in summer was 0.741042, with an approximate optimum range of 0.6178-0.8641(Figure 4).

Habitat Variation
Mapping the A. thazard HSI distribution in spring and summer using ArcGIS 10.8 revealed suitable habitats to be mainly patchy (Figure 5).Suitable habitat in spring occurs mainly in northeastern sea areas, and in patches southwest of Paracel Islands near Vietnam and east of Nansha Islands, such as Renai Reef and Mischief Reef (Figure 5).Suitable habitats in summer occur mainly around Nansha Islands in southeastern waters, with patchy habi also found in southwestern waters near Vietnam.Additionally, compared with that in spring, the suitable habitat completely moved from the northeast areas to the southeast areas in summer, and the area of suitable habitat increased significantly in summer, and the SI of the entire sea area increased (Table 6).Fishes 2022, 7, x FOR PEER REVIEW 11

Relationship between Significant Environmental Factors and Habitat Quality
Cross-correlation analysis revealed a significantly positive relationship bet average HSI and Chl-a (Figure 6 1.

Discussion
March-May (spring) is the breeding date of A. thazard [25,[39][40][41].The effect of w temperature on fish is especially obvious during pre-oviposition and oviposition per Generally speaking, suitable temperature ranges for spawning are narrower than t for survival, with temperatures below or above these ranges limiting gonad develop or maturation [42].Rising water temperatures during spring would theoretically re the reproductive capacity of A. thazard in southern areas.Changed currents and nutrient loadings in more northern waters might explain the migration of A. thaza them because of their increased suitability for larval and juvenile growth [43].During  1, Summer Voyage a, b, c, which are voyages 4, 6, and 8 in Table 1.

Discussion
March-May (spring) is the breeding date of A. thazard [25,[39][40][41].The effect of water temperature on fish is especially obvious during pre-oviposition and oviposition periods.Generally speaking, suitable temperature ranges for spawning are narrower than those for survival, with temperatures below or above these ranges limiting gonad development or maturation [42].Rising water temperatures during spring would theoretically reduce the reproductive capacity of A. thazard in southern areas.Changed currents and high nutrient loadings in more northern waters might explain the migration of A. thazard to them because of their increased suitability for larval and juvenile growth [43].During this Fishes 2022, 7, 218 13 of 17 time the southwest monsoon drives surface waters rich in nutrients and plankton northeast to the continental shelf and slope, providing suitable conditions for A. thazard larval growth and development [44,45].
From August-September (summer), the size of the A. thazard population experiencing spawning migration and reproduction increases rapidly, and larva and juveniles drifted to each distribution area with the current and monsoon [45].Adult fish migrate to areas with abundant food and gain conditions [40].Suitable habitat for A. thazard occurs mostly near the Nansha Islands from August to September (summer) (Figure 5).In the shallow waters of Nansha Islands, where the sun is high at noon, daylight hours are long, and water temperatures are high.Water currents in the SCS are complex from August to September.There are many islands and reefs around the Nansha Islands, and the water flow is affected by topographical obstacles.The Nansha Islands are located in the confluence of different flow systems and water masses, with the flow around the Nansha Islands also affected by terrain, creating many eddies [46,47].The eddy tends to cause strong vertical movements of seawater, resulting in the bottom water surging upwards and the surface water sinking downwards.Movements of seawater can promote the mixing of seawater and bring high-nutrient subsurface water to the surface layer, which can boost plankton growth and increase the aggregation of prey organisms [48].Thus, the Nansha Islands area is the important habitat of A. thazard (Figure 7).
Fishes 2022, 7, x FOR PEER REVIEW 14 of 18 time the southwest monsoon drives surface waters rich in nutrients and plankton northeast to the continental shelf and slope, providing suitable conditions for A. thazard larval growth and development [44,45].
From August-September (summer), the size of the A. thazard population experiencing spawning migration and reproduction increases rapidly, and larva and juveniles drifted to each distribution area with the current and monsoon [45].Adult fish migrate to areas with abundant food and gain conditions [40].Suitable habitat for A. thazard occurs mostly near the Nansha Islands from August to September (summer) (Figure 5).In the shallow waters of Nansha Islands, where the sun is high at noon, daylight hours are long, and water temperatures are high.Water currents in the SCS are complex from August to September.There are many islands and reefs around the Nansha Islands, and the water flow is affected by topographical obstacles.The Nansha Islands are located in the confluence of different flow systems and water masses, with the flow around the Nansha Islands also affected by terrain, creating many eddies [46,47].The eddy tends to cause strong vertical movements of seawater, resulting in the bottom water surging upwards and the surface water sinking downwards.Movements of seawater can promote the mixing of seawater and bring high-nutrient subsurface water to the surface layer, which can boost plankton growth and increase the aggregation of prey organisms [48].Thus, the Nansha Islands area is the important habitat of A. thazard (Figure 7).An accurate and precise species abundance index (AI) is generally considered to be most important for fisheries forecasting and stock assessment, and represents the most important input for a fishery model.Catch AI, CPUE AI, and Effort AI are common indicators of resource abundance when building fishery models [49,50].We use Catch AI An accurate and precise species abundance index (AI) is generally considered to be most important for fisheries forecasting and stock assessment, and represents the most important input for a fishery model.Catch AI, CPUE AI, and Effort AI are common indicators of resource abundance when building fishery models [49,50].We use Catch AI to construct an HSI model.From the SI curve, we determine that catch data are closely related to environmental conditions.High catch data indicate that A. thazard is concentrated within a region, indicating that environmental conditions are conducive to its survival.Another important reason to use these catch data is that they were obtained during scientific investigation rather than from commercial fisheries.Not only do catch levels not affect fishing behavior, but fishery data are collected using standard gear and protocols over a long time period, at defined spatial and/or temporal scales.Thus, these data theoretically provide an unbiased estimate of fish population parameters [51].In a study on the neon flying squid Ommastrephes bartramii, fishery models based on three AIs were compared, and results showed that for both simulated and real datasets, Catch AI was more accurate, consistent, and had better spatial representational ability at all spatial scales when compared with CPUE AI and Effort AI [52].Accordingly, Catch AI-based fishery models are important tools for assessing the impacts of environmental change on marine species habitat, and provide a scientific basis for resource management.
Obvious seasonal and spatial variations exist in environmental conditions and habitat patterns for A. thazard.Favorable environmental conditions play an important role in regulating the range and distribution of suitable habitats.In fact, for habitats, environmental factors affecting habitat distribution of different species, regions, and seasons all differ [53].When modeling fisheries, the fewer environmental variables used the better.Therefore, choosing fewer, more appropriate environmental variables can avoid overfitting and improve model interpretation.Common methods to statistically screen variables include stepwise regression (forward, backward, and forward-backward hybrid methods) and complete subset regression.In stepwise regression, each iteration can only follow the direction of the previous iteration [54].However, complete subset regression can fit the possible combination models of all variables, and select the optimal model with variable combinations according to R-squared, Mallows' Cp, Adjusted R-squared, AIC (Akaike Information Criterion), or other standards [55].For accurate prediction of habitat distribution, the selection of appropriate environmental factors is important.For example, through stepwise regression, the most important environmental factors affecting the distribution of bigeye tuna (Thunnus obesus) in southern Java-Bali waters were SST, SSC, and SSHD [56], but off western Guangdong, water depth had the greatest impact on the distribution of Japanese scad Decapterus maruadsi [38].Although variable screening through stepwise regression can establish a good model, it cannot guarantee that the constructed model is the best, because not every possible model is evaluated.To overcome this problem, we use complete subset regression to screen variables.Compared with stepwise regression, complete subset regression is computationally more complex, but it ensures that the selected combination of variables is optimal.According to Mallows' Cp and Adjusted R-squared, we selected the best model, environmental factors were statistically (objectively) screened, and we obtained a simple and effective model.
Statistically optimal combination of catch variables and significant environmental variables can be obtained through full subset regression methods.However, this does not mean that variables that are not screened out have no effect on habitat.For example, some studies showed that MLD had an important effect on tuna habitat distribution [57].In our study, the environmental data were all derived from satellite remote sensing, which is the secondary product of remote sensing data such as remote sensing images and radar data, and has certain limitations.Water depth, water turbidity, clouds, and reanalysis models will all affect the quality of secondary product, which in turn will affect the results of variable screening and model construction [58].Compared with traditional walk-around sampling and monitoring, remote sensing satellite have many irreplaceable advantages in constructing habitat models, such as quick access of information, consecutive data, and ability to cover large areas.Our study area is mainly located in the sea beyond the 200 m depth, where the water is relatively high-transparent and has relatively little influence on the inversion results of satellite remote sensing, enabling it to construct a credible habitat model.

Figure 1 .
Figure 1.Spatial distribution of survey stations for A. thazard in the South China Sea, spring to summer, 2015-2019.Inset with red rectangular box marks the study area.

Figure 2 .
Figure 2. Violin plot of lg (Catch) distribution of A. thazard for spring and summer.

Figure 3 .
Figure 3. QQ-plot of A. thazard lg (Catch) in spring and summer.

Figure 2 .
Figure 2. Violin plot of lg (Catch) distribution of A. thazard for spring and summer.

Figure 2 .
Figure 2. Violin plot of lg (Catch) distribution of A. thazard for spring and summer.

Figure 3 .
Figure 3. QQ-plot of A. thazard lg (Catch) in spring and summer.

Figure 3 .
Figure 3. QQ-plot of A. thazard lg (Catch) in spring and summer.

Fishes 2022, 7 , 18 Figure 5 .
Figure 5. Distribution of A. thazard HSI in spring and summer, South China Sea: (a) Distribution of A. thazard HSI in spring; (b) Distribution of A. thazard HSI in summer.

Figure 5 .
Figure 5. Distribution of A. thazard HSI in spring and summer, South China Sea: (a) Distribution of A. thazard HSI in spring; (b) Distribution of A. thazard HSI in summer.

Figure 6 .
Figure 6.Cross-correlation coefficients between average HSI and environmental factors for Auxis thazard in the South China Sea.Spring Voyage a, b, c, which are voyages 2, 3, and 5 inTable 1, Summer Voyage a, b, c, which are voyages 4, 6, and 8 in Table 1.

Figure 7 .
Figure 7. Environment impacts the habitat of A. thazard in summer.

Figure 7 .
Figure 7. Environment impacts the habitat of A. thazard in summer.
Figure 1.Spatial distribution of survey stations for A. thazard in the South China Sea, spring to summer, 2015-2019.Inset with red rectangular box marks the study area.

Table 2 .
Collinearity diagnostics tests of environmental variables.

Table 3 .
Generalized additive models (GAMs) fitted to A. thazard Catch and analysis of deviance.

Table 4 .
Significance test results of explanatory variables.

Table 5 .
SI model for each environmental factor that differed significantly between spring and summer.

Table 6 .
Areas of different habitat types in spring and summer.