Next Article in Journal
Effect of Water Area and Waterweed Coverage on the Growth of Pond-Reared Eriocheir sinensis
Next Article in Special Issue
Otoliths and Their Applications in Fishery Science
Previous Article in Journal
Decadal Changes in the Annual Growth of Two Commercial Ommastrephid Species in the Pacific Ocean
Previous Article in Special Issue
Timing of Increment Formation in Atlantic Bluefin Tuna (Thunnus thynnus) Otoliths
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Length–Weight Relationships, Growth Models of Two Croakers (Pennahia macrocephalus and Atrobucca nibe) off Taiwan and Growth Performance Indices of Related Species

1
Graduate Institute of Marine Affairs, National Sun Yat-sen University, Kaohsiung 804, Taiwan
2
Coastal and Offshore Resources Research Center, Fisheries Research Institute, Council of Agriculture, Kaohsiung 806, Taiwan
3
Department of Applied Mathematics, Tunghai University, Taichung 974, Taiwan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fishes 2022, 7(5), 281; https://doi.org/10.3390/fishes7050281
Submission received: 3 September 2022 / Revised: 2 October 2022 / Accepted: 7 October 2022 / Published: 11 October 2022
(This article belongs to the Special Issue Otoliths and Their Applications in Fishery Science)

Abstract

:
Information on age and growth is essential to modern stock assessment and the development of management plans for fish resources. To provide quality otolith-based estimates of growth parameters, this study performed five types of analyses on the two important croakers that were under high fishing pressure in southwestern Taiwan: Pennahia macrocephalus (big-head pennah croaker) and Atrobucca nibe (blackmouth croaker): (1) Estimation of length–weight relationships (LWR) with discussion on the differences with previous studies; (2) validation of the periodicity of ring formation using edge analysis; (3) examination of three age determination methods (integral, quartile and back-calculation methods) and selection of the most appropriate one using a k-fold cross-validation simulation; (4) determination of the representative growth models from four candidate models using a multimodel inference approach; and, (5) compilation of growth parameters for all Pennahia and Atrobucca species published globally for reviewing the clusters of estimates using auximetric plots of logged growth parameters. The study observed that features of samples affected the LWR estimates. Edge analysis supported the growth rings were formed annually, and the cross-validation study supported the quartile method (age was determined as the number of opaque bands on otolith plus the quartile of the width of the marginal translucent band) provided more appropriate estimates of age. The multimodel inference approach suggested the von Bertalanffy growth model as the optimal model for P. macrocephalus and logistic growth model for A. nibe, with asymptotic lengths and relative growth rates of 18.0 cm TL and 0.789 year−1 and 55.21 cm, 0.374 year−1, respectively. Auximetric plots of global estimates showed a downward trend with clusters by species. Growth rates of the two species were higher than in previous studies using the same aging structure (otolith) and from similar locations conducted a decade ago, suggesting a possible effect of increased fishing pressure and the need to establish a management framework. This study adds updated information to the global literature and provides an overview of growth parameters for the two important croakers.

1. Introduction

Information on age and growth is essential to modern stock assessment and the development of management plans for fish resources; it has many functions such as converting input catch estimates from biomass to numbers and converting output numbers into biomass [1,2] and can affect the estimation of natural mortality and the mean age at maturity [3,4,5]. Establishing appropriate aging procedures and selecting representative growth models are thus important steps in developing stock assessments [1,6]. Even if the growth curve has been estimated previously, reviewing the change in key parameters by a timely repeat of estimation is also beneficial to monitoring the stock status under fishing pressure and climate change [7,8,9,10,11].
Several approaches to estimate the growth curve have been used, including length-frequency analysis and growth bands analysis of fin rays, vertebrae, otoliths and scales. Otoliths are calcium carbonate structures formed inside the inner ears and grow with the deposition of continuous calcium carbonate layers, which can respond to both daily and seasonal changes. For teleost fishes, otoliths analysis is the most commonly used and reliable approach to age determination [12,13,14]. Age estimation using otoliths however requires high precision in laboratory techniques and estimating considerations, including validation of the periodicity of growth-rings formation, selection of age determination method (expression of age in integral count of growth rings or estimated value), and determination of representative growth curves.
Accuracy of age estimates is crucial in deriving age-based population parameters for making the right management decisions [15,16,17]. Edge analysis and marginal increment analysis (MIA) have been the most commonly used method to verify the increment periodicity of the otolith [18,19,20,21]. Edge analysis examines the timing of the translucent zone formation, and the ages are considered validated if only one opaque and one translucent band is formed annually [22,23,24,25]. The method is relatively easy to conduct compared to MIA, which plots the proportion of completion of the outermost increment against the month of capture to examine if it shows a yearly sinusoidal cycle [13,26,27] and can be used alone for the validation of annual periodicity of growth increment [18,24].
When determining the age under the microscope, several methods can be used for teleost fishes, including the commonly used integral counting method (integral method) where the age is the count of opaque bands, ignoring the growth within a year reflected in the width of the outmost translucent bands [28,29,30,31], and the back-calculate method where the age was determined as the proportion of the hypothesized hatching date to the sampling date plus the number of opaque bands [32,33]. Daily increments count has also been used [34]. Another method that can be explored is the quartile of marginal growth band method (quartile method) [35] to enhance the integral counting method by adding a quartile of the width of the marginal translucent band to the age. Ages estimates resulting from applying these methods are different and the accuracy needs to be investigated in advance.
The von Bertalanffy growth model (VBGM) has been widely used for fitting growth curves to the estimated age and the length data, but the standard VBGM is not always appropriate for all species since the ontogenetic changes in growth rates vary by species, and so alternative models should be explored [6,36]. The multimodel inference approach [37,38] is a useful way to select the most representative growth model from a suite of candidates for the species [27,39,40].
Plotting the estimated growth parameters together with those of other studies of the same species or related species can provide an overview of the similarities and differences between the current study and other global studies, and a chance to explore the differences [41,42,43]. Pauly and Munro [44] proposed a growth performance index based on logarithmized mean asymptotic length and relative growth rate obtained from the growth model, which allowed the identification of the bias in growth parameter estimates within a species while also revealing the difference between species. The auximetric plot created using these two parameters can help to visualize the clusters formed by the same species [43,45]. It has been applied widely to compare estimated growth parameters from different studies [27,46,47,48].
Sciaenidae includes 70 genera with around 282 species distributed in the Atlantic, Indian, and Pacific Oceans [49] and are important to coastal fisheries of many countries. There were 23 species in 12 genera of Sciaenidae recorded in Taiwan waters [50]. The high economic value of the Sciaenidae and its characteristics of mass migration and concentration during the breeding season make it easy to be caught in large quantities. Sciaenidae species are an important target to the trawl fishery of Taiwan. Among them, Pennahia macrocephalus ([51]; big-head pennah croaker) and Atrobucca nibe ([52]; blackmouth croaker) are two of the most abundant commercial species in southwestern Taiwan (Figure 1). They are widely distributed across Indo-Pacific waters with depth between 20 to 200 m, feeds on small invertebrates and crustaceans. Total landing of Pennahia sp. once reached a peak of 8.6 thousand tons in 1996 and declined to 852 tons in 2015, while that of A. nibe was 2.5 thousand tons in 1993 and declined to less than 500 tons after 2011 [50]. Lack of logbook from coastal fisheries due to exemption of submission responsibility [53] has hindered the assessment of the stock using catch data; studies of biological parameters are considered as alternatives for providing supplemental information on the stock.
As indicated above, we suppose that different age determination methods may affect the results of growth estimates, the most optimal growth model of these two species may not be the VBGM, and same fish species may have a relatively similar growth performance [44]. Therefore, five types of analyses were performed on the two species in this study: (1) estimation of length–weight relationships (LWR) which is also a key element in the examination of fish biology and contains valuable information to assess the general health condition of fish species [54,55,56]; (2) validation of the periodicity of ring formation using edge analysis; (3) examination of three age determination methods and selection of the most appropriate one using a k-fold cross-validation simulation study [57,58]; (4) determination of the best fitting growth models from four candidate models; and, (5) compilation of growth parameters for all Pennahia and Atrobucca species published globally for reviewing the clusters (groups) of estimates using auximetric plots. Most growth studies on the two species were conducted a decade ago and the majority of them were published in grey literature or in Chinese. This study adds updated information to the global literature and provides an overview of growth parameters for the two important croakers.

2. Materials and Methods

2.1. Data Sources and Estimation of Length-Weight Relationships (LWR)

Biological data: total length (TL) in cm, body weight (BW) in g, sex, and otolith samples of P. macrocephalus in southwestern Taiwan waters were collected monthly from July to October 2018 and February to August 2021 in the fish markets of Zihquan and Tungkang (Figure 2). Biological data of A. nibe were collected monthly from July 2017 to August 2018 at the fish markets of, or from fishing boats based in, Zihquan, Tungkang, Fangliao and Liuqiu (Figure 2). Excluding samples abandoned due to damage caused during otolith processing or disagreements on the number of ring marks between two age readings, a total of 359 P. macrocephalus and 386 A. nibe samples were used in this study.
The length–weight relationship (LWR) for each species was estimated using a power function of the form B W = a × T L b where a is the coefficient of the power function, and b is the allometric coefficient. Parameters a and b were estimated with a linear regression of ln(BW) against ln(TL). Size distribution and LWR by sex for each species were also examined. The difference of the parameters from the linear regression fits between sexes for both species were tested using ANOVA [59].

2.2. Otolith Processing, Periodicity of Increment and Age Determination

Otoliths (sagittae, lapilli, and asterisci, Supplementary Figure S1) were extracted from the fish samples. Sagittae, the largest of the three otoliths, are commonly used for estimating fish age [12] and were used for this study. They were firstly washed with water and air-dried for at least 24 h. Next, the otoliths were prepared for transverse sections by embedding the whole otolith in epoxy resin and epoxy hardener mixture (with the ratio 10:0.95), placing them at 60 °C for 2 h, and sectioning transversely with a low-speed diamond wheel saw (approximately 500 μm thick). Finally, the otolith sections were polished with silicone abrasive grinding paper (1200 to 4000 grit) until the rings became clear. The otoliths were examined using a microscope with transmitted light, by which the periods of fast growth accreted otolith material are translucent whereas it is opaque during slow growth periods, which is different from the observation viewed under a dissecting microscope with reflected light [60]. A translucent zone and subsequent opaque zone are referred to as an ‘increment’. Ring marks on the otoliths were counted and measured by at least two different readers or three times by the same reader. The results were accepted when the readings matched; otherwise, the results were abandoned and were not included in this study.
Edge analysis was used to examine the periodicity of increment formation (the timing of the translucent zone formation) and was calculated through the monthly proportion of the translucent zone on the outer edge of the sectioned otolith [23,24]. Subsamples for edge analysis were selected by month from the most abundant age groups, with efforts to restrict the number of age groups [13] (n = 296 for P. macrocephalus with age 2–4, TL 13.7–19.6 cm and 186 for A. nibe with age 4–6, TL 23.2–53.8 cm).
Three methods for determining the annual increments were conducted and evaluated for selecting the optimal one. (1) Integral counting method (integral method), where the opaque bands were counted from the nucleus to the edge of the otolith. This method was commonly used for many fish species, but since the marginal translucent band was ignored, this method may underestimate the actual age of the fish (Figure 3a). (2) Quartile of marginal growth band method (quartile method). The age was determined in the same way as the integral counting method when the edge of the otolith is an opaque band. If the edge of the otolith is a translucent band, the width of the marginal translucent band was compared to the width of the previous complete annulus to the approximate value of 0.25, 0.5 or 0.75; and the age was determined as the quartile plus the number of opaque bands (Figure 3b) [35]. (3) Back-calculation method. The age was determined as the proportion of the hypothesized hatching date to the sampling date plus the number of opaque bands [32].
A k-fold cross-validation simulation study [57,58] was performed to evaluate the goodness of fit of the growth model to the age-length data estimated by each of the three methods. The most commonly used VBGM was used as the reference model. The age-length data, with length data and age estimates from the three methods, were firstly split into four groups based on the results of the integral counting method (e.g., ages 0–1, age 2, age 3, and ages 4–5). Data within each age group were then randomly divided into five sub-groups, or folds, of approximately equal size. Each fold contains data from all age groups. The first fold (across the four age groups) was treated as a validation set, and the VBGM was fit on the remaining four folds (training set). The mean square error (MSE) was calculated from the observed length of the validation set, and the predicted length was estimated from the VBGM developed using the training set. Afterwards, the other folds were treated as validation sets, and the VBGM-developing and MSE-estimation were repeated. This loop was repeated 100 times to finally obtain the mean and standard deviation (sd) of the MSE for each type of age determination method.

2.3. Growth Parameters Estimation

An information-theoretic, multi-model inference approach was used to determine the optimal growth model for both species [37,38]. The candidate growth models examined include non-sigmoidal VBGM and three sigmoidal models: Gompertz, Logistic, and Richards models [61,62,63,64]. The first three are the most common three-parameter models [65], and the last is a four-parameter model.
Von   Bertalanffy   model :   L t = L 1 e K t t 0
Gompertz   model :   L t = L e e K t t 0
Logistic   model :   L t = L 1 + e K t t 0 1
Richards   model :   L t = L 1 + 1 p e K t t 0 p
where Lt is the length (TL in cm) at age t, L is the mean asymptotic length, p is a dimensionless parameter of Richards model, and t0 is the theoretical age when length equals zero for the VBGM or the age at the inflection point of the growth curve for the other models. Parameter k (year−1) is a relative growth rate, describing how quickly the asymptotic length is approached or the rate of exponential decrease of the relative growth rate with age [37,66].
The length-at-age data (combined and sex-specific) were fitted to the four candidate models using non-linear least squares in R version 4.1.2 [67]. The small-sample bias-corrected form (AICc) of the AIC (Akaike Information Criterion) was used for model selection [65,68,69]. The model with the smallest AIC value (AICc, min) was selected as the ‘best’ among the models tested. AIC differences Δ = AICc, min − AICc, i were computed overall candidate models i. The Akaike weight, wi, of each model was then calculated using these differences to quantify the plausibility, which is considered as the weight of evidence in favor of model i being the best of the available set of models [37,69].
w i = e 0.5 Δ i i = 1 4 e 0.5 Δ i

2.4. Global Estimate Compilation

A literature (e.g., books, peer-reviewed articles and project reports) search was conducted on the estimation of growth parameters for Pennahia and Atrobucca species; however, only four Pennahia and two Atrobucca species (i.e., P. macrocephalus, P. argentata, P. pawak, P. anea, A. nibe and A. alcocki) were available. Global estimates of growth parameters of the six species were compiled from the literature in Table 1 and Supplementary Table S1. Estimates of k were in year−1; estimates of L were in cm in either standard length (SL), fork length (FL), total length (TL), body length (BL) or not available (NA). Most of the L were in TL; thus, the estimates were converted to TL using the formulae listed below [70,71,72] or assumed to be TL for the cases of NA.
TL = 0.651 + 1.162 SL (R2 = 0.997, for P. macrocephalus)
TL = −0.37531 + 1.225 SL (R2 = 0.965, for P. argentata)
TL = FL (For P. anea)
Growth parameters were estimated using different methods: the aging of otoliths (sagittae) or scales, length-frequency data (LFD), or unknown method. Growth parameters for the two species in this study were compiled into the global estimates table (62 estimates in total) for classification analysis. All the parameters were from VBGM, except for the records coded as p_1 and nc_2, which were Logistic.
The two variables, log(L) and −log(k), of the growth performance index formula, φ = 2log(L) + log(k) [44], were used to identify the clusters of growth estimates from the global summary using the analysis of the auximetric plot [43,45,46]. Ellipses of 95% confidence intervals were also produced in the plot except for P. pawak with three records and A. alcocki with one record only.

3. Results

3.1. Samples and Length-Weight Relationships (LWR)

For P. macrocephalus, 359 individuals with total lengths and bodyweight ranging from 12.20 to 20.13 cm and from 21.70 to 104.82 g, respectively, were examined (Supplementary Table S2 with raw data and summary statistics). Males ranged from 12.20 to 20.13 cm TL and 21.70 to 104.70 g (n = 207) and females from 13.45 to 19.62 cm TL and 27.17 to 104.82 g (n = 152). The ANOVA test results indicated that there is no significant difference of the LWR parameters between sexes (p = 0.095, df = 355). Therefore, the sex combined LWR equation is presented as BW = 1.122 × 10−3·TL3.061 (R2 = 0.892) (Figure 4) (or BW = 2.770 × 10−2 SL2.937). To examine if the coefficients have been biased by possible overweighting by more samples from the breeding season (June to September), an additional LWR estimation was performed using limited sub-samples: all data in a month were used when the monthly sample size was <30, or 30 fish were randomly selected when the monthly sample size was >30. The resulted average was b = 3.122 (n = 10), suggesting that reducing sample size from the breeding season does not reduce the estimation of b and that the current estimates using all sample data are representative.
For A. nibe, 386 individuals with total lengths and bodyweight ranging from 17.02 to 56.6 cm and from 46.30 to 1941.60 g, respectively, were examined (Supplementary Table S3 with raw data and summary statistics). Males ranged from 17.01 to 52.70 cm TL and from 46.30 to 1235.07 g (n = 152) and females from 18.65 to 56.60 cm TL and from 78.37 to 1941.60 g (n = 234). The ANOVA test results indicated that there is no significant difference of the LWR parameters between sexes (p = 0.156, df = 382). The length–weight relationship equation is thus presented as BW = 1.448 × 10−2·TL2.886 (R2 = 0.985) for sex combined.

3.2. Periodicity of Increment Formation—Edge Analysis

Monthly percent occurrence of the translucent band at the edge for the two croakers is depicted in Figure 5. For P. macrocephalus, we were not able to obtain samples for November to January since it did not coincide with the active fishing season. However, the trend from the lowest proportion of translucent bands in February to the high level from June to October could still be observed. Although there was data deficiency during the winter period, the apparent trend suggested that the translucent bands started to form once a year near the end of winter and the early spring.
On the other hand, A. nibe can be caught throughout the year. The result of the edge analysis showed that the proportion of translucent bands started to increase from January to September (100%), with a drop in May (likely caused by sampling issues) followed by a sharp decline to 0% in December. The clear trend from edge analysis suggested that the translucent bands started to form once a year near the end of winter and during the spring.

3.3. Age Determination Methods and Growth Parameters Estimation

For P. macrocephalus, mean ± sd of MSE from the k-folds cross-validation runs for the three age determination methods was 1.722 ± 0.010 for the integral method, 1.572 ± 0.048 for the quartile method, and 1.730 ± 0.011 for the back-calculation method. The quartile method has the lowest mean MSE compared to the other two methods and was considered the optimal method. A similar conclusion was obtained for A. nibe that the quartile method has the lowest mean MSE, with mean ± sd of MSE of 31.233 ± 5.005 for the integral method, 30.538 ± 4.011 for the quartile method, and 33.783 ± 1.760 for the back-calculation method. The quartile method was thus selected as the age determination method for both species.
For P. macrocephalus, the Richards model could not converge in the growth fitting analysis. The growth model that fits the length and age data best was VBGM with AICc differences ∆ < 2 (indicating substantial support as the best model) [38] and Akaike weight w (the expected weight of evidence in favor of the model being the best among the four models) of 47% with sex combined (Figure 6, Table 2). The growth parameters of VBGM (L, k, t0) were: 18.0 cm TL, 0.789 year−1, and −0.872 year for sex-combined; 18.673 cm, 0.495 year−1, −2.095 year for female; and, 17.864 cm, 0.892 year−1, −0.571 year for male.
For A. nibe, the best growth model was logistic with AICc differences ∆ < 2 and Akaike weight w of 75% with sex combined (Figure 7, Table 2). The growth parameters of logistic (L, k, t0) were: 53.11 cm, 0.39 year−1, and 2.96 year for sex-combined; 52.14 cm, 0.45 year−1, 2.80 year for female; and, 55.31 cm, 0.32 year−1, 3.37 year for male.
The selection of the best model for each species was based on the quality of statistical fit assuming that the data could represent the potential species-specific growth pattern. Table 2 provides estimates of all the tested models for optional choices when considering the representativeness of the data set and the theoretical growth pattern from biological details of the fish [38,100].

3.4. Global Estimates Compilation

Few studies on P. macrocephalus and A. nibe were found in the literature. Especially for A. nibe, many were reported in grey literature and were cited in academic publications without detailed information (e.g., many were cited in Pauly [92], and information on the data sources are not available online). Due to the lack of studies in the two targeted species, studies from the same genus documented previously were also included in the analyses. A total of 62 growth records were obtained covering six species from ten regions with three estimation methods, including otolith, scale and length frequency (Table 1).
For Pennahia species, the age ranged from 0 to 10 years, and the growth parameters of VBGM ranged from 15.15–45.20 cm TL for L and 0.181–1.94 year−1 for k. The growth performance index was 2.10 to 3.16 (mean = 2.61), and for P. macrocephalus was between 2.10 to 2.77 (mean = 2.43). For Atrobucca species, the age ranged from 1 to 13 years, and the growth parameters ranged from 37.10–68.15 cm TL for L and 0.116–0.523 year−1 for k. The growth performance index was 2.58–3.05 (mean = 2.75), and for A. nibe was in the same range (mean = 2.76).
The two variables, log(L) and log(k), of the growth performance index formula and the results of the growth performance index, were shown in the auximetric plots (Figure 8). The same species tend to be in the same cluster with similar φ.

4. Discussion

4.1. Length-Weight Relationships

In order to compare results from this study with past studies, the LWR in TL was transformed to a LWR in SL for P. macrocephalus using the length-length relationship [70]. Divergencies were noted among the LWR estimates in Table 3. Although the reasons for the divergencies have not been fully explored, in addition to possible sampling errors in data collection methods [101], the following differences in data handling may, in part, account for the observed differences. (1) Different logarithmic functions are used for the LWR parameter estimation. Both base-10 logarithmic (Log10) [54,70] and natural logarithmic (Ln) functions [59,102,103,104] have been used to linearize the LWR equation for regression estimation, and most of the studies did not specify the logarithmic function they used. The two functions will result in different parameter values while the b parameter remains the same. (2) Uneven distribution of the samples by fish size, such as oversampling of medium-sized fish, can affect the regression fitness. For example, large fish were less fitted in the LWR of Tasi (Figure 19 in reference [34]) (Table 3), which might cause an overestimation of the allometric coefficient (b value). (3) Uneven distribution of sampling size among months, e.g., over-weighted samples were from the months when the fish grow faster or slower, could cause estimation differences since LWR may change seasonally with different levels of gonad maturity and other parameters [105]. Nhuận and Khuong [73] (Table 3) showed that b values estimated in the breeding season (third season) tended to be higher than in other seasons. For this study, most of the samples were from the breeding season, and thus, the b value was likely overestimated to some extent.
For A. nibe, differences in LWR were also observed (Table 3), which might result from several possible reasons, including seasonal or environmental changes, as proposed in Schneider [111] and De Giosa [112] and the abovementioned data handling issues. Sample distributions by fish size and month were evenly distributed in this study, so the LWR estimation was considered representative.

4.2. Periodicity of Increment Formation—Edge Analysis

Both edge analysis and marginal increment analysis (MIA) are commonly used methods to verify the increment periodicity of the otolith [20,21], with edge analysis being a qualitative approach and MIA a quantitative approach [40]. While MIA is generally considered more robust than edge analysis, edge analysis (or marginal analysis) only requires recording the otolith margin type, thus is relatively affordable both in terms of equipment and time [26] and was chosen in this study to examine the periodicity of increment formation. For P. macrocephalus, the occurrence of the translucent bands was lowest in February, increased in March and April, and then jumped to a high level from June to October (Figure 5), suggesting that the translucent bands were formed once a year near the end of winter and the early spring. This pattern was supported by an additional MIA (Supplementary Figure S2), which showed a continuously increasing trend of marginal increment from February to October.
The spawning season of this fish is from April to October (also the main fishing season), and the gonadosomatic index (GSI) value declines significantly from September to October [113]. The fish migrate to Taiwan mainly for spawning, and many large spawning fish species inhabit deeper depths after spawning (e.g., white-mouth croaker [114]) and thus are more difficult to catch by fisheries. Though sampling could not be performed year-round, this study assumed that most fish were hatched before or in October and that the lowest proportion of the translucent bands might be in the month without fish samples (November to January). The MIA of Attaqi [32] also showed a single sinusoidal cycle in which the highest appeared from August to September and decreased sharply in October. Results from this study and Attaqi [32] both support ring marks being formed once a year.
On the other hand, consecutive monthly samples were available for A. nibe, and the edge analysis showed a continuous increase in the occurrence of the translucent band from the lowest month in December to the highest month in September (Figure 5). The single sinusoidal cycle pattern also supported the annual formation of the ring mark.
The edge analysis and MIA are valuable tools for the validation of age estimates, but both have their limitations (refer to Campana [13] for a detailed discussion on the topic). Therefore, although the current analyses supported annual formation of the ring mark for the two species, further analyses are still needed for definite validations; parameters that would possibly affect ring formation (e.g., fluctuations of environmental factors and timing of high feeding sources) are also worth further investigation.

4.3. Age Determination Methods

A growth ring is a part of an otolith that consists of one translucent zone and one opaque zone [115]. When the annual formation of the growth ring is validated, the age is usually determined by the number of ring counts [20,116,117]. However, this method ignores valuable information contained within a year’s growth. Although age estimates for older fish may be best approximated by integer ages because the edge proportion is usually hard to identify for older fish (as seen in Figure 7), information on decimal ages can be gathered in younger fish. The quartile method could retain the lost information and be almost as fast as the integral method when conducting. Cross-validation simulation demonstrated that the quartile method could provide more appropriate estimates of age than the integral method.
The back-calculation method has been applied in many studies [32,33]. It needs an assumption of the birthdate of the fish. Vanderkooy [33] described the age from this method as the biological age, derived from the information of capture month (sampling month), number of annuli, and accepted birthdate estimate. This method could provide better distribution of age estimates when developing a growth model than the integral method if the study has representative samples of the sizes closest to the breeding and juvenile stages to determine the correct birthdate. In this regard, this method needs beforehand studies on fish reproduction or daily-rings counts to provide birthdate information. For P. macrocephalus, Attaqi [32] set 1st August as the hypothesized hatching date based on Wang [118], but the study from He [113] suggested that the spawning season could be as long as seven months since April. In this case, the difficulty in determining the birthdate would cause uncertainty in the estimate of age. Comparatively, the quartile method is more straightforward in interpreting the age of fish otoliths with both efficiency and accuracy.
Only Tsai [34] estimated the growth parameters of P. macrocephalus based on “daily rings”, while others used annuli. The study suggested one ring formed every two days (“two-day ring”) in the otolith based on the comparisons of the back-calculated birthdate and the GSI trend, and so the maximum age of this species was estimated to be 2.2 years. Besides lacking validations of the periodicity of increment formation and the first increment [27], the study also needs to consider the possibility that nondaily deposition can occur in suboptimal or extreme environmental conditions [119]. With these issues being addressed, daily rings for growth studies could provide more accurate estimates than the three methods described here, although the high cost of time and human resources need to be considered.

4.4. Growth Parameter Estimates—Global and Local

The global estimates of growth parameters for the Sciaenidae genera, together with those of this study, were collated and plotted for an overview (Figure 8) for the first time. The overall trend is downwards, from P. anea in the top-right corner to P. macrocephalus, P. pawak, P. argentata, and A. nibe in the bottom-left corner with comparatively higher L and smaller k.
Estimates of each species were clustered as a group, demonstrating that each species has a relatively similar growth performance [44], which is also noted for tuna species [43,45]. The elongated clusters were also found in catfish (Ariidae) and sardines (Clupeidae), despite the different sampling regions or aging methods [41,42]. Estimates of the same species but outside the group might have specific reasons or need further investigation [27,47,48]. For example, estimates of a_1 and a_3 of P. argentata are far away from the cluster of the species, and further examining the estimates could find that the two estimates were specifically on smaller fish (SL ≤ 26.59 cm for a_1 and SL ≤ 24.47 cm for a_3) that have smaller L but much higher k. As demonstrated in other studies, growth parameter estimates would be significantly affected by the lack of very young or old individuals [120], and the lack of older fish might underestimate asymptotic length [121].
The type of aging structure has been observed to significantly affect the estimation results [27,48]. However, although otolith, scale or length-frequency data were used for growth studies on the species in Table 1, differences in the estimates using different structures were not obvious for unknown reasons.
There were not many growth studies for P. macrocephalus. As explained earlier, Tsai’s [34] study (m_4) was likely biased. Excluding this estimate, the estimates from the South China Sea (m_1 to m_3) departed from those from Taiwan waters (m_5 and mc_1–mc_3), suggesting possible regional effects (different environmental productivity or fishing pressure) on the estimation (Figure 8). Diverse results using different age determination methods in this study were observed in the figure. The recommended estimate with higher precision (using quartile method, mc_2) departs from Attaqi’s [32] estimate that studied the fish in a similar location with smaller L and higher k. In addition to many possible technical reasons in the estimation [27], the phenomenon might affect high fishing pressure. The data used in Attaqi [32] was from 2010–2011, about a decade before the data used for this study (2018–2021). Total landings of the category of “white mouth croaker” (P. macrocephalus composed the majority of the catch) in the market had declined drastically from 800 tons in 2011 to <200 tons after 2018 [122]. High fishing pressure can lead to less competition for food and space in the fish population and thus produce higher growth coefficients [7,8,9,123]. On the other hand, high fishing pressure can cause the loss of larger individuals and a decrease of relative abundance of target species which have the capacity to grow to large size, and thus is attributed to a smaller estimate of asymptotic length [7,124]. The same inference was drawn in Yi [80] for P. pawak in which high fishing pressure contributed to the decline of L from 24.15 cm in 2008–2009 to 22.05 cm in 2018–2019, a decade later. These findings suggested the necessity to establish a management framework on the species to reduce fishing pressures.
For A. nibe, the growth estimates of VBGM (nc_2) from this study and from the previous studies (n_1 to n_18) formed a cluster, and all were within the 95% confidence interval. The best estimate in this study was the logistic model (nc_1); however, this estimate deviated from the 95% confidence interval ellipse, which was likely due to the difference in growth models. Except for one of the estimates from a study in the Oman Sea (n_10), the rest were all from the northeast Pacific Ocean (Taiwan, Japan, Korea). However, since this species is not a highly migrating species, those fish in different waters might belong to different stocks, so the growth parameters were not comparable.
Otolith, scale and length-frequency data have been used to estimate the parameters. Estimates using length-frequency data tend to have lower k than estimates using hard structures for dolphinfish and flyingfishes [27,48]; however, this pattern has not been observed for A. nibe. Compared to the estimates from Taiwan using the same aging material of otolith (n_11, n_12, n_15–n_18), the k estimates of this study were much higher. The obvious difference between these studies was that the data of previous studies were conducted about a decade ago (1993 and 2008). Based on market data, the annual landing of the species has dropped 75% from the early 2010s to the late 2020s and hence fishing pressure might have contributed to the increase in k estimates. The other difference was that this study included more samples of older fish (ages 9 and 10), which could make the estimates of L higher and k lower [121]. However, care needs to be taken that the quality of this study might be affected by the lack of small fish samples (age 1), as suggested in the following.
Although the VBGM was selected as the best model for P. macrocephalus because of the lowest ∆AICc, the differences among models were not substantial (Table 2). The growth curves at Figure 6 also showed little differences in fit among the models where data are present; the differences occur mainly beyond the data (e.g., age 0). In addition, all models underestimate lengths of the oldest ages (4+). For A. nibe, the logistic model fits the age groups 2–8 better but diverted in age 0–1 and age 9+. Meanwhile, L estimates for both species fall beyond the range of observed data. These indicated the limitation of the lack of old and very young (0–1) age group samples in this study. To refine the growth estimates (especially the L estimates), further sample collections targeted at larger fish and very young (age 0–1) fish are needed in the future.

5. Conclusions

P. macrocephalus and A. nibe are among the most abundant commercial Sciaenidae species in southwestern Taiwan. Landings of the two species have dropped rapidly since 1998, presumably due to overfishing. Most studies on the estimation of growth parameters for the two species were conducted a decade ago, and the majority of them were published in the grey literature or in Chinese. To fill the information gap, this study provides reliable growth estimates and length–weight relationships (in Section 3.1 and Table 2) as updated information to the global literature and provides an overview of growth parameters for the two important croakers (in Table 1), which could be used for the development of a sound stock assessment on the resources. Comparison of the resulted growth parameters with those estimated a decade ago suggested that these two species have undergone high fishing pressures and induced consequent higher growth coefficients and smaller asymptotic length, and thus suggested the need for establishing a management framework to reduce fishing pressures.
Methodology analyses of this study suggested that the quartile of marginal growth band method could provide more appropriate estimates of age from otoliths than the other two methods. The results suggested the VBGM as the optimal growth model for P. macrocephalus and the logistic model for A. nibe, based on the samples collected. However, both results need further refinements by collecting more samples of large fish and very young fish. Finally, the auximetric plot of this study supported the argument that the same fish species has a similar growth performance and could be used as a quick tool for identifying outliers of growth estimates for further investigations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/fishes7050281/s1, Table S1: Detail information on the estimates of growth parameters for Pennahia and Atrobucca species from studies published globally; Table S2: Length and weight raw data for Pennahia macrocephalus; Table S3: Length and weight raw data for Atrobucca nibe; Figure S1: Images of otoliths (sagitta, lapillus, and asteriscus) of P. macrocephalus and A. nibe; Figure S2: The marginal increment ratio (MIR) of P. macrocephalus.

Author Contributions

Conceptualization, S.-C.H. and S.-K.C.; Data curation, S.-C.H. and J.-S.H.; Formal analysis, S.-C.H. and T.-L.Y.; Funding acquisition, S.-K.C. and J.-S.W.; Investigation, S.-C.H. and C.-C.L.; Methodology, S.-C.H., S.-K.C. and T.-L.Y.; Project administration, S.-K.C. and C.-C.L.; Resources, S.-K.C.; Software, S.-C.H. and T.-L.Y.; Supervision, S.-K.C. and J.-S.W.; Validation, J.-S.H.; Visualization, S.-C.H.; Writing—original draft, S.-C.H.; Writing—review & editing, S.-K.C. and C.-C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology (MOST), Taiwan, grant number MOST 106-2611-M-110-006, MOST 108-2611-M-110-002, MOST 109-2611-M-110-004.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data of this study is included in the report, the Supplementary Materials or may be made available by the authors.

Acknowledgments

The authors appreciate the technical support from Ying-Chu Chen in the otolith ring reading and the financial support of the Ministry of Science and Technology for the open access publication fee.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maunder, M.N.; Piner, K.R. Contemporary Fisheries Stock Assessment: Many Issues Still Remain. ICES J. Mar. Sci. 2015, 72, 7–18. [Google Scholar] [CrossRef] [Green Version]
  2. Francis, R.I.C.C.; Aires-da-Silva, A.M.; Maunder, M.N.; Schaefer, K.M.; Fuller, D.W. Estimating Fish Growth for Stock Assessments Using Both Age–Length and Tagging-Increment Data. Fish. Res. 2016, 180, 113–118. [Google Scholar] [CrossRef]
  3. Then, A.Y.; Hoenig, J.M.; Hall, N.G.; Hewitt, D.A. Handling editor: Ernesto Jardim Evaluating the Predictive Performance of Empirical Estimators of Natural Mortality Rate Using Information on over 200 Fish Species. ICES J. Mar. Sci. 2015, 72, 82–92. [Google Scholar] [CrossRef]
  4. Farley, J.; Eveson, P.; Krusic-Golub, K.; Sanchez, C.; Roupsard, F.; Nicol, S.; Leroy, B.; Smith, N.; Chang, S.-K. Project 35: Age, Growth and Maturity of Bigeye Tuna in the Western and Central Pacific Ocean; CSIRO: Rarotonga, Cook Islands, 2017; p. 51. [Google Scholar]
  5. McKechnie, S.; Pilling, G.; Hampton, J. Stock Assessment of Bigeye Tuna in the Western and Central Pacific Ocean. In Proceedings of the Thirteenth Regular Session of the Scientific Committee, Rarotonga, Cook Islands, 9–17 August 2017; p. 149. [Google Scholar]
  6. Maunder, M.N.; Crone, P.R.; Punt, A.E.; Valero, J.L.; Semmens, B.X. Growth: Theory, Estimation, and Application in Fishery Stock Assessment Models. Fish. Res. 2016, 180, 1–3. [Google Scholar] [CrossRef]
  7. Jennings, S.; Greenstreet, S.P.R.; Reynolds, J.D. Structural Change in an Exploited Fish Community: A Consequence of Differential Fishing Effects on Species with Contrasting Life Histories. J. Anim. Ecol. 1999, 68, 617–627. [Google Scholar] [CrossRef]
  8. Polacheck, T.; Eveson, J.P.; Laslett, G.M. Increase in Growth Rates of Southern Bluefin Tuna (Thunnus maccoyii) over Four Decades: 1960 to 2000. Can. J. Fish. Aquat. Sci. 2004, 61, 307–322. [Google Scholar] [CrossRef]
  9. Reznick, D.N.; Ghalambor, C.K. Can Commercial Fishing Cause Evolution? Answers from Guppies (Poecilia reticulata). Can. J. Fish. Aquat. Sci. 2005, 62, 791–801. [Google Scholar] [CrossRef]
  10. Kimura, D.K. Extending the von Bertalanffy Growth Model Using Explanatory Variables. Can. J. Fish. Aquat. Sci. 2008, 65, 13. [Google Scholar] [CrossRef]
  11. Lavin, C.P.; Gordó-Vilaseca, C.; Stephenson, F.; Shi, Z.; Costello, M.J. Warmer Temperature Decreases the Maximum Length of Six Species of Marine Fishes, Crustacean, and Squid in New Zealand. Environ. Biol. Fishes 2022. [Google Scholar] [CrossRef]
  12. Campana, S.E.; Neilson, J.D. Microstructure of Fish Otoliths. Can. J. Fish. Aquat. Sci. 1985, 42, 1014–1032. [Google Scholar] [CrossRef]
  13. Campana, S.E. Accuracy, Precision and Quality Control in Age Determination, Including a Review of the Use and Abuse of Age Validation Methods. J. Fish Biol. 2001, 59, 197–242. [Google Scholar] [CrossRef]
  14. Morat, F.; Wicquart, J.; Schiettekatte, N.M.D. Individual Back-Calculated Size-at-Age Based on Otoliths from Pacific Coral Reef Fish Species. Sci. Data 2020, 7, 370. [Google Scholar] [CrossRef] [PubMed]
  15. Lai, H.L.; Gunderson, D.R. Effects of Ageing Errors on Estimates of Growth, Mortality, and Yield per Recruit for Walleye Pollock (Theragra chlacogramma). Fish. Res. 1987, 5, 287–302. [Google Scholar]
  16. Beamish, R.J.; McFarlane, G.A. A. A Discussion of the Importance of Ageing Errors, and an Application to Walleye Pollock: The World’s Largest Fishery. In Recent Developments in Fish Otolith Research; Secor, D.H., Dean, J.M., Campana, S.E., Eds.; University of South Carolina Press: Columbia, SC, USA, 1995; pp. 545–565. [Google Scholar]
  17. Porta, M.J.; Snow, R.A. Validation of Annulus Formation in White Perch Otoliths, Including Characteristics of an Invasive Population. J. Freshw. Ecol. 2017, 32, 489–498. [Google Scholar] [CrossRef] [Green Version]
  18. Froeschke, B.; Allen, L.G.; Pondella, D.J. Life History and Courtship Behavior of Black Perch, Embiotoca jacksoni (Teleostomi: Embiotocidae), from Southern California. Pac. Sci. 2007, 61, 521–531. [Google Scholar] [CrossRef] [Green Version]
  19. Okamura, H.; Punt, A.E.; Semba, Y.; Ichinokawa, M. Marginal Increment Analysis: A New Statistical Approach of Testing for Temporal Periodicity in Fish Age Verification. J. Fish Biol. 2013, 82, 1239–1249. [Google Scholar] [CrossRef]
  20. Smith, J. Age Validation of Lemon Sole (Microstomus kitt), Using Marginal Increment Analysis. Fish. Res. 2014, 157, 41–46. [Google Scholar] [CrossRef]
  21. Hidalgo-de-la-Toba, J.A.; Morales-Bojórquez, E.; González-Peláez, S.S.; Bautista-Romero, J.J.; Lluch-Cota, D.B. Modeling the Temporal Periodicity of Growth Increments Based on Harmonic Functions. PLoS ONE 2018, 13, 0196189. [Google Scholar] [CrossRef] [Green Version]
  22. Prince, E.D.; Lee, D.W.; Berkeley, S.A. Use of Marginal Increment Analysis to Validate the Anal Spine Method for Ageing Atlantic Swordfish and Other Alternatives for Age Determination. Col. Vol. Sci. Pap. ICCAT 1988, 27, 194–201. [Google Scholar]
  23. Pearson, D.E. Timing of Hyaline-Zone Formation as Related to Sex, Location, and Year of Capture in Otoliths of the Widow Rockfish, Sebastes entomelas. Fish. Bull. 1996, 94, 190–197. [Google Scholar]
  24. Laidig, T.E.; Pearson, D.E.; Sinclair, L.L. Age and Growth of Blue Rockfish (Sebastes mystinus) from Central and Northern California. Fish Bull 2003, 101, 800–808. [Google Scholar]
  25. Young, J.; Drake, A.; Groisson, A.-L. Age and Growth of Broadbill Swordfish (Xiphias glades) from Eastern Australian Waters—Preliminary Results; CSIRO, Division of Marine Research: Hobart, Australia, 2003. [Google Scholar]
  26. Chang, S.-K.; Chou, Y.-T.; Hoyle, S.D. Length-Weight Relationships and Otolith-Based Growth Curves for Brushtooth Lizardfish off Taiwan with Observations of Region and Aging-Material Effects on Global Growth Estimates. Front. Mar. Sci. 2022, 9, 921594. [Google Scholar] [CrossRef]
  27. Chang, S.-K.; Yuan, T.-L.; Hoyle, S.D.; Farley, J.H.; Shiao, J.-C. Growth Parameters and Spawning Season Estimation of Four Important Flyingfishes in the Kuroshio Current off Taiwan and Implications from Comparisons with Global Studies. Front. Mar. Sci. 2022, 8, 747382. [Google Scholar] [CrossRef]
  28. Yamaguchi, A.; Kume, G.; Takita, T. Geographic Variation in the Growth of White Croaker, Pennahia argentata, off the Coast of Northwest Kyushu, Japan. Environ. Biol. Fishes 2004, 71, 179–188. [Google Scholar] [CrossRef]
  29. Yan, Y.R.; Hou, G.; Lu, H.S.; Li, Z.L. Growth Characteristics and Population Composition of Big-Head Pennah Croaker, Pennahia macrocephalus in the Beibu Gulf. J. Fish. Sci. China 2010, 40, 61–68. [Google Scholar]
  30. Ju, P.L.; Yang, L.; Lu, Z.B.; Yang, S.Y.; Du, J.G.; Zhong, H.Q.; Chen, J.; Xiao, J.M.; Chen, M.R.; Zhang, C.Y. Age, Growth, Mortality and Population Structure of Silver Croaker Pennahia argentata (Houttuyn, 1782) and Red Bigeye Priacanthus macracanthus Cuvier, 1829 in the North-Central Taiwan Strait. J. Appl. Ichthyol. 2016, 32, 652–660. [Google Scholar] [CrossRef]
  31. Jeon, B.S.; Choi, J.H.; Kim, D.N.; Im, Y.J.; Lee, H.W. Age and Growth of White Croaker Pennahia argentata in the Southern Sea of Korea by Otolith Analysis. Korean J. Fish. Aquat. Sci. 2021, 54, 53–63. [Google Scholar]
  32. Attaqi, A.N. Age Structure and Growth of Big-Head Pannah Croaker, Pennahia macrocephalus in the Southwestern Waters off Taiwan: An Approach Using Both Thin-Section and Weight of Otoliths. Master’s Thesis, College of Ocean Science and Resource, National Taiwan Ocean University,, Keelung, Taiwan, 2018. [Google Scholar]
  33. Vanderkooy, S.; Carroll, J.; Elzey, S.; Gilmore, J.; Kipp, J. (Eds.) A Practical Handbook for Determining the Ages of Gulf of Mexico and Atlantic Coast Fishes, 3rd ed.; Gulf States Marine Fisheries Commission and Atlantic States Marine Fisheries Commission: Ocean Springs, MS, USA, 2020. [Google Scholar]
  34. Tsai, T.W. Age and Growth of Big Head Pennah Croaker (Pennahia macrocephalus) Sampled from the Southwestern Waters off Taiwan. Master’s Thesis, Department of Fisheries Production and Management of National Kaohsiung University of Science and Technology, Kaohsiung, Taiwan, 2009. [Google Scholar]
  35. Farley, J.; Krusic-Golub, K.; Eveson, P.; Clear, N.; Roupsard, F.; Sanchez, C.; Nicol, S.; Hampton, J. Age and Growth of Yellowfin and Bigeye Tuna in the Western and Central Pacific Ocean from Otoliths (No. WCPFC-SC16-2020/SA-WP-02). In Proceedings of the WCPFC Scientific Committee 16th Regular Session, Electronic Meeting, 11–20 August 2020. [Google Scholar]
  36. Roff, D.A. A Motion for the Retirement of the Von Bertalanffy Function. Can. J. Fish. Aquat. Sci. 1980, 37, 127–129. [Google Scholar] [CrossRef]
  37. Katsanevakis, S. Modelling Fish Growth: Model Selection, Multi-Model Inference and Model Selection Uncertainty. Fish. Res. 2006, 81, 229–235. [Google Scholar] [CrossRef]
  38. Katsanevakis, S.; Maravelias, C.D. Modelling Fish Growth Multi-model Inference as a Better Alternative to a Priori Using von Bertalanffy Equation. Fish Fish. 2008, 9, 178–187. [Google Scholar] [CrossRef]
  39. Williams, A.J.; Farley, J.H.; Hoyle, S.D.; Davies, C.R.; Nicol, S.J. Spatial and Sex-Specific Variation in Growth of Albacore Tuna (Thunnus alalunga) across the South Pacific Ocean. PLoS ONE 2012, 7, e39318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Carbonara, P.; Intini, S.; Kolitari, J.; Joksimović, A.; Milone, N.; Lembo, G.; Casciaro, L.; Bitetto, I.; Zupa, W.; Spedicato, M.T.; et al. A Holistic Approach to the Age Validation of Mullus barbatus L., 1758 in the Southern Adriatic Sea (Central Mediterranean). Sci. Rep. 2018, 8, 13219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Velasco, G.; Oddone, M.C. Growth Parameters and Growth Performance Indexes for Some Populations of Marine Catfishes (Actinopterygii, Siluriformes, Ariidae). Acta Biol. Leopoldensia 2004, 26, 307–313. [Google Scholar]
  42. Ba, K.; Thiaw, M.; Lazar, N.; Sarr, A.; Brochier, T.; Ndiaye, I.; Faye, A.; Sadio, O.; Panfili, J.; Thiaw, O.T.; et al. Resilience of Key Biological Parameters of the Senegalese Flat Sardinella to Overfishing and Climate Change. PLoS ONE 2016, 11, e0156143. [Google Scholar] [CrossRef] [Green Version]
  43. Murua, H.; Rodríguez-Marín, E.; Neilson, J.; Farley, J.; Juan-Jordá, M.J. Fast versus Slow Growing Tuna Species: Age, Growth, and Implications for Population Dynamics and Fisheries Management. Rev. Fish Biol. Fish. 2017, 27, 733–773. [Google Scholar] [CrossRef]
  44. Pauly, D.; Munro, J.L. Once More on the Comparison of Growth in Fish and Invertebrates. Fishbyte 1984, 2, 21. [Google Scholar]
  45. Juan-Jordá, M.J.; Mosqueira, I.; Dulvy, N. The Conservation and Management of Tunas and Their Relatives: Setting Life History Research Priorities. PLoS ONE 2013, 8, e70405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Pauly, D.; Binohlan, C. FishBase and AUXIM as Tools for Comparing the Life-History Patterns, Growth and Natural Mortality of Fish: Applications to Snappers and Groupers. ICLARM Conf. Proc. 1996, 48, 223–247. [Google Scholar]
  47. Chang, S.-K.; DiNardo, G.; Farley, J.; Brodziak, J.; Yuan, Z.-L. Possible Stock Structure of Dolphinfish (Coryphaena hippurus) in Taiwan Coastal Waters and Globally Based on Reviews of Growth Parameters. Fish. Res. 2013, 147, 127–136. [Google Scholar] [CrossRef]
  48. Chang, S.-K.; Maunder, M. Aging Material Matters in the Estimation of von Bertalanffy Growth Parameters for Dolphinfish (Coryphaena hippurus). Fish. Res. 2012, 119–120, 147–153. [Google Scholar] [CrossRef]
  49. Nelson, J.S. Fishes of the World, 4th ed.; John Wiley & Sons Inc: Hoboken, NJ, USA, 2006. [Google Scholar]
  50. Shao, K.T. Taiwan Fish Database. WWW Web Electronic Publication. 2022. Available online: http://fishdb.sinica.edu.tw (accessed on 21 February 2022).
  51. Tang, D.-S. A study of sciaenoid fishes of China. Amoy Mar. Biol. Bull. 1937, 2, 47–88. [Google Scholar]
  52. Jordan, D.S.; Thompson, W.F. A Review of the Sciaenoid Fishes of Japan. Proc. United States Natl. Mus. 1911, 39, 241–261. [Google Scholar] [CrossRef]
  53. Chang, S.-K.; Yuan, T.L.; Wang, S.P.; Chang, Y.J.; DiNardo, G.; Chen, Y.; Mello, L. Deriving a Statistically Reliable Abundance Index from Landings Data: An Application to the Taiwanese Coastal Dolphinfish Fishery with a Multispecies Feature. Trans. Am. Fish. Soc. 2019, 148, 106–122. [Google Scholar] [CrossRef] [Green Version]
  54. Froese, R. Cube Law, Condition Factor and Weight–Length Relationships: History, Meta-Analysis and Recommendations. J. Appl. Ichthyol. 2006, 22, 241–253. [Google Scholar] [CrossRef]
  55. Zakeyudin, M.-S.; Mat, M.; Salmah, M.; Md Sah, A.S.R.; Hassan, A. Assessment of Suitability of Kerian River Tributaries Using Length-Weight Relationship and Relative Condition Factor of Six Freshwater Fish Species. J. Environ. Earth Sci. 2012, 2, 52–60. [Google Scholar]
  56. Jisr, N.; Younes, G.; Sukhn, C.; El-Dakdouki, M.H. Length-Weight Relationships and Relative Condition Factor of Fish Inhabiting the Marine Area of the Eastern Mediterranean City, Tripoli-Lebanon. Egypt. J. Aquat. Res. 2018, 44, 299–305. [Google Scholar] [CrossRef]
  57. James, G.; Witten, D.; Hastie, T.; Tibshirani, R. (Eds.) An Introduction to Statistical Learning: With Applications in R; Springer Texts in Statistics; Springer: New York, NY, USA, 2013; ISBN 978-1-4614-7137-0. [Google Scholar]
  58. Ashworth, E.C.; Hesp, S.A.; Hall, N.G. A New Proportionality-Based Back-Calculation Approach, Which Employs Traditional Forms of Growth Equations, Improves Estimates of Length at Age. Can. J. Fish. Aquat. Sci. 2017, 74, 1088–1099. [Google Scholar] [CrossRef] [Green Version]
  59. Ogle, D.H. Introductory Fisheries Analyses with R, 1st ed.; CRC Press: New York, NY, USA, 2016; ISBN 978-1-315-37198-6. [Google Scholar]
  60. Karlson, S.; Michalsen, K.; Folkvord, A.; Karlson, S.; Age, A. Age Determination of Atlantic Halibut (Hippoglossus hippoglossus L.) along the Coast of Norway: Status and Improvements. ICES J. Mar. Sci. 2012, 70, 50–55. [Google Scholar] [CrossRef]
  61. Gompertz, B. On the Nature of the Function Expressive of the Law of Human Mortality, and on a New Mode of Determining the Value of Life Contingencies. Philos. Trans. R. Soc. 1825, 115, 513–583. [Google Scholar] [CrossRef]
  62. Von Bertalanffy, L. A Quantitative Theory of Organic Growth. Hum. Biol. 1938, 10, 181–213. [Google Scholar]
  63. Richards, F.J. A Flexible Growth Function for Empirical Use. J. Exp. Bot. 1959, 10, 290–301. [Google Scholar] [CrossRef]
  64. Ricker, W.E. Computation and Interpretation of Biological Statistics of Fish Populations. Bull. Fish. Res. Board Can. 1975, 191, 1–382. [Google Scholar]
  65. Flinn, S.A.; Midway, S.R. Trends in Growth Modeling in Fisheries Science. Fishes 2021, 6, 1. [Google Scholar] [CrossRef]
  66. Quist, M.C.; Pegg, M.A.; Devries, D.R. Age and Growth. Fisheries Techniques. In Age and Growth. Fisheries Techniques; Zale, A.V., Parrish, D.L., Sutton, T.M., Eds.; American Fisheries Society: Bethesda, MD, USA, 2012; pp. 677–731. [Google Scholar]
  67. R-Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  68. Akaike, H. Information Theory and an Extension of the Maximum Likelihood Principle. In Proceedings of the 2nd International Symposium on Information Theory, Tsahkadsor, Armenia, 2–8 September 1973; Volume I. pp. 267–281. [Google Scholar]
  69. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, 2nd ed.; Springer: New York, NY, USA, 2002; ISBN 978-0-387-95364-9. [Google Scholar]
  70. Wang, X.; Qiu, Y.; Du, F.; Lin, Z.; Sun, D.; Huang, S. Population Parameters and Dynamic Pool Models of Commercial Fishes in the Beibu Gulf, Northern South China Sea. Chin. J. Oceanol. Limnol. 2012, 30, 105–117. [Google Scholar] [CrossRef]
  71. Chen, Z.; Huang, Z.; Qiu, Y. 南海北部白姑鱼生长和死亡参数的估算 (Estimation of Growth and Mortality Parameters of Argyrosomus argentatus in Northern South China Sea). Chin. J. Appl. Ecol. 2005, 16, 712–716. (In Chinese). Available online: http://www.cqvip.com/qk/90626a/200504/15404417.html (accessed on 1 October 2021).
  72. FishBase. World Wide Web Electronic Publication. 2022. Available online: https://www.fishbase.de/summary/Pennahia-anea (accessed on 1 October 2021).
  73. Nhuận, M.C.; Khương, Đ.V. Stock Biological Characteristics of Big-Head Pennah Croacer (Pennahia macrocephalus, Tang 1937) in the Shared Fishing Zone in the Gulf of Tonkin, during 2006–2010. J. Fish. Sci. Technol. Vietnam. 2013, 2, 83–88. [Google Scholar]
  74. Liu, H.C.; Tzeng, W.N. Age and growth of the white croaker, Argyrosomus argentatus (HOUTTUYN), in the Southern part of the East China Sea and Taiwan Strait. J. Fish. Soc. Taiwan 1972, 1, 21–38. [Google Scholar]
  75. Lu, Z.; Dai, Q.; Yan, Y. Population dynamics of major demersal fish species offshore Fujian. J. Oceanogr. Taiwan Strait Chin. 1999, 18, 100–105. [Google Scholar]
  76. Hu, Y.; Qian, S. A study on the age and growth of the white croaker. Mar. Fish. 1989, 4, 158–162. [Google Scholar]
  77. Higuchi, T.; Yamaguchi, A.; Takita, T. Age and Growth of White Croaker, Pennahia argentata, in the Ariake Sound, Japan. Bull. Fac. Fish. Nagasaki Univ. Jpn. 2003, 47–51. [Google Scholar]
  78. Kakuda, S.; Matsumoto, K. On the Age and Growth of the White Croaker Argyrosomus argentatus. Hiroshima Univ. 1977, 16, 115–122. [Google Scholar]
  79. Yan, Y.; Hou, G.; Lu, H.; Yin, Q. Age and Growth of Pennahia pawak in the Beibu Gulf. Chin. Fish. Sci. 2011, 145–155. [Google Scholar]
  80. Yi, X.; Qiu, K.; Zhou, X.; Zhao, C.; Deng, Y.; He, X.; Yan, Y. Analysis of Fishery Biology of Pennahia pawak in Beibu Gulf. J. Shanghai Ocean Univ. 2021, 30, 515–524. [Google Scholar] [CrossRef]
  81. Bhuyan, S.K.; Kumar, P.; Jena, J.K.; Pillai, B.R.; Chakraborty, S.K. Studies on the Growth of Otolithes Ruber (Bloch & Schneider, 1801), Johnius Carutta Bloch, 1793 and Pennahia macrophthalmus (Bleeker, 1850) from Paradeep Coast, Orissa, India. Indian J. Fish 2012, 59, 89–93. [Google Scholar]
  82. Ingles, J.; Pauly, D. An Atlas of the Growth, Mortality and Recruitment of Philippine Fishes. Available online: https://www.worldfishcenter.org/publication/atlas-growth-mortality-and-recruitment-philippine-fishes (accessed on 2 May 2022).
  83. Ziegler, B. Growth and Mortality Rates of Some Fishes of Manila Bay, Philippines as Estimated from the Analysis of Length Frequencies. Master’s Thesis, Kiel University, Kiel, Germany, 1979. [Google Scholar]
  84. Abu Talib, A. Population Dynamics of Big-Eye Croaker (Pennahia macrophthalmus, Sciaenidae) off Kedah, Penang and Perak States, Malaysia; FAO Fisheries Report; FAO: Rome, Italy, 1988. [Google Scholar]
  85. Chakraborty, S.K.; Deshmukh, V.D.; Khan, M.Z.; Vidyasagar, K.; Raje, S.G. Estimates of Growth, Mortality, Recruitment Pattern and MSY of Important Resources from the Maharashtra Coast. J. Indian Fish. Assoc. 1994, 24, 1–39. [Google Scholar]
  86. Chakraborty, S.K. Stock assessment of big-eye croaker, Pennahia macropthalmus (Bleeker) (Pisces/Perciformes/Sciaenidae) from Bombay waters. Indian J. Mar. Sci. 1996, 25, 316–319. [Google Scholar]
  87. Jayasankar, P. Population Dynamics of Big-Eye Croaker Pennahia macrophthalmus and Blotched Croaker Nibea maculata (Pisces/Perciformes/Sciaenidae) in the Trawling Grounds off Rameswaram Island, East Coast of India. Indian J. Mar. Sci. 1995, 24, 153–157. [Google Scholar]
  88. Chakraborty, S.K.; Devadoss, P.; Manojkumar, P.P.; Feroz Khan, M.; Jayasankar, P.; Sivakami, S.; Gandi, V.; Appanna Sastry, Y.; Raju, A.; Livingston, P.; et al. Fishery, Biology and Stock Assessment of Jew Fish Resources of India. In Marine Fisheries Research and Management; Pillai, V.N., Menon, N.G., Eds.; Central Marine Fisheries Research Institute: Cochin, India, 2000; pp. 604–616. [Google Scholar]
  89. Chakraborty, S.K. Growth Studies of Sciaenids from Mumbai Waters Using Bhattarcharya’s Method. Naga ICLARM Q. 2001, 24, 40–41. [Google Scholar]
  90. Menon, M.; Maheswarudu, G.; Rohit, P.; Laxmilatha, P.; Madhumita, D.; Rao, K.N. Biology and Stock Assessment of the Bigeye Croaker Pennahia Anea (Bloch, 1793) Landed along Andhra Pradesh, North-East Coast of India. Indian J. Fish 2015, 62, 46–51. [Google Scholar]
  91. Wagiyo, K.; Tirtadanu, T.; Chodriyah, U.; Wang, X.; Qiu, Y.; Du, F.; Lin, Z.; Sun, D.; Huang, S. Biology Characteristic, Abundance Index and Fishing Aspect of Donkey Croaker Pennahia anea (Bloch, 1793) in the Tangerang Waters. E3S Web Conf. 2020, 30, 105–117. [Google Scholar] [CrossRef] [Green Version]
  92. Pauly, D. A Preliminary Compilation of Fish Length Growth Parameters. Ber. Inst. Meereskd. Univ. Kiel 1978, 55, 200. [Google Scholar]
  93. Kao, P. Stock Discrimination of Blackmouth Croaker Atrobucca Nibe in Taiwan. Master’s Thesis, The Institute of Fisheries Science, National Taiwan University, Taipei, Taiwan, 2019. (In Chinese). [Google Scholar]
  94. Sato, T. Fishery Biology of Black Croaker, Argyrosomus nibe (Jordan et Thompson). I. On the Age and Growth of the Black Croaker in the Central and Southern Parts of the East China Sea. Bull. Seikai Reg. Fish Res. Lab. 1963, 29, 75–96. [Google Scholar]
  95. Salarpouri, A.; Kaymaram, F.; Valinassab, T.; Behzadi, S.; Darvishi, M.; Kamali, E.; Rezwani, S.; Memarzadeh, M.; Karami, N. A Survey on Black Mouth Croaker (Atrobucca nibe) Resources in the North-West of Oman Sea; Iranian Fisheries Science Research Institute: Tehran, Iran, 2015. [Google Scholar]
  96. Chang, J.J. Effects of Fishing on Life History Parameters of Atrobucca Nibe from Surrounding Waters of Guei-Shan Island, Northeastern Taiwan. Master’s Thesis, College of Ocean Science and Resource, National Taiwan Ocean University, Keelung, Taiwan, 2008. [Google Scholar]
  97. Tsai, C.N. Age and growth of black croaker, Atrobucca nibe, in adjacent waters of Tungkang, Taiwan. Master’s Thesis, National Taiwan Ocean Univeristy, Keelung, Taiwan, 1993. [Google Scholar]
  98. Hwang, G.M.; Chen, C.T. Fishery Biology of Black Croaker, Atrobucca Nibe (Jordan et Thompson) from the Surrounding Waters of Guei Shan Island, Taiwan. J. Fish. Soc. Taiwan 1984, 11, 35–52. [Google Scholar]
  99. Memon, K.; Liu, Q.; Kalhoro, M.; Nazir, K.; Waryani, B.; Chang, M.S.; Mirjatt, A.N. Estimation of Growth and Mortality Parameters of Croaker Atrobucca alcocki in Pakistani Waters. J. Agric. Sci. Technol. 2016, 18, 669–679. [Google Scholar]
  100. Zhu, L.; Li, L.; Liang, Z. Comparison of Six Statistical Approaches in the Selection of Appropriate Fish Growth Models. Chin. J. Oceanol. Limnol. 2009, 27, 457–467. [Google Scholar] [CrossRef]
  101. Kimmerer, W.; Avent, S.R.; Bollens, S.M.; Feyrer, F.; Grimaldo, L.F.; Moyle, P.B.; Nobriga, M.; Visintainer, T. Variability in Length–Weight Relationships Used to Estimate Biomass of Estuarine Fish from Survey Data. Trans. Am. Fish. Soc. 2005, 134, 481–495. [Google Scholar] [CrossRef]
  102. Morey, G.; Moranta, J.; Massutí, E.; Grau, A.; Linde, M.; Riera, F.; Morales-Nin, B. Weight–Length Relationships of Littoral to Lower Slope Fishes from the Western Mediterranean. Fish. Res. 2003, 62, 89–96. [Google Scholar] [CrossRef]
  103. CMFRI/BOBP-IGO/GoI. Training Manual on Stock Assessment of Tropical Fishes; Central Marine Fisheries Research Institute (CMFRI): Kochi, India, 2016; Volume 156. [Google Scholar]
  104. Mredul, M.M.H.; Alam, M.R.; Akkas, A.B.; Sharmin, S.; Pattadar, S.N.; Ali, M.L. Some Reproductive and Biometric Features of the Endangered Gangetic Leaf Fish, Nandus nandus (Hamilton, 1822): Implication to the Baor Fisheries Management in Bangladesh. Aquac. Fish. 2021, 6, 634–641. [Google Scholar] [CrossRef]
  105. Tesch, F.W. Age and growth. In Methods for Assessment of Fish Production in Freshwaters; Ricker, W.E., Ed.; Blackwell Scientific Publications: Oxford, UK, 1968; pp. 101–113. [Google Scholar]
  106. Lee, J.J. Reproductive Biology of Big Head Pennah Croaker (Pennahia macrocephalus) Sampled from the Southwestern Waters off Taiwan. Master’s Thesis, Department of Fisheries Production and Management of National Kaohsiung University of Science and Technology, Kaohsiung, Taiwan, 2010. (In Chinese). [Google Scholar]
  107. Matsui, I. 東海黄海に於ける底曳網漁場と底棲生物群聚との関係に就て (Relation between the Trawling Grounds and the Benthic Association in the East China Sea and the Yellow Sea). Nippon Suisan Gakkaishi Jpn. 1951, 16, 159–167. [Google Scholar] [CrossRef]
  108. Kao, P.; Rong, J.; Zhang, W.; Li, Y. A Review of the Research on the Biology of Atrobucca nibe Fishery. Fish. Ext. Natl. Taiwan Univ. 2018, 29, 46–57. [Google Scholar]
  109. Liu, F.; Chen, J. A survey of Nibea nibe Jordan et Thompson in northern Taiwan. China Aquat. Prod. 1954, 24, 23–45. [Google Scholar]
  110. Hsiao, L.; Chen, C.; Wu, C.; Chen, Y.; He, J. Research on Reproductive Biology of Heliconia in the Southwest Waters of Taiwan. J. Taiwan Fish. Res. 2017, 25, 15–25. [Google Scholar]
  111. Schneider, J.C.; Laarman, P.W.; Gowing, H. Length-Weight Relationships. Chapter 17. In Manual of Fisheries Survey Methods II: With Periodic Updates; Michigan Department of Natural Resources, Fisheries Division: Lansing, MI, USA, 2000; pp. 1–18. [Google Scholar]
  112. De Giosa, M.; Czerniejewski, P.; Rybczyk, A. Seasonal Changes in Condition Factor and Weight-Length Relationship of Invasive Carassius gibelio (Bloch, 1782) from Leszczynskie Lakeland, Poland. Adv. Zool. 2014, 2014, e678763. [Google Scholar] [CrossRef] [Green Version]
  113. He, J.S.; Huang, X.H.; Wu, Y.S.; Lai, C.C.; Weng, J.S. Reproductive Biology of the Females of the Pennahia macrocephalus in the Southwestern Waters of Taiwan. J. Taiwan Fish. Res. 2020, 28, 13–24. [Google Scholar]
  114. Costa, M.D.P.; Muelbert, J.H.; Moraes, L.E.; Vieira, J.; Castello, J.P. Estuarine Early Life Stage Habitat Occupancy Patterns of Whitemouth Croaker Micropogonias furnieri (Desmarest, 1830) from the Patos Lagoon, Brazil. Fish. Res. 2014, 160, 77–84. [Google Scholar] [CrossRef]
  115. Putnis, I.; Korņilovs, G. Manual for Age Determination of Baltic Herring. In Proceedings of the ICES 2008 Report of the Workshop on Age Reading of Baltic Herring (WKARBH), Riga, Latvia, 9–13 June 2008; p. 37. [Google Scholar]
  116. Newman, S.J.; Dunk, I.J. Age Validation, Growth, Mortality, and Additional Population Parameters of the Goldband Snapper (Pristipomoides multidens) off the Kimberley Coast of Northwestern Australia. Fish. Bull. 2003, 101, 116–128. [Google Scholar]
  117. Hara, K.; Furumitsu, K.; Aoshima, T.; Kanehara, H.; Yamaguchi, A. Age, Growth, and Age at Sexual Maturity of the Commercially Landed Skate Species, Dipturus chinensis (Basilewsky, 1855), in the Northern East China Sea. J. Appl. Ichthyol. 2018, 34, 66–72. [Google Scholar] [CrossRef]
  118. Wang, S.B. Survey of Fisheries Economic Activities in the Nearby Area of Mai-Liao, Yunlin County, Taiwan; The Final Report for Capture Fisheries in 2012; Century Engineering Consultant Co., Ltd.: New Taipei, Taiwan, 2011. (In Chinese) [Google Scholar]
  119. Jones, C. Determining Age of Larval Fish with the Otolith Increment Technique. Fish. Bull. 1986, 84, 91–102. [Google Scholar]
  120. Cailliet, G.M.; Goldman, K.J.; Carrier, J.; Musick, J.A.; Heithaus, M. (Eds.) Age Determination and Validation in Chondrichthyan Fishes; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  121. Chen, Y.; Mello, L.G.S. Growth and Maturation of Cod (Gadus morhua) of Different Year-Classes in the Northwest Atlantic, NAFO Subdivision 3Ps. Fish. Res. 1999, 42, 87–101. [Google Scholar] [CrossRef]
  122. Fisheries Agency. Fisheries Statistical Yearbook—Taiwan, Kinmen and Matsu Area; Fisheries Agency, Council of Agriculture, Executive Yuan: Taipei, Taiwan, 2021. (In Chinese) [Google Scholar]
  123. Haimovici, M.; Cavole, L.M.; Cope, J.M.; Cardoso, L.G. Long-Term Changes in Population Dynamics and Life History Contribute to Explain the Resilience of a Stock of Micropogonias furnieri (Sciaenidae, Teleostei) in the SW Atlantic. Fish. Res. 2021, 237, 105878. [Google Scholar] [CrossRef]
  124. Rice, J.; Gislason, H. Patterns of Change in the Size Spectra of Numbers and Diversity of the North Sea Fish Assemblage, as Reflected in Surveys and Models. ICES J. Mar. Sci. 1996, 53, 1214–1225. [Google Scholar] [CrossRef]
Figure 1. Two important croaker species occurred in Taiwan. Pennahia macrocephalus (big-head pennah croaker, age = 3.5) and Atrobucca nibe (blackmouth croaker, age = 2).
Figure 1. Two important croaker species occurred in Taiwan. Pennahia macrocephalus (big-head pennah croaker, age = 3.5) and Atrobucca nibe (blackmouth croaker, age = 2).
Fishes 07 00281 g001
Figure 2. Black circles indicate the major fishing areas of P. macrocephalus and A. nibe from the trawlers based in Zihquan, Tungkang, Fangliao and other ports in southwestern Taiwan. The market samples used in this study came from these fishing areas.
Figure 2. Black circles indicate the major fishing areas of P. macrocephalus and A. nibe from the trawlers based in Zihquan, Tungkang, Fangliao and other ports in southwestern Taiwan. The market samples used in this study came from these fishing areas.
Fishes 07 00281 g002
Figure 3. Sample images of sagittal otolith (transverse sections) showing the determination of ages by different methods. (a) Integral counting method (sample of A. nibe; total length 41.5 cm): rn indicate the ring marks on the outer margin of the opaque bands used for determining age class n. (N: nucleus). For this case, the age is six. (b) Quartile of marginal growth band method (quantile method, sample of P. macrocephalus; total length 16.8 cm): r1 and r2 indicate the ring marks on the outer margin of the opaque bands. a is the width of the most recent formed annulus, and b is the width of the hyaline band extended from the margin of the outmost opaque band to the edge. Ages were determined as the number of opaque bands plus approximate quartile value of b against a. For this case, the age is 2.25.
Figure 3. Sample images of sagittal otolith (transverse sections) showing the determination of ages by different methods. (a) Integral counting method (sample of A. nibe; total length 41.5 cm): rn indicate the ring marks on the outer margin of the opaque bands used for determining age class n. (N: nucleus). For this case, the age is six. (b) Quartile of marginal growth band method (quantile method, sample of P. macrocephalus; total length 16.8 cm): r1 and r2 indicate the ring marks on the outer margin of the opaque bands. a is the width of the most recent formed annulus, and b is the width of the hyaline band extended from the margin of the outmost opaque band to the edge. Ages were determined as the number of opaque bands plus approximate quartile value of b against a. For this case, the age is 2.25.
Fishes 07 00281 g003aFishes 07 00281 g003b
Figure 4. Total length (TL) and body weight (BW) data and the estimated length-weight equation of P. macrocephalus and A. nibe. The blue lines and circles are male, and the red lines and circles are female.
Figure 4. Total length (TL) and body weight (BW) data and the estimated length-weight equation of P. macrocephalus and A. nibe. The blue lines and circles are male, and the red lines and circles are female.
Fishes 07 00281 g004
Figure 5. Monthly percent occurrence of the translucent band at the edge of samples from P. macrocephalus and A. nibe. Numbers indicate the sample sizes.
Figure 5. Monthly percent occurrence of the translucent band at the edge of samples from P. macrocephalus and A. nibe. Numbers indicate the sample sizes.
Fishes 07 00281 g005
Figure 6. Age and length of sex-combined and sex-separated data and fitted growth curves of the three growth models for P. macrocephalus off southwestern Taiwan. (a,d,e) are the quartile of the marginal growth band method; (b) is the integral counting method; (c) is the back-calculation method.
Figure 6. Age and length of sex-combined and sex-separated data and fitted growth curves of the three growth models for P. macrocephalus off southwestern Taiwan. (a,d,e) are the quartile of the marginal growth band method; (b) is the integral counting method; (c) is the back-calculation method.
Fishes 07 00281 g006
Figure 7. Age and length of sex-combined and sex-separated data and fitted growth curves of the three growth models for A. nibe off southwestern Taiwan.
Figure 7. Age and length of sex-combined and sex-separated data and fitted growth curves of the three growth models for A. nibe off southwestern Taiwan.
Fishes 07 00281 g007
Figure 8. Auximetric plot of the two variables, log(L) and log(k), of the growth performance index formula obtained from the previous studies globally of Sciaenidae family. Ellipses with a 95% confidence interval. Code referred to Table 1.
Figure 8. Auximetric plot of the two variables, log(L) and log(k), of the growth performance index formula obtained from the previous studies globally of Sciaenidae family. Ellipses with a 95% confidence interval. Code referred to Table 1.
Fishes 07 00281 g008
Table 1. Global estimates of growth parameters of Pennahia and Atrobucca spp. from the literature. Lengths are in the unit of cm and k in year−1. Refer to Table S1 for detailed information.
Table 1. Global estimates of growth parameters of Pennahia and Atrobucca spp. from the literature. Lengths are in the unit of cm and k in year−1. Refer to Table S1 for detailed information.
CodeSpeciesRegionMaterial (Method)L (TL)k (cm y−1)φRef
m_1P. macrocephalusBeibu Gulf, SCSOtolith (I)27.3710.4082.485[29]
m_2P. macrocephalusBeibu Gulf, SCSLen. freq.28.8990.5202.638[70]
m_3P. macrocephalusBeibu Gulf, SCSScale31.5020.5902.768[73]
m_4P. macrocephalusSouthwest, TWOtolith (D)27.0280.5962.639[34]
m_5P. macrocephalusYunlin, TWOtolith (B)21.3580.3712.228[32]
mc_1P. macrocephalusSouthwest, TWOtolith (I)19.5920.3282.099This study
mc_2P. macrocephalusSouthwest, TWOtolith (Q)18.0040.7892.408This study
mc_3P. macrocephalusSouthwest, TWOtolith (B)18.2580.6172.313This study
a_1P. argentataTW straitScale30.8870.9782.970[74]
a_2P. argentataTW straitScale35.7130.3752.680[74]
a_3P. argentataSouthern, ECSScale28.8781.2073.003[74]
a_4P. argentataSouthern, ECSScale34.7330.3152.580[74]
a_5P. argentataBeibu Gulf, SCSScale30.5000.3502.513[75]
a_6P. argentataBeibu Gulf, SCSScale28.2300.5002.600[76]
a_7P. argentataNorthern, SCSLen. freq.38.2000.4202.787[71]
a_8P. argentataBeibu Gulf, SCSLen. freq.31.5000.3502.541[71]
a_9P. argentataTW straitOtolith36.3010.3772.696[30]
a_10P. argentataSouthern Sea, SKOtolith35.5290.3802.681[31]
a_11P. argentataAriake Sound, JPOtolith31.2000.3842.573[77]
a_12P. argentataAriake Sound, JPOtolith29.4000.3602.493[77]
a_13P. argentata Tachibana Bay, JPOtolith35.4000.2722.533[28]
a_14P. argentata Tachibana Bay, JPOtolith29.0000.3422.459[28]
a_15P. argentataOmura Bay, JPOtolith37.0000.2562.545[28]
a_16P. argentataOmura Bay, JPOtolith37.7000.1812.410[28]
a_17P. argentataThe Sea of Goto, JPOtolith44.5000.2412.679[28]
a_18P. argentataThe Sea of Goto, JPOtolith45.2000.1952.600[28]
a_19P. argentataSeto Inland SeaScale39.560.3075.68[78]
p_1P. pawakBeibu Gulf, SCSScale (Logistic)22.0300.5802.449[79]
p_2P. pawakBeibu Gulf, SCSLen. freq.24.1500.3902.357[80]
p_3P. pawakBeibu Gulf, SCSLen. freq.22.0500.3202.192[80]
an_1P. aneaParadeep, INLen. freq.30.3000.8602.897[81]
an_2P. aneaSan Miguel Bay, PHILen. freq.20.0000.6002.380[82]
an_3P. aneaManila Bay, PHILen. freq.26.5001.4002.993[83]
an_4P. aneaPenang/Perak, MANA34.2000.4002.670[84]
an_5P. aneaBombay, INLen. freq.24.5000.6402.585[85]
an_6P. aneaBombay, INLen. freq.26.0001.2002.909[86]
an_7P. aneaRameswaram, INLen. freq.23.3001.2602.835[87]
an_8P. aneaMandapam, INLen. freq.26.0000.9802.821[88]
an_9P. aneaBombay, INLen. freq.27.3001.9403.160[89]
an_xP. aneaAndhra Pradesh, INLen. freq.33.0000.7002.882[90]
an_yP. aneaIndonesiaLen. freq.23.8900.8402.681[91]
n_1A. nibeFormosa StraitNA58.5000.1162.599[92]
n_2A. nibeNorth, TWLen. freq.58.7000.1452.699[92]
n_3A. nibeNorth, TWLen. freq.37.1000.3542.688[92]
n_4A. nibeNorth, TWLen. freq.37.1000.5232.857[92]
n_5A. nibeNorth, TWLen. freq.57.3000.1772.764[92]
n_6A. nibeSouth, SKLen. freq.50.9000.2272.769[92]
n_7A. nibeSouth, SKLen. freq.46.9000.2972.815[92]
n_8A. nibeTWScale56.2900.1202.580[93]
n_9A. nibeECSScale43.2000.2522.672[94]
n_10A. nibeOman SeaOtolith50.0000.2002.699[95]
n_11A. nibeGuei-Shan Island, TWOtolith48.0600.2752.803[96]
n_12A. nibeGuei-Shan Island, TWOtolith64.9100.1472.792[96]
n_13A. nibeNALen. freq.45.0000.2252.659[92]
n_14A. nibeNALen. freq.40.1000.2382.583[92]
n_15A. nibe Southwest, TWOtolith53.3840.1832.717[97]
n_16A. nibe Southwest, TWOtolith52.5770.1872.713[97]
n_17A. nibe Northeast, TWOtolith49.5270.2272.746[98]
n_18A. nibeNortheast, TWOtolith48.4170.2012.673[98]
nc_1A. nibeSouthwest, TWOtolith (VBGM)68.1490.130−0.919This study
nc_2A. nibeSouthwest, TWOtolith (Logistic)53.1110.3942.958This study
alA. alcockiPakistani watersLen. freq.47.2500.1802.604[99]
Note: L∞ (TL), the L∞ value in total length (TL, might be converted from the original value). For Sex: M, male; F, female; C, sex-combined. For Region: SCS, the South China Sea; TW, Taiwan; ECS, the East China Sea; SK, South Korea; PHI, Philippines; MA, Malaysia; IN, India; JP, Japan. For method: I, integral counting method; Q, quartile method; B, back-calculated method; D, daily growth ring method. NA, information not available.
Table 2. Parameter estimates (±standard error) from four candidate growth models, by sex, for P. macrocephalus and A. nibe off Southwestern Taiwan. The best-fit models are shown in bold.
Table 2. Parameter estimates (±standard error) from four candidate growth models, by sex, for P. macrocephalus and A. nibe off Southwestern Taiwan. The best-fit models are shown in bold.
SpeciesModelSexL (cm TL)k (Year−1)t0 (Year)pAICcAICcw
P. macrocephalusVBGMC18.004 (0.315)0.789 (0.157)−0.872 (0.368)-1029.0400.468
P. macrocephalusVBGMF18.673 (0.982)0.495 (0.239)−2.095 (1.280)-445.0800.386
P. macrocephalusVBGMM17.864 (0.352)0.892 (0.202)−0.571 (0.367)-585.8400.390
P. macrocephalusGompertzC17.921 (0.285)0.879 (0.162)−0.564 (0.294)-1029.810.7740.318
P. macrocephalusGompertzF18.647 (0.955)0.532 (0.243)−1.693 (1.015)-445.390.3110.330
P. macrocephalusGompertzM17.777 (0.316)1.002 (0.210)−0.286 (0.290)-586.170.3280.331
P. macrocephalusLogisticC17.850 (0.262)0.970 (0.168)−0.311 (0.241)-1030.591.5560.215
P. macrocephalusLogisticF18.643 (0.947)0.565 (0.246)−1.354 (0.808)-445.690.6120.284
P. macrocephalusLogisticM17.704 (0.288)1.112 (0.219)−0.055 (0.236)-586.520.6780.278
A. nibeVBGMC68.149 (8.182)0.130 (0.033)−0.919 (0.413)-2373.406.3800.029
A. nibeVBGMF60.548 (5.835)0.186 (0.047)−0.266 (0.439)-1443.834.0900.075
A. nibeVBGMM90.094 (37.291)0.070 (0.046)−1.926 (0.870)-911.611.7360.180
A. nibeGompertzC57.162 (3.512)0.263 (0.036)1.919 (0.194)-2369.532.5140.202
A. nibeGompertzF54.732 (3.226)0.318 (0.052)1.853 (0.160)-1441.391.6480.255
A. nibeGompertzM62.594 (9.428)0.197 (0.050)2.313 (0.720)-910.560.6790.306
A. nibeLogisticC53.111 (2.294)0.394 (0.040)2.958 (0.219)-2367.0200.709
A. nibeLogisticF52.142 (2.306)0.448 (0.052)2.804 (0.160)-1439.7400. 582
A. nibeLogisticM55.313 (5.285)0.322 (0.054)3.372 (0.615)-909.8800.429
A. nibeRichardsC55.345 (2.904)0.282 (0.035)1.843 (0.162)−3.661 × 106 (3.477 × 106)2371.954.9270.060
A. nibeRichardsF54.349 (3.260)0.322 (0.053)1.843 (0.183)−1.067 × 105 (3.060 × 106)1443.533.7880.088
A. nibeRichardsM57.794 (6.277)0.224 (0.049)1.944 (0.409)−6.994 × 106 (6.283 × 106)913.123.2480.085
Note: For sex: M, male; F, female; C, sex-combined.
Table 3. The coefficients a and b of the length–weight relationship equation (LWR) of P. macrocephalus (standard length) and A. nibe (total length) of the current and previous studies.
Table 3. The coefficients a and b of the length–weight relationship equation (LWR) of P. macrocephalus (standard length) and A. nibe (total length) of the current and previous studies.
SpeciesRegionabReference
P. macrocephalusBeibu Gulf, SCS2.37 × 10−53.248[29]
P. macrocephalusBeibu Gulf, SCS2.16 × 10−23.032[70]
P. macrocephalusBeibu Gulf, SCS0.01 × 10−6~0.01 × 10−72.95~3.57[73]
P. macrocephalusYunlin, TW4.30 × 10−22.747[32]
P. macrocephalusSouthwest, TW ♀4.70 × 10−22.730[106]
P. macrocephalusSouthwest, TW ♂4.08 × 10−22.786[106]
P. macrocephalusSouthwest, TW9.90 × 10−33.337[34]
P. macrocephalusSouthwest, TW2.77 × 1022.937This study
A. nibeYellow Sea and ECS2.03 × 10−33.463[107]; cited in [108]
A. nibeECS2.84 × 10−33.220[94]; cited in [108]
A. nibeNorth, TW6.39 × 10−43.207[109]; cited in [108]
A. nibeNortheast, TW1.41 × 10−32.934[98]; cited in [108]
A. nibeNortheast, TW6.92 × 10−33.100[96]
A. nibeSouthwest, TW1.21 × 10−32.968[97]
A. nibeSouthwest, TW2.00 × 10−32.907[110]
A. nibeOman Sea ♀1.21 × 10−22.939[95]
A. nibeOman Sea ♂7.5 × 10−33.074[95]
A. nibeSouthwest, TW1.45 × 10−22.886This study
Note: For Region: SCS, the South China Sea; TW, Taiwan; ECS, the East China Sea. ♀: female, ♂ male.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, S.-C.; Chang, S.-K.; Lai, C.-C.; Yuan, T.-L.; Weng, J.-S.; He, J.-S. Length–Weight Relationships, Growth Models of Two Croakers (Pennahia macrocephalus and Atrobucca nibe) off Taiwan and Growth Performance Indices of Related Species. Fishes 2022, 7, 281. https://doi.org/10.3390/fishes7050281

AMA Style

Huang S-C, Chang S-K, Lai C-C, Yuan T-L, Weng J-S, He J-S. Length–Weight Relationships, Growth Models of Two Croakers (Pennahia macrocephalus and Atrobucca nibe) off Taiwan and Growth Performance Indices of Related Species. Fishes. 2022; 7(5):281. https://doi.org/10.3390/fishes7050281

Chicago/Turabian Style

Huang, Shu-Chiang, Shui-Kai Chang, Chi-Chang Lai, Tzu-Lun Yuan, Jinn-Shing Weng, and Jia-Sin He. 2022. "Length–Weight Relationships, Growth Models of Two Croakers (Pennahia macrocephalus and Atrobucca nibe) off Taiwan and Growth Performance Indices of Related Species" Fishes 7, no. 5: 281. https://doi.org/10.3390/fishes7050281

Article Metrics

Back to TopTop