Combining ecosystem indicators and life cycle assessment for environmental assessment of demersal trawling in Tunisia

Purpose The present study assesses environmental performance of seafood production by demersal trawling in Tunisia (Gulf of Gabes) in order to analyze the contribution of each production stage to environmental impacts and to understand drivers of the impacts using life cycle assessment (LCA). Then a set of ecosystem quality indicators were determined using an ecosystem modeling tool, Ecopath with Ecosim (EwE), and were combined with LCA to increase the relevance of both tools ’ assessments when applied to fisheries. Methods The approach consisted of conducting LCA and calculating ecosystem indicators to provide a complete assessment of trawling ’ s environmental impacts and the ecosystem characteristics associated with seafood production. The functional unit for the LCAwas set to 1 t of landed seafood, and system boundaries included several operational stages related to demersal trawling. Several ecosystem indicators from EwE were calculated. Demersal trawling in the exploited ecosystem of the Gulf of Gabes (southern Tunisia) was used as a case study to illustrate the applicability of the approach. Several management plans were simulated and their influence on environmental performance was assessed. Ecospace, the spatial module of EwE, was used to simulate management scenarios: establishment of marine protected areas, extension of the biological rest period, and decrease in the number of demersal trawlers. Results and discussion LCA revealed that fuel consumption by fishing vessels, fuel production, and paint and antifouling production contributed most to environmental impacts. All management plans simulated decreased environmental impacts compared with the baseline scenario. The most effective management plan is extending the rest period, which increases demersal trawler yield and greatly decreases the PPR/catch of demersal trawlers. Conclusions The method developed in this study is relevant for supplementing LCA of fisheries and potentially that of seafood production systems. It provides policy makers with practical information to help implement effective management plans in the context of an ecosystem approach to


Introduction
Seafood products represent more than 9% of the economic value of total agricultural exports and 1% of the world merchandise trade (FAO 2016). The seafood trade has substantially expanded over the past few decades and is fueled by the increase in world demand for seafood, from 9.9 kg per capita in the 1960s to 19.7 kg in 2013, with preliminary estimates rising above 20 kg per capita in 2015 (FAO 2016). Seafood represents about 17% of the global population's intake of animal protein and 6.7% of all protein consumed (FAO 2016). Fishing activity has expanded considerably in terms of fishing effort and the number of fishing units and fishing grounds, because of the increased demand and consumption of seafood and the technological development in fishing technologies, in addition to the fact that fisheries operate in an increasingly globalized environment (Pauly et al. 2002;Swartz et al. 2010). On a global scale, marine catches have remained relatively stagnant since 1996 (FAO 2016). Compared with agriculture, the global seafood economic system is relatively small; however, the increased supply of seafood carries the risk of ecological degradation of ecosystems and has considerable impact on the world's environment (Kaiser and de Groot 2000). Therefore, it is imperative to apply sustainability principles to seafood production systems.
Long-term sustainability of fishing is a major concern from an environmental and ecological viewpoint. Fishing activity carries the risk of negative direct and indirect impacts on marine ecosystems (Kaiser and de Groot 2000), mainly because it relies entirely on extracting organisms from the ecosystem (Christensen et al. 2003). Consequently, environmental impacts of seafood production have been intensively studied in recent decades (World Bank 2017; Worm et al. 2009). However, most environmental analysis focused on the immediate impact of fisheries on targeted stocks (Costello et al. 2016), by-catch and discards (Glass 2000), benthic communities (Guyonnet et al. 2008), seabed damage (Kaiser et al. 2006) and changes in trophic dynamics, and the structure and functioning of the ecosystem (Jackson et al. 2001;Pauly et al. 2002;Tremblay-Boyer et al. 2011).
An integrated, science-based approach is important for impact assessment, and life cycle assessment (LCA) is an effective method that considers the entire supply chain to estimate potential environmental impacts associated with a production or service, such as seafood production (Pelletier et al. 2007). Since the 2000s, LCA has been applied to fisheries and seafood production systems around the world (Ziegler et al. , 2011Hospido and Tyedmers 2005;Iribarren et al. 2010;Vázquez-Rowe et al. 2010Ramos et al. 2011;Fréon et al. 2014;Avadí et al. 2015;Abdou et al. 2017Abdou et al. , 2018b. Although LCA is an extensive approach, given the range of impacts it is able to assess, authors of most seafood LCA studies have strongly recommended including ecosystem components (e.g., biotic resources, primary production, trophic interactions) to supplement the impact assessment (Avadí and Fréon 2013), especially with the current overexploitation of marine resources and increased disturbances in marine ecosystems caused by human activities (Halpern et al. 2008). Several LCA studies included and discussed direct fishery-specific impacts, most of which were calculated outside the LCA framework (Avadí and Fréon 2013). For example, a biotic resource use impact category was developed and applied to quantify the amount of biotic biomass removed by fishing activity (Boyd 2008;Driscoll et al. 2015;Langlois et al. 2014a;Parker and Tyedmers 2012a, b). Biotic resource use includes estimates of the primary production required (PPR). The sea-use impact category was also introduced in LCA to represent physical impacts due to occupation or transformation of marine areas (Langlois et al. 2014b). Nilsson and Ziegler (2007), Ramos et al. (2011), andHansson (2003) applied seafloor area impact to different fisheries. The last potential yield impact category was developed to quantify overfishing and the depletion of exploited fish stocks (Emanuelsson et al. 2014). Specific discard indexes in LCA were developed to characterize and standardize estimates of discards in fisheries (Vázquez-Rowe et al. 2012b).
Ecosystem models are used in ecosystem-based fishery management to assess the sustainability, health, productivity, and resilience of ecosystems. This concept is used to incorporate ecosystem considerations for the sustainable utilization of marine resources and not just single-species management plans. Ecopath with Ecosim (EwE, Christensen and Walters 2004) is one of the most frequently used models in the world to model marine and aquatic ecosystems (Plagányi 2007). The model provides better understanding of impacts of fishing on target and non-target species and assesses interactions between ecosystem components. The EwE approach explicitly considers trophic interactions and helps in studying fishing activities within an ecosystem context. EwE has three main modules: (i) Ecopath, a static mass-balanced snapshot of the system; (ii) Ecosim, a dynamic simulation module; and (iii) Ecospace, a spatial and temporal dynamic module.
EwE provides valuable information about ecosystem functioning; however, it considers only extraction of organisms from the ecosystem by fisheries and overlooks relevant aspects related to the performance of fishing. For instance, fuel consumption was found to be an important driver behind environmental impact of fisheries (Thrane 2004;Tyedmers et al. 2005). The model does not consider effects of the energy and material used to construct and maintain fishing vessels and gear (Hayman et al. 2000;) and the use of paint and antifouling paint (Hospido and Tyedmers 2005). LCA considers these aspects, which are crucial to establish a robust and complete impact assessment of fishing.

Study area: the Gulf of Gabes
The Gulf of Gabes is located in southern Tunisia and in the southern Mediterranean Sea and it covers approximately 35,900 km 2 (Fig. 1). The gulf is characterized by its wide continental shelf (a depth of 200 m is not reached until 400 km offshore), resulting in high sensitivity to atmospheric changes (Natale et al. 2006). It has the highest tidal amplitude in the Mediterranean Sea, reaching 1.8 m in height (Sammari et al. 2006). The gulf contains a large bed of Posidonia oceanica, an endemic Mediterranean seagrass (Batisse and Jeudy de Grissac 1998) that provides an important nursery, feeding, and breeding ground for many marine species (Hattour 1991). The Gulf of Gabes is under multiple natural and anthropogenic threats Lamon et al. 2014). It is a major fishing ground in Tunisia due to its richness in benthic fauna (i.e., shrimp, mullets, soles) and the presence of soft bottom habitats that facilitate access to resources (Missaoui et al. 2000). It is considered one of the most productive areas in the Mediterranean Sea in terms of catches Papaconstantinou and Farrugio 2000). The Gulf of Gabes supports 60% of fleets in the country and provides more than 40% (more than 40,000 t) of the annual national fish production (DGPA 2015). Catches are dominated by Sparidae (e.g., Diplodus annularis), mullet (Mullus barbatus and M. surmuletus), round sardinella (Sardinella aurita), European pilchard (Sardina pilchardus), and several benthic cephalopods (e.g. Sepia officinalis, Octopus vulgaris). Demersal trawling is the predominant fishing activity in the gulf Mosbah et al. 2013). Based on Tunisian Fisheries and Aquaculture Department (DGPA) statistics, the number of demersal trawlers increased to 285 from 1980 to 1991 and decreased to 229 in 1997. The number of trawlers has fluctuated around 250 since 1997 (DGPA 2015). In 2015, 226 demersal trawlers were operating in the gulf, and most (184) were wooden (DGPA, 2015) targeting shrimp and demersal finfish (Sparidae (D. annularis, Sparus aurata), mullets, rays (Raja clavata), and sharks (Mustelus mustelus). Demersal trawlers in the Gulf of Gabes produced 10,332 t of seafood in 2015, which generated approximately 41 million € (DGPA, 2015); however, it is considered the most destructive fishing gear worldwide. Fuel consumption by demersal trawlers in the Gulf of Gabes varies between 35 and 180 metric tons of diesel per year, depending on trawler's size, number and length of fishing trips, main species targeted, and other technical Fig. 1 Geographic location of the study area in the Gulf of Gabes ecosystem characteristics. It damages bottom habitats and harms benthic communities, in addition to its non-selectivity (Kumar and Deepthi 2006). The Gulf of Gabes is an archetypal ecosystem in which the effects of fisheries are the most pronounced, and according to stock assessments, is considered a highly exploited ecosystem . Indeed, the expansion of fisheries resulted in overexploitation of several stocks (Fiorentino et al. 2008). One sign of overexploitation in the Gulf of Gabes is a decrease in hourly yield from 75 kg h −1 in the 1970s to 30 kg h −1 in the 1990s. Therefore, it is necessary to establish adequate management practices to facilitate the recovery of marine resources. To address these issues, the EwE model was applied to the Gulf of Gabes. The massbalance model was developed to better understand ecosystem structure and functioning and to study impacts of fisheries ). In addition, the Ecospace module was used to investigate potential ecosystem responses to spatial and temporal management plans . LCA of demersal trawlers and EwE in the Gulf of Gabes were developed separately and published previously (Abdou et al. 2018a;Hattab et al. 2013). Combining EwE with fish supplychain models has been previously attempted . The ultimate goal of this study is to develop a new set of potential indexes from EwE to supplement seafood LCA studies and place them in an ecosystem context. The method was applied to the Gulf of Gabes ecosystem. In addition, EwE was used to assess potential ecosystem responses to different management scenarios.

Life cycle assessment
The present study was based on a seafood LCA of demersal trawlers in Tunisia as case study (Abdou et al. 2018a) which quantified environmental impacts associated with landing seafood caught by wooden demersal trawlers in the Gulf of Gabes. It was performed following the four steps suggested by the International Reference Life Cycle Data System (European Commission 2010). The functional unit (FU) was set as 1 t of landed seafood. Of the 226 demersal trawlers in the Gulf of Gabes, the 184 wooden vessels were selected for this study because of the availability of reliable data for this group. The studied system covers a wide range of operational stages related to demersal trawling such as fuel consumption by fishing vessels, other inputs to fishing (reflecting onboard activities, release of paint and antifouling substances, etc.), fuel production, lubricating oil production, paint and antifouling paint production, trawler and trawl net construction and maintenance, and transport of materials upstream. However, end-of-life stages (e.g., material recycling, disposal of certain materials) and post-landing stages (e.g., sorting and packaging, sale, use, and processing) were excluded from the analysis because those stages occur beyond the scope of fishery management decisions as examined in study. Although refrigerants have been shown to contribute heavily to certain environmental impact of fisheries, they were not included in this study because they are not used by any of the wooden demersal trawlers in Tunisia. Therefore, the LCA is considered as a "cradle-to-gate" because the final product is 1 t of seafood landed at the port gate. The main objective of the LCA was to analyze the contribution of each production stage to environmental impacts and to understand drivers of the impacts associated with demersal trawling. Thus, impacts were not allocated between coproducts.
The LCI was based on data collected from official DGPA records (e.g., landings, length, tonnage, engine type, engine power, lightship weight). Additional data were obtained from surveys of demersal trawler skippers and fishermen in the port of Sfax, the main fishing port in the Gulf of Gabes, and of demersal trawler builders in the shipyard. Gathered data included trawlers' operational aspects (e.g., fuel consumption, number of fishing trips, and number of days at sea) and information about vessel construction (e.g., material used for construction, paint and antifouling paint quantities, dimensions of vessels, life span). Collected data supplement and validate official statistics.
Mean lifespan of wooden vessels in the Gulf of Gabes is estimated at 40 years and maintenance and repair of vessels are taken into account by adding 25% of the total amount of wood used for vessel construction (Tyedmers 2000). The composition of paint and antifouling paint was obtained from material safety data sheets and validated by managers of the principal paint-producing company in Tunisia, and their emissions to water were estimated based on the assumption that two-thirds of the paint and antifouling paint used is released into the water (Hospido and Tyedmers 2005). Required background data were extracted from the ecoinvent 3.0 database (Weidema et al. 2013). See Abdou et al. (2018a) for further details about the method.
To estimate impacts, data were aggregated into 11 midpoint impact categories that reflect the environmental issues associated to fisheries and seafood production system. The selection was based on previous guidelines in fishery LCA studies (Avadí and Fréon 2013). Impacts were calculated per ton of seafood produced using SimaPro® 8.0 software (Goedkoop et al. 2008). The following impact categories were included: & Abiotic depletion potential: reflects the decrease in nonrenewable and renewable abiotic resources available for human use; expressed in kg Sb equivalent (eq). & Acidification potential: reflects negative acidic effects to the atmosphere generated by production; expressed in kg SO 2 eq. & Eutrophication potential: represents negative effects of discharging nitrogen and phosphorus into the environment; expressed in kg PO 4 eq. & Global warming potential: represents potential contribution to radiative forcing via greenhouse gas emissions; expressed in kg CO 2 eq. & Ozone depletion potential: represents potential damage to the ozone layer caused by chlorinated and brominated chemicals; expressed in kg CFC-11. & Photochemical oxidant formation potential: represents negative effects of chemical substances caused by sunlight reacting with emissions from reactive substances (fossil fuel); expressed in kg ethylene (C 2 H 4 ) eq. & Human toxicity potential: represents potential harm of a unit of chemical released into water, air or soil; expressed in kg 1,4-dichlorobenzene (1,4-DCB) eq. & Marine ecotoxicity potential: represents negative effects of toxic substances on marine ecosystems; expressed in kg 1,4-DCB eq. & Terrestrial ecotoxicity potential: represents the negative impacts of toxic substances on terrestrial ecosystems; expressed in kg 1,4-DCB eq. & Land occupation potential: represents the land area necessary to produce the FU (1 t of seafood); expressed in square meters per year (m 2 year −1 ). & Total cumulative energy demand: represents the amount of energy (e.g. fossil fuels, electricity) necessary to produce the FU (1 t of seafood); expressed in megajoules (MJ).
A full description of the LCI and impact categories is found in Abdou et al. (2018a). All impact categories were calculated for an average demersal trawler in the Gulf of Gabes.

Ecopath with Ecosim modeling
Principles, basic concepts, and assumptions of the EwE modeling approach are described in detail in Christensen et al. (2008); Christensen and Walters (2004); Walters et al. (1997). EwE has three main components: & Ecopath: a static snapshot of interactions among functional groups in an ecosystem (Christensen and Walters 2004). Ecopath is built from two linear equations: one for mass balance and one for energy balance. & Ecosim: a dynamic simulation model based on balanced Ecopath parameters. It enables exploring effects of fishing options and changes in ecosystem functioning (Christensen and Walters 2004). & Ecospace: a spatially explicit dynamic module of EwE (www.ecopath.org; Christensen and Walters 2004;Walters et al. 1999). It integrates trophic and temporal dynamics of Ecopath and Ecosim in a two-dimensional space. After spatial grid cells are defined, each cell is assigned a habitat type, which has a relative primary production (Christensen et al. 2008). Each functional group is assigned to its preferred habitat type, and each type of fishery is assigned to its allowed fishing zones (Walters et al. 1999).
Ecopath includes indicators that describe the ecosystem based on information about the food web, trophic flow, thermodynamic concept, information theory, and network analysis (Christensen and Walters 2004;Coll et al. 2006). Fishing activities can influence the maturity, stability, and complexity of ecosystems on several levels. Indictors in EwE models indicate the state of the ecosystem and how it changes over time (Christensen and Walters 2004). In the present study, the following indicators from the EwE model were chosen to supplement LCA to provide a complete assessment of trawling's environmental impacts and place its results in an ecosystem context:

Pressure indicators
& PPR: represents the primary production required to sustain the catches and the consumption by the trophic groups. Ecopath estimates PPR, by removing all cycles from the diet composition (DC') and identifying all paths in the flow network using the method of Ulanowicz (1995). PPR of yield Y of a given group of species is quantified by summing all pathways leading to the group: the consumption:production ratio. PPR, which is equivalent to net primary production, enables ecosystems, including terrestrial ecosystems to be compared. Many authors have developed methods to calculate it, and the one most frequently used is (Pauly and Christensen 1995): CR is the conversion rate of wet weight to carbon (a ratio of 9:1), TE is the transfer efficiency between trophic levels and was assumed to be 10%, TL i is the trophic level of species i, and n is the number of species caught. In this study, we used the PPR calculated by EwE, and the values expressed using EwE are different than those reported in other LCA studies because of the difference in calculation methods and the units. We also calculated PPR/catch, PPR of demersal trawlers, and PPR/catches of demersal trawlers. & Yield: yields in Ecopath are expressed in t year −1 km −2 .

Exploited resources indicators
& Mean trophic index (MTI): describes direct and indirect trophic interactions among functional groups (Pauly and Watson 2005). It equals ∑TL i :B i ∑B i , where B is the biomass of species i. For this indicator to be sensitive to fishing pressure, only species with a trophic level higher than 3.25 (a standard threshold) are considered; species of low trophic level are excluded because their biomass tends to vary greatly in response to environmental factors (Pauly and Watson 2005). & Mean trophic level of the catch (TLc): reflects effects of fishing on the food web (Pauly et al. 1998). It is calculated as: TLc

Trophic chain indicators
& Mean trophic transfer efficiency (TE): the percentage of prey production that is transferred to predator production. It equals predator production divided by prey production. & Apex predator indicator (API): the percentage of predators with a trophic level higher than 4 (top predators) out of all predators, excluding planktivores (i.e., trophic level < 3.25) (Bourdaud et al. 2016). Fishing decreases API by removing large individuals of high trophic level.
This study is based on the EwE model of Hattab et al. (2013) and the Ecospace model of Abdou et al. (2016). The EwE model for the Gulf of Gabes was created with EwE version 6.2 (www.ecopath.org; Christensen and Walters 2004;Walters et al. 1999). The model includes 62 species divided into 41 functional groups based on ecological and taxonomic similarities. It includes the main fleets operating in the area: demersal trawling, small seines, tuna purse seines, purse seines using lights, coastal motorized fishing, and sponge fishing. Landing statistics were obtained from D G PA . F u r t h e r d e t a i l s o f d a t a r e s o u r c e s a n d parameterization are found in Hattab et al. (2013).
The Ecosim approach was used to carry out dynamic simulation using the parameters from the balanced Ecopath model. Landing data, time series of fishing effort by fishing gear, and stock assessment estimates were collected by the Tunisian National Institute of Sciences and Technologies of the Sea. Ecosim was fitted to landings for the period 1995-2008. Primary production in the study area for the period 1997-2007 (data from SeaWIFS project, http://oceancolor.gsfc.nasa.gov/ SeaWiFS/) was used to calibrate Ecosim ).
The Ecospace model of the Gulf of Gabes covers approximately 25,000 square cells, each covering 3.2 km 2 . Each cell represents land or water and is assigned to a specific habitat type (nine habitat types were included). Each of the 41 functional groups was assigned to its preferred habitat (i.e., where the feeding rate was high, based on the available information on the ecology and biology of each species and expert advice), and each fishery fleet was assigned to the fishing zones authorized by Tunisian fisheries regulations. The modeled area extended from the 20-m contour depth to approximately the 200-m isobaths, because of lack of reliable data. Ecospace was used mainly to assess scenarios of establishment of marine protected areas (MPA). More details on Ecospace implementation and parameterization are available in Abdou et al. (2016).

Simulations of fishing management scenarios
Ecospace was used to simulate scenarios over a 15-year period (1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010). The baseline scenario reflects the state of the Gulf of Gabes ecosystem after 15 years under the current situation, which has no MPA and a 3-month biological rest period of demersal trawlers. In addition to the baseline scenario, seven fishing management scenarios were simulated. Most scenarios were assessed previously Halouani et al. 2016) and do not reflect desired or planned management plans, but rather represent contrasting scenarios for the spatial and temporal management plans. Scenarios consisted of implementing temporal and spatial measures to explore potential ecosystem response (Fig. 2 When using Ecospace to simulate the spatial management scenarios, the model redistributes the fishing effort according to the new regulations instead of reducing it. For each simulation, the model predicts changes in functional group biomass and the changes in yield (per group and per fishery). Then, the new catch and biomass values (after a 15-year run) predicted from the Ecopsace model were used to re-estimate the different environmental impacts (using LCA) and recalculate the ecosystem indicators. All results were then compared with the baseline scenario. The EwE model was run separately and the results were used to supplement the results of the LCA for demersal trawlers to place this activity in an ecosystem context.

LCA results for the current situation
Estimated environmental impacts of an average demersal trawler under the current situation are presented in Table 1. Fuel production is responsible for around 96% of abiotic depletion potential (total = 93.9 kg Sb eq t −1 ) and ozone depletion potential (total = 0.003 kg CFC-11 eq t −1 ). It is also responsible for around 55% of photochemical oxidant formation (total = 2.4 kg C 2 H 4 eq t −1 ), marine ecotoxicity (total = 1,583,942 kg 1,4-DCB eq t −1 ), land occupation potential (total = 119.2 m 2 year t −1 ), and total cumulative energy demand (total = 212,637 MJ eq t −1 ). Fuel consumption by fishing vessels is the largest contributor to acidification (84% of 155.8 kg SO 2 eq t −1 ), global warming (81% of 13,773 kg CO2 eq t −1 ), and eutrophication (57% of 9.32 kg PO 4 eq t −1 ). Paint and antifouling production contributes less than the other processes and is responsible for about 13% of marine ecotoxicity, land occupation, and total energy demand. Construction of the trawler and trawling net is responsible for 84% of terrestrial ecotoxicity (total = 88.8 kg 1,4-DCB eq t −1 ) and human toxicity (57% of 4.094 kg 1,4-DCB eq t −1 ), and also contributes to marine toxicity, land occupation, and total energy demand (30%). Transport of materials upstream does not contribute greatly to impacts (Table 1).

EwE indicators of the current situation
The pressure indicators reveal a PPR of yield of 240.6 kg of wet weight km −2 (126.9 and 113.8 kg of wet weight km −2 from primary producers and detritus, respectively). PPR/ catch is 139.7 kg of wet weight t −1 km −2 . PPR and PPR/ catch of demersal trawlers are 313.5 kg of wet weight km −2 and 522.9 kg of wet weight t −1 km −2 , respectively. Total yield in the ecosystem is 1.72 t km −2 year −1 , and demersal trawler yield is 0.85 t km −2 year −1 . The exploited resource indicators reveal MTI of 3.85 and mean TLc of 3.44. The trophic chain indicators reveal mean trophic TE of 20.3% and API of 44% (Table 2).

Baseline scenario
In the 15-year simulation of the current situation, Ecospace predicts that total yield decreases by 11% and demersal trawler yield decreases by 37%. Total biomass of functional groups decreases by 16%. LCA impacts after the 15-year simulation (calculated using the final catch values) show a 59% increase in all impacts. PPR increases by 31% to 314.61 kg of wet weight km −2 , while PPR/catch increases by 100% to 280 kg of wet weight t −1 km −2 . PPR of demersal trawlers decreases by 44%; however, PPR/catch of demersal trawlers decreases to 177 kg of wet weight t −1 km −2 . Total yield decreases by 20%, and demersal trawler yield decreases by 60%. MTI does not change (3.87). Although TL C increases by 10%, API increases by 12% (reaching 50%) ( Table 2).

MPA scenarios
Simulated scenarios of coastal MPA establishment have lower environmental impacts than the baseline scenario. All impacts increase by 22%, compared with 59% in the baseline scenario (Fig. 3). Total PPR increases to 305.64 kg of wet weight km −2 , which is a 27% increase compared with 31% in the baseline scenario; however, PPR of demersal trawlers decreases as much as in the baseline scenario (44%). Total PPR/catch Fig. 2 Demersal trawling management plans in the Gulf of Gabes simulated with the Ecospace model. MPA = marine protected area. Orange represents the MPA, blue is the area included in the model and gray represents land increases by 102%, compared with a 100% increase in the baseline scenario. PPR/catch of demersal trawlers decreases by 1.6% compared with 0.8% in the baseline scenario. Total yield decreases by 37%, and demersal trawler yield decreases by 60%. MTI remains the same. TL C increases by 6% to 3.6 compared with 3.8, and the API increases to reach 50% instead of 44% (Fig. 4).
Simulated scenarios of offshore MPA establishment have lower environmental impacts than the baseline scenario. All impacts increase by 27% with a small offshore MPA and 26% with a large offshore MPA, compared with 59% in the baseline scenario (Fig. 3). Total PPR, PPR/catch, PPR of demersal trawlers, and PPR/catch of demersal trawlers experience similar changes as the baseline scenario for both MPAs. Total yield decreases by 34% and demersal trawler yield decreases by 60% for both MPAs. MTI does not change. TL C increases to 3.8, and API increases to 50% (same as the baseline scenario) (Fig. 4).

Rest period scenarios
According to the 15-year simulations, establishment of 4-month and 6-month rest periods decrease all environmental impacts by 7% and 15%, respectively (Fig. 3), which is better than the 59% increase in impacts in the Table 1 Contribution of operational stages to mean environmental impacts per ton of seafood produced by an average demersal trawler in the Gulf of Gabes. ADP, abiotic depletion potential; AP, acidification potential; EP, eutrophication potential; GWP, global warming potential; ODP, ozone depletion potential; POFP, photochemical oxidant formation potential; HTP, human toxicity potential; METP, marine ecotoxicity potential; TETP, terrestrial ecotoxicity potential; LOP, land occupation potential; TCED, total cumulative energy demand baseline scenario. In both cases, total PPR increases by 38%, PPR/catch increases to 85%, and PPR/catch of demersal trawlers decreases by 6%; however, PPR of demersal trawlers decreases by 21% and 15% with 4-month and 6-month rest periods, respectively. Total yield decreases by 25% (1.3 t km −2 year −1 ), and demersal trawler yield decreases by 40% and 37% with 4-month and 6-month rest periods, respectively. MTI remains the same. TL C increases to 4.3, and API increases to reach 50% (Fig. 4).

Scenarios for decreasing the number of demersal trawlers
Effects of decreasing the number of demersal trawlers vary by impact category. According to the 15-year simulations when 20% of demersal trawlers are eliminated, abiotic depletion, acidification, eutrophication, global warming, ozone layer depletion, photochemical oxidant formation, and total cumulative energy demand decrease by approximately 31% compared with the baseline scenario. On the other hand, a 50% reduction in demersal trawler number results in the decrease of those impact categories by 55% compared with the baseline scenario. Marine aquatic ecotoxicity and land occupation decrease by 46% and 66% compared with the baseline scenario, when eliminating 20% and 50% of demersal trawlers, respectively. Human toxicity decreases by 56% compared with the baseline case when 20% of demersal trawlers were eliminated and decreases by 73% compared with the baseline scenario when 50% of demersal trawlers were eliminated. Terrestrial ecotoxicity decreases by 66% and 81% compared with the baseline scenario when 20% or 50% of demersal trawlers were eliminated, respectively (Fig. 3). When 20% of demersal trawlers were eliminated, total PPR increases by 30% and PPR/catch increases by 92%. PPR and PPR/catch of demersal trawlers decreases by 70% and 3%, respectively. Total yield decreases by 33% and demersal trawler yield decreases by 79% when 20% of demersal trawlers were eliminated; however, when 50% of the trawlers were eliminated, total yield decreases by 35% and demersal trawler yield decreases by 84%. MTI remains the same. TL C increases by 16% when 20% of trawlers were eliminated and increases by 10% when 50% of trawlers were eliminated. API increases by 1.5% in both scenarios (Fig. 4).

Discussion
LCA impacts calculated for 1 t of landed seafood showed that fuel production and fuel consumption by fishing vessels contribute most to environmental impacts. This large contribution is due to the use of demersal trawling, which has the most fuelintensive gear (Avadí and Fréon 2013;Schau et al. 2009). Avadí and Fréon (2013) stated that fuel combustion and production contributes most to environmental impacts of fishing vessels, and several fishery LCAs demonstrated that fuel consumption is a major contributor to many impact categories Hospido and Tyedmers 2005;Pelletier et al. 2007;Schau et al. 2009;Vázquez-Rowe et al. 2012a). Demersal trawlers in the Gulf of Gabes consume 4841 l of diesel to produce one metric ton of seafood, which is twice the global average for bottom trawl fisheries (around 2000 l per metric ton of seafood) estimated by Parker and Tyedmers (2015) and almost 10 times higher than the global average for fisheries (489 l of diesel per metric ton of seafood) estimated by Parker et al. (2018). Paint and antifouling production contributed the most to marine ecotoxicity, and also to land occupation and total cumulative energy demand due to emissions of copper, xylene, lead, tributyltin, and zinc oxides. In many fishery-related LCAs, trawler and fishing gear construction were not considered due to its supposedly small contribution to environmental impacts (Avadí and Fréon 2013). However, other studies found large impacts of the construction phase Svanes et al. 2011). In this study, construction contributed more than 80% of terrestrial ecotoxicity, more than 50% of human toxicity, and more than 30% of marine toxicity, land occupation, and total cumulative energy demand. This can be explained by the fact that only wooden vessels were considered in this study while most other fishery LCA studies considered vessels with steel hull. Thus, it seems necessary to include the construction phase in fishery LCA studies due to its large contribution to environmental impacts, especially toxicity impacts. Transport of materials upstream contributes much less to impacts than to the other processes, which suggests that it is not a key subsystem of demersal trawling in the Gulf of Gabes.
Although LCA provides environmental assessment from "cradle-to-grave", several key impacts are still lacking. To ensure the sustainability of fishing, it is important to consider ecosystem state. EwE provides ecosystem indicators to describe ecosystem state and impacts of fisheries on the ecosystem. PPR, one of the most common pressure indicators to describe ecosystem state, is widely used in LCA studies because it provides a measure of biotic resource use. Although PPR is commonly applied to seafood products, its fundamental assumptions can be challenged. For example, PPR is calculated by assuming a TE of 10% per trophic level. However, Libralato et al. (2008) estimated variability in TE ranging from 5 to 14% depending on the type of ecosystem and fish species variability. Luong et al. (2015) showed that standard estimates of PPR are 3.9-5.0 times as high as those estimated when adopting the food chain theory. We chose to calculate PPR using EwE, which considers the entire ecosystem and interactions between its compartments and fisheries.
Official statistics and stock assessments indicate that the Gulf of Gabes is a highly fished area. Most targeted species are overexploited or fully exploited (e.g., Pomatomus saltatrix (Dhieb et al. 2007), M. surmuletus (Ben Meriem et al. 1994), M. barbatus (Gharbi et al. 2004), Pagellus erythrinus (Jarboui et al. 1998)). Indicators from EwE indicate that fishing is unsustainable in the Gulf of Gabes (high MTI and TLc). This finding is in line with what was found in Hattab et al. (2013). LCA results demonstrated that demersal trawling in the gulf has high environmental impacts in all categories. This finding was supported by a comparison with other fishery LCA studies in Abdou et al. (2018a). Consequently, it is necessary to establish fishery management measures, such as MPAs, rest periods, and a substantial reduction in fishing effort.
We used the Ecospace module as a decision support tool to assess fishery management scenarios in the context of ecosystem-based fishery management. Impacts in all categories increased by 59% after a 15-year run of the baseline scenario (no management plan). All management measures had a positive impact on the environmental performance and decreased all impacts compared with the baseline scenario. MPA scenario results were better when a small coastal MPA or a large offshore MPA was established, due to the type of habitat it covered. Previous studies demonstrated the importance of bottom characteristics in establishing a successful MPA (Guizien et al. 2012). The rest period scenarios provided the best results for environmental impacts and ecosystem indicators (except for PPR and PPR of demersal trawlers) among all the scenarios assessed. Extending the rest period decreased all impacts compared to those of the baseline scenario (Table 3). Impacts of demersal trawlers decreased as the duration of the rest period increased: 4-month and 6-month rest periods decreased impacts by 7% and 15%, respectively, compared with those in the baseline scenario. Decreasing the number of demersal trawlers by 20% or 50% had better results for terrestrial, marine, and human toxicity than the MPA scenarios. For the other impact categories, the 20% decrease provided the worst results among all scenarios (33% increase), and the 50% decrease provided better results than the MPA scenarios. Total PPR and PPR/catch did not differ greatly among scenarios and, as expected, PPR of demersal trawlers was higher than that of the baseline when the rest period was extended and was lower than that of the baseline when the number of trawlers decreased. For the seven simulated management scenarios, total yield decreased compared with that of the baseline scenario; however, demersal trawler yield increased when the rest period was extended and decreased when the number of trawlers decreased. Simulation results indicate that the management measures increased TLc, except for the establishment of MPAs. API increased with the coastal MPA and the extent of the rest period, and decreased when the number of demersal trawlers decreased. The management scenarios were intended only to help understand ecosystem functioning and are not for direct application. The scenario results would interest stakeholders because they could help to identify management scenarios that would maintain or increase landings while also providing adequate environmental performance and not compromising ecosystem structure.
The overall results indicate that EwE indicators are able to supplement LCA results to provide a complete assessment of fisheries and place them in an ecosystem-based management context. However, this study could improve on many levels. The LCA results would be more accurate if the boundaries were extended to include post-harvesting processes related to demersal trawler yield and establish a "cradle-to-grave" LCA instead of just a "cradle-to-gate" LCA. In addition, conducting LCA to include all fisheries in the EwE model may provide better insight into the fishing activity and state of the Gulf of Gabes ecosystem. Thus, it may help in developing effective management plans to ensure sustainability. The EwE model would also improve if more updated ecosystem data were available, especially due to the high uncertainty of parameters in Ecopath. The uncertainty in Ecospace and LCA should be considered when interpreting the results. We focused on trends when comparing the management scenarios to reduce the uncertainties related with data reliability and the complexity of the ecosystem.

Conclusions
This study developed a framework to supplement the LCA of seafood production with ecosystem indicators from Ecopath with Ecosim. The Gulf of Gabes was used as a case study to conduct the LCA, which was combined with the EwE model to provide a complete assessment of the environmental performance and ecosystem characteristics associated with production of seafood landed by demersal trawlers. LCA results showed that fuel consumption by fishing vessels, fuel production, and paint and antifouling production were the main contributors to environmental impacts. Ecosystem indicators from the EwE model provide valuable information to conduct environmental analysis using LCA and place it in the context of an ecosystem approach to fisheries. The Ecospace module of EwE was used to simulate management scenarios. All management scenarios decreased environmental impacts compared with those of the baseline scenario; however, ecosystem indicators varied more among scenarios. Among the scenarios, extending the rest period to 6 months is the most effective management plan, which increases total yield and demersal trawler yield compared with those of the baseline scenario. Total primary production required of demersal trawlers increased in this scenario, but PPR/catch of demersal trawlers greatly decreased. Apex predator index and trophic level of catch increased with the implementation of this measure. Results of this study provide stakeholders and policy makers with practical information that can help identify the most effective management plan, since ecosystems may respond differently to management measures depending on their characteristics.
Funding information The authors would like to acknowledge valuable financial support from the "Institut de Recherche pour le Développement" (JEAI GAMBAS project). This study was partially funded by the "ICV pêche" project and the "LabexMer" (ANR-10-LABX-19).
Open Access This article is distributed under the terms of the Creative Comm ons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.