Distribution Patterns of Microplastics in Seawater Surface at a Portuguese Estuary and Marine Park

Measuring local levels of marine pollution by microplastics (MP) and identifying potential sources in coastal areas is essential to evaluate the associated impacts to environment and biota. The accumulation of floating MP at the sea surface is of great concern as the neustonic habitat consists of a feeding ground for primary consumers (including filter-feeders) and active predators, which makes these organisms a relevant via of MP input into the marine trophic chain. Here, a baseline evaluation of MP accumulation at the sea surface was conducted with a neuston net (335 μm mesh) at the Arrábida coastal area, in Portugal. The study site encompasses a marine protected area and an estuary, both under strong anthropogenic pressures due to multiple activities taking place. A short-term investigation on local spatiotemporal distribution, concentration and composition of MP was performed for the first time, through the monthly collection (summer 2018 to winter 2019) of samples at 6 stations. All the neuston samples contained MP and their mean concentration was 0.45 ± 0.52 items m−3 (mean ± SD). Both the averaged MP:neuston and MP:ichthyoplankton ratios were higher in December, when concentrations of organisms decreased. Temporal distribution patterns followed expected trends, as MP concentration was clearly higher in winter months due to precipitation and runoff. Although mean MP concentrations did not vary significantly between sampling stations, there was a spatial distribution of MP in relation to particle shape and size. Fragments were the most abundant shape and MP belonging to 1–2 mm size class were dominant. Amongst a diversity of 10 polymers identified by FTIR analysis, polyethylene (PE), polypropylene (PP) and copolymer PP/PE were the most abundant. Potential links between local sources/activities and the different polymers were suggested. Altogether, the information provided in this study aims to raise awareness among the identified sectors and consequently to act toward the prevention of MP inputs in the region.

Measuring local levels of marine pollution by microplastics (MP) and identifying potential sources in coastal areas is essential to evaluate the associated impacts to environment and biota. The accumulation of floating MP at the sea surface is of great concern as the neustonic habitat consists of a feeding ground for primary consumers (including filter-feeders) and active predators, which makes these organisms a relevant via of MP input into the marine trophic chain. Here, a baseline evaluation of MP accumulation at the sea surface was conducted with a neuston net (335 µm mesh) at the Arrábida coastal area, in Portugal. The study site encompasses a marine protected area and an estuary, both under strong anthropogenic pressures due to multiple activities taking place. A short-term investigation on local spatiotemporal distribution, concentration and composition of MP was performed for the first time, through the monthly collection (summer 2018 to winter 2019) of samples at 6 stations. All the neuston samples contained MP and their mean concentration was 0.45 ± 0.52 items m −3 (mean ± SD). Both the averaged MP:neuston and MP:ichthyoplankton ratios were higher in December, when concentrations of organisms decreased. Temporal distribution patterns followed expected trends, as MP concentration was clearly higher in winter months due to precipitation and runoff. Although mean MP concentrations did not vary significantly between sampling stations, there was a spatial distribution of MP in relation to particle shape and size. Fragments were the most abundant shape and MP belonging to 1-2 mm size class were dominant. Amongst a diversity of 10 polymers identified by FTIR analysis, polyethylene (PE), polypropylene (PP) and copolymer PP/PE were the most abundant. Potential links between local sources/activities and the different polymers were suggested. Altogether, the information provided in this study aims to raise awareness among the identified sectors and consequently to act toward the prevention of MP inputs in the region.
Keywords: microplastics, distribution, Sado estuary, marine park, MP:neuston ratio, MP:ichthyoplankton ratio, plastic polymers, Portugal INTRODUCTION Tackling marine plastic pollution became a major planetary challenge of the 21st century. Besides the worldwide scientific contribution to the topic for more than one decade and the increasing public awareness, governments have proven their commitment by implementing more sustainable measures and encouraging both initiatives and changes (European Commission, 2018;UNEP, 2018). Yet, although plastic production has recently decreased in Europe from 64.4 MT in 2017 to 61.8 MT in 2018, it has continued to grow at a global level, from 348 MT in 2017 to 359 MT in 2018(PlasticsEurope, 2019. This tendency largely relies on the persistent high demand for such a low-cost, lightweight, versatile, and durable material (Barnes et al., 2009). Consequently, and adding to the excessive consumption of disposable items (Napper et al., 2015) and poor waste management (Frias et al., 2014), plastic pollution represents a significant threat to the marine environment (Laskar and Kumar, 2019). The latest estimations pointed out between 1.1 and 8.8 MT of mismanaged plastic waste being generated annually by land-based human activities at each country (Jambeck et al., 2015). From this waste amount, a considerable part ends up in the marine environment, mainly through wastewater treatment plants discharges, land runoff or transported by the wind, rivers and tides (Andrady, 2011;Jambeck et al., 2015). Despite the greater relevance of terrestrial sources, there are several sea-based activities, such as fishing, aquaculture, maritime traffic, offshore platforms and recreational uses, which may also be considered as additional sources of plastic pollution (Browne et al., 2011;Jambeck et al., 2015;UNEP, 2016;Gewert et al., 2017).
Pioneer studies focused on plastic debris abundance and distribution in the marine environment inevitably verified that plastic pollution could act at a wide size range (from macro to nanoplastics), at a broad spectrum of impacts, as skin injuries or smothering from entanglement, gastrointestinal tract lesions or blockage from ingestion, and event act as vectors of pathogens and chemicals (Laist, 1987;Teuten et al., 2009;Kühn et al., 2015;Bowley et al., 2020). Indeed, the potential of smaller plastics to be ingested by marine biota (Barnes et al., 2009) and to be transferred throughout the trophic chain (Eriksson and Burton, 2003;Farrell and Nelson, 2013;Setälä et al., 2014), was rapidly recognized. This perception shifted the focus of investigation onto microplastics [hereafter MP; defined as particles between 1 µm and 5 mm (Arthur et al., 2009)] which developed into a new research topic addressed worldwide. In addition, the critical concern about the potential impacts of MP in human health through oral, dermal and inhalation exposure has triggered an increase of investigation on this subject although it remains poorly understood Galloway, 2015;Revel et al., 2018).
Regarding the origin of MP, it was considered to be either primary, if manufactured in microscopic size ranges (as industrial pellets and abrasives or microbeads from personal care products); or secondary, if resulting from fragmentation of larger objects (fishing gear, packaging, fibers from synthetic textile washing, paint flakes from nautical coating and dust from vehicle tire) (Cole et al., 2011;Hidalgo-Ruz et al., 2012;GESAMP, 2016;Rochman et al., 2019). The fragmentation of plastic may occur by photo-degradation, mechanical, chemical, and biological action (Barnes et al., 2009;Andrady, 2011;van Sebille et al., 2015). Regardless of its origin, an evident spatial distribution of MP in the water column occurs vertically, from the water surface to the seabed (Thompson, 2004). This mainly relies on polymers density and biofouling level, as both affect particles buoyancy (Gregory, 2009;Kaiser et al., 2017). Horizontal distribution of MP is also known to occur as a result of hydrodynamic forces, mainly by wind (Kukulka et al., 2012), tides, waves and thermohaline gradients (Zhang, 2017). In addition, both vertical (Choy et al., 2019) and horizontal distribution may be influenced by biota, through ingestion and egestion of MP in different compartments of the water column and different locations.
Regarding the impacts on marine biota upon MP ingestion, besides physical harm [e.g., damage in the gastrointestinal tract with inflammatory responses (von Moos et al., 2012) or false sense of satiation (Kühn et al., 2015)], toxicological effects have also been reported (Rochman et al., 2013;Wright et al., 2013). These rely on potential load of harmful chemicals adsorbed from seawater onto plastic and on the toxic additives incorporated during manufacture (Teuten et al., 2009). As a result, MP are suggested to act as trophic vectors of contaminants (Teuten et al., 2007;Garcia-Garin et al., 2020), although their contribution for bioaccumulation (and bioamplification) in organisms tissues may not be as relevant as other contamination pathways, such as prey ingestion or dermal uptake (Koelmans et al., 2016).
Understanding the exposure of primary consumers to MP became essential to evaluate the consequent implications in the marine trophic chain (including eventual detrimental impacts on human health due to seafood contamination). This triggered an increase in research aiming at calculating encounter rates between MP and primary consumers, based on their concentrations and ratio (Collignon et al., 2012(Collignon et al., , 2014Hitchcock and Mitrovic, 2019). Yet, such research has been scarcely conducted in Portuguese waters (Frias et al., 2014;Rodrigues et al., 2019), being insufficient for a country where fisheries have a large cultural and social importance (FAO, 2017) and where seafood constitutes a very important diet component (FAO, 2010;Almeida et al., 2015;EUMOFA, 2020).
Both the Sado estuary and Professor Luiz Saldanha Marine Park, located at the Portuguese west coast, are important nursery areas for fish larvae (Borges et al., 2007(Borges et al., , 2009) and constitute valuable artisanal fishing grounds (Horta e Costa et al., 2013a,b;Batista et al., 2015). However, there are multiple anthropogenic activities taking place at this coastal zone, potentially contributing to local and regional MP pollution and thus posing a threat to this hotspot of biodiversity (Cunha et al., 2014). In this context, this study aims to contribute with baseline data on MP pollution at a Portuguese estuary and marine park by assessing temporal and spatial variations in concentration, distribution, and composition of MP particles. Two hypotheses are tested: (i) MP concentration decreases at sampling stations far away from the metropolitan area of Setúbal; and (ii) MP concentration increase in winter months when compared to summer and autumn months. Additionally, since these are important nursery areas for fish larvae, the ratio of MP to ichthyoplankton was calculated separately from the MP to neuston ratio, aiming to be useful either for comparing with other regions or as a simple and clear take-home message at science outreach activities.

Study Area
The study area, located on the west coast of Portugal, encompasses the south-facing coastal area between the city of Setúbal and the village of Sesimbra (Figure 1). It comprises both the mouth of Sado estuary [designated as the transitional water body Sado-WB1 (ARH Alentejo, 2012)] and the Professor Luiz Saldanha Marine Park (from its eastern side -Figueirinha beach -until the buffer area contiguous to Sesimbra). The mesotidal homogeneous Sado estuary has a mainly tidally driven flow (Martins et al., 2002) with an annual average flow of 40 m 3 s −1 (Vale et al., 1993). It is under considerable anthropogenic pressure due to numerous activities (mostly occurring in its northern margin), from urban and industrial (including maritime traffic), to agriculture and animal production, fisheries and tourism sectors (APA, 2016). Nevertheless, this estuary (Ramsar site no. 826) encompasses a Nature Reserve [RNES (ICNF Reserva Natural do Estuário do Sado)], where birds and habitats are, respectively, protected by a Special Protection Area (PTZPE0011; Birds Directive) and a site classified under the Habitats Directive (PTCON0011), both belonging to Natura 2000 network. Located outside these conservation areas, the Sado-WB1 is adjacent to the city of Setúbal (ca 119.000 inhabitants: Statistics Portugal 2013), being close to a multipurpose terminal port, ship repair yard and to the submarine outfall of Setúbal wastewater treatment plant. Along the Sado-WB1 margins, where two streams discharge (Comenda and Livramento), there are diverse nautical and tourist facilities, an important cement industry, an Orthopedic Hospital and beaches of high demand.
Established westward from the estuary, is Professor Luiz Saldanha Marine Park (hereafter marine park), a sheltered coastline from the prevailing north and north-west winds by the Arrábida mountain chain (Henriques et al., 1999). The different protective measures established in this marine park aim to minimize the impacts of nautical, recreational and fishing activities on its biological and ecological patrimony (Henriques et al., 1999). The erosion of the adjacent cliffs of Arrábida (Gonçalves et al., 2002;Costa et al., 2013) contribute to the complex substratum found in this subtidal rocky reef, which is expected to export MP to the adjacent marine environment due to the breakdown of larger items through physical abrasion on rocks (Eriksson and Burton, 2003;Cheshire et al., 2009). One of 2 submarine outfalls regularly used for effluent discharges of wastewater treatment plants of Sesimbra Municipality (ca. 50.000 inhabitants; Statistics Portugal 2013) is located within the study area.

Sampling Methods
Six sampling campaigns were conducted from August 2018 to February 2019 (summer to winter), at 6 stations ( Table 1). These were located at the 5 m isobaths and distributed 5 km apart from each other, from the mouth of Sado estuary through the marine park (Figure 1). In each station, a 30 min neuston trawl (following Directive, 2013) was performed in the E-W direction, at a constant speed of 1-3 knots. Initial and final GPS positions were registered and enabled trawl length and area calculations to allow posterior standardization of MP data (following Law et al., 2014). Sampling campaigns were specifically scheduled to days with calm weather conditions (Beaufort wind scale ≤3) and tows were performed out of the vessel wake zone (ca 25 m behind the vessel). Precautions intended to reduce vertical mixing of buoyant plastic particles and consequently increase the efficiency of the selected equipment (neuston net). The 3 m long neuston net (Aquatic Biotechnology) had a stainless steel 0.8 × 0.3 m (width × height) rectangular opening and a 335 µm polyamide mesh. Its floatation system assured that only half of the opening frame was submerged (therefore collecting MP floating in the top 15 cm of the water column). The flowmeter (Hydro-bios) attached to the lower third of the net opening enabled the calculation of the volume of filtered water. As only half of the net opening is submerged, the volume was calculated with the following formulae: Tow length = flowmeter revolutions × hydraulic pitch Following each tow, the content in the cod end container was thoroughly poured into a 250 µm stainless steel mesh sieve (where larger pieces of biological material as sticks, seagrass leaves and algae, were rinsed with filtered seawater before being discarded) and then stored in glass jars. A small aliquot (ca 50 ml) per sample was collected and preserved separately, in 100 ml of 70% ethanol, to allow the identification of neustonic organisms and the calculation of the MP:neuston and MP:ichthyoplankton ratios. The neuston samples (n = 36) were transported in ice coolers to the laboratory and then frozen at −20 • C.

Sample Processing and Microplastics Characterization
Due to the considerable volume of biological material present, samples were processed according to Gago et al. (2018). After thawing, the sample was transferred to a 2 L glass beaker where the biovolume was measured after 1 h of sedimentation. Then, the organic content digestion was performed by adding a 10% KOH solution, with volume equivalent to at least 3 times the sample biovolume. Following the 48 h of digestion at room temperature, density separation was conducted by adding 1 L of a hypersaturated NaCl solution (1.2 g cm −3 ). After manual stirring, it was left to settle for 1 h before filtration of the supernatant with a vacuum filtration system. After filtration of every 500 ml (approx.), the sample was stirred and allowed to settle again before the next filtration. Each filter (MFV2 glass fiber filter with 47 mm Ø and 1,0 µm pore; FILTER-LAB) was stored in a covered Petri dish until observation under a stereoscopic microscope (Leica MZ12.5) equipped with a camera (MOTICAM 10+). Particles were measured with the Motic Images Plus 3.0 software, considering the 0.335-5 mm size range (the lower limit corresponds to mesh size of the neuston net) and then attributed to one of the following size classes: 0.335-1, 1-2, 2-3, 3-4, and 4-5 mm. Characterization consisted of registering both color and shape. Particles were assigned to one of six shapes: fragment, film, foam, fiber, filament, and bead [ Table 2; adapted from Lusher et al. (2017)]. The particles selected to follow polymer identification were isolated in covered concave slides. MP concentration was reported as items m −3 and items km −2 to enable comparisons with similar studies.

Polymer Identification
Selection of particles for polymer identification, from all shapes ( Table 3), was based on the best expert judgment according to similarity, texture, thickness, shine and reaction to touch (following Lusher et al., 2017). Polymer identification was achieved by Fourier Transformed Infrared Spectroscopy (FTIR). The majority of the particles (mainly between 1 and 5 mm) were analyzed in attenuated total reflectance (ATR) mode. Spectra were acquired using an Agilent Handheld 4300 FTIR Spectrometer with a DTGS detector, with controlled temperature and a diamond ATR sample interface; the analyses were performed at the sample surface. Spectra were acquired with a resolution of 4 cm −1 and 32 scans. For fibers and smaller particles (mainly at the 0.335-1 mm size range), analyses were carried out in a Nicolet Nexus spectrophotometer coupled to a Continuµm microscope (15x objective) with an MCT detector. Spectra were collected in transmission mode, with a resolution of 8 cm −1 and 128 scans. The spectra are shown here as acquired, without corrections or any further manipulations, except for the occasional removal of the CO 2 absorption at ca. 2,300-2,400 cm −1 . The identification of polymers was first made by searching in the extensive polymer spectral database of the Department of Conservation and Restoration (FCT NOVA) and the assignments were confirmed by analysis of the polymers characteristic bands (Hummel, 2002).

Quality Assurance (QA) and Quality Control (QC)
The airborne contamination was analyzed by exposing wet filters to the air (procedural controls; blanks), both during field (inside a hanging open glass jar, at the boat deck, one per sampling campaign) and lab work (inside Petri dishes, one at the left and one at the right of the working area, per group of 3 samples). All the fibers extracted from a sample which were similar to those found in the respective blanks (from field and lab work) were excluded from results. Sources of contamination were also minimized both during field and lab work by using glass, stainless steel and aluminum materials. Samples were kept covered at all times, both cotton lab coat and nitrile gloves were always worn, and benches and equipment were rinsed before use with Milli-Q filtered water and ethanol.

MP:neuston and MP:ichthyoplankton Ratios
The biovolume of neuston aliquots was registered after 1 h of sedimentation in the graduated cylinders and then homogenized (manual stirring). Three subsamples of 2 ml each were analyzed under a stereomicroscope using a Bogorov counting chamber. Apart from insects, neuston organisms mainly consisted of zooplankton. Dominant groups (fish larvae and eggs, Mysidacea, Polychaeta, Chaetognata, Apendiculata, Bivalvia larvae, zoea and megalopa of Brachyura, Cladocera, naupli of Cirripedia, Copepoda, Echinodermata larvae, Amphipoda, Isopoda and Insecta), rather than individual species or genera (Di Mauro et al., 2017), were counted with the support of a hand tally counter, enabling the calculation of each group abundance. Mean counts (all dominant groups were considered for MP:neuston ratio, whereas only fish larvae and eggs were considered for MP:ichthyoplankton ratio calculation) were extrapolated according to the aliquot and sample biovolume and then converted to individuals m −3 .

Statistical Analysis
To evaluate how the MP:neuston ratio varied temporally (along 6 months) and spatially (between the 6 stations), a Kruskal-Wallis test was performed. This non-parametric test, conducted after the invalidation of parametric assumptions, was followed by posthoc multiple comparisons with the Dunn's test. The same tests were applied for MP:ichthyoplankton ratio.
A two-way ANOVA without replication was performed to assess whether temporal (6 campaigns) and spatial (6 stations) variation occurred in MP concentration (dependent variable). This parametric test was used after Box-Cox transformation of original data to meet normality (Shapiro-Wilk test) and homogeneity of variances (Levene test) assumptions. Posthoc Bonferroni's test (p < 0.05) were used to identify the sources of significant differences. Analysis were conducted in Statistica 13 (Statsoft) software.
The effect of campaigns and stations (fixed factors; with 6 levels each) in MP concentration of each particle shape (multivariate data) was tested by a permutational multivariate analysis of variance (PERMANOVA), with 999 permutations. Data were square-root transformed and the resemblance matrix between samples was calculated based on Bray-Curtis similarities. When differences were statistically significant, pair-wise comparisons among levels were analyzed. Then, to determine which particle shape most contributed to explain the dissimilarity amongst each pair of samples, the similarity percentages routine (SIMPER; with a cut-off percentage of 90% for low contributions) was conducted. These statistical procedures, which were conducted in the Primer 6 software with the Permanova+ add-on (Clarke and Gorley, 2006;Anderson et al., 2008), were similarly applied to understand the response of MP concentration of each size class to both factors (campaigns and stations).

Presence and Concentration of Microplastics
From the total of particles (3,317) extracted from the 36 neuston samples, 353 (11%) were discarded for being considered airborne contamination fibers and 44 (1%) were excluded after  being identified as non-plastic particles by Fourier Transformed Infrared Spectroscopy. Therefore, the assessment of the temporal and spatial distribution of MP (size range 0.335 to 5 mm) was based on a total of 2,920 particles. All samples contained MP, with a mean concentration of 0.45 ± 0.52 items m −3 (mean ± SD) and 40,822.58 ± 43,578.63 items km −2 . While the highest concentration per cubic meter was found in February at Figueirinha beach (St2; 2.06 items m −3 ), the highest concentration per square kilometer was verified at Setúbal (St1; 203558.50 items km −2 ). Conversely, the lowest concentration (0.04 items m −3 or 2,068.85 items km −2 ) was observed at Mijona beach (St6) in October (Figure 2). The number of MP ranged from 405 at St1 in February to 5 MP at St6 in October.

Fourier Transformed Infrared Spectroscopy (FTIR) Analysis
Among 265 particles confirmed as microplastics by FTIR analysis, a total of 10 polymers were identified (

Temporal and Spatial Distribution
MP concentration in February was significantly higher than those found in all other campaigns (Figure 6A), except for January (Bonferroni test, p < 0.05). In October, concentration was the lowest and significantly different from January. MP concentration did not vary significantly between stations (p = 0.06; Figure 6B).
The PERMANOVA results showed significant differences in the  MP concentration of each particle shape, between sampling campaigns (Pseudo-F = 6.57, P(perm) = 0.001) and stations (Pseudo-F = 2.11, P(perm) = 0.008). MP concentrations per particle shape differed mainly between October and February, but also between each of these 2 months and all the other campaigns. The combination of the 3 predominant shapes: fragments, films and foams contributed with more than 71% (cumulative percentage) for the dissimilarities between all pairs, with concentrations being always higher in February ( Figure 8A). An additional result from pair-wise comparisons concerned the dissimilarities between November and January campaigns, which were based on the higher concentration of fragments, foams and beads (with contributions of 27.87, 24.54, and 14.28%, respectively), found in January. Regarding spatial variation, the relatively higher foam and bead concentrations at St1 (estuary) explained dissimilarities found between this station and both St5 (Sesimbra; foam-24.82% and bead-14.00% contributions) and St6 (Mijona beach; foam-27.44% and bead-13.26% contributions) ( Figure 8B). Moreover, at St6 concentrations of fragments were significantly lower than St3 (Portinho da Arrábida). No plastic pellets were collected in this study and all beads belonged to the smaller size class (0.335-1 mm).

Distribution Variations According to Size Class
By decreasing order, the relative abundance of each size class (mm) was: 1-2 (36%) > 2-3 (24%) > 3-4 (16%) > 0.335-1 (15%) > 4-5 (9%). According to PERMANOVA results, MP concentration varied according to size class between campaigns (Pseudo-F = 7.69, P(perm) = 0.001) and stations (Pseudo-F = 2.55, P(perm) = 0.005). MP belonging to the 1-2 and 2-3 mm size classes explained (with more than 46% of cumulative contribution) the dissimilarities found between February and all the other campaigns and also between January and both November and October months ( Figure 9A). In addition, while the 0.335-1 mm size class largely contributed (ca. 29%) to distinguish August from October (being more represented in August), the higher concentration of MP at the 1-2 and 3-4 mm class sizes in December, compared to October, contributed more than 48% for their differences. The particle size range at St6 (the furthest station from the estuary) was distinct from all the others, mainly due to its low concentration of MP belonging to the 3-4 mm size (contributions between 21 and 25%) and particularly different from St1 and St4 due to the smaller concentration of MP at the 4-5 and 1-2 mm size ranges, respectively ( Figure 9B).

Presence and Mean Concentration of Microplastics -Comparison With Other Studies
The presence of MP in all coastal samples collected in this study is in accordance to reported MP pollution levels close to shore and to estuaries, either at Portuguese (Frias et al., 2014), European (Pedrotti et al., 2016;Frère et al., 2017) In addition, if compared with surface waters of estuaries and contiguous coastal areas from other countries, our study area presents higher MP concentrations than those quantified by Lima et al. (2014) at the Goiana estuary in Brazil (0.26 items m −3 ). Conversely, mean MP concentration at Arrábida was more than one order of magnitude lower than the 3 estuaries in Australia east-coast investigated by Hitchcock and Mitrovic (2019) (with a range of 23-198 items m −3 at the Clyde estuary, the one with the lowest MP concentration). However, comparisons must be performed cautiously as local environmental conditions, levels of anthropogenic pressure, and methodologies applied may differ among studies (Lima et al., 2014). In fact, the lack of methodologies standardization has been often  highlighted (Gago et al., 2018;GESAMP, 2019) and remains a current challenge.

MP:neuston and MP:ichthyoplankton Ratios
The accumulation of floating MP at the seawater surface layer leads to concerns about the exposure of neustonic organisms, such as zooplankton (including ichthyoplankton), to these synthetic particles and, consequently, of their active predators and filter-feeding biota (Collignon et al., 2012).
As expected, the increasing tendency of MP concentration observed in winter months and the simultaneous decline of zooplankton and larval fish abundance (Cunha, 1993;Primo et al., 2011) increased both MP:neuston and MP:ichthyoplankton ratios in this time of the year. Regarding the MP:ichthyoplankton ratio, although MP have never exceeded ichthyoplankton in number at any sample, their similar proportions suggest a higher potential for MP to be ingested either by fish larvae or by ichthyoplankton's predators [crustaceans: crabs, shrimps, euphausiids, amphipods, and copepods; ctenophores; fishes; medusae (Bailey and Houde, 1989;Paradis et al., 1996)]. As a consequence, it would be expected to find critical variations in the following spring at these important nursery areas. Further studies focused on MP:ichtyoplankton ratio and on MP ingestion by wild fish larvae would be essential to confirm possible impacts, as their survival largely influence fish recruitment success and population fluctuations (Houde, 1987).
The average MP:neuston ratio verified in this work (0.0009) was low when compared to other studies: 0.002 at the Bay of Calvi (Collignon et al., 2014) and 0.2 at the Ligurian Sea (Pedrotti et al., 2014). To our knowledge, besides the MP to fish larvae ratio (1.5:1.0; fish eggs excluded) found in Douro river, Portugal (Rodrigues et al., 2019), no studies have assessed the proportion between MP and ichthyoplankton alone (0.091 in average; 1:1.294 maximum), as it has been pooled together with all other zooplankton organisms.

Temporal and Spatial Distribution
Both temporal and spatial distribution variations were verified for MP concentration in our study site. As expected, MP concentrations increased significantly in winter months, achieving a maximum in February. This is in agreement with the reported increase of MP concentrations in marine coastal waters after storms and heavy rainfall, typically frequent in winter season for Mediterranean-type climatic conditions (Santos et al., 2005), which induces frequent floods and increase river discharges (Veerasingam et al., 2016;Gündogdu et al., 2018;Hitchcock, 2020). Regarding the spatial distribution, it was anticipated a clear seaward decrease in MP concentration at stations further away from the metropolitan area of Setúbal (Sado estuary), with an eventual increase at the station close to Sesimbra. Instead, MP pollution level found at stations located between Setúbal and Sesimbra municipalities kept similar orders of magnitude, although with a slight decrease tendency. Such retention of MP, which might be related with the shelter provided by Arrábida mountain chain against the prevailing north and north-west winds, may impact the high biodiversity of this Marine Park. Therefore, the continuous input of MP in the estuary (at St1, the closest station to the urban area of Setúbal) is suggested to partially accumulate in the sheltered Arrábida nearshore area.
Further explanations could rely on the hydrodynamics at the Arrábida rocky reef which may potentially enhance fragmentation of both MP or even larger items, by mechanical action against rocks (Eriksson and Burton, 2003;Cheshire et al., 2009) contributing for the increase of secondary MP. Subsequently, the continuous exportation of these MP by local currents could explain the considerable concentration of particles at st5, despite being distant from the estuary. Concentrations calculated at this station may also result from Sesimbra village input of MP yet, due to the fragmentation potential at the sandy surf zone of this sheltered bay, particles may easily achieve sizes which are not retained by the neuston net.
Lastly, fragmentation enhanced during retention at the Arrábida nearshore may also contribute to export MP in the coastal drift, explaining the unexpected high concentration of MP reported further south by Frias et al. (2014) at Costa Vicentina.

Distribution Variations According to Particle Shape
Bead and foam shapes presented distinct patterns in their distribution at the study area, unlike the other MP shapes. Both were predominantly collected in station 1 (Setúbal), contrasting with station 5 and 6 (Sesimbra and Mijona beach), with concentrations being higher in the January and February campaigns. The preponderance of foam shape (expanded polystyrene) in the estuary is potentially related to fisheries activities, consisting of secondary MP from the breakdown of buoys and cooler boxes for bait and catches, which despite the decrease of fishery activities during winter (DGRM, 2018) are frequently kept close to the seashore and left exposed to adverse weather conditions till the next fishing season. Conversely, beads (primary MP) are suggested to enter in the marine ecosystem by wastewater treatment plant (WWTP) effluents after domestic use (Fendall and Sewell, 2009), as these particles may not be retained in the treatment processes. The predominance of fragments in this study is in line with results from a similar study performed in Portugal, in the Douro estuary (Rodrigues et al., 2019) and in Australia, at Clyde, Bega and Hunter estuaries (Hitchcock and Mitrovic, 2019). This suggests that secondary sources of MP prevail, rather than primary sources, and are related with the diverse activities taking place in the nearby urban area, including littering. Our findings differ from studies reporting fibers as the predominant shape detected (Beer et al., 2018;Bessa et al., 2018), usually attributed to fishing ropes degradation (Ramos et al., 2012) and to the inefficient retention of fibers from textile laundry by the WWTP (Browne et al., 2011). In fact, fibers represented only 4% of the total of MP found in our study, after the exclusion of airborne contamination (11%) from the original MP amount. The small abundance of fibers reported here may rely on the retention efficiency of treatment processes of WWTW (Gies et al., 2018) at both Setúbal (advanced secondary treatment) and Sesimbra (tertiary treatment) Municipalities, or be related to the sampling method applied, as neuston nets are suggested to underestimate the concentrations of fibers when compared with other methods (Barrows et al., 2017;Green et al., 2018).

Distribution Variations According to Size Class
Distribution patterns of MP according to their size were noticed both in time and space. In fact, the predominance of bigger sized MP (3-4 and 4-5 mm size) inside the estuary, the abundance increase of MP in December (beginning of winter), particularly MP belonging to the 3-4 mm size class and the high concentration of MP from intermediate size classes (1-2 and 2-3 mm) at january and mostly in February, suggest that MP inputs in this Portuguese region occur mostly close to Setúbal and mainly consists of larger particles which undergo fragmentation over time.
The preponderance of 1-2 mm sized particles among the 5 size classes, instead of the expected smallest size class (0.335-1 mm), according to Norén (2007) and Kang et al. (2015) findings, may be essentially related with the sampling method used here. As mentioned before, the use of neuston nets may underestimate fibers concentrations, which are more malleable and easier to escape through the net mesh, explaining the low concentrations of fibers collected here and in particular those belonging to our smallest size class (0.335-1 mm). Secondly, the reduced concentration of MP belonging to the 0.335-1 mm size class could be related with the retention time spent at sheltered stations, which could enhance biofouling levels and consequently cause smaller particles to sink (Kaiser et al., 2017) or to be ingested, as biofilms are suggested to increase MP palatability (Vroom et al., 2017).
As several studies have already highlighted (Song et al., 2014;Lenz et al., 2016), MP concentrations at surface waters are potentially underestimated due to the lower size limit of the range considered for monitoring, usually ca. 330 µm (net mesh used). Consequently, as there is a tendency over time for continuous fragmentation of plastic and permanent input into the marine ecosystem, studies are missing the size fraction which is potentially more abundant and easily ingested by primary consumers (Cole et al., 2013). Therefore, the selection of a sampling method that efficiently collects smaller MP in further studies would be required to clarify the abundance patterns found at this coastal area.

Polymer Diversity
Polymer identification of particles in plastic pollution studies is essential to confirm visual identification processes (Löder and Gerdts, 2019), to characterize the diversity of polymers available and to assist in identifying potential local sources, as it will empower authorities and stakeholders to tackle this global concern by implementing efficient prevention measures.
The high polymer diversity (10 polymers) detected mirrors the diverse activities performed in the area, both on land (domestic, commercial, industrial and tourism) and at sea (fishing and recreational activities, intense maritime traffic to shipyards). As expected, polyethylene (PE) and polypropylene (PP) showed higher percentages, since they are widely used in many applications (mainly packaging of consumer goods and single-use items). Nevertheless, there was also a considerable amount of particles identified as copolymer PP/PE, which occurs as an industrial way to recycle both PE and PP by giving origin to other high demanding applications and expanding market options (Graziano et al., 2019), as containers, outdoor decking or sack bags (Aumnate et al., 2019). Polystyrene (PS) particles, can be related to fragmentation of disposable cutlery, cups and Styrofoam R items (expanded and extruded PS, EPS and XPS), which are currently used in fishing activities, in food trays and other disposable items (Farrelly and Shaw, 2017). A note should be mentioned regarding the presence of kaolin in PS and copolymer PP/PE particles, which is used as a filler to improve the strength of the plastic material. Particles of polyvinyl alcohol polymer (PVA), considered of low environmental impact, may have been originated from medical and sanitary devices, as well as from food packaging. In fact, this polymer is considered appropriate for orthopedic applications (Baker et al., 2012) potentially linked to the Orthopedic Hospital located close to St2. Fibers were identified as Rayon, a cellulose-based semi-synthetic fiber frequently found in similar studies (Comnea-Stancu et al., 2017); polyester, widely used in packaging, textile, automotive, medical, electronic, and construction sectors (Camlibel, 2018); and polyamide (PA), predominantly used in fishing nets, but also used in the automotive sector and as a bone tissue scaffold in the medical sector (Winnacker, 2017;Atayeter and Atar, 2018). Polyurethane (PUR) is widely used in coating epoxy resins to protect boat hulls from deterioration and used as rigid foams to insulate boats from extreme temperatures and noise, besides biomedical, construction, and automotive applications (Akindoyo et al., 2016). Finally, polyacrylic acid (PAA) is used in the manufacture of household cleaning products, but also for enhancing the mechanical properties of hydrogels used as biological glues in the medical and tissue engineering sector (GVR -Grand View Reasearch, 2017).

CONCLUSION
As expected, MP pollution in this study was higher during the winter months, co-occurring with the usual decrease of primary consumers abundance in this season. The consequent increase of both MP:neuston and MP:ichthyoplankton ratios suggests therefore a critical time period for marine biota feeding in the neustonic habitat. Regarding MP spatial distribution, instead of a clear decreasing gradient from the estuary (area with higher human impact) to further coastal stations, a slight decline in concentrations was observed, suggesting a retention effect close to the Arrábida shore. Although fragments were the dominant shape, only foam and beads presented distinct variation in space, according to the location of their potential sources (fishing harbor and WWTP submarine outfall). The predominance of particles at the 1-2 mm size range instead of the smaller size range (0.335-1 mm), is suggested to be related with the sampling method used, although further studies would be required to clarify this hypothesis. The diversity of polymers reflects the multiple activities occurring in the estuary and in the marine park, highlighting the urgent need to disseminate findings locally, namely on fishing communities and in tourism, industrial and marine traffic sectors. Sharing scientific findings with society aims to increase public awareness about MP pollution and to inspire actions toward the prevention and reduction of plastic entering the marine environment.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
DR conducted fieldwork sampling, laboratory procedures (MP extraction and characterization, neuston identification), statistical analysis, and wrote the manuscript. JA performed FTIR analysis, collaborated in the discussion and selection of the best method for MP extraction from neuston samples, provided assistance with laboratory procedures and at reporting FTIR analysis, and results. VO performed micro -FTIR analysis, assisted in the interpretation of spectra, and at reporting FTIR analysis and results. PS co-ordinated the study, discussed results, gave important contributions to the writing and to the English review of the text. MC reviewed and made important contributions to the text. All authors contributed to the article and approved the submitted version.