Maximum sustainable yield of Greater lizardfish Saurida tumbil fishery in Pakistan using the CEDA and ASPIC packages

The catch and effort data analysis（CEDA） and ASPIC（a stock assessment production model incorporating covariates） computer software packages were used to estimate the maximum sustainable yield（MSY） from the catch and effort data of Greater lizardfish Saurida tumbil fishery of Pakistan from 1986 to 2009. In CEDA three surplus production models of Fox, Schaefer and Pella-Tomlinson were used. Here initial proportion（IP） of 0.5 was used because the starting catch was roughly 50% of the maximum catch. With IP = 0.5, the estimated MSY from Fox model were 20.59 mt and 38.16 mt for normal and log-normal error assumptions, while the MSY from Schaefer and Pella-Tomlinson were 60.40, 60.40 and 60.40 mt, for normal, log-normal and gamma error assumptions respectively. The MSY values from Schaefer and Pella-Tomlinson models of three error assumptions were the same. The R2 values from those three models were above 0.6. When IP = 0.5, the MSY values estimated from ASPIC from Fox were 132 mt, and from logistic model were 69.4 mt, with R2 value above 0.8. Therefore we suggest the MSY of S. tumbil fishery from Pakistan to be 60–70 mt, which is higher than the latest catch, thus we would recommend that the fishing efforts for this fishery may be kept at the current level.

Population dynamics of Japanese threadfin bream Nemipterus japonicus from Pakistani waters

Japanese threadfin bream Nemipterus japonicus（Bloch, 1791） is among the most abundant and commercially important species in Pakistan. From the coast of Pakistan, four demersal trawl surveys in October–November 2009 and May–June, August, October and November in 2010 were carried out. The purpose of this study is to estimate the population dynamics and status of the stock of the N. japonicus from Pakistani waters based on the research trawl surveys from the research area. The data consist of n=784 length-weight pairs and n=7 530 length frequency with the maximum length and weight of 29 cm and 358 g respectively. The length frequency data were analyzed using ELEFAN method in FiSAT computer package. The parameters of length and weight relationship were b=2.778, a=0.032 and R2=0.973. The estimated von Bertalanffy growth function parameters were L∞=30.45 cm, K=0.270 year-1. Based on length-converted catch curve analysis the total mortality（Z） during this study was estimated at 0.960 year-1. The natural mortality coefficient（M） was 0.74 year-1 using Pauly＇s equation（the annual average sea surface temperature was 27°C）, therefore, the fishing mortality coefficients（F） were 0.22 year-1. The yield per recruit analysis indicated that when tc was 2, Fmax was estimated at 1.2 and F0.1 at 1.1. When tc was 1, Fmax was estimated at 0.95 and F0.1 at 0.8. Because current age at first capture is about 1 year and Fcurrent was 0.22, Fcurrent is smaller than F0.1 and Fmax, which indicated that the fishery is about in a safe condition. When using Gulland（1971） biological reference point, Fopt was equals to M（0.74）. The current fishing mortality rate of 0.22 was smaller than the target biological reference point.