# An introduction to the EcoTroph R Package: analyzing aquatic ecosystem trophic networks

Mathieu Colleter, Jérôme Guitton, Didier D. Gascuel

## ► To cite this version:

Mathieu Colleter, Jérôme Guitton, Didier D. Gascuel. An introduction to the EcoTroph R Package: analyzing aquatic ecosystem trophic networks. The R Journal, 2013, 5 (1), pp.98-107. doi:10.3180/ijte-01210175

HAL Id: hal-01210175

https://hal.science/hal-01210175v1

Submitted on 29 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# An Introduction to the EcoTroph R Package: Analyzing Aquatic Ecosystem Trophic Networks

by Mathieu Colléter, Jérôme Guitton and Didier Gascuel

**Abstract** Recent advances in aquatic ecosystem modelling have particularly focused on trophic network analysis through trophodynamic models. We present here a R package devoted to a recently developed model, EcoTroph. This model enables the analysis of aquatic ecological networks and the related impacts of fisheries. It was available through a plug-in in the well-known Ecopath with Ecosim software or through implementations in Excel sheets. The R package we developed simplifies the access to the EcoTroph model and offers a new interfacing between two widely used software, Ecopath and R.

## Introduction

In the face of the global overexploitation of marine resources and the fast degradation of ecosystems integrity, scientists developed new modelling approaches at the scale of the ecosystem. In fact, the main tool used for fisheries regulation is a stock approach which does not account for the trophic network linking marine ecological components. An important challenge is to analyze the aquatic ecological networks and the related impacts of fishery. There are several ways to define and represent an ecosystem. One widely used approach is based on the trophic level concept. Elton (1927) and Lindeman (1942) introduced this concept for describing aquatic ecosystems by assigning integer trophic levels (TLs) to the individual numbers, to the biomass or to the biological production by its component species. This approach differentiated between primary producers and detritus (TL = 1), first-order consumers (TL = 2), second-order consumers (TL = 3) ... The ecosystem is so represented as a pyramid of number, biomass or production, from low to high TLs. Odum and Heald (1975) developed this concept by implementing fractional trophic levels resulting from the diet of the individual and the trophic level of its preys. The emergence of Ecopath as a widely used approach and software for modelling aquatic ecosystems (Polovina, 1984; Christensen and Pauly, 1992) contributed in a major way to the prominence of TLs, especially as they were not an input, but an output of the model (i.e. estimated parameters). As the use of Ecopath spread worldwide with hundreds of application cases, so did the trophic level concept.

EcoTroph (ET) is an approach and software for modelling marine and freshwater ecosystems, entirely articulated around the TL concept (Gascuel, 2005; Gascuel and Pauly, 2009). It has been developed at the same time as the Ecopath worldwide expansion happened and incorporated into the Ecopath plug-in family (Gascuel et al., 2009). The first key idea of ET is that it deals with the continuous distribution of the biomass in an ecosystem as a function of continuous TL. The biomass enters the foodweb at TL = 1, generated by the photosynthetic activity of primary producers, or recycled from the detritus by the microbial loop. Between TL = 1 and TL = 2, the biomass is composed of mixotrophs only, and is usually low. If any, it is conventionally split between biomasses at TL = 1 and 2. Then, at TLs ≥ 2, the biomass is composed by heterotrophic organisms with mixed diet and fractional TLs resulting in a continuous distribution of biomass along TLs. The second key feature of ET is that the trophic functioning of aquatic ecosystems is modelled as a continuous flow of biomass surging up the foodweb, from lower to higher TLs, through predation and ontogenic processes. All the equations of the model are detailed in Gascuel et al. (2011). Such an approach, wherein species as such disappear, may be viewed as the final stage in the use of the TL metric for ecosystem modelling. It provides a simplified but useful representation of ecosystem functioning and impact of fishing. Thus, ET has been used both in theoretical contexts based on virtual ecosystems (Gascuel and Pauly, 2009; Gascuel et al., 2011), or in specific case studies to assess the current fishing impacts at the ecosystem scale (Gasche et al., 2012, in the South African Benguela ecosystem, Lassalle et al., 2012, in the Bay of Biscay, Tremblay-Boyer et al., 2011, for a worldwide analysis), or to analyze the effects of marine protected areas on the whole food web (Colléter et al., 2012 in Sénégal, Valls et al., 2012 in the Mediterranean Sea). Furthermore, ET enables the construction of a unique comparison framework for Ecopath models, the trophic spectrum. This display, based on ET key ideas, is a graphical representation of the ecosystem parameters, such as biomass, production, catch, fishing mortality, etc., along the trophic level (Gascuel et al., 2005). Examples of the use of the ET model and the associated trophic spectra analyses are provided in Figure 1.

This package is the first attempt to offer interfacing between R and Ecopath through its plug-in
<figure><img src="image_2.png" /></figure>

(a) Fishing impact on biomass in the Guinean ecosystem. Simulations refer to various multipliers (mE) of the current (i.e. 2004) fishing mortalities, from mE = 0 (i.e. virgin state) to mE= 5. Results highlight the strong impact of the current fishing effort (mE = 1, dashed lines) on the high TLs biomass: biomass is divided by 2 compared to the virgin state for all trophic levels higher than 3.8. In contrast, TLs around 3 benefit from a release of predation which approximately counterbalances the fishing impact.

<figure><img src="image_3.png" /></figure>

(c) Diagnosis on the effect of restricting fishery in the Port-Cros Marine Protected Area (France, Mediterranean Sea): trophic spectra of the relative biomass for a total closure of the fishery and 3 other hypothetical fishing scenarios. Results show that the MPA fulfills its conservation objective with current biomass very close to the virgin state for all trophic levels. Conversely, the other fishing scenarios lead to a significant decrease in the biomass of high trophic levels. It was concluded that the MPA mostly benefits predators and the functional trophic biodiversity.

<figure><img src="image_4.png" /></figure>

(b) Diagnosis on catch made in the Guinean ecosystem. Relative values of catch simulated for some trophic classes are expressed as a function of multipliers of the current fishing mortalities. Results show that high TLs (TL = 4 and 4.5) are fully exploited as an increase in the fishing pressure leads to a decrease in catch. In contrast, yields could be increased for the lower TLs. In other words, increasing fishing pressure would lead to a decrease in the mean trophic level of catch.

<figure><img src="image_5.png" /></figure>

(d) Analysis of the Marine Protected Area effect in the Bolong de Bamboung (Sénégal). Using the ET-Transpose routine, biomass trophic spectra are built from the 2003 Ecopath model (before the closure of the fishery) and the 2006-2008 model (enforced closure). The ET-Diagnosis routine is used to simulate the MPA effect, starting from the 2003 model and applying a null fishing mortality. Results show that the closure of the fishery can explain, at least partially, the increase observed in the biomass of high TLs (≥ 3.3). The decrease observed for intermediate TLs may not only results from a release in predation. Behavioral (refuge, flee) and environmental effects, that are not included in the model, may thus explain the observed differences.

**Figure 1:** Examples of the use of the EcoTroph model: diagnosis on the fishing impact on biomass (a) and catches (b) in the Guinean Ecosystem (from Gascuel et al, 2011); assessment of a Marine Protected Area (MPA) in Port-Cros, France (c) (from Valls et al, 2012) and in the Bolong de Bamboung, Sénégal (d) (from Colléter et al., 2012).
EcoTroph. The plug-in will use R and the EcoTroph (Guitton et al., 2013) package within the free Ecopath with Ecosim software. This enables the use of ET for a large panel through the plug-in for inexperienced R users, or directly the R package for the more expert ones. On the developer side, this way to link a rich user interface (developed in Microsoft VB) and a well known software in the fishery scientists group such as R is a way to set up a community. These researchers can so focus on the model improvements without taking into account the rich user interface which is time consuming. Along with several functions, we include an example dataset on the Guinean marine ecosystem (Gascuel et al., 2011) within the package. This article introduces, using the example dataset, the three main components of EcoTroph:

1. The data import and validation

2. The ET-Transpose tool

3. The ET-Diagnosis tool

A great deal of documentation, both introductory and advanced, is available on the ET website (http://sirs.agrocampus-ouest.fr/EcoTroph). The EcoTroph package requires the XML package (Lang, 2012) in order to load the model input parameters coming from the Ecopath software. The read.ecopath.model function will parse the data.

## Data import and validation

The EcoTroph package requires an input data table to run. The user has to load this dataset under the different possible formats ('.xls', '.csv', '.txt') with the functions read.table, read.csv... The ecopath_guinee dataset (Table 1) is an example of a suitable input table constructed for use with the EcoTroph package.

The variable names have to be specified and strictly the same as above: group_name (name of the group representing one or several species gathered together), TL (the trophic level of the group), biomass (the biomass of the group), prod (the production on biomass ratio or P/B) and accessibility (the ratio of the group biomass that would be caught assuming an infinite fishing pressure). These parameters generally come from Ecopath inputs or outputs, but can also be independent. The entry catch.1, catch.2, catch.whatyouwant is necessary if several fisheries do exist. The 0I column (the omnivory index, an Ecopath output parameter) is optional, it is used in the create.smooth function for an alternative smooth form ($\sigma_{\text{LN}} = \frac{\text{OI}}{\text{TL}_j}$, see below).

The check.table function was developed to check the compatibility of the input data table with the EcoTroph package:

check.table(ecopath_guinee)

In the example, no warning message appears as no error is made. If a message appears, the user has to correct the dataset in view of the comments. No missing values (NAs) are accepted as input, a yield column has to be entered with 0-values if no catches are registered. As well as the yield, the P/B of detritus groups (not entered in Ecopath) has to be set to 0. The check.table function converts the input dataset into a "data.frame" object if it is not already one.

A read.ecopath.model function was also implemented. It allows users to import data in '.xml' format exported from the Ecopath/EcoTroph plug-in or distributed by a web service (a database of Ecopath models has been set up to allow a meta-analysis at a world scale, http://sirs.agrocampus-ouest.fr/EcoTroph). This function formats the inputs so they meet the EcoTroph package requirements. Nevertheless, it is advisable to check the data consistency using the check.table function.

## ET-Transpose tool

### Creation of the Smooth function

The Smooth function returns a table allowing the conversion of data referring to specific taxons or functional groups (Ecopath trophic groups for example) into data referring to trophic classes. The major assumption of this function is that distributions of the trophic groups' biomass, yield, etc., around their mean trophic level follow a lognormal curve using the equation:

$$P_{ij} = \frac{1}{(TL_i - \text{shift})\sigma_j\sqrt{2\pi}} \exp \left[ -\frac{(\ln(TL_i - \text{shift}) - \ln(\overline{TL}_j - \text{shift}))^2}{2\sigma_j^2} \right]$$
<table><caption>**Table 1:** Data from the Ecopath model of the Guinean ecosystem (extracts).</caption><thead><tr><th></th><th>group_name</th><th>TL</th><th>biomass</th><th>prod</th><th>catch.1</th><th>catch.2</th><th>accessibility</th><th>OI</th></tr></thead><tbody><tr><td>1</td><td>Whales</td><td>4.01</td><td>0.0309</td><td>0.020</td><td>0.000</td><td>0.000</td><td>0.0</td><td>0.059</td></tr><tr><td>2</td><td>Dolphins</td><td>4.48</td><td>0.0433</td><td>0.070</td><td>0.000</td><td>0.000</td><td>0.0</td><td>0.331</td></tr><tr><td>3</td><td>Turtles</td><td>2.19</td><td>0.0296</td><td>0.150</td><td>0.000</td><td>0.000</td><td>0.0</td><td>0.338</td></tr><tr><td>4</td><td>Sea birds</td><td>3.81</td><td>0.0013</td><td>0.300</td><td>0.000</td><td>0.000</td><td>0.0</td><td>0.353</td></tr><tr><td>5</td><td>Rays+</td><td>3.97</td><td>0.3860</td><td>0.363</td><td>0.012</td><td>0.024</td><td>0.9</td><td>0.329</td></tr><tr><td>6</td><td>Sharks+</td><td>4.31</td><td>0.1050</td><td>0.410</td><td>0.007</td><td>0.003</td><td>0.8</td><td>0.633</td></tr><tr><td>7</td><td>Large pelagics</td><td>4.21</td><td>0.3840</td><td>0.850</td><td>0.025</td><td>0.069</td><td>0.8</td><td>0.263</td></tr><tr><td>8</td><td>Barracudas+</td><td>4.12</td><td>0.0583</td><td>0.920</td><td>0.009</td><td>0.022</td><td>0.9</td><td>0.259</td></tr><tr><td>9</td><td>Carangids</td><td>4.16</td><td>0.0627</td><td>1.000</td><td>0.010</td><td>0.024</td><td>0.8</td><td>0.139</td></tr><tr><td>10</td><td>Horse Makerels+</td><td>3.13</td><td>2.3330</td><td>0.700</td><td>0.000</td><td>0.115</td><td>0.8</td><td>0.366</td></tr><tr><td>34</td><td>Primary producers</td><td>1.00</td><td>69.0000</td><td>84.000</td><td>0.000</td><td>0.000</td><td>0.0</td><td>0.000</td></tr><tr><td>35</td><td>Detritus</td><td>1.00</td><td>290.0000</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.0</td><td>0.193</td></tr></tbody></table>

The lognormal distribution is defined by: a mean (the mean trophic level of the group, $TL_j$), a standard deviation ($\sigma_j$ denoted as sigmaLN in the R code) which is a measure of the trophic level variability within the group, and a shift parameter defining the theoretical trophic level characterised by a null variability in TL within group.

The create.smooth function enables the creation of this Smooth function using several input parameters. The parameter ecopath corresponds to the input data table (ecopath_guinee in the example). The parameter pas defining the splitting of trophic classes has by default a value of 0.1. The parameter smooth_type defines the form of the standard deviation (sigmaLN) wanted for the lognormal distribution. Three options are implemented:

1. If smooth_type=1 (choice by default), sigmaLN is constant. This constant sigmaLN is equal to the parameter sigmaLN_cst specified in the function, and has by default a value of 0.12. The shift parameter is set equal to 1.8 by default.

2. If smooth_type=2, this is equivalent to sigmaLN=smooth_param*ln(TL-0.05). The parameter smooth_param of this formula (also a parameter of the create.smooth function) defines the slope of the log-linear trophic level variability increase around the mean trophic level of the group. Based on our experience gained partially through observations, default parameters have been defined as follows: smooth_param=0.07 and shift=0.95. (No need to change the shift value in the function, let shift=NULL, it will be automatically set to 0.95 for smooth_type=2. Same thing for the parameter smooth_param.)

3. If smooth_type=3, sigmaLN for each group is equal to the omnivory index calculated by Ecopath divided by the mean trophic level of the group. A warning message will appear if OIs are equal to 0, they will be automatically changed to a value of 0.01. The parameter shift is by default equal to 0. (No need to change the shift value in the function, let shift=NULL, it will be automatically set to 0 for smooth_type=3.)

The create.smooth function returns a table of the distribution of each mean trophic level within trophic classes (i.e. how a given species or ecological group, characterized by a given and known mean trophic level, is distributed around this trophic level). This table will be used in the next step of the analysis for the construction of trophic spectra.

```r
# default choice, constant sigmaLN
create.smooth(ecopath_guinee)
# sigmaLN = smooth_param*ln(TL-0.05)
create.smooth(ecopath_guinee, smooth_type = 2)
```

A graphic function, plot(smooth), was developed in order to display this Smooth function. The input parameter is the table returned by the create.smooth function. It returns a plot with the lognormal curve for each present trophic class (see Figure 2).

```r
plotcreate.smooth(ecopath_guinee))
```
<figure><img src="image_6.png" /><figcaption>**Figure 2:** Output of the plot(smooth) graphic function, applied to the Guinean ecosystem example. Each curve represents the distribution across trophic levels for all groups whose mean trophic level is equal to the mean value of the distribution. Only curves related to existing groups are displayed.</figcaption></figure>

## Data transposition

The Transpose function enables the conversion of data referring to specific taxons or functional groups (Ecopath trophic groups for example) into data referring to trophic classes. This function uses the table returned by the create.smooth function. The concerned variables are the biomasses, or the catches or others... Using the Transpose function, these variables are distributed continuously along the trophic classes for each group. This function will be reused in the create.ETmain function to build a summary table with all the variables calculated by trophic class.

Transpose takes as input parameters the table returned by the create.smooth function (tab.smooth), the input data table (ecopath), and the name of the column the user wants to distribute by trophic class (column):

```
A <- create.smooth(ecopath_guinee)
# Transpose of the biomass column
T_biomass <- Transpose(A, ecopath_guinee, "biomass")
# Transpose of the catch.1 column
Transpose(A, ecopath_guinee, "catch.1")
```

Results can be displayed graphically using the plot(Transpose) function. It takes as input parameter the table returned by the Transpose function (tab_Trans). The user has the possibility to use a log scale for the y-axis (scale=log, the minimum value considered on the graph is conventionally set up at 1/10000 of the total biomass), and to enter a title (title):

```
plot tab_Trans, title = NULL, scale = NULL) returns the principal plots according to the selected column, in particular a plot by group and the associated trophic spectra (see Figure 3).
```

```
# title and log scale for the biomass
plot(T_biomass, title = "biomass", log)
```

## ET Main table creation

The create.ETmain function enables the creation of a summary table, ET Main, containing the principal variables by trophic class. This function can be used directly, there is no need of the previous steps Smooth and Transpose (however these steps are necessary to a good understanding of the EcoTroph model). It takes as input parameters the input data table (ecopath) and the parameters of the create.smooth function:

```
# constant sigmaLN
create.ETmain(ecopath_guinee)
# sigmaLN = smooth_param*ln(TL-0.05)
create.ETmain(ecopath_guinee, smooth_type = 2)
```
<figure><img src="image_7.png" /></figure>

**Figure 3:** Output of the plot (Transpose) graphical function applied to the biomass data per ecological group in the Guinean ecosystem (ecopath_guinee dataset). The top panel displays the distributions for each group, while the bottom panel displays the biomass trophic spectrum (BTS, i.e. the sum of all groups).
<table><thead><tr><th>TL</th><th>B</th><th>B_ACC</th><th>P</th><th>P_ACC</th><th>Kin</th><th>Kin_ACC</th><th>Y_tot</th><th>F_loss</th></tr></thead><tbody><tr><td>1</td><td>359.00</td><td>0.00</td><td>5796.00</td><td>0.00</td><td>16.14</td><td>1.00</td><td>0.00</td><td>0.00</td></tr><tr><td>2</td><td>25.73</td><td>0.00</td><td>1029.20</td><td>0.00</td><td>40.00</td><td>0.68</td><td>0.00</td><td>0.00</td></tr><tr><td>2.1</td><td>6.26</td><td>0.00</td><td>13.40</td><td>0.00</td><td>2.14</td><td>0.68</td><td>0.00</td><td>0.00</td></tr><tr><td>2.2</td><td>82.06</td><td>0.06</td><td>147.64</td><td>0.04</td><td>1.80</td><td>0.68</td><td>0.01</td><td>0.00</td></tr><tr><td>2.3</td><td>11.94</td><td>0.27</td><td>21.17</td><td>0.20</td><td>1.77</td><td>0.74</td><td>0.05</td><td>0.00</td></tr><tr><td>2.4</td><td>0.62</td><td>0.44</td><td>0.90</td><td>0.57</td><td>1.44</td><td>1.31</td><td>0.24</td><td>0.00</td></tr><tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr></tbody></table>

<table><caption>**Table 2:** create.ETmain(ecopath_guinee)$ET_Main results (extracts).</caption><thead><tr><th>TL</th><th>F_loss_ACC</th><th>N_loss</th><th>Fish_mort</th><th>Fish_mort_ACC</th><th>Selec</th><th>Time</th><th>N_loss_ACC</th></tr></thead><tbody><tr><td>1</td><td>0.00</td><td>1.73</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>NaN</td></tr><tr><td>2</td><td>0.22</td><td>43.41</td><td>0.00</td><td>0.15</td><td>0.00</td><td>0.06</td><td>-197.08</td></tr><tr><td>2.1</td><td>0.22</td><td>-24.00</td><td>0.00</td><td>0.15</td><td>0.00</td><td>0.06</td><td>-70.66</td></tr><tr><td>2.2</td><td>0.22</td><td>19.42</td><td>0.00</td><td>0.15</td><td>0.00</td><td>0.11</td><td>-16.80</td></tr><tr><td>2.3</td><td>0.25</td><td>31.57</td><td>0.00</td><td>0.18</td><td>0.02</td><td>0.17</td><td>-10.82</td></tr><tr><td>2.4</td><td>0.42</td><td>-1.89</td><td>0.38</td><td>0.55</td><td>0.70</td><td>0.22</td><td>-6.45</td></tr><tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr></tbody></table>

In the ET_Main output table (Table 2), some parameters directly come from the input table (biomass B, catch Y_tot), while some are calculated based on the definitions of EcoTroph parameters (e.g. Fishing mortality $Fish_mort = Y/B$, fishing loss rate $F_loss = Y/P$, see details in Gascuel et al., 2009).

```R
# constant sigmaLN
ET_Main <- create.ETmain(ecopath_guinee) \ $ET_Main
ET_Main
```

The create.ETmain function also returns intermediate tables, i.e. the following tables are also contained in the returned list object:

* `biomass`, the table returned by the Transpose function for the column biomass.

* `biomass_ACC`, the table biomass multiplied by the accessibility parameter.

* `prod`, the table corresponding to the production by trophic class.

* `prod_ACC`, the table prod multiplied by the accessibility parameter.

* `tab.smooth`, the table returned by the create.smooth function.

* Y, the table(s) returned by the Transpose function for the column catch. (x).

```R
> names創成(創成.ETmain(eco path_guinee))
[1] "ET_Main" "biomass" "biomass_ACC" "prod" "prod_ACC" "tab.smooth" "Y"
```

As previously, we developed a graphic function to display the main results: plot(ETmain). Different plots are created, e.g. the biomass trophic spectrum (BTS), the accessible biomass trophic spectrum (ABTS), the catch trophic spectrum (CTS). It takes as input parameter the list object returned by the create.ETmain function.

```R
plot(創成.ETmain(eco path_guinee), log)
```

Naturally all the returned graphics are not exhaustive. The user can construct other ones using the returned list object. Moreover the plot(ETmain) function is implemented with a log scale parameter for the different trophic spectra. However, this does not always provide a good representation. We really encourage users to test different scales for the y-axis. Some plots could be falsely interpreted with no awareness of the y-axis form importance.

## ET-Diagnosis simulation tool

ET-Diagnosis is used to simulate the effect of different fishery mortality scenarios on trophic spectra coming from ET-Transpose (see above). Thus different effort multipliers (Mul_eff variable), ranging conventionally from 0 to 5, are applied to the initial fishing mortalities $F_T$. Mul_eff=0 corresponds
<figure><img src="image_8.png" /><figcaption>**Figure 4:** Simulated biomass trophic spectra (BTS) of the Guinean ecosystem for various fishing effort multipliers ranging from 0 (no fishing) to 5. Such a graph highlights the fishing effects on trophic levels around 2.5 or higher than 3.4, with larger impact for the top predators (highest TLs). In contrast, the low and intermediate ones are less damaged.</figcaption></figure>

**Bibliography**

V. Christensen and D. Pauly. ECOPATH II—a software for balancing steady-state ecosystem models and calculating network characteristics. *Ecological Modelling*, 61(3–4):169–185, 1992. [p98]

M. Colléter, D. Gascuel, J. Ecoutin, and L. Tito de Morais. Modelling trophic flows in ecosystems to assess the efficiency of marine protected area (MPA), a case study on the coast of Sénégal. *Ecological Modelling*, 232:1–13, 2012. [p98]

C. Elton. *Animal Ecology*. The Macmillan Company, New York, 1927. [p98]

L. Gasche, D. Gascuel, L. Shannon, and Y. Shin. Global assessment of the fishing impacts on the Southern Benguela ecosystem using an EcoTroph modelling approach. *Journal of Marine Systems*, 90 (1):1–12, 2012. [p98]

D. Gascuel. The trophic-level based model: A theoretical approach of fishing effects on marine ecosystems. *Ecological Modelling*, 189(3–4):315–332, 2005. [p98]

D. Gascuel and D. Pauly. EcoTroph: Modelling marine ecosystem functioning and impact of fishing. *Ecological Modelling*, 220(21):2885–2898, 2009. [p98]

D. Gascuel, Y.-M. Bozec, E. Chassot, A. Colomb, and M. Laurans. The trophic spectrum: Theory and application as an ecosystem indicator. *ICES Journal of Marine Science*, 62(3):443–452, 2005. [p98]

D. Gascuel, L. Tremblay-Boyer, and D. Pauly. EcoTroph (ET): A trophic level based software for assessing the impacts of fishing on aquatic ecosystems. *Fisheries Centre Research Reports*, 17(1):1–82, 2009. [p98, 104]

D. Gascuel, S. Guénette, and D. Pauly. The trophic-level-based ecosystem modelling approach: Theoretical overview and practical uses. *ICES Journal of Marine Science*, 68(7):1403–1416, 2011. [p98, 100]

J. Guitton, M. Colléter, D. Gascuel, and P. Gatti. *EcoTroph: EcoTroph modelling support*, 2013. URL http://sirs.agrocampus-ouest.fr/EcoTroph/. R package version 1.5-1. [p100]

D. T. Lang. XML: *Tools for parsing and generating XML within R and S-Plus.*, 2012. URL http://CRAN.R-project.org/package=XML. R package version 3.9-4. [p100]
G. Lassalle, D. Gascuel, F. Le Loc'h, J. Lobry, G. J. Pierce, V. Ridoux, M. B. Santos, J. Spitz, and N. Niquil. An ecosystem approach for the assessment of fisheries impacts on marine top predators: The Bay of Biscay case study. *ICES Journal of Marine Science*, 69(6):925-938, 2012. [p98]

R. L. Lindeman. The trophic-dynamic aspect of ecology. *Ecology*, 23(4):399-417, 1942. [p98]

W. E. Odum and E. J. Heald. The detritus-based food web of an estuarine mangrove community. In L. E. Cronin, editor, *Estuarine Research*, volume 1, pages 265-286. Academic Press, New York, 1975. [p98]

J. J. Polovina. An overview of the ECOPATH model. *Fishbyte*, 2(2):5-7, 1984. [p98]

L. Tremblay-Boyer, D. Gascuel, R. Watson, V. Christensen, and D. Pauly. Modelling the effects of fishing on the biomass of the world's oceans from 1950 to 2006. *Marine Ecology Progress Series*, 442: 169-185, 2011. [p98]

A. Valls, D. Gascuel, S. Guénette, and P. Francour. Modeling trophic interactions to assess the effects of a marine protected area: Case study in the NW Mediterranean Sea. *Marine Ecology Progress Series*, 456:201-214, 2012. [p98]

Mathieu Colléter

Université Européenne de Bretagne, Agrocampus Ouest, UMR985 ESE, INRA, Agrocampus Ouest
Rennes

France
and

The Fisheries Center, University of British Columbia, Vancouver
Canada

mathieu.colleter@agrocampus-ouest.fr

Jérôme Guitton

Université Européenne de Bretagne, Agrocampus Ouest, UMR985 ESE, INRA, Agrocampus Ouest
Rennes

France

jerome.guitton@agrocampus-ouest.fr

Didier Gascuel

Université Européenne de Bretagne, Agrocampus Ouest, UMR985 ESE, INRA, Agrocampus Ouest
Rennes

France

didier.gascuel@agrocampus-ouest.fr