Skip to main content
Advertisement
Browse Subject Areas
?

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

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Modelling microplastic bioaccumulation and biomagnification potential in the Galápagos penguin ecosystem using Ecopath and Ecosim (EwE) with Ecotracer

  • Karly McMullen ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Visualization, Writing – original draft, Writing – review & editing

    k.mcmullen@oceans.ubc.ca (KM); j.alava@oceans.ubc.ca (JJA)

    Affiliation Ocean Pollution Research Unit, Institute for the Oceans and Fisheries, University of British Columbia, Vancouver, British Columbia, Canada

  • Félix Hernán Vargas,

    Roles Supervision, Writing – review & editing

    Affiliation The Peregrine Fund, Isla Santa Cruz, Galápagos Islands, Ecuador

  • Paola Calle,

    Roles Supervision, Writing – review & editing

    Affiliation Escuela Superior Politécnica del Litoral, ESPOL, Facultad de Ciencias de la Vida, ESPOL Polytechnic University, Guayaquil, Ecuador

  • Omar Alavarado-Cadena,

    Roles Writing – review & editing

    Affiliation Escuela Superior Politécnica del Litoral, ESPOL, Facultad de Ciencias de la Vida, ESPOL Polytechnic University, Guayaquil, Ecuador

  • Evgeny A. Pakhomov,

    Roles Conceptualization, Funding acquisition, Resources, Supervision, Writing – review & editing

    Affiliations Ocean Pollution Research Unit, Institute for the Oceans and Fisheries, University of British Columbia, Vancouver, British Columbia, Canada, Department of Earth, Ocean and Atmospheric Sciences, University of British Columbia, Vancouver, British Columbia, Canada

  • Juan José Alava

    Roles Conceptualization, Formal analysis, Methodology, Resources, Supervision, Writing – review & editing

    k.mcmullen@oceans.ubc.ca (KM); j.alava@oceans.ubc.ca (JJA)

    Affiliation Ocean Pollution Research Unit, Institute for the Oceans and Fisheries, University of British Columbia, Vancouver, British Columbia, Canada

Abstract

Bioaccumulation and biomagnification of anthropogenic particles are crucial factors in assessing microplastic impacts to marine ecosystems. Microplastic pollution poses a significant threat to iconic and often endangered species but examining their tissues and gut contents for contaminant analysis via lethal sampling is challenging due to ethical concerns and animal care restrictions. Incorporating empirical data from prey items and fecal matter into models can help trace microplastic movement through food webs. In this study, the Galápagos penguin food web served as an indicator species to assess microplastic bioaccumulation and biomagnification potential using trophodynamic Ecopath with Ecosim (EwE) modelling with Ecotracer. Empirical data collected from surface seawater near Galápagos penguin colonies, zooplankton, penguin prey, and penguin scat in October 2021 were used to inform the ecosystem model. Multiple scenarios, including a 99% elimination rate, were employed to assess model sensitivity. Model predictions revealed that microplastics can bioaccumulate in all predator-prey relationships, but biomagnification is highly dependent on the elimination rate. It establishes the need for more research into elimination rates of different plastics, which is a critical missing gap in current microplastic ecotoxicological and bioaccumulation science. Compared to empirical data, modelling efforts underpredicted microplastic concentrations in zooplankton and over-predicted concentrations in fish. Ultimately, the ecosystem modelling provides novel insights into potential microplastics’ bioaccumulation and biomagnification risks. These findings can support regional marine plastic pollution management efforts to conserve native and endemic species of the Galápagos Islands and the Galápagos Marine Reserve.

1. Introduction

In the Plasticene era, escalating plastic production has outpaced plastic elimination efforts, leading to its surge in the environment [1, 2]. Over 170 trillion plastic particles now float in the global ocean, with this estimate consistently rising [35]. These particles, many of which are considered microplastics (plastic particles 1000 nm to < 5 mm in size), become bioavailable to marine organisms, and have been widely documented in fish [6, 7] and invertebrates [8], including zooplankton [9], and, to a lesser extent, seabirds [1013], marine mammals [1417], and sea turtles [14, 1820]. Microplastics, due to their size and pervasive nature, enter biota through various pathways, including direct intake via water [2124] and air [21, 25], as well as indirectly through contaminated prey [16, 21, 2629]. The latter may represent a major pathway of microplastic ingestion for top predators [16]. These particles can cause physical damage to the gut and digestive tract, endocrine disruption due to sorption and leaching of toxic chemicals, and can impair feeding ability and predator avoidance [3034]. These sublethal but chronic health impacts may have cascading effects at the species or ecosystem level [35].

Bioaccumulation, bioconcentration, and biomagnification are common toxicokinetic and bioaccumulation science concepts often used in ecotoxicological and environmental risk assessments for water soluble (hydrophilic) chemicals such as pesticides, metals, pharmaceuticals, and personal care products [36], hydrocarbons, and ionic or neutral chemicals such as per- and polyfluoroalkyl substances (PFAS), as well as for hydrophobic or lipophilic substances, including persistent organic pollutants (POPs) [3739]. Bioaccumulation refers to the increase of a pollutant, i.e. microplastics, in an organism over time, or the gradual net uptake from all environmental compartments, including the surrounding environment and the food items [21, 28, 29, 39, 40]. Bioconcentration is a subcategory of bioaccumulation and refers to the gradual buildup of a contaminant in an organism from the water alone, excluding uptake from prey [29, 37]. Conversely, biomagnification implies an increase of the contaminant concentration at each trophic level of the food web. The contaminant concentration amplifies through the food web and organisms at higher trophic levels or apex predators exhibit elevated contaminant concentrations compared to organisms at lower trophic levels [39, 41].

These concepts have recently been applied to microplastic science [21, 26, 28, 29, 36, 42]. A recent literature review across multiple empirical studies shows evidence for bioaccumulation of microplastics within trophic levels [42]. Studies have also established trophic transfer of microplastics in natural-like laboratory settings [16, 27]. However, while species-specific bioaccumulation is likely to occur in marine species and it is a function of the elimination rate [21], evidence for biomagnification of microplastics between trophic levels is lacking [29, 42, 43]. The bioaccumulation and biomagnification capacity of microplastics are a cause for concern due to potential impacts on high trophic level species, and therefore require further research [40, 42].

Ecosystem and bioaccumulation modelling in tandem with the application of several ecotoxicological and bioaccumulation metrics can help to determine exposure levels and bioaccumulation potential, without needing to conduct invasive and lethal sampling or dissection of live individual organisms of threatened, endemic marine species. Within this rationale, food web bioaccumulation modelling is a feasible and useful tool to predict and estimate the bioaccumulation potential and ecotoxicological risks of microplastics to higher trophic level species [21].

The Galápagos Islands offer a unique opportunity to study microplastic bioaccumulation and biomagnification in relatively isolated and simplistic food webs. Recent studies have highlighted presence of microplastics in the archipelago [44, 45], but little is known about the trophic transfer and accumulation of microplastics in endemic and endangered marine species in the unique food webs of the islands. Galápagos seabirds, including the Galápagos penguin (Spheniscus mendiculus), may serve as sentinel species for biomonitoring microplastic pollution [46]. The Galápagos penguin is the only tropical, non-migrating penguin species with a decreasing population currently estimated at 1200 individuals [47]. This species faces several challenges from intense El Niño events and ocean warming [48, 49] to ocean plastics [44].

Using the Galápagos penguin as “the canary in the coal mine,” this study predicted the bioaccumulation and biomagnification of microplastics in their food web. The objectives of this study were to: first, understand the bioaccumulation behaviour of microplastics in the food web of the Galápagos penguin; second, predict the biomagnification potential of microplastics in the food web of the Galápagos penguin; and third, compare these findings to empirical data gathered in the Galápagos. Trophodynamic ecosystem modelling, including basic and advanced Ecopath and Ecosim (EwE) models along with the Ecotracer routine [50], were applied. Empirical seawater and biotic data gathered in the Galápagos were used as uptake and environmental input data and later were compared against the predicted data to assess model bias and corroborate the modelling performance. The research was designed to provide essential insights into the bioaccumulation potential and ecotoxicological risks of microplastics to higher trophic level species and inform regional policies to combat marine plastic pollution in the Galápagos Islands.

2. Methods

2.1 Ecosystem modelling theory

Ecopath with Ecosim (EwE) is a trophodynamic ecosystem modelling tool used to understand how different species in an ecosystem interact with one another, how they respond to changes in the environment, and how human activities might impact the ecosystem [50]. Ecopath is central to the software suite and provides a mass-balance snapshot of the ecosystem in question, while Ecosim offers a dynamic approach for temporal simulations based on predator-prey interactions using Lotka-Volterra and foraging arena theories [5052]. The EwE model together integrates biotic and abiotic components in an ecosystem by incorporating the principles of mass balance as well as a set of linear equations that describe and track the average flow of mass and energy between functional groups according to a diet composition matrix. Functional groups can be species or groups of species that have similar life-history or characteristics which are combined into biomass pools. The diet composition matrix is used to represent the flow of mass and energy within the ecosystem. It also accounts for energy losses over time through processes such as respiration, emigration, and decomposition [5052]. The core principles and mathematical equations of EwE are described in the user guide, which is accessible at http://www.ecopath.org [53].

2.2 Model descriptions

Two models were used to simulate microplastic movements through the Galápagos penguin food web, namely, a novel and simplistic food web model for the Galápagos penguin based on the species diet; and a trophic model of the Bolivar Channel Ecosystem (BCE), based on Ruiz and Wolff model [54] which includes an advanced snapshot of the BCE, with functional groups including seabirds, a proportion of which are Galápagos penguins (i.e., ~26% of the functional group ‘pool’). The latter incorporates a wider variety of groups and energy flows for a holistic ecosystem view, while the former isolates the Galápagos penguin food web to closely track microplastic movements within the specific food web. Both models were run with the Ecotracer routine to assess and compare bioaccumulation and biomagnification.

2.2.1 The Galápagos Penguin (GP) food web model.

The Galápagos Penguin (GP) food web model was constructed using EwE 6.6.8 [53]. The model was developed using a top-down approach by first analyzing the diet of the Galápagos penguin (Table 1). Based on limited available dietary data, the Galápagos penguin is suggested to feed primarily on small planktivorous fish, including sardines (Sardinops sagax), piquitingas (Lile stolifera), and salema (Xenocys jessiae), as well as fry or juveniles of mullets (Mugil spp.) [48, 5557; F. H. Vargas, pers. comm., 16 August 2023]. Closely related to the Galápagos penguin, the Humboldt Penguin (Spheniscus humboldti) also feeds on anchovies (Engraulis ringens), Araucanian herring (Strangomera bentincki), and silverside (Odontesthes regia), and to a lesser extent, cephalopods (Patagonian squid, Doryteuthis gahi and Humboldt squid, Dosidiscus gigas) and crustaceans (stomatopods and isopods) [58]. The diet composition of the Magellanic penguin (Spheniscus magellanicus), which is also closely related to the Galápagos penguin, was found to be similar [59]. Assessed penguins primarily ate anchovy (Engraulis anchoita) and thornfish (Bovichtus argentinus) [59]. This information provided a foundation for creating the GP model based on plausible diet matrix assumptions.

thumbnail
Table 1. Diet composition of Galápagos penguin (GP) EwE model.

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

The finalized diet matrix was compared and adjusted based on an existing EwE model available on Ecobase (http://ecobase.ecopath.org/) [60, 61], the Floreana island rocky reef ecosystem model, which includes the Galápagos penguin in the seabird functional group along with other marine species [62], as well as the BCE model [54]. Initially, cephalopods were part of the GP model; however, this species functional group was eventually excluded due to the unavailability of sufficient biomass data and the fact that they were not accounted for in the Floreana and BCE models. Conversely, barracuda (Barracuda pelicano) were added to the GP model because they were recognized as seabird prey in the Floreana and BCE models [54, 62]. It is reasonable to conclude that juvenile barracuda are likely to be preyed upon, given the preferred prey size for penguins [58, 63].

The final diet matrix consisted of ten functional groups, as shown in Table 1. Species selection and identification were done using scuba diving field guides [64, 65] as well as using the BCE model [54]. Diet information for penguin prey and lower trophic level species were accessed through FishBase (https://www.fishbase.se/) and SeaLifeBase (https://www.sealifebase.ca/) [6668]. Where data was not available, the diet information of the closest related species was used. Functional group biomasses (t/km2) were calculated based on existing models and reasonable estimations of biomass in g/m2. Production per biomass (P/B) was calculated based on the average lifespan of the species (i.e., the inverse of mortality). Galápagos penguins, for example, live up to 15 to 20 years [69], therefore a P/B of 0.07 is reasonable when assuming a mortality of y-1. Consumption per biomass (Q/B) was obtained from FishBase and SeaLifeBase [6668]. The Q/B parameter from FishBase was adjusted and recalculated for barracuda and mullets, given the penguins’ preference for smaller-sized prey; therefore, assuming the barracuda and mullet functional groups are comprised of juveniles, the size i.e., max total length (TL), was set to 30 cm for barracuda (S. idiates) and 25 cm for mullet (M. galapagensis), then FishBase calculations were rerun. Ecotrophic efficiency was left to be calculated by the EwE model. The sources of information and input data can be found in Table 2 & S1 Table in S1 File [48, 54, 57, 62, 7073].

thumbnail
Table 2. List of species and functional groups of the Galápagos penguin (GP) EwE model.

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

After inputting all parameter estimates, the model was found to be unbalanced, therefore parameters were manually adjusted until balance was achieved. The Q/B was reduced, and the P/B was increased for fish so that the P/Q ratio was near 0.2–0.3 based on EwE best practices [74, 75]. The biomass of herbivorous and predatory zooplankton was increased to account for their predation. Given the intentionally limited scope of the model and specific focus on the Galápagos penguin, diet import was reasonably assumed for most functional groups given that not all prey were included (e.g., not all barracuda prey).

2.2.2 Bolivar Channel Ecosystem (BCE) model.

The Bolivar Channel Ecosystem (BCE) model is a more advanced EwE model consisting of thirty functional groups known to inhabit the highly productive, upwelling zone between Isabela and Fernandina islands [54]. The original model was provided by courtesy of the authors to support the assessment of microplastic bioaccumulation in the area. Ruiz and Wolff [54] estimated biomass from observed sightings and censuses where applicable. The model includes one functional group entitled ‘seabirds’, comprising 74% of flightless cormorants (Phalacrocorax harrisi) and 26% Galápagos penguins (S. mendiculus). Input parameters underwent a series of resampling through the Ecoranger resampling routine to select a random set of input values from normal distributions of the input parameters. The model is fully described in Ruiz and Wolff [54].

2.3 Modelling bioaccumulation of microplastics with EwE Ecotracer routine

The Ecotracer routine is an EwE module used to track and assess the bioaccumulation potential of pollutants in marine food webs over time [50, 76, 77]. Ecotracer uses Ecosim temporal simulations to predict the flow of contaminants through the food web [7678]. The contaminant concentration over time in each functional group is based on the flow rates from Ecosim, as well as decay, elimination, and physical exchange rates. The linear dynamical equation for changes in contaminant concentration over time for a given functional group (pool) or species and the dynamic changes in contaminant concentration can be found in Ecopath with Ecosim: A User’s Guide Ecosystem effects of invertebrate fisheries View project Ecosystem Fmsy estimations for data rich stocks [52].

The contaminant concentrations in immigrating and emigrating organisms were considered to be negligible for the purposes of this modelling work to simplify the model and because the Galápagos penguins are endemic and residents of the Galápagos Islands and do not undertake migration.

Input parameters were gathered from observed data from the field [79], laboratory studies, and existing modelling work outlined in Table 3. Microplastic concentrations (particle size > 10 μm) in the environment were obtained [79] and were entered as input data for the initial and inflow environmental concentrations of microplastics in the Ecotracer module for both EwE models. The direct uptake rates were derived from the observed zooplankton ingestion rates of anthropogenic particles [79]. This information was utilized as the exclusive source for the microplastic uptake, based on the premise that zooplankton occupy a pivotal position at the bottom of aquatic food webs and therefore represent the initial point of entry for the bioaccumulation of microplastics in the ecosystem [21]. Considering the purpose of the modelling work, which aimed to evaluate the levels and magnitude of bioaccumulation and biomagnification of microplastics in the Galápagos penguin food web, alternative sources of microplastic exposure such as inhalation of airborne particles or direct uptake from the environment were assumed to be negligible and thus not considered.

thumbnail
Table 3. Ecotracer environmental parameters, input data, and respective sources used in the Galápagos Penguin (GP) EwE model.

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

Microplastic elimination or egestion rates were adopted from the food web bioaccumulation model developed by Alava [21] or collected from the best available data (See Table 4). Briefly, the egestion or elimination rate of microplastic is computed from the retention time (τ) as the elimination or egestion rate is inversely related to the retention or residency time: τ = 1/kE; thus, solving for kE:kE = 1/τ [21]. Finally, an average decay rate was set at ~0.0283 per year (2.825% per year), based on plastic particles decay estimates (ranging 0.65% to 5% per year), accounting for weight loss due to solar radiation and oxygenation, documented by Everaert et al. [80]. Ecotracer was run for a simulation period of 100 years.

thumbnail
Table 4. Data values of microplastic retention times to calculate the egestion (elimination) rates with reference sources used for the Galápagos Penguin (GP) EwE model.

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

2.4 Model scenarios

To explore the microplastic bioaccumulation potential and biomagnification capacity in the food web, the model was run using four different scenarios: (1) A baseline scenario using observed environmental anthropogenic particle concentrations [79] and egestion rates from literature; (2) a high environmental concentration scenario, which assumed a higher abundance of microplastics in the environment, based on the upper limit of the standard deviation of the aforementioned observed microplastic abundance data; (3) a low environmental concentration scenario, which assumed a lower abundance of microplastics in the environment, based on the lower limit of the standard deviation of the microplastic abundance data; and finally, (4) a 99% egestion rate scenario which assumes all microplastics are excreted in under 24 hours, using the environmental concentration data from baseline scenario (S2 Table in S1 File). See S1-S4 Tables in S1 File for additional Ecotracer input data.

2.5 Bioaccumulation and biomagnification metrics

The application of bioaccumulation and biomagnification criteria and metrics for microplastics were based on Alava [21, 89, 93].

Bioaccumulation factor (BAF). The bioaccumulation factor (BAF) was calculated by comparing the amount of microplastics accumulated by aquatic biota to the total concentration of microplastics present in both the aquatic environment and diet (e.g., the predator’s prey). A BAF greater than 1 indicates plausible bioaccumulation under steady state (e.g., BAF >1).

(1)

Where Ci is the microplastic concentration in the predator (in units of g/kg), Co is the concentration in the environment and Cj is the concentration in the diet (prey) in units of g/kg.

Bioconcentration factor (BCF). The bioconcentration factor (BCF) was calculated by comparing the amount of microplastics accumulated by the aquatic organism or biota (in units of g/kg) to the total concentration of microplastics present in the water or aquatic environment (in units of g/L). A BCF greater than 1 indicates bioconcentration in the organism (e.g., BCF >1).

(2)

Where Ci is the microplastic concentration in biota (g/kg) and Cw is the concentration in the water (g/L).

Predator-prey biomagnification factor (BMFTL). The predator-prey biomagnification factor (BMFTL) was calculated by comparing the amount of microplastics accumulated by the predator in units of g/kg to the total concentration of microplastics (g/kg) in the prey, divided by the difference in trophic levels. A BMF greater than 1 indicates plausible trophic biomagnification (e.g., BMFTL >1).

(3)

Where Ci is the microplastic concentration in the predator, Cj is the concentration in prey, TL is the trophic level for the predator (i) and prey (j).

Trophic magnification factor (TMF). The trophic magnification factor (TMF) is a commonly employed metric in food web analysis that measures the biomagnification of pollutants at various trophic levels [21, 9092]. The TMF is calculated as the antilog of the regression slope of the linear regression between the logarithmic-transformed concentrations of microplastics (Log MPs) predicted in the GI tract of organisms of the food web and their trophic level, TL [21]. The linear equation is as follows.

(4)

Where b is the slope, TL is the trophic level, and a is the y-intercept. The TMF can be expressed as the equivalent exponential mathematical terms expressed as TMF = 10b, where b is the slope.

(5)

The TMF (slope, b) is statistically evaluated using a significance level (α) of 0.05. A TMF > 1 (b > 0) indicates that the contaminant biomagnifies in the food web. A TMF < 1 (b < 0) indicates trophic dilution of the contaminants, while a TMF = 1 (b = 0) indicates no change in contaminant concentrations among organisms of a food web [21, 90].

2.6 Sensitivity assessment and model bias

The sensitivity of the model was assessed by testing changes in the environmental concentrations and the functional group elimination rates. This was conducted by comparing the outcomes of the model through four of the different scenarios (S2 Table in S1 File), including high and low environmental concentrations as well as high (e.g., > 24 hours) and low (e.g., < 24 hours) retention times, equivalent to slow (high retention) or fast (low retention) elimination rates, respectively. The most sensitive parameters were determined by assessing the variance in output data according to changes in environmental concentrations and elimination rate.

2.7 Model bias

A model bias (MB) approach was applied to assess the performance of the food web model and corroborate the projections of microplastics under the three abiotic concentrations’ scenarios (scenario 1–3) and one conservative egestion rate scenario (scenario 4). The performance of the model was analyzed in terms of the model bias ratio: (6) where CBP,iMP and CBO,iMP are the model calculated (predicted) and observed microplastic concentrations in species i, respectively. This analysis was done by comparing the projected microplastics concentrations in biota (zooplankton, anchovy, and mullet) yielded by the model simulations to the empirical data measured in free-ranging zooplankton, wild-caught anchovy, and mullets in waters of the Galápagos National Park [79]. The microplastic and anthropogenic particle mean concentrations for zooplankton, anchovy, and mullet were 7.38x10-7, 7.69x10-10, and 1.48x10-8 g/kg, respectively [79].

3. Results

3.1 Galápagos penguin food web model in Ecopath

A successfully balanced GP model was constructed with a total of ten functional groups of varying trophic levels (Fig 1; Table 5). Analysis of the energy flow diagram revealed the highest energy flows from microalgae/phytoplankton to herbivorous zooplankton, as well as from anchovy, herring, sardines, and salema to barracuda and to the Galápagos penguin. A high biomass of macroalgae and detritus was necessary to maintain the balance of the model, and matched output of the BCE model [54]. Detritus, microalgae/phytoplankton, and macroalgae made up trophic level 1, followed by herbivorous zooplankton, mullet, predatory zooplankton, and anchovy, herring, sardines, salema as trophic level 2, in order of lowest to highest. Predatory zooplankton had a higher trophic level than mullets, which is expected given mullet’s preference for detritus compared to predatory zooplankton consuming other zooplankton. Likewise, anchovy, herring, sardines, and salema had a similar trophic level compared to predatory zooplankton, which is not surprising given similar food preferences, namely other zooplankton. The Galápagos penguin and barracuda were the high trophic level species included in the model. Ecotrophic efficiencies (EE) were comparable to Ruiz & Wolff [54] and ranged from 0 to 0.68 where the highest EE resulted in predatory zooplankton, herbivorous zooplankton, and microalgae/phytoplankton.

thumbnail
Fig 1. Galápagos penguin EwE food web model energy flow diagram.

Circles represent the amount of biomass and labels indicate the respective functional group. Thicker lines between functional groups reflect higher energy flows from one functional group to another, and trophic levels are indicated on the y-axis.

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

thumbnail
Table 5. Parameter output from the Galápagos penguin EwE food web model.

https://doi.org/10.1371/journal.pone.0296788.t005

3.2 Microplastic bioaccumulation and biomagnification via Ecotracer

3.2.1 Galápagos penguin (GP) food web model.

Utilizing the baseline Ecotracer scenario for the GP food web model, output data revealed a rapid increase in contaminant concentration (g) per biomass (kg) until around year 5, after which it shifted to a more gradual increase and eventually plateaued, reaching steady state around year 60 (Fig 2A). The microplastic concentration in the Galápagos Penguin plateaued at 7.6x10-5 g/kg, while the microplastic concentration in their primary prey item, planktivorous fish, reached a steady state at 4.93x10-8 g/kg. The Galápagos penguin was found to have the highest level of microplastics per biomass, followed by barracuda, anchovy, sardine, herring, and salema, and then predatory zooplankton (Fig 2A). In contrast, species or functional groups at lower trophic levels had lower concentrations of microplastics per unit of biomass.

thumbnail
Fig 2. Galápagos penguin model microplastic concentration over time and trophic level.

(A) Bioaccumulation simulations using the EwE Ecotracer routine, showing the projections of microplastics (MPs) bioaccumulation in the Galápagos penguins food web model under the baseline scenario with seawater initial concentrations = 0.00411 t/km2. The simulations for the bioaccumulation include the elimination rates based on the literature. For zooplankton, as the key trophic level for the initial uptake of microplastics, the ingestion rate of microplastics and anthropogenic particles by zooplankton (i.e., 7.38x10-7 g/kg/day) was used, based on [79]. (B) Linear regressions showing the significant relationship between predicted concentrations of microplastics (log-transformed data) and trophic levels in the GP model at year 100. The antilog of the regression slope was used to determine TMF.

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

The baseline simulation results predicted the plausible bioaccumulation and biomagnification of microplastics in the functional groups of species within this model. The relationship between the logarithmic-transformed concentration of microplastics and the species’ trophic levels shows that the TMF was significantly greater than 1 (i.e., TMF > 1; when the slope is statistically different or greater than zero [(b > 0)] in a positive, significant linear regression). The TMF increased from 7.08 at year 1 to 61.7 at year 100 (Fig 2; S5 Table in S1 File). By year 100, there was a substantial difference between the Galápagos penguin and Barracuda (trophic level 3) compared to lower trophic level groups, indicating biomagnification in the model. The regression between the logarithmic microplastic concentration (g/kg) revealed a slope of 1.79 with a TMF = 61.7 (p < 0.05), further indicating a statistically significant relationship for biomagnification in the model (Fig 2B).

The Galápagos penguin had the highest BAF, regardless of what scenario was assessed, while the mullet had the lowest BAF and BCF (Fig 2; S6 Table in S1 File). The BMF increased with each trophic level and was high (BMFTL = 1239 and BMFTL = 1796) for the Galápagos penguin, with planktivorous fish and mullet being considered as prey, respectively, though it was substantial for all prey items from the Galápagos penguin. The BMFTL was the lowest for planktivorous fish (BMFTL = 3.06) as well as for predatory zooplankton with herbivorous zooplankton as prey (BMFTL = 3.05).

3.2.2 Bolivar Channel (BCE) ecosystem model.

The BCE model with the Baseline Ecotracer scenario yielded similar results to the GP food web model in terms of overall trend, but with lower microplastic concentrations per functional group. Output data once again revealed a rapid increase in contaminant concentration (g) per biomass (kg) until around year 5, and diverging slightly from the GP model, concentrations plateaued earlier, reaching steady state around year 35 as shown in Fig 3A. Microplastic concentrations in the seabirds plateaued at around 5.9x10-7 g/kg, two orders of magnitude lower than results from the GP model. Microplastic concentrations in the penguins’ primary prey items, planktivorous fish, reached steady states at 6.63x10-8 g/kg (Fig 3A), which is within the same order of magnitude as yielded in the GP model. The functional groups with the highest concentrations of microplastics per biomass displayed a pattern consistent with the one established in the GP model. The functional groups of interest with elevated microplastic concentrations included seabirds as the highest, followed by barracuda, predatory zooplankton, and finally, small planktivorous reef fish (Fig 3A). Unlike the GP model, predatory zooplankton had slightly higher microplastic concentrations compared to small planktivorous fish and herbivorous zooplankton had higher concentrations of microplastics than both detritivorous species, mullets, and lobsters.

thumbnail
Fig 3. BCE model microplastic concentration over time and trophic level.

(A) Bioaccumulation simulations using the EwE Ecotracer routine, showing the projections of microplastics (MPs) bioaccumulation in the BCE model under the baseline scenario where seawater initial concentrations = 0.00411 t/km2. The simulations for the bioaccumulation include the elimination rates based on the literature. For zooplankton, as the key trophic level for the initial uptake of microplastics, the ingestion rate of microplastics and anthropogenic particles by zooplankton (i.e., 7.38x10-7 g/kg/day) was used, based on [79]. (B) Linear regressions showing the significant relationship between predicted concentrations of microplastics (log-transformed data) and trophic levels in the BCE model at year 100. The antilog of the regression slope was used to determine TMF.

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

Similar to the GP model simulations, BCE simulation resulted in plausible bioaccumulation and biomagnification of microplastics in the functional groups of species in the model. The relationship between the concentration of microplastics and the species’ trophic levels shows that the TMF was once again significantly greater than 1 (slope > 0; Fig 3A). The TMF ranged from 6.17 at year 1 to 16.6 at year 100 (S7 Table in S1 File). At year 100, the TMF was lower than estimated by the GP model. Seabirds and barracuda displayed a similar trend, yielding higher bioaccumulation and biomagnification potentials than other functional groups of interest. Of particular interest were the projections for predatory marine mammals and sharks, having the highest biomagnification capacity of all functional groups in the model. Overall, the regression indicated a high biomagnification potential as the logarithmic microplastic concentration (g/kg) versus the trophic level revealed a statistically significant slope of 1.22, i.e., TMF = 16.7 (p < 0.001; Fig 3; S7 Table in S1 File).

As in the GP model, seabirds had the highest BAF and BCF potential, regardless of what prey scenario was assessed, and the mullet once again had the lowest BAF and BCF (S8 Table in S1 File). The resulting BMFTL was different than the GP model. Foremost, the planktivorous reef fish revealed BMFTL < 1, indicating lack of biomagnification. The Galápagos penguin revealed positive BMFTL across all prey items, where the highest was once again with prey item lobster. That said, BMFTL was substantially lower than the GP model; the highest BMFTL in the BCE model was a factor of 76 compared to a factor of 2308 in the GP model.

3.3 Model sensitivity

To assess the sensitivity of the model parameters, four different scenarios were evaluated for each model, resulting in a total of eight simulations to examine microplastic bioaccumulation and biomagnification. The simulation outcomes showed that the elimination rate was the most sensitive parameter. This was evidenced by the fact that the results from this scenario displayed the greatest variability in terms of BAF and BCF, as well as smaller ranges in BMFTL across both the GP model and the BCE model, as shown in Table 6.

thumbnail
Table 6. Bioaccumulation, bioconcentration, and biomagnification factors in the Galápagos penguin and BCE models across scenarios.

https://doi.org/10.1371/journal.pone.0296788.t006

In the 99% elimination rate scenario for both the GP and BCE models, the BAF and BCFs were substantially lower, with values that were two to four orders of magnitude lower than those estimated in the high and low abundance scenarios as well as the baseline scenarios (Table 6). The BAF and BCF for the high and low abundance simulations were comparable to the baseline scenario for the GP and BCE models.

For the GP model using the 99% elimination rate scenario, the BMFTL ranged from 0.02 to 49.6 depending on the predator and prey (Table 6). Conversely, the BMFTL had a much broader range (e.g., 3.06 to 2308) across the other three scenarios. When comparing the BMFTL for the Galápagos penguins and planktivorous fish, the baseline scenario resulted in the highest BMFTL (BMFTL = 1239). This value was comparable to the high and low abundance scenarios but was substantially higher than the BMFTL of 27.7 as estimated by the 99% elimination scenario for the same species combination.

For the BCE model, the BMFTL ranged from -0.04 to 76.71 in the 99% elimination rate scenario, which was more comparable to the other three scenarios yielding -2.73 to around 76 BMFTL, depending on the predator-prey combination (Table 6). Once again, the BMFTL for the Galápagos penguins and planktivorous fish was much higher in the high and low abundance scenarios as well as the baseline scenarios (BMFTL = ~8.5) compared to the 99% elimination rate scenario (BMFTL = 0.10).

Simulation outcomes of the different scenarios yielded interesting results for TMF calculations, where all scenarios, except for the 99% elimination rate, predicted the plausible trophic biomagnification of microplastics in the biomass species’ functional groups (Table 7; Fig 4A and 4B). For the 99% elimination rate, some slopes were negative (i.e., year 1 in the GP model and years 1 to 100 in the BCE model), indicating lack of trophic magnification, however regressions were not significant (p > 0.05) in the GP modelled scenario (Fig 4A) or the BCE (Fig 4B). A negative slope, if significant, would indicate trophic dilution of microplastics (i.e., a decline in microplastic concentrations as the trophic level increases, TMF < 1 [slope < 0]).

thumbnail
Fig 4. Microplastic concentrations and trophic level linear regressions.

Linear regressions showing the relationship between predicted concentrations of microplastics (log-transformed data) and trophic levels for: A) the Galápagos penguin (GP) model; and, B) the Bolivar Channel Ecosystem (BCE) model, for year 1, 25, 50, and 100, under four different scenarios including low microplastic abundance in seawater, high microplastic abundance in seawater, 99% elimination rates for all functional groups, and the baseline scenario which includes elimination rates based on available literature. The antilog of the regression slope was used to determine TMF. This panel was made to visualize the overall trends, not to display the fine detail of each graph.

https://doi.org/10.1371/journal.pone.0296788.g004

thumbnail
Table 7. Apparent trophic magnification factors (TMF) and regression statistics.

https://doi.org/10.1371/journal.pone.0296788.t007

The highest TMF value was predicted at year 100 in the baseline and high abundance scenarios running in the GP model, i.e., TMF = 61.7 (Table 7). The highest TMF predicted in the BCE model was 16.7 in the high abundance and baseline scenario. Across both models under the 99% elimination rate scenario, there was very low TMF values from year one to year 100. TMF ranged from 0.48 to 2.57 in the GP model, while the lowest TMF values, ranging from 0.19 to 0.37 were predicted from year one to 100 in the BCE model (Table 7). TMF rose substantially from year 1 to year 25 and then gradual increases in TMF through to year 100 were predicted for the other scenarios tested (Fig 4A and 4B).

3.4 Model bias

The projected concentration of MPs in the environment, and biota including zooplankton, anchovies, and mullet were compared to empirically collected data from similar abiotic (environmental) and biotic compartments. The outcomes of the MB ratio analysis revealed systematic under prediction (MB < 1) of microplastic concentrations in zooplankton, across all scenarios, with a MB ranging from 0.002 to 0.1 at low abundance scenario and high abundance scenario, respectively while using the GP model (Fig 5A and S1A, S1B Fig in S1 File). Similarly, under prediction for zooplankton occurred in the BCE model, with MB ranging from 0.002 to 0.08 in the low and baseline scenarios, respectively (Fig 5B and S1C, S1D Fig in S1 File). Conversely, the baseline, low, and high abundance (1, 2, 3) scenarios revealed over-prediction of microplastic concentrations (MB > 1) in fish ranging from 5.4 (low abundance) to 105.1 (high abundance) for anchovies in the GP model and 9.2 (low abundance) to 72.8 (baseline) for anchovies in the BCE model, as well as, 13.0 (low abundance) to 252.3 (high abundance) for mullets in the GP model and 11.8 (high abundance) to 65.5 (baseline) for mullets in the BCE model (Fig 5B; S1 Fig in S1 File). Results for fish are different in the 99% elimination rate scenario, namely, fish had the closest concentration values to the observed data (MB = ~1) (S1B and S1D Fig in S1 File). In the GP model, the 99% elimination rate scenario yielded MB of 0.29 and 0.65 and the BCE yielded MB of 0.38 and 0.48 for the anchovies and mullets, respectively.

thumbnail
Fig 5. Model bias.

Assessment of the model bias (MB = CBP,iMP/CBO, iMP; where CBP,iMP and CBO,iMP are the model predicted and observed microplastic concentrations, respectively, in zooplankton, anchovies, and mullets) and performance of the (A) Galápagos penguin (GP) and (B) Bolivar Channel Ecosystem (BCE) models by comparing the simulations of microplastic concentrations averaged from year 1–100. Four different scenarios are compared for model bias, including low microplastic abundance in seawater, high microplastic abundance in seawater, 99% elimination rates for all functional groups, and the baseline scenario which includes elimination rates based on available literature. The dotted line indicates MB = 1 (i.e., predicted microplastics data is equivalent to the observed microplastic concentration).

https://doi.org/10.1371/journal.pone.0296788.g005

4. Discussion

Ecosystem modelling simulation results in this study revealed bioaccumulation potential in all predator-prey combinations in both GP and BCE models across all scenarios. Biomagnification of microplastics was apparent in all simulations, except for the 99% elimination rate scenario, which showed a lack of trophic magnification with insignificant negative slopes (p > 0.05) in both GP and BCE models. Biomagnification was highly dependent on microplastic egestion rates. As well, the microplastic concentration in zooplankton was systematically under-predicted in both GP and BCE models, whereas microplastic concentrations in fish were over-predicted, with the closest concentration values to the observed data being seen in the 99% elimination rate scenario and the baseline scenario for fish and zooplankton respectively (Fig 5 & S1 Fig in S1 File).

These results are comparable to limited existing microplastic bioaccumulation and biomagnification simulations using EwE, which is not surprising given similar methods and parameters used [93, 94]. Ma and You’s [94] EwE simulations found microplastics bioaccumulate quickly in fish food webs of Baiyangdian Lake, China. Likewise, Alava et al. [93] ran Ecotracer through an extensive list of 20 marine ecosystem models available through EcoBase. The projections of microplastic concentrations in biota revealed that top predators are likely exposed to higher levels of microplastics accumulated through their diet (e.g., prey items). However, the elimination or egestion rate was likewise a key factor in determining the net bioaccumulation behaviour of microplastics. Emphasis was placed on better understanding of the role of retention times and elimination rates of microplastics in different functional groups. Likewise, microplastic bioaccumulation modelling in a cetacean food web [21] was comparable to EwE results, indicating that species-specific bioaccumulation of microplastics is likely, while biomagnification is highly dependent on species-specific elimination rates.

According to a comprehensive meta-analysis [42] evidence has been found for bioaccumulation of microplastics in marine species; however, biomagnification of these pollutants in the food web has yet to be confirmed by field data. A select number of field studies have observed trophic transfer of microplastics [16, 27, 95], but there is a disagreement as to whether microplastics are eliminated or retained in the GI tract or gut [95]. Microplastic dilution, for example, has been found in mussels and fish in the Persian Gulf [26]. To better understand this phenomenon, Miller et al. [29] explored microplastic movement through a coral reef food web and found bioconcentration evident in zooplankton, crustaceans, and fish, but no bioaccumulation or biomagnification.

Bioaccumulation and biomagnification currently represent an important debate in the field of microplastic science and uncertainty remains. While conducting research to investigate these ecotoxicological factors, it is important to consider (1) microplastics are a class of contaminants and may move differently within a food web depending on their physical and chemical characteristics [96], (2) microplastic input into the ocean is increasing year over year, and (3) evidence suggest ingestion of microplastics provides no advantage to marine organisms in a changing ocean. Thus, though microplastic biomagnification is not agreed upon, it is still important to reduce the anthropogenic emissions of microplastics to the sea.

First, it is critical to consider microplastics as a class of contaminants [96]. It is reasonable that their behaviour in the gastrointestinal tract is varied based on their characteristics (e.g., size, shape, microbial biosphere, additives, dyes), as different microplastics and associated chemicals will react differently in the gut [97, 98]. There is evidence to suggest that small microplastics and nanoplastics (<1 μm), not measured in this study, can pass through tissue membranes, translocate to tissues, and enter the blood stream [99102]. Chemical additives or harmful bacteria on microplastics can adsorb within the biota [103, 104] and may or may not bioaccumulate, biomagnify, or multiply within the organism [22, 105]. Specific shapes and sizes may be more toxic [98] and prone to lower or higher GI tract retention times, influencing microplastic bioaccumulation [21].

Secondly, plastic and microplastic pollution inputs into the ocean is increasing in the Plasticene [32, 80, 106109]. The increase in these anthropogenic particles will inevitably lead to microplastics accumulation in ecosystems, whether it will be a marine ecosystem accessible to researchers and the public, or deep-sea ecosystems that are not commonly sampled.

Lastly, numerous studies note prominent ingestion of microplastics when organisms are exposed to plastic particles in the environment or in laboratory conditions. Microplastic ingestion has been linked with several primarily sublethal health effects [31, 34, 110, 111] and, therefore, it is plausible that ingested microplastics do not offer advantages in stressful changing oceanic conditions further harmed by overfished resources and climate change effects [112]. Instead, the presence of microplastics or nanoplastics may independently pose a significant concern and may be magnified by other environmental challenges [113].

Overall, future research should prioritize laboratory assessments of microplastic accumulation across trophic levels, with particular attention on egestion rates and GI tract retention times, physical-chemical characteristics of retained microplastics, and ecotoxicological health effects. Future modelling work should explore different interactions with microplastic at the primary consumer level. For example, microplastics can be added as a functional group in EwE and mediation applied to increase plastic consumption when phytoplankton abundance is low and vice versa [74, 93]. It would also be prudent to create model scenarios with varying rates of increase for microplastics in the specific ecosystem. Future modelling should likewise consider multiple-anthropogenic stressors, including the impact of climate change forcing and El Niño Southern Oscillation (ENSO) events in the Galápagos penguin food web when there are drastic changes on sea surface temperature (negative anomalies), disruption of primary production and pronounced food shortages when fewer fish serving as prey for penguins during El Niño are available [48, 114]. These multifactorial stressors may well affect the exposure pathways, bioaccumulation, and elimination of microplastics in marine top predators, including seabirds like the Galápagos penguins.

5. Conclusion

It is imperative to continue prioritizing efforts to reduce the input of microplastics into vulnerable ecosystems and food webs, such as that of the endangered Galápagos penguin. Despite ongoing research, the biomagnification of microplastics remains unclear, and additional studies are needed to fully understand this phenomenon. This study identified a knowledge gap in microplastic elimination rates, which are needed to determine biomagnification potential.

As microplastic research continues, it is key to continue efforts raising awareness and mitigating microplastic pollution entering waterways, while not losing sight of other pressing threats to the world’s oceans. It is critical to adopt a balanced, holistic approach when considering oceanic threats, including plastics, but also overfishing, and climate change, in order to effectively protect and preserve iconic species like the Galápagos penguin and our precious ocean environment.

Ultimately, this trophodynamic ecosystem modelling work provided predictions and insights on the potential bioaccumulation and biomagnification risks of microplastics as a global pollutant of concern to support regional marine plastic pollution management efforts for the conservation of native and endemic species of the Galápagos Islands and the Galápagos Marine Reserve.

Please refer to the Supplementary Information for additional figures and tables, including model input values, scenario details, functional group details, linear regression data, bioaccumulation factors (BAF), bioconcentration factors (BCF), predator-prey biomagnification factors (BMFTL), and model bias figures. All data are accessible online at SEANOE (https://doi.org/10.17882/97201) [115].

Supporting information

S1 File. Supplementary information includes model input values, scenario details, functional group details, linear regression data, bioaccumulation factors (BAF), bioconcentration factors (BCF), predator-prey biomagnification factors (BMFTL), and model bias figures.

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

(DOCX)

References

  1. 1. Borrelle SB, Ringma J, Law KL, Monnahan CC, Lebreton L, Mcgivern A, et al. Predicted growth in plastic waste exceeds efforts to mitigate plastic pollution. Science. 2020 Sep 18;369(6510): 1515–8. pmid:32943526
  2. 2. Cordier M, Uehara T. How much innovation is needed to protect the ocean from plastic contamination? Sci, Total Environ. 2019 Jun 20;670:789–99. pmid:30921712
  3. 3. Eriksen M, Cowger W, Erdle LM, Coffin S, Villarrubia-Gómez P, Moore CJ, et al. A growing plastic smog, now estimated to be over 170 trillion plastic particles afloat in the world’s oceans—Urgent solutions required. PLOS One. 2023 Mar 8;18(3):e0281596. pmid:36888681
  4. 4. Van Sebille E, Wilcox C, Lebreton L, Maximenko N, Hardesty BD, Van Franeker JA, et al. A global inventory of small floating plastic debris. Environmental Research Letters. 2015 Dec 8;10(12).
  5. 5. Eriksen M, Lebreton LCM, Carson HS, Thiel M, Moore CJ, Borerro JC, et al. Plastic Pollution in the World’s Oceans: More than 5 Trillion Plastic Pieces Weighing over 250,000 Tons Afloat at Sea. PLOS One. 2014 Dec 10;9(12).
  6. 6. Azevedo-Santos VM, Gonçalves GRL, Manoel PS, Andrade MC, Lima FP, Pelicice FM. Plastic ingestion by fish: A global assessment. Environmental Pollution (Barking, Essex: 1987). 2019 Dec;255(Pt 1):112994. pmid:31541837
  7. 7. Jovanović B, Gökdağ K, Güven O, Emre Y, Whitley EM, Kideys AE. Virgin microplastics are not causing imminent harm to fish after dietary exposure. Mar. Pollut. Bull. 2018 May 1;130:123–31. pmid:29866538
  8. 8. Ugwu K, Herrera A, Gómez M. Microplastics in marine biota: A review. Mar. Pollut. Bull. 2021 Aug 1;169:112540. pmid:34087664
  9. 9. Botterell ZLR, Beaumont N, Dorrington T, Steinke M, Thompson RC, Lindeque PK. Bioavailability and effects of microplastics on marine zooplankton: A review. Environmental Pollution. 2019 Feb 1;245:98–110.
  10. 10. Ivar Do Sul JA, Costa MF. The present and future of microplastic pollution in the marine environment. Environmental pollution. 2014 Feb 1;185:352–64. pmid:24275078
  11. 11. Kühn S, van Franeker JA. Quantitative overview of marine debris ingested by marine megafauna. Mar. Pollut. Bull. 2020 Feb 1;151:110858. pmid:32056640
  12. 12. Roman L, Hardesty BD, Hindell MA, Wilcox C. A quantitative analysis linking seabird mortality and marine debris ingestion. Sci Rep. 2019 Dec 1;9(1). pmid:30824751
  13. 13. Wilcox C, Van Sebille E, Hardesty BD, Estes JA. Threat of plastic pollution to seabirds is global, pervasive, and increasing. Proc. Natl. Acad. Sci. USA. 2015 Sep 22;112(38):11899–904. pmid:26324886
  14. 14. López-Martínez S, Morales-Caselles C, Kadar J, Rivas ML. Overview of global status of plastic presence in marine vertebrates. Global Change Biology. 2021 Feb;27(4):728–37. pmid:33111371
  15. 15. Lusher AL, Hernandez-Milian G. Microplastic Extraction from Marine Vertebrate Digestive Tracts, Regurgitates and Scats: A Protocol for Researchers from All Experience Levels. Bio. Protoc. 2018 Nov;8(22). pmid:35036474
  16. 16. Nelms SE, Galloway TS, Godley BJ, Jarvis DS, Lindeque PK. Investigating microplastic trophic transfer in marine top predators. Environmental Pollution. 2018 Jul 1;238:999–1007. pmid:29477242
  17. 17. Zantis LJ, Carroll EL, Nelms SE, Bosker T. Marine mammals and microplastics: A systematic review and call for standardisation. Environmental Pollution. 2021 Jan 15;269:116142. pmid:33288297
  18. 18. Lynch JM. Quantities of Marine Debris Ingested by Sea Turtles: Global Meta-Analysis Highlights Need for Standardized Data Reporting Methods and Reveals Relative Risk. Environ. Sci. Technol. 2018 Nov 6;52(21):12026–38. pmid:30251846
  19. 19. Yaghmour F, Al Bousi M, Al Naqbi H, Samara F, Ghalayini T. Junk food: A preliminary analysis of ingested marine debris by hawksbill Eretmochelys imbricata and olive ridley Lepidochelys olivacea sea turtles (Testudines: Cheloniidae) from the eastern coast of the United Arab Emirates. Mar. Pollut. Bull. 2021 Dec 1;173. pmid:34741925
  20. 20. Schuyler Q, Hardesty BD, Wilcox C, Townsend K. Global Analysis of Anthropogenic Debris Ingestion by Sea Turtles. Conservation Biology. 2014 Feb;28(1):129–39. pmid:23914794
  21. 21. Alava JJ. Modeling the Bioaccumulation and Biomagnification Potential of Microplastics in a Cetacean Foodweb of the Northeastern Pacific: A Prospective Tool to Assess the Risk Exposure to Plastic Particles. Front. Mar. Sci. 2020 Sep 22;7.
  22. 22. Koelmans AA. Modeling the role of microplastics in bioaccumulation of organic chemicals to marine aquatic organisms. A critical review. In: Bergmann M, Gutow L, Klages M, editors. Marine Anthropogenic Litter. Springer International Publishing; 2015. p. 309–324.
  23. 23. Roch S, Friedrich C, Brinker A. Uptake routes of microplastics in fishes: practical and theoretical approaches to test existing theories. Sci. Rep. 2020 Dec 1;10(1). pmid:32127589
  24. 24. Setälä O, Fleming-Lehtinen V, Lehtiniemi M. Ingestion and transfer of microplastics in the planktonic food web. Environmental Pollution. 2014;185:77–83. pmid:24220023
  25. 25. Amato-Lourenço LF, Carvalho-Oliveira R, Júnior GR, dos Santos Galvão L, Ando RA, Mauad T. Presence of airborne microplastics in human lung tissue. J. Hazard. Mater. 2021 Aug 15;416:126124. pmid:34492918
  26. 26. Akhbarizadeh R, Moore F, Keshavarzi B. Investigating microplastics bioaccumulation and biomagnification in seafood from the Persian Gulf: a threat to human health? Food Additives & Contaminants: Part A. 2019 Nov 2;36(11):1696–708. pmid:31393815
  27. 27. Farrell P, Nelson K. Trophic level transfer of microplastic: Mytilus edulis (L.) to Carcinus maenas (L.). Environmental Pollution. 2013 Jun 1;177:1–3. pmid:23434827
  28. 28. Goswami P, Vinithkumar NV, Dharani G. First evidence of microplastics bioaccumulation by marine organisms in the Port Blair Bay, Andaman Islands. Mar. Pollut. Bull. 2020 Jun 1;155. pmid:32469778
  29. 29. Miller ME, Motti CA, Hamann M, Kroon FJ. Assessment of microplastic bioconcentration, bioaccumulation and biomagnification in a simple coral reef food web. Science of the Total Environment. 2023 Feb 1;858:159615. pmid:36309288
  30. 30. Bergmann M, Gutow L, Klages M. Marine Anthropogenic Litter. Springer Nature; 2015.
  31. 31. Cole M, Lindeque P, Fileman E, Halsband C, Goodhead R, Moger J, et al. Microplastic ingestion by zooplankton. Environ. Sci. Technol. 2013 Jun 18;47(12):6646–55. pmid:23692270
  32. 32. Everaert G, De Rijcke M, Lonneville B, Janssen CR, Backhaus T, Mees J, et al. Risks of floating microplastic in the global ocean. Environmental Pollution. 2020 Dec 1;267:115499. pmid:33254632
  33. 33. Kühn S, Bravo Rebolledo EL, Van Franeker JA. Deleterious effects of litter on marine life. In: Bergmann M, Gutow L, Klages M, editors. Marine Anthropogenic Litter. Springer International Publishing; 2015. p. 75–116.
  34. 34. Wright SL, Thompson RC, Galloway TS. The physical impacts of microplastics on marine organisms: a review. Environmental pollution. 2013 Jul 1;178:483–92. pmid:23545014
  35. 35. Galloway TS, Cole M, Lewis C. Interactions of microplastic debris throughout the marine ecosystem. Nature ecology & evolution. 2017 Apr 20;1(5):0116. pmid:28812686
  36. 36. Chormare R, Kumar MA. Environmental health and risk assessment metrics with special mention to biotransfer, bioaccumulation and biomagnification of environmental pollutants. Chemosphere. 2022 Sep 1;302:134836. pmid:35525441
  37. 37. Arnot JA, Gobas FAPC. A review of bioconcentration factor (BCF) and bioaccumulation factor (BAF) assessments for organic chemicals in aquatic organisms. Environmental Reviews. 2006 Dec 1;14(4):257–97.
  38. 38. Gobas FAPC, De Wolf W, Burkhard LP, Verbruggen E, Plotzke K. Revisiting bioaccumulation criteria for POPs and PBT assessments. Integrated Environmental Assessment and Management: An International Journal. 2009 Oct;5(4):624–37. pmid:19552497
  39. 39. Gobas FAPC Morrison HA. Bioconcentration and bioaccumulation in the aquatic environment. In: Boethling R, Mackay D, editors. Handbook of property estimation methods for chemicals: Environmental health sciences. CRC press; 2000 Mar 29. p. 189–231.
  40. 40. Provencher JF, Ammendolia J, Rochman CM, Mallory ML. Assessing plastic debris in aquatic food webs: What we know and don’t know about uptake and trophic transfer. Environmental Reviews. 2019;27(3):304–17.
  41. 41. Kelly BC, Ikonomou MG, Blair JD, Morin AE, Gobas FAPC. Food web-specific biomagnification of persistent organic pollutants. Science. 2007 Jul 13;317(5835):236–9.
  42. 42. Miller ME, Hamann M, Kroon FJ. Bioaccumulation and biomagnification of microplastics in marine organisms: A review and meta-analysis of current data. PLOS One. 2020 Oct 16;15(10):e0240792. pmid:33064755
  43. 43. Covernton GA, Cox KD, Fleming WL, Buirs BM, Davies HL, Juanes F, et al. Large size (>100-μm) microplastics are not biomagnifying in coastal marine food webs of British Columbia, Canada. Ecological Applications. 2022 Oct;32(7):e2654.
  44. 44. Jones JS, Porter A, Muñoz-Pérez JP, Alarcón-Ruales D, Galloway TS, Godley BJ, et al. Plastic contamination of a Galapagos Island (Ecuador) and the relative risks to native marine species. Science of the Total Environment. 2021 Oct 1;789:147704. pmid:34049146
  45. 45. Alfaro-Núñez A, Astorga D, Cáceres-Farías L, Bastidas L, Soto Villegas C, Macay K, et al. Microplastic pollution in seawater and marine organisms across the Tropical Eastern Pacific and Galápagos. Sci. Rep. 2021 Dec 1;11(1).
  46. 46. Boersma PD. Penguins as marine sentinels. Bioscience. 2008 Jul 1;58(7):597–607.
  47. 47. International BirdLife. Spheniscus mendiculus. 2020 [cited 14 November 2023]. In: The IUCN Red List of Threatened Species 2020 [Internet]. Available from: https://dx.doi.org/10.2305/IUCN.UK.2020-3.RLTS.T22697825A182729677.en.
  48. 48. Vargas H, Harrison S, Rea S, Macdonald DW. Biological effects of El Niño on the Galápagos penguin. Biol. Conserv. 2006 Jan;127(1):107–14.
  49. 49. Wolff M, Ruiz DJ, Taylor M. El Niño induced changes to the Bolivar Channel ecosystem (Galapagos): Comparing model simulations with historical biomass time series. Mar. Ecol. Prog Ser. 2012 Feb 23;448:7–22.
  50. 50. Christensen V, Walters CJ. Ecopath with Ecosim: Methods, capabilities and limitations. Ecol. Modell. 2004 Mar 1;172(2–4):109–39.
  51. 51. Christensen V, Pauly D. ECOPATH II—a software for balancing steady-state ecosystem models and calculating network characteristics. Ecological modelling. 1992 Jun 1;61(3–4):169–85.
  52. 52. Christensen V, Walters CJ, Pauly D. Ecopath with Ecosim: A User’s Guide. 2005 [cited April 2023]. In: Fisheries Centre, University of British Columbia, Vancouver, Canada and ICLARM, Penang, Malaysia [Internet]. Available from: https://www.researchgate.net/publication/267193103
  53. 53. Christensen V, Walters C, Pauly D, Forrest R. Lenfest Ocean Futures Project. Ecopath with Ecosim Version 6: User Guide. 2008 [cited April 2023]. In: Lenfest Ocean Futures Project 2008 [Internet]. Available from: https://repository.oceanbestpractices.org
  54. 54. Ruiz DJ, Wolff M. The Bolivar Channel Ecosystem of the Galapagos Marine Reserve: Energy flow structure and role of keystone groups. J. Sea. Res. 2011 Aug;66(2):123–34.
  55. 55. Mills KL. Diving behaviour of two Galápagos Penguins Spheniscus mendiculus. Marine Ornithology. 2000 Apr 15;28:75–9.
  56. 56. Vargas H, Lougheed C, Snell H. Population size and trends of the Galápagos Penguin Spheniscus mendiculus. Ibis. 2005 Apr;147(2):367–74.
  57. 57. Wilson RP, Wilson MPT. Foraging Ecology of Breeding Spheniscus Penguins. Penguin biology. 1990;181–206.
  58. 58. Herling C, Culik BM, Hennicke JC. Diet of the Humboldt penguin (Spheniscus humboldti) in northern and southern Chile. Mar. Biol. 2005 May;147(1):13–25.
  59. 59. Fernandez SJ, Yorio P, Ciancio JE. Diet composition of expanding breeding populations of the Magellanic Penguin. Marine Biology Research. 2019 Jan 2;15(1):84–96.
  60. 60. Colleter M, Valls A, Guitton J, Gascuel D, Pauly D, Christensen V. Global overview of the applications of the Ecopath with Ecosim modeling approach using the EcoBase models repository. Ecol. Modell. 2015;302:42–53.
  61. 61. Colléter M, Valls A, Guitton J, Lyne M, Arreguín-Sánchez F, Christensen V, et al. EcoBase: a repository solution to gather and communicate information from EwE models. 2013 [cited April 2023] [Internet]. Available from: https://institut-agro-rennes-angers.hal.science/hal-01103442
  62. 62. Okey TA, Banks S, Born AF, Bustamante RH, Calvopiña M, Edgar GJ, et al. A trophic model of a Galápagos subtidal rocky reef for evaluating fisheries and conservation strategies. Ecol. Modell. 2004 Mar 1;172(2–4):383–401.
  63. 63. Fernández Severini MD, Villagran DM, Buzzi NS, Sartor GC. Microplastics in oysters (Crassostrea gigas) and water at the Bahía Blanca Estuary (Southwestern Atlantic): An emerging issue of global concern. Reg. Stud. Mar. Sci. 2019 Nov 1;32.
  64. 64. Constant P. Marine life of the Galápagos: the diver’s guide to fishes, whales, dolphins and marine invertebrates. Odyssey; 2007.
  65. 65. Humann P, DeLoach N. Reef fish identification: Galápagos. New World Publications Incorporated; 2003.
  66. 66. Froese R, Pauly D. FishBase 2000: Concepts, design and data sources. ICLARM, Los Baños, Laguna, Philippines; 2000.
  67. 67. Froese R, Pauly D. FishBase. World Wide Web electronic publication; 2022.
  68. 68. Palomares MLD, Pauly D. SeaLifeBase. World Wide Web electronic publication; 2022.
  69. 69. Rafferty JP. Galapagos penguin. 2020 [cited April 2023]. In: Encyclopedia Britannica [Internet]. Available from: https://www.britannica.com/animal/Galapagos-penguin. 2020.
  70. 70. Karpouzi VS. Modelling and mapping trophic overlap between fisheries and the world’s seabirds. M.Sc. Thesis. The University of British Columbia. 2005. Available from: https://open.library.ubc.ca/soa/cIRcle/collections/ubctheses/831/items/1.0074880
  71. 71. Grove JS, Lavenberg RJ. The fishes of the Galápagos Islands. Stanford: Stanford University Press; 1997. p. 863.
  72. 72. Muck P, Rojas de Mendiola B, Antonietti E. Comparative studies on feeding in larval anchoveta (Engraulis ringens) and sardine (Sardinops sagax). In: Pauly D, Mendo J, Tsukayama I, editors. The Peruvian upwelling ecosystem: dynamics and interactions. ICLARM Conference Proceedings 18, 438 p. IMARPE, Callao, Peru; GTZ, Eschborn, Germany; and ICLARM, Manila, Philippines; 1989. p. 86–96.
  73. 73. Serra RT, Tsukayama I. Sinopsis de datos biologicos y pesqueros de la sardina Sardinops sagax sagax en el Pacifico sororiental. FAO Sinop. Pesca. 1988;1:60.
  74. 74. Christensen V. Fish 501 Mass Balance. 2021 [cited April 2023]. In: University of British Columbia, Institute for Oceans and Fisheries [Internet]. Available from: https://sites.google.com/view/fish501-2021/lectures
  75. 75. Heymans JJ, Coll M, Link JS, Mackinson S, Steenbeek J, Walters C, et al. Best practice in Ecopath with Ecosim food-web models for ecosystem-based management. Ecol. Modell. 2016 Jul 10;331:173–84.
  76. 76. Booth S. Pollutants in the seas around us. In: Pauly D, Zeller D, editors. Global Atlas of Marine Fisheries A critical appraisal of catches and ecosystem impacts. Washington, Covelo, London: Island Press; 2016. p. 152–70.
  77. 77. Coombs AP. Marine mammals and human health in the eastern Bering Sea: Using an ecosystem-based food web model to track PCBs. M.Sc. Thesis. University of British Columbia. 2004. Available from: https://open.library.ubc.ca/soa/cIRcle/collections/ubctheses/831/items/1.0074886
  78. 78. Alava JJ, Cisneros-Montemayor AM, Sumaila UR, Cheung WWL. Projected amplification of food web bioaccumulation of MeHg and PCBs under climate change in the Northeastern Pacific. Sci. Rep. 2018 Dec 1;8(1). pmid:30194394
  79. 79. McMullen K. The Galápagos penguin as the “canary in the coal mine” for microplastics research in the Galápagos Marine Reserve & plastic pollution perceptions in Ecuadorian mangrove communities. M.Sc. Thesis. University of British Columbia. 2023. Available from: http://hdl.handle.net/2429/84664
  80. 80. Everaert G, Van Cauwenberghe L, De Rijcke M, Koelmans AA, Mees J, Vandegehuchte M, et al. Risk assessment of microplastics in the ocean: Modelling approach and first conclusions. Environmental Pollution. 2018 Nov 1;242:1930–8. pmid:30061084
  81. 81. Lusher A. Microplastics in the marine environment: distribution, interactions and effects. In: Bergmann M, Gutow L, Klages M, editors. Marine Anthropogenic Litter. Springer; 2015. p. 245–307.
  82. 82. Lebreton L, Egger M, Slat B. A global mass budget for positively buoyant macroplastic debris in the ocean. Sci. Rep. 2019 Dec 1;9(1). pmid:31515537
  83. 83. Jackson S. Do seabird gut sizes and mean retention times reflect adaptation to diet and foraging method?. Physiological Zoology. 1992 May 1;65(3):674–97.
  84. 84. Laugksch RC, Duffy DC. Food transit rates in cape gannets and jackass penguins. Condor. 1986 Feb 1:117–9.
  85. 85. Critchell K, Hoogenboom MO. Effects of microplastic exposure on the body condition and behaviour of planktivorous reef fish (Acanthochromis polyacanthus). PLOS One. 2018 Mar 1;13(3). pmid:29494635
  86. 86. Ory NC, Gallardo C, Lenz M, Thiel M. Capture, swallowing, and egestion of microplastics by a planktivorous juvenile fish. Environmental Pollution. 2018 Sep 1;240:566–73. pmid:29758531
  87. 87. Ohkubo N, Yoneda M, Ito M, Hano T, Kono K. Microplastic uptake and gut retention time in Japanese anchovy (Engraulis japonicus) under laboratory conditions. Mar. Pollut. Bull. 2022 Mar 1;176. pmid:35202882
  88. 88. Murray F, Cowie PR. Plastic contamination in the decapod crustacean Nephrops norvegicus (Linnaeus, 1758). Mar. Pollut. Bull. 2011 Jun;62(6):1207–17. pmid:21497854
  89. 89. Alava JJ. Proposing a Bioaccumulation Metric Criteria Framework for Plastic Particles in Marine Biota and Foodwebs. In Abstracts Book: Society of Environmental Toxicology and Chemistry North America SETAC. In: (Society for Environmental Toxicology and Chemistry) (SETAC) (SETAC SciCon4) Solutions with Respect for Our Community and Environment [Internet]. Portland, Oregon; 2021. p. 25. Available from: https://scicon4.setac.org/wp-content/uploads/2021/11/SciCon4-abstract-book.pdf
  90. 90. Borgå K, Kidd KA, Muir DCG, Berglund O, Conder JM, Gobas FAPC, et al. Trophic magnification factors: Considerations of ecology, ecosystems, and study design. Integr. Environ. Assess. Manag. 2012 Jan;8(1):64–84. pmid:21674770
  91. 91. Conder JM, Gobas FAPC, Borgå K, Muir DCG, Powell DE. Use of trophic magnification factors and related measures to characterize bioaccumulation potential of chemicals. Integr. Environ. Assess. Manag. 2012 Jan;8(1):85–97. pmid:21538835
  92. 92. Walters DM, Jardine TD, Cade BS, Kidd KA, Muir DCG, Leipzig-Scott P. Trophic Magnification of Organic Chemicals: A Global Synthesis. Environ. Sci. Technol. 2016 May 3;50(9):4650–8. pmid:27014905
  93. 93. Alava JJ, Moreno-Báez M, McMullen K, Tekman MB, Barrows APW, Bergmann M, et al. Ecological Impacts of Marine Plastic Pollution, Microplastics’ Foodweb Bioaccumulation Modelling and Global Ocean Footprint: Insights into the Problems, the Management Implications and Coastal Communities Inequities. The Nippon Foundation-Ocean Litter Project (2019–2023). IOF Working Papers 2023 (01), 81 pp., Institute for the Oceans and Fisheries, University of British Columbia, Vancouver, Canada; 2023. Available from: https://oceans.ubc.ca/research/publications/working-papers/
  94. 94. Ma Y fei, You X yi Modelling the accumulation of microplastics through food webs with the example Baiyangdian Lake, China. Science of the Total Environment. 2021 Mar 25;762.
  95. 95. Chagnon C, Thiel M, Antunes J, Ferreira JL, Sobral P, Ory NC. Plastic ingestion and trophic transfer between Easter Island flying fish (Cheilopogon rapanouiensis) and yellowfin tuna (Thunnus albacares) from Rapa Nui (Easter Island). Environmental Pollution. 2018 Dec 1;243:127–33. pmid:30172118
  96. 96. Rochman CM, Brookson C, Bikker J, Djuric N, Earn A, Bucci K, et al. Rethinking microplastics as a diverse contaminant suite. Vol. 38, Environmental Toxicology and Chemistry. Wiley Blackwell; 2019. p. 703–11.
  97. 97. Monclús L, McCann Smith E, Ciesielski TM, Wagner M, Jaspers VLB. Microplastic Ingestion Induces Size-Specific Effects in Japanese Quail. Environmental Science & Technology. 2022 Oct 27;56(22):15902–11. pmid:36302260
  98. 98. Thornton Hampton LM, Brander SM, Coffin S, Cole M, Hermabessiere L, Koelmans AA, et al. Characterizing microplastic hazards: which concentration metrics and particle characteristics are most informative for understanding toxicity in aquatic organisms? Microplastics and Nanoplastics. 2022 Aug 6;2(1).
  99. 99. Ramsperger AFRM, Bergamaschi E, Panizzolo M, Fenoglio I, Barbero F, Peters R, et al. Nano- and microplastics: a comprehensive review on their exposure routes, translocation, and fate in humans. NanoImpact. 2023 Jan 1;29. pmid:36427812
  100. 100. Leslie HA, van Velzen MJM, Brandsma SH, Vethaak AD, Garcia-Vallejo JJ, Lamoree MH. Discovery and quantification of plastic particle pollution in human blood. Environ Int. 2022 May 1;163. pmid:35367073
  101. 101. Ragusa A, Svelato A, Santacroce C, Catalano P, Notarstefano V, Carnevali O, et al. Plasticenta: First evidence of microplastics in human placenta. Environment international. 2021 Jan 1;146:106274. pmid:33395930
  102. 102. Ragusa A, Notarstefano V, Svelato A, Belloni A, Gioacchini G, Blondeel C, et al. Raman Microspectroscopy Detection and Characterisation of Microplastics in Human Breastmilk. Polymers (Basel). 2022 Jul 1;14(13). pmid:35808745
  103. 103. Bakir A, O’Connor IA, Rowland SJ, Hendriks AJ, Thompson RC. Relative importance of microplastics as a pathway for the transfer of hydrophobic organic chemicals to marine life. Environmental Pollution. 2016 Dec;219:56–65. pmid:27661728
  104. 104. Naik RK, Naik MM, D’Costa PM, Shaikh F. Microplastics in ballast water as an emerging source and vector for harmful chemicals, antibiotics, metals, bacterial pathogens and HAB species: A potential risk to the marine environment and human health. Marine Pollution Bulletin. 2019 Dec 1;149:110525. pmid:31470206
  105. 105. Beiras R, Muniategui-Lorenzo S, Rodil R, Tato T, Montes R, López-Ibáñez S, et al. Polyethylene microplastics do not increase bioaccumulation or toxicity of nonylphenol and 4-MBC to marine zooplankton. Science of the Total Environment. 2019 Nov 20;692:1–9. pmid:31326592
  106. 106. Andrady AL. Plastics in the Anthropocene. In: Plastics and the Ocean. Wiley; 2022. p. 1–42.
  107. 107. Haram LE, Carlton JT, Ruiz GM, Maximenko NA. A Plasticene Lexicon. Mar. Pollut. Bull. 2020 Jan 1;150. pmid:31753559
  108. 108. Thompson RC, Olson Y, Mitchell RP, Davis A, Rowland SJ, John AWG, et al. Lost at Sea: Where Is All the Plastic? Science. 2004 May 7;304(5672):838. pmid:15131299
  109. 109. United Nations Environment Programme. UNEP year book 2011: emerging issues in our global environment. 2011 [cited April 2023] [Internet]. Available from: https://wedocs.unep.org/20.500.11822/8276
  110. 110. Anbumani S, Kakkar P. Ecotoxicological effects of microplastics on biota: a review. Vol. 25, Environmental Science and Pollution Research. Springer Verlag; 2018. p. 14373–96.
  111. 111. Yin K, Wang Y, Zhao H, Wang D, Guo M, Mu M, et al. A comparative review of microplastics and nanoplastics: Toxicity hazards on digestive, reproductive and nervous system. Science of The Total Environment. 2021 Jun 20;774:145758.
  112. 112. Alava JJ, McMullen K, Jones J, Barragán-Paladines MJ, Hobbs C, Tirapé A, et al. Multiple anthropogenic stressors in the Galápagos Islands’ complex social–ecological system: Interactions of marine pollution, fishing pressure, and climate change with management recommendations. Integrated Environmental Assessment and Management. John Wiley and Sons Inc; 2022.
  113. 113. Rowlands E, Galloway T, Cole M, Lewis C, Peck V, Thorpe S, et al. The Effects of Combined Ocean Acidification and Nanoplastic Exposures on the Embryonic Development of Antarctic Krill. Front. Mar. Sci. 2021 Aug 3;8.
  114. 114. Vargas H, Lacy RC, Johnson PJ, Steinfurth A, Crawford RJM, Dee Boersma P, et al. Modelling the effect of El Niño on the persistence of small populations: The Galápagos penguin as a case study. Biol. Conserv. 2007 Jun;137(1):138–48.
  115. 115. McMullen K, Vargas FH, Calle P, Alavarado-Cadena O, Pakhomov EA, Alava JJ (2024) Modelling microplastic bioaccumulation and biomagnification potential in the Galápagos penguin ecosystem using Ecopath and Ecosim (EwE) with Ecotracer. PLoS ONE 19(1): e0296788.