Printer Friendly

Analysis of sex-specific spawning biomass per recruit of the sailfish (Istiophorus platypterus) in the waters off eastern Taiwan.

Abstract--Analyses of sex-specific yield per recruit and spawning stock biomass per recruit were conducted to evaluate the current status of the sailfish sailfish, common name for a marine game and food fish belonging to the family Istiophoridae and related to the swordfish and the marlin. It is named for its high, wide dorsal fin, colored deep blue with black spots.  (Istiophorus platypterus) fishery in the waters off eastern Taiwan. Natural mortality rates estimated from Pauly's empirical equation were 0.26/yr for females and 0.27/yr for males. The current fishing mortality rates were estimated as 0.24/yr and 0.43/yr for females and males, respectively, which are much lower than the estimated [F.sub.0.1] (0.62/yr and 0.79/yr for females and males, respectively) and [F.sub.SSB SSB Statistisk Sentralbyrå (Statistics Norway)
SSB Super Smash Bros (video game)
SSB Space Studies Board
SSB Single Side Band
SSB Single Stranded DNA-Binding Protein
SSB Salomon Smith Barney
40] (0.46/yr for females) which are commonly used as target reference points in fisheries management. The effects of the fishing mortality, natural mortality, and age at first capture on the estimates of biological reference points were evaluated by using the Monte Carlo simulation Monte Carlo Simulation

A problem solving technique used to approximate the probability of certain outcomes by running multiple trial runs, called simulations, using random variables.
. The results indicate that failure to consider the uncertainty in parameters such as natural mortality or age at first capture may lead to the improper estimation of biological reference points. This study indicates the possibility of current fishing mortality exceeding the target biological reference points may be negligible for sailfish in the waters off eastern Taiwan. However, in view of the recent rapid increase in fishing effort, it is evident that the stock status and development of the fishery need to be closely monitored.

**********

Sailfish (Istiophorus platypterus) is a circumtropically distributed species (Hoolihan, 2005). Sailfish is a member of the billfish billfish

Any of several long-jawed fishes, especially those in the family Istiophoridae, including marlins, spearfishes, and sailfishes. The name is also applied to the gar, needlefish, and sauries (family Scomberesocidae), as well as to the swordfish (family Xiphiidae).
 family, Istiophoridae, which also includes marlins and spearfishes, and is considered a bycatch species in commercial fisheries. Off the eastern coast of Taiwan, sailfish are economically important and seasonally abundant from April to October (abundance peaks from May to July). Sailfish are mainly caught by drift gill nets, although some are also caught by set nets, harpoons, and as incidental Contingent upon or pertaining to something that is more important; that which is necessary, appertaining to, or depending upon another known as the principal.

Under Workers' Compensation statutes, a risk is deemed incidental to employment when it is related to whatever a
 bycatch in inshore in·shore  
adv. & adj.
1. Close to a shore.

2. Toward or coming toward a shore.


inshore
Adjective

in or on the water, but close to the shore:
 long-line fisheries (Chiang, 2004). There are virtually no discards of sailfish in Taiwan. For the past decade, the annual landings of sailfish off Taiwan waters have fluctuated between 500 and 1000 metric tons, of which over 50% have come from waters off Taitung (eastern Taiwan).

Globally, large predatory fish species, including billfish, are declining at alarming rates because of excessive exploitation (Myers and Worm, 2003). Recent increases in the exploitation of billfish stocks by both commercial and recreational fisheries clearly point to the need for accurate assessments, if the goal is to develop sustainable billfish fisheries (Uozumi, 2003). However, because few fisheries target sailfish, assessments have not been conducted, resulting in few or no effective management measures. For sailfish in the Pacific Ocean and Indian Ocean, no assessments have been conducted, and stock status remains unknown. Historical sailfish catch data from the Pacific Ocean are scant, as are data on the length and age composition of the catch. This lack of data precludes the use of most stock assessment tools such as production models and age-structured models (Punt Punt (pnt), ancient land S of Egypt accessible by way of the Red Sea. Its exact location has not been identified, but it probably included the Somali coast. , 1997; Prager and Goodyear, 2001; Liu et al., 2006). However, recent biological studies on sailfish in the waters off eastern Taiwan (Chiang et al., 2004, 2006) have provided an opportunity to apply yield per recruit (Y/R) and spawning biomass per recruit (SSB/R) models (Govender, 1995; Griffiths, 1997) to estimate biological reference points including the fishing mortality rate corresponding to the point where the slope of the yield-per-recruit curve equals 10% of the slope at the origin ([F.sub.0.1]; Gulland and Boerema, 1973) and the fishing mortality rate corresponding to a specific percentage (x%) of the spawning biomass per recruit at the unfished level ([F.sub.SSBx]) for the eastern Taiwan sailfish stock. The status of the fishery for this stock could be examined by comparing the fishing mortality at the current level ([F.sub.CUR]) with the biological reference points.

A per-recruit analysis requires information on growth, mortalities, and selectivity selectivity /se·lec·tiv·i·ty/ (se-lek-tiv´i-te) in pharmacology, the degree to which a dose of a drug produces the desired effect in relation to adverse effects.

selectivity

1.
 of fishing gear. Ca analysis (Ricker, 1975) is the most common method employed for estimating total mortality when data on the age composition of catch are available. For a specified natural mortality, [F.sub.CUR] could be computed simply by subtracting the natural mortality from the total mortality. However, in most cases large uncertainty is associated with the estimation of natural mortality and other life history parameters, which can lead to large uncertainty in the estimation of [F.sub.CUR] and biological reference points.

The objective of this study was to evaluate the current status of the sailfish fishery in waters off eastern Taiwan by comparing the current fishing mortality rate (estimated from analyzing length composition data collected from the fishery) with the biological reference points derived from the per-recruit analyses (Butter-worth et al., 1989; Sun et al., 2002, 2005). In addition, a Monte Carlo simulation study was conducted for evaluating the influence of uncertainty associated with mortalities and the age at first catch ([t.sub.c]) on the estimation of biological reference points. This study provides an approach that can be used to assess the status of fisheries for which limited information does not allow us to conduct a full stock assessment.

Materials and methods

Length and age composition of the catch

Length composition data were obtained by measuring sailfish landed at the Shinkang fish market in eastern Taiwan (Fig. 1) during the period from July 1998 to July 2005. Specimens were randomly selected from the landings and measured for their lengths and weights. The sex of each specimen was identified from the appearance of its gonads. Samples of the first dorsal dorsal /dor·sal/ (dor´s'l)
1. pertaining to the back or to any dorsum.

2. denoting a position more toward the back surface than some other object of reference; a synonym of posterior
 fin were taken from 1166 of the sampled individuals for which lengths were measured and used to age the sailfish (Chiang et al. 2004). These subsampled fish were used to construct sex-specific age-length keys, which in turn were used to convert the length-frequency data into age-composition data.

Estimating mortality rates

For each sex, the dynamics of a simulated year class can be projected forward from one year to another by using the exponential 1. (mathematics) exponential - A function which raises some given constant (the "base") to the power of its argument. I.e.

f x = b^x

If no base is specified, e, the base of natural logarthims, is assumed.
2.
 survival equation (Ricker, 1975):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII ASCII or American Standard Code for Information Interchange, a set of codes used to represent letters, numbers, a few symbols, and control characters. Originally designed for teletype operations, it has found wide application in computers. ] (1)

where [N.sub.t] = the number of fish at the beginning of age t;

M = the instantaneous natural mortality rate;

F = t the fishing mortality of fully-recruited fish; and

[S.sub.t] = the fishing gear selectivity of fish at age t.

[FIGURE 1 OMITTED]

Selectivity is the relative vulnerability of different age or size classes to the fishing gear. In this study, we assumed that the selectivity follows a dome-shaped distribution because our length-frequency data were mostly collected from gill nets. This dome-shaped selectivity can be quantified with the following normal distribution density function:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

where [mu] = the age at the mode of the dome-shaped selectivity; and

[sigma] = the standard deviation of the dome-shaped selectivity.

The expected catch ([[??].sub.t]) offish off·ish  
adj.
Inclined to be distant and reserved; aloof.



offish·ly adv.

off
 at age t can be estimated based on the catch equation (Ricker, 1975):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (3)

The parameters of [N.sub.0], F, [mu], and [sigma] can be estimated simultaneously by minimizing the following composite objective function:

[summation summation n. the final argument of an attorney at the close of a trial in which he/she attempts to convince the judge and/or jury of the virtues of the client's case. (See: closing argument)  over t] [([C.sub.t] - [[??].sub.t]).sup.2] + [(max([S.sub.t]) - 1).sup.2], (4)

where [C.sub.t] = the observed numbers of catch at age t.

The F estimated above was considered as the current fishing mortality ([F.sub.CUR]) in this study. The approach outlined above for the catch-curve analysis is similar to the one described in Rudershausen et al. (2008).

A total of 1000 independent bootstrap See boot.

(operating system, compiler) bootstrap - To load and initialise the operating system on a computer. Normally abbreviated to "boot". From the curious expression "to pull oneself up by one's bootstraps", one of the legendary feats of Baron von Munchhausen.
 samples of F were derived from 1000 sets of length-frequency data drawn randomly with replacement from the individuals of original length-frequency data.

Pauly's (1980) empirical equation was used to estimate M for each sex, and the mean sea surface temperature around eastern Taiwan waters fitted to the equation was about 26[degrees]C.

Per-recruit analyses

Yield per recruit (Y/R) of sailfish in the waters off eastern Taiwan was estimated from the following model:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

where [bar.W.sub.t] = the mean weight of fish at age t; and [t.sub.c] = the age at first capture.

Mean weight at age was computed as a power function of midyear mid·year  
n.
1. The middle of the calendar or academic year.

2.
a. An examination given in the middle of a school year.

b. midyears A series of such examinations.
 lower jaw fork length ([L.sub.t]):

[W.sub.t] = a x [L.sup.b.sub.t], (6)

and midyear lower jaw fork length was estimated from the von Bertalanffy growth function

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (7)

where K = the growth parameter;

[L.sub.[infinity] = the average asymptotic length; and

[t.sub.[infinity]] = hypothetic age at length of 0 (Ricker, 1975; see Table 1).

The maximum lifespan ([t.sub.max]) of sailfish in the waters off eastern Taiwan was unknown but was estimated by using the empirical relationship In science, an empirical relationship is one based solely on observation rather than theory. An empirical relationship requires only confirmatory data irrespective of theoretical basis.  of Taylor (1958):

[t.sub.max] = [t.sub.0] + [2.996/K], (8)

The equation for spawning stock biomass per recruit (SSB/R) is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9)

where [W.sub.t] = the mean weight at age t that was calculated from the von Bertalanffy function and length-weight relationship for female sailfish; and

fr = the fraction of female sailfish that are mature.

In this case fr is represented by a logistic lo·gis·tic   also lo·gis·ti·cal
adj.
1. Of or relating to symbolic logic.

2. Of or relating to logistics.



[Medieval Latin logisticus, of calculation
 equation fitted to maturity data collected from sailfish caught in the eastern waters off Taiwan (Chiang et al., 2004, 2006). The logistic equation can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (10)

where [r.sub.m] = the slope of the logistic curve; and [t.sub.m] = age at which 50% of fish are mature.

Biological reference points

The following biological reference points were estimated in order to determine the current status of the sailfish fishery: [F.sub.0.1], [F.sub.SSB25], and [F.sub.SSB40] x [F.sub.SSB25], and [F.sub.SSB40] are fishing mortality rates corresponding to the 25% and 40% of the spawning biomass per recruit at unfished level. The choice of 25% or 40% was relatively arbitrary for the fishery, but these values have been used as different levels of reference points for other relatively long-lived marine fishes (e.g., Griffiths, 1997; Kirchner, 2001; Sun et al., 2002, 2005). The spawning potential ratio (SPR spr Spring
SPR Strategic Petroleum Reserve
SPR Surface Plasmon Resonance
SPR Suomen Punainen Risti
SpR Specialist Registrar (UK doctor who supports a consultant)
SPR Society for Psychical Research
SPR Stop Prisoner Rape
) is the SSB/R at a given fishing mortality divided by the SSB/R without fishing (Gabriel et al., 1989; Goodyear, 1993; Katsukawa et al., 1999; Watanabe et al., 2000; Sun et al., 2002, 2005) and can be calculated as

SPR = [SSB/R/[SSB/R|.sub.F=0]]. (11)

Several authors have advocated designating [F.sub.0.1] or [F.sub.SSB40] as target reference points and [F.sub.SSB25] as a threshold reference point in order to obtain near optimal yields while minimizing the likelihood of stock collapse (Gulland and Boerema, 1973; Deriso, 1987; Hilden, 1993; Sun et al., 2002, 2005). We adopted these target and threshold reference points in this study.

Simulation study

The Monte Carlo simulation approach was applied to evaluate the sensitivity of estimating biological reference points with respect to parameters F, M, and [t.sub.c]. To quantify the uncertainty of F, F was assumed to follow a normal distribution with a mean and standard deviation ([[sigma].sub.F])--the latter estimated from the bootstrapped estimation of F. However, there was no information on the distribution of M and [t.sub.c]. We assumed a uniform distribution for M and [t.sub.c] by referring to the estimation of Pauly's empirical equation and the age at full recruitment from the age composition of sailfish in the waters off eastern Taiwan. The values of M and [t.sub.c] were sampled randomly from the corresponding uniform distributions defined in Table 2. Ten scenarios were designed to examine the effects of different combinations of the uncertain in parameters F, M, and [t.sub.c] on the estimation of biological reference points (Table 2). Scenarios A-D were used to evaluate the effects of changes in a single parameter, scenarios E-H were used to evaluate the results of changes in combinations of two parameters, and scenarios I and J were used to evaluate the results when three parameters were subject to uncertainty. For each scenario, 100 replicates of biological reference points were estimated by using the parameters of F, M, or [t.sub.c] randomly drawn from their assumed distributions. The median and the interquartile range In descriptive statistics, the interquartile range (IQR), also called the midspread, middle fifty and middle of the #s, is a measure of statistical dispersion, being equal to the difference between the third and first quartiles.  were used to quantify the central tendency and variation for the distributions of estimated biological reference points.

Results

Age composition

Length data were obtained for 12,323 sailfish (3532 females and 8791 males), and age data for 1166 of these fish (446 females and 720 males). The range of lower jaw fork length was 80-239 cm for females and 78-227 cm for males (Fig. 2). Age compositions of samples collected during the entire study period indicated that most sailfish caught off eastern Taiwan are larger than the age at 50% maturity, and peak in length at 5 years (Fig. 3). Accordingly, the estimates of age at full recruitment ([t.sub.c]) and the sample size during studying periods are shown in Table 3. The estimates of [t.sub.c] varied from 5 to 7 years for different time periods. Few sampled fishes were older than 11 years for both sexes and hence age 12 and higher were combined into the 12+ group (Fig. 3). The empirical estimates of maximum lifespan ([t.sub.max]) were 13 years for males and 21 years for females.

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

Mortality rates

The estimates of F were 0.24/yr for females and 0.43/yr for males based on the samples collected during the entire study period. Based on bootstrap analysis, the standard deviation of F was 0.046/yr and 0.045/yr for females and males, respectively. In addition, 79% of bootstrap replicates of [t.sub.c]. were 5 years old and few were 6 or 7 years old. The values of M estimated from the Pauly's empirical equation were 0.26/yr for females and 0.27/yr for males. In this study, therefore, F of 0.24/yr for females and 0.43/yr for males, [t.sub.c] of 5 years and M of 0.26/yr for females and 0.27/yr for males were set as the base case values for the subsequent analyses. Sensitivity analyses were also conducted for examining the results of Y/R and SPR analyses by assuming values of 0.20, 0.30, and 0.35 for M.

Y/R and SPR models and biological reference points

In this study, selectivities for females and males were assumed to be dome-shaped, and the estimated selectivity curves are shown in Figure 4. The estimates of [F.sub.CUR], [F.sub.0.1], and Y/R under various values of M are summarized in Table 4. For the base case, [F.sub.CUR] (0.24/yr for females and 0.43/yr for males) were substantially lower than the corresponding biological reference points [F.sub.0.1] (0.62/yr for females and 0.79/yr for males). The estimates of Y/[R.sub.CUR] were 3.37 kg for females and 3.72 kg for males and the estimates of Y/[R.sub.0.1] were 5.11 kg for females and 4.68 kg for males (Fig. 5). Even in the most conservative case, when M was assumed to be 0.2, [F.sub.CUR] was still lower than [F.sub.0.1] for both females and males.

The effect of varying [t.sub.c] on Y/R is shown in Figures 6 and 7. At low levels of F, Y/R generally increased rapidly over the range of [t.sub.c] values tested. The values of [t.sub.c] that maximized the yield per recruit decreased with the magnitude of M and increased with the level of F but typically ranged between 2 and 5 years for females and males.

The rapid growth and relatively late [t.sub.c] (5 years) of sailfish in the waters off eastern Taiwan produced low Y/R and high SPR at current fishing mortalities compared with reference points (Figs. 4 and 7). The estimates of [F.sub.SSB25], [F.sub.SSB40], SPR and Y/R under the various values of M are summarized in Table 5 and Figure 8. The base-case estimate of [F.sub.CUR] for females (0.24/ yr) was lower than the corresponding reference points [F.sub.SSB40] (0.46/yr) and [F.sub.SSB25] (0.94/yr); the current SPR was estimated to be about 57.20% of its unfished level. Under the low value of M (0.2/yr), the estimate of [F.sub.CUR] for females (0.26/yr) was lower than [F.sub.SSB40] (0.36/yr) and substantially lower than [F.sub.SSB25] (0.67/yr). Increasing [t.sub.c] to older than six years of age would ensure that the SPR was maintained at a value higher than the threshold level Noun 1. threshold level - the intensity level that is just barely perceptible
intensity, intensity level, strength - the amount of energy transmitted (as by acoustic or electromagnetic radiation); "he adjusted the intensity of the sound"; "they measured the
 (i.e., 25% of its unfished level) at almost any level of fishing mortality (Fig. 9) with relatively little effect on Y/R.

[FIGURE 4 OMITTED]

Simulation scenarios

The box plots of the estimates of [F.sub.0.1] and [F.sub.CUR]/[F.sub.0.1] are shown in Figure 10 for scenarios A-J with the assumption of uncertainty in F, M, and [t.sub.c] defined in Table 2. Although the medians of these two quantities were close to those of the base case for most scenarios, the variations of these quantities were diverse depending on the assumptions of uncertainty of the parameters. For scenarios A-D with only one parameter subject to uncertainty, the estimates of [F.sub.0.1] were independent of the changes in F but were very sensitive to the uncertainty in M. Adding the uncertainty in [t.sub.c]. resulted in slightly higher estimates of [F.sub.0.1] than that for the base case. Higher estimates of [F.sub.0.1] could have resulted from [t.sub.c] with values larger than 5 years (base case) selected in this scenario. In the case of this study, higher [t.sub.c] produced higher estimates of [F.sub.0.1] for sailfish in the water of eastern Taiwan. For scenarios E-J that incorporated the combinations of uncertainties in F, M, and [t.sub.c], large variations were observed for estimates of [F.sub.0.1] when a higher level of uncertainty was assumed for M (scenarios F, H, and J). Even though combinations of parameters with higher uncertainties were considered, the upper bounds (the third quartile Quartile

A statistical term describing a division of observations into four defined intervals based upon the values of the data and how they compare to the entire set of observations.

Notes:
Each quartile contains 25% of the total observations.
) of [F.sub.CUR]/[F.sub.0.1] were lower than 1 for all scenarios. This result implied that it was highly unlike that [F.sub.CUR] could exceed [F.sub.0.1]. In contrast, the effects of adding uncertainties in F and [t.sub.c] were relatively minor.

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

Figure 11 shows the box plots of the estimates of [F.sub.SSB40], [F.sub.SSB25], [F.sub.CUR]/[F.sub.1SSB40], and [F.sub.CUR]/[F.sub.0.1SSB25] for scenarios A-J (with uncertainty). Similarly, higher levels of uncertainties in M resulted in higher variations for the estimates of [F.sub.SSB40] and [F.sub.SSB25] (scenarios C, F, H, and J). Moreover, the estimates of [F.sub.SSB40] and [F.sub.SSB25] were obviously higher than those of the base case when the uncertainty in [t.sub.c] was considered (scenarios D, G, H, I, and J), which might result from female fish younger than [t.sub.c], being less vulnerable to fishing gear. Therefore, higher values of SSB/R would be obtained when larger values of [t.sub.c] were selected and higher estimates of [F.sub.SSB40] and [F.sub.SSB25] were revealed for these scenarios. For all scenarios, the medians of [F.sub.CUR]/[F.sub.SSB40] and [F.sub.CUR]/[F.sub.SSB25] were 0.71 and 0.37, respectively, which were similar to those of the base case (Table 5). The upper boundaries of [F.sub.CUR]/[F.sub.SSB40] and [F.sub.CUR]/[F.sub.SSB25] were substantially lower than 1 for all the scenarios and no simulation runs were observed to have the values higher than 1.

Discussion

Mortality rates

The Y/R curve is so flat-topped that maximum Y/R cannot be well defined and may not be attained under any practical fishing mortality rate (Fig. 5). Reference points such as [F.sub.0.1], [F.sub.SSB25], and [F.sub.SSB40] have often been used to develop fishery management strategies. Previous studies have indicated that Y/R and SPR are sensitive to values of M (Griffiths, 1997; Kirchner, 2001; Sun et al., 2005). In this study, however, the estimates of [F.sub.CUR] were much lower than the target levels of [F.sub.0.1] or [F.sub.SSB40], except for the scenario when M was assumed at the lower bound of 0.2. Moreover, the estimates of [F.sub.CUR] were never larger than the threshold level of [F.sub.SSB25] even when M was assumed to be 0.2. This result would indicate that the stock of sailfish in the waters off eastern Taiwan appears to be moderately exploited and has relative low risk of being overfishing Overfishing occurs when fishing activities reduce fish stocks below an acceptable level. This can occur in any body of water from a pond to the oceans. More precise biological and bioeconomic terms define 'acceptable level'. .

[FIGURE 7 OMITTED]

Sex-specific per-recruit analyses and sensitivity analyses

Most assessment methods require historic information on catch, effort, and catch-at-age (or catch-at-length). The lack of long-term fishery statistics usually makes it difficult to evaluate the status of populations exploited by small-scale fisheries or taken incidentally (e.g., Govender, 1995; Barbieri et al., 1997; Griffiths, 1997; Jones and Wells, 2001; Sun et al., 2002, 2005). Per-recruit analyses, which require only parameters related to life span and mortality, combined with an analysis of catch curves, can become an alternative method for evaluating the status of a fishery. Because per-recruit analyses can provide the estimates of biological reference points and the catch curve analysis can yield the estimate of current fishing mortality rate, the status of a fishery can be readily determined by comparing current fishing mortality with the biological reference points. Nevertheless, few results of previous studies have been discussed as to how the biological reference points based on per-recruit analyses were influenced by the uncertainties of biological parameters.

[FIGURE 8 OMITTED]

[FIGURE 9 OMITTED]

In this study, the effects of F, M, and [t.sub.c] on the estimates of [F.sub.0.1], [F.sub.SSB40], and [F.sub.SSB25] were evaluated by using the Monte Carlo simulation method. Although other parameters (e.g., growth) are essential inputs for per-recruit analyses, we focused on the effects of F, M, and [t.sub.c] which are generally more difficult to estimate owing to owing to
prep.
Because of; on account of: I couldn't attend, owing to illness.

owing to prepdebido a, por causa de 
 a lack of enough auxiliary information (Chen et al., 2007). In addition, the assumption of selectivity of fishing gear could influence the results of per-recruit analyses. Generally, longline long·line  
n.
A heavy fishing line usually several miles long and having a series of baited hooks.



long
 selectivity is assumed to be asymptotic and gillnet gill·net  
tr.v. gill·net·ted, gill·net·ting, gill·nets
To catch (fish) by means of a gill net.
 selectivity is assumed to be dome shaped. In this study, a dome-shaped selectivity was assumed for incorporating into the per-recruit analyses. Sailfish in the waters off eastern Taiwan were exploited by various fishing gears although large proportion of the catch was made by gillnet. However, insufficient length-frequency data recorded by fishing gear lead to difficulty in estimating the selectivity for different fishing gear. Therefore, collecting the information from fishing gear for length-frequency data is necessary to evaluate the influence of different selectivity assumptions on the results of assessment.

Implications of sex-specific assessment and management of the species

Uozumi (2003) indicated that the problems in the stock assessment of sailfish are the inability to obtain reliable biological parameters, standardization of catch per unit of effort, and a mechanism to develop reliable abundance indices. Sailfish are known to be sexually dimorphic dimorphic

see dimorphic fungus.
 (females grow faster and to a larger size than males (Chiang et al., 2004; Hoolihan, 2007), females become mature later than males, and sex ratio varies with length (Chiang, 2004). These attributes indicate that sex-specific assessments should be done to evaluate the status of sailfish. The sex-specific per-recruit analyses with the consideration of parameter uncertainty used in this study provide a method for evaluating management strategies.

[FIGURE 10 OMITTED]

There is a possibility that the low F estimates may be an artifact A distortion in an image or sound caused by a limitation or malfunction in the hardware or software. Artifacts may or may not be easily detectable. Under intense inspection, one might find artifacts all the time, but a few pixels out of balance or a few milliseconds of abnormal sound  of an influx of fish from outside the study area. Such a bias could occur if the probability that a fish will move from distant waters to within range of the eastern Taiwanese fleet is substantial and increases with the age of the fish (Sun et al., 2005). In principle, estimates for M could be obtained from research (such as tagging studies to determine M; Hampton, 2000). However, in the short to medium term, the values for M will have to be obtained from the results of studies for other stocks of sailfish. They have been similar to those obtained by Pauly's (1980) method. This difference should be evaluated in relation to environmental factors. Although recent analyses of molecular markers do not support recognition of separate Atlantic and Indo-Pacific species of sailfish (Graves and McDowell, 2003), there is no evidence to indicate that sailfish become increasingly likely to migrate to Taiwan with increasing age. Although Prince et al. (2006) examined the tagging results in the eastern Pacific Ocean, they considered this species a single stock. It is unclear whether the sailfish population in the Pacific Ocean comprises a single or multiple stocks, and their regional or global abundance is unknown (Ehrhardt and Fitchett, 2006). A tagging program Noun 1. tagging program - a computer program that attaches labels to the grammatical constituents of textual matter
tagger

computer program, computer programme, programme, program - (computer science) a sequence of instructions that a computer can interpret
 in which electronic and conventional tags are used to examine the spatial movement patterns and stock structure of sailfish in this geographical region would prove beneficial for the sustainable management of the species.

There are no management measures for sailfish in the waters off eastern Taiwan at present. In the waters off eastern Taiwan, sailfish are targeted by the gillnet fishery and caught incidentally in the longline, harpoon harpoon (härpn`), weapon used for spearing whales and large fish. The early type was a flat triangular piece of metal with barbed edges and a socket for attaching a wooden handle, to the , and set net fisheries. This makes it difficult to effectively control the fishing effort of these fisheries for sailfish. In this case, the isopleths An isopleth is a feature of meteorological charts, connecting points which have an equal value of some variable at a given time and spatial area. The particular variables shown on a given chart may include values such as pressure, temperature, wind speed, etc.  of SPR indicate that increasing [t.sub.c] to an age between six and seven years old would likely result in modest gains in terms of SPR (on the order of 40%) and also hedge against recruitment overfishing. For example, the target level of 40% SPR could be achieved at even double values of [F.sub.CUR] if [t.sub.c] for female sailfish was larger than five years (the age at sexual maturity). Of course the efficacy of increasing [t.sub.c] would be mitigated by any substantial release mortality. At present little is known about the mortality rates of fish released from gillnet and longline vessels operating off eastern Taiwan, and further study is needed before we can confidently recommend increasing [t.sub.c] as a measure to prevent overfishing.

[FIGURE 11 OMITTED]

In summary, sex-specific per-recruit modeling coupled with Monte Carlo simulation analyses are effective in evaluating the stock status of billfish because of the sexual dimorphism Sexual dimorphism

Any difference, morphological or behavioral, between males and females of the same species. In many animals, the sex of an individual can be determined at a glance.
 and uncertainty in key life history and fishery parameters of these species. This study reveals that sailfish in the waters off eastern Taiwan appear to be moderately exploited and have relative low risk of being overfished. However, in view of the recent rapid increase in fishing effort, it is evident that the stock status and development of the fishery need to be closely monitored.

Acknowledgments

We would like to thank two reviewers, Michael H. Prager, NMFS NMFS National Marine Fisheries Service
NMFS National Mortality Followback Survey
NMFS Network Multimedia File System
NMFS Nested Mount File System
 Southeast Fisheries Science Center, and Gerard DiNardo, NMFS Pacific Islands Fisheries Science Center, and one anonymous reviewer for their valuable and constructive comments. We thank Michael Prager for proposing the modified catch-curve method. This study was partially financially supported by the Fisheries Agency, Council of Agriculture, Taiwan, through the grant 93AS-9.1.1-FA-F1(2) to Chi-Lu Sun.

Manuscript submitted 17 July 2008. Manuscript accepted 28 January 2009. Fish. Bull. 107:265-277.

The views and opinions expressed or implied in this article are those of the author and do not necessarily reflect the position of the National Marine Fisheries Service, NOAA NOAA
abbr.
National Oceanic and Atmospheric Administration

Noun 1. NOAA - an agency in the Department of Commerce that maps the oceans and conserves their living resources; predicts changes to the earth's environment;
.

Literature cited

Barbieri, L. R., M. E. Chittenden Jr., and C. M. Jones. 1997. Yield-per-recruit analysis and management strategies for Atlantic croaker Atlantic croaker
n.
A small silvery food fish (Micropogonias undulatus) common in Atlantic waters south of Massachusetts.

Noun 1.
, Micropogonias undulates, in the Middle Atlantic Adj. 1. middle Atlantic - of a region of the United States generally including Delaware; Maryland; Virginia; and usually New York; Pennsylvania; New Jersey; "mid-Atlantic states"
mid-Atlantic
 Bight bight, broad bend or curve in a coastline, forming a large open bay. The New York bight, for example, is the curve in the coast described by the southern shore of Long Island and the eastern shore of New Jersey. The term bight may also refer to the bay so formed. . Fish. Bull. 95: 637-645.

Butterworth, D. S., A. E. Punt, D. L. Borches, J. G.. Pugh, and G. S. Hughes. 1989. A manual of mathematical techniques for linefish assessment (incorporating report of the SANCOR Marine Linefish Programme's Workshop on Population Dynamics Population dynamics is the study of marginal and long-term changes in the numbers, individual weights and age composition of individuals in one or several populations, and biological and environmental processes influencing those changes. , February 4-6, 1987, Cape Town Cape Town or Capetown, city (1991 pop. 854,616), legislative capital of South Africa and capital of Western Cape, a port on the Atlantic Ocean. It was the capital of Cape Province before that province's subdivision in 1994. ). S. Afr. Natl. Sci. Prog. Rep. 160, 89 p.

Chen, Y., L. Xu, X. Chen, and X. Dai. 2007. A simulation study of impacts of at-sea discarding and bycatch on the estimation of biological reference points [F.sub.0.1] and [F.sub.max]. Fish. Res. 85:14-22.

Chiang, W. C. 2004. Population dynamics and stock assessment of the sailfish (Istiophorus platypterus) in waters off eastern Taiwan. Ph.D. diss diss  
v.
Variant of dis.


diss
Verb

Slang, chiefly US to treat (a person) with contempt [from disrespect]

Verb 1.
., 171 p. National Taiwan University, Taipei, Taiwan.

Chiang, W. C., C. L. Sun, S. Z. Yeh, and W. C. Su. 2004. Age and growth of sailfish (Istiophorus platypterus) in waters off eastern Taiwan. Fish. Bull. 102:251-263.

Chiang, W. C., C. L. Sun, S. Z. Yeh, W. C. Su, D. C. Liu, and W. Y. Chen. 2006. Sex ratios, size at sexual maturity, and spawning seasonality of sailfish Istiophorus platypterus from eastern Taiwan. Bull. Mar. Sci. 79:727-737.

Deriso, R. B. 1987. Optimal F0.1 criteria and their relationship to maximum sustainable yield In population ecology, maximum sustainable yield or MSY is, theoretically, the largest yield/catch that can be taken from a species' stock over an indefinite period. . Can. J. Fish. Aquat. Sci. 44:339-348.

Ehrhardt, N. M., and M. D. Fitchett. 2006. On the seasonal dynamic characteristics of the sailfish, Istiophorous platypterus, in the eastern Pacific off central America Central America, narrow, southernmost region (c.202,200 sq mi/523,698 sq km) of North America, linked to South America at Colombia. It separates the Caribbean from the Pacific. . Bull. Mar. Sci. 79:589-606.

Gabriel, W. L., M. P. Sissenwine, and W. J. Overhaltz. 1989. Analysis of spawning stock biomass per recruit: an example for Georges Bank Georges Bank

Submerged sandbank in the Atlantic Ocean east of Massachusetts, U.S. It has long been an important fishing ground, with scallops harvested in its northeastern portion. Navigation is made dangerous by crosscurrents and fog.
 Haddock haddock: see cod.
haddock

Valuable North American food fish (Melanogrammus aeglefinus, family Gadidae). A bottom-dweller that feeds on invertebrates and fishes, it resembles the cod, with its chin barbel (fleshy feeler) and two anal and three dorsal
. N. Am. J. Fish. Manag. 9:383-391.

Goodyear, C. P. 1993. Spawning stock biomass per recruit in fisheries management: foundation and current use. Can. Spec. Publ. Fish. Aquat. Sci. 120:67-81.

Govender, A. 1995. Mortality and biological reference points for the king mackerel mackerel, common name for members of the family Scombridae, 60 species of open-sea fishes, including the albacore, bonito, and tuna. They are characterized by deeply forked tails that narrow greatly where they join the body; small finlets behind both the dorsal and  (Scomberomorus commerson) fishery off Natal Natal, city, Brazil
Natal (nətäl`), city (1991 pop. 606,887), capital of Rio Grande do Norte state, NE Brazil, just above the mouth of the Potengi River.
, South Africa South Africa, Afrikaans Suid-Afrika, officially Republic of South Africa, republic (2005 est. pop. 44,344,000), 471,442 sq mi (1,221,037 sq km), S Africa.  (based on a per-recruit assessment). Fish. Res. 23:195-208.

Graves J. E., and J. R. McDowell. 2003. Stock structure of the world's istiophorid billfishes: a genetic perspective. Mar. Freshw. Res. 54:287-298.

Griffiths, M. H. 1997. The application of per-recruit models to Argyrosomus inodorus, an important South African sciaenid fish. Fish. Res. 30:103-115.

Gulland, J. A., and L. K. Boerema. 1973. Scientific advice on catch levels. Fish. Bull. 71:325-335.

Hampton, J. 2000. Natural mortality in tropical tunas: size really does matter. Can. J. Aquat. Sci. 55:1002-1010.

Hilden, M. 1993. Reference points for fisheries management: the ICES experience. Can. Spec. Publ. Fish. Aquat. Sci. 120:59-66.

Hoolihan, J. P. 2005. Horizontal and vertical movements of sailfish (Istiophorus platypterus) in the Arabian Gulf Arabian Gulf: see Persian Gulf. , determined by ultrasonic ultrasonic /ul·tra·son·ic/ (-son´ik) beyond the upper limit of perception by the human ear; relating to sound waves having a frequency of more than 20,000 Hz.

ul·tra·son·ic
adj.
1.
 and pop-up satellite tagging. Mar. Biol. 146:1015-1029.

2006. Age and growth of Indo-Pacific sailfish, Istiophorus platypterus, from the Arabian Gulf. Fish. Res. 78:218-226.

Jones, C. M., and B. K. Wells. 2001. Yield-per-recruit analysis for black drum The black drum (Pogonias cromis) is a saltwater fish similar to its cousin, the red drum. It is the only species in the genus Pogonias. Though most specimens are generally found in the 5-30lb range, the black drum is well known as the largest of all the drum , Pogonias cromis, along the East Coast of the United States The "Eastern Seaboard," or "Atlantic Seaboard" are terms referring to the easternmost coastal states in the United States. They touch the Atlantic Ocean and stretch up to Canada.  and management strategies for Chesapeake Bay Chesapeake Bay, inlet of the Atlantic Ocean, c.200 mi (320 km) long, from 3 to 30 mi (4.8–48 km) wide, and 3,237 sq mi (8,384 sq km), separating the Delmarva Peninsula from mainland Maryland. and Virginia. . Fish. Bull. 99:328-337.

Katsukawa, T., Y. Lee, and Y. Matsumiya. 1999. Spawning per recruit analysis for female snow crab Chionoecetes opilio in the Sea off Kyoto prefecture. Nippon Suisan Gakk. 65:288-293.

Kirchner, C. H. 2001. Fisheries regulations based on yield-per-recruit analysis for the linefish silver kob Argyrosomus inodorus in Namibian waters. Fish. Res. 52:155-167.

Leaman, B. M. 1993. Reference points for fisheries management: the western Canadian experience. Can. Spec. Publ. Fish. Aquat. Sci. 120:15-30.

Liu, K. M., Y. T. Chang, I. H. Ni, and C. B. Jin. 2006. Spawning per recruit analysis of the pelagic thresher The pelagic thresher, Alopias pelagicus, is a thresher shark of the genus Alopias, found in tropical oceans, at depths down to 500 m. It was first described in 1935 by Japanese ichthyologist Hiroshi Nakamura.  shark shark, member of a group of almost exclusively marine and predaceous fishes. There are about 250 species of sharks, ranging from the 2-ft (60-cm) pygmy shark to 50-ft (15-m) giants. They are found in all seas, but are most abundant in warm waters. , Alopias pelagicus, in the eastern Taiwan waters. Fish. Res. 82:56-64.

Myers, R. A., and B. Worm. 2003. Rapid worldwide depletion of predatory fish communities. Nature 423:280-283.

Pauly, D. 1980. On the interrelationships between natural morality, growth parameters and mean environmental temperatures in 175 fish stocks. J. Cons. Int. Mer 39:175-192.

Prager, M. H., and C. P. Goodyear. 2001. Effects of mixed-metric data on production model estimation: simulation study of a blue-marlin-like stock. T. Am. Fish. Soc. 130:927-939.

Prince, E. D., D. B. Holts, D. Snodgrass, E. C. Orbesen, J. Luo, M. Domeier and J. E. Serafy. 2006. Transboundary and movement of sailfish, Istiophorus platypterus, off the Pacific coast of Central America. Bull. Mar. Sci. 79:827-838.

Punt, A. E. 1997. The performance of VPA-based management. Fish. Res. 29:217-243.

Ricker, W. E. 1975. Computation and interpretation of biological statistics offish population. Bull. Fish. Res. Board Can. 191,382 p.

Rudershausen, P.J., E.H. Williams, J.A. Buckel, J.C. Potts, and C.S. Manooch III. 2008. Comparison of reef fish catch per unit effort and total mortality between the 1970s and 2005-2006 in Onslow Bay Onslow Bay is an indentation of the North Carolina coast, between Cape Fear in the south and Cape Lookout in the north. Thirteen barrier islands form the shore of the bay. It is part of the open Atlantic Ocean and affords no protection from oceanic swells. , North Carolina North Carolina, state in the SE United States. It is bordered by the Atlantic Ocean (E), South Carolina and Georgia (S), Tennessee (W), and Virginia (N). Facts and Figures


Area, 52,586 sq mi (136,198 sq km). Pop.
. Tran. Am. Fish. Soc. 137:1389-1405.

Sun, C.L., N.M. Ehrhardt, C.E. Porch, and S.Z. Yeh. 2002. Analyses of yield and spawning stock biomass per recruit for the South Atlantic albacore albacore: see tuna.
albacore

Large oceanic tuna (Thunnus alalunga) that is noted for its fine flesh. The streamlined bodies of these voracious predators are adapted to fast and continuous swimming.
 (Thunnus alalunga). Fish. Res. 56:193-204.

Sun, C. L., S. P. Wang, C. E. Porch, and S. Z. Yeh. 2005. Sex-specific yield per recruit and spawning stock biomass per recruit for the swordfish, Xiphias gladius, in the waters around Taiwan. Fish. Res. 71:61-69.

Taylor, C. C. 1958. Cod growth and temperature. J. Cons. Int. Exp. Mer 23:360-370.

Uozumi, Y. 2003. Historical perspective of' global billfish stock assessment. Mar. Freshw. Res. 54:555-565.

Watanabe, K., T. Hesho, K. Saiura, T. Okazake, and Y. Matumiya. 2000. Fisheries management by spawning per recruit analysis for the threeline grunt Parapristipoma trilineatum around Mugi-oshima Island on the coastal area of Tokushima prefecture. Nippon Suisan Gakk. 66:690-696.

Wei-Chuan Chiang (1)

Chi-Lu Sun (contact author) (2)

Sheng-Ping Wang (3)

Su-Zan Yeh (2)

Yong Chen (4)

Wei-Cheng Su (5)

Don-Chung Liu (5)

Wen-Yie Chen (1)

Email address See Internet address.  for contact author: chilu@ntu.edu.tw

(1) Eastern Marine Biology marine biology, study of ocean plants and animals and their ecological relationships. Marine organisms may be classified (according to their mode of life) as nektonic, planktonic, or benthic. Nektonic animals are those that swim and migrate freely, e.g.  Research Center of Fisheries Research Institute No. 22 Wuchuan Rd. Chenkung, Taitung, Taiwan 961

(2) Institute of Oceanography National Taiwan University No. 1, Sec. 4, Roosevelt Rd. Taipei, Taiwan 106

(3) Department of Environmental Biology and Fisheries Science National Taiwan Ocean University No. 2, Beining Rd. Keelung, Taiwan 202

(4) School of Marine Science University of Maine "UMO" redirects here, but this abbreviation is also used informally to mean the Mozilla Add-ons website, formerly Mozilla Update

Should not be confused with Université du Maine, in Le Mans, France
The University of Maine
 218 Libby Hall Orono, Maine 04469

(5) Fisheries Research Institute No. 199, Ho-lh Rd. Keelung, Taiwan 202
Table 1

Biological parameters used in the per-recruit analysis for the
sailfish (Istiophorus platypterus) in the waters off eastern
Taiwan during the period from July 1998 to July 2005. VBGF is
the von Bertalanffy growth function, [L.sub.[infinity]] = the
asymptotic length, K =the growth coefficient, [t.sub.0] = the
hypothetical age at length zero; Length-weight relationship is
W = A x [L.sup.B] , where W = rought weight (in kg) and L = lower
jaw fork length (in cm); Maturity fraction parameters [r.sub.m] =
the slope of logistic equation fitted to the maturity data
collected, and [t.sub.m] = the age at 50% sexual
maturity.

Parameter                              Female

VBGF
  [L.sub.[infinity]]           250.29 cm
  K                              0.138/yr
  [t.sub.0]                     -2.99

Length-weight relationship
  A                              2.3234 x [10.sup.-6]
  B                              3.1013

Maturity fraction
  [r.sub.m]                      1.525
  [t.sub.m]                      5

Parameter                               Male

VBGF
  [L.sub.[infinity]]           240.539 cm
  K                              0.145/yr
  [t.sub.0]                     -2.781

Length-weight relationship
  A                              1.1933 x [10.sup.-5]
  B                              2.7828

Maturity fraction
  [r.sub.m]
  [t.sub.m]

Table 2

Scenarios designed to examine the effects of uncertainty of F, M, and
[t.sub.c] on the estimates of biological reference points of sailfish
(Istiophorus platypterus) in the waters off eastern Taiwan. Scenarios
A-D were used to evaluate the effects of changes in a single
parameter, scenarios E-H were used to evaluate the results of changes
in combinations of two parameters, and scenarios I and J were used
to evaluate the results when three parameters were subject to
uncertainty. [[sigma].sub.F] = standard deviation from the
boot-strapped estimation; M = natural mortality per year;
[t.sub.c] = age at first catch)

                             Parameters

            [[sigma].sub.F]              M

Scenario    Female     Male     Female        Male       [t.sub.c]

Base        0         0        0.26         0.27            5
A           0.046     0.045    0.26         0.27            5
B           0         0        0.2-0.3      0.2-0.3         5
C           0         0        0.15-0.35    0.15-0.35       5
D           0         0        0.26         0.27            5-7
E           0.046     0.045    0.2-0.3      0.2-0.3         5
F           0.046     0.045    0.15-0.35    0.15-0.35       5
G           0         0        0.2-0.3      0.2-0.3         5-7
H           0         0        0.15-0.35    0.15-0.35       5-7
I           0.046     0.045    0.2-0.3      0.2-0.3         5-7
J           0.046     0.045    0.15-0.35    0.15-0.35       5-7

Table 3

Sample sizes and estimates of age at full recruitment
for female and male sailfish (Istiophorus platypterus)
by year in the waters off eastern Taiwan during the
period from July 1998 to July 2005.

              Sample          Age at
               size      full recruitment

Female
  1998          131             6
  1999          702             5
  2000          387             7
  2001          358             6
  2002          349             7
  2004          484             5
  2005         1121             6
  Overall      3532             5

Male
  1998          525             6
  1999         1943             5
  2000          578             7
  2001          730             5
  2002          576             5
  2004         1123             5
  2005         3316             5
  Overall      8791             5

Table 4

Estimates of current fishing mortality ([F.sub.CUR]), current yield
per recruit (Y/[R.sub.CUR]), and the reference points of [F.sub.0.1]
and Y/[R.sub.0.l] at different levels of natural mortality (M) for
female and male sailfish (Istiophorus platypterus) in the waters
off eastern Taiwan during the period from July 1998 to July 2005.
[F.sub.0.1] = the fishing mortality rate corresponding to
Y/[R.sub.0.1]; Y/[R.sub.0.1] = the point of a yield-per-recruit
curve where the slope equals 10%n of the slope at the origin.

M        [F.sub.CUR]   [F.sub.0.1]   Y/[R.sub.CUR]   Y/[R.sub.0.1]
(1/yr)     (1/yr)        (1/yr)          (kg)             (kg)

Female
  0.20      0.26          0.54           5.36             7.08
  0.26      0.24          0.62           3.37             5.11
  0.30      0.22          0.69           2.45             4.13
  0.35      0.20          0.80           1.63             3.19

Male
  0.20      0.46          0.70           5.79             6.66
  0.27      0.43          0.79           3.72             4.68
  0.30      0.41          0.83           3.06             4.04
  0.35      0.38          0.91           2.21             3.17

Table 5

Estimates of current fishing mortality ([F.sub.CUR]) and the
reference points of SSB/[R.sub.SSBx] and SPR at different levels
of natural mortality (M) for female sailfish (Istiophorus
platypterus) in the waters off eastern Taiwan during the period
from July 1998 to July 2005. [F.sub.SSBx] = the fishing mortality
rate corresponding to SSB/[R.sub.SSBx]; SSB/[R.sub.SSBx] = a
specific percentage (x%) of the spawning biomass per recruit at
the unfished level; SPR (spawning potential ratio) = the SSB/R
at a given fishing mortality divided by the
SSB/R without fishing.

M            [F.sub.CUR]         [F.sub.SSB40]      [F.sub.SSB25]
(1/yr)         (1/yr)               (1/yr)             (1/yr)

0.20            0.26                 0.36               0.67
0.26            0.24                 0.46               0.94
0.30            0.22                 0.54               1.22
0.35            0.20                 0.70               1.73

M          SSB/[R.sub.CUR]     SSB/[R.sub.SSB40]
(1/yr)          (kg)                 (kg)

0.20            32.50                26.56
0.26            22.22                15.53
0.30            17.42                11.14
0.35            12.97                 7.52

M         SSB/[R.sub.SSB25]      [SPR.sub.CUR]
(1/yr)          (kg)                  (%)

0.20            16.60                48.96
0.26             9.71                57.20
0.30             6.96                62.57
0.35             4.70                68.95
COPYRIGHT 2009 National Marine Fisheries Service
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2009 Gale, Cengage Learning. All rights reserved.

 Reader Opinion

Title:

Comment:



 

Article Details
Printer friendly Cite/link Email Feedback
Author:Chiang, Wei-Chuan; Sun, Chi-Lu; Wang, Sheng-Ping; Yeh, Su-Zan; Chen, Yong; Su, Wei-Cheng; Liu, Don-C
Publication:Fishery Bulletin
Article Type:Report
Geographic Code:9TAIW
Date:Jul 1, 2009
Words:6684
Previous Article:Age- and length-at-maturity of female arrowtooth flounder (Atheresthes stomias) in the Gulf of Alaska.
Next Article:Abundance, condition, and diet of juvenile Pacific ocean perch (Sebastes alutus) in the Aleutian Islands.
Topics:

Terms of use | Copyright © 2014 Farlex, Inc. | Feedback | For webmasters