First Assessment of the Thryssa vitrirostris ( Engraulidae ) Beach Seine Fishery in Northeastern Mozambique

Monthly length-frequency data, from 2009 to 2014, was used to estimate the growth parameters, mortality, and spawning season, and to assess the exploitation status of Thryssa vitrirostris (Engraulidae) fisheries in Pebane. The von Bertalanffy asymptotic length (L∞) and growth rate (K) were 25.1 cm (standard length) and 0.41 per year (standard length), respectively. Two proxy recruitment peaks were found: the first peak of recruitment occurs from April to July, and the second recruitment peak from September to October. The total estimated mortality rates (Z), natural mortality (M), and fishing mortality (F) were Z = 1.31, M = 0.92, and F = 0.39. For the beach seine gear, the size at first catch of T. vitrirostris was Lc25 = 4.43 cm. The 50% retention size of the catch was Lc50 = 5.39 cm. The retention probability analyses revealed a large rate of juvenile fishing mortality (54.2%). The estimated exploitation rate (0.30) was below the maximum exploitation rate (0.48), and above the optimal sustainable exploitation rate (E50 = 0.28), evidencing a sustainable fishery. However, under such an exploitation regime, it is advised that a continuously monitoring-survey of T. vitrirostris is maintained. An increase in migration of fishermen has been recently recorded in Pebane, due to its rich fisheries, which can increase the fishing effort and the risk of overexploitation if management measures (such as mesh size increase) are not taken in advance.


Introduction
The Mozambique's fisheries are divided in three subsectors: industrial, artisanal, and subsistence.The artisanal and subsistence fisheries contribute to the large amount of landings, and more than 2/3 of landings in the economic exclusive zone [1].Small-scale and subsistence small-medium pelagic fisheries have a major socio-economic role for coastal communities, comprising an important direct protein food source [2].The small-scale and subsistence fishing takes place both from shore and from canoes and dhow-type planked boats, mostly propelled by sails [3]), and they almost exclusively target the nearshore waters, up to 40 m depth or less [4].Doherty et al. [5] reconstructed Mozambique catches and estimated them to be between 55,000 and 64,000 t/year in 1950s, and between 120,000 and 130,000 t/year by the late 2000s.According to the National Institute of Fisheries Research (IIP) statistics, on average, the annual catch for the period of 2003-2006 of small-scale fishes was around 40,000 t/year [6,7]).Given the great potential available in terms of small pelagic fish, it is believed that this estimated value is only what has been reported by the fisheries authorities, since most of them are not reported [1].In fact, Mozambique has a history of high under-reporting of catches, mainly due to the illegal, unreported, and undeclared (IUU) catches from the artisanal and subsistence fishing sector.Artisanal fishermen often fish during the closed season, and in protected areas when it is not allowed [8].In the country, there are around 120,000 fishers and 658 small-scale coastal landing sites.National fish stocks are under great pressure, and many fisheries show signs of overexploitation and the shallow coastal waters are severely overfished [1].
Most marine species, including small pelagic fish, are distributed along to the Sofala bank, characterized by an extensive and wide shallow continental shelf with a high marine biodiversity [9].In Sofala bank, the Engraulidae catches are very abundant, and Thryssa vitrirostris is considered a particularly important fishery resource [6].Based on spatial factors, seasonal catch rate, and demographic data from the beach-seine fishery, it is considered that the large fishery of the species in Mozambique occurs in the Sofala bank region.The core of the exploitable population is determined in Pebane, and the species seems to migrate towards the coast at the onset of the rainy season (November-February), and then northwards [6].T. vitrirostris is a small pelagic fish that belongs to the Engraulidae family, inhabiting coastal and estuarine waters [10] and forming large schools or shoals at depths ranging from 0 to 50 m (www.fishbase.org,9 January 2018).It distributes throughout the Indian Ocean, Madagascar and the coast of Africa from the Alfredo port to the north of the Persian Gulf, and along the coast of Pakistan and India ( [11], www.fishbase.org,9 January 2018).
According to [5,12], T. vitrirostris represents the largest source of income for most of the population that lives along the coast, as well as for the country's economy, providing crucial social benefits to these communities.The few studies carried out in Mozambique related to this species have reported greater representation in artisanal catches, as well as in terms of bycatch in industrial and semi-industrial demersal shrimp fisheries, which represent the most frequent and most important fisheries in Mozambique, and so this aspect make this species an important source of animal protein for coastal communities, playing an important role for the diet and local economy [7,13,14].
The few studies focused on T. vitrirostris in Mozambique are: (i) the study of [14,15], which is based on data that has been collected through scientific cruises to address aspects of the biology of the species in general; (ii) the study of [16], which addresses aspects relating to temporal distribution, growth, and the reproductive biology of T. vitrirostris on Zalala beach; (iii) and lastly, the study developed by [6] that addresses the distribution and biology of T. vitrirostris and other Engraulidae along the coast of Sofala bank, western Indian ocean.This latter study is one of the main references currently, in terms of information about this species, namely concerning the analysis of the influence of fishing mortality (F), on the dynamics of the explored stock along the Sofala bank.Nevertheless, due to continuously inconsistent time series records, such as demographic data, most studies that assess T. vitrirostris are based on short-term and/or episodic data.T. vitrirostris fishing activity is socio-economically important and it assures daily livelihood incomes and food.The monthly length frequency data obtained between 2009 and 2014 by the Institute of Fisheries Research were used to determine the biological growth parameters, mortality rates, and recruitment season, and to assess the fishery exploitation status of T. vitrirostris.

Growth Parameters
The minimum and maximum length-classes were 4 and 25 cm.The asymptotic length (L ∞ ) values required to seed the Von Bertalanffy Plot (VBP) estimated by the different techniques are given in Table 1.
The L ∞ , used for "seed" in the VBF equation, was estimated in Electronic Length-Frequency Analysis (ELEFAN) using the response surface routine.As given in the Table 1, the estimated (L ∞ ) was L ∞ = 25.1 cm.The instantaneous growth rate (K) was K = 0.41 year −1 (goodness of fit: Rn = 0.15).
The monthly von Bertalanffy growth curve of T. vitrirostris was superimposed over the normal length frequency histograms and over the restructured length frequency histograms (Figure 1).The L∞, used for "seed" in the VBF equation, was estimated in Electronic Length-Frequency Analysis (ELEFAN) using the response surface routine.As given in the Table 1, the estimated (L∞) was L∞ = 25.1 cm.The instantaneous growth rate (K) was K = 0.41 year −1 (goodness of fit: Rn = 0.15).The monthly von Bertalanffy growth curve of T. vitrirostris was superimposed over the normal length frequency histograms and over the restructured length frequency histograms (Figure 1).The analyses of the VBF histogram revealed several cohorts along the year.Overall, six cohorts were observed across years (Figure 2).The length-classes between 12 and 15 cm comprised most of the catches (32.52%), with the length-modal class set at 13 cm (Figure 2).The analyses of the VBF histogram revealed several cohorts along the year.Overall, six cohorts were observed across years (Figure 2).The length-classes between 12 and 15 cm comprised most of the catches (32.52%), with the length-modal class set at 13 cm (Figure 2).

Recruitment Patterns
According to the recruitment pattern analyses, the species had two peaks of recruitment over the year.The analyses of the recruitment pattern (Figure 3A), together with length frequency histogram observations (black slices in Figure 1B) revealed that the first peak of recruitment occurred from April to July, and the second recruitment peak from August to October.The analyses of the percentage of recruits of young of the year (4 and 5 cm) revealed a peak in May and in June (Figure 3B,C).

Recruitment Patterns
According to the recruitment pattern analyses, the species had two peaks of recruitment over the year.The analyses of the recruitment pattern (Figure 3A), together with length frequency histogram observations (black slices in Figure 1B) revealed that the first peak of recruitment occurred from April to July, and the second recruitment peak from August to October.The analyses of the percentage of recruits of young of the year (4 and 5 cm) revealed a peak in May and in June (Figure 3B,C).

Recruitment Patterns
According to the recruitment pattern analyses, the species had two peaks of recruitment over the year.The analyses of the recruitment pattern (Figure 3A), together with length frequency histogram observations (black slices in Figure 1B) revealed that the first peak of recruitment occurred from April to July, and the second recruitment peak from August to October.The analyses of the percentage of recruits of young of the year (4 and 5 cm) revealed a peak in May and in June (Figure 3B,C).

Mortality Parameters
The instantaneous total mortality (Z) rate, determined with the length-converted catch curve, estimated a mortality rate of Z = 1.31 year −1 (Figure 4).The current fishing exploitation rate estimated (E est.), derived from the analysis of mortality rates, was E est.= 0.30.The estimated instantaneous natural mortality (M) coefficient was M = 0.92.The determination of the instantaneous fishing mortality coefficient (F) was F = 0.39.
J 2018, 2, x FOR PEER REVIEW 5 of 16

Mortality Parameters
The instantaneous total mortality (Z) rate, determined with the length-converted catch curve, estimated a mortality rate of Z = 1.31 year −1 (Figure 4).The current fishing exploitation rate estimated (Eest.), derived from the analysis of mortality rates, was Eest.= 0.30.The estimated instantaneous natural mortality (M) coefficient was M = 0.92.The determination of the instantaneous fishing mortality coefficient (F) was F = 0.39.

Probabilities of Capture
The length at the first catch (Lc25) and the length at which the probably of 50% (Lc50) and 75% (Lc75) of individuals are captured by gear, determined through the cumulative probability of capture, (Figure 5) were the following: Lc25 = 4.43 cm, Lc50 = 5.93 cm, and Lc75 = 6.35 cm.The percentage of juveniles (length-classes < 13 cm) comprised 54.2% of the catches, while adults were 45.2%.No statistically significant differences in yearly percentage in number average, were observed (t-test: t = 0.95; df = 5; p = 0.384) between juveniles and adults.

Probabilities of Capture
The length at the first catch (Lc 25 ) and the length at which the probably of 50% (Lc 50 ) and 75% (Lc 75 ) of individuals are captured by gear, determined through the cumulative probability of capture, (Figure 5) were the following: Lc 25 = 4.43 cm, Lc 50 = 5.93 cm, and Lc 75 = 6.35 cm.The percentage of juveniles (length-classes < 13 cm) comprised 54.2% of the catches, while adults were 45.2%.No statistically significant differences in yearly percentage in number average, were observed (t-test: t = 0.95; df = 5; p = 0.384) between juveniles and adults.
J 2018, 2, x FOR PEER REVIEW 5 of 16

Mortality Parameters
The instantaneous total mortality (Z) rate, determined with the length-converted catch curve, estimated a mortality rate of Z = 1.31 year −1 (Figure 4).The current fishing exploitation rate estimated (Eest.), derived from the analysis of mortality rates, was Eest.= 0.30.The estimated instantaneous natural mortality (M) coefficient was M = 0.92.The determination of the instantaneous fishing mortality coefficient (F) was F = 0.39.

Probabilities of Capture
The length at the first catch (Lc25) and the length at which the probably of 50% (Lc50) and 75% (Lc75) of individuals are captured by gear, determined through the cumulative probability of capture, (Figure 5) were the following: Lc25 = 4.43 cm, Lc50 = 5.93 cm, and Lc75 = 6.35 cm.The percentage of juveniles (length-classes < 13 cm) comprised 54.2% of the catches, while adults were 45.2%.No statistically significant differences in yearly percentage in number average, were observed (t-test: t = 0.95; df = 5; p = 0.384) between juveniles and adults.

Yield Isopleths
The estimated Lc/L∞ ratio for the currently estimated exploitation regime (Eest.) was Lc/L∞ = 0.24.

Discussion
The Mozambique coast is subject to intense fishing pressure by artisanal fishing [1].Generally, no scientific information is available to support fishery management or ecosystem ecological-based management, with consequent socio-economic and biodiversity impacts.This increases the impact of a rise in poverty, under the current scenario of an increase in population along estuarine and coastal areas [17].The analysis of length-frequency data is a reliable way to obtaining information regarding fishery status in tropical areas, and in underdeveloped poor countries where continuously monitoring programs often fail.Therefore, this study can contribute to support fishery management and sustainability of one of the most important small-pelagic subsistence/artisanal coastal fisheries,  As fishing exploitation increased relative yield per recruit (Y'/R) decreased (Figure 6).The Beverton and Holt relative yield per recruit (model Figure 6) estimated: E 10 = 0.40, E 50 = 0.28, E max.= 0.48 (Figure 6A).

Yield Isopleths
The estimated Lc/L∞ ratio for the currently estimated exploitation regime (Eest.) was Lc/L∞ = 0.24.

Discussion
The Mozambique coast is subject to intense fishing pressure by artisanal fishing [1].Generally, no scientific information is available to support fishery management or ecosystem ecological-based management, with consequent socio-economic and biodiversity impacts.This increases the impact of a rise in poverty, under the current scenario of an increase in population along estuarine and coastal areas [17].The analysis of length-frequency data is a reliable way to obtaining information regarding fishery status in tropical areas, and in underdeveloped poor countries where continuously monitoring programs often fail.Therefore, this study can contribute to support fishery management and sustainability of one of the most important small-pelagic subsistence/artisanal coastal fisheries,

Yield Isopleths
The estimated Lc/L ∞ ratio for the currently estimated exploitation regime (E est. ) was Lc/L ∞ = 0.24.

Discussion
The Mozambique coast is subject to intense fishing pressure by artisanal fishing [1].Generally, no scientific information is available to support fishery management or ecosystem ecological-based management, with consequent socio-economic and biodiversity impacts.This increases the impact of a rise in poverty, under the current scenario of an increase in population along estuarine and coastal areas [17].The analysis of length-frequency data is a reliable way to obtaining information regarding fishery status in tropical areas, and in underdeveloped poor countries where continuously monitoring programs often fail.Therefore, this study can contribute to support fishery management and sustainability of one of the most important small-pelagic subsistence/artisanal coastal fisheries, and consequently quality of life (socioeconomic) in the region.To the authors' knowledge, the monthly time-series used in this study was one of the largest datasets available for Thryssa vitrirostris (Engraulidae) for stock status assessment.As no discards were observed, we considered that the landings were representative of total catches.
To obtain a reliable estimate of growth and biological parameters of the population, the suitability of length-frequency data must be ascertained.The raw length-frequency data should exhibit peaks with apparent shifts in modal length over time [18,19].As revealed by yearly length-frequency monthly plots, the data met these criteria.The estimated growth parameters (L ∞ = 25.1 cm; K = 0.41) differed from those found in a nearby beach at Zalala (L ∞ = 22.26;K = 0.44; [19]), and in grouped data of Sofala bank (L ∞ =19; K = 0.66; [6]), with the asymptotic length being markedly higher, while the growth rate was lower.According to [20,21], the asymptotic length (L ∞ ) was basically influenced by food availability and population density, while the growth rate (K) was a parameter that was dependent on genetic and physiological factors, which can vary according to environmental fluctuations.
Environmental conditions, and namely, sea surface temperature (SST) can affect fish growth and recruitment, namely in small pelagic species [15,22,23]).According to [24], oceanographic conditions, such as SST, upwelling, current regimes, and wind conditions at Pebane and nearby areas (present study) were similar to environmental conditions recorded in other fisheries studies in the region [19,24].Thus, biological parameter differences (L ∞ and K) estimated among studies should not be related to putative environmental causes.
Previous studies on the biology of T. vitrirostris [19,24] were based on episodic length frequency data collected through discontinuous years and/or areas.In the present study, the observed maximum length-class sampled was 25 cm, which was close to the L ∞ estimations.Thus, we believe that the growth parameter estimated in the present work, using continuous monthly surveys, can contribute to enhancing the knowledge of the species biology and support Mozambique fisheries managers.The L ∞ value estimated corresponded to an increase in previous maximum L ∞ off 2.5 cm in northeastern Mozambique [6,19].However, comparing to [13], the L ∞ herein estimated was 0.7 cm below the L ∞ value estimated in the Maputo region (Southeaster Mozambique).According to [21], the asymptotic growth (L ∞ ) and growth rate (K) parameters were inversely related, which means that the greater the asymptotic growth (L ∞ ), the lower the growth rate (K), and the higher the growth rate (K), the smaller the asymptotic growth.Overall the L ∞ value was higher than in other studies conducted in northeastern Mozambique, while the growth rate value (K) was lower [7,19].The scientific cruise of [13] was conducted with a pelagic trawl.The range of the length-size in [25] did not include length classes below 8 cm, nor higher than 22 cm.The size of the trawl mesh was not mentioned in the report, although we can assume (as standard in a scientific survey) that a small mesh size was used in the cod end.Thus, it is a little speculative to debate whether the growth parameter differences among studies are due to fishing gear, since gear selectivity and fish size retained are linked.However, environmental conditions were conservative.Therefore, we cannot exclude the probability that gear or fishing technique differences could explain the regional differences in parameter estimations.
The estimated mortality rates in northeastern Mozambique (Z = 1.31,M = 0.92 and F = 0.39) were also different from those found by [26] (Z = 2.30, M = 1, and M = 1.5), with the M value depending on the method used by the authors, respectively.Consequently, fishing mortality (F) estimated based on the difference between F = Z − M would be 1 and 0.8 for each M value estimated by [26].Therefore, it can be questioned why growth rate (K) and mortality rates in the present study are lower than recorded elsewhere in other close areas (for instance, mortality rate values were almost at half, compared to similar studies).The observed differences may be associated with the methods that were used to analyze the size-classes frequency data and the range size of the length across studies.For instance, [26] estimated mortality using Monte Carlos simulations (to cope with scarce data), while in the present work, a six year monthly size-class time-series dataset was used.Another reason can be also linked to the data acquisition areas/surveys.[26] sampled offshore areas and pooled demographic information (inshore and offshore) for the estimation of biological parameters and stock assessment models, including mortality rate.Nevertheless, the observed maximum length found by these authors and [19] were lower than those found in this study.In fact, the demography differences found among studies are considerable, with modal length-class values of the present study being higher comparatively to former studies.Therefore, new estimates of mortality achieved are more accurate, because the monthly data have captured a more complete population size distribution.
VBF analysis clear showed several cohorts over the year.In fact, the monthly modal progression analysis represented by histograms reveals a clear upward trend in the growth of the observed cohorts.This is related to the growth rate profile in tropical areas, which is continuous throughout the year, whilst in temperate areas, the growth rate shows a seasonal trend [27][28][29].According to [26], the main spawning activity of T. vitrirostris extends from November to January, and July to August.However, latter study does not cover all of the months in the year.In the present study, the recruitment analysis based on species demography matches the results obtained by [26] regarding the species recruitment season.The recruitment analysis showed two clear proxies recruitment peaks: the first peak of recruitment occurs from April to July, and the second recruitment peak from September to October.The ontogeny period of small pelagic larvae species is considered to be short [30], with the warm temperatures contributing to reduce the larvae stage duration [31].Considering a lag among spawning periods, larvae growth, and the fishing recruitment of young of the year fish (length-classes of 4-5 cm), a time lag of 1.5 to 2.5 months between spawning and fish gear recruitment can be determined.
The natural mortality rate (M) was higher than the fishing mortality (F), when it was partly expected to be the opposite trend if the schooling social behavior of species that form groups, making them more vulnerable to fishing mortality, was accounted for.However, natural mortality is related to the biological characteristics of the species [32].Small pelagic fish species usually present a fast growth and a relatively short life cycle and higher natural mortality.The tropical fish tend to have high natural mortality rates (M) in relation to fishing mortality (F), with the natural mortality (M) not evidencing any relationship with the asymptotic size (L ∞ ) or the growth rate (K) [27,33].This can be related to the fact that, unlike fishing mortality, natural mortality is associated with predation and diseases, two factors that are not related to the age of individuals [33].
The probability of capture analysis recorded a standard selectivity logistic curve for purse seine fisheries.In fact, more than 50% of undersized fish are retained above the 6 cm size-class.Considering the size of first maturity is 13 cm [26], our results showed a large juvenile mortality for beach seine gear.The percentage of juvenile fish in catch (size-classes) comprises 54.2% of the catches.According to the current fishing legislation the allowable minimum mesh size to beach seining net used to capture this species is 30 mm.For this mesh size the estimated first capture length was Lc = 4.43 cm and the retention size of 50% of the capture is Lc 50 = 5.39 cm.Thus, the mesh size used may be inadvisable to exploit the resource if we consider that the first maturation size is 13 cm, which is much larger than the size of the Lc 50 , indicating that it is necessary to adapt the biology of the species (maturity size) to the selectivity of the fishing gear.
The estimated fishing exploitation rate value (E est.= 0.30) was below the E max.value (E max.= 0.48) but above the optimal exploitation rate value (E 50 = 0.28).This means that the exploitation pressure is above the value of E under which the stock has been reduced to 50% of its unexploited biomass.According to [32,33], a fishing stock is considered to be at a sustainable exploitation level when the exploitation rate does not exceed 50% (E 50 ), the point at which natural mortality (M) and fishing mortality (F) are at equilibrium.Recent studies confirmed significant increases in the number of fishing gear and fishermen into Zambézia Province, Pebane, due to higher fishing yields expected to be obtained in this area [6].So, some concerns arise from the future of T. vitrirostris artisanal fishing exploitation regimes in Sofala bank, namely in Pebane, an area where species forms a single demographic population and where the core of the exploitable population is found [6].
Along Mozambique coastal fishing communities, T. vitrirostris is captured with beach seine, and no discards occurs.Thus landings are representative of the catch.The fishing gear has low selectivity, as shown by the analysis of the capture probability.It is generally recognized that in areas where food security is poor, fisheries tend to be overexploited and too high a proportion of juveniles are caught, placing populations at risk.Yet T. vitrirostris population seems healthy and sustainably produces recruits and catches.Such results are related with fishing exploitation regime.However, considering the E est.= 0.3 and Lc/L ∞ = 0.24 values, the production rate falls in a exploitation regime that requires some concern on the part of the fishery managers [34], namely, continuously monitoring of the fishery to ensure maintained fishing mortality rate at a steady state.Another concern is the amount of catches under small mesh sizes that can increase the risk of overexploitation [34].If fishing managers considered maintaining the fishing exploitation rate as it stands (E max.> E est.> E 50 ), an increase in mesh size can be enforced in advance as a precaution management approach aiming directly at reduce fishing mortality.In several parts of the world, small pelagic species under continuous stock assessment and enforced regulations have been reported to be at biological risk, due to high exploitation levels [26,35].Taken together, and considering the socio-economic importance of T. vitrirostris, the currently exploitation regime might require the implementation of management measures to avoid biomass reduction to unsustainable levels.At the current exploitation stage, it is essential to continue with monitoring surveys and evaluate the risk associated with fishing effort increases as fishing precautionary approaches.

Study Area
The Mozambique Channel is located between Madagascar and the African continent (11−27 • S).The coastline of Mozambique extends over 2700 km in length, and it is oriented from north to south.
The study area has located to the north of the Sofala bank, in Pebane (Figure 7).The continental shelf in the northern region of Mozambique is shallow and narrow, and it extends only a few kilometers offshore until the 200 m isobath that is found close to the coast [23].Along the central region (Pebane region) the Mozambican coast is characterized by a wide continental shelf with the 200 m isobath much further offshore [23].The climate in the region is tropical humid, with two seasons throughout the year; the summer-hot and rainy season, which varies from November to April, and the winter-dry and cold period, which varies from May to September [36].Over the coast, the equatorial current flows from the north, and it follows towards the south forming cyclonic and anticyclonic vortices [36,37].The tides are semi-diurnal with a diurnal inequality and the amplitude pattern does not exceed 4 m [37].The data for the study was obtained along Pebane fishing villages/communities.The sampling sites include several estuarine (Molocue-Cuassiane; Malaua) and coastal fishing villages (Morremone; Maverane, Macuacuane; Cuassiane; Therrepuane; Sacone).

Data Collection
Monthly size frequency data was obtained by the technicians of Instituto Nacional de Investigação Pesqueira (IIP-Mozambique).Due to the lack of fisheries statistics on the species, the main fishing centers and/or fishing villages in Pebane (N = 8) were initially randomly chosen and visited on a monthly basis from February 2009 to November 2014 (six years).The data acquired during the monthly field sampling surveys in all communities over mentioned were stored in the PESCART database, used by the [38].The sampling process targeted fishes that had been caught by beach seine, an artisanal fishery, which was the main gear used by local fishermen.From the beach seine catches, a sample of Thryssa vitrirostris was measured, and the individual total length recorded (measured from the upper jaw to the tip of the caudal fin), allowing a size frequency distribution by length-classes to be obtained.Data was grouped into 1 cm length-classes rounded to the nearest mm.The sampling of the specimens was conducted in the beach, and fish were afterwards commercialized.All fishes caught in beach seine were used size-independently for food consumption.Therefore, the technicians did not need to consider sampling before discards, as landings were equivalent to the total catch.Samples of Thryssa spp.from within a single catch were randomly collected.

Data Collection
Monthly size frequency data was obtained by the technicians of Instituto Nacional de Investigação Pesqueira (IIP-Mozambique).Due to the lack of fisheries statistics on the species, the main fishing centers and/or fishing villages in Pebane (N = 8) were initially randomly chosen and visited on a monthly basis from February 2009 to November 2014 (six years).The data acquired during the monthly field sampling surveys in all communities over mentioned were stored in the PESCART database, used by the [38].The sampling process targeted fishes that had been caught by beach seine, an artisanal fishery, which was the main gear used by local fishermen.From the beach seine catches, a sample of Thryssa vitrirostris was measured, and the individual total length recorded (measured from the upper jaw to the tip of the caudal fin), allowing a size frequency distribution by length-classes to be obtained.Data was grouped into 1 cm length-classes rounded to the nearest mm.The sampling of the specimens was conducted in the beach, and fish were afterwards commercialized.All fishes caught in beach seine were used size-independently for food consumption.Therefore, the technicians did not need to consider sampling before discards, as landings were equivalent to the total catch.Samples of Thryssa spp.from within a single catch were randomly collected.
The size of the gear varied from 100 to 600 m in length.In the central part between the ends, the gear varied from 6 to 20 m in height, depending on the total length; the height of the net varied from 2 to 10 m.Overall, the difference between height of the center and height at the ends of the net, allowed for the formation of a catch bag or bag, where the captured fish was retained.Currently the mesh size fixed in the regulation of this fishery was 3 cm.The modus operandi was the following: (i) the net is thrown into the sea with the help of a canoe at one end, while the other end is secured or otherwise held by one of the crew members of the canoe that lies on the ground; (ii) the launching of the net must follow the orientation of an arc towards the margin with a certain distance between the cables, in the sense of having a considerable opening; (iii) the opening of the network is related to the slope of the beach where it operates; this means that the greater the slope of the area, the greater the difference.The size of the gear varied from 100 to 600 m in length.In the central part between the ends, the gear varied from 6 to 20 m in height, depending on the total length; the height of the net varied from 2 to 10 m.Overall, the difference between height of the center and height at the ends of the net, allowed for the formation of a catch bag or bag, where the captured fish was retained.Currently the mesh size fixed in the regulation of this fishery was 3 cm.The modus operandi was the following: (i) the net is thrown into the sea with the help of a canoe at one end, while the other end is secured or otherwise held by one of the crew members of the canoe that lies on the ground; (ii) the launching of the net must follow the orientation of an arc towards the margin with a certain distance between the cables, in the sense of having a considerable opening; (iii) the opening of the network is related to the slope of the beach where it operates; this means that the greater the slope of the area, the greater the difference.

Data Analysis
T. vitrirostris monthly length-frequency data were analyzed using the computational package FISAT II, FAO-ICLARM Stock Assessment Tools [39].The monthly size data were analyzed to determine important life history parameters that would inform the biology and exploitation status of T. vitrirostris.

Growth Parameters
The T. vitrirostris growth parameters were estimated by means of ELEFAN (Electronic Length Frequency Analysis) [40].This method allowed the growth curves from fish length-frequency data to be followed, sequentially arranged in time.The method allowed an objective estimation of growth parameters whenever representatively samples were available from a fish population.The method was suited both to the analysis of species growth in tropical (such as Quelimane area, in Northeastern Mozambique) and temperate stocks [40].With ELEFAN I, data were reconstructed to generate "peaks" and "troughs", and the best statistical estimation of the growth parameters (K and L ∞ ) was selected based on the Rn value, that measure the goodness of fit (the higher Rn, the better the fitting).The distribution of length frequencies in histograms was adjusted and analyzed, based on the seasonal growth curves of von Bertalanffy, since this technique serves to adjust the observed mean length at age for different combined cohorts.
The asymptotic growth (L ∞ ) and growth rate (K) were estimated from the seasonal equation of von Bertalanffy, developed by [41], and later modified by Somers [26,27] as the following: where: L t -is the length at time t, L ∞ -asymptotic length, K-growth rate, t-time, t o -the age at zero length, C-amplitude of growth oscillations, t s -the time between birth and onset of the first growth oscillations.
According to the von Bertalanffy function (VBF), fish grow (length in time L t ) on average towards the asymptotic length (L ∞ ) at an instantaneous growth rate (K) [42].This means that growth is based on the modal displacement in time sequences of the length samples [31,43].There are several techniques that are available, that were tested to estimate the L ∞ required for the "seed" L ∞ in the VBF Equation (1) when the K-Scan procedure is used: (i) The Powell-Wetherall plot method was used to obtain an initial estimate of (L ∞ = 27.11cm).
This method takes into consideration the distribution of length frequency, pooled against the midlength [26,44], allowing for verification of the behavior or variation of the individual sizes in relation to L ∞ .The mathematical model of this method is described below: L middle -is the mean length of all fish (≥L ) and was calculated as: L is the cutoff length, the length that once reached, the individual has maximum probability of being caught by fishing, L ∞ = −a/b and Z/K = −(1 + b)/b.(ii) L ∞ required for "seed" L ∞ in VBF Equation (1) was determined based on the maximum length ratio observed (L max./0.95) in the original length frequency data (maximum observed length-class = 25 cm) according to [45].(iii) At last, the L ∞ value from length frequency data was estimated using the response value routine in ELEFAN.Using the latter approach (iii), and removing the two higher length-classes, that showed an irregular distribution (frequency occurrence very low), we were able to get a good fit for the VBF growth Equation ( 1) and for both mortality rates and catch plot curves.

Growth Parameters (L ∞ ) and Model Selection
According to previous works, the asymptotic length of the species ranged between 19 [6,7] and 25.8 cm [15].Raw data showed a maximum length-class of 25 cm.Therefore, initially the L ∞ was estimated considering all length-classes, from 4 to 25 cm.The L ∞ values estimated using (i), (ii), and (iii) allowed for a good adjustment of the VBF growth Equation (1).L ∞ values estimated using in (i) and (ii) with all length-classes did not perform well for subsequent analyses, namely, the length converted-catch curves and the mortality rates estimations.This can be related with the low frequency of occurrence and the low number of individuals, for large length-class size distribution [39,43].In fact, data exploitation analyses reveal that the two higher-length classes (24 and 25 cm) showed an irregular distribution (frequency occurrence that was very low).By removing the 24 and 25 cm length classes, it was possible to achieve a good adjustment of the VBF growth equation (Rn and r-coefficient of correlation), and both mortality rates and catch plot goodness of fit in all analyses.Taken together, the best model or the best estimated L ∞ value (statistical fit) necessary to seed the VBF equation to carry out subsequent analyses, was the one estimated by the response value routine in ELEFAN excluding 24-25 length-classes (Table 1).

Recruitment Patterns
The species recruitment patterns were estimated, allowing for the reconstruction of recruitment pulses from a time series of length-frequency data; determine the number of pulses per year and the relative strength of each pulse.The recruitment pattern was estimated based on the back projection of the restructured length frequency data estimated by the VBF [46].There were two recruitment pulses in data (such as preliminary observed in histogram plots).Therefore, the composite data was decomposed using the approach implemented in NORMSEP.In this case, starting estimates of the mean (in months) were automatically estimated above by the ELEFAN routine.The average yearly percentage of adults and juveniles were compared by the t-test paired sample test after an arcsine of square root transformation, as advised for percentage data [47].

Mortality Parameters Natural Mortality
The natural instantaneous mortality rate (M) was estimated based on the empirical model suggested by [48], incorporated in the FISAT II tool [43,49] as the following: where: L ∞ and K are the parameters obtained from ELEFAN I, and T ( • C) is the mean yearly surface water temperature in the study area (average SST = 21 • C, www.fishbase.orgstudies;[11]).This model takes into account the instantaneous mortality rate, and is very suitable model for small pelagic fish [27,43].

Total Mortality
The instantaneous total mortality rate (Z) was estimated using the length-converted catch curve method, and the respective analysis was based on the middle length of the individual in the catch, as shown by the following expression: N i is the initial number of fish in length class i, and N t is the number at time t (t is the age corresponding to the mid-length class i).Z with seasonality was then computed from the regression equation.The instantaneous fishing mortality (F) coefficient was estimated as: Length-converted catch curves are built around assumptions that are similar to those involved in age-structured catch curves, and they require as an input the growth parameter estimations obtained from the previous routines.Therefore, the input parameters required were L ∞ , K, and t o , and also the identification of the smallest length class that is recruited (L ).The best model can be selected based on statistical values of length-converted catch curves regression analyses.By selecting or deselecting the first and last data observed points, regression analyses are computed automatically, and based on the goodness of fit, the best model can be selected.The extrapolated or selected points will be used to approximate the probability of capture of length-converted catch curves (for details see: [48]).

Probabilities of Capture
The probability of the individuals being captured from the moment they reach the reference cohort length was estimated using the length-converted catch curve, in order to assess the degree of vulnerability of the recruits to fishing, and to the exploitation process in general.On the other hand, since this species shows a fast growth, we used the same length frequency data to identify possible recruitment peaks throughout the year.The length at first catch (Lc) was determined from the logistic curve adjusted to the cumulative probability of capture at each size class.This can be applied also to identify the fishing regime [27,50].For T. vitrirostris, Lc was estimated for cumulative probabilities of capture of 25, 50, and 75% (Lc 25 , Lc 50 and Lc 75 ).The size values for 25, 50, and 75% of capture (Lc 25 , Lc 50 , and Lc 75 ) were estimated from the logistic curve adjusted to the cumulative probability of capture of T. vitrirostris [43].The estimation of relative yield and relative biomass per recruit allows for the prediction of changes in Lc, and in the rate of exploitation (Ec) [27,49].
The stock exploitation rate (E) of T. vitrirostris, was analyzed through relative yield and relative biomass per recruit using the [32], as shown by Equations ( 7) and ( 8), respectively: where: U = 1 − (Lc/L ∞ )-the growth to be completed after entry into the exploitation phase; m = (1 − E)/(M/K) = (K/Z), and E = F/Z, the exploitation rate, i.e., the fraction of mortality of the T. vitrirostris caused by the fishery, F-the instantaneous fishing mortality coefficient, and Lc is the length at first capture.The relative biomass per recruit (B'/R) was estimated as follows: This model allows an estimate of: (i) the optimum yield (E 10 ), the exploitation rate at which the marginal increase of the relative yield-per-recruit is 1/10th of its value at E = ¨0; (ii) the optimum sustainable yield (E 50 ), or the yield for the reduction of 50% of the stock or value of E under which the stock has been reduced to 50% of its unexploited biomass, and (iii) the maximum sustainable yield (MSY or E max. ), which is the exploitation rate that produces the maximum yield.

Yield Isopleths
Population parameters-growth, mortality, age at recruitment, can be used to explain and understand the history of the fishery, and to predict its future, namely, to predict the effect of changes in the pattern of fishing due to any proposed conservation or regulatory measure.In stock assessment, fish availability can be estimated in terms of the yield of a single year-class of fish throughout its life, which will be (considering a steady state situation) the same as the yield in one year from all year-classes present in the fishery.The life span of a species can be divided into relatively short periods, and therefore the number of fish alive, the number caught, the number of fish dying due to natural) causes, and the number or yield in the weight of fish surviving to the beginning of the next period can be calculated.The yield contours were plotted to identify the impact on the yield of changes in the maximum sustainable yield (E max. ) and the critical length ratio Lc/L ∞ [42].

Table 1 .
Summary of the stock assessment parameters for all models tested: L ∞ = Asymptotic length (cm); K = Growth rate; Rn = Goodness of fit; Z/K = The ratio of total mortality to growth rate; Z = Total mortality; M = Natural mortality; F = Fishing mortality (cm); E est.= Estimated exploitation rate; E max.= Maximum Sustainable Exploitation rate; M/K = ratio of natural mortality to growth rate; Lc/L ∞ = ratio of middle-length to asymptotic length; r = coefficient of correlation.

Figure 1 .
Figure 1.Monthly von Bertalanffy growth curve of Thryssa vitrirostris superimposed over the normal length frequency histograms (A) and superimposed over the restructured length frequency (LF) histograms (B).LF (A) allows for the covering of the life span of several cohorts, and reconstructs an average cohort, while reconstructed LF (B) allows for study of temporal fluctuations of the recruitment.

Figure 1 .
Figure 1.Monthly von Bertalanffy growth curve of Thryssa vitrirostris superimposed over the normal length frequency histograms (A) and superimposed over the restructured length frequency (LF) histograms (B).LF (A) allows for the covering of the life span of several cohorts, and reconstructs an average cohort, while reconstructed LF (B) allows for study of temporal fluctuations of the recruitment.

Figure 3 .
Figure 3. Monthly recruitment pattern of Thryssa vitrirostris by the year (A), monthly percentage of 4 cm of length class registered along the years (B), monthly percentage of 5 cm of length class registered along the years (C).

Figure 3 .
Figure 3. Monthly recruitment pattern of Thryssa vitrirostris by the year (A), monthly percentage of 4 cm of length class registered along the years (B), monthly percentage of 5 cm of length class registered along the years (C).

Figure 3 .
Figure 3. Monthly recruitment pattern of Thryssa vitrirostris by the year (A), monthly percentage of 4 cm of length class registered along the years (B), monthly percentage of 5 cm of length class registered along the years (C).

Figure 4 .
Figure 4. Length converted catch curves applied to length frequencies data of Thryssa vitrirostris.The slope of the right descending arm (black dots) of the curve allow the estimation of total mortality (Z).White dots (expected values) are the expected number of fish, and grey dots (empirical data) the numbers of fish actually sampled (observed).

Figure 4 .
Figure 4. Length converted catch curves applied to length frequencies data of Thryssa vitrirostris.The slope of the right descending arm (black dots) of the curve allow the estimation of total mortality (Z).White dots (expected values) are the expected number of fish, and grey dots (empirical data) the numbers of fish actually sampled (observed).

Figure 4 .
Figure 4. Length converted catch curves applied to length frequencies data of Thryssa vitrirostris.The slope of the right descending arm (black dots) of the curve allow the estimation of total mortality (Z).White dots (expected values) are the expected number of fish, and grey dots (empirical data) the numbers of fish actually sampled (observed).

Figure 5 .
Figure 5. Probability of capture to each size class of Thryssa vitrirostris obtained from the length-converted probability catch curve (A) and percentage of juveniles (length-class < 13 cm) and adults (B) for the six year period (2009-2014).

Figure 5 .
Figure 5. Probability of capture to each size class of Thryssa vitrirostris obtained from the length-converted probability catch curve (A) and percentage of juveniles (length-class < 13 cm) and adults (B) for the six year period (2009-2014).

J 16 Figure 5 .
Figure 5. Probability of capture to each size class of Thryssa vitrirostris obtained from the length-converted probability catch curve (A) and percentage of juveniles (length-class < 13 cm) and adults (B) for the six year period (2009-2014).

Figure 6 .
Figure 6.Beverton and Holt's relative yield per recruit and average biomass per recruit models for Thryssa vitrirostris (A): E 10 -green line; E 50 -optimum sustainable yield, redline, and the (E max.-maximum sustainable yield, yellow line.Yield isopleths plot (B) with exploitation rate variability across the critical length ratio, Lc/L ∞ .

J 16 Figure 7 .
Figure 7. Location of study area (Adapted from Google Maps).

Figure 7 .
Figure 7. Location of study area (Adapted from Google Maps).

Table 1 .
Summary of the stock assessment parameters for all models tested: L∞ = Asymptotic length (cm); K = Growth rate; Rn = Goodness of fit; Z/K = The ratio of total mortality to growth rate; Z = Total mortality; M = Natural mortality; F = Fishing mortality (cm); Eest.= Estimated exploitation rate; Emax.= Maximum Sustainable Exploitation rate; M/K = ratio of natural mortality to growth rate; Lc/L∞ = ratio of middle-length to asymptotic length; r = coefficient of correlation.