Assessing the Fishing Impact on the Marine Ecosystem of Guishan Island in the Northeastern Waters of Taiwan Using Ecopath and Ecosim

: The northeastern waters of Guishan Island constitute one of the crucial fishing grounds for coastal trawl fishery in Taiwan and have been exploited for many decades. To construct the marine ecosystem and to examine the interactions among trophic levels of fisheries resources in the waters of Guishan Island, historical catch, catch composition, biological information, fishing effort, environmental data such as sea surface temperature, salinity


Introduction
Taiwan is located in the shared boundary of the East and the South China Sea; the complicated habitats and current systems contribute to the high biodiversity of its surrounding waters.The Guishan Island is located in the northeastern waters of Taiwan which the Kuroshio Current passes, and the surrounding upwelling brings rich nutrients.The high biodiversity of fishery resources, including shrimps, crabs, cattle fishes, squids, sea breams, and tile fishes, is due to the high primary production and rich nutrients in this area.Thus, this area is one of the most important fishing grounds for coastal trawl fishery in Taiwan.Additionally, this area is one of the best dolphin-watching areas due to the high occurrence of cetaceans.The trawl fishery is based in the Tahsi fishing port, while dolphin-watching vessels is based in Wushi port (Figure 1). the high occurrence of cetaceans.The trawl fishery is based in the Tahsi fishing port, while dolphin-watching vessels is based in Wushi port (Figure 1).According to the Taiwan Fisheries Year Book [1], the major fisheries in this area include the bottom trawl, torch light, gill net, hand line, mid-water trawl, and longline.High biodiversity is found in this region, and the major economic species include the black pomfret, Parastromateus niger, skipjack, Katsuwonus pelamis, yellowfin tuna, Thunnus albacares, sergestid shrimp, Sergia lucens, crimson seabream, Dentex tumifrons, black seabream, Spondyliosoma cantharus, bigeye, Priacanthus macracanths, hairtail, Trichiurus lepturus, pomfrets, narrowbarred mackerel, Scomberomorus commerson, cuttle fishes, squids, etc.However, due to the increase in fishing pressure in recent decades, the coupling of the catch of economic species decreased, with smaller sized fish at catch provoking coastal fishery management by the local government [2].The Fisheries Agency, Taiwan, has implemented management measures which included prohibiting trawl and gillnet fishing within 3 miles from the coast and seasonal closure for trawl fishery that reduced the fishing effort.The yield of these areas decreased from 42,990 tons in 2000 to 28,476 tons in 2008.
Conventional single species' stock assessments, such as production models [3,4] and yield per recruit models [5], did not take account of the interactions among species.Hence, several multispecies' stock assessment methods taking account of the interactions among marine species based on single species models have been developed since late 1970s [6], i.e., multispecies production models (MSP, [7]), multispecies virtual population analysis [8], and multispecies bioenergetics models (MSBE, [9][10][11]).In addition, ecosystem-based assessment packages which can simulate the interactions among species in an ecosystem such as Ecopath with Ecosim [12] and NETWRK [13] have also been developed.Many more input parameters required for these packages resulted in higher uncertainty of the output [14].According to the Taiwan Fisheries Year Book [1], the major fisheries in this area include the bottom trawl, torch light, gill net, hand line, mid-water trawl, and longline.High biodiversity is found in this region, and the major economic species include the black pomfret, Parastromateus niger, skipjack, Katsuwonus pelamis, yellowfin tuna, Thunnus albacares, sergestid shrimp, Sergia lucens, crimson seabream, Dentex tumifrons, black seabream, Spondyliosoma cantharus, bigeye, Priacanthus macracanths, hairtail, Trichiurus lepturus, pomfrets, narrowbarred mackerel, Scomberomorus commerson, cuttle fishes, squids, etc.However, due to the increase in fishing pressure in recent decades, the coupling of the catch of economic species decreased, with smaller sized fish at catch provoking coastal fishery management by the local government [2].The Fisheries Agency, Taiwan, has implemented management measures which included prohibiting trawl and gillnet fishing within 3 miles from the coast and seasonal closure for trawl fishery that reduced the fishing effort.The yield of these areas decreased from 42,990 tons in 2000 to 28,476 tons in 2008.
Conventional single species' stock assessments, such as production models [3,4] and yield per recruit models [5], did not take account of the interactions among species.Hence, several multispecies' stock assessment methods taking account of the interactions among marine species based on single species models have been developed since late 1970s [6], i.e., multispecies production models (MSP, [7]), multispecies virtual population analysis [8], and multispecies bioenergetics models (MSBE, [9][10][11]).In addition, ecosystem-based assessment packages which can simulate the interactions among species in an ecosystem such as Ecopath with Ecosim [12] and NETWRK [13] have also been developed.Many more input parameters required for these packages resulted in higher uncertainty of the output [14].
Ecopath with Ecosim (EwE) is a package which is used to simulate the marine ecosystem mechanism and population dynamics based on the concept of predator-prey interactions [12,15].Regarding the energy equilibrium condition, Ecosim can be used to simulate the variation of biological parameters at different fishing pressures.In addition, Ecospace can be further used to conduct the simulations of spatial dynamics, which is useful for marine protected area planning.Ecopath with Ecosim has been used by several authors for ecosystem-based assessments and fishery management [14].For example, it has been used for simulations on the impact of removing large sharks in marine ecosystems [16], and simulating the impact of reducing the fishing effort of large sharks on other species in the pelagic ecosystem of the North Pacific [17].Resource management based on an ecosystem approach has been used in the Northeastern Atlantic [18].However, despite the work by Lin et al. [19,20] and Liu et al. [21] on small-scale coral reefs or estuaries in the coastal waters of Taiwan, this approach has never been used to assess the fishery impact on a large-scale marine ecosystem in Taiwan waters.Thus, how the Taiwanese fishery affects the marine ecosystem and its structure is still unknown.
In the waters of Guishan Island, a single species' stock assessment approach (virtual population analysis) has been applied on the bigeye, Priacanthus macrocanthus [2].However, due to the complex and multispecies' targeting nature of the coastal and offshore fisheries in this area, conventional single-species' stock assessment and management may not be appropriate for describing the impact of these fisheries [22].Thus, this study aims to assess the fishery impact on the marine ecosystem by using the Ecopath and Ecosim approach by constructing the marine ecosystem and examining the interactions among trophic levels of the resources of fisheries in the waters of Guishan Island.It is hoped that the results derived from this study can provide useful information for better managing the fishery resources in this region.

Data Source
Fishery data and biological information were obtained from the sampling vessels operating in the waters of Guishan Island (Figure 1).Catch of commercial species and trash fish were collected from sales records at Tahsi fish market, northeastern Taiwan.The species identification, species composition, condition factor, stomach content, maturation condition, and feeding cycle followed the methods mentioned by Wang et al. [23].

Ecopath Model Structure
Ecopath with Ecosim 6.20 (EwE 6.20) [24] was used to analyze the biological and biomass data of marine organisms in Guishan Island waters.EwE is a modelling package that is used to estimate the energy transfer in a freshwater or marine ecosystem and can simulate the fishery impact on marine ecosystems.The model is based on the massbalanced principle assumption, where the production of an animal species/group (P i ) and its immigration (I i ) in year i equals predatory loss (M2 i ), non-predatory mortality (M0 i ), harvest (Y i ), and the sum of emigration (E i ) [25].
Production + immigration = predatory losses + non-predatory mortality + harvest + emigration where B i is the biomass of species i in the study area, P i /B i is the production per unit biomass of species i, Ecopath assumes a mass balanced ecosystem, and P i /B i equals the total mortality of species i.
, where M0 i is other natural mortality of species i and EE i is utility rate of species i, which can be expressed as follows: , where EE i is the proportion utilized by other species in the ecosystem or fisheries for specie i.
, where M2 i is the mortality from predation for species i.
is the predation rate of predator j.I i and E i are the immigration and emigration rate of species i; the net immigration rate is the difference of these two values.I i and E i are assumed to be the same for all scenarios in simulations; thus, the net immigration rate was 0.

Estimation of Trophic Levels
Pauly et al. [26] mentioned that the decreases in trophic levels in global marine ecosystems were due to overfishing.To understand the variation of population structure of marine ecosystem of Guishan Island waters, the trophic level of each species/group was estimated by weighting the catch of each species/group [27]: where DC i,j represents the proportion in weight of prey for species i in the stomach of predator j, TL i represents the trophic level of predator i, the trophic level of primary producer was set as 1, and TL j represents the trophic level of prey j for predator i.
Omnivory index (OI) represents the variation of trophic level of the prey for predator i.When the trophic level of prey was the same OI = 0, OI increased with the diversity of trophic levels of prey.OI can be expressed as .

Mixed Trophic Impact (MTI)
Ulanowicz and Puccia's [29] matrix was used to estimate the mixed trophic impact (MTI) in Ecopath.The MTI represents the direct impact of species i on species j through predation minus the overall impact from other species on species j by predation.MTI was expressed as follows: MTI i,j = DC i,j − FC i,j where MTI i,j is the MTI of species i on species j, DC i,j is the direct impact from species i on species j by consumption, and FC i,j is the overall impact from other species on species j by consumption.

Keystone Species Identification
Keystone species is the species in an ecosystem that plays an important role in the food web of an ecosystem through predator-prey relation [30,31].The overall mixed trophic impact and keystone index derived from Ecopath were used to identify the keystone species in the ecosystem.
The overall mixed trophic impact (ε i ) can be expressed as follows: where m i,j is the relative impact of a slight increase in biomass of impacting group i on biomass of impacted group j.
The keystone index was estimated based on the following method [31]: , KS i is the keystone species index 1.
, where P i is the proportion of the biomass of species i to the summary of biomass of all species.

Simulations of Management Measures and Fishing Pressures Using Ecosim
After construing the ecosystem of Guishan Island by using Ecopath, Ecosim was used to simulate the impact on the ecosystem from various management measures.We used Ecosim to simulate the impact of management measures with the reduction of 10%, 20%, and 30% fishing effort for the subsequent 30 years.The core formula of Ecosim is as follows: , where g i is the net growth efficiency of species i, which equals the ratio of the production and predatory rate of species i.
F i is the fishing mortality of species i.

Model Construction
Based on the literature review, sampling vessel records, daily sales records of fish market, and Taiwan FishBase [32], a total of 28 species/functional groups were included in the model building.The ecosystem model's constructions are as follows: 2.7.1.Cetaceans Small-size dolphins were the major cetaceans in the study area, including the spinner dolphin (Stenella longirostris), common bottlenose dolphin (Tursiops truncates), striped dolphin (Stenella coeruleoalba), etc. [33].The biomass (B) of this group was set as 3.125 × 10 −3 ton/km 2 [33].Production per biomass (P/B) was set as 0.32 year −1 ; predation rate (Q/B) was set as 18.37 year −1 [34].The major prey of this group, based on their stomach analysis, were small fishes and cephalopods [35][36][37].

Small Skates and Rays
The small demersal skates and rays, such as the sepia stingray (Urolophus aurantiacus) and sharpspine skate (Okamejei acutispina), comprised the majority in this group.The P/B was set as 0.50, Q/B was 2.50, and EE was 0.67 for this group.The major prey of this group were demersal seashells, demersal benthivores, and shrimps.

Small Sharks
The sharks with a maximum size less than 100 cm TL are defined as small sharks, including the whitespotted bamboo shark (Chiloscyllium plagiosum) and Japanese spurdog (Squalus japonicas), which are demersal species.The P/B was set as 0.45, Q/B was 2.25, and EE was 0.58 for this group.Their prey are similar to that of skates and rays and mainly comprised demersal benthivores and shrimps.The skates and rays, and small sharks were mainly caught by bottom trawl or bottom longline fishery.

Dolphinfish
Two dolphinfish species were found in this marine ecosystem: common dolphinfish (Coryphaena hippurus) and Pompano dolphinfish (C.equiselis).The former comprised the majority in the study area.The P/B was set as 1.681, Q/B was 8.48, and EE was 0.85.Smallpelagic fishes, zooplanktons, shrimps, and scombrids were their prey.

Scombrids
The spotted chub mackerel (Scomber australasicus) and chub mackerel (S. japonicus) were the dominant species, and the former comprises 60~80% of the catch of the Taiwanese mackerel purse seiner fishery [38].The P/B was set as 3.37, Q/B was 32.57, and EE was 0.99 for this group.The major prey for this group were cephalopods, zooplanktons, and small fishes.

Shrimps
The majority of this group comprised sergestid shrimp (Sergia lucens) and other shrimps, which were caught by middle-water or bottom trawlers.The P/B was set as 10, Q/B was 40, and EE was 0.99.The major prey for this group was phytoplankton.

Lobsters and Crabs
This group included the spiny lobster (Panulirus penicillatus), Japanese spiny lobster (P.japonicas), and swimming crab (Portunus sanguinolentus), which were caught by the trap net.The P/B was set as 1.09, Q/B was 6.54, and EE was 0.70.

Oil Fish/Escolar
The two major species were oil fish (Ruvettus pretiosus) and escolar (Lepidocybium flavobrunneum), which were the bycatch of demersal longline fishery.These species inhabited at the depth of 200-400 m with high swimming ability and mainly fed on crustaceans, cephalopods, and small fishes.The P/B was set as 2.0, Q/B was 10.0, and EE was 0.90.

Carangids
This group included the bluefin trevally (Caranx melampygus), greater yellowtail (Seriola dumerili) and other species in Carangidae.The major prey items were demersal benthivores, shrimps, and zooplanktons for this group.The P/B was set as 1.79, Q/B was 7.14, EE was 0.95.

Pomfret
This group included the silver pomfret (Pampus argenteus), black pomfret, and other butterfish.The P/B was set as 0.57, Q/B was 2.85, and EE was 0.95.Their major prey item was zooplankton.

Sciaenids
This group included the large yellow croaker (Larimichthys crocea), small yellow croaker (L.polyactis), and other species of Sciaenidae.The P/B value was set as 1.5, Q/B was 7.5, and EE was 0.95.

Cephalopods
This group included the species Octopodidae, Sepiidae, and Loliginidae.The major prey of this group was phytoplankton, and this group was the major prey of several high-trophic level species.The P/B was set as 2.5, Q/B was set as 25, and EE was set as 0.99.

Sparids
This group includes the red seabream (Pagrus major), Crimson seabream (Dentex tumifrons), and longfinned bullseye (Cookeolus japonicas), etc.The species in this group have high commercial value and are carnivorous or omnivorous.Their major prey items were crustaceans and demersal benthivores.The P/B was set as 0.80, Q/B was 4.00, and EE was 0.95.

Clupeids
The major species of this group included the red-eye round herring (Etrumeus teres), spotted sardinella (Amblygaster sirm), and other pilchard species.These small coastal species fed on shrimps, sea algae, and larvae.The P/B was set as 2.00, Q/B was 10.00, and EE was 0.99.

Congers
This group included congers such as the daggertooth pike conger (Muraenesox cinereus).The major prey of this group was other fishes, crustaceans, and demersal benthivores.P/B was set as 1.00, Q/B was set as 5.00, and EE was set as 0.85.

Hairtail
Hairtail, including Trichiurus lepturus and T. japonicas, are the important commercial species in this area.They have diel movement behavior and mainly feed on small fishes and crustaceans.The P/B was set as 2.50, Q/B was 12, and EE was 0.95.

Blackbelly Triggerfish
This group included mudbank filefish (Paramonacanthus sulcatus) and other file fish that mainly fed on zooplankton.The P/B was set as 2.5, Q/B was 12, and EE was 0.95.

Anchovy
This group included Engraulis japonicas and Encrasicholina spp. in Engraulidae.These species mainly feed on crustaceans, zooplankton, and phytoplankton.The P/B was set as 3, Q/B was 15, and EE was 0.99.
The remaining species/group in this marine ecosystem included other pelagic fishes, demersal fishes, zooplanktons, phytoplankton, sea weeds, demersal benthivores, shellfish, and detritus.After referencing parameters used in similar marine ecosystems, the P/B ratio for phytoplankton was set as 400, while the P/B and Q/B ratios for zooplankton were set as 50 and 200, respectively [11,24].The information regarding the biomass, production (catch), and consumption (stomach contents) of 28 species/functional groups collected from the literature were used as the input parameters in Ecopath with Ecosim for analysis (Table 1).The pedigree analysis in Ecopath was used to estimate the uncertainty of input parameters.

Results and Discussion
The output of Ecopath and Ecosim for the ecosystem model in the waters of Guishan Island developed in this study was described as follows:

Trophic Position
In the Guishan Island marine ecosystem model developed in this study using Ecopath, trophic positions (TP) were determined based on the analysis of stomach contents and the literature review.Cetaceans occupied the highest trophic position with a value of 4.59, while demersal elasmobranchs, including small sharks (3.75) and skates and rays (3.73), held intermediate positions.Among teleost species, hairtail and congers exhibited the highest TP at 4.05, followed by dolphin fish at 4.0.The trophic positions of other teleost species ranged from 3.06 for anchovy to 3.98 for oil fish/escorla.Crustaceans displayed varying TP values, ranging from 2.06 for seashells to 3.53 for lobster (Table 1 and Figure 2).

Biomass Estimate
Based on the biological and stomach content information, the biomass of skates and rays, and small sharks estimated from Ecopath were 0.032 ton/km 2 and 0.038 ton/km 2 , respectively.The estimated biomass of pelagic species ranged from 0.001 ton/km 2 for oil fish

Biomass Estimate
Based on the biological and stomach content information, the biomass of skates and rays, and small sharks estimated from Ecopath were 0.032 ton/km 2 and 0.038 ton/km 2 , respectively.The estimated biomass of pelagic species ranged from 0.001 ton/km 2 for oil fish to 0.012 ton/km 2 for dolphinfish.While the estimated biomass of other teleost species ranged from 0.001 ton/km 2 for flat fish to 5.12 ton/km 2 for clupeids, the biomass of shrimp was 3.935 ton/km 2 and that of crabs and lobsters was 0.127 ton/km 2 , while the biomass of cephalopods was 0.114 ton/km 2 and that of demersal benthivores was 2.645 ton/km 2 .The biomass estimates of zooplankton, phytoplankton, and sea grass were 4.55 ton/km 2 4.740 ton/km 2 and 0.842 ton/km 2 , respectively (Table 1).

Fishing Mortality, Exploitation Rate, and Predation Mortality
Although all cetaceans have been included on the list of conservation animals in Taiwan since 1993 [39], incidental catches still occur through longline and drift net fisheries, resulting in a low fishing mortality (F) and exploitation rate (E) of 0.010 and 0.031 year −1 , respectively.Compared with fishes, cetaceans have life history characteristics of late maturity and small numbers of offspring.Although the F and E were low compared with other animals in the ecosystem, a continuous monitoring of the abundance is needed.
Fishermen mentioned that cetaceans will not only prey on fish with a high economic value in this area, but also cause a decline in fish populations.However, Ecosim results in this study showed that after the introduction of strict bans on fishing and bycatch of cetaceans, the biomass of cetaceans will increase significantly, but it will not have a great impact on other fish populations due to its small number of abundance compared with other animals.The variations of marine mammal resources in this ecosystem may not have a great impact on the balance of the ecosystem, but the economic impact on fishermen still requires further study.
The F and E were estimated to be 0.308 year −1 and 0.617 for skates and rays, and 0.261 year −1 and 0.580 for small sharks.Most species of this group had a low economic impact and were treated as trash fish in the past; however, some fishing vessels started to target these species due to the increase in their sales price in recent years.The F and E of dolphinfish, a high-commercial-value species, were estimated to be 0.411 year −1 and 0.245.The hairtail and flat fish had the highest F and E of 1.640 year −1 and 0.656, and 1.433 year −1 and 0.573, respectively, while crabs and lobsters had the lowest F and E of 0.172 year −1 and 0.158 (Table 2).The estimated inter-species/group predatory mortality in the Guishan Island marine ecosystem is showed in Table 3.The predatory mortality (PR) of dolphinfish mainly resulted from cetaceans, and PR was estimated to be 1.018 year −1 , that of scombrids was 1.989 year −1 , oil fish/escorla was 1.265 year −1 , and it was 0.262 year −1 for mackerels.The PR of other teleost species ranged from 0.372 year −1 for sparids to 0.779 year −1 for carangids.The highest PR occurred in anchovy at 2.084 year −1 , shrimps at 8.898 year −1 , and 1.711 year −1 for demersal benthivores (Table 2).

Trophic Impact and Keystone Species
Negative trophic impact on keystone species from fisheries was identified in the analysis (Figure 3).Pelagic fisheries such as longline and drift net fisheries had a strong negative trophic impact on cetaceans.On the other hand, the demersal fisheries such as bottom longline and bottom trawl fisheries had a strong negative trophic impact on demersal skates and rays, and small sharks.In addition, fishing also causes a great trophic impact on mackerels, hairtail, sparids, and carangids.However, positive trophic impacts resulting from the removal of predators were observed for dolphinfish and oil fish/escorla (Figure 3).Cetaceans exhibited a negative impact on prey species such as dolphinfish, mackerel, and oil fish/escorla but had a mild impact on other teleost species.The increase in shrimp biomass was noted to benefit predators like dolphinfish and mackerels.Phytoplankton and sea grass showed positive impacts due to a decrease in their predators, while zooplankton displayed similar effects (Figure 3).The increase in predators had a negative impact on shrimps but benefitted other species with higher trophic levels.The increase in the biomass of anchovy can provide more food for other teleost fishes and thus have a positive impact on these species and result in the increase of catch (Figure 3).
Phytoplankton, despite its low trophic position, emerged as the top keystone species in the Guishan Island marine ecosystem, being the major source of primary production with an abundant biomass.Zooplankton and anchovy were identified as the second keystone species, contributing significantly to the ecosystem with their ample biomass.Carangids and hairtail, despite their lower biomasses, exerted a higher influence due to their wide trophic breadth and substantial food intake (Table 4).

Primary Production Required for Fisheries
Harvesting high-trophic-level species such as cetaceans, dolphinfish, oil fish/escorla, and mackerels not only reduces the top-down control of the other species in the ecosystem but also necessitates more primary production to sustain fisheries.Simulations indicated that harvesting 1 kg cetacean requires 10,604 kg of primary production.Although the incidental catch of cetaceans was rare, their primary production requirement constituted

Primary Production Required for Fisheries
Harvesting high-trophic-level species such as cetaceans, dolphinfish, oil fish/escorla, and mackerels not only reduces the top-down control of the other species in the ecosystem but also necessitates more primary production to sustain fisheries.Simulations indicated that harvesting 1 kg cetacean requires kg of primary production.Although the incidental catch of cetaceans was rare, their primary production requirement constituted only 0.01% of the gross primary production.Dolphinfish, hairtail, and mackerels required 302, 317, and 272 kg primary production, respectively, comprising 0.05%, 0.18%, and 0.07% of the gross primary production.Anchovy and clupeids had the highest requirements at 2.00% and 2.51% of the gross primary production (Table 5).PPR: primary production, PPR/catch: primary production required per catch (kg), PPR/TotPP (%): primary production required per gross primary production (%).

Simulations with Reduction of Fishing Effort
Ecosim results with 10%, 20%, and 30% reduction of fishing effort for subsequent 30 years indicated similar results but the former two management measures did not lead to a significant increase in the resource levels of various species within this ecosystem, while simulations with a 30% reduction in fishing effort showed increased biomass for most species, particularly a 47% increase for hairtail and a 30% increase for small sharks, skates, and rays.Conversely, species at lower trophic levels, such as zooplankton, demersal benthivores, and shrimps, experienced decreased biomass due to the increase in their predator's biomass (Figures 4 and 5).In our 30-year simulations, the high trophic species, such as cetaceans, small skates and rays, small sharks, and dolphinfishs, were benefited from reduced fishing effort.The biomass of these functional groups increased greatly; however, the biomass of their main prey-scombrids, anchovy, and cephalopods-increased only to a very limited extent.The result suggested that a strong top-down control in Guishan Island marine ecosystem.Predators often play a crucial role in marine ecosystems due to predation.If large amount predators are removed, it can lead to changes in the biomass, and habitat of all related species in the ecosystem due to trophic cascade effect [15,16].

Model Uncertainty
Dame and Christian [14] suggested that the ecosystem-based analyses should consider the following uncertainties: (1) natural variations of input parameters, such as biomass of species/group or stomach content (consumption); (2) sampling methods, such as different selectivity of fishing gears may cause errors in biomass estimation; (3) model construction, such as that single species or aggregated functional group may cause its influence to be enhanced or dispersed [40]; (4) basic assumptions of ecological network analysis (ecological network analysis; ENA).For example, ecological network analysis research usually assumes that the ecosystem is in a steady state [41], which is not realistic.However, only a few studies of ENA have considered the uncertainty in model construction and analysis (Baird et al., 1998;Christian and Luczkovich, 1999) [42,43].The pedigree analysis indicated that the biomass had the highest uncertainty, and a lower uncertainty was found for the input parameters of higher-trophic-level species (Table 6).Although input parameter uncertainties were not fully addressed in this study, the use of Ecosim to simulate the impact of biomass and energy changes provides a dynamic perspective.This approach can make up for the insufficiencies of Ecopath, which assumes a steady-state ecosystem (Walters et al., 1997) [22].Pinnegar et al. [44] and Coll et al. (2006) [45] used principal component and cluster analysis to analyze the biological parameters of each species in a marine ecosystem and aggregated them into groups, and removed species that had a small impact on the entire ecosystem.This method can effectively reduce the uncertainty in the construction method of the ecological network model.However, in this study, due to the lack of species-specific biological information, mackerels, cephalopods, other surface fish, and other demersal fish were aggregated as functional groups without being based on multivariate analysis.In addition, the input parameters for some species were obtained from Fishbase or the literature, which might not represent the actual values in the study area.Future studies should focus on collecting more biological information on each species/group, particularly those of lower trophic levels, to improve the results of analyses.

Conclusions
This study aimed to construct the ecosystem structure of the waters of Guishan Island, incorporating 28 major species/functional groups using Ecopath with Ecosim.Fisheries were found to have a negative impact on major commercial species, with the exception of dolphinfish and oil fish, Lepidocybium spp., which benefitted from the reduction in their predators or competitors.Keystone species in this ecosystem were identified as phytoplankton, zooplankton, and anchovy, owing to their abundant biomass and significance as crucial food sources of many species in the ecosystem.It was observed that reducing the fishing effort could potentially increase the biomass of species occupying higher trophic positions in the ecosystem.
It is important to note that the analysis did not include trash fish caught by trawl fishery as bycatch, which comprises under-size individuals of many commercial species and non-commercial species.The growing proportion of trash fish in the area in recent years [23] underscores the need for a further investigation.Subsequent studies should focus on examining the composition of trash fish, estimating their abundance, and incorporating these data into simulations to enhance the accuracy of results.

J 20 Figure 2 .
Figure 2. The flow diagram of the food web derived from the marine ecosystem model of the Guishan Island, northeastern Taiwan.

Figure 2 .
Figure 2. The flow diagram of the food web derived from the marine ecosystem model of the Guishan Island, northeastern Taiwan.

Figure 3 .
Figure 3.The trophic and fishery impacts on species/functional group derived from the marine ecosystem model of the Guishan Island.Larger size of circles indicates larger impact.

Figure 3 .
Figure 3.The trophic and fishery impacts on species/functional group derived from the marine ecosystem model of the Guishan Island.Larger size of circles indicates larger impact.

Figure 4 .
Figure 4. 30-year trend of relative biomass of species/functional group in the marine ecosystem of the Guishan Island derived from simulations of 30% deduction in fishing effort.

Figure 4 . 20 Figure 5 .Figure 5 .
Figure 4.The 30-year trend of relative biomass of species/functional group in the marine ecosystem of the Guishan Island derived from simulations of 30% deduction in fishing effort.J. Mar.Sci.Eng.2023, 11, x FOR PEER REVIEW 16 of 20

Table 1 .
The input parameters and output of 28 species/functional groups from Ecopath for the Guishan Island marine ecosystem modeling.

Table 2 .
Estimated species/functional group specific total mortality, fishing mortality, predatory mortality, other mortality, and exploitation rate in the Guishan Island marine ecosystem.

Table 3 .
Estimated inter-species/functional group predatory mortality in the Guishan Island marine ecosystem.

Table 4 .
The keystone index and relative total impact of the species/functional group in the Guishan Island marine ecosystem.

Table 5 .
The species/functional group specific primary productivity and trophic position (TP) required to support the sustainability of Guishan Island marine ecosystem.

Table 6 .
Pedigree of biomass (B), production (P/B), predation (Q/B), consumption, and catch of Guishan Island marine ecosystem derived from Ecopath.The larger value indicates lower uncertainty.