The shrimp trawling fishery is the most important one in Mexico in value terms and given its putative environmental, societal and economical implications, it is also the most difficult to manage. Although this fishery was restricted from national bays and estuaries since the 1970’s, local fisheries cooperatives recently claimed access to shrimp stocks within La Paz Bay by using an artisanal fleet and a low impact trawling net. This study is aiming at simulating some ecosystem level effects of such a potential fishing effort release. We explored the response of three ecosystem indicators under two different exploitation scenarios: 30% and 80% of shrimp biomass removal. The indicators were: relative ecosystem biomass distribution as function of trophic level, trophic replacement and interaction strength, all computed from the outputs of a mass balance dynamic model (Ecopath with Ecosim) of this ecosystem. Our results suggest that moderate fishing scenario (30%) would not cause major changes in either indicator whilst the scenario of strong fishing pressure (80%) seems to increase not only the fish resources variability at the population level but also the variability of the overall biomass, hence potentially reducing ecosystem stability.
In 2006, the economic income of the shrimp fishery in the Gulf of California was in excess of $411 million, ranking the first place as for the total value of marine resources. In the Mexican Pacific, roughly 1500 vessels target marine shrimp stocks; more than 70% of this fleet operates in Sonora, Sinaloa, Baja California and Baja California Sur [
Although the shrimp fishery is the most lucrative one in Mexico, it is also the most difficult to manage. The fact that this activity has been a cornerstone for the regional development during the last 50 years, counteracts with the environmental effects attributable to the fleet operation (habitat loss, fishing mortality of non-target animal and plant species), the recent decreasing trend of the catch per unit of effort, overcapitalization of the industry, federal subsidies and the concomitant societal and economic impacts. In light of this situation the academic, governmental and social sectors are deeply concerned over the long term ecological sustainability of the Mexican shrimp fishery. As a part of the government’s response to this problem, since 1993 a portion of Sonora and Sinaloa’s shrimp trawling fleet that had operated on sandy bottoms in La Paz Bay, was officially banned [
However in 2005, local fishing cooperatives claimed access to shrimp stocks within the bay (
This study is aiming to address some potential ecosystem level effects if fishing effort of artisanal shrimp trawling was unleashed in La Paz Bay, by applying three ecological indicators: relative ecosystem biomass distribution
as function of trophic level (RBD), trophic replacement (TR) and interaction strength (IS) computed from the outputs of a mass balance dynamic model Ecopath with Eocsim [
The present study is based on an Ecopath mass balance model of La Paz Bay [
Time simulations were made by using Ecopath with Ecosim (EwE) software, which provides a tropho-dynamic
simulation with ecosystem level considerations. EwE, through a system of differential equations, estimates energy flow rates among groups as a function of harvest rates and biomass changes through time for all functional groups in the model [6-8]. The instantaneous (Ecopath) mass balance model is defined by the following Equation:
where Bi and Bj are biomasses (the latter pertaining to j, the consumers of i), P/B their production/biomass ratio, equivalent to total mortality rate, EEi the fraction of production that is consumed within, or caught from the system, Yi is the fisheries catch (Y = FB; F is fishing mortality rate), Q/Bj the food consumption per unit biomass of j, and DCji the contribution of prey i to the diet of predator j.
The dynamic form of the model (Ecosim) is defined as:
where dB/dt is rate of change in biomass, g is growth efficiency, F is fishing mortality rate, M is natural mortality rate (excluding predation), e is emigration rate, I is immigration rate, and Cij (Cji) is the consumption rate of type i (j) biomass by type j (i) organisms.
We explored two different fishing scenarios affecting the group denominated as “brown shrimp” (Farfantepeneaus californiensis), and their ecosystem level effects through three indicators: relative ecosystem biomass distribution by trophic level (RBD), trophic replacement (TR) [
ER is proportional to the total stock biomass and can have values from 0 to 1, being 1 the exploitation rate equivalent to extracting the whole population. For each ER value (at each simulation year) we iteratively estimated the correspondent fishing mortality. The first scenario was one of moderate fishing exerted on shrimp stock, beginning in year 20 such as at the end of the simulation ER = 0.30 (i.e. 30% of the virgin biomass is removed). The correspondent F values ranged from F = 0 up to F = 5. In the second scenario, one of stronger fishing intensity, the exploitation began at year cero and ended with an ER = 0.80. In this case F values ranged from F = 0 up to F = 15. For both scenarios the simulation time-frame was 70 years.
Trophic replacement TR indicator quantifies the trophic replacement of a species that is removed from a system or drastically reduced in biomass (e.g. collapsed through overfishing) by other species or group in the system. It quantifies the extent to which trophic niche replacement occurs following a collapse of, or increases in, a target stock. The TR indicator is calculated for each group j (and year) belonging to the set of groups (K) that show a change in biomass of the opposite sign to that of the target group i (e.g. if target group i collapses, K is a group undergoing an increase following the collapse of group i:
where ΔB is the change in biomass of any group between consecutive years. The indicator ranges between 0 (no replacement) and 1 (total replacement).
IS indicator assumes that a change in biomass of a strong interactor causes large changes in the biomass of other groups. The relative change in the biomass of a species or functional group is expressed as a proportion of the sum of the relative changes of all groups in the system:
where B is biomass in year 0, i is the group being tested, j is another functional group in the ecosystem, ΔB is the change in biomass of a group and n is the total number of groups in the model. Values for the indicator lie between 1 and 0, with large values indicating that species (or group) i has a stronger impact on other groups than species with smaller values. IS and TR were derived from the trophic model outputs for each group and year of simulation.
After the simulation, biomass data were standardized considering the highest value of each individual group at each year as the unit (1). In this way, differences in magnitude of biomass among groups become proportional (RBD). Groups in all cases were sorted by trophic level. Resulting matrices (groups/species sorted by trophic level in the ordinate, simulation years in the abscissa and indices values in the vertical axis) were plotted as surface diagrams.
Results are shown in
The second scenario, as it involves extracting 80% of the shrimp population biomass, notably modifies the behavior of all indicators over time, particularly RBD. Analyzing the RBD in different time windows, it is evident that the degree of variation is constantly increasing over time until a similar pattern to that of the last stage of the biomass distribution in the first scenario is reached, though exacerbated and of much less duration (10 years); i.e. small (high) values in scenario 1 are smaller (higher) in scenario 2. The IS index shows comparable changes between the two scenarios. In the first one, the TR level is considerable high in all functional groups, and in the second one it displays higher variation but of lesser intensity.
Previous population and ecosystem dynamic models applied to the shrimp fishery in La Paz Bay indicate that if trawling activity operates without any depth restriction, then extracting 38% of the whole population represents a target reference point. At this level of extraction, ecosystem impacts would not be significant in terms of the trophic structure and the mean trophic level of the catch. If the fishery is restricted to depths in excess of 40 m, then the fraction of the population that can be extracted increases to 50% without significantly affecting those ecosystem indices (F. Arreguín-Sánchez, unpublished data).
These results are compatible with ours. As the first scenario indicates, a moderate fishing pressure causes slight changes in ecosystem biomass and biomass stability, even after removing 30% of shrimp biomass. Likewise, it was found that trawling fishery in La Paz Bay positively impacts shrimp populations [
Other indicators show that brown shrimp in La Paz Bay naturally have low TR and IS values. Once fishing pressure starts affecting shrimp populations, both indicators increase (Figures 2(c) and (d)). This suggests that fishing pressure has small effects on either indicator; perhaps shrimps and groupers display the more noticeable changes (Figures 2(c), (e) and (f)). In the case of shrimps, being trophic generalists they are more prone to be replaced by similar species. Furthermore, faunal assemblages in La Paz Bay may be sufficiently diverse as to allow a dynamic functional replacement of fish populations decimated by exploitation [12,13].
On the other hand, IS evolution suggests that brown shrimp, being a foraging species, has an important trophic role in the ecosystem, tightening the interaction with its many predators as it is increasingly fished down (
by trawling gears in La Paz Bay, representing up to 90% of total catch, and are highly predated.
Moderate fishing scenario would not cause mayor changes in RBD, TR nor IS. The more noticeable modification is the biomass reversion across all trophic levels, although, without risking the ecosystem stability. Conversely, a scenario of strong fishing pressure seems to increase fish resources variability at the population level and, as suggested [
Actual harvest rate of the trawling fishery on shrimp populations within La Paz Bay is negligible and even if it started at moderately exploitation levels (F £ 5) would cause no significant change in other fish stocks (i.e. ecologically sustainable). However, a heavy fishing regime has the potential of disrupting ecosystem stability and resilience. More research is needed in order to quantitatively delimitate critical ecosystem threshold levels between these two conditions.
LASG thanks PIFI-IPN Program and CONACyT schoolarships. Authors thank partial support received through projects SIP-IPN (20121172, 20121444) and SEPCONACyT (104974, as well as COFAA and EDI programs of the Instituto Politécnico Nacional (IPN).