Response to Environmental Factors of Spawning Ground in the Pearl River Estuary, China

: Spawning grounds are important areas for ﬁsh survival and reproduction, and play a key role in the supplement of ﬁshery resources. This study investigated environmental effects on the spatiotemporal variability of spawning ground in the Pearl River Estuary (PRE), China, using the generalized additive model (GAM), based on satellite remote sensing (sea surface temperature ( SST ), chlorophyll-a concentration ( Chl-a ), sea surface salinity ( SSS ), depth), and in situ observations. Results showed that 39.8% of the total variation in ﬁsh egg density was explained by these factors. Among them, the most important factor was SST , accounting for 14.3%, followed by Depth , SSS , and Chl-a , with contributions of 9.7%, 8.5%, and 7.3%, respectively. Spawning grounds in the PRE were mainly distributed in the waters with SST of 22 ◦ C, depth of 30–50 m, SSS of 16–35‰, and Chl-a of 6–15 mg/m 3 . From spring to summer, the spawning ground moved from the outlet of the PRE to the east. The distribution of the spawning ground in the PRE was mainly affected by the Pearl River Plume (PRP), Guangdong Coastal Current (GCC), and monsoons in this area.


Introduction
Spawning grounds are the areas where fish, shrimps, and shellfish mate, spawn, hatch, and breed, and are critical for the survival and reproduction of aquatic animals [1]. Fish eggs, as an important biological group, are the foundation for the replenishment of fishery resources and sustainable utilization of fishery resources [2]. Therefore, it is of great significance to understand the distribution of fish eggs for maintaining the balance of marine ecosystems [3]. As a critical and vulnerable period in the life cycle of fish populations, the distribution of fish eggs is influenced not only by the spawning behavior of fish (including spawning locations, timing, etc.) and ocean dynamics phenomena (currents, eddies, upwelling and stratification of the water column, etc.) [4][5][6][7][8], but also by water temperature, salinity, water depth, and topography [9][10][11]. Studies have shown that the timing of spawning population emergence in fish distributed in the English Channel was related to water temperature. In warmer years, the spawning season of the summer spawning population came earlier; while in colder years, the spawning season of the winter spawning population came earlier [11]. The distribution of Lutjanus campechanus in the East Seto Inland Sea was related to the substrate, water depth, and temperature, and the sandy substrate area was more suitable for spawning [12]. The distribution of fish eggs in the Mediterranean Sea was related to depth, wind direction, coastal eddies current direction, and topographic features of the coastline, and the nearshore ecosystem of the rocky coasts resulted in the accumulation of eggs [5,7,13]. The Baltic Sea had a strong thermocline and The data of fish eggs were obtained from spawning ground surveys from 2014-2017 (April-July) ( Table 1). The study area is located at 112-116 • E, 20-23 • N ( Figure 1). The plankton nets (shallow water type I) were used in the waters less than 30 m, and macroplankton nets were used in the waters more than 30 m. The trawling depth from the bottom to the surface and different water layers were not sampled separately. In the area of water depth less than 50 m, sampling every 3 h, samplings were performed a total of 9 times; in the area of water depth greater than 50 m, sampling every 4 h, samplings were performed a total of 7 times. During the trawling, the wire rope inclination angle was less than 45 • . Fish egg samples were preserved in 5% formaldehyde solution and taken back to the laboratory for identification and counting. Fish eggs were identified and classified based on their morphological characteristics [28,29]. The major species of eggs in this study were Carangidae, Trichiurus haumela, Thunnidae, Sardinella aurita, and Stolephorus commersoni. The egg data were grouped by 0.5 × 0.5 • grid cells, and the formula for calculating the density of individual fish eggs is as follows.
where C is the fish-egg density, in ind/(10m) 3 . N is the number of fish eggs, in ind. V is the volume of filtered water, in (10m) 3 . water depth less than 50 m, sampling every 3 h, samplings were performed a total of 9 times; in the area of water depth greater than 50 m, sampling every 4 h, samplings were performed a total of 7 times. During the trawling, the wire rope inclination angle was less than 45°. Fish egg samples were preserved in 5% formaldehyde solution and taken back to the laboratory for identification and counting. Fish eggs were identified and classified based on their morphological characteristics [28,29]. The major species of eggs in this study were Carangidae, Trichiurus haumela, Thunnidae, Sardinella aurita, and Stolephorus commersoni. The egg data were grouped by 0.5 × 0.5° grid cells, and the formula for calculating the density of individual fish eggs is as follows.
where C is the fish-egg density, in ind/(10m) 3 . N is the number of fish eggs, in ind. V is the volume of filtered water, in (10m) 3 .

Environmental Data
Satellite data include sea surface temperature (SST), sea surface salinity (SSS), sea surface chlorophyll a concentration (Chl-a), water depth (Depth), and currents (Current) ( Table 2). SST and Chl-a data are from MODIS Aqua products, acquired from April 2014 to July 2017 (http://oceancolor.gsfc.nasa.gov). SSS and Current data come from the Global

Environmental Data
Satellite data include sea surface temperature (SST), sea surface salinity (SSS), sea surface chlorophyll a concentration (Chl-a), water depth (Depth), and currents (Current) ( Table 2). SST and Chl-a data are from MODIS Aqua products, acquired from April 2014 to July 2017 (http://oceancolor.gsfc.nasa.gov (accessed on 31 May 2021)). SSS and Current data come from the Global Ocean Physical Reanalysis Product of Copernicus Marine Environment Management Service (CMEMS), acquired from April 2014 to July 2017 (http://marine.copernicus.eu/ (accessed on 31 May 2021)). Depth data come from Google Earth elevation data with an elevation class of 18. Water depth and stations are one-by-one correspondence. R v.4.0.0 (R Development Core Team, 2020) [30] software was applied to remove invalid values, month-by-month averaging, and data fusion for SST, SSS, and Chl-a. The same month, the same grid of fish egg data, and environmental data are matched together for subsequent statistical analysis. ArcGIS 10.3 (Esri, Redlands, CA, USA) was used to map the spatial distribution of environmental factors and fish eggs [31][32][33].

GAM
The generalized additive model (GAM) is a nonparametric extension of the generalized linear model (GLM), which can explain the nonlinear relationship between the dependent variable and multiple independent variables [21]. The basic GAM is as follows: where Y is the egg density, and log transformation is performed for Y + 1 to incorporate zero values in the response variable; s(x) is spline-smoothing function of covariate x; SST is the sea surface temperature; SSS is the sea surface salinity; Chl-a is the chlorophyll-a concentration; Depth is the water depth; and ε is the modelling error. In the rug plots of GAMs, the x-axis indicates the relative density of data points. As the y-axis scale is relative, the mean effect of each predictor on the response variable is indicated by a y-value of zero, with positive values indicating a positive effect and negative values indicating a negative effect. The "mgcv" package in software R v.4.0.0 was used to construct and test the GAM model [34].

Model Testing
Akaike information criterion (AIC) was used to test the fit of the model after adding the influence factors. The smaller is the AIC, the better is the fit [35]. Generalized cross validation (GCV) was applied to evaluate predictor variables. The smaller the GCV value is, the stronger the generalization is [36,37]. A stepwise method was applied to test the influencing effects of selecting factors including SST, SSS, Chl-a, and Depth. The influencing factors were added sequentially to the model to select the best model according to the reduction of GCV and AIC values [37]. The significance of these influencing factors was assessed by F-test and chi-square test [38][39][40].
The calculation formula for the AIC value is as follows: where θ is the deviation, df is the degree of freedom, and ϕ is the variance.

GAMs Analysis
In GAMs, the contribution of each selected factor to the egg density indicates the degree of influence of each factor on the egg density (Table 3). Among them, SST was the most influential factor affecting egg density, with a contribution of 14.3%, followed by Depth, SSS, and Chl-a, with contributions of 9.7%, 8.5%, and 7.3%, respectively. The chi-square test and F-test showed that all factors selected in the model were significantly correlated with egg density (p < 0.001; Table 3).
GAM analysis of the effects of environmental and spatial factors on fish egg density in the PRE is shown in Figure 2. In the ranges of 22-24 • C and 29-30 • C, egg density decreased with increasing SST with narrow confidence intervals and high confidence levels. In the range of 24-29 • C, egg density was stable with increasing SST with narrow confidence intervals and high confidence levels ( Figure 2a). The contribution of Depth to egg density in GAMs was 9.7%. When Depth was in the range of 0-50 m, the egg density increased with increasing Depth, with a narrow confidence interval and high confidence level. While in the range of 50-280 m, the egg density decreased with increasing Depth, the confidence interval increased, and the confidence level decreased (Figure 2b). When SSS was in the range of 10-17‰, the egg density increased with the increasing SSS, with a wide confidence interval and low confidence level. While in the range of 17-35‰, the egg density fluctuated with the increase of SSS, with a narrow confidence interval and high confidence level (Figure 2c). In the range of 1-10 mg/m 3 , the egg density increased with increasing Chl-a, the confidence interval increased with increasing Chl-a, and the confidence level decreased (Figure 2d).  GAM analysis of the effects of environmental and spatial factors on fish egg density in the PRE is shown in Figure 2. In the ranges of 22-24 °C and 29-30 °C, egg density decreased with increasing SST with narrow confidence intervals and high confidence levels. In the range of 24-29 °C, egg density was stable with increasing SST with narrow confidence intervals and high confidence levels ( Figure 2a). The contribution of Depth to egg density in GAMs was 9.7%. When Depth was in the range of 0-50 m, the egg density increased with increasing Depth, with a narrow confidence interval and high confidence level. While in the range of 50-280 m, the egg density decreased with increasing Depth, the confidence interval increased, and the confidence level decreased (Figure 2b). When SSS was in the range of 10-17‰, the egg density increased with the increasing SSS, with a wide confidence interval and low confidence level. While in the range of 17-35‰, the egg density fluctuated with the increase of SSS, with a narrow confidence interval and high confidence level (Figure 2c). In the range of 1-10 mg/m 3 , the egg density increased with increasing Chl-a, the confidence interval increased with increasing Chl-a, and the confidence level decreased (Figure 2d). The cumulative explained deviations of spatiotemporal and environmental factors on fish egg density were obtained using GAMs to fit and predict the effects of the added variables on the model ( Table 4). The influencing factors selected based on AIC and GCV values included SST, SSS, Chl-a, and Depth. The cumulative explanatory bias of these factors to fish egg density was 39.8%, with an adjusted R 2 of 0.35 (Table 4). According to the model checking plots for GAMs, the distribution of the residuals conforms to a normal distribution. Therefore, the prediction results are relatively satisfactory (Figure 3).  The cumulative explained deviations of spatiotemporal and environmental factors on fish egg density were obtained using GAMs to fit and predict the effects of the added variables on the model ( Table 4). The influencing factors selected based on AIC and GCV values included SST, SSS, Chl-a, and Depth. The cumulative explanatory bias of these factors to fish egg density was 39.8%, with an adjusted R 2 of 0.35 (Table 4). According to the model checking plots for GAMs, the distribution of the residuals conforms to a normal distribution. Therefore, the prediction results are relatively satisfactory (Figure 3).

Spatiotemporal Distribution of Egg Density and SST, SSS, Chl-a
Higher egg density (10,000-25,000 ind/(10m) 3   SSS in the PRE was 12-35‰ in spring and 5-35‰ in summer (Figure 4c,d). Waters with low SSS (12-20‰) spread from the PRE to the southwest in spring and to the northeast in summer. High egg density (10,000-25,000 ind/(10m) 3 ) was found in waters with SSS 12-20‰ in spring, and in waters with SSS 5-15‰ in summer.

Spatiotemporal Distribution of Fish Egg Density and Current
Egg density and flow field superimposed maps were obtained by kriging interpolation in the whole area, based on the egg density data in each station ( Figure 5). The Pearl River Plume (PRP) fresh water at the PRE and the coastal area forms an approximate plume current under the joint effects of buoyancy and external forcing [41]. In spring, the direction of the current in coastal waters west of the mouth of PRE is southwestward, while in the offshore waters, the direction of current is northeastward. As a result, a cyclonic vortex formed between the Pearl River runoff and the high salinity seawater off the coast, and the density of fish eggs was higher (4000-8000 ind/(10m) 3 ) (Figure 5a). In addition, in the waters near the PRE, multiple small-and medium-scale eddies formed due to the confluence of salt and freshwater, and the egg density was also higher in these areas (8000-18,000 ind/(10m) 3 ) (Figure 5a). In coastal waters east of the PRE, as the direction of the current is northeastward, the egg density was higher in waters at the mouth of Daya Bay and Dapeng Bay (4000-10,000 ind/(10m) 3 ) (Figure 5a). The egg density was lower in offshore waters (0-1000 ind/(10m) 3  Bay and Dapeng Bay (4000-10,000 ind/(10m) 3 ) (Figure 5a). The egg density was lower in offshore waters (0-1000 ind/(10m) 3 ) than in coastal waters (4000-18,000 ind/(10m) 3 ) (Figure 5a). In summer, the direction of the current in the PRE was northeastward, and the higher egg density (4000-18,000 ind/(10m) 3 ) was mainly distributed in the waters east of the PRE. Among them, the highest egg density (8000-18,000 ind/(10m) 3 ) was found in Dapeng Bay (Figure 5b). The egg density in the waters west of the PRE was about 0-4000 ind/(10m) 3 (Figure 5a).

Relationship between Egg Density and SST, SSS, Chl-a, Depth
GAM analysis of the environmental and spatial factors selected in this study on egg density showed that SST was the most important factor affecting egg density, with a contribution of 14.3% (Table 3). Fish eggs, which are in the early stage of fish life history, are sensitive to changes in their habitat. Studies have shown that water temperature plays an important role in fish distribution and spawning [42]. It is an important factor affecting fish maturation and egg development [9], controlling the timing of fish spawning and the number of adults entering the spawning grounds [4,9,43]. Higher water temperatures promote fish gonad development and increase spawning activity [44]. This study showed that eggs of major economic fishes in the PRE were mainly distributed in waters with SST 22-30 °C, and the most suitable SST was 22 °C, followed by 25-28 °C (Figure 2a). Major economic fishes in the PRE are Trichiurus haumela, Trachinidae, Nemipteras virgatus, Pneumatophorus japonicus, etc. [45][46][47]. Among them, the suitable spawning water temperature for Trichiurus haumela is 25-28 °C [45]. The suitable spawning water temperature for Pneu- In summer, the direction of the current in the PRE was northeastward, and the higher egg density (4000-18,000 ind/(10m) 3 ) was mainly distributed in the waters east of the PRE. Among them, the highest egg density (8000-18,000 ind/(10m) 3 ) was found in Dapeng Bay (Figure 5b). The egg density in the waters west of the PRE was about 0-4000 ind/(10m) 3 (Figure 5a).

Relationship between Egg Density and SST, SSS, Chl-a, Depth
GAM analysis of the environmental and spatial factors selected in this study on egg density showed that SST was the most important factor affecting egg density, with a contribution of 14.3% (Table 3). Fish eggs, which are in the early stage of fish life history, are sensitive to changes in their habitat. Studies have shown that water temperature plays an important role in fish distribution and spawning [42]. It is an important factor affecting fish maturation and egg development [9], controlling the timing of fish spawning and the number of adults entering the spawning grounds [4,9,43]. Higher water temperatures promote fish gonad development and increase spawning activity [44]. This study showed that eggs of major economic fishes in the PRE were mainly distributed in waters with SST 22-30 • C, and the most suitable SST was 22 • C, followed by 25-28 • C (Figure 2a). Major economic fishes in the PRE are Trichiurus haumela, Trachinidae, Nemipteras virgatus, Pneumatophorus japonicus, etc. [45][46][47]. Among them, the suitable spawning water temperature for Trichiurus haumela is 25-28 • C [45]. The suitable spawning water temperature for Pneumatophorus japonicus is 22-29 • C [46] and is 18.79-25.98°C for Nemipteras virgatus [47], which is basically consistent with the results of this study. In this suitable temperature range, the overall egg density is negatively correlated with water temperature (Figure 2a). For fish eggs with a higher suitable temperature, higher water temperature promotes shorter egg development time, and accelerates the speed of fish egg hatching, which also leads to a decrease in the fish egg density [13,48,49].
In GAMs, SSS contributed 8.5% to the variation in egg density (Table 3). Salinity has long been considered a dominant factor affecting the distribution, survival, and community composition of estuarine fishes [50,51]. Studies have shown that lower or higher salinities can interfere with egg development and the hatching processes, and reduce the survival rate of larvae and juveniles [48,52]. Drastic changes in salinity can cause large-scale changes in the spawning area [53]. In the PRE, the optimum salinity for egg density was 17-35‰, and the egg density decreased when SSS was lower than 17‰ (Figure 2c). The suitable salinity range for spawning grounds of major economic fishes (Pneumatophorus japonicus, Nemipteras virgatus, and Trichiurus haumela) in the PRE is 31-35‰ [45][46][47], which is basically consistent with the results of this study. In addition, salinity also affects the buoyancy of fish eggs, which determines their vertical distribution in the waters. In waters with low salinity, fish eggs tend to accumulate and cannot get enough oxygen, which hinders egg development. In waters with higher salinity, fish eggs can be suspended or float in the water, which facilitates the absorption of oxygen from the surrounding water, and thus, improves the survival rate of fish eggs. Moreover, fish fertilization activity also increases with high salinity, leading to increased egg density [10,54]. In the PRE, there were multiple peaks in egg density in the range of SSS 17-35‰ (Figure 2c). On the one hand, it may be caused by the noise in the model fitting process, which attempted to match fish egg and salinity data from a limited number of records. On the other hand, it may be related to the characteristics of salty and freshwater confluence in the estuary that influence the abundance of eggs. The fishery resources in the area are diverse and complex in composition, forming three ecological groups: semi-saline, inshore, and offshore species [19,55]. Suitable salinity provides spawning grounds for fish with different reproductive habits [24,56].
The contribution of Chl-a to egg density was 8.07% (Table 3). The spawning goal of fish is to provide a favorable foraging and feeding environment for their larvae [6,57]. Spawning grounds usually coincide with bait distribution, which contributes to reproductive success and the increase of population replenishment [9,58]. Chlorophyll concentration reflects the standing stock of phytoplankton in the waters, which is the main feeding source of zooplankton and some marine organisms [59]. This study showed that the egg density generally increased with increasing Chl-a (Figure 2d), which is consistent with preferred conditions for spawning grounds. In addition, the effect of higher chlorophyll a concentration on fish populations usually has a lag effect. In the Gulf of Cadiz, anchovy adults required high chlorophyll a concentrations 3-4 days before the spawning process occurred to increase the food supply of anchovy larvae, so the chlorophyll a concentration with a time lag of 3 days was the environmental variable that best characterizes their spawning grounds [60]. The hysteresis effect was also evident in different stages of other organisms; for example, in the South China Sea, the hysteresis effect of chlorophyll a concentration on protozoa was about 2 months [61]. This phenomenon mainly depends on the trophic level and life stages of different organisms [60][61][62]. However, the egg density increased when Chl-a was less than 1 mg/m 3 (Figure 2d). This may be due to the increased transparency of seawater at lower Chl-a, which has a positive impact on the spawning activities of fish [63].
In the PRE, the contribution of Depth to the egg density was 5.1%, and the most suitable depth in the area was 30-50 m (Figure 2b). Water depth affects the spawning activity of fish by influencing the level of dissolved oxygen in the water [64]. Studies have shown that spawning fish populations have a clear preference for different depths [12]. The area of 30-50 m water depth is the most concentrated spawning ground in the PRE, where the major spawning grounds of Trichiurus haumela, Decapterus maruadsi, Pneumatophorus japonicus, Nemipteras virgatus, Upeneus sulphureus, and Priacanthus tayenus are distributed [45][46][47]65], with the highest egg density. Fewer fish spawned in areas with water depths less than 30 m and greater than 100 m [65], and egg densities were lower (Figure 2b). In addition, the egg density in the deep-water area was significantly lower than that in the shallow-water area in the PRE (Figure 5a,b). Studies have shown that eggs of major economic fishes in the western Guangdong waters were usually distributed in shallow waters, forming main spawning grounds [18], which is consistent with the distribution characteristics of spawning grounds in this study area. Coastal habitats are generally considered to provide more suitable conditions for the survival of fish eggs, larvae, and juveniles than other offshore waters [5,66], which may be related to the ability of shallow marine habitats to provide higher water mass stability and more abundant bait biomass. Another possible explanation is that this may be due to its unique recruitment strategy. For example, in Biscay Bay, although in that region the food availability offshore is approximately half that on the shelf, the predatory risk is also lower. Anchovies spawn along the coast, larval advection offshore, juvenile, and then swim to the coast, which is a way of moving that makes them the least vulnerable to predation. This interpretation is based on the possibility of the population exploiting the off the shelf waters for recruitment through a loophole of lower predation [67].

Spatiotemporal Variation of Fish Eggs
Although the spawning site preferences may be the biggest factor affecting the density distribution of fish eggs, the distribution of fish eggs is also subject to wind fields and currents to some extent because fish eggs are not capable of swimming [5,68]. The current in the PRE is controlled by monsoons, tides, river discharges, coastal currents, and mesoscale eddies [69][70][71]. Of those, the monsoon plays the most important role in the surface shape and extension of the Pearl River Plume [41,69,71]. In spring, the sea level in the PRE is higher than that in the western Guangdong waters, but lower than that in the eastern Guangdong waters. The surface fresh-water discharge mainly diffuses westward [72], forming a sea area with low SSS and high Chl-a in the western PRE (Figure 4c,e). While the fresh-water plume expands outward from the estuary, highly saline water between the estuary is transported to the shore during high tides. In the sea area beyond the front formed by the diffusion of freshwater (roughly bounded by the salinity 33 contour), seawater is driven by the southwest monsoon and flows in a northeastward direction. As a result, the formation of cyclonic eddies is caused by the shear effects of estuarine discharged water, intrusion of highly saline water from the inter-estuarine shelf, and northeastward flow of the shelf [73] (Figure 5a). The formation of such eddies is thought to be a key spawning ground for pelagic fish, as these features increase the feeding opportunities of larvae and juveniles [6]. Analysis based on the currents and egg density distribution showed that higher egg density (about 4000-6000 ind/(10m) 3 ) was in the center of this vortex, and the highest egg density (about 14,000-18,000 ind/(10m) 3 ) was in the central waters of the PRE (Figure 5a). This is partly due to the fact that in this area, the flow field is more complex, forming multiple eddies, which is conducive to the bait organisms. On the other hand, there are many islands here and complex topography are advantageous to the aggregation of fish eggs [18]. In addition, this area is where freshwater and seawater meet, with nutrient-rich waters and abundant bait organisms (Figure 4c,e). These topographic, current and bait conditions are beneficial to the formation of spawning grounds [74][75][76]. In summer, the Pearl River runoff enters the northern South China Sea, forming a freshwater hypersaline tongue near the mouth of the Pearl River, and expanding out to the sea or coast under the action of tides, monsoons, and coastal currents [77]. The surface-flushing freshwater expands in a "symmetrical expansion-like pattern" to both eastern and western Guangdong waters [72]. Under the influence of the southwest monsoon, the runoff from four outlets (Humen, Jiaomen, Hongqili, and Hengmen) and the east mouth of Modaomen in Lingding Bay spread eastward [78], and the area with high egg density also transfers from the mouth of PRE to the eastern Guangdong waters, forming the most suitable spawning ground near the Dapeng Bay (Figure 5b). Studies have shown that spawning in the bay and downstream of the cape is effective in reducing egg loss and maintaining juveniles in the spawning area, and therefore, can be a suitable spawning ground [42,79]. Under the influence of the southwest monsoon, the direction of the coastal current is mainly eastward [78] (Figure 5b), and the fish eggs near the mouth of PRE follow the current to spread eastward as well as southeastward (Figure 5b).
In this study, survey data cover spring and summer, which are major spawning seasons in the PRE. Fishery in the area has characteristics of multi-species fishery resources, and single species have large varieties and small quantities [18,78]. Therefore, fish eggs are analyzed as a whole in this study. Environmental factors such as SSS, SST, and Chla available from satellite remote sensing are considered, and other influencing factors including fishing vessel parameters, egg biological characteristics, current fields, wind fields, eddies, water turbidity, and predation, which may improve the accuracy of the model will be considered in a future study.

Conclusions
For the first time, the environmental effects of the spatiotemporal distribution of fish eggs were quantified in the PRE. SST was the most important environmental factor affecting fish egg density (contribution of 14.3%), followed by depth, SSS, and Chl-a. The spawning grounds in the PRE are mainly distributed in the waters with SST of 22 • C, water depth of 30-50 m, SSS of 16-35‰, and Chl-a of 4-10 mg/m 3 . The egg density along the coast is higher than that of the outer sea, and the Pearl River plume is favorable for the aggregation of fish eggs. Major spawning grounds are found in the PRE in spring and in the eastern PRE (Dapeng Bay) in summer. In the next study, a long term survey and analysis on fish eggs will be conducted, and more oceanographic parameters will be considered to clarify the early replenishment mechanism for fishery resources in the PRE.
Author Contributions: L.Y. and D.W. designed the study. J.Y. and L.Y. collected the fishery data. D.W. and P.C. analyzed the data. L.Y. and R.H. helped with data collection and analysis. J.Y. and D.W. wrote the article. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: Data available on request due to restrictions. 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.