Uncertainties regarding the natural
mortality of fish can increase due global
                                            climate change
                                            Caroline Pereira Campos1, Sandro Dimy Barbosa Bitar2 and
                                            Carlos Freitas3
                                        1
                                          Instituto Federal de Roraima, Caracarai, Roraima, Brazil
                                        2
                                          Departamento de Matemática, Universidade Federal do Amazonas, Manaus, Amazonas, Brazil
                                        3
                                          Departamento de Ciências Pesqueiras, Universidade Federal do Amazonas, Manaus, Amazonas,
                                          Brazil
                                                ABSTRACT
                                                The increase in temperature resulting from global climate change can directly affect
                                                the survival of fish and therefore population parameters such as natural mortality
                                                (M). The estimation of this parameter and the understanding of the uncertainties in
                                                its estimates are enormous challenges for studies that evaluate fish stocks.
                                                In addition, the effects of increases in temperature may be associated with life
                                                strategies. Therefore, the fuzzy set theory was used to evaluate the effects of
                                                temperature increase on the natural mortality of fish, considering different life
                                                strategies. The model showed that the increase in temperature increased the
                                                uncertainties in M estimates for all species, regardless of the life strategy. However,
                                                opportunistic species present greater uncertainties in estimates of M compared to
                                                equilibrium species. The patterns found in uncertainties of M associated with species
                                                groupings by life strategies can be used in holistic approaches for the assessment and
                                                management of recently exploited fisheries resources or for those with limited
                                                biological data.
                                            Subjects Aquaculture, Fisheries and Fish Science, Ecology, Mathematical Biology, Zoology,
                                            Climate Change Biology
Submitted 24 August 2022
                                            Keywords Climate change, Natural mortality of fish, Fuzzy theory, Uncertainty, Life strategies of
Accepted 12 February 2023
Published 7 March 2023                      fish
Corresponding author
Carlos Freitas, [email protected]
                                            INTRODUCTION
                                            Climate change is considered as one of the biggest and most complex environmental and
Academic editor
Matteo Zucchetta                            socioeconomic problem of the world at the present time (Pörtner et al., 2022). Global
Additional Information and                  warming is the result of natural alterations in the climate system and human activities,
Declarations can be found on                which contribute to rising levels of greenhouse gas (GHG) emissions (Pörtner et al., 2022).
page 16
                                            Scientific evidence indicates that anthropogenic emissions of GHGs are the main cause of
DOI 10.7717/peerj.14989                     the increase in temperature and the biggest culprit in relation to the current environmental
   Copyright                                imbalance, causing more than half of the observed increase in global average surface
2023 Pereira Campos et al.
                                            temperatures of the planet (Cleugh et al., 2011). The temperature increases on the surface
Distributed under
Creative Commons CC-BY 4.0
                                            of the planet are among the phenomena with greater predictability and with wide-ranging
                                            effects (Pörtner et al., 2022).
                                            How to cite this article Pereira Campos C, Bitar SDB, Freitas C. 2023. Uncertainties regarding the natural mortality of fish can increase due
                                            global climate change. PeerJ 11:e14989 DOI 10.7717/peerj.14989
                                       Temperature is an important environmental variable, and it is directly related to the
                                    survival of many ecothermic animals, especially fish (Baudron et al., 2014). The effects on
                                    fish populations can be manifested through changes in abundance, migration patterns and
                                    geographical distribution (Baudron et al., 2014; Ficke, Myrick & Hansen, 2007). Oremus
                                    et al. (2020) made projections using global climate change scenarios in order to verify
                                    changes in the migration and distribution of 779 species of fish from each exclusive
                                    economic zone (EEZ) and concluded that tropical countries are more vulnerable to a
                                    potential loss of stocks of these species, since these migrate to cooler waters to maintain
                                    their thermal comfort.
                                       There are patterns that correspond to different life strategies of freshwater and marine
                                    fish species (Winemiller, 1989; Winemiller & Rose, 1992; Röpke et al., 2017). The plasticity
                                    of species with different life strategies that depend on environmental variations has also
                                    been widely studied (Winemiller & Rose, 1992). Researchers have been indicating that
                                    species with life strategies characterized as having a large size, high longevity, low growth
                                    and low mortality rate are at a disadvantage in terms of temperature increase, due to their
                                    higher metabolic rate, and will likely face local extinction (Freitas et al., 2013). On the other
                                    hand, smaller species in lower trophic levels can benefit from this potential reduction of
                                    predators (Woodward et al., 2016). This shows that the effects of climate change on fish
                                    can be of a greater or lesser degree and can have distinct signs, depending on the life
                                    strategy of each species (Blanck & Lamouroux, 2007; Röpke et al., 2017).
                                       The natural mortality M is one of the most important parameters in the assessment of
                                    fish stocks and it is essential in order to estimate the productivity of a fish stock, but it is
                                    also one of the most difficult to estimate, due to limitations of data (Punt et al., 2021).
                                    The magnitude of error, bias and variance in the estimation of this parameter can be
                                    substantial and can be affected by temporal and spatial effects (Kenchington, 2014).
                                    Incorporating the uncertainties in the calculation of M estimates is essential to make
                                    robust assessments of the fish stock (Brodziak et al., 2011).
                                       Most estimation methods for M consider only intrinsic factors of the species: length,
                                    maturation size, growth rate and age, without considering the influences of extrinsic
                                    factors, such as environmental variables (e.g., temperature) (Kenchington, 2014), which can
                                    be determinants in population dynamics and are intrinsically variable (Sæther & Bakke,
                                    2000). The method proposed by Pauly (1980) to estimate M includes intrinsic and extrinsic
                                    factors, since it estimates M as a function of the von Bertalanffy (1938) growth parameters
                                    and the average annual water temperature. It has been shown that the higher the
                                    temperatures, the higher the natural mortality rates in fish (Sparre & Venema, 2007).
                                       However, the effects of temperature increase, due to climate change, incorporate high
                                    variability, which is in part due to the uncertainties existing in the climate models and
                                    underlying scenarios (Jeppesen et al., 2014). Thus, it is necessary to apply methodological
                                    approaches that incorporate variations and uncertainties in the estimation of parameters.
                                    In this sense, fuzzy logic and fuzzy set theory are presented as appropriate tools in the
                                    modeling of biological phenomena (Barros, Bassanezi & Lodwick, 2017). This theory was
                                    introduced by Zadeh (1965), with the idea of contrasting deterministic models to more
                                    flexible models, which enable mathematical modeling of uncertainty, besides
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                            2/22
                                    mathematically formulating the subjectivity inherent to natural phenomena. Fuzzy logic
                                    has been used as a tool in several issues related to the fisheries, such as assessing the
                                    vulnerability of marine species (Cheung, Pitcher & Pauly, 2005), assessment of pelagic
                                    fishery ecosystems (Paterson et al., 2007), classification of fishing areas (Sylaios,
                                    Koutroumanidis & Tsikliras, 2010), spatial management of marine fisheries (Teh & Teh,
                                    2011), behavior of fishers in pelagic fisheries (Wise et al., 2012), stock-recruitment relations
                                    (Gutiérrez-Yurrita, 2014), parameter estimates of the relationship between weight and
                                    length (Bitar, Campos & Freitas, 2016), dynamic of competition among species taking into
                                    account abiotic factors and fishing effort (Souza & Bassanezi, 2018), and the fuzzy
                                    interactivity in the generalized model of von Bertalanffy (Souza & Prata, 2019).
                                       In this study, we assumed two points: (i) the estimate of the average annual water
                                    temperature, calculated as a fuzzy number, could, therefore, allow an estimate of natural
                                    mortality that incorporates uncertainties in the face of global climate changes; and (ii) the
                                    amplitude and intensity of the effects of environmental phenomena could be
                                    species-specific (Röpke et al., 2017). We then answered three questions: (i) Will the
                                    temperature increases, in the scenarios proposed by the IPCC, increase uncertainties
                                    regarding the natural mortality in fish? (ii) Are these uncertainties, if noticeable, the same
                                    for marine and freshwater fish species? (iii) Are these uncertainties, if noticeable, the same
                                    for fish species with different life strategies?
                                    MATERIALS AND METHODS
                                    Study area
                                    The area covered in this study was the Brazilian tropical region, including continental and
                                    maritime areas (Fig. 1). The main feature of the climate of the tropics is the high average
                                    monthly temperatures. However, factors such as latitude, altitude, circulation of air masses
                                    and global climate change can cause climate differentiation.
                                    Classification of life strategies of the target species
                                    As an alternative to the concepts of r and k strategists (Pianka, 1970) and based on
                                    mechanisms of trade-offs between reproduction, growth and survival of freshwater fish,
                                    Winemiller (1989) proposed a classification using three categories of life strategies:
                                    equilibrium, opportunistic and seasonal. The classification proposed by Winemiller (1989)
                                    gave rise to the triangular continuum of strategies model, which resulted from adaptive
                                    responses to environmental predictability and variability. Subsequently, Winemiller &
                                    Rose (1992), while studying freshwater and marine fish, confirmed this classification.
                                    Basically, these patterns of life strategies suggest that:
                                       (1) Equilibrium strategists—These are fish with low growth rates and natural mortality,
                                    with a large body size, late first maturation, split spawning, large oocytes, low fecundity,
                                    broad reproductive period, relatively high longevity, parental care, which produce larger
                                    offspring in low quantities. Since they produce relatively few offspring, survival from the
                                    early stages must be high for these populations to persist. These are species that live in
                                    relatively stable environments.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                         3/22
  Figure 1 Map of the Brazilian tropical zone. The zones indicates watersheds and temperature collection points (river monitoring stations).
                                                                                                    Full-size  DOI: 10.7717/peerj.14989/fig-1
                                           (2) Seasonal strategists—These are fish that exhibit migratory behavior, absence of
                                        parental care, growth rates and natural mortality from intermediate to high, with
                                        intermediate to large body size, first maturation is intermediate to late, small oocytes, high
                                        fecundity and short periods of reproduction. This strategy would be associated with both
                                        temporal and seasonal variations and, if these conditions are favorable, a synchronization
                                        in reproduction will occur, guaranteeing juvenile individuals greater chances of growth
                                        and survival, in addition to providing a rapid population replacement.
                                           (3) Opportunistic strategists—These are fish with a small body size, early first
                                        maturation, small oocytes, low fecundity, prolonged breeding period, short life span and
                                        with little or no parental care. These characteristics are related to the ability of rapid
                                        colonization of environments marked by unpredictable temporal variations.
                                        Species selected for the study
                                        Species representing the different life strategies were selected according to the classification
                                        of Winemiller (1989) and Winemiller & Rose (1992). In addition, we sought to select species
                                        from freshwater and marine environments that have been fished for commercial purposes.
                                        Based on these criteria, three freshwater and three marine species were selected, each
                                        representing a particular type of life strategy. The growth parameters (L∞ and k), of these
                                        species can be observed in Table 1.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                                       4/22
                                       Table 1 Population parameters of the species selected for the study, where: L∞ = theoretical
                                       maximum length; k = individual growth coefficient; L50 = first maturation; Tmáx = longevity.
                                     Species                     L∞ (cm)     k (year−1)    L50 (cm)     Tmáx (year)     Strategy         Environment
                                     Cichla temensis             68.05       0.20          31.11        14.00           Equilibrum       Freshwater
                                     Prochilodus nigricans       34.60       0.44          23.08        6.80            Seasonal
                                     Triportheus angulatus       26.78       0.77          16.30        3.89            Opportunist
                                     Cynoscion acoupa            142.90      0.13          42.07        10.00           Equilibrum       Marine
                                     Lutjanus synagris           56.00       0.22          23.00        18.00           Seasonal
                                     Sardinella brasiliensis     23.68       0.26          15.80        10.73           Opportunist
                                    Note:
                                       Sources : 1—Campos, Freitas & Amadio (2015), Campos, Catarino & Freitas (2020); 2—Catarino et al. (2014), Camargo &
                                       Lima (2008); 3—Prestes et al. (2010); 4—Oliveira et al. (2020); 5—Aschenbrenner et al. (2017); 6—Costa, Tubino &
                                       Monteiro-Neto (2018).
                                       Species of freshwater environments:
                                       Equilibrium Strategist: Cichla temensis (Humboldt, 1821)
                                       Order: Perciformes. Family: Cichlidae
                                       The peacock bass, Cichla temensis, has its distribution restricted to the blackwater rivers
                                    and their tributaries in the basins of the Amazon River and Orinoco River. It has great
                                    importance in the diets of the riparian population, as well as in the ornamental fish trade
                                    and, especially, in sport fishing (Soares et al., 2007). In 2003, the peacock bass represented
                                    2.58% of the fishing production that was landed in the main ports of the state of Amazonas
                                    (Ruffino et al., 2006) and in the middle Negro River region, and it was the third most
                                    landed fish in 2013 (Inomata & Freitas, 2015). It is characterized by a preference for lentic
                                    environments, edges of lakes and sandbanks in the main channel of rivers (Winemiller,
                                    Taphorn & Barbarino-Duque, 1997). Cichla temensis is a predator at the top of the chain,
                                    with a piscivorous habit and it differs among the 15 species of the genus Cichla, since it
                                    reaches larger sizes, which can be up to 80 cm and more than 11 Kg (Winemiller, Taphorn
                                    & Barbarino-Duque, 1997; Kullander & Ferreira, 2006). The couple builds the nest and,
                                    one or both, take care of the offspring (Zaret, 1980; Jepsen, Winemiller & Taphorn, 1999).
                                    It presents split spawning, with a long reproductive period and low fecundity (Jepsen,
                                    Winemiller & Taphorn, 1999). We used these natural history data to approximate model
                                    parameters for major life history characteristics for each of the target species.
                                    We considered C. temensis to have low growth rates (k) and natural mortality (M), high
                                    values of theoretical maximum lengths (L∞) and first maturation (L50), as well as relatively
                                    high longevity (Tmax).
                                       Seasonal Strategist: Prochilodus nigricans (Spix & Agassiz, 1829)
                                       Order: Characiformes. Family: Prochilodontidae
                                       The curimatã, Prochilodus nigricans, is found in the Amazon River and Tocantins River
                                    basins (Castro & Vari, 2003). It has significant commercial value and is one of the most
                                    caught species in Brazilian continental waters. It is predominant in the landings in
                                    Amazonian cities (Batista et al., 2012) in addition to its importance in subsistence fishing
                                    (Isaac et al., 2015; Zacarkim et al., 2015). This species inhabits large rivers, flooded forests,
                                    lakes and streams, and can be found at any time of the year, even if its abundance varies,
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                                                5/22
                                    due to the influence of seasonal fluctuations in the level of the rivers in the Amazon (Silva
                                    & Stewart, 2017). It performs trophic, reproductive and dispersal migrations throughout
                                    the year and feeds, basically, on organic debris and periphyton (Silva & Stewart, 2017).
                                    Like other migratory fish that feed on detritus, it has access to an abundant source of
                                    energy, which contributes to the modulation of carbon flow and productivity of ecosystems
                                    (Silva & Stewart, 2017). It is a medium-sized species, between 35 and 50 cm long and is
                                    around 3 Kg in weight (Santos, Ferreira & Zuanon, 2006). It has high fecundity, small
                                    oocytes and short reproductive period (at the beginning of the rising water period) and
                                    does not provide parental care (Silva & Stewart, 2017). In addition, it presents high
                                    intermediate values of population parameters, such as growth constant (k), theoretical
                                    maximum length (l∞), natural mortality (M) and first maturation length (L50). Röpke et al.
                                    (2017) also classified P. nigricans as a seasonal strategist.
                                       Opportunistic Strategist: Triportheus angulatus (Spix & Agassiz, 1829)
                                       Order: Characiformes. Family: Characidae
                                       A total of 16 species have been described for the genus Triportheus, in which three are
                                    well known, Triportheus albus (Cope, 1872), T. auritus (Valenciennes, 1850) and
                                    T. angulatus (Spix & Agassiz, 1829), and their occurrences are recorded in the basins of the
                                    Amazon, Tocantins and Orinoco Rivers (Malabarba, 2004). In recent years, the demand
                                    and commercial interest for Amazonian fish species that are numerous, but of smaller size,
                                    has increased significantly. For example, in the 1970s, sardines were infrequently
                                    commercialized (Petrere Jr., 1978), but in the past decade, there has been an increase in the
                                    volume of landings in the ports of Manaus and Manacapuru, Amazonas state, which has
                                    risen from 2% to 12% (Gonçalves & Batista, 2008). Triportheus angulatus is a pelagic fish
                                    that commonly lives near the surface and close to the banks of rivers and lakes. It has an
                                    omnivorous diet that includes insects, zooplankton, fruits and seeds, and it has been
                                    verified that, even under conditions of low oxygen concentrations, its feeding activity is not
                                    influenced as it has the ability to modify its lower lips to absorb surface oxygen from water
                                    (Yamamoto, Soares & Freitas, 2004).
                                       Although T. angulatus does not meet some characteristics to be assumed as a typical
                                    opportunistic species, it is one of the exploited species by fishing fleets in freshwater
                                    systems that shows more characteristics of this type of life strategy, such as small body size,
                                    ranging between 20 and 24 cm (Malabarba, 2004); low fecundity and a single spawning
                                    (Araújo et al., 2012); no parental care and high growth (k) and natural mortality (M) rates,
                                    low value of theoretical maximum lengths (L∞), small first maturation lengths (L50) and
                                    low life span (Tmax).
                                       Species of marine environments:
                                       Equilibrium Strategist: Cynoscion acoupa (Lacepède, 1801)
                                       Order: Perciformes. Family: Sciaenidae
                                       The red mullet, Cynoscion acoupa, has a wide distribution and occupies coastal waters
                                    from Panama, in Central America, to Argentina, in southern Latin America as well as most
                                    of the Brazilian coast (except in the extreme south of the country), both in marine and
                                    estuarine environments (Cervigón, 1993). Cynoscion acoupa is the third most landed fish
                                    species in Brazil and has a high commercial value due to the quality of its meat and for its
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                          6/22
                                    swim bladder, which is used as a raw material in the production of emulsifiers and clarifiers
                                    (Barletta, Barletta-Bergan & Saint-Paul, 1998). Cynoscion acoupa is a species with
                                    nectonic, demersal and coastal habits that lives in shallow and brackish waters under the
                                    influence of estuaries and mangroves, which, for this species, are places of refuge, feeding
                                    and reproduction (Szpilman, 2011). It feeds mainly on fish and crustaceans and is a large
                                    species that can reach up to 130 cm in length and 20 Kg in weight (Ferreira et al., 2016).
                                    It has batch spawning, a long reproductive period and high fecundity (a common
                                    characteristic for marine species in equilibrium) (Almeida et al., 2016). It has low growth
                                    rates (k) and high values of theoretical maximum lengths (L∞) and first maturation
                                    lengths (L50). Its longevity is considered as intermediate to high (Tmáx).
                                       Seasonal Strategist: Lutjanus synagris (Linnaeus, 1758)
                                       Order: Perciformes. Family: Lutjanidae
                                       The lane snapper, Lutjanus synagris, has an area of occurrence that extends from North
                                    Carolina (USA) to southeastern Brazil (Manooch & Mason, 1984). It is an important
                                    fishing resource of tropical marine waters, due to its abundance and quality of meat, which
                                    gives it high commercial value and, as a result, it is the target of commercial, artisanal and
                                    sports fisheries (Rezende, Ferreira & Frédou, 2003). Lutjanus synagris is a species of
                                    demersal habits that lives in warm waters, and is often associated with rocky and coral
                                    bottoms, between the coastal zone and out to depths of about 400 m (Pimentel & Joyeux,
                                    2010). It feeds mainly on crustaceans and fish. Considered a medium-sized species, it can
                                    reach lengths of between 40 and 50 cm (Pimentel & Joyeux, 2010). It has batch spawning
                                    and high fecundity (Grimes, 1987; Freitas et al., 2014). It presents intermediate to slow
                                    growth (k), as well as intermediate values of theoretical maximum lengths (L∞) and first
                                    maturation (L50), and is relatively long-lived (Tmax). King & McFarlane (2003) concluded
                                    that species of the family Lutjanidae are seasonal strategists.
                                       Opportunistic Strategist: Sardinella brasiliensis (Steindachner, 1879)
                                       Order: Clupeiformes. Family: Clupeidae
                                       The Brazilian sardinella, Sardinella brasiliensis, is a species that is endemic to the
                                    southeastern coast of Brazil, and is found along the continental shelf, between Cabo de São
                                    Tomé, RJ and Cabo de Santa Marta Grande, SC (Paiva & Motta, 2000). In terms of volume
                                    of production, it is the most important marine fishing resource in Brazil since it is the most
                                    commercialized and consumed fish in the country (Paiva, 1997). Shoals occur near the
                                    surface, in water up to 80 m deep, and reduce their frequency as the depth increases (Paiva
                                    & Motta, 2000). In the region where it occurs, specific oceanographic characteristics are
                                    observed, due to the periodic entry of infiltrations from the Central South Atlantic water.
                                    These waters act as a fertilization mechanism and provide a greater amount of plankton
                                    that favors the survival of the larvae of S. brasiliensis (Matsuura, 1998). It is a species of
                                    planktophagous feeding habits (Kurtz & Matsuura, 2001), of small size, ranging from 9 to
                                    27 cm, with a high fertility rate, sequenced spawning and without parental care (Costa,
                                    Tubino & Monteiro-Neto, 2018). In addition, it also presents a high natural mortality rate
                                    (M) and early maturation (L50), with low values for theoretical maximum lengths (L∞).
                                    Costa, Tubino & Monteiro-Neto (2018) also classified S. brasiliensis as an opportunistic
                                    strategist.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                         7/22
                                       Table 2 Area of study of the target species and their respective basins of origin.
                                     Study No                                 Study area                                          Basin
                                     1                                        Negro River                                         1—Amazon River
                                     2                                        Solimões River                                      1—Amazon River
                                     3                                        Solimões River                                      1—Amazon River
                                     4                                        Baía de São Marcos                                  3—N/NE Atlantic
                                     5                                        Banco dos Abrolhos                                  5—Eastern Atlantic
                                     6                                        Costa do Itaipu                                     5—Eastern Atlantic
                                    Note:
                                       Sources : 1—Campos, Freitas & Amadio (2015), Campos, Catarino & Freitas (2020); 2—Catarino et al. (2014), Camargo &
                                       Lima (2008); 3—Prestes et al. (2010); 4—Oliveira et al. (2020); 5—Aschenbrenner et al. (2017); 6—Costa, Tubino &
                                       Monteiro-Neto (2018).
                                    Data collection
                                    To estimate the fuzzy set of natural mortality (M) we use: (i) input variables for the natural
                                    mortality equation (M) proposed by Pauly (1980): von Bertalanffy growth parameters
                                    (Table 1) of species and annual average water surface temperature (National Water Agency
                                    of Brazil (ANA), 2020); and (ii) environmental variable in the modeling: temperature
                                    estimates proposed by the Intergovernmental Panel on Climate Change (IPCC, 2014).
                                    Average annual temperature
                                    The National Water Resources Information System (SNIRH), whose organization,
                                    implementation and management are the responsibility of the National Water Agency,
                                    divides the national territory into eight large basins: 1—Amazon River, 2—Tocantins
                                    River, 3—Atlantic (northern/northeastern stretch), 4—São Francisco River, 5—Atlantic
                                    (eastern stretch), 6—Paraná, 7—Uruguay, and 8—Atlantic (southeastern stretch)
                                    (National Water Agency of Brazil (ANA), 2020). The temperature data concerning the
                                    basins 1, 3 and 5 (Fig. 1), where the studies of estimation of growth parameters of the
                                    species were carried out (Table 2), were obtained and used in the modeling.
                                    IPCC scenarios
                                    In 2014, the Intergovernmental Panel on Climate Change (IPCC, 2014) proposed four
                                    more simplified scenarios, but using a more complete system, known as RCPs
                                    (representative concentration pathways), which use the amount of energy absorbed
                                    (radiative forcing, in W/m2) by greenhouse gases (GHG), during or near the end of the 21st
                                    century. RCP 2.6: mitigation scenario, which leads to a very low level of absorption; RCP
                                    4.5 and RCP 6.0: stabilization scenarios and RCP 8.5: scenario with very high greenhouse
                                    gas emissions (Table 3) (IPCC, 2014).
                                    Data analysis
                                    Foundations of the fuzzy set theory
                                    To understand the modeling process presented in this study, it is necessary to define
                                    concepts related to the theory of fuzzy sets. All concepts were defined according to Barros,
                                    Bassanezi & Lodwick (2017) and are detailed in Appendix S1.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                                                8/22
  Table 3 Features of the RCP scenario of the fifth IPCC report (AR5).
 RCP Radioactive forcing           Concentration CO2-         Behavior                         Rise in sea level Elevation of T C on the planet
                                   equiv. (ppm)                                                (cm)              in 2,100
 2.6    Peak of <3 W.m−2 in 2,100 490                         Rising with a peak in 2,040 and 26–55              0.3–1.7
                                                               declining
 4.5    Additional storage of ≈4.5 650                        Rising until 2,060 and stabilizing 32–63           1.1–2.6
         W.m−2
 6.0    Additional storage of ≈6   850                        Rising until 2,100 and stabilizing 33–63           1.4–3.1
         W.m−2
 8.5    Storage of ≈8.5 W.m−2      >1370                      Rising until 2,100               45–82             2.6–4.8
Note:
   Source: IPCC (2014).
                                           Estimation of fuzzy natural mortality
                                           The fuzzy set for natural mortality of each species was performed in four stages:
                                              (i) Fuzzy set estimate for average annual temperature:
                                              For the fuzzification of average annual temperature, the normality of the data
                                           distribution was first verified using the Shapiro–Wilk normality test and its variance. Then,
                                           the average annual temperature was modeled as a fuzzy number, from a set of confidence
                                           intervals, based on the methodology of Buckley (2005).
                                              Buckley (2005) suggests finding the confidence interval 100 (1 − β)% for all 0.01 ≤ β < 1.
                                           Each of these intervals can be denoted by: [θ1(β), θ2(β)].
                                              The range is considered [θ , θ ] for β = 1, where θ it is the point estimate for the
                                           parameter θ, thus, there are intervals                                                       ^
                                                                            h ia for 0.01 ≤ β ≤ 1. These ranges define the fuzzy set h
                                           through its a-levels as
                                                                 h follows:
                                                                   ia        h^ ¼ ½u1 ðaÞ; u2 ðaÞ.
                                              For 0 ≤ a < 0.01, h^ ¼ ½u1 ð0:01Þ; u2 ð0:01Þ was defined. With this, from the
                                           confidence intervals for a parameter θ, a fuzzy set for the average annual temperature was
                                           constructed whose a-cuts are those intervals. The fuzzy average annual temperature was
                                           estimated for each set of temperature data, according to the basin of origin in the study
                                           area for each species.
                                              (ii) Fuzzy modeling for the temperature variations suggested by the IPCC (Di)    ~
                                              For each IPCC scenario, a fuzzy average annual temperature was projected.
                                           The membership functions, used for each IPCC scenario, were assumed to be the
                                           characteristic function of the set of each temperature range, since it is believed that all
                                           values of this range have the same importance, since there is no evidence of differences.
                                              For each interval Di ¼ ½Dmi ; DM i , with ði ¼ 1; 2; 3; 4Þ and where Dmi is the minimum
                                           value and DMi is the maximum value in the IPPC interval, associated with the temperature
                                           variation suggested by the IPCC, a fuzzy set D   ~ i is assigned. Therefore, a membership
                                           function is defined to characterize this fuzzy set. We will use the characteristic function (or
                                           indicator function) of each set Di . With this configuration, the a-levels of such sets are
                                           ½Aa ¼ Di , 0  a  1. Algebraic operations on these fuzzy sets are performed on their
                                           a-levels, based on interval arithmetic (Moore, 2009). These fuzzy sets will be used to
                                           estimate the average temperature variation under the influence of variations predicted by
                                           the IPCC.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                                         9/22
  Figure 2 Fuzzy annual average temperature. (A) Amazon River, (B) Northern—Eastern Atlantic and (C) Eastern Atlantic. Each curve was
  obtained as fuzzy set to represent the interval of values of temperature.                  Full-size  DOI: 10.7717/peerj.14989/fig-2
                                          (iii) Membership function for average annual temperature values under IPCC scenarios
                                       e
                                      (T i , i ¼ 1; 2; 3; 4)
                                          Let D ~ i the diffuse variation, associated with scenario i, obtained with the methodology
                                      presented above. The fuzzy set estimate for average temperature T~i , relative to scenario i,
                                      will be calculated with the sum of two fuzzy numbers, acocording to the following
                                      equation:
                                          Tei ¼ T ~ D ~ i , where  denotes sum of fuzzy numbers (Barros, Bassanezi & Lodwick,
                                      2017).
                                          It is noteworthy that, (Theorem 2.4, Barros, Bassanezi & Lodwick, 2017)
                                                  a h ia
                                      ½f
                                       Ti a ¼ T   ~ þ D     ~i
                                         and, by Theorem 1.4 in Barros, Bassanezi & Lodwick (2017), it is possible to characterize
                                      ~
                                      Ti through its a-levels.
                                         (iv) Fuzzy set for natural mortality
                                         In each scenario, the fuzzy natural mortality was estimated by applying the Zadeh’s
                                      Extensions Principle (Zadeh, 1965) of the classical function proposed by Pauly (1980):
                                      lnðM Þ ¼ 0:0152  0:279 lnðL1Þ þ 0:6543 lnðkÞ þ 0:463 lnðT Þ                                  (1)
                                      where L∞, k and T denote theoretical maximum length, individual growth coefficient and
                                      average annual water surface temperature, respectively.
                                         The procedures for calculating the Zadeh extension are presented in the pseudo code,
                                      detailed in Appendix S2. We coded our simulations in Python programming language,
                                      version 3.10.6, along with matplotlib, scipy, numpy and pandas packages, on a computer
                                      equipped with an Intel(R) Core (TM) i5-2450M and 2.50 GHz CPU. The four stages are
                                      contained in pseudocode, and documented in Appendix S2.
                                      RESULTS
                                      The fuzzy sets for the average annual temperature for the Amazon River basin have a
                                      temperature range between 24  C and 32  C (Fig. 2A), while the northern/northeastern
                                      Atlantic basin has a temperature range of between 26  C and 33  C (Fig. 2B) and the
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                              10/22
  Figure 3 Distribution of possibilities of natural mortality (M) of freshwater species. (A) Cichla temensis, (B) Prochilodus nigricans and
  (C) Triportheus angulatus. The vertical bar on the right indicates the degree of membership for each IPCC (Intergovernmental Panel on Climate
  Change) scenario, after the application of the Zadeh extension. The vertical axis on the left presents the possibilities of M. The bar on the right, in
  grayscale, presents the degree of membership, and the closer to 1 (one), the greater the pertinence of M.C1 = RCP 2.6; C2 = RCP 4.5; C3 = RCP 6.0;
  C4 = RCP 8.5 (IPCC, 2014).                                                                               Full-size  DOI: 10.7717/peerj.14989/fig-3
                                           eastern Atlantic basin ranges between 19  C and 30  C (Fig. 2C). Taking into account the
                                           highest degree of membership, that is, 1 (one), the average temperature is 27.3  C, 28.8  C
                                           and 25.3  C, for the Amazon River, northern/northeastern and eastern basins, respectively.
                                               For the three freshwater species, a pattern of increased uncertainty of M with increasing
                                           temperature was observed. This uncertainty is most evident between scenarios RCP 2.6
                                           and RCP 8.5. Cichla temensis, which presents strategy in equilibrium, presented the lowest
                                           degree of uncertainty of M among the freshwater species (Fig. 3A). Among these species,
                                           Triportheus angulatus, that presents some characteristics of opportunistic species,
                                           presented the highest degree of uncertainty of M for all scenarios (Fig. 3B). Prochilodus
                                           nigricans, a seasonal species, presented an intermediate degree of uncertainty when
                                           compared to other freshwater species (Fig. 3C).
                                               In the same manner as was presented for the freshwater species, the observed pattern is
                                           an increase in the uncertainty of M with the increase in temperature for the three saltwater
                                           species, and this is more evident between scenarios RCP 2.6 and RCP 8.5. When the life
                                           strategies were compared, Cynoscion acoupa, an equilibrium strategist, presented the
                                           lowest degree of uncertainty between the species (Fig. 4A). This increased for Lutjanus
                                           synagris (Fig. 4B), which is a seasonal strategist, and Sardinella brasiliensis, an opportunist,
                                           presented the highest degree of uncertainty of M, among the species (Fig. 4C). Considering
                                           a degree of relevance greater than or equal to 0.05, the a-cut associated with the
                                           distribution of the possibility of M for each species and for each scenario of the IPCC can
                                           be seen in Table 4.
                                               It was possible to identify a general pattern associated with estimates of M, which
                                           resulted from IPCC scenarios, and two patterns related to the interactions between the life
                                           strategy and the type of environment inhabited by the species. The uncertainty of M
                                           increases gradually from one scenario to another, and is consistently higher in RCP 8.5 for
                                           all species. However, the variation in the uncertainty is not constant between species with
                                           different life strategies. Species in equilibrium have lower uncertainty in estimates of M in
                                           all scenarios. Seasonal species present estimates of M with intermediate values of
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                                                 11/22
  Figure 4 Distribution of possibilities of natural mortality (M) of marine species. (A) Cynoscion acoupa, (B) Lutjanus synagris and (C) Sardinella
  brasiliensis. The vertical bar on the right indicates the degree of membership for each IPCC (Intergovernmental Panel on Climate Change) scenario,
  after the application of the Zadeh extension. The vertical axis on the left presents the possibilities of M. The bar on the right, in grayscale, presents the
  degree of membership, and the closer to 1 (one), the greater the pertinence of M.C1 = RCP 2.6; C2 = RCP 4.5; C3 = RCP 6.0; C4 = RCP 8.5 (IPCC,
  2014).                                                                                                        Full-size  DOI: 10.7717/peerj.14989/fig-4
  Table 4 Natural mortality possibilities (M) for each target species in the study and IPCC (Intergovernmental Panel on Climate Change)
  scenarios, considering a degree of relevance greater than or equal to 0.05.
 Species                        Current           RCP 2.6            RCP 4.5           RCP 6.0           RCP 8.5           Strategy              Environment
 Cichla temensis                0.48              0.49–0.50          0.51              0.49–0.51         0.50–0.53         Equilibrium           Freshwater
 Prochilodus nigricans          0.98–0.99         0.99–1.02          1.00–1.03         1.01–1.04         1.03–1.07         Seasonal
 Triportheus angulatus          1.52–1.54         1.53–1.58          1.55–1.60         1.56–1.62         1.59–1.66         Opportunist
 Cynoscion acoupa               0.30              0.31               0.31              0.31              0.33              Equilibrium           Marine
 Lutjanus synagris              0.52–0.53         0.52–0.55          0.53–0.56         0.54–0.56         0.55–0.58         Seasonal
 Sardinella brasiliensis        0.74–0.76         0.75–0.78          0.76–0.79         0.76–0.80         0.78–0.82         Opportunist
                                            uncertainty, while opportunistic species or with characteristics that come close to the
                                            opportunistic strategy, such as T. angulatus, exhibit the greatest uncertainties in the
                                            estimated M values for all scenarios, and the variation in the increase in natural mortality is
                                            slightly greater for these species. Apparently, species that inhabit freshwater systems have a
                                            higher degree of uncertainty of M when compared to marine species.
                                            DISCUSSION
                                            The difficulty in estimating natural mortality is associated with the number of factors that
                                            influence it and the complexity of taking into account all or most of these factors in the
                                            estimation process. To minimize this problem, Brodziak et al. (2011) suggested considering
                                            the relative influence of intrinsic vs. extrinsic factors and using the estimated distribution
                                            of M to provide a more accurate approximation of parametric uncertainties. Our model
                                            used temperature as a key factor, which can influence the ratio of metabolic rate and body
                                            mass as an intrinsic factor, as well as an environmental and extrinsic factor. Therefore, this
                                            can be considered a more robust model compared to the models that have already been
                                            used in the attempt to incorporate the uncertainties about M, but that do not combine the
                                            two groups of factors that affect M or do not use a force variable that can influence the
                                            intrinsic and extrinsic factors (MacCall, 2011; Punt et al., 2021).
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                                                       12/22
                                       In a recent study, Dahlke et al. (2020) analyzed 694 freshwater and marine fish species
                                    from all climate zones and concluded that 60% of fish species may not survive if climate
                                    warming reaches the worst-case scenario, i.e., increasing to around 4  C. However, for a
                                    more optimistic scenario, with an increase of up to 1.5  C, only 10% of the species surveyed
                                    would be at risk in the next 80 years. The physiology of fish is directly related to
                                    temperature, which has an inverse relationship with the amount of dissolved oxygen in the
                                    aquatic environment (Ficke, Myrick & Hansen, 2007). This inverse relationship (increase
                                    in temperature and the reduction of oxygen level) led to an increase in metabolic rate, an
                                    increase in physiological stress, and an imbalance in biochemical functions, as well as an
                                    increase in the mortality rate of eggs, larvae and adult fish, thus changing the structure,
                                    composition and population dynamics of the fish species (Freitas et al., 2013; Barros &
                                    Albernaz, 2013). Moreover, the natural mortality is related to body length, asymptotic
                                    length, and growth rate (Gislason et al., 2010).
                                       The changes in fish populations that are induced by climate change seem to be more
                                    noticeable in smaller species with a higher growth rate (Levangie, Blanchfield & Hutchings,
                                    2021). These authors investigated the effects on M of a 10% reduction in the asymptotic
                                    length of marine fish species of Canada’s Scotian shelf, associated with a 1  C increase in
                                    ocean temperature, and found that species with smaller body size suffered greater increases
                                    in M when compared to species with a larger body size. Additionally, considering the 1  C
                                    increase in ocean temperature, Wang et al. (2020) showed that small-sized, fast-growing
                                    species were more likely to experience declines in population growth than larger,
                                    slow-growing species. These studies corroborate our results, which show that
                                    opportunistic species presented the highest M values, in addition to the greatest increases
                                    in uncertainty of M (Table 4).
                                       The effects of climate change due to an increase in temperature may vary by region and
                                    cause adverse environmental events, such as hydrological, oceanographic alterations, as
                                    well as extreme climatic events (Marengo et al., 2009). One reason that would explain the
                                    additional natural mortality of freshwater species, which in this study are Amazonian
                                    species (Cichla temensis, Prochilodus nigricans and Triportheus angulatus), would be the
                                    environmental changes associated with extreme climatic events in the Amazon region,
                                    such as extreme flooding and droughts (Freitas et al., 2013). The marine species (Cynoscion
                                    acoupa, Lutjanus synagris and Sardinella brasiliensis) may have been affected by impacts
                                    on coastal circulation patterns, such as coastal resurgences and South Atlantic central
                                    water infiltrations (Garcia et al., 2004; Diaz & Rosenberg, 2008; Gruber, 2011). Compared
                                    to marine species, freshwater species are more vulnerable to habitat loss, which can result
                                    from geographical barriers and environmental degradation induced by anthropogenic
                                    actions, such as hydroelectric dams and pollution (Comte & Olden, 2017). In addition, the
                                    thermal tolerance ranges of freshwater species are on average 1  C greater than for marine
                                    species, which probably reflects the greater temperature variability in lakes and rivers and
                                    greater stability of the oceans (Dahlke et al., 2020), and which could explain the lower
                                    number of uncertainties about M for marine species.
                                       The uncertainties and variability of estimates of M observed in this study are related to
                                    the life strategies of each species. Many aspects of the life history of fish are related to
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                       13/22
                                    climate on an evolutionary scale, which gives them high or low plasticity, i.e., the ability to
                                    adapt to natural changes in their habitats (Val & Almeida-Val, 1995). However, the climate
                                    changes that are currently underway are occurring at a higher rate than the adaptation
                                    rates, and it is necessary to understand, from the aspects of the life history of a species, if
                                    the species has the potential to deal with these changes in the environment (Dahlke et al.,
                                    2020; Sunday, 2020), since life history traits are the underlying determinants for population
                                    responses to environmental forces (King & McFarlane, 2003).
                                        Several studies have been carried out that relate the characteristics of the life history of
                                    fish species to temperature in the face of climate change. Studies comparing lake fish from
                                    temperate and tropical regions along latitudinal gradients have shown that fish species
                                    from lower latitudes are often smaller, grow faster, mature earlier, have a shorter life
                                    expectancy, and allocate less energy for reproduction than species from higher latitudes
                                    (Blanck & Lamouroux, 2007; Jeppesen et al., 2014; Meerhoff et al., 2012). Jeppesen et al.
                                    (2014) conducted a literature review and synthesized the expected changes in the main
                                    characteristics of fish life history of inland water systems. In a temperature gradient,
                                    species with opportunistic life history characteristics showed greater plasticity compared to
                                    species in equilibrium. It has also been shown that the influence of temperature is greater
                                    and more significant for fish populations with opportunistic characteristics (faster growth,
                                    early age of sexual maturity and shorter life expectancy), than fish with slower life histories,
                                    such as species in equilibrium (Free et al., 2019), which was corroborated by our results.
                                        King & McFarlane (2003) analyzed the life history characteristics of 42 commercially
                                    important fish species from the west coast of Canada. The authors grouped the species by
                                    life strategies and made inferences in responses to environmental conditions for fisheries
                                    management purposes. They concluded that equilibrium strategists can withstand only
                                    low or moderate rates of capture, since they have low fecundity and late maturation and,
                                    therefore, are not able to recover as quickly as other fish after population reduction by
                                    fishing and, because they have low growth rates, their population dynamics have a very low
                                    variability. As seasonal strategists are relatively long-lived, they could benefit from this by
                                    ensuring a relatively long reproductive cycle, which minimizes the risk that periods of
                                    unfavorable environmental conditions might result in the loss of a stock. However, these
                                    species exhibit lower variability in abundance and were classified as having a steady-state
                                    population pattern. Opportunistic species occupy habitats not only with a high degree of
                                    variability, but also with large energy resources, so their high population abundance
                                    variability.
                                        These studies corroborate our results about increased variability and uncertainty in
                                    natural mortality estimates (M) and, consequently, in regard to population dynamics.
                                    The increase in temperature has less influence and, therefore, a lower degree of uncertainty
                                    of M for strategists in equilibrium (Cichla temensis and Cynoscion acoupa) regarding
                                    climate change. The seasonal strategists (Prochilodus nigricans and Lutjanus synagris)
                                    presented an intermediate and/or stationary degree of uncertainty of M and the
                                    opportunists or with characteristics that come close to the opportunistic strategy
                                    (Triportheus angulatus and Sardinella brasiliensis) presented the highest degree of
                                    uncertainty of M. This suggests that the adaptation and plasticity of fish species in order to
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                          14/22
                                      Figure 5 Expected behavior according to the characteristics of the life history of fish according to a
                                      gradient of temperature ( C) and natural mortality (M) (adapted from Jeppesen et al., 2014). Each
                                      groups of species are described by key life history traits. Full-size  DOI: 10.7717/peerj.14989/fig-5
                                    deal with current climate changes, associated with life strategies, follow a gradient of
                                    temperature and natural mortality (M), in addition to a general behavior in the variability
                                    of the population dynamics (Fig. 5).
                                       The natural mortality (M) is directly related to the productivity of the stock and the
                                    yields that can be obtained, in addition to being one of the most sensitive parameters of fish
                                    stock assessment models. Therefore, it is essential to use M estimates that incorporate
                                    variations and uncertainties in the face of global climate change. In addition, the patterns
                                    found in uncertainties about M associated with species groupings by life strategies can be
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                                   15/22
                                    used in more holistic approaches to the assessment and management of recently exploited
                                    fishery resources or those with limited biological data.
                                     ADDITIONAL INFORMATION AND DECLARATIONS
                                    Funding
                                    This work was supported by Coordenação de Apoio de Pessoal de Nível Superior—
                                    CAPES. The funders had no role in study design, data collection and analysis, decision to
                                    publish, or preparation of the manuscript.
                                    Grant Disclosures
                                    The following grant information was disclosed by the authors:
                                    Coordenação de Apoio de Pessoal de Nível Superior—CAPES.
                                    Competing Interests
                                    The authors declare that they have no competing interests.
                                    Author Contributions
                                     Caroline Pereira Campos conceived and designed the experiments, performed the
                                      experiments, analyzed the data, prepared figures and/or tables, authored or reviewed
                                      drafts of the article, and approved the final draft.
                                     Sandro Dimy Barbosa Bitar conceived and designed the experiments, analyzed the data,
                                      prepared figures and/or tables, authored or reviewed drafts of the article, and approved
                                      the final draft.
                                     Carlos Freitas conceived and designed the experiments, authored or reviewed drafts of
                                      the article, and approved the final draft.
                                    Data Availability
                                    The following information was supplied regarding data availability:
                                      The mathematical concepts and definitions of Fuzzy Set Theory are available in the
                                    Supplemental File.
                                    Supplemental Information
                                    Supplemental information for this article can be found online at http://dx.doi.org/10.7717/
                                    peerj.14989#supplemental-information.
                                    REFERENCES
                                    Almeida ZS, Santos NB, Sousa HL, Neta RNFC, Andrade TSOM. 2016. Biologia reprodutiva da
                                      pescada amarela (Cynoscion acoupa) capturada na baía de São Marcos, Maranhão, Brasil. Biota
                                      Amazônia 6(1):46–54 DOI 10.18561/2179-5746/biotaamazonia.v6n1p46-54.
                                    Araújo AS, Lima LTB, Nascimento WS, Yamamoto ME, Chellappa S. 2012. Características
                                      morfométricas-merísticas e aspectos reprodutivos da sardinha de água doce, Triportheus
                                      angulatus (Osteichthyes: Characiformes) do rio Acauã do bioma Caatinga. Biota Amazônia
                                      2(1):59–73 DOI 10.18561/2179-5746/biotaamazonia.v2n1p59-73.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                       16/22
                                    Aschenbrenner A, Freitas MO, Rocha GRA, de Moura RL, Francini-Filho RB, Minte-Vera C,
                                      Ferreira BP. 2017. Age, growth parameters and fisheries indices for the lane snapper in the
                                      Abrolhos Bank, SW Atlantic. Fisheries Research 194(Suppl. A):155–163
                                      DOI 10.1016/j.fishres.2017.06.004.
                                    Barletta M, Barletta-Bergan A, Saint-Paul U. 1998. Description of the fisheries structure in the
                                      mangrove-dominated region of Bragança (State of Pará, North Brazil). Ecotropica 4:41–53.
                                    Barros DF, Albernaz ALM. 2013. Possible impacts of climate change on wetlands and its biota in
                                      the Brazilian Amazon. Brazilian Journal of Biology 74(4):810–820
                                      DOI 10.1590/1519-6984.04013.
                                    Barros LC, Bassanezi RC, Lodwick WA. 2017. A first course in fuzzy logic, fuzzy dynamical
                                      systems, and biomathematics. Studies in Fuzziness and Soft Computing. Vol. 347. Berlin:
                                      Springer.
                                    Batista VS, Isaac VJ, Fabré NN. 2012. A Produção Desembarcada por Espécie e sua Variação por
                                      Macrorregião Amazônica. In: Batista VS, Isaac VJ, eds. Peixes e pesca no Solimões Amazonas:
                                      uma avaliação integrada. Brasília: IBAMA/ProVárzea, 107–133.
                                    Baudron AR, Needle C, Rijnsdorp AD, Marshall CT. 2014. Warming temperatures and smaller
                                      body sizes: synchronous changes in growth of North Sea fishes. Global Change Biology
                                      20(4):1023–1031 DOI 10.1111/gcb.12514.
                                    Bitar SD, Campos CP, Freitas CEC. 2016. Applying fuzzy logic to estimate the parameters of the
                                       length-weight relationship. Brazilian Journal of Biology 76(3):611–618
                                       DOI 10.1590/1519-6984.20014.
                                    Blanck A, Lamouroux N. 2007. Large-scale intraspecific variation in life-history traits of European
                                       freshwater fish. Journal of Biogeography 34(5):862–875 DOI 10.1111/j.1365-2699.2006.01654.x.
                                    Brodziak J, Ianelli J, Lorenzen K, Methot RD. 2011. Estimating natural mortality in stock
                                      assessment applications. Seattle: Alaska Fisheries Science Center, U.S. Department of Commerce,
                                      National Oceanic and Atmospheric Administration, National Marine Fisheries Service, NOAA
                                      Technical Memorandum, 119.
                                    Buckley JJ. 2005. Fuzzy statistics: hypothesis testing. Soft Computing 9(7):512–518
                                      DOI 10.1007/s00500-004-0368-5.
                                    Camargo M, Lima WMA Jr. 2008. Aspectos da biologia reprodutiva de seis espécies de peixes de
                                      importância comercial do médio Rio Xingu-bases para seu manejo. Scientific Magazine UAKARI
                                      3(1):64–77 DOI 10.31420/uakari.v3i1.20.
                                    Campos CP, Catarino MF, Freitas CEC. 2020. Stock assessment of the peacock bass Cichla
                                      temensis (Humboldt, 1821), an important fishing resource from the middle Negro River,
                                      Amazonas, Brazil. Brazilian Journal of Biology 80(3):506–510 DOI 10.1590/1519-6984.203124.
                                    Campos CP, Freitas CEC, Amadio S. 2015. Growth of the Cichla temensis Humboldt, 1821
                                      (Perciformes: Cichlidae) from the middle rio Negro, Amazonas, Brazil. Neotropical Ichthyology
                                      13(2):413–420 DOI 10.1590/1982-0224-20140090.
                                    Castro RMC, Vari RP. 2003. Family Prochilodontidae. In: Reis RE, Kullander SO, Ferraris JR, eds.
                                      Check List of the Freshwater Fishes of South and Central America. Porto Alegre: EDIPUCRS,
                                      65–70.
                                    Catarino MF, Campos CP, Garcez R, Freitas CEC. 2014. Population dynamics of Prochilodus
                                      nigricans caught in Manacapuru Lake (Amazon basin, Brazil). Boletim do Instituto de Pesca
                                      40(4):589–595.
                                    Cervigón F. 1993. Los peces marinos de Venezuela. Second Edition. Caracas: Fundación Científica
                                      Los Roques.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                            17/22
                                    Cheung WWL, Pitcher TJ, Pauly D. 2005. A fuzzy logic expert system to estimate intrinsic
                                      extinction vulnerabilities of marine fishes to fishing. Biological Conservation 124(1):97–111
                                      DOI 10.1016/j.biocon.2005.01.017.
                                    Cleugh H, Smith MS, Battaglia M, Graham P. 2011. Climate change: science and solutions for
                                       Australia. Campbell: Commonwealth Scientific and Industrial Research Organisation (CSIRO)
                                       Publishing.
                                    Comte L, Olden JD. 2017. Climatic vulnerability of the world’s freshwater and marine fishes.
                                       Nature Climate Change 7(10):718–722 DOI 10.1038/nclimate3382.
                                    Costa MR, Tubino RA, Monteiro-Neto C. 2018. Length-based estimates of growth parameters
                                       and mortality rates of fish populations from a coastal zone in the Southeastern Brazil. Zoologia
                                       35:e22235 DOI 10.3897/zoologia.35.e22235.
                                    Dahlke FT, Wohlrab S, Butzin M, Pörtner HO. 2020. Thermal bottlenecks in the life cycle define
                                      climate vulnerability of fish. Science 369(6499):65–70 DOI 10.1126/science.aaz3658.
                                    Diaz RJ, Rosenberg R. 2008. Spreading dead zones and consequences for marine ecosystems.
                                      Science 321(5891):926–929 DOI 10.1126/science.1156401.
                                    Ferreira GVB, Barletta M, Lima ARA, Dantas D, Justino AKS, Costa MF. 2016. Plastic debris
                                       contamination in the life cycle of Acoupa weakfish (Cynoscion acoupa) in a tropical estuary.
                                       ICES Journal of Marine Sciences 73(10):2695–2707 DOI 10.1093/icesjms/fsw108.
                                    Ficke AD, Myrick CA, Hansen LJ. 2007. Potential impacts of global climate change on freshwater
                                       fisheries. Reviews in Fish Biology and Fisheries 17(4):581–613 DOI 10.1007/s11160-007-9059-5.
                                    Free CM, Thorson JT, Pinsky ML, Oken KL, Wiedenmann J, Jensen O. 2019. Impacts of
                                       historical warming on marine fisheries production. Science 363:979–983
                                       DOI 10.1126/science.aau1758.
                                    Freitas CEC, Rivas AA, Campos CP, Sant’ana I, Kahn JR, Correa MAA. et al. 2013. The potential
                                       impacts of global climatic changes and dams on Amazonian fish and their fisheries.
                                       In: Turker H, ed. New Advances and Contributions to Fish Biology. London: IntechOpen,
                                       175–195.
                                    Freitas MO, Rocha GRA, Chaves PDTDC, De Moura RL. 2014. Reproductive biology of the lane
                                       snapper, Lutjanus synagris, and recommendations for its management on the Abrolhos Shelf,
                                       Brazil. Journal of the Marine Biological Association of the United Kingdon 94(8):1711–1720
                                       DOI 10.1017/S0025315414001088.
                                    Garcia AM, Vieira JP, Winemiller KO, Grimm AM. 2004. Comparison of the 1982–1983 and
                                       1997–1998 El Nino effects on the shallow-water fish assemblage of the Patos Lagoon estuary,
                                       Brazil. Estuaries 27(6):905–914 DOI 10.1007/BF02803417.
                                    Gislason H, Daan N, Rice JC, Pope JG. 2010. Size, growth, temperature and the natural mortality
                                       of marine fish. Fish and Fisheries 11(2):149–158 DOI 10.1111/j.1467-2979.2009.00350.x.
                                    Gonçalves C, Batista VS. 2008. Avaliação do desembarque pesqueiro efetuado em Manacapuru,
                                       Amazonas, Brasil. Acta Amazonica 38(1):135–144 DOI 10.1590/S0044-59672008000100015.
                                    Grimes CB. 1987. Reproductive biology of the Lutjanidae: a review. In: Polovina JJ, Ralston S, eds.
                                       Tropical Snappers and Groupers: Biology and Fisheries Management. Boulder: Westview Press,
                                       239–294.
                                    Gruber N. 2011. Warming up, turning sour, losing breath: ocean biogeochemistry under global
                                      change. Philosophical Transactions of the Royal Society A: Mathematical, Physical and
                                      Engineering Sciences 369(1943):1980–1996 DOI 10.1098/rsta.2011.0003.
                                    Gutiérrez-Yurrita PJ. 2014. Fuzzy logic applied into a holistic model to manage seasonal basins in
                                      México. EnviroGeoChimica Acta 1(6):352–374.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                             18/22
                                    Inomata SO, Freitas CEC. 2015. Fishing in the Middle rio Negro: socioeconomic aspects and
                                      operational structure. Boletim do Instituto de Pesca 41(1):79–87.
                                    IPCC. 2014. In: Core Writing Team, Pachauri RK, Meyer LA, ed. Synthesis report. Contribution of
                                      working groups I, II and III to the fifth assessment report of the intergovernmental panel on
                                      climate change. Geneva, Switzerland: IPCC, 151.
                                    Isaac VJ, Almeida MC, Cruz REA, Nunes LG. 2015. Artisanal fisheries of the Xingu River basin in
                                       Brazilian Amazon. Brazilian Journal of Biology 75(3 suppl 1):125–137
                                       DOI 10.1590/1519-6984.00314BM.
                                    Jeppesen E, Meerhoff M, Holmgren K, Gonzalez-Bergonzoni I, Teixeira-de Mello F,
                                       Declerck SAJ, De Meester L, Søndergaard M, Lauridsen TL, Bjerring R, Conde-Porcuna JM,
                                       Mazzeo N, Iglesias C, Reizenstein M, Malmquist HJ, Liu ZW, Balayla D, Lazzaro X. 2014.
                                       Impacts of climate warming on lake fish community structure and potential effects on ecosystem
                                       function. Hydrobiologia 646(1):73–90 DOI 10.1007/s10750-010-0171-5.
                                    Jepsen DB, Winemiller KO, Taphorn DC. 1999. Age structure and growth of peacock cichlids
                                       from rivers and reservoirs of Venezuela. Journal of Fish Biology 55(2):433–450
                                       DOI 10.1111/j.1095-8649.1999.tb00689.x.
                                    Kenchington TJ. 2014. Natural mortality estimators for information-limited fisheries. Fish and
                                      Fisheries 15(4):533–562 DOI 10.1111/faf.12027.
                                    King JR, McFarlane GA. 2003. Marine fish life history strategies: applications to fishery
                                      management. Fisheries Management and Ecology 10(4):249–264
                                      DOI 10.1046/j.1365-2400.2003.00359.x.
                                    Kullander SO, Ferreira EJG. 2006. A review of the South American cichlid genus Cichla, with
                                      descriptions of nine new species (Teleostei: Cichlidae). Ichthyological Explorations of Freshwater
                                      17(4):289–398.
                                    Kurtz FW, Matsuura Y. 2001. Food and feeding ecology of Brazilian sardine (Sardinella
                                      brasiliensis) from the southeastern Brazilian Bight. Brazilian Journal of Oceanography
                                      49(1–2):61–74 DOI 10.1590/S1413-77392001000100006.
                                    Levangie PEL, Blanchfield PJ, Hutchings JA. 2021. The influence of ocean warming on the
                                      natural mortality of marine fishes. Environmental Biology of Fishes 105(10):1447–1461
                                      DOI 10.1007/s10641-021-01161-0.
                                    MacCall A. 2011. What is the precision of M? In: Brodziak J, Ianelli J, Lorenzen K, Methot RD, eds.
                                      Estimating Natural Mortality in Stock Assessment Applications. Seattle: Alaska Fisheries Science
                                      Center, 21–22.
                                    Malabarba MCSL. 2004. Revision of the Neotropical genus Triportheus Cope, 1872
                                      (Characiformes: Characidae). Neotropical Ichthyology 2(4):167–204
                                      DOI 10.1590/S1679-62252004000400001.
                                    Manooch CS III, Mason DL. 1984. Age, growth and mortality of lane snapper from Southern
                                      Florida. Northeast Gulf Science 7(1):109–115 DOI 10.18785/negs.0701.08.
                                    Marengo JA, Jones R, Alves LM, Valverde MC. 2009. Future change of temperature and
                                      precipitation extremes in South America as derived from the PRECIS regional climate modeling
                                      system. International Journal of Climatology 29(15):2241–2255 DOI 10.1002/joc.1863.
                                    Matsuura Y. 1998. Brazilian sardine (Sardinella brasiliensis) spawning in the southeast Brazilian
                                      Bight over the period 1976–1993. Brazilian Journal of Oceanography 46(1):33–43
                                      DOI 10.1590/S1413-77391998000100003.
                                    Meerhoff M, Teixeira-de Mello F, Kruk C, Alonso C, González-Bergonzoni I, Pacheco JP,
                                      Lacerot G, Arim M, Beklioğlu M, Brucet S, Goyenola G, Iglesias C, Mazzeo N, Kosten S,
                                      Jeppesen E. 2012. Environmental warming in shallow lakes: a review of potential changes in
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                              19/22
                                       community structure as evidenced from space-for-time substitution approaches. Advances in
                                       Ecological Research 46(2):259–349 DOI 10.1016/B978-0-12-396992-7.00004-6.
                                    Moore RE. 2009. Introduction to interval analysis. Philadelphia: Society for Industrial and Applied
                                       Mathematics.
                                    National Water Agency of Brazil (ANA). 2020. Sistema de informações sobre recursos hídricos/
                                       Agência Nacional de Águas. Brasília: ANA.
                                    Oliveira CD, Lessa R, Almeida Z, Santana FM. 2020. Biology and fishery of Acoupa Weakfish
                                       Cynoscion acoupa (Lacepède, 1801): a review. Neotropical Biology and Conservation
                                       15(3):333–349 DOI 10.3897/neotropical.15.e55563.
                                    Oremus KL, Bone J, Costello C, García Molinos J, Lee A, Mangin T, Salzman J. 2020.
                                      Governance challenges for tropical nations losing fish species due to climate change. Nature
                                      Sustainability 3(4):277–280 DOI 10.1038/s41893-020-0476-y.
                                    Paiva MP. 1997. Recursos pesqueiros estuarinos e marinhos do Brasil. Fortaleza: Edições UFC.
                                    Paiva MP, Motta PCS. 2000. Cardumes da sardinha-verdadeira, Sardinella brasiliensis
                                      (Steindachner), em águas costeiras do estado do Rio de Janeiro, Brasil. Revista Brasileira de
                                      Zoologia 17(2):339–346 DOI 10.1590/S0101-81752000000200004.
                                    Paterson B, Jarre A, Moloney CL, Fairweather TP, van der Lingen CD, Shannon LJ, Field JG.
                                      2007. A fuzzy-logic tool for multi-criteria decision making in fisheries: the case of the South
                                      African pelagic fishery. Marine and Freshwater Research 58(11):1056–1068
                                      DOI 10.1071/MF07060.
                                    Pauly D. 1980. On the interrelationships between natural mortality, growth parameters, and mean
                                      environmental temperature in 175 fish stocks. ICES Journal of Marine Sciences 39(2):175–192
                                      DOI 10.1093/icesjms/39.2.175.
                                    Petrere Jr M. 1978. Pesca e esforço de pesca no estado do Amazonas. II – Locais e aparelhos de
                                       captura e estatística de desembarque. Acta Amazonica 8(2):1–54
                                       DOI 10.1590/1809-439219780832s005.
                                    Pianka ER. 1970. On r- and K-Selection. The American Naturalist 104(940):592–597
                                       DOI 10.1086/282697.
                                    Pimentel C, Joyeux JC. 2010. Diet and food partitioning between juveniles of mutton Lutjanus
                                       analis, dog Lutjanus jocu and lane Lutjanus synagris snappers (Perciformes: Lutjanidae) in a
                                       mangrove fringed estuarine environment. Journal of Fish Biology 76(10):2299–2317
                                       DOI 10.1111/j.1095-8649.2010.02586.x.
                                    Prestes L, Soares MGM, Silva FR, Bittencourt MM. 2010. Dinâmica populacional de Triportheus
                                       albus, T. angulatus e T. auritus (CHARACIFORMES: CHARACIDAE) em lagos da Amazônia
                                       Central. Biota Neotropica 10(3):177–181 DOI 10.1590/S1676-06032010000300020.
                                    Punt AE, Castillo-Jordán C, Hamel OS, Cope JM, Maunder MN, Ianelli JN. 2021. Consequences
                                      of error in natural mortality and its estimation in stock assessment models. Fisheries Research
                                      233:105759 DOI 10.1016/j.fishres.2020.105759.
                                    Pörtner H-O, Roberts DC, Adams H, Adelekan I, Adler C, Adrian R, Aldunce P, Ali E,
                                      Ara Begum R, Bednar-Friedl B, Bezner Kerr R, Biesbroek R, Birkmann J, Bowen K,
                                      Caretta MA, Carnicer J, Castellanos E, Cheong TS, Chow W, Cissé G, Clayton S,
                                      Constable A, Cooley SR, Costello MJ, Craig M, Cramer W, Dawson R, Dodman D, Efitre J,
                                      Garschagen M, Gilmore EA, Glavovic BC, Gutzler D, Haasnoot M, Harper S, Hasegawa T,
                                      Hayward B, Hicke JA, Hirabayashi Y, Huang C, Kalaba K, Kiessling W, Kitoh A, Lasco R,
                                      Lawrence J, Lemos MF, Lempert R, Lennard C, Ley D, Lissner T, Liu Q, Liwenga E,
                                      Lluch-Cota S, Löschke S, Lucatello S, Luo Y, Mackey B, Mintenbeck K, Mirzabaev A,
                                      Möller V, Moncassim Vale M, Morecroft MD, Mortsch L, Mukherji A, Mustonen T,
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                             20/22
                                       Mycoo M, Nalau J, New M, Okem A, Ometto JP, O’Neill B, Pandey R, Parmesan C,
                                       Pelling M, Pinho PF, Pinnegar J, Poloczanska ES, Prakash A, Preston B, Racault MF,
                                       Reckien D, Revi A, Rose SK, Schipper ELF, Schmidt DN, Schoeman D, Shaw R, Simpson NP,
                                       Singh C, Solecki W, Stringer L, Totin E, Trisos CH, Trisurat Y, van Aalst M, Viner D,
                                       Wairiu M, Warren R, Wester P, Wrathall D, Zaiton Ibrahim Z. 2022. Technical summary.
                                       In: Pörtner HO, Roberts DC, Tignor M, Poloczanska ES, Mintenbeck K, Alegría A, Craig M,
                                       Langsdorf S, Löschke S, Möller V, Okem B, Rama B, eds. Climate Change 2022: Impacts,
                                       Adaptation and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report
                                       of the Intergovernmental Panel on Climate Change. Cambridge, UK and New York, NY, USA:
                                       Cambridge University Press, 37–118.
                                    Rezende SM, Ferreira BP, Frédou T. 2003. A pesca de lutjanídeos no Nordeste do Brasil: histórico
                                      das pescarias, características das espécies e relevância para o manejo. Boletim Técnico Científico
                                      da CEPENE 11(1):257–270.
                                    Ruffino ML, Soares ECS, Lopes-Junior U, Estupinãn G, Fonseca SN, Pinto WH, Oliveira C,
                                      Barthem RB, Batista V. 2006. Estatística Pesqueira do Amazonas e Pará 2003. Manaus: Ibama/
                                      ProVárzea.
                                    Röpke CP, Amadio S, Zuanon J, Ferreira EJ, De Deus CP, Pires TH, Winemiller K. 2017.
                                       Simultaneous abrupt shifts in hydrology and fish assemblage structure in a floodplain lake in the
                                       central Amazon. Scientific Reports 7(1):40170 DOI 10.1038/srep40170.
                                    Santos GM, Ferreira E, Zuanon J. 2006. Peixes comerciais de Manaus. Manaus: Ibama/AM, Pró-
                                       Várzea.
                                    Silva EA, Stewart DJ. 2017. Reproduction, feeding and migration patterns of Prochilodus nigricans
                                       (Characiformes: Prochilodontidae) in northeastern Ecuador. Neotropical Ichthyology
                                       15(3):e160171 DOI 10.1590/1982-0224-20160171.
                                    Soares MGM, Costa EL, Siqueira-Souza FK, Anjos HDB, Yamamoto KC, Freitas CEC. 2007.
                                      Peixes de lagos do Médio rio Solimões. Manaus: EDUA.
                                    Souza SD, Bassanezi RC. 2018. Modelo de competição de Gause com esforço de pesca e fatores
                                      abióticos. Biomatemática 28:87–100.
                                    Souza SD, Prata RAC. 2019. Interatividade fuzzy no modelo generalizado de von Bertalanffy.
                                      Biomatemática 29:49–58.
                                    Sparre P, Venema SC. 2007. Introdução a avaliação de mananciais de peixes tropicais. Roma:
                                      FAO-Organização das Nações Unidas para a Alimentação e a Agricultura.
                                    Sunday J. 2020. When do fish succumb to heat? Science 369(6499):35–36
                                      DOI 10.1126/science.abd1272.
                                    Sylaios GK, Koutroumanidis T, Tsikliras AC. 2010. Ranking and classification of fishing areas
                                       using fuzzy models and techniques. Fisheries Management and Ecology 17(3):240–253
                                       DOI 10.1111/j.1365-2400.2009.00714.x.
                                    Szpilman M. 2011. Peixes marinhos do Brasil: guia prático de identificação. Rio de Janeiro: MAUA
                                       Editora Ltda.
                                    Sæther BE, Bakke O. 2000. Avian life history variation and contribution of demographic traits to
                                      the population growth rate. Ecology 81(3):642–653
                                      DOI 10.1890/0012-9658(2000)081[0642:ALHVAC]2.0.CO;2.
                                    Teh LCL, Teh LSL. 2011. A fuzzy logic approach to marine spatial management. Environmental
                                      Management 47(4):536–545 DOI 10.1007/s00267-011-9645-0.
                                    Val AL, Almeida-Val VMF. 1995. Fishes of the Amazon and their environments. Physiological and
                                      biochemical features. Heidelberg: Springer Verlag.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                             21/22
                                    von Bertalanffy L. 1938. A quantitative theory of organic growth (inquiries on growth laws II).
                                      Human Biology 10(2):181–213.
                                    Wang HY, Shen SF, Chen YS, Kiang YK, Heino M. 2020. Life histories determine divergent
                                      population trends for fishes under climate warming. Nature Communications 11(1):4088
                                      DOI 10.1038/s41467-020-17937-4.
                                    Winemiller KO. 1989. Patterns of variation in life history among South American fishes in
                                      seasonal environments. Oecologia 81(2):225–241 DOI 10.1007/BF00379810.
                                    Winemiller KO, Rose KA. 1992. Patterns of life-history diversification in North American fishes:
                                      implications for population regulation. Canadian Journal of Fisheries and Aquatic Sciences
                                      49(10):2196–2218 DOI 10.1139/f92-242.
                                    Winemiller KO, Taphorn DC, Barbarino-Duque A. 1997. Ecology of Cichla (Cichlidae) in two
                                      blackwater rivers of Southern Venezuela. Copeia 4:690–696 DOI 10.2307/1447287.
                                    Wise L, Murta AG, Carvalho JP, Mesquita M. 2012. Qualitative modelling of fishermen’s
                                      behaviour in a pelagic fishery. Ecological Modelling 228:112–122
                                      DOI 10.1016/j.ecolmodel.2011.12.008.
                                    Woodward G, Bonada N, Brown LE, Death RG, Durance I, Gray C, Hladyz S, Ledger ME,
                                     Milner AM, Ormerod SJ, Thompson RM, Pawar S. 2016. The effects of climatic fluctuations
                                     and extreme events on running water ecosystems. Philosophical Translation of the Royal Society
                                     of London B 371(20150274):1–15 DOI 10.1098/rstb.2015.0274.
                                    Yamamoto KC, Soares MGM, Freitas CEC. 2004. Alimentação de Triportheus angulatus (Spix &
                                      Agassiz, 1829) no lago Camaleão, Manaus, AM, Brasil. Acta Amazonica 34(4):653–659
                                      DOI 10.1590/S0044-59672004000400017.
                                    Zacarkim CE, Piana PA, Baumgartner G, Aranha JMR. 2015. The panorama of artisanal fisheries
                                      of the Araguaia River, Brazil. Fisheries Sciences 81(3):409–416 DOI 10.1007/s12562-015-0853-z.
                                    Zadeh LA. 1965. Fuzzy sets. Information and Control 8(3):338–353
                                      DOI 10.1016/S0019-9958(65)90241-X.
                                    Zaret TM. 1980. Life history and growth relationships of Cichla ocellaris, a predatory South
                                      América Ciclids. Biotropica 12(2):144–157 DOI 10.2307/2387730.
Campos et al. (2023), PeerJ, DOI 10.7717/peerj.14989                                                                          22/22