<<

Nova Southeastern University NSUWorks

HCNSO Student Capstones HCNSO Student Work

12-6-2019

A Risk Analysis of Microplastic Consumption in Filter Feeders

Sheri Rahman Nova Southeastern University, [email protected]

This document is a product of extensive research conducted at the Nova Southeastern University . For more information on research and degree programs at the NSU , please click here.

Follow this and additional works at: https://nsuworks.nova.edu/cnso_stucap

Part of the Marine Biology Commons, and the Oceanography and Atmospheric Sciences and Meteorology Commons

Share Feedback About This Item

NSUWorks Citation Sheri Rahman. 2019. A Risk Analysis of Microplastic Consumption in Filter Feeders. Capstone. Nova Southeastern University. Retrieved from NSUWorks, . (347) https://nsuworks.nova.edu/cnso_stucap/347.

This Capstone is brought to you by the HCNSO Student Work at NSUWorks. It has been accepted for inclusion in HCNSO Student Capstones by an authorized administrator of NSUWorks. For more information, please contact [email protected]. Capstone of Sheri Rahman

Submitted in Partial Fulfillment of the Requirements for the Degree of

Master of Science M.S. Marine Biology

Nova Southeastern University Halmos College of Natural Sciences and Oceanography

December 2019

Approved: Capstone Committee

Major Professor: Dr. Abigail Renegar

Committee Member: Dr. Bernhard Riegl

This capstone is available at NSUWorks: https://nsuworks.nova.edu/cnso_stucap/347

HALMOS COLLEGE OF NATURAL SCIENCES AND OCEANOGRAPHY

A Risk Analysis of Microplastic Consumption in Filter Feeders

By

Sheri Rahman

Submitted to the faculty of Halmos College of Natural Sciences and Oceanography in partial fulfillment of the requirements for the degree of Master of Science with a specialty in:

Marine Science Concentration: Marine Biology

Nova Southeastern University

1 December 2019 Submitted in partial fulfillment of the requirements for the degree of

Master of Science: Marine Biology

Sheri Rahman Nova Southeastern University Halmos College of Natural Sciences and Oceanography

December 2019

Capstone Committee Approval

______Dr. Abigail Renegar, Major Professor

______Dr. Bernhard Riegl, Committee Member

2

Table of Contents

Abstract 3

I. Introduction 5

II. Statement of Purpose & Objectives 7

III. Materials and Methods

3.1 Data Acquisition 9

3.2 Data Analysis 13

IV. Results and Review

4.1 Microplastic Abundance 17

4.2 Rates 21

4.3 Microplastic Consumption Rates 23

4.4 Characteristics 28

V. Summary and Conclusions

5.1 Overall Risk Assessment 52

5.2 Future Considerations 56

VI. References 57

3 Abstract

Microplastics (plastic particles < 5 mm) pose a serious threat to marine organisms, as researchers have documented such particles in the gut contents of numerous species. In particular, filter feeders are at risk of consuming microplastics because they may accidentally consume the particulates when feeding or they may prey on species that have already consumed them. The goals of this research were to evaluate the risks that different filter feeders face in regards to microplastic consumption through the analysis of the calculated Microplastic Consumption Rates for numerous species of filter feeders. Factors that could potentially affect this risk were also considered, including ocean basin, environment type, salinity, life stage, IUCN status, and filtration technique. Initial analysis showed that body size greatly impacted a species’ risk of microplastic consumption and further tests were completed to evaluate overall microplastic contamination for each species. Microplastic consumption and microplastic contamination values were evaluated and analyzed to determine which filter feeding species were most at risk of experiencing ecological effects from microplastic pollution. From a resource management perspective, this research highlights the filter feeding species most at risk, contributing to the development of more effective plastic waste management policies.

Keywords: microplastics, plastics, filtration, microplastic consumption, microplastic contamination, filter feeding species

4 I. Introduction

More than nine million tons of plastic fibers are produced every year, and microplastics (plastics < 5 mm) are now found in aquatic environments around the globe (Barrows et al. 2018). Plastics were first produced in the 1950s and became popular very quickly due to their durability and low production costs (Lusher et al. 2017). Although they offer many benefits to the average consumer, including lower prices and convenience, plastic materials have become a danger to the environment. When improperly managed, plastic waste is often allowed to reach freshwater and marine environments. There, the material is exposed to the sun’s ultraviolet rays, causing it to degrade slowly (Lusher et al. 2017). This leads to the breakdown of the material and formation of small, microplastic particles, which have become such a prevalent problem today that they are now considered one of the greatest threats to the health of ecosystems and biodiversity on land and in marine and freshwater regions (Barrows et al. 2018, Lusher et al. 2017). Microplastics can generally be categorized as either primary or secondary. Primary microplastics are fibers and beads manufactured to a small size, which are often used in the cosmetic industry. These particles might be used in soaps, shampoos, toothpastes, shaving cream, makeup, bubble bath, and other cosmetic products around the world (Leslie 2014). When consumers rinse off the product and wash it down the drain, these plastics find their way into wastewater. And while effective management facilities will retain a small portion of these microplastics, the rest flow into freshwater or marine environments (Leslie 2014). Secondary microplastics, on the other hand, are produced from the degradation of larger items (Lusher et al. 2017), such as plastic bottles, bags, and other forms of waste. This degradation occurs as a result of exposure to saltwater and ultraviolet sunlight (Lusher et al. 2017). Plastics are known to include a variety of toxins, as they are often comprised of toxic chemicals and various additives that can have adverse effects on the health of marine organisms (Gallo et al. 2018). A variety of chemicals, such as monomers, plasticizers, and flame-retardants, are added to plastics during production (Lusher et al. 2017). The material can also adsorb contaminants like polychlorinated biphenyls (PCBS), polycyclic aromatic hydrocarbons (PAH), and persistent bioaccumulative toxic substances (PBTs) from the

5 surrounding environment. Contaminants accumulate through predator-prey relationships and trophic transfers, potentially leading to adverse health effects, such as increased immune responses, decreased growth, and decreased fecundity (Gallo et al. 2018, Lusher et al. 2017). Due to their popularity, long lifespan, process of degradation, and potential for toxicity, microplastics have become ubiquitous and a persistent pollutant. As such, it is increasingly important to understand their distribution and concentration around the globe (Barrows et al. 2018). In recent years, new research has expanded knowledge in this area, with much of the work being completed by citizen science initiatives (Barrows et al. 2018). A great example is the Global and Gallatin Microplastics Initiative, which launched a massive project that called for environmentally minded citizens who spend time on the water to take water samples and send it to their facilities for processing. The response was enormous, with samples collected from around the globe, encompassing marine and freshwater environments; this initiative has produced a large microplastic concentration dataset that can be used to bridge knowledge gaps (Global & Gallatin Microplastic Initiatives 2018). It is widely known that many species, including filter feeders, consume microplastics as previous studies have found such particles in the stomachs and guts of various organisms (Cole et al. 2013, Taylor et al. 2016, Wieczorek et al. 2018). Even some of the smallest species, like , bivalve larvae, and decapod larvae, ingest microplastics although the ability to uptake these particles may depend on size (Cole et al. 2013). Species that are larger in size or at higher trophic levels have also been documented interacting with microplastic pollution, whether directly or indirectly (Lusher et al. 2017). Although the direct ingestion of plastic particulates is more commonly studied, trophic transfer might also occur when an organism ingests a prey species that has already consumed the microplastics (Cole et al. 2013, Moore et al. 2001). Evidence even suggests that organisms in the deep have been exposed, as they frequently ingest microplastic fibers (Taylor et al. 2016) Like most other marine species, filter feeding organisms ranging in size and complexity from and jellyfish to whale are also known to consume microplastics either directly if mistaken for food or indirectly as a result of prey consumption of plastic particles or fibers (Cole et al. 2013, Moore et al. 2001). Because filter feeders must filter small food items from the water, such as and , they cannot

6 always be selective and avoid the consumption of other particulates that may also be present (Cole et al. 2013). Some organisms have developed adaptations prevent the consumption of unwanted materials, such as the mesh size of rakers and other anatomical components that can prevent consumption of items larger than a specific size (Roesch et al. 2013). Microplastics can still easily be consumed, however, even if the filter feeder has such adaptations to prevent it. After all, these adaptations were developed over thousands of generations, but microplastics have only been an issue within our oceans for less than a century (Roesch et al. 2013, Lusher et al. 2017). In order to assess the risks that microplastics may pose to filter feeding organisms, it is thus necessary to determine how likely it is that a filter feeder might consume microplastics by considering their filtration rate and the concentration of microplastics in the water. Most recent studies involving the interactions between living organisms and microplastics rely on the use of molluscs or crustaceans, though some may also focus on various species, both in the laboratory and in field observations (Lusher et al. 2017). In almost every niche environment, whether at the sea surface, on beaches, within the water column, or in the deep sea, microplastic uptake occurs among the organisms living there. Seabirds and marine mammals ingest microplastics regularly – an occurrence that can have significant consequences for both the organism and human health (Lusher et al. 2017, Taylor et al. 2016). However, little research has been done to better understand the ecological consequences of this phenomenon, particularly among filter feeders. Though some researchers believe the effects of microplastic consumption would not extend beyond the level of the individual, others have demonstrated that the trend might reduce primary productivity, either directly or indirectly (Lusher et al. 2017). In this study, the risk of microplastic consumption among filter feeders was assessed to bridge such knowledge gaps.

II. Statement of Purpose and Objectives

The goal of this research was to quantitatively assess the risks faced by different filter feeding organisms with regards to the consumption of microplastics based on three primary factors: the abundance of marine plastic debris across geographic locations, as demonstrated

7 by recent studies (Global & Gallatin Microplastic Initiatives 2018, Barrows et al. 2018, Woodall et al. 2014); the location that filter feeding species primarily live and feed; and the different filtration rates utilized by filter feeders. Such information can be used to determine the likelihood filter feeders might consume microplastic particles. In addition to the quantification of the microplastics consumed by these species while feeding, the study also determined if various factors had a significant impact on the estimated consumption of such particles. Perhaps most importantly, the study considered the impact of feeding location on these risks, potentially allowing conservation and waste managers in different areas to fully understand the risks filter feeding species face in their region. Feeding locations – including specific ocean basins, regions, types of environments, and whether the species feeds in marine or freshwater – provided insight into whether the specific variables could potentially impact a species’ risk of microplastic consumption. Some species – such as basking sharks, jellyfish, and others – are globally distributed (Priede et al. 2008, Sims et al. 2003), leading to the expectation that they might be more likely to consume microplastics in areas with greater abundance of these particles than in those with less abundance. Other species are specific to smaller regions. The blue , for instance, is generally found in the North Atlantic, in both the east and west regions of the basin (Boström & Bonsdorff 1997, Wildish & Miyares 1990). The study also considered the vulnerability of each species by considering IUCN Red List status labels (IUCN 2019), as well as the effect of organism age. The filtration technique used by these species was also considered, as distinctive strategies result in differing filtration rates that affect the quantity of microplastics potentially consumed. Filtration technique was expected to have an effect on the quantity of microplastics potentially consumed by filter feeders. Most filter feeders rely on at least one of four primary techniques: ram filtration, suspension feeding, water pumping, and lunge feeding. Ram filtration occurs when a species, such as the whale Rhincodon typus), swims forward slowly with an open mouth to capture food-laden water (Motta et al. 2010). Suspension feeding, however, occurs when an organism like the Pacific (Crassostrea gigas) can capture and extract food items out of the surrounding water as it flows over the (Harris 2008). Water pumping occurs when an organism actively pumps water through the mouth to capture food (Wildish & Miyares 1990), while lunge feeding is frequently seen in

8 large species, such as whales, to capture large quantities of food in one mouthful (Simon et al. 2012). Involving a comparison of multiple representative filter feeders, this study hypothesized that: 1) filter feeders searching for food and feeding in geographic locations with higher microplastic abundance would be more likely to consume plastic; and 2) specific factors, such as filtration technique, could have a significant effect on the risk of microplastic consumption. This study aimed to fill knowledge gaps by analyzing relevant datasets, including filtration rates and microplastic abundance worldwide. Altogether, this valuable information will enable managers to make informed environmental decisions and may aid in the development of more effective resource and waste management policies. Until now, little research has been done to attempt to quantify to what extent different species might consume such particles. After an extensive literature review was performed, a new database of 50 different species of filter feeding organisms was created to facilitate the evaluation of a wide range of filter feeders, from sea worms and bryozoans to whale sharks and fin whales.

II. Materials & Methods

Data Acquisition

This research study required a metadata analysis approach and a risk analysis framework, necessitating the use of various datasets to effectively characterize the risks associated with microplastics (Lusher et al. 2017). To accurately assess these risks, data was collated from a variety of sources, spanning decades of research.

a. Microplastic Abundance Data

The Global & Gallatin Microplastics Initiatives of Adventure Scientists conducted microplastic pollution surveys in aquatic environments around the globe from 2013 to 2017 (2018). This made it possible to assess and analyze where microplastics typically accumulate

9 geographically. After the collection of 2,677 surface water samples in four years, this dataset demonstrates the ubiquity of microplastics in marine and freshwater environments worldwide (Global & Gallatin Microplastics Initiative, 2018). The datasets provided from this research project included 1,394 samples of marine water and 1,009 samples of freshwater (Global & Gallatin Microplastics Initiative, 2018). Samples were taken from a broad range of water sources including coastal regions and open ocean areas of all ocean basins within the marine water dataset. In general, these data points only include surface water because all samples were obtained within the first 50 meters (Global & Gallatin Microplastics Initiative, 2018). Microplastic abundance data was also confirmed with a study conducted by Kanhai et al. (2017). The researchers collated data from previously conducted studies to review microplastic abundance in various locations. They included data for each of the ocean basins, including the region from which samples were taken and the method used to collect water samples (Kanhai et al. 2017). Although this dataset was not directly used in the statistical analysis and calculations within this paper, it was useful in confirming the validity of the mean microplastic abundances determined in the Global & Gallatin Microplastics Initiative project.

b. Filter Feeder Species Selection

Next, datasets illustrating filtration rates for specific representative species was acquired. Because these studies typically focus on one species at a time, data points were gathered individually and collated for further analysis. It was necessary to acquire data for a large variety of filter feeding species, including cnidarians, sponges, bivalves, whales, and fish, to accurately represent the diversity of such organisms. Because no filter feeder database currently exists in an easily accessible manner, one had to be created. Filtration rates for 50 species were collected from 44 published research papers (Table 1).

10 c. Species Characteristics

Other types of characteristic data were also collected for each filter feeding species because this information was necessary to determine which factors have a potentially significant effect on microplastic consumption in different species. This required a more in- depth review of literature for each of the 50 species. The information was included in the filter feeder database to allow for the tracking and analysis of each characteristic. These traits were: feeding locations and distribution, IUCN Red List Status, filtration techniques, whether the species lives in marine or freshwater areas, and whether the species tends to feed in coastal or open ocean areas. To obtain data for all these characteristics, the process entailed a review of an additional 190 papers (Table 1).

Table 1. The different species reviewed in this paper, as well as all the sources from which filtration data and other characteristics were drawn.

Species Sources of Data Motta et al. (2010); Duffy (2002); Heyman et al. (2001); de la Parra Venegas et al. (Rhincondon typus) (2011); Taylor (2006); Graham et al. (2005) (Cetorhinus maximus) Sims (1999); Skomal et al. (2004); Sims et al. (2003); Priede & Miller (2008) Blue Wildish et al. (1990); Bostrom & Bonsdorf (1996); Kotta & Orav (2001); Riisgard (Mytilus edulis) (1991) Jellyfish () Linnaeus (1758); Segura-Puertas et al. (2009); Oleson (1995) Bowhead whales Simon et al. (2009); Goldbogen et al. (2017); Laidre et al. (2007); Wursig et al. (1989); (Balaena mysticetus) Moore et al. (2010); Ashjian et al. (2010); Schick & Urban (2000) Humpback whale Simon et al. (2012); Clapham (2018); D'vincent (1985); Goldbogen et al. (2008); Hain (Megaptera novaeangliae) et al (1982) Blue whales Doniol-Valcroze et al. (2011); Goldbogen et al. (2011); Acevedo-Gutierrez (2002); (Balaenoptera musculus) Watkins & Schevill (1979); Fiedler et al. (1998); Gill et al. (2011); Gill (2002); Fuller & Clark (1936); Prokopchuk & Sentyabox (2006); Speirs et al. (2006); Aksnes & (Calanus finmarchicus) Magnusen (1979); Marshall & Nicholls (1934) Atlantic (Scomber scombrus) Sutherland et al. (1995) Langoy et al. (2012); Overholtz & Keith (2011) Antarctic minke whale Friedlaender et al. (2014); Thiele et al. (2004); Ohsumi et al. (1970); Goldbogen et al. (Balaenoptera bonaerensis) (2017); Tamura & Konishi (2009); (Crassostrea gigas) Qiu et al. (2015); Gerdes (1982); Harris (2008); Fey et al. (2010); Cognie et al. (2006) Bochdansky & Deibel (1998); Gorsky et al. (1982); Tomita et al. (2019); Sato et al. (Oikopleura dioica) (2001); Shelbourne (1953); Hopcroft & Roff (1995) Silver (Hypophthalmichthys molitrix) Zhao et al. (2011); Lazarro (1987) Divi et al. (2018); Paig-Tran et al. (2013); Paig-Tran et al. (2011); Dewar et al. (2008); (Manta birostis) Braun et al. (2014); Stewart et al. (2016) Pelagic Tunicate (Pegea confederata ) Harbison & Gilmer (1976); Harbison & Campenot (1979); Sutherland et al. (2010) Fin whales Goldbogen et al. (2010); Vikingsson et al. (2009); Mizroch et al. (1984); Monestiez et (Balaenoptera physalus) al. (2004); Panigada et al. (1999)

11 Glass Leys et al. (2011); Kahn et al. (2015); Yahel et al. (2007); Austin et al. (2007); Buhl- (Aphrocallistes vastus) Mortensen (2009) Cockle (Cardium edule) Riisgard et al. (2002); Richardson et al. (1993); Kater et al. (2006) Soft-shell Riisgard et al. (2002); Strasser (1999); Snelgrove et al. (1999); Seitz et al. (2001); (Mya arenia) Armonies & Reise (2003) Atlantic (Brevoortis tyrannus) Durbin & Durbin (1975); Love et al. (2006); Buchheister et al. (2016) Mysid (Rhopalophthalmus Jerling & Wooldridge (1994); Webb et al. (1997); Wooldridge (1986); Shlachler & terranatalis) Wooldridge (1995) Mysid shrimp Jerling & Wooldridge (1994); Webb et al. (1997); Paul & Calliari (2017); Froneman (Mesopodopsis wooldridgei) (2001) Burrowing shrimp (Upogebia deltaura) Lindahl & Baden (1997); Christiansen (2000); Tunberg (1985); Howe et al. (2004) Antarctic (Euphausia Boyd et al. (1984); Atkinson et al. (2008); Hill et al. (2013); Clarke & Tyler (2008); superba) Schmidt et al. (2014) Porcelain Crab (Porcellana longicornis) Achituv & Pedrotti (1999); Lance (1964); Werding et al. (2003) Ocean Quahog () Winter (1969); Cargnelli et al. (1999); Witbaard & Bergman (2003) Wrinkled Rockborer (Hiatella arctica) Ali (1970); Gordillo (2001); Sejr et al. (2002); Wlodarska-Kowalczuk (2007) Bay (Pecten irradians) Chipman & Hawkins (1954); MacKenzie (2008); Smith et al. (1988) Orange (Ptilosarcus gurneyi) Best (1988); Stone (2006) Feather star (Oligometra serripinna) Leonard et al. (1988)); Holland et al. (1991); Tay et al. (2016); Hellal (2012) Manila Clam (Ruditapes philippinarum) Nakamura (2001); Velez et al. (2015); Dang et al. (2010); Lewis et al. (2007) Yesso scallop (Patinopecten yessoensis) Yamamoto (1968); Sato et al. (2004); Silina (1996) Spaghetti Bryozoan Bullivant (1967); Minchin (2012); Amat & Tempera (2009); McCann et al. (2015); (Zoobotryon verticillatum) Jebakumar et al. (2017) Bryozoan (Electra pylosa) Riisgard & Manriquez (1997); Nikulina et al. (2007); Hermansen et al. (2007) Bryozoan (Conopeum reticulum) Riisgard & Manriquez (1997) Bryozoan (Celleporella hyalina) Riisgard & Manriquez (1997); Hermansen et al. (2007) Sea vase Randlov & Riisgard (1979); Runnstrom (1936); Havenhand (1991); Therriault & (Ciona intestinalis) Herborg (2008) Sea squirt (Ascidella aspersa) Randlov & Riisgard (1979); Schmidt (1983); Chebbi et al. (2010); Mastrotaro (2008) Polychaete worm (Myxicola infundibulum) Dales (1957); Gotshall (2005); Greathead et al. (2011) Peacock worm (Sabella pavonina) Dales (1957); Greathead et al. (2011); Murray et al. (2011) Keel worm (Pomatoceros Dales (1957); Kupriyanova & Badyaev (1998); Ponti et al. (2002); Southward (1957); triqueter) Ekaratne et al. (1982) Polychaete worm (Hydroides norvegica) Dales (1957); Moen (2006); Southward (1957) Sinistral spiral tubeworm (Spirorbis borealis) Dales (1957); O'Connor & Lamont (1978) Polychaete worm (Salmacina dysteri) Dales (1957); Isaac (1974); Eldredge & Smith (2001); Nishi (1992); Parnell (2001) Breadcrumb sponge Riisgard et al. (1993); Hansen et al. (1995); Vethaak et al. (1982); Forester (1979); (Halichondria panicea) Peattie & Hoare (1981) (Abramis brama) van den Berg (1993); Lammens (1986); Kuparinen et al. (2014); Lyons & Lucas (2002)

12 White Bream (Blicca bjoerkna) van den Berg (1993); Lammens (1986) Roach (Rutilus rutilus) van den Berg (1993) Gizzard shad (Dorosoma cepedianum) van den Berg (1993); Drenner et al. (1984); Wuellner et al. (2008) North Atlantic van der Hoop et al. (2019); Baumgartner & Mate (2005); Baumgartner et al. (2003); (Eubalaena glacialis) Baumgartner & Mate (2003)

Data Analysis

The sources reviewed to obtain these data points provided a more in-depth look at the risk each filter feeder faced regarding their consumption of microplastics. With the creation of the database, a risk assessment framework was used to evaluate the likelihood of adverse ecological effects as a result of filter feeder exposure to microplastics. More than 2,000 data points illustrated global microplastic abundance. To simplify calculations, the values were categorized based on larger regions, encompassing specific locations as well as surrounding areas (Global & Gallatin Microplastics Initiative 2018). The relationship between microplastic abundance and geographic location were assessed based on two distinct factors: ocean basin and environment. The ocean basin variable had five fixed levels. Microplastic abundance in various ocean basins were not normally distributed (p<0.05, Shapiro-Wilkes) or homoschedastic (p=0.01, Bartlett’s). A fixed factor One Way ANOVA of log-transformed data was thus used in the assessment. The environment variable had only two fixed levels, and the data were not normally distributed (p<0.05, Shapiro- Wilkes) or homoschedastic (p=0.005, Bartlett’s). A two-tailed, two sample t-test of log- transformed data was used to assess any significant differences between microplastic abundance and environment. Using the collated filtration rates, the Microplastic Consumption Rate (MCR) was then calculated in order to quantify how many microplastic particles are likely to be consumed by each filter feeding species. Calculation of the MCR required that filtration rates for each species be converted to mL s-1. Additionally, each filter feeder needed an assigned estimated feeding location based on its known geographic distribution. Once filter feeders were assigned at least one location, the corresponding mean microplastic abundance for that region was multiplied by the species’ filtration rate, as in the following equation:

13 Equation 1. Determination of the Microplastic Consumption Rate

Filtration Rate (mL/s) * Mean Microplastic Abundance (particles/mL) = Microplastic Consumption Rate (particles/s)

After calculating MCR values, the mean, median, and mode of microplastic consumption were determined. In total, 68 data points were considered for the 50 different filter feeding species, as some species were assessed at multiple feeding locations or life stages. The significance of various factors, including salinity, IUCN status, filtration technique, life stage, ocean basin, and environment, in relation to MCR were assessed using R software. In determining the influence of salinity on MCR values, the factor was defined as a categorical variable, indicating whether each species feeds in marine or freshwater areas. Raw and transformed data were not normally distributed (p<0.05, Shapiro-Wilkes) or homoschedastic (p<0.05, Bartlett’s), so the non-parametric two-tailed Mann-Whitney Wilcoxon test was used to compare MCR values between salinity levels. Also defined as a categorical variable, the IUCN Red List Status included five levels at which the different filter feeding species were labeled: Not Evaluated, Least Concern, Near Threatened, Vulnerable, or Endangered (IUCN 2019). Raw and transformed data were not normally distributed (p<0.05, Shapiro-Wilkes) or homoschedastic (p<0.05, Bartlett’s), so the non-parametric Kruskal-Wallis test was used to determine the significance of the relationship because the categorical factor had more than three levels in this assessment. Filtration technique was defined as another categorical variable with four levels: lunging, suspension, pumping, and ram. The levels were determined through a review of literature, which indicated the typical techniques used by study species. Raw and transformed data were not normally distributed (p<0.05, Shapiro-Wilkes) or homoschedastic (p<0.05, Bartlett’s), so the non-parametric Kruskal-Wallis test was used to assess the significance of filtration technique. This factor was also further reviewed to consider which species is most likely to experience higher MCR values at each of the four techniques by only analyzing each level at a time.

14 To determine whether filter feeders experienced significant differences in MCR as adults or juveniles, only those species that included data at different life stages were considered. In this case, life stage was defined as simply being Adult or Juvenile. Data were normally distributed (p>0.05, Shapiro-Wilkes) and homoschedastic (p=0.075, Bartlett’s). Thus, a two-tailed, two sample t-test was used because this factor only had two levels. To consider if the ocean basin influenced the MCR, the feeding locations for filter feeding species were estimated. For example, whale sharks are known to feed in coastal areas near Mexico (Motta et al. 2010, de la Parra Venegas et al. 2011). For this reason, the mean microplastic abundance value for Pacific Central America Coastal was used to calculate the whale shark’s microplastic consumption rate. Like a few other species, whale sharks were assessed at multiple locations. Because they might also feed near the of New Zealand (Duffy 2002), they were also assessed using the mean microplastic abundance values from the Pacific West Coastal category.

Although basking sharks are known to be a global species, they are often found in waters near Scotland and thus their feeding location was estimated to be around the Atlantic East Coastal category for the purpose of this research study (Priede & Miller 2008, Sims et al. 2003, Skomal et al. 2004). Blue mussels were analyzed in both Atlantic NW Coastal and Atlantic NE Coastal regions (Bostrom & Bonsdorf 1996, Kotta & Orav 2001, Riisgard 1991, Wildish & Miyares 1990), and bowhead whales were also estimated to feed in multiple locations: Atlantic NW Coastal and Pacific SE Alaska Coastal (Ashjian et al. 2010, Goldbogen et al. 2017, Laidre et al. 2007, Moore et al. 2010, Schick & Urban 2000, Simon et al. 2009). Continuing through the database of 50 species, feeding locations for all filter feeders were estimated, and some relied on the analysis of more than one region. Raw and transformed data for the ocean basin variable were not normally distributed (p<0.05, Shapiro-Wilkes) or homoschedastic (p<0.05, Bartlett’s), so the non-parametric Kruskal-Wallis test was used to assess the significance of feeding location. Additionally, an unbalanced Two Way ANOVA was completed to analyze both filtration technique and ocean basin to determine if any significant interactions occurred between these variables. Because the same dataset was used, parametric assumptions were once again not met and a non- parametric Kruskal-Wallis test was used.

15 Finally, the significance of the environment in relation to MCR was assessed. Defined as another categorical variable, the environment indicated whether species fed in coastal or open ocean locations. Data were not normally distributed (p<0.05, Shapiro-Wilkes) or homoschedastic (p<0.05, Bartlett’s). Thus, a two-tailed, two sample t-test of log-transformed data was used. Analysis of the MCR values for each species indicated the possibility that organism size played a key role in a species’ risk of microplastic consumption. To determine the nature of this relationship, further data regarding average bodyweight for each review species was collated. MCR values were then normalized as an MCR-to-bodyweight ratio with the values reported in units of particles/s/kg. After the data was normalized, analytical tests were run once again to determine if bodyweight affected the significant differences in MCR values for each of the six factors considered. In the consideration of Normalized Microplastic Consumption Rates (or NMCR), data for salinity was once again not normal (p<0.05, Shapiro-Wilkes) or homoschedastic (p<0.05, Bartlett’s). To analyze this factor, the non-parametric Mann-Whitney Wilcoxon test was used. Data for IUCN status was also found to be not normal (p<0.05, Shapiro-Wilkes) or homoschedastic (p<0.05, Bartlett’s). The normalized data for this variable, then, required a non-parametric Kruskal-Wallis test for analysis. Analysis of the normalized data for filtration technique indicated that data was still not normal (p<0.05, Shapiro-Wilkes) or homoschedastic (p<0.05, Bartlett’s). The non- parametric Kruskal-Wallis ANOVA, thus, was used for analysis. The life stage factor once again required analysis of only data from relevant species. Normalized data were found to be normal (p>0.05, Shapiro-Wilkes) and homoschedastic (p=0.614, Barlett’s), so analysis required a two-tailed two sample t-test. After the normalized data for the ocean basin variable was transformed, however, the data was found to be normal (p>0.05, Shapiro-Wilkes) and homoschedastic (p=0.08, Bartlett’s). For this variable, a One Way ANOVA could be used for the analysis. Similarly normalized data for the environment variable was found to be normal (p>0.05, Shapiro- Wilkes) and homoschedastic (p=0.497, Bartlett’s) after a log transformation. Thus, analysis required the use of a two-tailed, two sample t-test.

16 The analysis also showed that there was a possible interaction between the two factors, filtration technique and ocean basin. To investigate further, an unbalanced two-way ANOVA was run to determine if interactions between the two factors had any significant effect on MCR. As previously noted, data was not normal or homoschedastic for either variable, so non-parametric tests were used in the analysis. Additionally, NMCR data was also considered and the test was run a second time to determine if taking body weight into consideration impacted the results. As a final step in this project, the different filter feeding species were then divided into groups based on one factor found to be significant in the analysis: filtration technique. Once they were grouped as such, mean MCR and NMCR values were graphed to determine which species within each sub-category were most at risk of microplastic consumption or contamination. While it would be useful to determine if filtration technique had a significant impact on MCR and NMCR values for each of the subcategories, it was not possible to test with a One Way ANOVA because there were not enough data points for each species.

IV. Results and Discussion

Microplastic Abundance

Because several of the ocean basins are so large, spanning across different nations and localities, the microplastic abundance data were first categorized to make cross- referencing with filtration rates simpler. (Table 2). For example, coastal samples from the Atlantic Ocean were considered part of Caribbean, Gulf of Mexico, Mediterranean, Northwest (including North America, Bermuda, and Canada), Northeast (including United Kingdom, Europe, and Africa), and the South Atlantic regions. The Pacific coastal data points were also categorically divided into Central America, Gulf of Alaska, SE Alaska, SE Asia, West (including Australia, New Zealand, Niue, and Beveridge), and East (North and South America, Hawaii, Mexico, and Canada) regions (Figure 1). To analyze the factors affecting microplastic abundance, data was then further grouped by ocean basin.

17 The raw data shows that samples from the open ocean typically contain the greatest abundance of particles compared to coastal water samples (Table 2). The mean value of 54.57 ± 16.07 particles/L was found for open ocean samples from the Arctic basin, while coastal values of the same basin were 23.87 ± 6.46 particles/L (Table 2). Of the open ocean samples, highest mean values of microplastic abundance were found for the Arctic, Pacific (18.42 ± 3.47 particles/L), Atlantic (17.96 ± 1.22 particles/L), Southern (17.5 ± 1.22 particles/L), and Indian (16.87 ± 10.22 particles/L) oceans, respectively (Table 2) (Global & Gallatin Microplastic Initiative, 2018).

Table 2. Microplastic abundance data calculated at ocean basins and environments. Ocean Coastal or Mean Microplastic Abundance Basin Regional Sea Open Ocean (particles / L)(± SE) Arctic Coastal 23.8708 (± 6.4608) Arctic Open Ocean 54.5680 (± 16.0698) Atlantic Caribbean Coastal 9.9372 (± 3.5674) Atlantic Gulf of Mexico Coastal 3.0120 (± 1.4593) Atlantic Mediterranean Coastal 2.1180 (± 0.8149) Atlantic NW (America, bermuda, canada) Coastal 5.8342 (± 0.6392) Atlantic NE (UK, Europe, Africa) Coastal 1.9975 (± 0.4088) Atlantic South Coastal 2.2262 (± 0.6452) Atlantic Caribbean Open Ocean 5.840 (± 3.7176) Atlantic Mediterranean Open Ocean 9.0476 (± 1.3877) Atlantic Open Ocean 18.0176 (± 1.2235) Indian Coastal 2.9480 (± 0.5434) Indian Open Ocean 16.8722 (± 10.2184) Pacific Central America Coastal 4.3898 (± 0.7205) Pacific Gulf of Alaska Coastal 8.1858 (± 1.8316) Pacific SE Alaska Coastal 5.6129 (± 1.2623) Pacific SE Asia Coastal 5.3268 (± 1.4358) West (Australia, New Zealand, Niue, Pacific Beveridge) Coastal 1.0850 (± 0.2545) East (America, Mexico, Canada, Pacific Hawaii, S. America) Coastal 2.7056 (± 0.8773) Pacific Central America Open Ocean 3.1231 (± 1.7898) Pacific SE Asia Open Ocean 19.0741 (± 12.2426) Pacific Open Ocean 18.4176 (± 3.4670) Southern Coastal 15.29 (± 8.7241) Southern Open Ocean 17.5 Freshwater North America 1.1493 (± 0.0858) Freshwater Europe 1.5720 (± 0.3808)

18 The analysis showed that ocean basins experience significant difference in mean microplastic abundance (p = 0.0432, F4,19=0.0323, One-Way ANOVA). The Arctic had the highest mean and the Indian had the lowest mean compared to other sample locations (Figure 2). Furthermore, post-hoc analysis (Multiple Comparisons) indicated that microplastic abundance in the Arctic Ocean was significantly higher than abundance data in the Atlantic Ocean, while the Indian, Pacific, and Southern Oceans were not significantly different from each other. The Open Ocean and coastal environments were also compared to determine if this factor affected microplastic abundance. Open ocean environments had a significantly higher mean microplastic abundance compared to coastal samples (p = 0.005, t = -3.22, two-tailed two sample t-test) (Figure 3). When considering the results from these analyses, it is important to know that the raw data was not evenly distributed throughout the global ocean. Rather, very few samples were taken from the Arctic and Southern Oceans, likely because data was collected on a volunteer basis and fewer individuals were able to visit these locations. The Atlantic and Pacific Oceans, however, had a far greater quantity of data points available. Such an unbalanced distribution could affect the reliability of these results.

Microplastic Abundance in Categorized Sample Locations 80 70 60 ) 50

SE 40 ± 30 20 10

0

(particles/L (particles/L Mean Mean Microplastic Abundance

Location

Figure 1. Mean microplastic abundance (particles/L ± SE) for each of the sample locations as categorized for further analysis.

19 Microplastic Abundance in Different Ocean Basins 40 35 a

) 30 SE

± 25 ab 20

15 b

(particles/L (particles/L ab 10 ab

5 Mean Mean Microplastic Abundance 0 Arctic Atlantic Indian Pacific Basin or Water Source

Figure 2. Mean microplastic abundance (particles/L ± SE) for each of the marine sample locations.

Microplastic Abundance in Each Environment 20 18

16 )

SE 14

± 12 10 8

6 (particles/L (particles/L 4

2 Mean Mean Microplastic Abundance 0 Coastal Open Ocean Type of Environment

Figure 3. Mean microplastic abundance (particles/L ± SE) found in the two different types of environment, Coastal and Open Ocean.

20 Filtration Rates

In general, species of the smallest sizes, such as oyster larvae, bryozoans, copepods, seaworms, and , filter the least amount of water (Table 3). Larvae of Pacific , for example, filter 1.39x10-6 mL water/second. Much larger – and therefore, stronger and faster – species, however, tend to filter greater quantities of water. Fin whales, for example, can filter volumes of water as large as 9.75x106 mL/second (Table 3). Following this filtration rate would be that of bowhead whales (3.02x10-6 mL/s), then North Atlantic right whales (1.39x106 mL/s), humpback whales (7.0x105 mL/s), basking sharks (1.20x10-5 mL/s), and whale sharks (9.06x104 mL/s).

Table 3. Mean minimum and maximum filtration rates for each species (mL/s). Filtration Rate Filtration Rate Maximum Species Name (Scientific) Minimum (mL/s) (mL/s) Source Whale Shark (Rhincondon typus) 9.06E+04 1.71E+05 Motta et al. (2010) Basking Shark (Cetorhinus maximus) 1.20E+05 1.20E+05 Sims (1999) Blue mussels (Mytilus edulis) 6 38 Wildish et al. (1990); Riisgard et al. (2002) Jellyfish (Aurelia aurita) 2.17E-03 7.56E-02 Oleson (1995) Bowhead whales (Balaena mysticetus) 3.20E+06 3.20E+06 Simon et al. (2009); Goldbogen et al. (2017) Humpback whale (Megaptera novaeangliae) 7.00E+05 7.00E+05 Simon et al. (2012) Blue whales (Balaenoptera Doniol-Valcroze et al. (2011); Goldbogen et musculus) 5.85E+02 5.85E+02 al. (2011) Copepod (Calanus finmarchicus) 5.21E-05 5.21E-05 Fuller & Clark (1936) Atlantic mackerel (Scomber scombrus) 26.67 51.67 Sutherland et al. (1995) Antarctic minke whale (Balaenoptera bonaerensis) 5.36E+04 5.36E+04 Friedlaender et al. (2014) Oyster LARVAE (Crassostrea gigas) 1.39E-06 1.39E-06 Qiu et al. (2015) Pacific Oyster ADULT (C. gigas) - smaller size 0.108 0.108 Gerdes (1982) Pacific Oyster ADULT (C. gigas) - larger size 0.288 0.288 Gerdes (1982) Tunicate (Oikopleura dioica) 2.31E-04 2.31E-04 Bochdansky & Deibel (1998) (Hypophthalmichthys molitrix Val.) 9.58 10.42 Zhao et al. 2011 Manta Ray (Manta birostis) 1.51E+04 1.51E+04 Divi et al. (2018); Paig-Tran et al. (2013) Pelagic Tunicate (Pegea confederata ) 6.17E-03 7.77E-02 Harbison & Gilmer (1976) Fin whales (Balaenoptera physalus) 9.75E+06 9.75E+06 Goldbogen et al. (2010) Glass sponge (Aphrocallistes vastus) 17.25 1.73E+01 Leys et al. (2011)

21 Cockle (Cardium edule) 0.111 1.03 Riisgard et al. (2002) Soft-shell clam (Mya arenia) 0.333 1.056 Riisgard et al. (2002) (Brevoortis tyrannus) 41.67 87.83 Durbin & Durbin (1975) Mysid shrimp (Rhopalophthalmus terranatalis) ADULTS 2.36E-03 2.36E-03 Jerling & Wooldridge (1994) (Rhopalophthalmus terranatalis) JUVENILES 3.75E-03 3.75E-03 Jerling & Wooldridge (1994) Mysid shrimp (Mesopodopsis wooldridgei) ADULTS 9.50E-03 9.50E-03 Jerling & Wooldridge (1994) (Mesopodopsis wooldridgei) JUVENILES 5.17E-03 5.17E-03 Jerling & Wooldridge (1994) Burrowing shrimp (Upogebia deltaura) 0.972 9.72E-01 Lindahl & Baden (1997) (Euphausia superba) 0.125 0.125 Boyd et al. (1984) Porcelain Crab (Porcellana longicornis) 3.94E-02 7.42E-02 Achituv & Pedrotti (1999) Ocean Quahog (Arctica islandica) 0.555 1.14E+00 Winter (1969) Wrinkled Rockborer (Hiatella arctica) 1.53E-03 9.47E-03 Ali (1970) Bay Scallop (Pecten irradians) 0.906 4.089 Chipman & Hawkins (1954) Orange Sea Pen (Ptilosarcus gurneyi) 100 1000 Best (1988) Feather star (Oligometra serripinna) 68 111.6 Leonard et al. (1988) Manila Clam (Tapes philippinarum) 2.78E-02 0.278 Nakamura (2001); Hosokawa (1988) Yesso scallop (Patinopecten yessoensis) 0.694 1.1 Yamamoto (1968) Spaghetti Bryozoan (Zoobotryon verticillatum) 4.22E-05 2.92E-04 Bullivant (1967) Bryozoan (Electra pylosa) 6.94E-05 7.78E-05 Riisgard & Manriquez (1997) Bryozoan (Conopeum reticulum) 4.72E-05 5.56E-05 Riisgard & Manriquez (1997) Bryozoan (Celleporella hyalina) 3.33E-05 4.17E-05 Riisgard & Manriquez (1997) Sea vase (Ciona intestinalis) 0.05 0.2 Randlov & Riisgard (1979) Sea squirt (Ascidella aspersa) 0.067 0.333 Randlov & Riisgard (1979) Polychaete worm (Myxicola infundibulum) 7.90E-02 7.94E-02 Dales (1957) Peacock worm (Sabella pavonina) 2.03E-02 2.03E-02 Dales (1957) Keel worm () 7.50E-03 7.50E-03 Dales (1957) Polychaete worm (Hydroides norvegica) 3.10E-03 3.10E-03 Dales (1957) Sinistral spiral tubeworm (Spirorbis borealis) 6.39E-05 6.39E-05 Dales (1957) Polychaete worm (Salmacina dysteri) 8.06E-04 8.06E-04 Dales (1957) Breadcrumb sponge (Halichondria panicea) 7.17E-02 7.17E-02 Riisgard et al. (1993) Common Bream (Abramis brama) 7.6389 7.6389 van den Berg (1993) White Bream (Blicca bjoerkna) 6.389 6.389 van den Berg (1993) Roach (Rutilus rutilus) 9.833 9.833 van den Berg (1993) Gizzard shad (Dorosoma van den Berg (1993) & Drenner et al. cepedianum) 20.833 20.833 (1984) North Atlantic Right Whale (Eubalaena glacialis) 1.39E+06 1.39E+06 van der Hoop et al. (2019)

22 Microplastic Consumption Rates

With filtration rates and microplastic abundance data collated, the estimated microplastic abundance was determined for each species’ location categories and the MCR (in particles/s) was calculated (Table 4). The species with the highest maximum microplastic consumption rate (MCR) while feeding was the in the Pacific Open Ocean (1.79x105 particles/s). Among the species reviewed in this paper, the lowest maximum MCR occurred in larvae of Pacific Oysters (1.51x10-9 particles/s). However, among only the adults (and thus, excluding juveniles), the lowest maximum MCR occurred in the bryozoan, E. pylosa (8.48x10-8 particles/s). The minimum mean MCR was found to be 1.51x10 -09 particles/s, while the maximum mean MCR was found to be 6.235x103. The variance and standard deviation values, 9.21x108 and 3.03x104, respectively, further indicated that the data was very spread out. The MCR data indicated a strong increasing trend with increasing body weight (Figure 4), suggesting that an organism’s size played a significant role in MCR values. For this reason, the data was further analyzed to create a new dataset of with these values reported as a MCR-to-bodyweight ratio in units of particles/s/kg (Table 4). While MCR values provide information regarding a species’ risk of microplastic consumption, these Normalized MCR (or NMCR) values provide information regarding a species’ risk of microplastic contamination because the values are reported in terms of body weight. Analysis of the data showed that the pelagic tunicate (P. confederata) actually experiences the highest risk of microplastic contamination, as it had a NMCR value of 5.17x104 particles/s/kg. The bryozoan (E. pylosa), however, experiences the lowest risk of microplastic contamination with a NMCR value of 1.88x10 -07 particles/s/kg (Figure 5). Analysis of the relationship between body size and MCR showed that bodyweight does have a significant relationship with MCR values (p=1.92x10-11, z=6.71, tau=0.563, Non-parametric correlation). This relationship indicates that smaller organisms are more at risk of microplastic contamination, as they appear to consume larger quantities of microplastics on a per-kg basis.

23 Table 4. Mean minimum and maximum MCR (particles/s) for each species. Mean Normalized Microplastic Abundance MCR Species Sample Locations Mean MCR (particles/s) (particles/s/kg) Antarctic Krill (Euphausia superba) Southern, Coastal 1.91E-03 3.35E+00 Antarctic minke whale (Balaenoptera 9.90E-02 bonaerensis) Southern, Open Ocean 9.38E+02 Atlantic mackerel (Scomber scombrus) Atlantic, NW, Coastal 3.01E-01 9.77E-01 Atlantic mackerel (Scomber scombrus) Atlantic, NE, Coastal 1.03E-01 3.34E-01 Atlantic, Mediterranean, 1.10E-01 Atlantic mackerel (Scomber scombrus) Coastal 3.39E-02 Atlantic menhaden (Brevoortis 8.53E-01 tyrannus) Atlantic, NW, Coastal 5.12E-01 Basking Shark (Cetorhinus maximus) Pacific, East, Coastal 3.25E+02 8.13E-02 Bay Scallop (Pecten irradians) Atlantic, NW, Coastal 2.38E-02 9.52E-02 Blue mussels (Mytilus edulis) Atlantic, NW, Coastal 2.22E-01 3.36E+01 Blue mussels (Mytilus edulis) Atlantic, NE, Coastal 7.59E-02 1.15E+01 Blue whales (Balaenoptera musculus) Pacific, East, Coastal 1.59E+00 2.00E-05 Bowhead whales (Balaena mysticetus) Atlantic, NW, Coastal 1.87E+04 2.49E-01 Bowhead whales (Balaena mysticetus) Pacific, SE Alaska, Coastal 1.80E+04 2.40E-01 Breadcrumb sponge (Halichondria 2.20E-03 panicea) Atlantic, NE, Coastal 1.43E-04 Breadcrumb sponge (Halichondria 1.20E-03 panicea) Pacific, West, Coastal 7.82E-05 Bryozoan (Celleporella hyalina) Atlantic, NW, Coastal 2.43E-07 5.40E-07 Bryozoan (Celleporella hyalina) Pacific, East, Coastal 1.13E-07 2.41E-07 Bryozoan (Conopeum reticulum) Atlantic, NE, Coastal 1.11E-07 2.46E-07 Bryozoan (Electra pylosa) Pacific, West, Coastal 8.48E-08 1.88E-07 Atlantic, Mediterranean, 1.03E+00 Burrowing shrimp (Upogebia deltaura) Coastal 2.06E-03 Cockle (Cardium edule) Atlantic, NE, Coastal 2.06E-03 2.58E-01 Common Bream (Abramis brama) Freshwater, Europe 1.17E-02 1.95E-03 Copepod (Calanus finmarchicus) Atlantic Open Ocean (surface) 9.37E-07 2.86E+00 Feather star (Oligometra serripinna) Pacific, SE Asia, Coastal 2.13E+00 7.10E+02 Fin whales (Balaenoptera physalus) Pacific Open Ocean 1.79E+05 3.58E+00 Fin whales (Balaenoptera physalus) Atlantic, Open Ocean 1.76E+05 3.52E+00 Gizzard shad (Dorosoma cepedianum) Freshwater, North America 2.40E-02 1.26E-02 Pacific, Gulf of Alaska, 1.57E-01 Glass sponge (Aphrocallistes vastus) Coastal 1.41E-01 Humpback whale (Megaptera 1.41E-01 novaeangliae) Atlantic, NW, Coastal 4.08E+03 Humpback whale (Megaptera 2.63E-02 novaeangliae) Pacific, West, Coastal 7.63E+02 Jellyfish (Aurelia aurita) Atlantic, NW, Coastal 4.41E-04 6.35E-03 Keel worm (Pomatoceros triqueter) Atlantic, NE, Coastal 1.50E-05 3.75E+00 Keel worm (Pomatoceros triqueter) Arctic, Coastal 1.79E-04 4.48E+01 Manila Clam (Tapes philippinarum) Indian, Coastal 8.20E-04 7.13E-02

24 Manta Ray (Manta birostis) Atlantic, Caribbean, Coastal 1.50E+02 9.10E-02 Manta Ray (Manta birostis) Pacific, SE Asia, Coastal 8.05E+01 4.88E-02 Mysid shrimp (Mesopodopsis 1.18E+02 wooldridgei) ADULTS Indian, Coastal 2.80E-05 Mysid shrimp (Mesopodopsis 3.19E+02 wooldridgei) JUVENILES Indian, Coastal 1.53E-05 Mysid shrimp (Rhopalophthalmus 2.32E+00 terranatalis) ADULTS Indian, Coastal 6.96E-06 Mysid shrimp (Rhopalophthalmus 1.11E+01 terranatalis) JUVENILES Indian, Coastal 1.11E-05 North Atlantic Right Whale (Eubalaena 1.09E+00 glacialis) Atlantic, Open Ocean 2.50E+04 Ocean Quahog (Arctica islandica) Atlantic, NW, Coastal 6.64E-03 2.92E-02 Orange Sea Pen (Ptilosarcus gurneyi) Pacific, East, Coastal 2.71E+00 1.81E+01 Pacific Oyster ADULT (C. gigas) - 7.85E-04 larger size Pacific, West, Coastal 3.14E-04 Pacific Oyster ADULT (C. gigas) - 5.10E-04 smaller size Pacific, West, Coastal 1.18E-04 Pacific Oyster LARVAE (Crassostrea 7.55E-02 gigas) Pacific, West, Coastal 1.51E-09 Atlantic, Mediterranean, 1.08E+01 Peacock worm (Sabella pavonina) Coastal 4.30E-05 Pelagic Tunicate (Pegea confederata ) Atlantic, NE, Coastal 1.55E-04 5.17E+04 Atlantic, Mediterranean, 1.64E+00 Polychaete worm (Hydroides norvegica) Coastal 6.57E-06 Polychaete worm (Myxicola 3.58E+02 infundibulum) Atlantic, Open Ocean 1.43E-03 Polychaete worm (Myxicola 3.63E+02 infundibulum) Pacific Open Ocean 1.45E-03 Polychaete worm (Salmacina dysteri) Pacific Open Ocean 1.48E-05 3.70E+00 Porcelain Crab (Porcellana longicornis) Atlantic, NE, Coastal 1.48E-04 5.92E+00 Roach (Rutilus rutilus) Freshwater, Europe 1.50E-02 8.33E-03 Sea squirt (Ascidella aspersa) Atlantic, NE, Coastal 6.65E-04 4.43E-03 Sea vase (Ciona intestinalis) Atlantic, Open Ocean 3.60E-03 1.76E-04 Silver Carp (Hypophthalmichthys 4.68E-04 molitrix Val.) Freshwater, Asia 2.34E-02 Sinistral spiral tubeworm (Spirorbis 3.20E-02 borealis) Atlantic, NE, Coastal 1.28E-07 Sinistral spiral tubeworm (Spirorbis 4.33E-02 borealis) Pacific, East, Coastal 1.73E-07 Soft-shell clam (Mya arenia) Atlantic, NW, Coastal 6.16E-03 1.81E-01 Soft-shell clam (Mya arenia) Atlantic, NE, Coastal 2.11E-03 6.20E-02 Spaghetti Bryozoan (Zoobotryon 4.83E-04 verticillatum) Atlantic, Caribbean, Coastal 2.90E-06 Tunicate (Oikopleura dioica) Atlantic, NE, Coastal 4.62E-07 1.54E+02 Pacific, Central America, 2.20E-02 Whale Shark (Rhincondon typus) Coastal 7.51E+02 Whale Shark (Rhincondon typus) Pacific, West, Coastal 1.86E+02 5.50E-03 White Bream (Blicca bjoerkna) Freshwater, Europe 9.78E-03 9.78E-03 Wrinkled Rockborer (Hiatella arctica) Atlantic, South, Coastal 2.13E-05 1.42E-02 Yesso scallop (Patinopecten yessoensis) Pacific, West, Coastal 1.20E-03 9.23E-04

25

Figure 4. Mean microplastic consumption rate (particles/s ± SE) for each filter feeding species in order of bodysize. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

26 Figure 5. Mean normalized microplastic consumption rate (particles/s/kg ± SE) for each filter feeder in order of body size. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

27 Filter Feeder Characteristics

a. Salinity

Salinity was analyzed to allow for comparison between marine and freshwater species. Only a few freshwater species were considered in this review, including silver carp (Hypophthalamichthys molitrix), common bream (Abramis brama), white bream (Blicca bjoerkna), roach (Rutilus rutilus), and gizzard shad (Dorosoma cepedianum). The remaining

49 were marine species (Table 4).

Table 4. The salinity type (Marine or Freshwater) to which each species belongs.

Species Salinity Whale Shark (Rhincondon typus) Marine Basking Shark (Cetorhinus maximus) Marine Blue mussels (Mytilus edulis) Marine Jellyfish (Aurelia aurita) Marine Bowhead whales (Balaena mysticetus) Marine Humpback whale (Megaptera novaeangliae) Marine Blue whales (Balaenoptera musculus) Marine Copepod (Calanus finmarchicus) Marine Atlantic mackerel (Scomber scombrus) Marine Antarctic minke whale (Balaenoptera bonaerensis) Marine Pacific Oyster LARVAE (Crassostrea gigas) Marine Pacific Oyster ADULT (C. gigas) - smaller size Marine Pacific Oyster ADULT (C. gigas) - larger size Marine Tunicate (Oikopleura dioica) Marine Silver Carp (Hypophthalmichthys molitrix Val.) Freshwater Manta Ray (Manta birostis) Marine Pelagic Tunicate (Pegea confederata ) Marine Fin whales (Balaenoptera physalus) Marine Glass sponge (Aphrocallistes vastus) Marine Cockle (Cardium edule) Marine Soft-shell clam (Mya arenia) Marine Atlantic menhaden (Brevoortis tyrannus) Marine Mysid shrimp (Rhopalophthalmus terranatalis) ADULTS Marine Mysid shrimp (Rhopalophthalmus terranatalis) JUVENILES Marine Mysid shrimp (Mesopodopsis wooldridgei) ADULTS Marine Mysid shrimp (Mesopodopsis wooldridgei) JUVENILES Marine

28 Burrowing shrimp (Upogebia deltaura) Marine Antarctic Krill (Euphausia superba) Marine Porcelain Crab (Porcellana longicornis) Marine Ocean Quahog (Arctica islandica) Marine Wrinkled Rockborer (Hiatella arctica) Marine Bay Scallop (Pecten irradians) Marine Orange Sea Pen (Ptilosarcus gurneyi) Marine Feather star (Oligometra serripinna) Marine Manila Clam (Tapes philippinarum) Marine Yesso scallop (Patinopecten yessoensis) Marine Spaghetti Bryozoan (Zoobotryon verticillatum) Marine Bryozoan (Electra pylosa) Marine Bryozoan (Conopeum reticulum) Marine Bryozoan (Celleporella hyalina) Marine Sea vase (Ciona intestinalis) Marine Sea squirt (Ascidella aspersa) Marine Polychaete worm (Myxicola infundibulum) Marine Peacock worm (Sabella pavonina) Marine Keel worm (Pomatoceros triqueter) Marine Polychaete worm (Hydroides norvegica) Marine Sinistral spiral tubeworm (Spirorbis borealis) Marine Polychaete worm (Salmacina dysteri) Marine Breadcrumb sponge (Halichondria panicea) Marine Common Bream (Abramis brama) Freshwater White Bream (Blicca bjoerkna) Freshwater Roach (Rutilus rutilus) Freshwater Gizzard shad (Dorosoma cepedianum) Freshwater North Atlantic Right Whale (Eubalaena glacialis) Marine

Salinity was assessed to have no significant relationship with microplastic consumption rates (p = 0.3719, w = 196, Mann-Whitney Wilcoxon test). Although marine species had a higher mean MCR than freshwater species (Figure 4), the difference was not significant. The differences seen are likely due to the Marine outliers, which are above 1.5x105 particles/s. However, it is important to note that the differences are likely results of the few data points collected for freshwater species. Only five species out of the 50 live in freshwater environments, and the sample size can easily impact the reliability and precision of the non-parametric test used to analyze the relationship.

29 When taking bodyweight into consideration, however, significant differences did occur in NMCR values between salinity levels (p = 0.026, w=62, Mann-Whitney Wilcoxon test). Marine species had a significantly higher NMCR than freshwater species (Figure 7). This suggests that species in marine water would experience higher risks of microplastic contamination.

MCR & Salinity 1.00E+05 9.00E+04 8.00E+04

) 7.00E+04 SE

± 6.00E+04 5.00E+04

4.00E+04 Mean Mean MCR

3.00E+04 (particles/s (particles/s 2.00E+04 1.00E+04 8.97E+04 1.68E-02 0.00E+00 Marine Freshwater Salinity

Figure 6. The calculated MCR (particles/s ± SE) at both types of salinity, freshwater and marine.

NMCR & Salinity

1.00E+03 )

SE 5.00E+02 ±

6.63E-03 0.00E+00

Marine Freshwater Mean Mean NMCR -5.00E+02

(particles/s/kg (particles/s/kg 7.11E-02

-1.00E+03 Salinity

Figure 7. The calculated NMCR (particles/s/kg ± SE) at both types of salinity, freshwater and marine.

30 b. IUCN Red List Status

To effectively determine which species are most at risk of experiencing harmful ecological impacts from microplastics, the IUCN Red List status for each species was also collected (Table 5) (IUCN 2019). The IUCN generally categorizes species based on the vulnerability status, including labels that range from Least Concern to Vulnerable, Endangered, Critically Endangered, Extinct in the Wild, and Extinct. For those species on which very little data has been collected, the organization generates the default label, Not Evaluated (IUCN 2019). Of the 50 different filter feeding species, most have yet to be evaluated and are thus given the label “NE.” Of the evaluated species, only three were considered endangered (EN) –whale sharks (Rhindondon typus), blue whales (Balaenoptera musculus), and North Atlantic right whales (Eubalaena glacialis) (IUCN 2019). Those labeled Least Concern (LC) include bowhead whales (Balaena mysticetus), humpback whales (Megaptera novaeangliae), Atlantic mackerel (Scomber scombrus), Atlantic menhaden (Brevoortis tyrannus), Antarctic krill (Euphausia superba), common bream (Abramis brama), white bream (Blicca bjoerkna), roach (Rutilus rutilus), and gizzard shad (Dorosoma cepedianum) (IUCN 2019). A few were considered to be vulnerable (VU) species, including the basking shark (Cetorhinus maximus), the manta ray (Manta birostis), and fin whales (Balaenoptera physalus). Antarctic minke whales (Balaenoptera bonaerensis) and silver carp (Hypophthalmichthys molitrix) were considered near threatened (NT), while only the whale sharks (Rhincondon typus), blue whales (Balaenoptera musculus), and North Atlantic right whales (Eubalaena glacialis) were considered endangered (EN) species (IUCN 2019).

31 Table 5. IUCN Red List Status of study species (IUCN 2019).

Species IUCN Red List Status Whale Shark (Rhincondon typus) EN Basking Shark (Cetorhinus maximus) VU Blue mussels (Mytilus edulis) Not Evaluated Jellyfish (Aurelia aurita) Not Evaluated Bowhead whales (Balaena mysticetus) LC Humpback whale (Megaptera novaeangliae) LC Blue whales (Balaenoptera musculus) EN Copepod (Calanus finmarchicus) Not Evaluated Atlantic mackerel (Scomber scombrus) LC Antarctic minke whale (Balaenoptera bonaerensis) NT Pacific Oyster JUVENILES (Crassostrea gigas) Not Evaluated Pacific Oyster ADULT (C. gigas) - smaller size Not Evaluated Pacific Oyster ADULT (C. gigas) - larger size Not Evaluated Tunicate (Oikopleura dioica) Not Evaluated Silver Carp (Hypophthalmichthys molitrix Val.) NT Manta Ray (Manta birostis) VU Pelagic Tunicate (Pegea confederata ) Not Evaluated Fin whales (Balaenoptera physalus) VU Glass sponge (Aphrocallistes vastus) Not Evaluated Cockle (Cardium edule) Not Evaluated Soft-shell clam (Mya arenia) Not Evaluated Atlantic menhaden (Brevoortis tyrannus) LC Mysid shrimp (Rhopalophthalmus terranatalis) ADULTS Not Evaluated Mysid shrimp (Rhopalophthalmus terranatalis) JUVENILES Not Evaluated Mysid shrimp (Mesopodopsis wooldridgei) ADULTS Not Evaluated Mysid shrimp (Mesopodopsis wooldridgei) JUVENILES Not Evaluated Burrowing shrimp (Upogebia deltaura) Not Evaluated Antarctic Krill (Euphausia superba) LC Porcelain Crab (Porcellana longicornis) Not Evaluated Ocean Quahog (Arctica islandica) Not Evaluated Wrinkled Rockborer (Hiatella arctica) Not Evaluated Bay Scallop (Pecten irradians) Not Evaluated Ornage Sea Pen (Ptilosarcus gurneyi) Not Evaluated Feather star (Oligometra serripinna) Not Evaluated Manila Clam (Tapes philippinarum) Not Evaluated Yesso scallop (Patinopecten yessoensis) Not Evaluated Spaghetti Bryozoan (Zoobotryon verticillatum) Not Evaluated Bryozoan (Electra pylosa) Not Evaluated Bryozoan (Conopeum reticulum) Not Evaluated

32 Bryozoan (Celleporella hyalina) Not Evaluated Sea vase (Ciona intestinalis) Not Evaluated Sea squirt (Ascidella aspersa) Not Evaluated Polychaete worm (Myxicola infundibulum) Not Evaluated Peacock worm (Sabella pavonina) Not Evaluated Keel worm (Pomatoceros triqueter) Not Evaluated Polychaete worm (Hydroides norvegica) Not Evaluated Sinistral spiral tubeworm (Spirorbis borealis) Not Evaluated Polychaete worm (Salmacina dysteri) Not Evaluated Breadcrumb sponge (Halichondria panicea) Not Evaluated Common Bream (Abramis brama) LC White Bream (Blicca bjoerkna) LC Roach (Rutilus rutilus) LC Gizzard shad (Dorosoma cepedianum) LC North Atlantic Right Whale (Eubalaena glacialis) EN

The quantitative analysis of the IUCN status for each species allowed determination of whether a significant relationship exists with the corresponding MCR for the species. VU-labeled species had a significantly higher mean MCR than any other status (p<<0.000, χ2= 38.195, df = 4, Kruskal-Wallis test). Not Evaluated species, however, had a significantly lower mean MCR than any other status (Figure 8). The significance of these results indicate that IUCN status could be used as a potential indicator of a species’ microplastic consumption risk. However, similar biological characteristics must be met when drawing similar conclusions for other species. Post-hoc analysis (Multiple Comparisons) found that species categorized as endangered, least concern, near threatened, and vulnerable were not significantly different from each other. Species categorized as endangered and near threatened were also not significantly different from each other. Once bodyweight was taken into account and NMCR values were calculated, the statistical tests were re-run to determine if this impacted the results. It was determined that no significant differences occurred in NMCR values between IUCN statuses (p = 0.51, χ2=3.29, df=4, Kruskal-Wallis Test). This suggests that IUCN status does not indicate whether a species is at more or less risk of microplastic contamination (Figure 9).

33 MCR & IUCN Status 1.40E+05 1.20E+05 a ) a

SE 1.00E+05 ± 8.00E+04 7.11E+04 6.00E+04 a a Mean Mean MCR 4.00E+04

(particles/s (particles/s b ab 2.00E+04 6.49E+03 3.19E+03 1.21E-01 4.69E+02 0.00E+00 EN LC NE NT VU IUCN Status

Figure 8. Mean MCR (particles/s ± SE) for each IUCN Red List Status labels, including EN, LC, NE, NT, and VU.

NMCR & IUCN Status 3.00E+03

) ) 2.40E+03

SE ± 1.80E+03 1.22E+03

1.20E+03 Mean Mean NMCR

6.00E+02 (particles/s/kg (particles/s/kg 4.86E-01 4.97E-02 0.00E+00 2.79E-01 1.46E+00 EN LC NE NT VU IUCN Status

Figure 9. Mean NMCR (particles/s/kg ± SE) for each IUCN Red List Status labels, including EN, LC, NE, NT, and VU.

c. Filtration Technique

Of the 50 marine species reviewed, most relied on at least one of four main techniques: ram filtration, suspension feeding, water pumping (or suction feeding), and lunge feeding. For the purposes of this study, whale sharks, basking sharks, bowhead whales,

34 Atlantic mackerel, manta rays, Atlantic menhaden, and North Atlantic right whales are primarily considered to use ram filtration techniques. Humpback whales, blue whales, Antarctic minke whales, and fin whales typically rely on lunge feeding methods, while blue mussels, copepods, tunicates, and pelagic tunicates use water-pumping methods (Table 6). The remaining species, including jellyfish, Pacific oysters, glass sponges, cockles, soft-shell , porcelain crabs, ocean quahogs, wrinkled rockborers, bay , orange sea pens, feather stars, Manila clams, Yesso scallops, bryozoans, sea vase, sea squirts, polychaete worms, peacock worms, and keel worms, are considered suspension feeders (Table 6). In some cases, a species might be known to use more than one technique, such as whale sharks. Although these gentle giants primarily rely on ram filtration techniques by swimming forward at slow speeds, they have also been documented using an active suction feeding method. To do this, they frequently position themselves vertically just below the water’s surface and use a powerful buccal pump to create a suction, trapping their prey in gill rakers (Heyman et al. 2001). Despite multiple techniques, the mean filtration rate obtained for this review corresponds with the whale shark’s use of ram filtration, and the species is considered a ram filter feeder. Although some species may use highly specialized methods to obtain prey, their overall technique is still considered to fall into one of these four categories. For example, tunicates (Pegea confederata and Oikopleura doica) are considered to use the water pumping technique, accomplishing filtration by creating a “house” and pumping water through it (Bochdansky & Deibel 1998, Tomita et al. 2019). Antarctic krill (Euphausia superba) offer another excellent example, as these suspension feeders frequently create a feeding apparatus with the use of their front legs (Boyd et al. 1984, Clark & Tyler 2014). Manta rays are also unique in that, although they use a ram filtration technique, their specific strategy is known as ricochet filtration (Divi et al. 2018). Despite the unique methods and adaptations these species use in water filtration, an overall assessment required the categorization of their techniques into one of the four primary methods.

35 Table 6. Filtration technique (ram, suspension, lunge, or pumping) used by study species.

Species Filtration Technique Whale Shark (Rhincondon typus) Ram filtration Basking Shark (Cetorhinus maximus) Ram filtration Blue mussels (Mytilus edulis) Suspension Jellyfish (Aurelia aurita) Suspension Bowhead whales (Balaena mysticetus) Ram filtration Humpback whale (Megaptera novaeangliae) Lunge feeding Blue whales (Balaenoptera musculus) Lunge feeding Copepod (Calanus finmarchicus) Suspension feeding Atlantic mackerel (Scomber scombrus) Ram Filtration Antarctic minke whale (Balaenoptera bonaerensis) Lunge feeding Oyster JUVENILES (Crassostrea gigas) Suspension feeding Pacific Oyster ADULT (C. gigas) - smaller size Suspension feeding Pacific Oyster ADULT (C. gigas) - larger size Suspension feeding Tunicate (Oikopleura dioica) Water pumping Silver Carp (Hypophthalmichthys molitrix Val.) Water pumping Manta Ray (Manta birostis) Ram Filtration Pelagic Tunicate (Pegea confederata ) Water pumping Fin whales (Balaenoptera physalus) Lunge feeding Glass sponge (Aphrocallistes vastus) Suspension feeding Cockle (Cardium edule) Suspension feeding Soft-shell clam (Mya arenia) Suspension feeding Atlantic menhaden (Brevoortis tyrannus) Ram Filtration Mysid shrimp (Rhopalophthalmus terranatalis) ADULTS Suspension feeding Mysid shrimp (Rhopalophthalmus terranatalis) JUVENILES Suspension feeding Mysid shrimp (Mesopodopsis wooldridgei) ADULTS Suspension feeding Mysid shrimp (Mesopodopsis wooldridgei) JUVENILES Suspension feeding Burrowing shrimp (Upogebia deltaura) Suspension feeding Antarctic Krill (Euphausia superba) Suspension feeding Porcelain Crab (Porcellana longicornis) Suspension feeding Ocean Quahog (Arctica islandica) Suspension feeding Wrinkled Rockborer (Hiatella arctica) Suspension feeding Bay Scallop (Pecten irradians) Suspension feeding Ornage Sea Pen (Ptilosarcus gurneyi) Suspension feeding Feather star (Oligometra serripinna) Suspension feeding Manila Clam (Tapes philippinarum) Suspension feeding Yesso scallop (Patinopecten yessoensis) Suspension feeding Spaghetti Bryozoan (Zoobotryon verticillatum) Suspension feeding Bryozoan (Electra pylosa) Suspension feeding Bryozoan (Conopeum reticulum) Suspension feeding

36 Bryozoan (Celleporella hyalina) Suspension feeding Sea vase (Ciona intestinalis) Suspension feeding Sea squirt (Ascidella aspersa) Suspension feeding Polychaete worm (Myxicola infundibulum) Suspension feeding Peacock worm (Sabella pavonina) Suspension feeding Keel worm (Pomatoceros triqueter) Suspension feeding Polychaete worm (Hydroides norvegica) Suspension feeding Sinistral spiral tubeworm (Spirorbis borealis) Suspension feeding Polychaete worm (Salmacina dysteri) Suspension feeding Breadcrumb sponge (Halichondria panicea) Suspension feeding Common Bream (Abramis brama) Ram or Suction White Bream (Blicca bjoerkna) Ram or Suction Roach (Rutilus rutilus) Ram or Suction Gizzard shad (Dorosoma cepedianum) Ram or Suction North Atlantic Right Whale (Eubalaena glacialis) Ram Filtration

In assessing the effect of filtration technique on Microplastic Consumption Rates, significant differences occurred in MCR values at the different levels of filtration technique. (p = 2.015e-08, χ2 = 38.694, df = 3, Kruskal-Wallis test This is likely due to the tendency that such species, including whales, are often much larger than other filter feeders and can therefore filter far greater quantities of particulates from water (Figure 10). Furthermore, a post-hoc analysis (Multiple Comparisons) determined that species using lunge and ram filtration techniques had significantly higher MCR values than the others. Species that relied on water pumping and ram filtration were not significantly different from each other, while those that relied on suspension feeding and water pumping were had significantly lower MCR values than the techniques. However, the results appeared to change when taking bodyweight into account. No significant differences occurred in NMCR values between filtration techniques (p = 0.185, χ2=4.83, df=3, Kruskal-Wallis ANOVA). Thus, the different filtration techniques used by filter feeders are not associated with higher risks of microplastic contamination (Figure 11). The filter feeding species were also further separated into groups based on their filtration techniques to determine which species of each category faced the greatest risks. When bodyweight was not taken into account, this analysis showed that of the lunge feeding species, fin whales experience the highest MCR values, while blue whales (B. musculus) experienced the lowest MCR values (Figure 12). Taking bodyweight into consideration did

37 not cause any change to this result in regards to NMCR values (Figure 13). Of the water pumping filter feeders, blue mussels (M. edulis) had the highest MCR values but copepods (C. finmarchicus) had the lowest MCR values when bodyweight was not considered (Figure 14). When considering bodyweight, however, the pelagic tunicate (P. confederata) had the highest NMCR values, while the silver carp (H. molitrix) had the lowest NMCR values (Figure 15). Of ram filter feeders, the North Atlantic right whale (E. glacialis) had the highest MCR values, while the white bream (B. bjoerkna) had the lowest MCR values (Figure 16). When taking bodyweight into account, the North Atlantic right whale still has the highest NMCR values, but the common bream (A. brama) has the lowest NMCR values (Figure 17). Finally, among suspension feeders, the orange sea pen (P. gurneyi) had the highest MCR values but Pacific oyster larvae (C. gigas) had the lowest MCR values when bodyweight was not considered (Figure 18). When organism size was considered, the feather star (O. serripinna) had the highest NMCR values but the bryozoan (E. pylosa) had the lowest NMCR values (Figure 19). While this analysis provides new insight into the filtration technique categories, it was not possible to test for significant differences because only one data point existed for each species.

MCR & Filtration Technique a 1.00E+05

8.00E+04 )

SE 6.01E+04 ± 6.00E+04

4.00E+04

Mean Mean MCR ab

(particles/s (particles/s bc c 2.00E+04 3.95E+03 5.35E-02 1.26E-01 0.00E+00 Ram Pumping Suspension Lunge Filtration Technique

Figure 10. Calculated MCR (particles/s ± SE) for each of the four types of filtration technique: lunge feeding, water pumping, ram filtration, and suspension feeding.

38 NMCR & Filtration Technique 2.00E+04

) 1.60E+04

SE ± 1.20E+04 8.64E+03

8.00E+03 Mean Mean NMCR

(particles/s/kg (particles/s/kg 4.00E+03

4.93E+01 0.00E+00 2.58E-01 1.23E+00 Ram Pumping Suspension Lunge Filtration Technique

Figure 11. Calculated NMCR (particles/s/kg ± SE) for each of the four types of filtration technique: lunge feeding, water pumping, ram filtration, and suspension feeding.

MCR for Lunge Feeders 2.00E+05

) 1.60E+05 SE ± 1.20E+05

MCR 8.00E+04

(particles/s (particles/s 4.00E+04

0.00E+00 Balaenoptera Balaenoptera Balaenoptera Balaenoptera Megaptera Megaptera bonaerensis musculus physalus physalus novaeangliae novaeangliae Lunge Feeding Species

Figure 12. Calculated MCR (particles/s ± SE) for lunge feeders. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

39 NMCR for Lunge Feeders 4

3.5 )

SE 3 ± 2.5 2 NMCR NMCR 1.5

1 (particles/s/kg (particles/s/kg 0.5 0 Balaenoptera Balaenoptera Balaenoptera Balaenoptera Megaptera Megaptera bonaerensis musculus physalus physalus novaeangliae novaeangliae Lunge Feeding Species

Figure 13. Calculated NMCR (particles/s/kg ± SE) for lunge feeders. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

MCR for Pumping Feeders 2.50E-01

) 2.00E-01

SE ± 1.50E-01

MCR 1.00E-01

5.00E-02 (particles/s (particles/s

0.00E+00

Water Pumping Species

Figure 14. Calculated MCR (particles/s ± SE) for pumping feeders. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

40 NMCR for Pumping Feeders

60000 )

SE 50000 ± 40000 30000 NMCR NMCR 20000

10000 (particles/s/kg (particles/s/kg 0

Water Pumping Species

Figure 15. Calculated NMCR (particles/s/kg ± SE) for pumping feeders. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

MCR for Ram Feeders

2.40E+04 )

SE 1.80E+04 ±

1.20E+04 MCR

6.00E+03 (particles/s (particles/s

0.00E+00

Rutilusrutilus

Mantabirostis Mantabirostis

Blicca bjoerkna

Abramis brama

Rhincondontypus Rhincondontypus

Scomber scombrus Scomber scombrus Scomber scombrus

Balaenamysticetus Balaenamysticetus

Eubalaenaglacialis

Brevoortistyrannus

Cetorhinusmaximus Dorosomacepedianum Ram Feeding Species

Figure 16. Calculated MCR (particles/s ± SE) for ram feeders. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

41

NMCR for Ram Feeders

1.2 )

SE 1 ± 0.8 0.6

NMCR NMCR 0.4 0.2

(particles/s/kg (particles/s/kg 0

Dorosoma

cepedianum

Rutilusrutilus

Mantabirostis Mantabirostis

Blicca bjoerkna

Abramis brama

Rhincondontypus Rhincondontypus

Scomber scombrus Scomber scombrus Scomber scombrus

Balaenamysticetus Balaenamysticetus

Eubalaenaglacialis

Brevoortistyrannus Cetorhinusmaximus Lunge Feeding Species

Figure 17. Calculated NMCR (particles/s/kg ± SE) for ram feeders. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

MCR for Suspension Feeders 3.00E+00

) 2.50E+00 SE

± 2.00E+00

1.50E+00 MCR 1.00E+00

5.00E-01 (particles/s (particles/s

0.00E+00

Mya arenia Mya arenia Mya

Mesopodopsis… Mesopodopsis…

Aurelia aurita Aurelia

Electra pylosa Electra

Cardium edule Cardium

Hiatella arctica Hiatella

Pecten irradians Pecten

Arctica islandica Arctica

Ascidella aspersa Ascidella

Sabella pavonina Sabella

Ciona intestinalis Ciona

Spirorbis borealis Spirorbis borealis Spirorbis

Salmacina dysteri Salmacina

Upogebia deltaura Upogebia

Rhopalophthalmus… Rhopalophthalmus…

Ptilosarcus gurneyi Ptilosarcus

Celleporella hyalina Celleporella hyalina Celleporella

Hydroides norvegica Hydroides

Conopeum reticulum Conopeum

Tapes philippinarum Tapes

Aphrocallistes vastus Aphrocallistes

Halichondria panicea Halichondria panicea Halichondria

Oligometra serripinna Oligometra

Porcellana longicornis Porcellana

Pomatoceros triqueter Pomatoceros

Myxicola infundibulum Myxicola infundibulum Myxicola

Patinopecten yessoensis Patinopecten

Zoobotryon verticillatum Zoobotryon

Crassostrea gigas_juvenile Crassostrea

Crassostrea gigas_adult_large Crassostrea Crassostrea gigas_adult_small Crassostrea Suspension Feeding Species

Figure 18. Calculated MCR (particles/s ± SE) for suspension feeders. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

42 NMCR for Suspension Feeders 800

) 700 SE

± 600 500 400

NMCR NMCR 300 200

100 (particles/s/kg (particles/s/kg

0

Mya arenia Mya Mya arenia Mya

Mesopodopsis…

Aurelia aurita Aurelia

Electra pylosa Electra

Cardium edule Cardium

Hiatella arctica Hiatella

Pecten irradians Pecten

Arctica islandica Arctica

Ascidella aspersa Ascidella

Sabella pavonina Sabella

Ciona intestinalis Ciona

Spirorbis borealis Spirorbis Spirorbis borealis Spirorbis

Salmacina dysteri Salmacina

Upogebia deltaura Upogebia

Rhopalophthalmus… Rhopalophthalmus…

Ptilosarcus gurneyi Ptilosarcus

Celleporella hyalina Celleporella hyalina Celleporella

Hydroides norvegica Hydroides

Conopeum reticulum Conopeum

Tapes philippinarum Tapes

Aphrocallistes vastus Aphrocallistes

Halichondria panicea Halichondria panicea Halichondria

Oligometra serripinna Oligometra

Porcellana longicornis Porcellana

Pomatoceros triqueter Pomatoceros

Myxicola infundibulum Myxicola infundibulum Myxicola

Patinopecten yessoensis Patinopecten

Zoobotryon verticillatum Zoobotryon

Crassostrea gigas_juvenile Crassostrea

Crassostrea gigas_adult_large Crassostrea

Crassostreat gigas_adult_small Crassostreat Mesopodopsis wooldridgei_adult Mesopodopsis Suspension Feeding Species

Figure 19. Calculated NMCR (particles/s/kg ± SE) for suspension feeders. Multiple columns indicate data at different geographic locations for a single species, as described in Table 8.

d. Life Stage

Filtration rates of three species were also considered at different life stages, as either Adults or Juveniles. Only data for species that included filtration rates at both life stages were considered in the review (Table 7). These were the , Crassostrea gigas, and two species of mysid shrimp, Rhopalophthalmus terranatalis and Mesopodopsis wooldridgei (Wildish et al. 1990, Riisgard et al. 2002, Jerling & Wooldridge 1994).

43 Table 7. The life stages (Adult or Juvenile) at which they each of the three species were considered and their MCR (particles/s) Mean Microplastic Consumption Species LifeStage Rate (particles/s) Blue mussel (Crassostrea gigas) Juvenile 1.51E-09 Blue mussel (Crassostrea gigas) – small size Adult 1.18E-04 Blue mussel (Crassostrea gigas) – large size Adult 3.14E-04 Mysid shrimp (Rhopalophthalmus terranatalis) Adult 2.36E-03 Mysid shrimp (Rhopalophthalmus terranatalis) Juvenile 3.75E-03 Mysid shrimp (Mesopodopsis wooldridgei) Adult 9.50E-03 Mysid shrimp (Mesopodopsis wooldridgei) Juvenile 5.17E-03

Life stage was assessed to have no significant relationship with microplastic consumption rates (p = 0.2209, t = 1.5382, Welch two-sample t-test). The adult group has a higher mean MCR (Figure 20), though the difference is not significant. This difference likely occurs because only a few species were considered at both life stages, constricting the sample size. It might also occur as a result of the juveniles being less efficient at water filtration. Similarly, significant differences still did not occur in NMCR values between the life stages when bodyweight was taken into account (p = 0.336, t = 1.06, df = 4.997, two-tailed two-sample t-test). This result suggests both adults and juveniles organisms experience equal risks of microplastic contamination (Figure 21).

MCR & Life Stage 1.20E+04

) 1.00E+04 SE

± 8.00E+03 6.52E+03 6.00E+03

Mean Mean MCR 4.00E+03

(particles/s (particles/s 2.00E+03 8.77E-06 0.00E+00 Adult Juvenile Life Stage

Figure 20. The calculated MCR value (particles/s ± SE) at both types of life stages, Adult or Juvenile.

44 NMCR & Life Stage 1.80E+03

) 1.50E+03 SE

± 1.20E+03

9.00E+02 8.23E+02

Mean Mean NMCR 6.00E+02

(particles/s/kg (particles/s/kg 3.00E+02 1.10E+02 0.00E+00 Adult Juvenile Life Stage

Figure 21. The calculated NMCR value (particles/s/kg ± SE) at both types of life stages, Adult or Juvenile.

e. Ocean Basin

Locations were generalized so that filter feeders could be placed into one of several categories (Table 8): Pacific, Central America, Coastal (PCC); Pacific, West, Coastal (PWC); Pacific, East, Coastal (PEC); Pacific, Southeast Alaska, Coastal (PSAC); Pacific, Southeast Asia, Coastal (PSC); Pacific, Open Ocean (PO); Pacific, Gulf of Alaska, Coastal (PGC); Atlantic, Northwest, Coastal (ANWC); Atlantic, Northeast, Coastal (ANEC); Atlantic, Open Ocean (AO); Atlantic, Mediterranean, Coastal (AMC); Atlantic, Caribbean, Coastal (ACC); Atlantic, South, Coastal (ASC); Southern, Open Ocean (SO); Indian, Coastal (IC); Arctic, Coastal (AC); Freshwater, Asia (FA); Freshwater, Europe (FE); and Freshwater, North America (FNA). Some species required than one category, particularly for those that are globally distributed. For example, the whale shark is commonly found in PCC waters as well as PWC waters; blue mussels can be found in ANWC waters and in ANEC waters; bowhead whales are found in ANWC and PSAC waters; humpback whales are found in ANWC and PWC locations; Atlantic mackerel can be found in ANWC, ANEC, and AMC waters; manta rays are found in ACC and PSC locations; fin whales are located in PO and AO waters; soft-shell

45 clams are located in ANWC and ANEC waters; bryozoans (Celleporella hyaline) are generally found in ANWC and PEC waters; polychaete worms (Myxicola infundibulum) are found in AO and PO waters; keel worms are in ANEC and AC waters; sinistral spinal tubeworms are found in ANEC and PEC locations; and, finally, breadcrumb sponges have been documented in ANEC and PWC waters.

Table 8. Estimated geographic distribution and sampling locations (indicating ocean basin and environment) for each species. Cross-referenced with the microplastic abundance data to be used in calculation of MCR. Species Estimated Geographic Distribution Ocean Basin Environment (Microplastic Abundance Sample Locations) Whale Shark (Rhincondon typus) Pacific, Central America, Coastal Pacific Coastal

Whale Shark (Rhincondon typus) Pacific, West, Coastal Pacific Coastal

Basking Shark (Cetorhinus maximus) Pacific, East, Coastal Pacific Coastal

Blue mussels (Mytilus edulis) Atlantic, NW, Coastal Atlantic Coastal

Blue mussels (Mytilus edulis) Atlantic, NE, Coastal Atlantic Coastal

Jellyfish (Aurelia aurita) Atlantic, NW, Coastal Atlantic Coastal Bowhead whales (Balaena mysticetus) Atlantic, NW, Coastal Atlantic Coastal

Bowhead whales (Balaena mysticetus) Pacific, SE Alaska, Coastal Pacific Coastal

Humpback whale (Megaptera Atlantic, NW, Coastal Atlantic Coastal novaeangliae) Humpback whale (Megaptera Pacific, West, Coastal Pacific Coastal novaeangliae) Blue whales (Balaenoptera musculus) Pacific, East, Coastal Pacific Coastal

Copepod (Calanus finmarchicus) Atlantic Open Ocean (surface) Atlantic Open Ocean

Atlantic mackerel (Scomber scombrus) Atlantic, NW, Coastal Atlantic Coastal

Atlantic mackerel (Scomber scombrus) Atlantic, NE, Coastal Atlantic Coastal

Atlantic mackerel (Scomber scombrus) Atlantic, Mediterranean, Coastal Atlantic Coastal

Antarctic minke whale (Balaenoptera Southern, Open Ocean Southern Open Ocean bonaerensis)

46 Pacific Oyster JUVENILES (Crassostrea Pacific, West, Coastal Pacific Coastal gigas) Pacific Oyster ADULT (C. gigas) - Pacific, West, Coastal Pacific Coastal smaller size Pacific Oyster ADULT (C. gigas) - Pacific, West, Coastal Pacific Coastal larger size

Tunicate (Oikopleura dioica) Atlantic, NE, Coastal Atlantic Coastal

Silver Carp (Hypophthalmichthys Freshwater, Asia Freshwater Coastal molitrix) Manta Ray (Manta birostis) Atlantic, Caribbean, Coastal Atlantic Coastal

Manta Ray (Manta birostis) Pacific, SE Asia, Coastal Pacific Coastal

Pelagic Tunicate (Pegea confederata ) Atlantic, NE, Coastal Atlantic Coastal

Fin whales (Balaenoptera physalus) Pacific Open Ocean Pacific Open Ocean

Fin whales (Balaenoptera physalus) Atlantic, Open Ocean Atlantic Open Ocean

Glass sponge (Aphrocallistes vastus) Pacific, Gulf of Alaska, Coastal Pacific Coastal

Cockle (Cardium edule) Atlantic, NE, Coastal Atlantic Coastal

Soft-shell clam (Mya arenia) Atlantic, NW, Coastal Atlantic Coastal

Soft-shell clam (Mya arenia) Atlantic, NE, Coastal Atlantic Coastal

Atlantic menhaden (Brevoortis tyrannus) Atlantic, NW, Coastal Atlantic Coastal

Mysid shrimp (Rhopalophthalmus Indian, Coastal Indian Coastal terranatalis) ADULTS

Mysid shrimp (Rhopalophthalmus Indian, Coastal Indian Coastal terranatalis) JUVENILES

Mysid shrimp (Mesopodopsis Indian, Coastal Indian Coastal wooldridgei) ADULTS

Mysid shrimp (Mesopodopsis Indian, Coastal Indian Coastal wooldridgei) JUVENILES Burrowing shrimp (Upogebia deltaura) Atlantic, Mediterranean, Coastal Atlantic Coastal

Antarctic Krill (Euphausia superba) Southern, Coastal Southern Coastal

Porcelain Crab (Porcellana longicornis) Atlantic, NE, Coastal Atlantic Coastal

Ocean Quahog (Arctica islandica) Atlantic, NW, Coastal Atlantic Coastal

Wrinkled Rockborer (Hiatella arctica) Atlantic, South, Coastal Atlantic Coastal

47 Bay Scallop (Pecten irradians) Atlantic, NW, Coastal Atlantic Coastal

Orange Sea Pen (Ptilosarcus gurneyi) Pacific, East, Coastal Pacific Coastal

Feather star (Oligometra serripinna) Pacific, SE Asia, Coastal Pacific Coastal Manila Clam (Tapes philippinarum) Indian, Coastal Indian Coastal

Yesso scallop (Patinopecten yessoensis) Pacific, West, Coastal Pacific Coastal

Spaghetti Bryozoan (Zoobotryon Atlantic, Caribbean, Coastal Atlantic Coastal verticillatum) Bryozoan (Electra pylosa) Pacific, West, Coastal Pacific Coastal

Bryozoan (Conopeum reticulum) Atlantic, NE, Coastal Atlantic Coastal

Bryozoan (Celleporella hyalina) Atlantic, NW, Coastal Atlantic Coastal

Bryozoan (Celleporella hyalina) Pacific, East, Coastal Pacific Coastal

Sea vase (Ciona intestinalis) Atlantic, Open Ocean Atlantic Open Ocean Sea squirt (Ascidella aspersa) Atlantic, NE, Coastal Atlantic Coastal

Polychaete worm (Myxicola Atlantic, Open Ocean Atlantic Open Ocean infundibulum) Polychaete worm (Myxicola Pacific Open Ocean Pacific Open Ocean infundibulum) Peacock worm (Sabella pavonina) Atlantic, Mediterranean, Coastal Atlantic Coastal

Keel worm (Pomatoceros triqueter) Atlantic, NE, Coastal Atlantic Coastal

Keel worm (Pomatoceros triqueter) Arctic, Coastal Arctic Coastal

Polychaete worm (Hydroides norvegica) Atlantic, Mediterranean, Coastal Atlantic Coastal

Sinistral spiral tubeworm (Spirorbis Atlantic, NE, Coastal Atlantic Coastal borealis) Sinistral spiral tubeworm (Spirorbis Pacific, East, Coastal Pacific Coastal borealis)

Polychaete worm (Salmacina dysteri) Pacific Open Ocean Pacific Open Ocean

Breadcrumb sponge (Halichondria Atlantic, NE, Coastal Atlantic Coastal panicea) Breadcrumb sponge (Halichondria Pacific, West, Coastal Pacific Coastal panicea) Common Bream (Abramis brama) Freshwater, Europe Freshwater Coastal

48 White Bream (Blicca bjoerkna) Freshwater, Europe Freshwater Coastal

Roach (Rutilus rutilus) Freshwater, Europe Freshwater Coastal

Gizzard shad (Dorosoma cepedianum) Freshwater, North America Freshwater Coastal

North Atlantic Right Whale (Eubalaena Atlantic, Open Ocean Atlantic Open Ocean glacialis)

No significant difference in microplastic consumption rates was found among species feeding in the different ocean basins (p = 0.1512, χ2 = 0.5, df = 1, Kruskal-Wallis test). Although no significant difference occurs, species feeding in the Pacific Ocean had the highest mean microplastic consumption rates compared to other ocean basins (Figure 14). Those feeding in freshwater, the Indian Ocean, and Arctic Ocean had the lowest mean microplastic consumption rates (Figure 22). When considering how bodyweight might affect these results, the analysis showed that still no significant differences occurred in NMCR values at the different ocean basins (p

= 0.09, F5,62=2.01, Kruskal-Wallis Test). Thus, it can be concluded that the ocean basin is not associated with higher risks of microplastic contamination in different filter feeders (Figure 23). Additionally, both ocean basin and filtration technique were further tested to determine if any interactions between these two variables significantly affected MCR or NMCR values. Analysis showed that when bodyweight was not considered, no significant interactions occurred between ocean basin and filtration technique to affect

MCR values (p = 0.1, F10,55=4.0, Kruskal-Wallis Test). When bodyweight was taken into consideration, the analysis determined that significant interactions still did not occur between the two variables to affect MCR values (p = 0.967, F10,55=4.0, Kruskal-Wallis Test).

49 MCR & Ocean Basin 2.00E+04

) 1.50E+04

SE ± 9.48E+03 1.00E+04

Mean Mean MCR 6.57E+03

(particles/s (particles/s 5.00E+03

4.69E+02 1.79E-04 1.68E-02 1.76E-04 0.00E+00 Arctic Atlantic Freshwater Indian Pacific Southern Ocean Basin

Figure 22. Calculated MCR value (particles/s ± SE) for the six different ocean basins/water sources: Arctic, Atlantic, Freshwater, Indian, Pacific, Southern.

NMCR & Ocean Basin 4.00E+03

) 3.00E+03

SE ±

2.00E+03

1.54E+03 Mean Mean NMCR

(particles/s/kg (particles/s/kg 1.00E+03

9.00E+01 4.48E+01 6.63E-03 5.23E+01 0.00E+00 1.72E+00 Arctic Atlantic Freshwater Indian Pacific Southern Ocean Basin

Figure 23. Calculated NMCR value (particles/s/kg ± SE) for the six different ocean basins/water sources: Arctic, Atlantic, Freshwater, Indian, Pacific, Southern.

50 f. Environment

Each species was determined to feed in one of two types of environment: coastal or open ocean (Table 8). Analysis showed that the different environments do not have significant differences in regards to MCR for study species (p = 0.173, t = -1.48, Welch two sample t-test). Despite the lack of a significant difference, species in the open ocean had a higher mean MCR compared to those in the coastal areas (Figure 24), which supports similar values found in previous studies (Barrows et al. 2018). Additional analysis considered bodyweight and showed that no significant differences in NMCR values occurred at the different environments (p = 0.173, t = -1.48, two-tailed two sample t-test). This result suggests that filter feeders in either environment experience equal risks of microplastic contamination (Figure 25).

MCR & Environment 8.00E+04

) 6.00E+04 SE

± 4.23E+04

4.00E+04 Mean Mean MCR

(particles/s 2.00E+04

7.28E+02 0.00E+00 Coastal Open Ocean Environment

Figure 24. Calculated MCR (particles/s ± SE) for each environment, coastal or open ocean.

51 NMCR & Environment 2.00E+03 1.80E+03

) 1.60E+03

SE 1.40E+03 ± 1.20E+03 1.00E+03 9.00E+02 8.00E+02 Mean Mean NMCR 6.00E+02

(particles/s/kg (particles/s/kg 4.00E+02 2.00E+02 0.00E+00 8.16E+01 Coastal Open Ocean Environment

Figure 25. Calculated NMCR (particles/s/kg ± SE) for each environment, coastal or open ocean.

V. Summary and Conclusions

Overall Risk Assessment

This study estimated the quantity of microplastics likely consumed by filter feeders and analyzed the factors that affected that statistic. When bodyweight was not taken into account, it was found that fin whales (Balaenoptera physalus) consume the highest mean quantity of microplastics per second of feeding. Given that this species can consume 10 kilograms of krill in 70,000 liters of water, this conclusion is well supported in the literature (Goldbogen et al. 2010). Relying on lunge techniques, feeding among whales (Balaenopteridae) is energetically costly (Goldbogen et al. 2008, Goldbogen et al. 2011) and inadvertently consuming microplastic particulates could potentially take a major toll on even these massive organisms. The larger species are also at risk of consuming other types of debris, including macroplastics, which could potentially block the digestive system if consumed. Globally, fin whale populations are on the rise. They are no longer considered

52 endangered, but are still labeled as vulnerable on the IUCN Red List (2019). Although this provides a better outlook than their estimated MCR values might suggest, caution must be taken to ensure that these organisms are exposed to plastics as minimally as possible. Pacific oyster larvae (Crassostrea gigas) and bryozoans (Electra pylosa), on the other hand, consumed the lowest mean quantity of microplastics per second of feeding when bodyweight was not considered. Unfortunately, neither of these species is evaluated by the IUCN (2019) and it is difficult to infer how microplastics might affect their overall population. Yet, their comparatively small MCR indicates that they likely experience lesser risk of microplastic consumption compared with most other filter feeding species, including fin whales. Similar to other species, feeding in E. pilosa and other bryozoans is expected to incur some energetic costs, as the organisms actively filter with the use of a mechanical laterofrontalfilter (Riisgard & Manriquez 1997). It is possible that some inorganic particles may be filtered out post-capture by these species, but more research is required to determine if they are actually capable of removing any sediment or debris as has been previously described (Riisgard & Manriquez, 1997). The factors found to have significant differences in mean MCR values were IUCN Red List status and filtration technique. The species with higher levels of vulnerability according to the IUCN Red List statuses (i.e. vulnerable and endangered species) had higher mean Microplastic Consumption Rates compared to those that were not evaluated or threatened. Species that had not yet been evaluated tended to be small and widely distributed, including crabs, scallops, bryozoans, sea worms, tunicates, and copepods (IUCN 2019). Such organisms are generally incapable of filtering massive quantities of particulates from the water regardless of the microplastic abundance in their location. This result can be beneficial to resource managers IUCN could potentially be used as a predictor, as it shows that vulnerable species are more likely to consume higher quantities of microplastics over time. Effective strategies, then, could be implemented to protect these species. It is important to note, however, that reasonable biological characteristics must be met to use this factor as a predictor for specific species. The variable has only been considered in terms of marine filter feeding species and thus, conclusions should only be drawn for similar organisms. Species that filter water with lunge feeding techniques, such as humpback whales and bowhead whales, had significantly higher mean MCR values compared with those that

53 rely on other techniques, like suspension feeding or water pumping. Lunge feeding is energetically expensive and, as such, it is a method frequently used by larger and stronger species, which are also capable of filtering greater quantities of particulates in water (Acevedo-Gutiérrez et al. 2002, Watkins & Schevill 1979). The remaining variables considered – ocean basin, environment, life stage, and salinity – were not found to have significant differences in Microplastic Consumption Rates. These factors, then, do not increase or decrease the risk that individuals will experience higher risks of microplastic consumption. Although fin whales – the species with highest mean MCR in this study – are known to feed in offshore, subpolar marine waters (Vikingsson et al. 2009), for example, it is impossible to conclude from this knowledge that they are at risk of consuming high quantities of microplastics. Instead, it is much more valuable to consider the population’s vulnerability and filtration technique. Similarly, Pacific oyster larvae and bryozoans are known to feed in coastal marine waters of the Atlantic and Pacific Oceans (Harris 2008, Fey et al. 2010, Cognie et al. 2006). Though these areas tend to have a lower abundance of microplastics, conclusions cannot be drawn without first considering vulnerability and filtration techniques. Both of these species are not yet evaluated by the IUCN (2019) and rely on suspension feeding techniques (Gerdes 1982, Harris 2008, Riisgard & Manriquez 1997), factors that support the conclusion that such species are not at great risks of microplastic consumption. When bodyweight is factored into the analysis, results showed that pelagic tunicates (P. confederata) had the highest NMCR values. As one of the smallest species studied in this review, this result is likely caused by the species’ incredible efficiency and high filtration rate in relation to its size. No other factors considered here would have had a significant effect on the NMCR, so it would be important for future studies to take this into account. The only factor that had a significant relationship with NMCR was salinity, while the remaining variables did not experience significant differences. Bryozoans (E. pylosa) still experienced the smallest NMCR values. Thus, it can be concluded that pelagic tunicates experience the highest risk of microplastic contamination, while bryozoans experience the lowest risk of contamination. Understanding the species most at risk of consuming microplastics – including fin whales (Balaenoptera physalus), North Atlantic right whales (Eubalaena glacialis), and

54 bowhead whales (Balaena mysticetus) – is critical because these particles are known to contain toxic chemicals and pose serious dangers to the species that consume them (Gallo et al. 2018). Chemicals commonly associated with microplastics include Persistent Organic Pollutants, polychlorinated biphenyls, and Persistent, Bioaccumulative, and Toxic Compounds, are found in marine plastic litter (Gallo et al. 2018, Lusher et al. 2017). Some of these chemicals and additives have endocrine disrupting properties (Lusher et al. 2017). And PBTs are known to bioaccumulate, leading to the dangerous hazards that plastics pose (Lusher et al. 2017). Toxins and chemicals frequently associated with microplastics are often either added during the manufacturing process or absorbed from the surrounding environment. These harmful additives are expected to have significant and detrimental effects on entire populations and ecosystems, as they can reduce an individual’s ability to survive in their environment (Gallo et al. 2018). The whale species found to be most at risk of consuming microplastic are thus more likely to be exposed to such toxins and chemicals, providing them with yet another human-caused challenge to overcome and recover from their statuses as endangered or vulnerable species. This study also considered factors that affect microplastic abundance. It was determined there are significant differences in microplastic abundance among the ocean basins and between the different environments. The open ocean had higher mean microplastic abundance in surface waters compared to coastal environments. Furthermore, the Arctic and Southern Oceans had significantly higher mean levels of microplastic abundance than other basins. This can pose a potentially substantial problem in the Arctic Ocean, because researchers expect that climate change may lead to the release of even greater quantities of microplastics from melting sea ice in the region (Lusher et al. 2017). When drawing conclusions from these results, however, caution must be taken because data was not equally distributed between the different oceans. Far fewer water samples existed in the Arctic and Southern Oceans than in the Pacific, Atlantic, and Indian Oceans, and this disparity could cause the results to be slightly unreliable. The presence of marine litter has been a problem for decades in the open ocean, as solid waste was frequently discarded from ships prior to the 1980s (Lusher et al. 2017), most likely due to ghost gear or shipping container losses. Yet, even as international regulations and conservation efforts attempt to reduce the quantity of microplastics in

55 offshore waters today, the findings in this study show that open ocean environments continue to harbor vast quantities of litter. Due to continuous ocean currents, improper waste disposal, and dramatic events, such as floods and cyclones, it can be extremely difficult to manage the levels of marine litter found ((Lusher et al. 2017).

Future Considerations

In future studies, it would be beneficial for researchers to focus on individual species and consider their specific and unique risks in terms of microplastic consumption. Here, it was necessary to make generalizations and estimates of geographic distribution for each species simply as a result of the quantity of species considered throughout the review. Although the mean filtration rates would remain the same for each species, geographic distribution greatly determines the quantity of microplastics to which filter feeders are exposed. It was extremely beneficial to take an overall assessment of the many different filter feeders to better understand which are most at-risk of consuming toxic particulates and which factors affect that risk. But focusing future studies on specific species – particularly those that are commercially and ecologically important – could further this understanding. Additionally, consumption of macroplastics is an important topic to highlight in future studies. Communities around the globe are familiar with widely publicized news articles concerning the occurrence of beached (Lusher et al. 2015). Many of the necropsies that result from these incidents indicate that macroplastics are consumed, particularly in whales, sharks, seabirds, and other species that are vital to the ecosystem (Lusher et al. 2015, Bråte et al. 2017). It would be beneficial to develop a broader understanding of the risks associated with macroplastic consumption in conjunction with the risks of microplastic consumption, as reviewed in this paper. Such an understanding could illuminate the different ecological impacts associated with plastics of varying sizes. Previous studies have also shown that mesh size and the size of microplastic particulates should be considered when evaluating microplastic consumption (Roesch et al. 2013, Zhao et al. 2014). Thus, it would be beneficial if future studies consider how specific mesh sizes of gill rakers in each species, as well as the average microplastic size, could potentially affect the quantity of microplastics consumed. This paper aimed to

56 determine how likely it is that different filter feeding species will consume microplastics, and while this complex problem was simplified to estimate risks of consumption and contamination for many different species for the purposes of this review, it did not provide concrete quantities of microplastics actually consumed. With the use of ever emerging technologies and techniques, it is expected that actual consumption data will be provided for many of these species, allowing researchers to consider these risks further and more accurately predict their ecologically and environmental impacts.

57 VI. References

Acevedo-Gutierrez A, Croll DA, Tershy BR. 2002. High feeding costs limit dive time in the largest whales. J Exp Biol. 205: 1747-1753.

Achituv Y, Pedrotti ML. 1999. Costs and gains of porcelain crab suspension feeding in different flow conditions. Mar Ecol Prog Ser. 184: 161-169. doi: 10.3354/meps184161.

Aksnes DL, Magnesen T. 1979. Distribution, development, and production of Calanus finmarcchicus (Gunnerus) of Lindåspollene, western Norway, 1979. Sarsia. 68: 195- 208. doi: 10.1080/00364827.1983.10420572.

Ali RM. 1970. The influence of suspension density and temperature on the filtration rate of Hiatella arctica. Mar Biol. 6(4): 391-302. doi: 10.1007/BF00353662.

Amat JN, Tempera F. 2009. Zoobotryon verticillatum Della Chiaje, 1822 (), a new occurrence in the archipelago of the Azores (North-Eastern Atlantic). Mar Poll Bul. 58: 761-764. doi: 10.1016/j.marpolbul.2009.02.019.

Armonies W, Reise K. 2003. Empty habitat in coastal sediments for populations of macrozoobenthos. Helgol Mar Res. 56:279-287. doi: 10.1007/s10152-002-0129-8.

Ashjian CJ, Braund SR, Campbell RG, George JC, Kruse J, Maslowski W, Moore SE, Nicolson CR, Okkonnen SR, Sherr BF, Sherr EB, Spitz YH. 2010. Climate, variability, oceanography, distribution, and Iñupiat subsistence whaling near Barrow, Alaska. Arctic. 63(2): 179-194.

Atkinson A, Siegel V, Pakhomov EA, Rothery P, Loeb V, Ross RM, Quetin LB, Schmidt K, Fretwell P, Murphy EJ, Tarlin GA, Fleming AH. 2008. Oceanic circumpolar habitats of Antarctic krill. Mar Ecol Prog Ser. 362:1-23. doi: 10.3354/meps07498.

Austin WC, Conway KW, Barrie JV, Krautter M. 2007. Growth and morphology of a reef- forming glass sponge, Aphrocallistes vastus (Hexactinellida), and implications for recovery from widespread trawl damage. Porifera Research: Biodiversity, Innovation & Sustainability, Serie Livros. 28: 139-145.

Barrows APW, Cathey SE, Petersen CW. 2018. Marine environment microfiber contamination: Global patterns and the diversity of microparticle origins. Environ Pollut 237:275-284. doi:10.1016/j.envpol.2018.02.062.

Baumgartner MF, Cole TVN, Clapham PJ, Mate BR. 2003. North Atlantic right whale habitat in the lower Bay of Fundy and on the SW Scotian Shelf during 1999–2001. Mar Ecol Prog Ser. 264:137-154. doi: 10.3354/meps264137.

58 Baumgartner MF, Mate BR. 2003. Summertime foraging ecology of North Atlantic right whales. Mar Ecol Prog Ser. 264:123-135. doi: 10.3354/meps264123.

Baumgartner MF, Mate BR. 2004. Summer and fall habitat of North Atlantic right whales (Eubalaena glacialis) inferred from satellite telemetry. Can J Fish Aquat Sci. 62:527- 543. doi: 10.1139/F04-238.

Best, BA. 1988. Passive suspension feeding in a sea pen: Effects of ambient flow on volume flow rate and filtering efficiency. Biol Bul. 175(3): 332-342. doi: 10.2307/1541723.

Bochdansky AB, Deibel D. 1998. Functional feeding response and behavioral ecology of Oikopleura vanhoeffeni (Appendicularia, Tunicata). J Exp Mar Biol Ecol. 233: 181- 211.

Boström C, Bonsdorff E. 1997. Community structure and spatial variation in benthic associated with Zostra marina (L.) beds in the northern . J of Sea Res. 37: 153-166. doi: 10.1016/S1385-1101(96)00007-X.

Boyd CM, Heyraud M, Boyd C. 1984. Feeding of the antarctic krill Euphausia superba. J of Crustacean Biol, 4: 123-141.

Bråte ILN, Huwer B, Thomas KV, Eidsvoll DP, Halsband C, Almroth BC, Lusher A. 2017. Micro-and-macro-plastics in marine species from Nordic waters. Nordic Council of Ministers. 549. doi: 10.6027/TN2017-549.

Braun CD, Skomal GB, Thorrold SR, Berumen ML. Diving behavior of the reef manta ray links coral reefs with adjacent deep pelagic habitats. PLoS ONE 9(2): e88170. doi: 10.1371/journal.pone.0088170.

Buchheister A, Miller TJ, Houde ED, Secor DH, Latour RJ. 2016. Spatial and temporal dynamics of Atlantic menhaden (Brevoortia tyrannus) recruitment in the Northwest Atlantic Ocean. ICES J Mar Sci. 73(4): 1147-1159. doi: 10.1093/icesjms/fsv260.

Buhl-Mortensen L, Vanreusel A, Gooday AJ, Levin LA, Priede IG, Buhl-Mortensen P, Gheerardyn H, King NJ, Raes M. 2009. Biological structures as a source of habitat heterogeneity and biodiversity on the deep ocean margins. Mar Ecol. 31: 21-50. doi: 10.1111/j.1439-0485.2010.00359.x.

Bullivant JS. 1968. The rate of feeding of the bryozoan, Zoobotryon verticillatum. New Zealand J Mar Freshwater Res, 2:111-134, doi: 10.1080/00288330.1968.9515230.

Cargnelli LM, Griesbach SJ, Packer DB, Weissberger E. 1999. Essential Fish Habitat Source Document: Ocean Quahog, Arctica islandica, Life History and Habitat Characteristics. NOAA Technical Memorandum NMFS-NE-148.

59 Chebbi N, Mastrototaro F, Missaoui H. Spatial distribution of ascidians in two Tunisian lagoons of the Mediterranean Sea. Cah Biol Mar. 51: 117-127.

Chipman WA, Hawkins JG. 1954. Water filtration by the bay scallop, pecten irradians, as observed with the use of radioactive . Biol Bul. 107: 80-91. doi: 10.2307/1538632.

Christiansen ME. 2000. On the occurrence of thalassinidea (decapoda) in norwegian waters. J Crustacean Biol. 20(2): 230-237. doi: 10.1163/1937240X-90000024.

Clarke A, Tyler PA. 2008. Adult Antarctic krill feeding at abyssal depths. Current Biol. 18: 282-285. doi: 10.1016/j.cub.2008.01.059.

Cognie B, Haure J, Barillé. 2006. Spatial distribution in a temperate coastal ecosystem of the wild stock of the farmed oyster Crassostrea gigas (Thunberg). Aquacult. 259:249-259. doi: 10.1016/j.aquaculture.2006.05.037

Cole M Lindeque P, Fileman E, Halsband C, Goodhead R, Moger J, Galloway TS. 2013. Microplastic ingestion by zooplankton. Environ Sci & Technol 47(12): 6646-6655. doi:10.1021/es400663f

Compagno LJV. 2001. Sharks of the world: An annotated and illustrated catalogue of shark species known to date. FAO Spec Cat Fish Purp 1(2): 269. FAO, Rome.

Dales RP. 1957. Some quantitative aspects of feeding in sabellid and serpulid fan worms. J Mar Biol Assoc UK. 36: 309-316.

Dang C, de Montaudouin X, Gam M, Paroissin C, Bru N, Caill-Milly N. 2010. The Manila clam population in Arcachon Bay (SW France): Can it be kept sustainable? J Sea Res. 63:108-118. 10.1016/j.seares.2009.11.003. de la Parra Venegas R, Hueter R, Cano JG, Tyminski J, Remolina JG, Maslanka M, Ormos A, Weigt L, Carlson B, Dove A. 2011. An unprecedented aggregation of whale sharks, Rhincodon typus, in Mexican coastal waters of the Caribbean Sea. PLoS ONE 6(4): e18994. doi:10.1371/journal.pone.0018994.

Dewar H, Mous P, Domeier M, Muljadi A, Pet J, Whitty J. 2008. Movements and site Wdelity of the giant manta ray, Manta birostris, in the Komodo Marine Park, Indonesia. Mar Biol. 155:121-133. doi: 10.1007/s00227-008-0988-x.

Divi RV, Strother JA, Paig-Tran EWM. 2018. Manta rays feed using ricochet separation, a novel nonclogging filtration mechanism. Sci Adv. 4(9). doi: 10.1126/sciadv.aat9533.

Doniol-Valcroze T, Lesage V, Giard J, Michaud R. (2011). Optimal foraging theory predicts diving and feeding strategies of the largest marine predator. Behav Ecol. 22:880-888. 10.1093/beheco/arr038.

60

Duffy CAJ. 2002. Distribution, seasonality, lengths, and feeding behaviour of whale sharks (Rhincodon typus) observed in New Zealand waters. New Zealand Journal of Marine and Freshwater Research 36:565-570. doi: 0028-8330/02/3603-0565.

Durbin AG, Durbin EG. 1975. rates of the Atlantic menhaden Brevoortia tyrannus as a function of particle size and concentration. Mar Biol. 33: 265-277. doi: 10.1007/BF00390931.

D’Vincent CG, Nilson RM, Hanna RE. 1985. Vocalization and coordinated feeding behavior of the humpback whale in southeastern Alaska. Sci Rep Whales Res Inst. 36: 41-47.

Drenner RW, Mummert JR, deNoyelles F, Kettle D. 1984. Selective particle ingestion by a filter-feeding fish and its impact on phytoplankton community structure. Limnol Oceanogr. 29(5): 941-948. doi: 10.4319/lo.1984.29.5.0941.

Ekaratne K, Burfitt AH, Flowerdew MW, Crisp DJ. 1982. Separation of the two Atlantic species: Pomatoceros, P. lamarckii, and P. triqueter (Annelida: Serpulidae) by means of biochemical genetics. Mar Biol. 1:257-264. doi: 10.1007/BF00397042.

Eldrige LC, Smith CM. 2001. A guidebook of introduced marine species in Hawaii. Bishop Museum Technical Report. 21.

Fey F, Dankers N, Steenbergen J, Goudswaard K. 2010. Development and distribution of the non-indigenous Pacific oyster (Crassostrea gigas) in the Dutch Wadden Sea. Aquacult Int. 18:45-59. doi: 10.1007/s10499-009-9268-0

Fiedler PC, Reilly SB, Hewitt RP, Demer D, Philbrick VA, Smith S, Armstrong W, Croll DA, Tershy BR, Mate BR. 1998. habitate and prey in the California Channel Islands. Deep Sea Res. 2(45): 1781-1801. doi: 10.1016/S0967- 0645(98)80017-9.

Forester AJ. 1979. The association between the sponge Halichondria panacea (pallas) and scallop Chlamys varza (l.): a commensal-protective mutualism. J Exp Mar Biol Ecol. 36:1-10. doi: 10.1016/0022-0981(79)90096-0.

Friedlaender AS, Goldbogen JA, Nowacek DP, Read AJ, Johnston D, Gales N. 2014. Feeding rates and under-ice foraging strategies of the smallest lunge filter feeder, the Antarctic minke whale (Balaenoptera bonaerensis). J Exp Biol. 217:2851-2854. doi: 10.1242/jeb.106682.

Froneman PW. 2001. Feeding ecology of the mysid, Mesopodopsis wooldridgei, in a temperate estuary along the eastern seaboard of South Africa. J Plankton Res. 23(9): 999-1008. doi: 10.1093/plankt/23.9.999.

61 Fuller JL, Clarke GL. Further experiments on the feeding of Calanus finmarchicus. Biol Bul. 70(2): 308-320. 10.2307/1537475.

Gallo F, Fossi C, Weber R, Santillo D, Sousa J, Ingram I, Nadal A, Romano D. 2018. Marine litter plastics and microplastics and their toxic chemicals components: The need for urgent preventive measures. Environ Sci Eur 30(1): 13. doi:10.1186/s12302-018- 0139-z

Gerdes D. 1983. The Pacific oyster Crassostrea gigas part I. Feeding behaviour of larvae and adults. Aquaculture. 31:195-219. doi: 10.1016/0044-8486(83)90313-7.

Gill P. (2002). A blue whale (Balaenoptera musculus) feeding ground in a southern Australian coastal zone. 4:179-184.

Gill PC, Morrice MG, Page B, Pirzl R, Levings AH, Coyne M. 2011. Blue whale habitat selection and within-season distribution in a regional upwelling system off southern Australia. Mar Ecol Prog Ser. 421: 243-263. doi: 10.3354/meps08914

Goldbogen JA. 2010. The ultimate mouthful: lunge feeding in rorqual whales. Am Sci. 98(2): 124-131. doi: 10.1511/2010.83.124.

Goldbogen JA, Cade DE, Calambokidis J, Friedlaender AS, Potvin J, Segre PS, Werth AJ. 2017. How baleen whales feed: the biomechanics of engulfment and filtration. Annu Rev Mar Sci 9:367-86. doi: 10.1146/annurev-marine-122414-033905

Goldbogen JA, Calambokidis J, Croll DA, Harvey JT, Newton KM, Oleson EM, Schorr G, Shadwick RE. 2008. Foraging behavior of humpback whales: kinematic and respiratory patterns suggest a high cost for a lunge. J Exp Biol. 211: 3712-3719. doi:10.1242/jeb.023366

Goldbogen JA, Calambokidis J, Oleson E, Potvin J, Pyenson ND, Schorr G, Shadwick RE. 2011. Mechanics, hydrodynamics and energetics of blue whale lunge feeding: efficiency dependence on krill density. J Exp Biol. 214: 131-146. doi: 10.1242/jeb.054726.

Gordillo S. 2001. Puzzling distribution of the fossil and living genus Hiatella (). Palaeogeography, Palaeoclimatology, Palaeoecology. 165: 231-249. doi: 10.1016/S0031-0182(00)00162-0.

Gorsky G, Fisher NS, Fowler SW. 1983. Biogenic debris from the pelagic tunicate, Oikoplewu dioicu, and its role in the vertical transport of a transuranium element. Estuar Coast Shelf Sci. 18:13-23. doi: 10.1016/0272-7714(84)90003-9.

Gotshall D. 2005. Guide to marine invertebrates —Alaska to Baja (2d ed.): Monterey." Calif., Sea Challengers.

62 Graham RT, Roberts CM, Smart JCR. 2005. Diving behaviour of whale sharks in relation to a predictable food pulse. J. R. Soc. Interface. 3:109-116. doi:10.1098/rsif.2005.0082.

Greathead C, Demain D, Dobby H, Allan L, Weetman A. 2011. Quantitative Assessment of the Distribution and Abundance of the Burrowing and Large Epifauna Community in the Fladen Fishing Ground, Northern North Sea. Scottish Mar Fresh Sci. 2(2).

Hain JHW, Carter GR, Kraus SD, Mayo CA, Winn HE. 1982. Feeding behavior of the humpback whale, Megaptera novaeangliae, in the western North Atlantic. Fish Bul. 80(2): 259-268.

Hansen IV, Weeks JM, Depledge MH. 1995. Accumulation of copper, zinc, cadmium and chromium by the marine sponge Halichondria panicea Pallas and the implications for biomonitoring. Mar Poll Bul. 31: 133-138.

Harbison GR, Campenot RB. 1979. Effects of temperature on the swimming of (Tunicata, ): Implications for vertical migration. Limnol Oceangr. 24(6): 1081-1091. doi: 10.4319/lo.1979.24.6.1081.

Harbison GR, Gilmer RW. 1976. Feeding rates of the pelagic tunicate Pegea confederata and two other salps. Limnol Oceanogr. 21(4): 517-528. doi: 10.4319/lo.1976.21.4.0517

Harris J. 2008. Pacific oyster, Crassostrea gigas (Thunberg, 1793). Aquat Inv Ecol.

Havenhand JN, Svane I. 1991. Roles of hydrodynamics and larval behaviour in determining spatial aggregation in the tunicate Ciona intestinalis. Mar Ecol Prog Ser. 68: 271- 276. doi: 10.3354/meps068271.

Hellal AM. 2012. Taxonomic Study on the Feather Stars (Crinoidea: Echinodermata) from Egyptian Red Sea and Suez Canal, Egypt. Open J Mar Sci. 2:51-57.

Herdis L, Nøttestad L, Skaret G, Broms C, Fernö A. 2012. Overlap in distribution and diets of Atlantic mackerel (Scomber scombrus), Norwegian spring-spawning ( harengus) and blue whiting (Micromesistius poutassou) in the Norwegian Sea during late summer. Mar Biol Res. 8:5-6, 442- 460, doi: 10.1080/17451000.2011.642803.

Hermansen P, Larsen PS, Riisgard HU. 2001. Colony growth rate of encrusting marine bryozoans (Electra pilosa and Celleporella hyalina). J Exp Mar Biol Ecol. 263:1-23. doi: 10.1016/S0022-0981(01)00243-X.

Heyman WD, Graham RT, Kjerfve B, Johannes RE. 2001. Whale sharks Rhincondon typus aggregate to feed on fish in Belize. Mar Eco Prog Ser 215: 275-282. doi: 10.3354/meps215275.

63 Hill SL, Phillips T, Atkinson A. 2013. Potential Climate Change Effects on the Habitat of Antarctic Krill in the Weddell Quadrant of the Southern Ocean. PLoS ONE 8(8): e72246. doi:10.1371/journal.pone.0072246.

Holland, ND, Leonard AB, Meyer DL. 1991. Digestive mechanics and gluttinous feeding in the feather star Oligometra serripinna (Echinodermata: Crinoidea). Mar Biol. 111:113-119. 10.1007/BF01986352.

Hopcroft RR, Roff JC. 1995. Zooplankton growth rates: extraordinary production by the larvacean Oikopleura dioica in tropical waters. J Plankton Res. 17(2): 205-220. doi: 10.1093/plankt/17.2.205.

Howe RL, Rees AP, Widdicombe S. 2004. The impact of two species of bioturbating shrimp (Callianassa subterranea and Upogebia deltaura) on sediment denitrification. J Mar Biol Assoc UK. 84: 629-632. doi: 10.1017/S002531540400966Xh.

Isaac MJ. 1974. Copepoda monstrilloida from South-west Britain including six new species. J Mar Biol Assoc UK. 54:126-140. doi: 10.1017/S0025315400022116.

IUCN 2019. The IUCN Red List of Threatened Species. Version 2019-2. http://www.iucnredlist.org. Downloaded on 18 July 2019.

Jebakumar JPP, Nandhagopel G, RajanBabu B, Ragumaran S, Ravichandran V, Marchini A, Minchin D. 2017. Bioinvasion Rec. 6. doi: 10.3391/bir.2017.6.3.05.

Jerling HL, Wooldridge TH. Comparative morphology of the feeding appendages of four mesozooplankton species in the Sundays River estuary. African Zool. 29(4).

Kahn AS, Yahel G, Chu JWF, Tunnicliffe V, Leys SP. 2015. Benthic grazing and carbon sequestration by deep-water glass sponge reefs. Limnol Oceanogr. 60: 78-88. doi: 10.1002/lno.10002.

Kanhai LDK, Officer R, Lyashevska O, Thompson RC, O'Connor I. 2017. Microplastic abundance, distribution and composition along a latitudinal gradient in the Atlantic Ocean. Mar Pollut Bull 115(1): 307-314. doi: 10.1016/j.marpolbul.2016.12.025.

Kater BJ, Geurts van Kessel AJM, Baars JJMD. 2006. Distribution of cockles Cerastoderma edule in the Eastern Scheldt: habitat mapping with abiotic variables. Mar Ecol Prog Ser. 318: 221-227.

Kotta J, Orav-Kotta H. 2001. Role of benthic macroalgae in regulating macrozoobenthic assemblages in the Väinameri (north-eastern Baltic Sea). Ann. Zool. Fennici. 38: 163- 171.

64 Kuparinen A, Vinni M, Teacher AGF, Kähkönen K, Merliä J. 2013. Mechanism of hybridization between bream Abramis brama and roach Rutilus rutilus in their native range. J Fish Biol. 84:237-242. doi: 10.1111/jfb.12272.

Kupriyanova EK, Badyaev AV. 1998. Ecological correlates of Arctic serpulidae (Annelida, Polychaeta) distributions. Ophelia. 49(3), 181-193. doi: 10.1080/00785326.1998.10409381.

Laidre KL, Heide-Jørgensen MP, Nielson TG. Role of the bowhead whale as a predator in West Greenland. Mar Ecol Prog Ser 346: 285-297. doi: 10.3354/meps06995.

Lammens E. 1986. Interactions between and the structure of fish communities in Dutch shallow eutrophic lakes. Landbouwuniversitet Waginingen

Lance J. 1964. The Salinity Tolerances of Some Estuarine Planktonic Crustaceans . Biol Bul. 127(1): 108-118. doi: 10.2307/1539348.

Lazarro X. 1987. A reviewof planktivorous fishes: Their evolution, feeding behaviours, selectivities, and impacts. Hydrobiologia. 146:97-167. doi: 10.1007/BF00008764.

Leonard AB, Strickler JR, Holland ND. Effects of current speed on filtration during suspension feeding in Oligometra serripinna (Echinodermata: Crinoidea). Mar Biol. 97: 111-125. doi: 10.1007/BF00391251.

Leslie HA. 2014. Review of microplastics in cosmetics: Scientific background on a potential source of plastic particulate marine litter to support decision-making. IVM Institute for Environmental Studies.

Lewis TL, Esler D, Boyd WS. Effects of by sea ducks on clam abundance in soft- bottom intertidal habitats. Mar Ecol Prog Ser. 329: 131-144. doi: 10.3354/meps329131.

Leys SP, Yahel G, Reidenbach MA, Tunnicliffe V, Shavit U, Reiswig HM. The Sponge Pump: The Role of Current Induced Flow in the Design of the Sponge Body Plan. PLoS ONE 6(12): e27787. doi:10.1371/journal.pone.0027787.

Linnaeus, C. (1758). Systema Naturae per regna tria naturae, secundum classes, ordines, genera, species, cum characteribus, differentiis, synonymis, locis. Editio decima, reformata [10th revised edition], vol. 1: 824.

Love JW, Johnson AK, May EB. 2006. Spatial and temporal differences of Atlantic menhaden (Brevoortia tyrannus) recruitment across major drainages (1966-2004) of the watershed. Estuar Coasts. 29(5): 794-801. doi: 10.1007/BF02786530.

65 Lusher AL, Hernandez-Milian G, O’Brien J, Berrow S, O’Connor I, Officer R. 2015. Microplastic and macroplastic ingestion by a deep diving, oceanic cetacean: The True’s beaked whale Mesoplodon mirus. Env Pol. 199: 185-191. doi: 10.1016/j.envpol.2015.01.023.

Lusher AL, Hollman P, Mendoza-Hill J. 2017. Microplastics in and aquaculture: Status of knowledge on their occurence and implications for aquatic organisms and food safety. FAO and Aquaculture Technical Paper. Rome. 615. doi:978-92- 5-109882-0.

Lyons J, Lucas MC. 2002. The combined use of acoustic tracking and echosounding to investigate the movement and distribution of common bream (Abramis brama) in the River Trent, England. Hydrobiologia. 483: 265-273. doi: 10.1023/A:1021364504129.

MacKenzie CL. 2008. The Bay Scallop, Argopecten irradians, Massachusetts Through North Carolina: Its Biology and the History of Its Habitats and Fisheries. Mar Fish Rev. 70:3-4.

Marshall SM, Nicholls AG, Orr AP. 1934. On the biology of Calanus finmarchicus. V. seasons distribution, size, weight and chemical composition in Loch Striven in 1933 and their relation to the phytoplankton. J Mar Biol Assoc of the UK 19: 793–828. doi: 10.1017/S0025315400046804.

Mastrototaro F, D’onghia G, Tursi A. Spatial and seasonal distribution of ascidians in a semi- enclosed basin of the Mediterranean Sea. J Mar Biol Assoc UK. 88(5): 1053-1061. doi: 10.1017/S0025315408001392.

Minchin D. 2012. Rapid assessment of the bryozoan, Zoobotryon verticillatum (Delle Chiaje, 1822) in marinas, Canary Islands. Mar Poll Bul. 64: 2146-2150. doi: 10.1016/j.marpolbul.2012.07.041.

McCann L, Keith I, Carlton JT, Ruiz GM, Dawson TP, Collins K. 2015. First record of the non-native bryozoan Amathia (= Zoobotryon) verticillata (delle Chiaje, 1822) (Ctenostomata) in the Galápagos Islands. Bioinvasions Rec. 4(4): 255-260. doi: 10.3391/bir.2015.4.4.04 .

Mizroch SA, Rice DW, Breiwick JM. 1984. The Fin Whale, Balaenoptera physalus. Mar Fish Rev. 46(4): 20-24.

Moen TL. 2006. A translation of Bishop Gunnerus’ description of the species Hydroides norvegicus with comments on his triqvetra. Scienta Marina. 70S3: 115- 123. doi: 10.3989/scimar.2006.70s3115.

Monestiez P, Dubroca L, Bonnin E, Durbec JP, Guinet C. 2006. Geostatistical modelling of spatial distribution of Balaenoptera physalus in the Northwestern Mediterranean Sea

66 from sparse count data and heterogeneous observation efforts. Ecol Model. 193: 615- 628. doi: 10.1016/j.ecolmodel.2005.08.042.

Moore CJ, Moore SL, Leecaster MK, Weisberg SB. 2001. A comparison of plastic and plankton in the North Pacific Central Gyre. Mar Pollut Bull 42(12): 1297-1300. doi:https://doi.org/10.1016/S0025-326X(01)00114-X

Moore SE, George JC, Sheffield G, Bacon J, Ashjian CJ. Bowhead whale distribution and feeding near Barrow, Alaska, in late summer 2005-06. Arctic. 63(2): 195-205.

Motta PJ, Maslanka M, Hueter RE, Davis RL, de la Parra R, Mulvany SL, Habegger ML, Strother JA, Mara KR, Gardiner JM, Tyminski JP, Zeigler, LD. 2010. Feeding anatomy, filter-feeding rate, and diet of whale sharks Rhincodon typus during surface ram filter feeding off the Yucatan Peninsula, Mexico. Zoology 113(4): 199-212. doi: https://doi.org/10.1016/j.zool.2009.12.001.

Murray JM, Watson GJ, Giangrande A, Bentley MG, Farrell P. 2011. Reproductive biology and population ecology of the marine fan worm Sabella pavonina (Savigny) (Polychaeta: Sabellidae). Invertebr Reprod Dev. 55(3): 183-196. doi: 10.1080/07924259.2011.555619.

Nakamura Y. 2004. Suspension feeding and growth of juvenile Manila clam Ruditapes philippinarum reared in the laboratory. J Exp Mar Biol Ecol. 266(2): 181-192. doi: 10.1016/S0022-0981(01)00354-9.

Nikulina EA, Hanel R, Schäfer P. 2007. Cryptic speciation and paraphyly in the cosmopolitan bryozoan Electra pilosa—Impact of the Tethys closing on species evolution. Mol Phylogenet Evol. 45:765-776. doi: 10.1016/j.ympev.2007.07.016.

Nishi E. 1993. Notes on reproductive biology of some serpulid polychaetes at Sesoko Island, Okinawa, with brief accounts of setal morphology of three species of Salmacina and Filograna implexa. 10(1): 11-16. doi: 10.4282/sosj1979.10.11.

O’Connor RJ, Lamont P. 1978. The spatial organization of an intertidal spirorbis community. J Exp Mar Biol Ecol. 32: 143-169. doi: 10.1016/0022-0981(78)90112-0.

Ohsumi S, Masaki Y, Kawamura A. 1970. Stock of the Antarctic Minke Whale. Sci Rep Whales Res Inst. 22: 5-125.

Oleson NJ. 1995. Clearance potential of jellyfish Aurelia aurita, and predaction impact on zooplankton in a shallow cove. Mar Ecol Pro Ser. 124: 63-72.

Overholtz WJ, Hare, JA, Keith, CM. 2011. Impacts of Interannual Environmental Forcing and Climate Change on the Distribution of Atlantic Mackerel on the U.S. Northeast Continental Shelf. Mar Coast Fish. 3: 219-232. doi: 10.1080/19425120.2011.578485.

67 Paig-Tran EWM, Bizzarro JJ, Strother JA, Summers AP. 2011. Bottles as models: predicting the effects of varying swimming speed and morphology on size selectivity and filtering efficiency in fishes. J Exp Biol. 214: 1643-1654. doi: 10.1242/jeb.048702.

Paig-Tran EWM, Kleinteich T, Summers AP. 2013. The Filter Pads and Filtration Mechanisms of the Devil Rays: Variation at Macro and Microscopic Scales. J Morphol. 274(9): 1026-43. doi: 10.1002/jmor.20160.

Panigada S, Zanardelli M, Canese S, Jahoda M. 1999. How deep can baleen whales dive? Mar Ecol Prog Ser. 187:309-311.

Parnell PE. 2001. The distribution of estuarine and oceanic water masses on the southern shore of O’ahu, Hawai’i: Ecological and coastal management implications, and novel methodology. Limnol Oceanogr. 46(6): 1468-1485. doi: 10.4319/lo.2001.46.6.1468.

Paul S, Calliari D. 2016. Salinity and temperature tolerances of Neomysis americana (Crustacea: ) sub adults: Perspectives on vulnerability, and distribution in South America. J Exp Mar Biol Ecol. 486: 373-378. doi: 10.1016/j.jembe.2016.10.027.

Peattie ME, Hoare R. 1980. The sublittoral ccology of the Menai Strait, II. The sponge Halichondria panicea (Pallas) and its associated fauna. Estuar Coast Shelf Sci. 13:621-635.

Ponti M, Abbiati M, Ceccherelli VU. 2002. Drilling platforms as artificial reefs: distribution of macrobenthic assemblages of the ‘‘Paguro’’ wreck (northern Adriatic Sea). ICES J Mar Sci. 59: S316–S323. doi: 10.1006/jmsc.2002.1225.

Priede IG, Miller PI. 2009. A basking shark (Cetorhinus maximus) tracked by satellite together with simultaneous remote sensing II: New analysis reveals orientation to a thermal front. Fisheries Research. 95: 370-372. doi: 10.1016/j.fishres.2008.09.038

Prokupchuck I, Senyabov E. 2006. Diets of herring, mackerel, and blue whiting in the Norwegian Sea in relation to Calanus finmarchicus distribution and temperature conditions. J Mar Sci. 63: 117-127. doi: 10.1016/j.icesjms.2005.08.005.

Qiu T, Liu Y, Zheng J, Zhang T, Qi J. 2015. A feeding model of oyster larvae (Crassostrea angulata). Physiol Behav. 147: 169-174. doi: 10.1016/j.physbeh.2015.04.043.

Randlov A, Riisgard HU. 1979. Efficiency of Particle Retention and Filtration Rate in Four Species of Ascidians. Mar Ecol Prog Ser. 1:55-59. doi: 10.3354/meps001055.

Richardson CA, Ibarrola I, Ingham RJ. 1993. Emergence pattern and spatial distribution of the Cerastoderma edule. Mar Ecol Prog Ser. 99: 71-81. doi: 10.3354/meps099071.

68 Riisgard HU. 1991. Filtration rate and growth in the blue mussel, Mytilus edulis, Linnaeaus 1758: Dependence on algal concentration. J Shellfish Res. 10(1): 29-35.

Riisgard HU, Kittner C, Sheerup DF. 2003. Regulation of opening state and filtration rate in filter-feeding bivalves (Cardium edule, Mytilus edulis, Mya arenaria) in response to low algal concentration. J Exp Mar Biol Ecol. 284: 105-127. doi: 10.1016/S0022- 0981(02)00496-3.

Riisgard HU, Manriquez P. 1997. Filter-feeding in fifteen marine ectoprocts (Bryozoa):particle capture and water pumping. Mar Ecol Prog Ser. 154: 223-239. doi: 10.3354/meps154223.

Riisgard HU, Thomassen S, Jakobsen H, Weeks JM, Larsen PS. 1993. Suspension-feeding in marine sponges halichondria-panicea and haliclona-urceolus - effects of temperature on filtration-rate and energy-cost of pumping. Mar Ecol Prog Ser. 96(2): 177-188.

Roesch C, Lundsgaard-Hansen B, Vonlanthen P, Taverna A, Seehausen O. 2013.

Experimental evidence for trait utility of number in adaptive radiation of a north temperate fish. J Evol Biol. 26:1578-1587. doi: 10.1111/jeb.12166.

Runnstrum S. 1936. A study on the life history and migrations of the Norwegian Spring- Herring based on the analysis of the winter rings and summer zones of the scale. Rep Norwegian Fish Mar Invest. 5(2): 1-107.

Sako R, Tanaka T, Ishimaru T. 2001. House production by Oikopleura dioica (Tunicata, Appendicularia) under laboratory conditions. J Plankton Res. 23(4): 415-423. doi: 10.1093/plankt/23.4.415

Sato M, Kawamata K, Zaslavskaya N, Nakamura A, Ohta T, Nishikiori T, Brykov V, Nagashima K. 2005. Development of microsatellite markers for Japanese scallop (Mizuhopecten yessoensis) and their application to a population genetic study. Mar Biotech. 7:713-728. doi: 10.1007/s10126-004-0127-8.

Schick RS, Urban DL. 2000. Spatial components of the bowhead whale (Balaena mysticetus) distribution in the Alaskan Beaufort Sea. Can J Fish Aquat Sci. 57: 2193-2200.

Schlachler TA, Wooldridge TH. 1995. Small-scale distribution and variability of demersal zooplankton in a shallow, temperate estuary: tidal and depth effects on species- specific heterogeneity. Cah Biol Mar. 36: 211-227.

Schmidt GH. 1983. The hydroid Tubularia larynx causing 'bloom' of the ascidians Ciona intestinalis and Ascidiella aspersa. Mar Ecol Prog Ser. 12: 103-105. doi: 10.3354/meps012103.

69 Schmidt K, Atkinson A, Pond DW, Ireland LC. 2014. Feeding and overwintering of Antarctic krill across its major habitats: The role of sea ice cover, water depth, and phytoplankton abundance. Limnol Oceangr. 59(1): 17-36. doi: 10.4319/lo.2014.59.1.0017.

Segura-Puertas, L, Celis L, Chiaverano L. 2009. Medusozoans (: Cubozoa, Scyphozoa, and Hydrozoa) of the Gulf of Mexico, Pp. 369–379 in Felder, DL, DK Camp (eds.), Gulf of Mexico–Origins, Waters, and Biota. Biodiversity. Texas A&M University Press, College Station, Texas.

Seitz RD, Lipcius RN, Hines AH, Eggleston DB. 2001. Density-dependent predation, habitat variation, and the persistence of marine bivalve prey. Ecol. 82(9): 2435-2451.

Sejr MK, Sand MK, Jensen KT, Petersen JK, Christensen PB, Rysgaard S. 2002. Growth and production of Hiatella arctica (Bivalvia) in a high-Arctic fjord (Young Sound, Northeast Greenland). Mar Ecol Prog Ser. 244: 163-169. doi: 10.3354/meps244163

Shelbourne JE. 1953. The feeding habits of plaice post-larvae in the Southern Bight. J Mar Biol Assoc U.K. 32(1): 149-159). doi: 10.1017/S0025315400011474.

Simon M, Johnson M, Madsen PT. 2012. Keeping momentum with a mouthful of water: behavior and kinematics of humpback whale lunge feeding. J Exp Biol. 215: 3786- 3798. doi:10.1242/jeb.071092

Simon M, Johnson M, Tyack P, Madsen PT. 2009. Behaviour and kinematics of continuous ram filtration in bowhead whales (Balaena mysticetus). Proc. R. Soc. B. 276: 3819- 3828. doi:10.1098/rspb.2009.1135

Sims DW. 1999. Threshold foraging behaviour of basking sharks on zooplankton: life on an energetic knife-edge? Proc R Soc Lond B. 266. 10.1098/rspb.1999.0798.

Sims DW, Southall EJ, Richardson AJ, Reid PC, Metcalfe JD. 2003. Seasonal movements and behaviour of basking sharks from archival tagging: no evidence of winter hibernation. Mar Ecol Prog Ser. 248: 187-196. doi: 10.3354/meps248187.

Skomal GB, Wood G, Caloyianis N. 2004. Archival tagging of a basking shark, Cetorhinus maximus, in the western North Atlantic. J. R. Biol. Ass. U.K. 84: 795-799. doi: 10.1017/S0025315404009968h.

Silina AV. 1995. Mortality of late juvenile and adult stages of the scallop Mizuhopecten yessoensis (Jay). Aquaculture. 141:97-105. doi: 10.1016/0044-8486(95)01207-9.

Smith I, Fonseca MS, Rivera JA, Rittmaster KA. 1968. Habitat value of natural versus recently transplanted eelgrass, Zostera marina, for the bay scallop, Argopecten irradians. Fish Bul. 87:189-196.

70 Snelgrove PVR, Grant J, Pilditch CA. 1999. Habitat selection and adult-larvae interactions in settling larvae of soft-shell clam Mya arenaria. Mar Ecol Prog Ser. 182: 149-159. doi: 10.3354/meps182149.

Southward EC. 1957. The distribution of polychaeta in offshore deposits in the Irish Sea. J Mar Biol Assoc UK. 36:49-75. doi: 10.1017/S0025315400017069.

Speirs DC, Gurney WSC, Heath MR, Horbelt W, Wood SN, de Cuevas BA. 2006. Ocean- scale modelling of the distribution, abundance, and seasonal dynamics of the copepod Calanus finmarchicus. Mar Ecol Prog Ser. 313: 173-192. doi: 10.3354/meps313173.

Stewart JD, Hoyos-Padilla EM, Kumli KR, Rubin RD. 2016. Deep-water feeding and behavioral plasticity in Manta birostris revealed by archival tags and submersible observations. Zoology. http://dx.doi.org/10.1016/j.zool.2016.05.010.

Stone RP. Coral habitat in the Aleutian Islands of Alaska: depth distribution, fine-scale species associations, and fisheries interactions. Coral Reefs. 25: 229-238. doi: 10.1007/s00338-006-0091-z.

Strasser M. 1999. Mya arenaria - an ancient invader of the North Sea coast. Helgoländer Meeresuntersuchungen. 52:309-324.

Sutherland KR, Madin LP, Stocker R, Koehl MAR. 2010. Filtration of submicrometer particles by pelagic tunicates. Proc Nat Acad Sci. 107(34): 15129-15134. doi: 10.1073/pnas.1003599107.

Sutherland SJ, Macy WK, Durbin AG. 1995. Effects of light intensityon the filter-feeding rates of mackerel (scomber scombrus). ICES. H:31.

Tamura T, Konishi K. 2009. Feeding habits and prey consumption of Antarctic minke whale (Balaenoptera bonaerensis) in the Southern Ocean. J Northw Atl Fish Sci, 42: 13–25. doi:10.2960/J.v42.m652

Tay TS, Low JKY. 2016. diversity in the subtidal non- habitats of Singapore. Raffles Bul Zool. 34: 659-665.

Taylor JG. 2007. Ram filter-feeding and nocturnal feeding of whale sharks (Rhincondon typus) at Ningaloo Reef, Western Australia. Fisheries Research. 84:65-70. doi: 10.1016/j.fishres.2006.11.014

Taylor ML, Gwinnett C, Robinson LF, Woodall LC. 2016. Plastic microfibre ingestion by deep-sea organisms. Sci Rep 6:33997. doi:10.1038/srep33997.

Therriault TW, Horberg LM. 2008. Predicting the potential distribution of the vase tunicate Ciona intestinalis in Canadian waters: informing a risk assessment. ICES J Mar Sci. 65: 788–794. doi: 10.1093/icesjms/fsn054.

71

Thiele D, Chester ET, Moore SE, Sirovic A, Hildebrand JA, Friedlaender AS. 2004. Seasonal variability in whale encounters in the Western Antarctic Peninsula. Deep Sea Res. 51:2311-2325. doi: 10.1016/j.dsr2.2004.07.007

Tomita M, Shiga N, Ikeda T. 2003. Seasonal occurrence and vertical distribution of appendicularians in Toyama Bay, southern Japan Sea. J Plankton Res. 25(6): 579- 589. doi: 10.1093/plankt/25.6.579.

Tunberg B. 1985. Studies on the population of ecology of Upogebia deltaura (Leach) (Crustacea, Thalassinidea). Estuar Coast Shelf Sci. 22: 753-765. doi: 10.1016/0272- 7714(86)90097-1.

Ulf L, Baden SP. 1997. Type three functional response in filter feeding of the burrowing shrimp Upogebia deltaura (Leach). Ophelia. 47:1, 33-41. doi: 10.1080/00785326.1997.10433388 van den Berg C. 1993. Filter-feeding in common bream (Abramis brama), white bream (Blicca bjoerkna) and roach (Rutilus rutilus); structures, functions and ecological significance. Landbouwuniversitet Waginingen. van der Hoop JM, Nousek-McGregor AE, Nowacek DP, Parks SE, Tyack PE, Madsen PT. 2019. Foraging rates of ram‐filtering North Atlantic right whales. Func Ecol. 33(7): 1290-1306. 10.1111/1365-2435.13357.

Velez C, Figueira E, Soares A, Freitas R. 2015. Spatial distribution and patterns in three clam populations from a low contaminated ecosystem. Estuar Coast Shelf Sci. 155: 114-125. doi: 10.1016/j.ecss.2015.01.004.

Vethaak AD, Cronie RJA, Van Soest RWM. 1982. Ecology and distribution of two sympatric, closely related sponge species, Halichondria panicea (Pallas,1766) and H. bowerbanki Burton, 1930(Porifera, Demospongiae), with remarks on their speciation. 52(2): 52-102.

Vikingsson GA, Pike DG, Desportes G, Øien N, Gunnlaugsson T, Bloche D. 2009. Distribution and abundance of fin whales (Balaenoptera physalus) in the Northeast and Central Atlantic as inferred from the North Atlantic Sightings Surveys 1987- 2001. NAMMCO Sci. Publ. 7:49-72. doi: 10.7557/3.2705.

Watkins WA, Schevill WE. 1979. Aerial observations of feeding behavior in four baleen whales: Eubalaena glacialis, Balaenoptera borealis, Megaptera novaeangliae, and Balaenopter physalus. J Mammal. 60(1): 155-163. doi: 10.2307/1379766.

Webb P, Wooldridge T, Schlacher T. 1997. Osmoregulation and Spatial Distribution in Four Species of Mysid . Comp Biochem Physiol. 117(4): 427-431. doi: 10.1016/S0300-9629(96)00235-6.

72

Werding B, Hiller A, Lemaitre R. 2003. Geographic and depth distributional patterns of western Atlantic Porcellanidae (Crustacea: Decapoda: Anomura), with an updated list of species. Memoris of Museum Victoria. 60:79-85.

Wieczorek AM, Morrison L, Croot PL, Allcock AL, MacLoughlin E, Savard O, Brownlow H, Doyle, TK. 2018. Frequency of microplastics in mesopelagic fishes from the Northwest Atlantic. Mar Pollut 5(39). doi:10.3389/fmars.2018.00039

Wildish DJ, Miyares MP. 1990. Filtration rate of blue mussels as a function of flow velocity: reliminary experiments. J Exp Marine Biol Ecol 142(3): 213-219. doi: 10.1016/0022- 0981(90)90092-Q.

Winter, J.E. 1969. Über den Einfluß der Nahrungskonzentration und anderer Faktoren auf Filtrierleistung und Nahrungsausnutzung der Muscheln Arctica islandica und Modiolus modiolus. Mar Biol. 4: 87. doi: 10.1007/BF00347037.

Witbaard R, Bergman MJN. 2003. The distribution and population structure of the bivalve Arctica islandica L. in the North Sea: what possible factors are involved? J Sea Res. 50:11-25. doi: 10.1016/S1385-1101(03)00039-X.

Wlodarska-Kowalczuk M. 2007. Molluscs in Kongsfjorden (Spitsbergen, Svalbard): a species list and patterns of distribution and diversity, Polar Research, 26:1, 48-63, doi: 10.1111/j.1751-8369.2007.00003.x.

Woodall LC, Sanchez-Vidal A, Canals M, Paterson GL, Coppock R, Sleight V, Calafat A, Rogers AD, Narayanaswamy BE, Thompson RC. 2014. The deep sea is a major sink for microplastic debris. R Soc Open Sci 1(4): 140317. doi:10.1098/rsos.140317.

Wooldrige TH. 1986. Distribution, population dynamics and estimates of production for the estuarine mysid, Rhopalophthalmus terranatalis. Estuar Coast Shelf Sci. 23(2): 205- 223. doi: 10.1016/0272-7714(86)90055-7.

Wuellner MR, Graeb BDS, Ward MJ, Willis DW. 2008. Review of gizzard shad population dynamics at the Northwestern ddge of its range. Am Fish Soc Symp. 62:637-653.

Würsig B, Dorsey EM, Richardson WJ, Wells RS. 1989. Feeding, aerial, and play behaviour of the bowhead whale, Balaena mysticetus, summering in the Beaufort Sea. Aquatic Mammals. 15(1): 27-37. doi: 10.1016/0006-3207(85)90111-9.

Yahel G, Whitney F, Reiswig HM, Eerkes-Medrano DI, Leys SP. 2007. In situ feeding and metabolism of glass sponges (Hexactinellida, Porifera) studied in a deep temperate fjord with a remotely operated submersible. Limnol Oceanogr. 52(1): 428-440. doi: 10.4319/lo.2007.52.1.0428.

Yamamoto, G. 1968. Studies on the propagation of the scallop, Patinopecten yessoensis(Jay),

73 in Mutsu Bay. Fish Res Board Can Transl Ser 1054: 68.

Zhao Z, Dong S, Wang F, Tian X, Gao Q. 2014. Effect of density on breathing and feeding of filter-feeding silver carp (Hypophthalmichthys molitrix Val.). Aquaculture. 433: 133-136. doi: 10.1016/j.aquaculture.2014.05.043.

74