Estimation of Maximum Sustainable Yield Using Production Modeling: A Stock Appraisal of Indian Oil Sardine (Sardinella longiceps) from Pakistani Waters
Estimation of Maximum Sustainable Yield Using Production Modeling: A Stock Appraisal of Indian Oil Sardine (Sardinella longiceps) from Pakistani Waters
Abdul Baset1, Qun Liu1,*, Muhammad Tariq Hanif1, 2, Baochao Liao1, Aamir Mahmood Memon1,3 and Muhammad Mohsin1
1College of Fisheries, Ocean University of China, Qingdao 266003, China
2Marine Fisheries Department, Fish Harbor West Wharf, Karachi 74000, Pakistan
3Sindh Fisheries Department, Hyderabad 71000, Sindh, Pakistan
Time series of catch and effort data were analyzed to estimate the Maximum Sustainable Yield (MSY) of Indian oil sardine, Sardinella longiceps (Valenciennes, 1847, Family: Clupeidae), fishery resources from Pakistani waters. The catch and effort data (1994 – 2009) were analyzed by CEDA (catch and effort data analysis) consists of Fox, Schaefer and Pella Tomlinson models and ASPIC (a surplus production model incorporating covariates) contain Fox and Logistic models. The average annual landings of S. longiceps stocks were 38 353 t while the highest and lowest 65 050 t in 1994 and 26 937 t in 2009 was recorded correspondingly. The MSY was estimated using the initial proportion (IP) value of 1 because the starting catch was approximately 100% of the maximum catch. The estimated values of MSY using CEDA with the Fox model under three error assumptions of normal, log-normal and gamma were 25 860.07 t (R2 = 0.864), 26 812.65 t (R2 = 0.836) and 26 665 t (R2 = 0.85), respectively. The MSY values with Schaefer and Pella Tomlinson models under two error assumptions normal and log-normal were 29 160.13 t (R2 =0.84) and 30 455.5 t (R2 =0.807), respectively. The values of Schaefer and Pella Tomlinson models were the same. The gamma error assumption showed minimization failures for Schaefer and Pella Tomlinson models. The estimated values of MSY using ASPIC computer software with two surplus production models of Fox and logistic were 28 180 t (R2 =0.915) and 33 260 t (R2 =0.897) respectively. The estimated values of MSY from CEDA were about 26 000 - 31 000 t and from ASPIC were 29 000 - 31 000 t. The estimated MSY values of CEDA and ASPIC were higher than the most recent catch. This shows that the stocks are in a safe condition.
Received 29 January 2016
Revised 26 April 2016
Accepted 03 August 2016
Available online 02 February 2017
QL conceived and designed the study. MTH arranged the data. BL and AMM helped in data analysis. MM prepared the Figures. AB analyzed the data and wrote the article. QL reviewed and edited the article.
Sardinella longiceps, Maximum sustainable yield, Catch and effort data analysis (CEDA), A surplus production model incorporating covariates (ASPIC), Pakistan.
* Corresponding author: firstname.lastname@example.org
0030-9923/2016/0004-1161 $ 8.00/0
Copyright 2016 Zoological Society of Pakistan
Fisheries resources are not only playing a significant role in the national economy, but also take part in the human development and welfare, such as in terms of providing a quality source of food and employment with the major aspect of self-renewable in nature. Therefore, if the fish stocks were managed in a sustainable manner then their effects can be limitless ().
The coastal belt of Pakistan is about 1,120 km long from the southeast Indian border to the northwest Iranian border () with an EEZ (exclusive economic zone) of 240,000 km2 with an additional continental shelf area of about 50,270 km2. The EEZ of the country is rich in fish and mineral resources (). There were about 250 commercially important demersal fishes (), 50 small pelagic, 15 medium-sized pelagic and 20 large pelagic fish species in Pakistani water (). Small pelagic fishes have a cosmopolitan distribution in the world seas. These species were caught, canned and consumed globally (). These resources mainly include species from family Clupeidae, Engraulidae, Scomboridae and Carangidae. Indian oil sardine (Sardinella longiceps Valenciennes, 1847) belongs to the family Clupeidae and local names are Tarli, Luar in Sindhi and Lugger, Luar in Baluchi language (). Body length has estimated about15-17 cm at the end of first year and 19 cm at the end of second year (). Its maximum length is 23 cm and commonly 16 cm in Pakistan ().
Oil sardine, S. longiceps is a valuable commercial fish, used for food as well as for the production of fishmeal and oil (). Several studies are related on the biological and ecological aspects of Indian oil sardine (; , ; ).
In the present paper we analyzed the catch and effort data of S. longiceps from Pakistan to estimate maximum sustainable yield (MSY). Several surplus production models have been used to appraise the fish stocks for the estimation of MSY (e.g. ; ; ; ; ; ). These models express the fishery yields as a function of catch and fishing effort. Some work has already been done on the estimation of MSY of various fish species from Pakistan and China (e.g. , ; ; ; ).
The aim of this article is to provide baseline information on the current stock status of Indian oil sardine (S. longiceps) in Pakistani waters. For this purpose, we estimated the MSY biological reference points from a variety of surplus production models (; ; ) with different error assumptions by catch and effort data analysis (CEDA) () and a non-equilibrium surplus production model ASPIC (a surplus-production model incorporating covariates, ) software. The estimated values may help to draw the desired fisheries regulations for the sustainable utilization of this fishery resource in the country.
Note: Catch is in the form of metric ton.
MATERIALS AND METHODS
The catch and effort data of Indian oil sardine (Sardinella longiceps) for the period of 1994 to 2009 (16 years) were from the handbook of Fisheries Statistics of Pakistan compiled by Marine Fisheries Department (MFD), Karachi, Pakistan (). Fishing effort is presented by the number of operational fishing boats in the maritime region of Pakistan, and the per annum total catch is presented in the form of catch weight (metric tons).
Surplus production models
The fishery statistics of S. longiceps were analyzed by CEDA and ASPIC. In order to minimize the biasness in the results we used both of them which were developed by fisheries scientists from UK and US. Both have some technical differences but they are based on the same theory. There were three surplus production models (SPMs) also called biomass dynamic models (BDMs) used in CEDA package i.e. Schaefer, Fox and Pella-Tomlinson. The most commonly used model is which is based on the logistic population growth model:
CEDA is a pc-based software package, designed by fishery scientists in UK for fish stock assessments, which was first developed in 1992 and has been adjusted and further developed on request of users.
CEDA, version 3.0.1 includes three non-equilibrium surplus production models (; ; ) with three error assumption (normal, log-normal and gamma). The computer package can calculate the following key parameters: MSY, q (catchability coefficient), K (carrying capacity), r (intrinsic growth rate), replacement yield, final biomass, whereas CV (coefficient of variation) of the estimated MSY values can also be calculated. The package requires an input value of initial biomass (B1) over the carrying capacity by the user. When the IP (initial proportion) value has set to zero, it points out that the fishery started from a virgin population and if the IP value is close to 1, it indicates that the fishery initiate from a heavily harvested population. However, in some cases starting biomass is fixed by programmer i.e.
ASPIC is a non-equilibrium application of the surplus production models which was developed by a fishery scientist in USA and started in 1994 and is freely available from the internet. The model is more formally described in and .
ASPIC package contains two surplus production models Logistic (Schaefer) and Fox (a special case of GENFIT). At the beginning, the program was based on the logistic model () but the most current version 5.0 allows the selection criterion between logistic and generalized production models. The main output parameters are: MSY, q, K, ratio of starting biomass over carrying capacity (B1/K), coefficient of determination (R2), coefficient of variation (CV), stock biomass giving MSY (BMSY), and fishing mortality rate at MSY (FMSY).
CEDA is a computer program which requires catch and effort data with three non-equilibrium surplus production models and three error assumptions (). The average catch of S. longiceps from Pakistani waters in the 16-year period of 1994 – 2009 was 38,353 t. shows that the observed catch in 1994 was 65,050 t which was the highest catch while in 2009 it was 26,937 t which was the lowest catch. In and , bootstrapping confidence limit method was used for obtaining the CV. The results from the CEDA package are highly responsive to the input IP (initial proportion) values (). The estimated values of MSY were different when IP from 0.1 to 1. showed the CEDA results with IP of 1 because the initial yield was about 100% of the maximum catch. The estimated values of MSY with CV from the Fox model with three error assumptions (normal, log normal and gamma) were produced at 25,860 t (CV = 0.060), 26,812 t (CV = 0.064) and 26,665 t (CV = 0.051), respectively. Similarly, the estimated values of MSY for Schaefer and Pella Tomlinson with two error assumptions normal and log normal were 29,160 t (CV = 0.101), (CV=0.088) and 30,455 t (CV = 0.084), (CV=0.092) respectively. The Fox model produced R2 = 0.864, 0.836 and 0.85, respectively, whereas Schaefer and Pella Tomlinson models produced R2 =0.84 and 0.807, respectively. The values of MSY were the same from Schaefer and Pella Tomlinson. Gamma error assumption showed minimization failure for Schaefer and Pella Tomlinson. There was a little difference among the estimated values of MSY from all models, i.e. Fox, Schaefer and Pella Tomlinson while estimated R2 values show a better fit to the data.
|Normal||Log- normal||Gamma||Nor mal||Log-normal||Gam ma||Nor mal||Log-normal||Gam ma|
|(0.048)||(2.34E- 05)||MF||(5688.9)||(4.74E- 05)||MF||(7473)||(6.7E- 05)||MF|
MF, showing minimization failures.
|Fox (Log normal)||322858||0.231573||27504.57||26812||0.0643||2.22E-05||0.836||93147.27|
|Schaefer (Log normal)||268659.6||0.453444||30455.5||23880||0.08494||2.70E-05||0.807||71915.17|
|Pella Tomlinson (Normal)||290645.5||0.401315||29160.13||22656.||0.0884||2.46E-05||0.84||76692.88|
|Pella Tomlinson (Log normal)||268659.6||0.453444||30455.5||2388||0.0921||2.70E-05||0.807||71915.17|
|Pella Tomlinson (Gamma)||MF||MF||MF||MF||MF||MF||MF||
Note: MF stands for minimization failure.
The parameters are: K, carrying capacity; r, intrinsic population growth rate; MSY, maximum sustainable yield; Ryield, replacement yield; CV, Coefficient of variation; q, catchability coefficient; R2, coefficient of determination and biomass.
ASPIC () was also used for the estimation of MSY and related parameters of the Indian oil sardine fishery from Pakistani waters. Results from the non-equilibrium surplus production models of Logistic and Fox in ASPIC computer package are shown in . The IP was 1, because the initial yield was about 100% of the maximum catch. The estimated values of MSY from Fox and logistic with CV were 28,180 t (CV=0.039) and 33,260 t (CV=0.084), respectively. The rate of fishing mortality at MSY (FMSY) and stock biomass giving MSY (BMSY) model was 0.268 and 105,100 t from the Fox and 0.329 and 100,900t from the Logistic model correspondingly, while the estimated values of K (carrying capacity) in the Fox model were a little higher than in the Logistic model. In addition to this, shows ASPIC results for the
MSY, maximum sustainable yield; B1/K, starting biomass over carrying capacity; K,carrying capacity; FMSY , fishing mortality rate at MSY; q, catchability coefficient; BMSY, stock biomass giving MSY; R2, coefficient of determination.
|0.1||52020||3.15E -0||2.23E +05||1.22E-04||6.33E-01||0.813||25.10636||82190|
|0.2||34800||5.57E -0||2.93E +05||4.51E-05||3.23E-01||0.907||0.229254||107700|
|0.3||28180||1.01E +0||2.86E +05||2.52E-05||2.68E-01||0.915||0.036071||105200|
|0.4||28180||1.01E +0||2.86E +05||2.52E-05||2.68E-01||0.915||0.035486||105200|
|0.7||28180||1.01E +0||2.86E +05||2.52E-05||2.68E-01||0.915||0.034121||
|0.8||28180||1.01E +0||2.86E +05||2.52E-05||2.68E-01||0.915||0.04855||105100|
|0.9||28180||1.01E +0||2.86E +05||2.52E-05||2.68E-01||0.915||0.039971||
|1||28180||1.01E +0||2.86E +05||2.52E-05||2.68E-01||0.915||0.039191||105100|
S. longiceps fishery in Pakistan with IP ranging from 0.1 to 1. The program does not show higher sensitivity to initial proportion values as compared to CEDA computer package. The MSY values of Logistic model were high when IP = 0.1 and 0.2, the remaining were between 28,000 t to 50,000 t from both Fox and Logistic models which were relatively close to the annual catch of an Indian oil sardine fishery. The Fox model produced R2 =0.915 while Logistic produced R2 =0.897, which were higher than those from the CEDA computer package and showed the excellent fitting to the data.
Decisions derived to manage the fisheries are commonly directed by the stock assessment results. Therefore, it is essential that scientists provide a reliable description of stock dynamics and stock status to the managers (). The main purpose of the paper was to estimate the MSY of S. longiceps by Surplus Production models. The results obtained from CEDA and ASPIC in this study were basically consistent among the three surplus production modeling approaches i.e. Fox, Schaefer and Pella Tomlinson. Logistic approach resulted in a higher value of MSY which showed that the Fox model is more conservative.
Catch-per-unit-effort (CPUE) is considered to be proportional to the fish population abundance and referred as the relative abundance index. Several population dynamic models were used to predict the relative abundance index in order to obtain the future values of predicted absolute abundance by multiplying with a constant (catchability coefficient, q) (). To do so, the SPMs or BDMs are the plausible approach due to their holistic or simple in nature. The surplus production models can determine the level of effort at which a fishery produces a maximum yield of a fish stock in a sustainable manner without changing the long term productivity termed as MSY ().
Assumptions, limitations and trends
The perception or understanding of MSY has played a significant role in the fishery science over the decades and often regarded as the target Biological Reference Point (BRPs) (), which serve as a standard of measurement of stock status from a biological point of view. These reference points are generally summarized under the three main fish stock dynamic models, i.e., stock-recruit, dynamic pool and production models ().
Fish stock sizes and distributions can fluctuate widely even in their unexploited state, due to the change in environmental factors and also the interrelationships with other species. Some researchers have indicated that if the surplus production is higher than catch it means the population size is increasing, whereas when the catch is greater than surplus production the population has a decreasing trend. Moreover, when the catch and production were equal, it showed a constant trend in fish population size (; ; ).
Small pelagic fishes accounted for 33.3% of the total marine fishery resources of Pakistan, including Herrings, Sardines, and Anchovies. During the year 1998, the total biomass and MSY of small pelagic fishes was estimated to be 700,000 and 300,000 t, respectively (). However, previously no work has been reported on the MSY of Indian oil sardine (S. longiceps) separately. To our knowledge, this is the first study to investigate the MSY of sardine fisheries from Pakistani water’s. We have found that the trajectory of stock biomass and exploitation rate remains steady for years 1999-2009 except for few minimization failures only in the gamma error assumption.
Tools of analysis
CEDA package is limited to three surplus production models, i.e., Fox, Schaefer, Pella-Tomlinson with three error assumptions, i.e., Normal, Log-Normal and Gamma respectively () from which we can estimate the MSY, and other related population parameters of fish stocks. It can be used as a supportive tool for fishery management. When we use an initial proportion (IP) ranging from 0.1–1 shows that CEDA package is sensitive to the IP values. When IP values were low, the estimated MSY values were higher and when IP values were high the estimated MSY values were lower (). When the IP was from 0.5 - 1, the estimated MSY values from Fox, Schaefer and Pella-Tomlinson production models were lower than that of observed catch data. To minimize the bias in the data we use IP= 1 because the starting catch was roughly 100% of the maximum catch. The estimated values from CEDA package were 25,860 t, 26,812 t and 26,665 t for the Fox model with three error assumptions normal, log normal and gamma respectively, while 29,160 t and 30,455 t for Schaefer and Pella Tomlinson with two error assumptions normal and log normal, respectively ().
Similarly, the ASPIC () was used to estimate the values of MSY and related parameters from time series catch and fishing effort data by non-equilibrium surplus production models. The non-equilibrium surplus production models were often used because the fish stocks were rarely in equilibrium state due to several factors affecting the population size i.e., biology, environment and fishing effort. The MSY values and other parameters of Indian oil sardine fishery from Pakistan were estimated and shown in which indicates that the ASPIC package is not sensible with IP values. The Fox and Logistic models produced R2 =0.915 and 0.897, respectively. In contrast to CEDA, the ASPIC showed a better fit to the data.
The estimated values of MSY from CEDA package were about 26,000 – 31,000 t and from ASPIC the values were 26,000 - 33,000 t .These values are higher than the recent catches of S. longiceps from Pakistani waters. From the outcome of our investigation, it is possible to conclude that the Indian oil sardine fishery in Pakistan were in managed condition.
Abdul Baset acknowledges the Chinese Scholarship Council (CSC) for funding his PhD. degree. This work is supported by the Fundamental Research Funds for the Central Universities (No. 201562030).
Statement of conflict of interest
Authors have declared no conflict of interest.
Bianchi, G., 1985. FAO species identification sheets for Fishery Purposes. Field guide to the commercial marine and brackish-water species of Pakistan. Prepared with the support of PAK/77/033 and FAO (FIRM) Regular Programme. Rome, FAO, 8p.
Gabriel, W.L. and Mace, P.M., 1999. A review of biological reference points in the context of the precautionary approach. Proceedings of the fifth national NMFS stock assessment workshop: providing scientific advice to implement the precautionary approach under the Magnuson-Stevens fishery conservation and management act. NOAA Tech Memo NMFS-F/SPO-40: 34-45.
Hoggarrth, D.D., Abeyasekera, S., Arthur, R.I., Beddington, J.R., Burn, R.W., Halls, A.S, Kirkwood, G.P., McAllister, M., Medley, P., Mees, C.C., Parkes, G.B., Pilling, G.M., Wakeford, R.C. and Welcomme, R.L., 2006. Stock assessment for fishery management-a framework guide to the stock assessment tools of the fisheries management science programme. FAO Fish. Tech. Pap. 487. FAO, Rome, Italy. pp. 261.
Kalhoro, M.A., Liu, Q., Memon, K.H., Chang, M.S. and Jatt, A.N., 2013. Estimation of maximum sustainable yield of Bombay duck, Harpodon nehereus fishery in Pakistan using the CEDA and ASPIC packages. Pakistan J. Zool., 45: 1757-1764.
Maunder, M.N., Sibert, J.R., Fonteneau, A., Hampton, J., Kleiber, P. and Harley, S.J., 2006. Interpreting catch per unit effort data assess the status of individual stocks and communities. ICES J. mar. Sci., 63: 1373-1385.
Panwar, S.K., Liu, Q., Khan, F. and Siddiqui, P.J.A., 2012a. Maximum sustainable yield estimates of Ladypees, Sillago sihama (Forsskål), fishery in Pakistan using the ASPIC and CEDA packages. J. Ocean Univ. China, 11: 93-98.
Panwer, S.K., Liu, Q., Khan, F. and Waryani, B., 2012b. Maximum sustainable yield estimates of spiny lobster fishery in Pakistan using non-equilibrium CEDA package. Russian J. mar. Biol., 38: 448-453.
Prager, M.H., 2005. A stock-production model incorporating covariates (Version 5) and auxiliary programs. CCFHR (NOAA). Miami Laboratory Document MIA-92/93-55, Beaufort Laboratory Document BL-2004-01.
Smith, T. and Punt, A., 2001. The gospel of maximum sustainable yield in fisheries management: birth, crucifixion and reincarnation. Conservation of exploited species (eds. J.D. Reynolds, G.M. Mace, K.H. Redford, and J.G. Robinson), Cambridge University Press, Cambridge, UK: pp. 41-66.