Evaluating the Effects Related to Restocking and Stock Replenishment of Penaeus penicillatus in the Xiamen Bay, China

: The quantitative evaluation of restocking and stock replenishment is essential for providing operational feedback and implementing adaptive management for future restoration projects. Since 2010, approximately 700 million juvenile shrimp ( Penaeus penicillatus ) have been released into Xiamen Bay, Fujian Province, China, each year, through stock replenishment programs. The recruited shrimp were sampled through three-year bottom trawl surveys from 2014 to 2017. The biological characteristics and catch equation were used to evaluate the effect of restocking and stock replenishment. The analysis uses the FAO-ICLARM Stock Assessment Tool (FISAT II) program. In general, there are two sources of recruitments—one from spawning brood stock and the other from released juvenile shrimp. We constructed an evaluation model for an effect evaluation based on Baranov’s catch equation to separate the initial recruitment volume using survey data. The relationship between body weight and total length was W = 1.638 TL 2.9307 . There is no statistically signiﬁcant difference between males and females. The von Bertalanffy growth parameters derived for prawns, using FiSAT II, were L ∞ = 209.6 mm and K = 0.51 per year. In spring 2014, the initial resource amount was 49,200, while the ratio of effective recruitment and parent amount was 3.92. The survival rate of the released shrimp larvae, 1.88 ‱ , seems to be very unsatisfactory. The resource amount in summer and autumn is higher than in winter and spring. Obviously, the restocking effect is lower and the programs need to be improved. To improve the restocking effect, the replenishment performance should be adjusted to reduce the mortality rate and increase its release effectiveness. Therefore, corresponding implementations are recommended, including standard extensive culture, reduction in stress during transportation, and temporary culture.


Introduction
In many offshore and pelagic fisheries in the world, nearly 90% of marine fish resources have been fully exploited, overexploited, or depleted, and are affected by the degradation of aquatic habitats [1]. Global capture fisheries production has stagnated, while seafood demand has steadily increased [1,2]. However, for some coastal fisheries, restocking, stock replenishment, and sea ranching via the application of aquaculture technology are also expected to help to restore the lost production and possibly increase harvests beyond the original level [3]. The improvement of fisheries based on aquaculture is a set of management methods that involve the release of cultured organisms to enhance, protect, or restore the and restocking. We have particularly constructed an assessment model to evaluate the effect based on Baranov's catch equation, which can separate the amount of initial recruitment using survey data from 2014 to 2017. The aim of this research, which analyzes the biological characteristics and restocking effects of P. penicillatus in Xiamen Bay, is to provide an idea regarding the effect of fisheries replenishment and restocking, provide suggestions on how to improve the effects, and serve as a reference for formulating redtail prawn management policies and the conservation of fisheries resources.

Study Area
The study area is in Xiamen Bay (117 • 50 -118 • 20 E and 24 • 14 -24 • 42 N) and has a total marine area of 1281 km 2 (Figure 1). It is a semi-enclosed bay, and the whole area includes the eastern areas (including Tong'an Bay, the Dadeng-Xiaodeng region, and Weitou Bay) with a sandy substrate and the western areas (including Western Harbor and Jiulong River Estuary), which are characterized by estuarine processes [25]. The water depth around the bay ranges from 6 to 25 m, and the bay has a deep-water coastal area of 30 km [26]. The bay contains many habitat types (including mangrove wetlands and sandy beaches) and diverse biological groups [27]. Xiamen Bay is an important economic region for Xiamen as well as the whole of Fujian Province. There are numerous ports, transport infrastructure, shipbuilding, and petrochemical industries. Therefore, the area is also affected by intensive development activities, including shipping, aquaculture, reclamation, and tourism. These activities have led to several issues, including a decline in marine biodiversity, habitat loss, and water pollution [28].
To manage shrimp fisheries resources and sustainability in Xiamen waters, biological information and population assessments are needed. This research aims to analyze the biological aspect of redtail prawn (P. penicillatus) in Xiamen waters, consisting of the sex ratio, length-weight analysis, growth parameter, level of exploitation (natural mortality rate, fishing mortality rate, and total mortality rate), and evaluation of stock replenishment and restocking. We have particularly constructed an assessment model to evaluate the effect based on Baranov's catch equation, which can separate the amount of initial recruitment using survey data from 2014 to 2017. The aim of this research, which analyzes the biological characteristics and restocking effects of P. penicillatus in Xiamen Bay, is to provide an idea regarding the effect of fisheries replenishment and restocking, provide suggestions on how to improve the effects, and serve as a reference for formulating redtail prawn management policies and the conservation of fisheries resources.

Study Area
The study area is in Xiamen Bay (117°50′-118°20′ E and 24°14′-24°42′ N) and has a total marine area of 1281 km 2 (Figure 1). It is a semi-enclosed bay, and the whole area includes the eastern areas (including Tong'an Bay, the Dadeng-Xiaodeng region, and Weitou Bay) with a sandy substrate and the western areas (including Western Harbor and Jiulong River Estuary), which are characterized by estuarine processes [25]. The water depth around the bay ranges from 6 to 25 m, and the bay has a deep-water coastal area of 30 km [26]. The bay contains many habitat types (including mangrove wetlands and sandy beaches) and diverse biological groups [27]. Xiamen Bay is an important economic region for Xiamen as well as the whole of Fujian Province. There are numerous ports, transport infrastructure, shipbuilding, and petrochemical industries. Therefore, the area is also affected by intensive development activities, including shipping, aquaculture, reclamation, and tourism. These activities have led to several issues, including a decline in marine biodiversity, habitat loss, and water pollution [28]. The biotic survey was conducted at six sample stations (Figure 1). These stations are briefly described as follows: XM01 includes the estuary of Jiulong River with larger freshwater exchange and mud sedimentation, where the seabed substrata are muddy because The biotic survey was conducted at six sample stations (Figure 1). These stations are briefly described as follows: XM01 includes the estuary of Jiulong River with larger freshwater exchange and mud sedimentation, where the seabed substrata are muddy because of the river's material deposition. XM02 and XM03 are located within the shallow waters, which face the open sea, with a substrate with mud-sand mixing and fluent water exchange. Compared to site XM02, site XM03 is better sheltered from monsoon wind waves. In addition, XM02 is at the front of the open zone with the greatest depth. Additionally, it stands where the wave current is the strongest, but also has the most drifted sand. Sites XM04, XM05, and XM06 are behind Kinmen Island and Xiaokinmen Island, especially site XM06, which is located at the mouth of Tongan Bay, a semi-enclosed water environment. A lot of sludge is deposited at Tongan Bay behind the bay mouth because of the sea water backflow.

Biological Sampling
From May 2014 to February 2017, field surveys were conducted to explore the effect of stock replenishment and restocking of the P. penicillatus population in Xiamen Bay, China. A total of 6 sites ( Figure 1) were set up to conduct biological sampling. The "Minlongyu 62,678 Fishing Boat" in Longhai City was used, and twelve cruises were performed as a seasonal sampling in February, May, August, and November. The vessel was a single-boat truss bottom trawler with a main engine power of 330 kilowatts. The height of the bottom trawl net mouth was 2.5 m, the length of the net was 24 m, the mesh of the bag net was 20 mm, the width of the truss was 27 m, and the average towing speed was 2-3 n mile/h. Surveys were conducted in accordance with the "Regulations for the Survey of Marine Fishery Resources (SC/T 9403-2012)". Towing took place approximately parallel to the coastline at a speed of 2.5 knots for 1 h.
After lifting the net, all the catch was poured onto the deck. All unwanted debris, plants, and garbage were first removed from catches; thereafter, the remainder was sorted into fish and shrimp categories. The shrimp and fish were put in marked plastic bags, and samples were frozen for species identification and further analysis in a laboratory. The crustacean samples were thawed in the laboratory, rinsed, and species identification, counting, weighing (g), and biological parameter determination at substations were performed; males and females were identified; the total length (mm) and carapace length (mm) of all samples of P. penicillatus were measured. A total of 312 shrimps (172 females and 140 males) were collected for the present study.

Biological Characteristics
The length-weight relationships were determined according to the allometric equation given by Sparre and Venema [29], as follows: where W = weight (g), L = total length (mm), a is the proportionality constant, and b is the isometric exponent or slope indicating isometric growth, also known as the power exponent coefficients. The weight-length relationships were also determined for each sex. The t test was used to analyze differences in total lengths and the weights of the sexes. Growth in length and weight was analyzed separately for each sex using the von Bertalanffy growth function (VBGF). By means of the input data from the length frequencies and the FiSAT II program, the asymptotic length (L ∞ ) and the growth coefficient (K) were estimated.
Data analysis was conducted with the most recent version of FiSAT statistical software [30,31]. The instantaneous rate of natural mortality (M) was obtained using Pauly's empirical formula [32]. The annual average surface water temperature T was 21.6 • C in Xiamen Bay.
The instantaneous rate of total mortality (Z) was found from the estimate of the growth parameters (K, L ∞ ), using the length-converted catch curve method (FiSAT II software) [31].
The instantaneous rates of fishing mortality (F) were calculated by the subtraction of the estimates of M from Z, as follows: The exploitation rate was calculated as follows: Then, the total mortality (A) was obtained using the follow equation: The derived estimates were then compared with other studies to allow further assessment.

Biomass Estimation
A method for assessment of stocking was developed and applied based on stock assessment theory and methods. Shrimp biomass was calculated using the swept area method [33]. The swept area (a, nm 2 ) or 'effective path swept' for each tow was calculated as follows: where a is the sweeping area of each tow, v is the towing speed, t is the duration of the towing, and B is the width of the path swept by the trawl (the width of the truss bar in this study is 27 m). Catch rates were calculated as catch (C, kg) divided by the time spent trawling (t, h) and converted to catch per unit area (CPUA, kg/nm 2 = biomass b per unit area) by dividing by the swept area ((C/t)/(a/t) = C/a).
The formula of the average abundance (N, kg), total biomass, was calculated from the following: where C/a is the CPUA of all tows (kg/square kilometers), ai is the sweeping area of the i-th station, ci is the number of catches at the i-th station, q is the proportion of catches (0.5 in this study; we assumed that all shrimps in the path of the tow would be captured), Ar is the overall area (154.18 square kilometers in this study), and n is the total number of tows, which is derived from the number of survey stations in all seasons. The survey frequency is quarterly in this study. Therefore, the average abundance in that year is the average value from the four quarters.

Model Building and Biomass Assessment
Generally, the spawning season of P. penicillatus is between April and May in Xiamen Bay. To facilitate calculation and setting conditions, we assume that all spawning will be completed on 1 May (ts: time of spawning). The release time is generally 15 June. This time is usually around the international ocean week (8 June). Therefore, we assume that P is the number of spawning shrimps in the natural sea area, te is the annual release time of juvenile shrimp (te: release time of stock replenishment), tr and NR are the recruitment time and the number of recruitments for juvenile shrimp to enter the fisheries, respectively. At this time, there are two sources of recruitments NR; one is from spawning brood stock NRP and the other is from released juvenile shrimp NRr. The schematic diagram of assessment model is shown in Figure 2.
respectively. At this time, there are two sources of recruitments NR; one is from spawning brood stock NRP and the other is from released juvenile shrimp NRr. The schematic diagram of assessment model is shown in Figure 2.

NR = NRP + NRr
According to Baranov's catch equation, the average abundance during the interval is as follows: where is average abundance and N × A = total deaths. Then, a simplified algorithm can be conducted for three years, which can estimate the number of initial recruitments.
The average abundance for the first year can be constructed in the form of the following equations: The average abundance for the second and the third year can be constructed in the form of the following equations: According to Baranov's catch equation, the average abundance during the interval is as follows: where N is average abundance and N × A = total deaths. Then, a simplified algorithm can be conducted for three years, which can estimate the number of initial recruitments.
The average abundance for the first year can be constructed in the form of the following equations: The average abundance for the second and the third year can be constructed in the form of the following equations: where, the following apply: P is the amount of female brood stock in natural seas. This value represents the amount of brood stock in the previous period, which will determine the number of next-generation juvenile shrimp. C r is the conversion ratio, which represents the ratio between the recruitment amount entering the fishing area in natural seas and the amount of brood stock P. t is the period of released shrimps within the year. R t is the number of shrimps released from artificial farming each year.
α t is the ratio of the number of females to the total resources at time t. S is the survival rate of released shrimps entering the fishing area. A t is the total mortality of shrimp in a natural sea area at time t. Z t is the total mortality coefficient of shrimp in a natural sea area at time t. N is the average abundance of shrimp in a natural sea area at time t.
The above continuous Equation (4) includes a total of eight parameters, where R t can be obtained from the website of Xiamen Municipal Bureau of Marine Fisheries, A t and Z t can be obtained by the body length transformation catch curve method in the FiSAT II software, and t can be calculated by the sea sweeping area method and substituted into Equation (12). Finally, we can calculate the value of N 0 , C r , and S by solving the continuous equations.
We assumed that (1) the population is a one-year-old species, that is, one year of sexual maturity, and the parent dies after laying eggs; (2) there is no difference between the released population and the natural one in the sea area, and they are fully mixed; (3) after the release activity, the fishing and natural mortality coefficients of this population generation are constant; (4) the reproduction rate is consistent every year; (5) there is poor mobility of species, and it is assumed that there is no population moving in and moving out in the sea area.

Biological Information and Growth Estimation
During the period from 2014 to 2017, a total of 312 shrimps (172 females and 140 males), with a total of 7573.18 gm, were collected from 12 cruises. Among different stations, the largest number was a total of 77 prawns, 1423.51 gm, accounting for 24.68% of the total catch number and 18.80% of the total catch weight, obtained at station XM05. In contrast, the smallest number was 20 tails, 436.6 gm, accounting for 6.41% of the total catch number and 5.77% of the total catch weight, obtained at station XM02.
The number of shrimps caught in summer and autumn is higher than that in winter and spring. Among them, 107 shrimps were caught in summer and autumn, accounting for 34.29% of the total catch, while 21 shrimps were caught in winter and 77 in spring, which accounted for 6.73% and 24.68% of the total catch, respectively. The highest catch weight in autumn and the lowest in spring, which are 3623.30 g and 754.13 g, accounted for 47.84% and 9.96% of the total catch weight, respectively.
The size frequency data are plotted in Figure 3. The total length of the shrimps was 62 to 191 mm, the main range was 100 to 150 mm (63.14%), the weight range was 2.8 to 79.6 gm, and the main range was 10.0 to 20.0 gm (31.41%). There are obvious seasonal changes here. The figure also shows the season when the recruitment group enters the sea and the state of seasonal growth.
For females, the smallest individuals (75 mm TL) were observed in the months of May, June, August, and September, while the largest sizes (>215 mm) were found in all the months (Figure 2). Most individuals, however, fell within the range of 145 to 195 mm TL. The smallest male individual (65 mm TL) was found in the month of June, while large sizes (>215 mm) were found in the month of May (Figure 3). Most female shrimps fell within the range of 115 to 155 mm TL. These results reveal that males are generally larger than females.
The weight-length equations obtained were as follows ( Figure 4):  For females, the smallest individuals (75 mm TL) were observed in the months of May, June, August, and September, while the largest sizes (>215 mm) were found in all the months (Figure 2). Most individuals, however, fell within the range of 145 to 195 mm TL. The smallest male individual (65 mm TL) was found in the month of June, while large sizes (>215 mm) were found in the month of May (Figure 3). Most female shrimps fell within the range of 115 to 155 mm TL. These results reveal that males are generally larger than females.
The weight-length equations obtained were as follows ( Figure 4):  In general, there are two kinds of allometric growth patterns, positive allometric (b > 3) and negative allometric (b < 3). Positive allometric indicates that the growth in weight is dominant compared to length, while negative allometric indicates that the growth in length is more dominant than the growth in weight. This seems to imply that growth in In general, there are two kinds of allometric growth patterns, positive allometric (b > 3) and negative allometric (b < 3). Positive allometric indicates that the growth in weight is dominant compared to length, while negative allometric indicates that the growth in length is more dominant than the growth in weight. This seems to imply that growth in length is more dominant than growth in weight.
The estimation of the growth parameters of P. penicillatus uses the von Bertalanffy equation, with length-frequency data as the input data in the FISAT II program ( Figure 5). The growth coefficient (K) value is 0.51 year −1 , and the asymptotic total length (L ∞ ) can reach 209 mm. In general, there are two kinds of allometric growth patterns, positive allometric (b > 3) and negative allometric (b < 3). Positive allometric indicates that the growth in weight is dominant compared to length, while negative allometric indicates that the growth in length is more dominant than the growth in weight. This seems to imply that growth in length is more dominant than growth in weight.
The estimation of the growth parameters of P. penicillatus uses the von Bertalanffy equation, with length-frequency data as the input data in the FISAT II program ( Figure  5). The growth coefficient (K) value is 0.51 year −1 , and the asymptotic total length (L∞) can reach 209 mm. Regarding the weight-length equations obtained in this paper, W = 1.638 TL 2.9307 . Compared with the parameter of weight-length from the study of Chen and Hong [34], the coefficients of a and b are very close (Table 1). The results of the statistical F test show that there is no difference between males and females. However, Liu and Zhong [35] showed, in 1986, that males are heavier than females with the same body length. In contrast, the results of Wang et al. showed no difference between males and females [36]. Regarding the weight-length equations obtained in this paper, W = 1.638 TL 2.9307 . Compared with the parameter of weight-length from the study of Chen and Hong [34], the coefficients of a and b are very close (Table 1). The results of the statistical F test show that there is no difference between males and females. However, Liu and Zhong [35] showed, in 1986, that males are heavier than females with the same body length. In contrast, the results of Wang et al. showed no difference between males and females [36]. Comparing the results of Liu and Zhong in 1986 with this 2014-2017 study, it is shown that both males and females are heavier in the 1986 study than they are in the results of this study [35]. This seems to imply that the weight of P. penicillatus is lighter in this study. The reason for the difference may come from different sea areas; Liu and Zhong conducted their study in the South China Sea [35], and this study (2021) was conducted in the waters of Xiamen, using the same ocean fishing method. However, comparing the different environment, aquaculture, and natural sea, the body size of the aquaculture [36] is smaller than the natural sea [35] and this study. Regarding the difference in growth between the marine environment and aquaculture environment, the growth of shrimp may be caused by environmental differences.

Estimation of Mortality
The results of the analysis showed that the instantaneous rate of the natural mortality (M) of P. penicillatus was 0.5919 year −1 . In this study, the temperature used was 21.6 • C, which was the average water temperature in the three years in Xiamen waters. Then, the instantaneous rate of total mortality (Z) was obtained using the length-converted catch curve method in the FiSAT II software, Z = 1. 15 year −1 (Figure 6). Hereafter, the total mortality (A) was 0.683, the instantaneous rates of fishing mortality (F) = 0.5581 year −1 , and the exploitation rate (E) = 0.485 year −1 , respectively. According to Sparre and Venema (1992), the optimum exploitation ratio Eopt is 0.5, implying that the stock of P. penicillatus in Xiamen waters is heavily exploited. The optimum level of exploitation (Eopt) is 0.5, beyond which the stock was said to be over-exploited [29]. This level implies that the stock of the P. penicillatus in Xiamen waters is exploited in a suitable state. Komi et al. reported a fishing mortality of F = 0.39 year −1 , natural mortality of M = 1.640 year −1 , total mortality of Z = 2.030 year −1 , and exploitation rate of E = 0.1906 of Metapenaeus elegans in the Andoni River, Nigeria, which indicated that the stock was under-exploited [37]. Another species of shrimp, Metapenaeus elegans of Segara Anakan Lagoon Cilacap, Central Java, showed overexploitation of this population, obtained from F = 6.760, M = 1.430, and E = 0.8830 [38]. Nwosu also reported an over-exploitation case of P. notialis, with E = 0.77 in 2007 and E = 0.69 in 2008 [39]. Most species of penaeid shrimps are over fished [40][41][42]. The exploitation rate could help understand the losses caused by fishing efforts and natural deaths, so that recommendations for management can be made. The catch deaths and natural deaths were almost equal, which also means that the biomass was not wasted in natural death. Therefore, the exploitation rate of P. penicillatus in Xiamen waters was conducive to management purposes. Such conditions could be harvested sustainably to maintain the socio-economic capacity and development of coastal communities. Assuming that the instantaneous rates of total mortality (Z) and total mortality (A) in three years are the same, the data obtained above are substituted into 1.15/year and 0.683 respectively. According to the survey data, the proportion of females in these three years are 0.5, 0.6512, and 0.5046, respectively. At the same time, using the method of sweeping the sea area, it can be calculated that the average volume in the past three years is 76,000, 103,100, and 64,600. According to Equation (4), it can be calculated that for the initial female resource of P. penicillatus in the waters of Xiamen, the estimate was 49.2 thousand in the spring of 2014, while the ratio between the effective recruitment amount The optimum level of exploitation (Eopt) is 0.5, beyond which the stock was said to be over-exploited [29]. This level implies that the stock of the P. penicillatus in Xiamen waters is exploited in a suitable state. Komi et al. reported a fishing mortality of F = 0.39 year −1 , natural mortality of M = 1.640 year −1 , total mortality of Z = 2.030 year −1 , and exploitation rate of E = 0.1906 of Metapenaeus elegans in the Andoni River, Nigeria, which indicated that the stock was under-exploited [37]. Another species of shrimp, Metapenaeus elegans of Segara Anakan Lagoon Cilacap, Central Java, showed overexploitation of this population, obtained from F = 6.760, M = 1.430, and E = 0.8830 [38]. Nwosu also reported an overexploitation case of P. notialis, with E = 0.77 in 2007 and E = 0.69 in 2008 [39]. Most species of penaeid shrimps are over fished [40][41][42]. The exploitation rate could help understand the losses caused by fishing efforts and natural deaths, so that recommendations for management can be made. The catch deaths and natural deaths were almost equal, which also means that the biomass was not wasted in natural death. Therefore, the exploitation rate of P. penicillatus in Xiamen waters was conducive to management purposes. Such conditions could be harvested sustainably to maintain the socio-economic capacity and development of coastal communities. Assuming that the instantaneous rates of total mortality (Z) and total mortality (A) in three years are the same, the data obtained above are substituted into 1.15/year and 0.683 respectively. According to the survey data, the proportion of females in these three years are 0.5, 0.6512, and 0.5046, respectively. At the same time, using the method of sweeping the sea area, it can be calculated that the average volume in the past three years is 76,000, 103,100, and 64,600. According to Equation (4), it can be calculated that for the initial female resource of P. penicillatus in the waters of Xiamen, the estimate was 49.2 thousand in the spring of 2014, while the ratio between the effective recruitment amount and the brood stock amount was 3.92, and the survival rate of released shrimp larvae was 1.88‱.

Estimation of Biomass and Replenishment Effect Assessment
Jiag et al. evaluated the recaptured yield and abundance of spawners originating from the hatchery-released juveniles of Nibea albiflora in Xiangshan Bay, with a mark-recapture experiment [11]. The natural and fishing mortalities of the hatchery-released individuals were 0.51/a and 1.31/a, respectively. The effectiveness of N. albiflora stock replenishment is strongly dependent on the level of fishing effort, and an appropriate reduction in the fishing effort would benefit both the recapture yield and the abundance of spawners originating from hatchery-released juvenile N. albiflora. In this study, the survival rate of the larvae of P. penicillatus stock replenishment released in Xiamen Bay from 2014 to 2016 was 1.88‱. The result showed that only a very small number of individuals could survive during the release process. This means that although billions of prawn larvae were released, the amount of resources in the sea was not well restored.
Wang used historical production to evaluate the release effect of P. penicillatus in Luoyuan Bay and believed that the recapture rate of release was between 5% and 7% [16]. An in vitro listing method evaluation of the reproduction and release effect of P. penicillatus in Sanduao waters was conducted [21], and the recapture rates of the two times release were 9.3% and 20.45%, respectively. Lin et al. conducted multiple releases of prawns in Daya Bay waters. In vitro plastic labeling was carried out, and the results showed that the recovery rate of P. penicillatus was 5.0% [43]. The abovementioned studies seem to have very different results from the results of this article. We speculate that the main reason may be twofold, as follows: the size of the released juveniles is different, and the release process is different. We discussed the two release procedures and conditions of Guan et al. [21]. The average body lengths of the shrimp larvae were 7 cm and 10.72 cm, respectively, and they were temporarily raised before release. Obviously, the use of larger juveniles and good temporary rearing is beneficial [21]. In contrast, regarding the release process and conditions of Lin et al., the body length of juvenile shrimp is between 1 and 5.7 cm, which also seems to represent a small juvenile shrimp [36]. Compared with our conditions, the average body length of the shrimp released in Xiamen Bay, from 2014 to 2016, was 1 cm. In short, the two conditions, including small individuals and no temporary rearing, may be the reason behind the higher mortality and lower fishing rate.
At that time, due to the small size of juvenile shrimp, it was not easy to use the method of marking in vitro, and the damage of the body surface was expected to cause many deaths, so this study did not use the method of marking and recapture evaluation. To overcome this problem, molecular labeling and isotope labeling seem to be more suitable methods for smaller release species. This is also an expected topic of ongoing research in the future. As we could not use the method of marking and recapture, we chose traditional fisheries resource evaluation. In this study, we used mathematical modeling to separate the two sources of recruitments. The highlight of this approach is that it confirms the effect of the release. Furthermore, the data used can basically be obtained through years of field surveys.

Conservation Suggestions
At this moment, the scientific basis for marine restocking, stock replenishment, and sea ranching continues to advance rapidly. The government assesses the possible contribution of these methods to fisheries management goals before designing enhancements and needs to be judged as having good potential, with the plan implemented effectively and responsibly [44]. It is widely recognized that there is a need to reduce fishing efforts and restore habitats to increase the resilience of capture fisheries [44]. However, the application of this technology still has a long way to go before an integrated management system that successfully solves all biological, ecological, social, cultural, and economic problems is in place [5]. The main challenges include determining when and where to use interventions to add value to management, integrating these initiatives with institutions and fisheries management regimes, monitoring the success of interventions, considering the cost effectiveness of culturing farmed juveniles, and releasing them into the wild, so that they survive in high proportions [5,44]. Marine ranching, which is considered a sustainable fisheries mode that has advantages for the ecosystem approach to fisheries [45], aquaculture, and capture-based aquaculture, is rapidly growing in China [45]. Habitat restoration and construction technology, stock replenishment, and the behavioral control of fisheries resources are some of the important methods of marine ranching [23,[44][45][46].
Due to the small body length of the released juvenile P. penicillatus, their ability to resist impact is low. In addition, there are many predators in the environment. Ultimately, the number that can survive is limited. Therefore, in order to improve the survival rate of juvenile shrimps, corresponding implementations should be taken, including standard crude culture, reduction in stress during transportation, and temporary breeding. These implementations can improve the resources of P. penicillatus for the purpose of restoration in Xiamen Bay.
In addition, the factor of the habitat environment is also very important. This factor will have a great impact on the survival and growth of the prawn. With the development of the economy in Xiamen Bay, there are a lot of human-made disturbances, such as estuary pollution, marine garbage, seabed sand mining, and artificial land reclamation, etc. Reducing the damage to the habitat of the prawns requires the relevant departments to formulate corresponding implementations. For this, strengthening sewage treatment and monitoring, resolutely combating illegal sand mining, and salvaging seabed garbage are recommended, among others [46]. Improving the living environment and optimizing the habitat of P. penicillatus are meaningful implementations to increase the survival rate of the released prawns.
Illegal fishing is also one of the influencing factors [46]. Jiag et al. mentioned that the effectiveness of N. albiflora stock replenishment is strongly dependent on the level of fishing effort, and an appropriate reduction in the fishing effort would benefit both the recapture yield and the abundance of spawners originating from hatchery-released juvenile N. albiflora [11]. Thus, an appropriate reduction in the fishing effort is essential to improve its effectiveness. Therefore, we make the following recommendations:

1.
Vigorously crack down on illegal fishing activities, including the confiscation and destruction of illegal fishing gear, and strengthen the punishment of offenders.

2.
Conduct legal and ecological education to make relevant fishermen understand the importance of ecology, the environment, and resources.

3.
Increase the number of patrols, regular patrols, and random inspections, so that illegal fishing personnel have nowhere to hide, and increase the cost of illegal fishing crimes.

4.
Organize professionals from relevant scientific research institutions to regularly carry out publicity and education on marine biological resources and environmental protection, enter primary and secondary schools or communities, and improve the ecological awareness of ordinary people regarding protecting and caring for the ocean.
In the future, after the restoration of P. penicillatus in Xiamen Bay, the genetic diversity of the species must be ensured on a sustainable basis. In order to prevent all the released individuals from being close relatives, the gene integrity of Penaeus prawns in natural seas is destroyed. Therefore, in the process of releasing, it is necessary to collect and hatch eggs of different gene groups. Finally, enhancements enter into complex fisheries systems. It is crucial to consider the fisheries system, broad objectives for management, and the full range of management options when assessing the potential for developing and using enhancements [44,45].

Conclusions
Since 2010, approximately 700 million juvenile shrimp (P. penicillatus) have been released into Xiamen Bay, Fujian Province, China, each year, through stock replenishment programs. The biological aspect consists of length-weight analysis, growth parameters, level of exploitation (natural mortality rate, fishing mortality rate, and total mortality rate), and an evaluation of stock replenishment and restocking was conducted. An assessment model for effect evaluation, according to Baranov's catch equation, to separate the amount of initial recruitment, was constructed using survey data from 2014 to 2017. The survival rate of the released shrimp larvae, 1.88‱, seems to be very unsatisfactory. Obviously, the restocking effect is lower and the programs need to be improved. To improve the restocking effect, the replenishment performance should be adjusted to reduce the mortality rate and increase its release effectiveness. Therefore, corresponding implementations are recommended, including standard extensive culture, reduction in stress during transportation, and temporary culture. Of course, the results of this evaluation are unsatisfactory. We believe that the main limitations of the created evaluation model are the consequence of certain assumptions that were made. Therefore, identifying how to improve our ability to carry out successful evaluations is expected to be the focus of efforts in the future.