Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

More than Anecdotes: Fishers’ Ecological Knowledge Can Fill Gaps for Ecosystem Modeling

Abstract

Background

Ecosystem modeling applied to fisheries remains hampered by a lack of local information. Fishers’ knowledge could fill this gap, improving participation in and the management of fisheries.

Methodology

The same fishing area was modeled using two approaches: based on fishers’ knowledge and based on scientific information. For the former, the data was collected by interviews through the Delphi methodology, and for the latter, the data was gathered from the literature. Agreement between the attributes generated by the fishers’ knowledge model and scientific model is discussed and explored, aiming to improve data availability, the ecosystem model, and fisheries management.

Principal Findings

The ecosystem attributes produced from the fishers’ knowledge model were consistent with the ecosystem attributes produced by the scientific model, and elaborated using only the scientific data from literature.

Conclusions/Significance

This study provides evidence that fishers’ knowledge may suitably complement scientific data, and may improve the modeling tools for the research and management of fisheries.

Introduction

Ecological models have been advocated by the Food and Agriculture Organization of the United Nations (FAO) as promising methodologies for ecosystem-based assessments aimed at managing fishing resources [1,2]. Ecosystem-based fisheries management (EBFM) focuses on a multi-species approach as part of an integrated vision of the ecosystem, rather than on single target populations or fishing fleets [3,4]. Therefore, it should include ecological, social and economic aspects [5], and simultaneously, consider fish, fishers, the maintenance of fishery resources, and the environment [6]. Ecosystem modeling may provide the analytical tool for achieving these goals by considering the ecosystem as a whole, and underpinning the direct and indirect effects of a specific fishery on target and non-target species, and on other fisheries [7].

Ecosystem modeling in fisheries is often hampered by a shortage of consistent and regular fishing statistics [8]. At the same time, the models are also susceptible to the lack of local and quantitative ecological data [9]. Fisheries, on the other hand, are often data limited, markedly so in developing countries, which have no current mandatory fishery statistics [10]. As a result, roughly 99% of fisheries or fish stocks worldwide are unassessed, and remain data-poor fishery scenarios. Beyond that, many of the most data-poor fishing grounds are also the world’s most vulnerable and rapidly changing ecosystems [11].

In contrast, fishers can provide broad and detailed information on fish species and ecology [12], habitats [13], and the environment as a whole [14], including its past and recent changes [15,16,17]. Fishers also retain important information on fishing activity, such as species composition, catches, and abundance [11,18]. This critical data availability with regard to fisheries has led to increasing interest in the place-based knowledge of resource users [19,20], increasing the use of fishers’ ecological knowledge (FEK) [11], or simply FK, as we refer to here.

FK is generated by a fisher’s observation of their fishing environment [21], and by interactions with and observations of other fishers [22]. Its use has been increasingly recognized as a good source of data [12], which has a low cost and is easily sampled by simple methods [11]. However, a fisher’s knowledge must be used with caution, since fishers, like everybody else, may miss or forget information [12]. Moreover, the costs of data collection vary under different circumstances, and data reliability should be tested whenever possible [10]. Aiming to appropriately collect, analyze, and use fishers’ knowledge, researches have recently focused on methodological issues for using FK, defining a suitable sample size, and using the proper variance in the samples, underscoring the awareness of the suitability of the data registered [11].

For modeling and management purposes, the use of FK associated with science could enhance knowledge and user participation, by sourcing the most experienced and, therefore, expert fishers in the fishery and the environment [15,22]. Notwithstanding the evidence on the validity of FK, and its growing employment to fill the gaps in the scientific knowledge [12,13,15,16,18,23], fishers’ knowledge is sometimes referred to as “anecdotal knowledge”, likely limiting its translation into scientific knowledge, the development of methods to properly sample this knowledge, and its application to management [21]. As a result, the inclusion of FK into fisheries management remains the exception, rather than the rule [14]. Following the slow pace for new ideas, FK has thus far not been included in data sampling aiming to model fisheries, despite the scientific recommendations [24].

As a first step to using FK in the ecosystem modeling approach, we interviewed fishers to estimate the parameter values for a model of a coastal fishing area. We used one of the most widely applied modeling approaches and software, Ecopath with Ecosim (EwE) [8], to analyze their contribution to the modeling process, while describing the marine fishing environment, its food web, and the fisheries. The main questions addressed were: Could an Ecopath model be elaborated using fishers’ knowledge? If so, which variables are fishers more capable of informing? Do the attributes generated in the FK model match those of a scientific-based model? Is it possible to improve a food web model sourcing fishers as experts to complement the information? The main expectation here is to compare the ecological role of the relevant species between the two models, investigate the ecosystem attributes, identify the lack of information that could be provided properly by fishers, and determine what should be focused on by the scientific community. We hope to underscore the role of fishers and science as partners to inform fisheries, enhancing managers’ understanding of the fisheries-related food web.

This paper is structured as follows: The FK and science-based (SC) information is discussed within the context of fishing ecosystem modeling. Then, the FK data collection is described, focusing on the methods for selecting interviewees (our experts), and for consulting them. Similarly, the approach used for collecting the science-based data is presented. Next, the methodology used to build the Ecopath model and the output analyses involved are summarized. The attributes of the FK model and SC model are compared, and the paper is concluded with a view on the potential of fishers providing information for modeling and the role of science in filling the knowledge gaps. The aim was to show the value of this new approach to model fishing ecosystems and to better contribute to the understanding, studies, and future management of fisheries.

Study Area

The study area included the small-scale fisheries at the continental shelf in the Rio Grande do Norte state, on the Brazilian northeast coast, corresponding to the FAO area “Northern-Natal division 41.1.2” (Fig 1). In this region, the continental shelf break occurs approximately 30–40 km from the coast, usually at a depth of 60–80 m [25]. The small-scale fisheries occur over an area of approximately 16,400 km2, with substrates characterized by reef coral sediments and calcareous algae (6 km from coast, and 20 m in depth), showing a high diversity and low abundance of fishes and invertebrates.

thumbnail
Fig 1. Study area of the fishery community (Baía Formosa) on the Brazilian northeast coast.

https://doi.org/10.1371/journal.pone.0155655.g001

Small-scale fisheries in Rio Grande do Norte catch around 11,500 tons per year, which corresponds to 80% of the total catch of this state [26]. The fishery activities are distributed over 95 coastal communities in Rio Grande do Norte, and directly provide employment for 12,300 fishers [25]. Hand lines and gillnets are the most commonly used gear, but the fishers also use traps and diving, as well as trawlers and diving for invertebrates. The fishing fleets operating in this state include sail or paddle powered (about 3,750, or 80% of the total number of boats), motorized (about 900, or 20% of the total), and industrial (about 50, or 1.5%) [27].

The major species in the official landing statistics include tuna and tuna-like species (Thunnus atlanticus and Euthynnus alletteratus), mackerel (Scomberomorus brasiliensis and S. cavalla), grouper (Mycteroperca bonaci), and snapper (Lutjanus analis and L. jocu). Invertebrate catches are also economically important, and include octopus (Octopus insulares), lobster (Panulirus argus and P. laevicauda), and shrimp (Litopenaeus schmitti, Xiphopenaeus kroyeri, and Penaeus brasiliensis) [27].

Material and Methods

Modeling and data used. Two ecosystem models were constructed for the same fishing area using the Ecopath with Ecosim software: one model based on the fishers’ knowledge and one model based on the available scientific information. For the first model, the data were collected by interviews, and for the second, the data were gathered from the scientific literature.

In both models, we used the same landing data provided by the Brazilian Environmental Agency (IBAMA). From 2001 to 2011, the artisanal catches for each species or functional group were selected. The catch data included in the model was restricted to the landings, and does not include the estimate of the discards or catches from other fishing activities (recreational, subsistence, research, or illegal). The landing data were expressed in tons (t) per square kilometer per year.

Due to the lack of information on the fishing area of focus, the biomass was predicted based on catches using the Vasconcellos and Cochrane model [28] for both Ecopath with Ecosystem (EwE) models. This model predicts the biomass dynamics according to the Schaeffer model: where B is the biomass, r is the intrinsic rate of the increase in biomass and K is the virgin biomass or stock at the arraying capacity [28]. The K parameters were obtained from fishbase.org, and the catch (Ct-1) data are described above. Finally, the analyses were carried out using R software, specifically, the COM-SIR package.

Fishers’ Knowledge Model (FK Model)

Selecting interviewees.

All of the interviewed fishers live in the fishing community of Baía Formosa, the most productive fishing community in the austral portion of the state. This community has roughly 200 fishers, and all of them land their catches at Beach Porto, where the landings were sampled.

Only those fishers with higher catches were approached to provide (mainly) the information unavailable in the scientific literature. The interviewee choice was guided by the landing data previously registered from February of 2013 through February of 2014. Based on the landing data registered by the research team during this one-year period, we identified the ten most landed fish species. The recorded landings further enabled us to identity the fishers with higher catches, when using one of the two main types of gear (hand lines or gillnets). For the purpose of this study, these were considered to be the expert fishers that should be consulted.

In addition to having higher fishing productivity, the expert fishers also needed to fulfill pre-established criteria to be considered in this survey: (1) willingness and availability to participate in the research, (2) fishing experience (in years), and (3) being an elder, assuming that older fishers had fished for a longer period of time and would have greater ecological knowledge [29].

Method for consulting expert fishers.

After the aforementioned selection of experts, 23 fishers were individually selected and interviewed face-to-face in the Baía Formosa community, between February and April of 2014. The same researcher conducted all of the interviews, and the expert fishers were asked for ecological information about the ten most often caught fish species.

The interviews were conducted following an adaptation of the Delphi methodology [30], seeking to add value to the local knowledge by achieving consensus among experts over several rounds of investigation. This method is driven by the assumption that combining the expertise of several individuals will provide more accurate results than what can be achieved from consulting a single individual only [14]. The fishers may have been interviewed under many rounds of consultation, aiming to achieve consensus on the information provided by them. Then, the results of all of the responses in each round were summarized and presented individually to each participant, with the purpose of finding consensus among them. Consensus occurred when at least 50% of the experts reached the same conclusion and agreed on the results. Then, the data were based on individual knowledge and consensus, avoiding bias due to the presence of a group (known as social desirability) or groupthink [30].

Here, were used two rounds of interviews (Fig 2): in the first one, the individual expert fishers could choose one fish species to address. After that, other fish species were presented (in color photographs) and the fisher could discuss them. At this step, the interviewer showed the information from the first round (from other fishers), and each fisher could review his responses according to the new general data. As stated before, the information was considered consensual when confirmed by more than 50% of the fishers.

thumbnail
Fig 2. Interviewing methodology: Steps for interviewing methodology applied to register fisher knowledge aiming to select expert fishers and to fill the gaps in the scientific data needed for modeling.

https://doi.org/10.1371/journal.pone.0155655.g002

Considering the fact that the fishers could not follow the goals of the scientists with regard to the research, the research team conducted two previous meetings with the fishers in the middle of 2013, which had theoretical and practical components. During the first meeting, the fishers were presented with scientific studies on fisheries ranging from unispecific to ecosystem-based evaluations. A few months later, they took part in a participatory approach to learning and evaluating a fishery resource, stock assessment, and management. Moreover, the research team has sampled landing data with the collaboration of local fishers since January of 2013. This scenario has enabled the interviews with the fishers, increasing their involvement with and understanding of the scale and scope of this study.

All of the fishers approached were first informed about the research purpose; then, they were asked about their willingness to participate in this study. Only those who gave verbal consent were included, and their participation was always witnessed by a third person from the university, who was not involved in the research. Despite sounding unusual, oral consent has been accepted and approved by ethics committees in previous research [16]. Oral consent is an instrument that allows the inclusion of illiterate interviewees and/or traditional communities who are not comfortable with tape recording. Actually, there was a high proportion of illiteracy among the interviewed fishers, and to register their consent to participate, the fishes were asked to mark a specific place on the consent form or write their name, if they were formerly trained in this for other purposes.

The data sampling was approved according to the guidelines of the Committee of Ethics at the Federal University of Rio Grande do Norte (Protocol of approval No. 48378415.9.0000.5537). The rules of approval state that the data and metadata related to personal interviews may be available under a specific request and justification, aiming to protect interviewees who may have provided private information.

Registering Fishers’ Knowledge and estimating parameters.

During the above-mentioned consultation of fishers, the information registered included population parameters like the (1) maximum weight (Wmax) and/or length (Lmax), (2) modal weight (Wmodal) or length (Lmodal), (3) growth rate, and (4) longevity, as well as trophic information like the (5) diet composition, (6) amount of food required per day, (7) main predators, (8) total catch per species during the last year, and (9) total stock size. The details about the data collection instruments are available in the Supporting Information (S1 Table).

The expert fishers were presented with pictures of food items in order to assist them when providing information about their diet items. Even when using the pictures, the interviewees added other species as diet items, providing detailed information about their diet compositions (S3 Table). In addition, the expert fishers could use pieces of graduated tape to show their perceptions of the estimation of length (Lmax and Lmodal in cm). In the cases in which the fishers provided the fish dimensions in weight (kilograms), the length (L∞) was estimated by weight-length relationships using an online database (fishbase.org). The fishers provided information on the Lmax, which was used to estimate the asymptotic length (L) for each species, based on the empirical equation of Froese and Binohlan [31].

The fishers’ knowledge data provided information to estimate the Natural Mortality through Pauly’s equation [31], which was assumed to be equal to the production/biomass (P/B). The weight estimates based on the fishers’ knowledge were used to calculate the consumption rate (Q/B), using the Palomares and Pauly equation [32].

During the interviews, in addition to being asked about the stock sizes for the ten most caught fish species, the fishers were provided with corn seeds that they could use to represent the quantity of each species they caught annually. They were encouraged to add as many seeds as they needed to coarsely represent the total size of each fish stock. However, the fishers were unable to provide accurate information on the annual catch or the stock sizes, by any means.

Science data-based Model (SC Model)

Aiming to provide information for all of the compartments in the scientific model, data was gathered from scientific papers, reports, theses, governmental reports, online databases, and other EwE models created in Brazil or in other tropical ecosystems [3336].

As in the FK model, the P/B for the fish species was calculated based on Pauly’s equation [37], and the Q/B was estimated using the method of Palomares and Pauly [32]. For the non-fish compartments, we used Gascuel et al. [38] for the P/B estimates, and multiplied them by three for the Q/B. All of the basic input parameters were expressed per unit of surface area, that is, in tons per square kilometer. The basic input estimates for each species were used to calculate an average value of that basic input for each functional group.

The diet matrix was summarized as the percentage of each group in terms of the total wet weight (or volume) in the diet of the predator. The diet composition matrix was determined using an online database (fishbase.org) and other studies of the same group or species and region. The main studies consisted of dolphins [3941], sharks [4245], dogfish [46], large pelagic fish [43,47], Thunnus atlanticus [48], grouper [49], Lutjanus spp. [50], Coryphaena hippurus [43], Euthynnus alletteratus [51], octopus and squid [5254], carnivorous zoobenthos [5557], detritivorous zoobenthos [57,58], and Scomberomorus brasiliensis, Seriola fasciata, Cynoscion jamaicensis, and S. cavalla (fishbase.org) (S3 Table).

Model Description

EwE version 6.4.3 was used to elaborate both models (FK and SC). This modeling tool was initially developed by Polovina [59], and has been further developed by Christensen and Walters [60]. EwE enables the integration of a large body of information into a coherent description of an aquatic food web, while accounting for human activities in an ecosystem context, and environmental changes [8].

The Ecopath model is built on a system of linear equations to describe the average flows of mass and energy between a series of functional groups that represent the organisms in a food web. The flow to and from each functional group is described by the main Ecopath equation representing the production of each group i: where Bi is the biomass of i, P/Bi is the production/biomass ratio, EEi is the ecotrophic efficiency, Bj is the biomass of the consumers or predators j, Q/Bj is the consumption per unit of biomass of j, DCji is the fraction of the prey i in the diet of predator j, Yi is the total fishery catch rate, Ei is the net migration rate, and BAi is the biomass accumulation rate [60].

In order to parameterize the model, three out of four basic parameters must be provided for each functional group: biomass (B), production per unit of biomass (P/B), consumption per unit of biomass (Q/B), and ecotrophic efficiency (EE). The algorithm then estimates the fourth parameter to ensure the mass balance. For each group, only one of the basic input data could be missing, which was then estimated by the model [61]. The diet composition and landings are also required.

Output Analysis.

For this analysis, we estimated the ranking correlation between the 10 main species in the landing data (in tons) and other catch species according to the fishers’ interviews, using the non-parametric Spearman index. The paired t test was used to compare the input parameters, P/B, and Q/B between the models, assuming that each interview could be treated as an independent sample.

A paired t test was also used to compare the trophic level (TL), omnivory index (OI), and keystone index (KS) to each of the main fish species in both EwE models. The TL provides an estimate of the trophic position of a particular species or functional group in the food web. The OI measures the distribution of the feeding interactions among the trophic levels by functional group [61], while the KI estimates the effects of each group on all of the other groups in the food web, including the indirect ones [62].

The attributes used to compare the outputs in both models (FK and SC) were the total production, total throughput system (total flow, including detritus), and total net primary production. The general attributes of the ratios used to evaluate the matches between the models were the total primary production/total respiration, total primary production/total biomass, and total biomass/total throughput. The total biomass (excluding detritus) and total catch were used to compare the biomasses. We also compared the mean trophic levels of the catch. The main indexes to measure the resilience and stability in each model were the connectance, system omnivory, ascendancy, overhead, Finn cycling index, and Finn mean path length [60].

Results

Between February and April of 2014, 23 fishers were interviewed (12 hand line users and 11 gillnet users). Their ages ranged from 26 to 57 years old (median = 41 years old and mode statistic = 35 years old), and their mean age when starting in the fishery was 14 years old. Their fishing experiences varied from 10 to 45 years, with an average of 24 years. Altogether, they represented 470 years of fishing experience.

The most commonly caught species according to the interviews were Lutjanus analis (mutton snapper, 17.1%), Scomberomorus brasiliensis (Spanish mackerel, 14.6%), and Thunnus atlanticus (blackfin tuna, 12.2%). Euthynnus alletteratus (Atlantic little tunny) and Rhizoprionodon lalandii (dogfish) were mentioned by 9.8% of the fishers. Each one of the following species Coryphaena hippurus (dolphinfish), Cynoscion jamaicensis (weakfish), Mycteroperca bonaci (black grouper), Scomberomorus cavalla (king mackerel), and Seriola fasciata (lesser amberjack) were referred to by 7.3% of the fishers. Moreover, the fishers reported all of the aforementioned species as the ten most caught fish species.

The correlation between the catches and frequency in which the species were chosen by each fisher during the interviews (Spearman index ρ = 0.83) shows how confident the fishers were when talking about the species with which they were most familiar and had greater knowledge, that is, species with higher landings.

Basic and Output Parameters

A total of 25 compartments represented the ecosystem in both models (FK and SC). Each compartment was expressed by either a single species or functional groups, and was comprised of ecologically related species, depending on their importance in the landing records (Table 1, Fig 3).

thumbnail
Fig 3. Modeling of trophic network by Ecopath with Ecosim for the fishing habitat of Baía Formosa (Brazil).

(A) Fisher knowledge model and (B) science-data model as input parameters. Black square: keystone groups.

https://doi.org/10.1371/journal.pone.0155655.g003

thumbnail
Table 1. Basic input for the EwE to the FK model and SC model.

B = biomass (wet weight–t*km-2); P/B = production/biomass (year-1); Q/B = consumption/biomass (year -1); EE = ecotrophic efficiency. Species aggregation in each compartment. FK model = fishers’ knowledge model; SC model = science-data model.

https://doi.org/10.1371/journal.pone.0155655.t001

The expert fishers were unable to answer questions about longevity, total catch, and stock size, or about the amount of food required per day, which they overestimated by more than 100 times that expected by the fish metabolism (Table 2). For example, according to the fishers, the piscivorous species would eat 3 times their own weight in a day. Therefore, this, last information was not used for the modeling, and the Wmax allowed us to calculate the W used to estimate the Q/B for the FK model [32].

thumbnail
Table 2. Average values of fishers’ knowledge for the maximum (Wmax) and modal (Wmodal) sizes, and amount of food required per biomass per year.

The weight-length relationships were used to estimate the length (L∞) from the related maximum weight. The values in the brackets refer to the standard deviation.

https://doi.org/10.1371/journal.pone.0155655.t002

The growth rates were only obtained from a few of the fishers, and only for some of the pelagic surface species, such as Thunnus atlanticus, Scomberomorus brasiliensis, Euthynnus alletteratus, and Scomberomorus cavalla. The information for the total catch in the previous year and stock size were obtained only qualitatively (i.e. higher, the same, or lower). Nevertheless, the expert fishers provided information on the maximum (Wmax) and modal weights (Wmodal) (Table 2).

The FK and SC models did not differ in the P/B (t = 1.28, df = 23, p-value = 0.211) and Q/B (t = -1.92, df = 21, p-value = 0.068) ratios. Once the other data used to estimate the P/B and Q/B were the same for both models, they indicated the likeness between the values of the weights observed by the fishers and sampled scientifically, with regard to the ten most caught species.

For the fishers, the most mentioned predators were sharks (36 times), medium pelagics (23), great barracuda (20), dolphins (13), S. cavalla (11), and groupers (9). This information was used to build the diet matrix in the fishers’ model. However, the overall diet composition was more detailed in the matrix based on the scientific data (S2 Table). Nonetheless, such itemized information did not generate differences in either the trophic level (t test = -0.10242, df = 20, p-value = 0.9194) or omnivory index (t test = -1.759, df = 20, p-value = 0.1555) between the FK and SC models (Table 3). The keystone species indicated the first five species in the scientific model: carnivorous reef fishes, small pelagics, Cynoscion jamaicensis, detritivorous zoobenthos, and zooplankton. In the FK model, the most important species were carnivorous reef fish, small pelagics, large pelagics, and sharks (Table 3).

thumbnail
Table 3. Trophic level, omnivory index, and keystone index position obtained by the FK and SC models.

The values did not differ between the models for the trophic levels (p = 0.92) and omnivory index (p = 0.15).

https://doi.org/10.1371/journal.pone.0155655.t003

In general, the ecosystem attributes (Table 4) were very similar in each model. However, the system omnivory index, throughput cycled (excluding detritus), Finn’s cycling index, and overhead had higher values in the FK model. These differences are related to the different diet compositions used in each model.

thumbnail
Table 4. The parameters of the FK model and SC model to the fishing system modeled using EwE.

https://doi.org/10.1371/journal.pone.0155655.t004

Discussion

Globally, the information provided by the fishers has been used broadly for wide assessments, such the evaluation of fishery resources from the users’ point of view [63], as well as to track changes in fish diversity [15], complement available information on fishing and ecological changes [23], define management units [10], assess ecosystem services through spatial modeling [8], and depict trophic relationships among fishes [52].

With regard to modeling, few studies use fishers’ knowledge, mostly because this knowledge is considered to be anecdotal, non-methodological, and of limited application among scientists and managers [12]. As a result, studies using FK are mostly conceptual and qualitative, overlooking the potential of such data to enhance outcomes from ecosystem modeling under distinct views: from the scientists and from the fishers’ observations. One exception, thus far, is the quantitative modeling approach, which indicates the ability of fishers to report biodiversity losses and food web changes from a historical perspective [15].

Our study is the first to use fishers’ knowledge to parameterize an ecosystem model based on the Ecopath with Ecosim approach. Here, the 19 ecosystem attributes provided by the model which were elaborated using FK were equivalent to the ecosystem attributes produced by the model elaborated using only the data from the literature. Clearly, an EwE model may (partly) be successfully elaborated using the basic parameters provided by fishers. However, the fishers were unable to, even roughly, inform on the growth rates, longevity, biomass, and stock sizes for all of the species. These are precisely the main gaps to modelers in the literature, mainly in developing countries, given the lack of funding and policies to provide these data [8]. Additionally, the information on the food required per day was entirely inconsistent between the two sources of data (the fishers and the literature). When the fishers did not observe the fish during their entire life cycles, they did not understand the rate of growth of those fish. In addition, since they did not see their target fish species feeding, and did not have information about their predation pressure on the prey stocks, they were unlikely to inform on their requirements of food per day.

Conversely, despite the greater detail in the scientific data, the expert fishers were more precise in their information about the maximum and average weights (Wmax and Wmodal). These parameters tend to be more easily answered by the fishers, since they represent exceptional experiences [18] in our case, in terms of fish size (Wmax), and the expected typical size of each species (Wmodal). Specifically, the Wmax reminds them on uncommon and pleasant moments, which turn into memories, that are easily recalled.

The fishers were also precise in their information on diet composition, mainly for the ten most caught fish species. Indeed, the FK for all of the species has proven to be a good source of information about important commercial fishes’ diet compositions [13,64,65]. In some cases, however, this information may be based on which prey species they used as bait, and the researcher must be aware of this, and not overemphasize the bait items as food resources [11]. Most importantly, the fishers acquire their knowledge of fish diets gradually, through daily observations of species biodiversity and ecology, or the stomach contents while cleaning the fish that were caught [19].

In this regard, during the interviews, the fishers upheld slightly higher predation rates on small and medium pelagic species, than those reported in the literature. Once it is established that the group of small pelagics is comprised of local species used for the fishers’ consumption, and the group of medium pelagics contains species of interest to recreational fishing, this information is easier to grasp. Beyond that, the small pelagics emerge as keystone species in both the FK and SC models, possibly supporting the fishers’ observations.

With regard to the keystone groups, is interesting to note that the models elaborated upon the observations of the fishers, and the science underscored the compartments of the intermediate trophic level, which includes small pelagics, C. jamaicensis, and carnivorous reef fishes. What can be said here is, in general, fishers do not bias their information towards fishes of higher economic value, but by simple observation, they guided their modeling to the same outcome that the science would come up with. On the other hand, according to the FK model, the keystone groups would also include top-predators, while the science-based model indicated microinvertebrates. At the current poor data level, these opposed outcomes are conceivable, but challenging, if not impossible, to confirm.

The few discrepancies found here are noteworthy, and may have important management implications. They may indicate minutiae not observed by the fishers (such as the food required per day), and therefore, assumed to be unreal or unimportant. On the other hand, as immediate observers of the sea [23], the fishers may consider the science to be inaccurate if they are confident of specific local evidence (like the higher predation on small and medium pelagics), which is not given by the scientific data. In both of the aforementioned situations, the mismatch between the knowledge revealed by the fishers and the information attained by the science may prompt the fishers to exhibit non-compliance with further management decisions [16].

Currently, there is also a mismatch and low connection between the modeling efforts and their real-world applications, mostly because fishery scientists are the main ecosystem modelers, while the application relies on managers. However, the disconnection EwE modeling has indicated the most threatening fishing gear, the most vulnerable fishing resources or compartments, and the effects of fishing pressure and changes in the environment on the biomass of exploited species [23,66]. All of this evidence guides management suggestions, and the managers may or may not implement them. Regardless of the application of the modeling outcomes, modeling using actual temporal data allows for the proposal of management actions, and its validation under, for instance, changes in the ecosystem [67] or in the fishing pressure.

The proof of consistency among the estimates achieved by both models (FK and SC) is not a new finding, and the agreement between the local fishers’ knowledge and scientific data has been shown previously [12,16,18,23]. The novelty here relies on using fishers’ observations for quantitative modeling, and on revealing that the modeling outcomes were not enhanced or worsened by the fishers’ information. Under this evidence emerges the question of whether only the fishers’ knowledge could be used to model ecosystems in areas that have been poorly studied.

Prompted by the results achieved here, we contend that under dataless scenarios, the use of fishers’ knowledge should not be disregarded. However, simple methodological caution, like carefully selecting the informants when seeking specialized advice, may increase the reliability of the information provided by the fishers [11]. In this case, avoiding random interviews and seeking older fishers is likely to minimize the level of uncertainty [16,23]. This methodology is a practical tool to access FK and conduct ecosystem modeling, and it can be used as a basis for EBFM. One fundamental aspect of EBFM is that it is an interdisciplinary approach, where the fishers’ participation is essential for success [4]. In this sense, including the fishers’ views and experiences may enhance commitment from fishers when establishing management strategies [22].

The use of fishers’ knowledge as a tool to model fisheries is not common, and it is still in its early steps. Fisher experience-based knowledge has been combined with literature and data sampling to understand the loss of biodiversity and food web changes [15], and fishers’ knowledge about the trophic relationships of fishes has been compared with scientific data for providing new hypotheses [13]. In addition, the total catch modeling using the catch per unit of effort estimated by the fishers and by official data has been applied to data evaluation and support management [63].

Despite the consistency of information described by our results, and shown in many articles, fishers and scientists may fail. Therefore, divergences should be treated with similar skepticism, and tested wherever practical, and when scientific and local knowledge diverge, both should be re-examined [11]. The main message remaining here is that some variables may be suitably informed by either the fishers or scientists. Maybe the time has come for science to focus on information that is usually completely unavailable for fish species, such as the life expectancy, growth rate, stock size, and total catch per year, mainly in developing countries. Thus, the main lack of learning relies on the lack of biological and statistics surveys. Improved statistics may, in many ways, enhance fisheries science [68], and here, combined with FK, it could fill the gaps in the knowledge of fishers and scientists, leading to better management of fisheries resources. Furthermore, combining the data from different sources may expand the results, increase the reliability, and be used independently for quick and low-cost assessments.

Supporting Information

S1 Table. Fishers’ questionnaire.

Questionnaire that was applied for experts fishers.

https://doi.org/10.1371/journal.pone.0155655.s001

(DOCX)

S2 Table. Scientific model diet matrix.

Diet matrix for the 25 compartments using Ecopath with Ecosim. The predators are set in columns and the prey is set in lines.

https://doi.org/10.1371/journal.pone.0155655.s002

(DOCX)

S3 Table. Fishers’ knowledge model diet matrix.

Diet matrix for the 25 compartments using Ecopath with Ecosim. The predators are set in columns and the prey is set in lines.

https://doi.org/10.1371/journal.pone.0155655.s003

(DOCX)

Acknowledgments

We would like to thank all of the fishers from the studied communities for their kind support and patience, and the two anonymous referees who improved the paper.

Author Contributions

Conceived and designed the experiments: AHB ARC. Performed the experiments: AHB. Analyzed the data: AHB RA VC. Contributed reagents/materials/analysis tools: AHB ARC RA VC. Wrote the paper: AHB ARC RA VC.

References

  1. 1. FAO. The State of World Fisheries and Aquaculture. FAO Fisheries Proceedings. 2014. 6–32 p. Available: http://www.fao.org/3/a-i3720e.pdf
  2. 2. Pikitch EK, Santora C, Babcock EA, Bakun A, Bonfil R, Conover DO, et al. Ecosystem-Based Fishery Management. Science. 2004. (80–91). Available from: http://www.sciencemag.org/content/307/5716/1725.short
  3. 3. Hilborn R. Future directions in ecosystem based fisheries management: A personal perspective. Fish Res. 2011.
  4. 4. Berkes F. Implementing ecosystem-based management: evolution or revolution? Fish Fish. 2012. doi.wiley.com/10.1111/j.1467-2979.2011.00452.x
  5. 5. FAO. Voluntary Guidelines for Securing Sustainable Small-Scale Fisheries. 2015. Available: http://www.fao.org/3/a-i4356e.pdf.
  6. 6. Hind EJ. Knowledge Research: a Challenge To Established Fisheries Science. 2015. 72:341–58.
  7. 7. Christie P, Pollnac RB, Fluharty DL, Hixon MA, Lowry GK, Mahon R, et al. Tropical Marine EBM Feasibility: A Synthesis of Case Studies and Comparative Analyses. Coast Manag. 2009;37(3–4): 374–85.
  8. 8. Coll M, Akoglu E, Arreguín-Sánchez F, Fulton EA, Gascuel D, Heymans JJ, et al. Modelling dynamic ecosystems: venturing beyond boundaries with the Ecopath approach. Rev Fish Biol Fish. 2015; 413–24.
  9. 9. Angelini R, de Morais RJ, Catella AC, Resende EK, Libralato S. Aquatic food webs of the oxbow lakes in the Pantanal: A new site for fisheries guaranteed by alternated control? Ecol Modell. 2013;253: 82–96.
  10. 10. Pitcher TJ, Heymans JJ, Pitcher T, Forrest R, Stanford R, Martell S, et al. Back to the future: advances in methodology for modelling and evaluating past ecosystems as future policy goals. Fish Cent Res Reports. 2004;12(1):158pp.
  11. 11. Johannes R, Freeman M, Hamilton R. Ignore fishers’ knowledge and miss the boat. Fish Fish. 2000. doi.wiley.com/10.1046/j.1467-2979.2000.00019.x\nhttp://onlinelibrary.wiley.com/doi/10.1111/j.1467-2979.2000.00019.x/full
  12. 12. Bender MG, Machado GR, Silva PJ de A, Floeter SR, Monteiro-Netto C, Luiz OJ, et al. Local Ecological Knowledge and scientific data reveal overexploitation by multigear artisanal fisheries in the Southwestern Atlantic. PLoS One. 2014;9(10):e110332. pmid:25333661
  13. 13. Ramires M, Clauzet M, Barrella W, Rotundo MM, Silvano RA, Begossi A. Fishers’ knowledge about fish trophic interactions in the southeastern Brazilian coast. J Ethnobiol Ethnomed. 2015;11(1): 19.
  14. 14. Leite MCF, Gasalla MA. A method for assessing fishers’ ecological knowledge as a practical tool for ecosystem-based fisheries management: Seeking consensus in Southeastern Brazil. Fish Res. 2013;145(0): 43–53.
  15. 15. Rosa R, Carvalho AR, Angelini R. Integrating fishermen knowledge and scientific analysis to assess changes in fish diversity and food web structure. Ocean Coast Manag. 2014;102: 258–68.
  16. 16. Damasio LM, Lopes PFM, Guariento RD, Carvalho AR. Matching Fishers’ Knowledge and Landing Data to Overcome Data Missing in Small-Scale Fisheries. PLoS One. 2015. journals.plos.org/plosone/article?id=10.1371/journal.pone.0133122
  17. 17. Begossi A, Lopes P, Silvano RAM. Co-Management of reef fisheries of the Snapper-Grouper complex in a human ecological context in Brazil. Glob Prog Ecosyst Fish Manag. 2012;353–74.
  18. 18. Tesfamichael D, Pitcher TJ, Pauly D. Assessing Changes in Fisheries Using Fishers ‘ Knowledge to Generate Long Time Series of Catch Rates: a Case Study from the Red Sea. Ecol Soc. 2014;19(1):18.
  19. 19. Johannes R.E. The case for data less marine resource management: examples from tropical nearshore fisheries. Trends Ecol Evol. 1998;13(98):pp. 243–6.
  20. 20. Consulting H, River E. Using Traditional Ecological Knowledge in Science: Methods and Applications. 2000;10(5):1270–4.
  21. 21. Silvano RAM, Valbo-Jørgensen J. Beyond fishermen’s tales: contributions of fishers’ local ecological knowledge to fish ecology and fisheries management. Environ Dev Sustain. 2008;10(5): 657–75. link.springer.com/10.1007/s10668-008-9149-0
  22. 22. Grant S, Berkes F. Fisher knowledge as expert system: A case from the longline fishery of Grenada, the Eastern Caribbean. Fish Res. 2007;84(2):162–70
  23. 23. Coll M, Carreras M, Ciércoles C, Cornax MJ, Gorelli G, Morote E, et al. Assessing fishing and marine biodiversity changes using fishers’ perceptions: The Spanish Mediterranean and Gulf of Cadiz case study. PLoS One. 2014;9(1):15.
  24. 24. Coll M, Libralato S. Contributions of food web modelling to the ecosystem approach to marine resource management in the Mediterranean Sea. Fish Fish. 2012;13(1): 60–88. doi.wiley.com/10.1111/j.1467-2979.2011.00420.x
  25. 25. Brazil. Boletim da Estatística da Pesca Maítima e Estuarina do Nordeste do Brasil. 2006;6.4.1.
  26. 26. Silva AC. A pesca de pequena escala nos litorais setentrionais e orientais do Rio Grande do Norte. J Chem Inf Model. 2010;1:201.
  27. 27. Lessa R. P., Nóbrega M. F., Bezerra JLJ. Dinâmica de Populações e Avaliação de Estoques dos Recursos Pesqueiros da Região Nordeste. 2004;246.
  28. 28. Vasconcellos M, Cochrane K. Overview of world status of data-limited fisheries: Inferences from landings statistics. Fish Assess Manag Data-Limited Situations. 2005;1–20. seagrant.uaf.edu/bookstore/pubs/item.php?id=11958
  29. 29. Silvano RM, MacCord PFL, Lima R V., Begossi A. When does this fish spawn? Fishermen’s local knowledge of migration and reproduction of Brazilian coastal fishes. Environ Biol Fishes. 2006;76: 371–86.
  30. 30. MacMillan DC, Marshall K. The Delphi process—An expert-based approach to ecological modelling in data-poor environments. Anim Conserv. 2006;9(1): 11–9.
  31. 31. Froese R, Binohlan C. Empirical relationships to estimate asymptotic length, length at first maturity and length at maximum yield per recruit in fishes, with a simple method to evaluate length frequency data. J Fish Biol. 2000;56(4): 758–73. doi.wiley.com/10.1006/jfbi.1999.1194
  32. 32. Palomares MLD, Pauly D. Predicting food consumption of fish populations as functions of mortality, food type, morphometrics, temperature and salinity. Mar Freshw Res. 1998;49(5): 447.
  33. 33. Freire K, Christensen V, Pauly D. Description of the Brazil large marine ecosystem using a tropic model. Sci Mar. 2008;72(3): 477–91. scientiamarina.revistas.csic.es/index.php/scientiamarina/article/view/1005/1046
  34. 34. Gasalla M a., Rossi-Wongtschowski CLDB. Contribution of ecosystem analysis to investigating the effects of changes in fishing strategies in the South Brazil Bight coastal ecosystem. Ecol Modell. 2004;172(2–4): 283–306.
  35. 35. Guénette S, Hill RL. A trophic model of the coral reef ecosystem of La Parguera, Puerto Rico: Synthesizing fisheries and ecological data. Caribb J Sci. 2009;45(2–3): 317–37.
  36. 36. Telles MD. Modelo trofodinâmico dos Recifes em Franja do Parque Nacional Marinho dos Abrolhos—Bahia. 1988;150.
  37. 37. Pauly D. On the interrelationships between natural mortality, growth parameters, and mean environmental temperature in 175 Fish Stocks. ICES J Mar Sci. 1980;39(2):175–92. icesjms.oxfordjournals.org/content/39/2/175.short
  38. 38. Gascuel D, Morissette L, Palomares MLD, Christensen V. Trophic flow kinetics in marine ecosystems: Toward a theoretical approach to ecosystem functioning. Ecol Modell. 2008;217(1–2): 33–47.
  39. 39. Cremer MJ, Pinheiro PC, Simões-Lopes PC. Prey consumed by Guiana dolphin Sotalia guianensis (Cetacea, Delphinidae) and franciscana dolphin Pontoporia blainvillei (Cetacea, Pontoporiidae) in an estuarine environment in southern Brazil. Iheringia Série Zool. 2012;102(2): 131–7.
  40. 40. Gannon DP, Waples DM. Diets of coastal bottlenose dolphins from the U. S. Mid-Atlantic coast differ by habitat. Mar Mammal Sci. 2004;20(3): 527–45.
  41. 41. Lopes XM, da Silva E, Bassoi M, dos Santos RA, de Oliveira Santos MC. Feeding habits of Guiana dolphins, Sotalia guianensis, from south-eastern Brazil: new items and a knowledge review. J Mar Biol Assoc United Kingdom. 2012;92(8): 1–11.
  42. 42. Bornatowski H, Costa L. Hábitos alimentares de tubarões-martelo jovens Sphyrna zygaena (Carcharhiniformes: Sphyrnida), no litoral sul do Brasil. Biota Neotrop. 2007;7(1): 213–6.
  43. 43. Satoh K, Yokawa K, Saito H, Matsunaga H, Okamoto H, Uozumi Y. Preliminary stomach contents analysis of pelagic fish collected by Shoyo-Maru 2002 research cruise. Collect Vol Sci Pap ICCAT. 2004;56(3): 1096–114.
  44. 44. Sebastián Lopez K, Patricio Barría M, Roberto Meléndez C. Feeding and trophic relationships of two highly migratory sharks in the eastern south Pacific Ocean. Panam J Aquat Sci. 2012;7(1): 50–6.
  45. 45. Vaske T Jr., Vooren CM, Lessa RP. Feeding Strategy of the Night Shark (Carcharhinus signatus) and Scalloped Hammerhead Shark (Sphyrna lewini) Near Seamounts Off Northeastern Brazil. 2009;57(2): 97–104.
  46. 46. Bornatowski H, Heithaus MR, Abilhoa V, Corrêa MFM. Feeding of the Brazilian sharpnose shark Rhizoprionodon lalandii (Müller & Henle, 1839) from southern Brazil. J Appl Ichthyol. 2012;28(4): 623–7.
  47. 47. Rudershausen PJ, Buckel J a., Edwards J, Gannon DP, Butler CM, Averett TW. Feeding ecology of Blue Marlins, Dolphinfish, Yellowfin Tuna, and Wahoos from the North Atlantic Ocean and comparisons with other oceans. Trans Am Fish Soc. 2010;139(5): 1335–59.
  48. 48. Headley BM, Oxenford H a., Peterson MS, Fanning P. Size related variability in the summer diet of the blackfin tuna (Thunnus atlanticus Lesson, 1831) from Tobago, the Lesser Antilles. J Appl Ichthyol. 2009;25(6): 669–75. Adoi.wiley.com/10.1111/j.1439-0426.2009.01327.x
  49. 49. Freitas MO. Auto-ecologia de Epinephelus morio e Mycteroperca bonaci: Epinefelídeos comercialmente importantes e ameaçados no banco dos Abralhos. PhD Thesis. 2014;160.
  50. 50. Duarte LO, García CB. Diet of the mutton snapper Lutjanus analis (Cuvier) from the Gulf of Salamanca, Colombia, Caribbean Sea. Bull Mar Sci. 1999;65(2):453–65.
  51. 51. Douba VN, Guessan KJN, Kouame EP, Bahou L, Kone T, Gouli B. Food composition and feeding habits of little tunny (Euthynnus alletteratus) in continental shelf Cost d’ Ivoire (West Africa) waters of Co. 2007;1044–52.
  52. 52. Coelho M, Domingues P, Balguerias E, Fernandez M, Andrade JP. A comparative study of the diet of Loligo vulgaris (Lamarck, 1799) (Mollusca: Cephalopoda) from the south coast of Portugal and the Saharan Bank (Central-East Atlantic). Fish Res. 1997;29(3):245–55.
  53. 53. Gasalla MA, Rodrigues AR, Postuma FA. The trophic role of the squid Loligo plei as a keystone species in the South Brazil Bight ecosystem. 2010;1979:2010.
  54. 54. Leite TS, Haimovici M, Mather J. Octopus insularis (Octopodidae), evidences of a specialized predator and a time-minimizing hunter. Mar Biol. 2009;156(11): 2355–67. link.springer.com/10.1007/s00227-009-1264-4
  55. 55. Domingos SS, Athiê AA, Rossi-Wongtschowski, CDBL. Diet of Chaceon notialis (Decapoda, Brachyura) off the coast of Rio Grande, RS, Brazil. Brazilian J Oceanogr. 2007;55(4): 327–9.
  56. 56. Ivo C, Pereira J. Sinopse das principais observações sobre as lagostas Panulirus argus (Latreille) e Panulirus laevicauda (Latreille), capturadas em águas costeiras do Brasil, entre os estados do Amapá e do Espirito Santo. Bol Técnico Cient do CEPENE. 1996;4: 7–94.
  57. 57. Tararam AS, Wakabara Y, Eqüi MB. Hábitos alimentares de onze espécies de magafauna bêntica da plataforma continental de Ubatuba, SP. Inst Ocean S Paulo. 1993;10(10): 159–67.
  58. 58. Yokoyama LQ, Amaral ACZ. Allometric growth of a common Nassariidae (Gastropoda) in south-east Brazil. J Mar Biol Assoc United Kingdom. 2011;91(05): 1095–105.
  59. 59. Polovina JJ. Coral Reefs Model of a Coral Reef Ecosystem. Coral Reefs. 1984;3: 1–11.
  60. 60. Christensen V, Walters CJ. Ecopath with Ecosim: methods, capabilities and limitations. Ecol Modell. 2004;172(2–4): 109–39.
  61. 61. Christensen V, Walters CJ, Pauly D. Ecopath with Ecosim: a user’s guide. Fish Cent Res Reports. 2005;154. Available: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.122.6467&rep=rep1&type=pdf
  62. 62. Libralato S, Christensen V, Pauly D. A method for identifying keystone species in food web models. Ecol Modell. 2006;195(3–4): 153–71.
  63. 63. Otero J, Rocha F, González ÁF, Gracia J, Guerra Á. Modelling artisanal coastal fisheries of Galicia (NW Spain) based on data obtained from fishers: the case of Octopus vulgaris. Sci Mar. 2005;69(4): 577–85.
  64. 64. Begossi A. Local knowledge and training towards management. Environ Dev Sustain. 2008;10(5): 591–603.
  65. 65. Silvano RAM, Begossi A. Fishermen’s local ecological knowledge on southeastern Brazilian coastal fishes: Contributions to research, conservation, and management. Neotrop Ichthyol. 2012;10(1): 133–47.
  66. 66. Coll M, Navarro J, Olson RJ, Christensen V. Assessing the trophic position and ecological role of squids in marine ecosystems by means of food-web models. Deep Sea Res Part II Top Stud Oceanogr. 2013;95: 21–36.
  67. 67. Angelini R, Vaz-Velho F. Ecosystem structure and trophic analysis of Angolan fishery landings. Sci Mar. 2011;75(2): 309–19.
  68. 68. Pauly D, Zeller D. Catch reconstructions reveal that global marine fisheries catches are higher than reported and declining. Nat Commun. 2016;7:10244. pmid:26784963