Floating Marine Litter in Eastern Mediterranean From Macro to Microplastics: The Lebanese Coastal Area as a Case Study

Despite emerging and increasing concerns related to marine micro and macroplastics, no systematic surveys have been undertaken yet in the Lebanese marine area. To understand the spatio-temporal variation of plastic litter (macro and microplastics) in the Lebanese marine environment and to determine the sources of pollution, this study investigated the characteristics of plastic pollution in sea surface waters during wet and dry seasons in 22 sites of Beirut and Tyre regions. A total of 23,023 items were identified and assessed according to the shape, color, and concentration; moreover, the risk of microplastics (MPs) contamination was explored based on a risk assessment model. The obtained results demonstrated that the average macroplastics concentration was 0.45 ± 0.6 items/m3. The average microplastics concentration was found to be 20.1 ± 21.8 and 3.78 ± 5.2 items/m3 in spring and fall respectively. During fall, MPs fragments were dominant in Beirut (97%) and Tyre (91%), and no pellets were observed. During spring, filaments were most encountered in Beirut (76.5%). The most dominant marine litter color was blue followed by black and white. The Pollution Load Index (PLI) values showed a moderate contamination of the Lebanese coast with MPs (PLI: 5.79 ± 3.93) except for several sites in Beirut that showed high values of PLI, highlighting the local influence of cities and rivers on MPs concentration. This study serves as an important baseline for understanding the characteristics of the seasonal variation of MPs along the Lebanese marine environment; it will help stakeholders and countries to take proactive and reactive actions to face plastic litter pollution in the Lebanese coastal area.


INTRODUCTION
The Mediterranean Sea is one of the most affected marine environments with regards to marine litter (UNEP, 2015). Even though the type of litter found is highly heterogeneous, plastics are the most abundant material recorded (Cincinelli et al., 2019). According to the Marine Strategy Framework Directive TG10, marine litter is defined as any persistent, manufactured or processed solid material discarded, disposed of or abandoned in the marine and coastal environment (Galgani et al., 2010). The range and scale of impacts from marine litter can target public safety, as well as, the socio-economic and environmental sectors (Ferronato and Torretta, 2019).
The global annual production of plastics has increased substantially in recent decades, paving the way for the occurrence of more plastic litter in the marine environment (Da Costa 2017). About 80% of marine litter generation is mainly continental and enters the marine environment via multiple pathways, such as, rivers and discharged wastewater (Alomar et al., 2016;Kazour et al., 2019).
Among many different types of marine litter, plastic is considered the most harmful to the marine environment and wildlife (Barnes et al., 2009). One of the most appreciated qualities of plastic products is their durability and corrosion resistance (Rios et al., 2010;Szeteiová, 2010). However, these qualities when combined with improper waste management leads to environmental contamination, and impacts terrestrial, freshwater and marine ecosystems. More than an estimated 80% of the total pressure in marine waters comes from land, while 20% comes from the sea (Jambeck et al., 2015). Single-use plastics account for around 50% of total plastics production (Xanthos and Walker, 2017); hence, the packaging sector produces the highest pressure followed by the construction (Boucher and Friot, 2017;Abalansa et al., 2020), and then the agricultural sector. In addition, the combined pressures from seabased fisheries and aquaculture represent the fourth most important economic sector in terms of pressure on the marine ecosystem (Boucher and Friot, 2017;Charter et al., 2018). Finally, the tourism economic sector, related to the use of beaches, has been reported to increase pressure by 40% in the Mediterranean region every summer (Galgani et al., 2010;Murani et al., 2016), taking into consideration that the surrounding countries annually attract about one third of the world tourism (WTO, 2018).
Once in the ocean, plastic tends to stay at or near the surface where the photo-chemical, mechanical and biological processes degrade large items into smaller fragments (less than 5 mm), forming what is called as "microplastics" (Thompson et al., 2004;Andrady, 2011;Isobe et al., 2014). Microplastics (MPs) can result either directly from synthetic fabric residues and cosmetic products (i.e., industrial pellets and scrubbers) or from the breakdown of larger plastic products generating fragments or fibers (Arthur et al., 2009;Browne et al., 2011;Hidalgo-Ruz et al., 2012;Cole et al., 2013). Particles with a density that exceeds that of seawater (1.02 g/cm 3 ) will sink and accumulate in the sediment (Woodall et al., 2015), while low-density particles tend to float on the sea surface or in suspension in the water column (Fossi et al., 2012;Suaria and Aliani, 2014). However, the main threat of microplastics floating on seawater is the possibility of ingestion of their ingestion by species of commercial importance for fisheries and aquaculture, which can obstruct the digestive tract and limit the entry of food, consequently leading to malnutrition for the animal (Wright et al., 2013;Ivar do Sul and Costa 2014;GESAMP, 2016;Karbalaei et al., 2018;Kazour et al., 2019;Xiong et al., 2019). In the marine environment MPs are referred to the 'cocktail of contaminants' due to their association with additives; hence, their ability to reach biota and human (Rochman, 2015;Hartmann et al., 2017). Two important seafood species (European anchovy, Engraulis encrasicolus, and the spiny oysters, Spondylus spinosus) that are wholly consumed by the Lebanese community have high presence of MPs which highlights the potential risks that MPs might have on the biota. The risk of MPs in the environment has been previously assessed using the pollution load index (PLI) Liu et al., 2019). But such assessment remains difficult to determine due to the environmental influences and the absence of a standardized method to monitor MPs in the environment (Hidalgo-Ruz et al., 2012). In addition, PLI has been not assessed in the Mediterranean Sea yet even though this semi-enclosed sea is considered to be a hotspot for plastic contamination (Llorca et al., 2020). The Lebanese coast is affected by a wide range of anthropogenic pressures. Urban sprawl and pollution from wastewater effluents, including river discharges, appear to be the most severe and permanent problems of the coastal zone (UNEP ERML, 2012). In addition, the presence of numerous mismanaged open-air dumpsites and landfills along the coast point out their potential contribution in MPs entry into the marine environment (UNDP, 2017;Kazour et al., 2019). In the Eastern Mediterranean, studies related to marine litter and microplastics in marine waters are relatively scarce for their environmental and social costs (Özden et al., 2021). Such studies are deficient along Lebanese coast, with the exception of a recent study that highlighted the presence of high levels of microplastics in surface waters and sediments compared to other Mediterranean regions (Kazour et al., 2019). Therefore, a baseline study assessing a significant number of samples and covering more seasons is necessary to have an overview of plastic debris and its variations in surface waters along the Lebanese coast. Hence, the present study aims 1) to quantify marine litter in coastal areas under a gradient of anthropogenic pressures from populated coastal zones to Marine Protected Areas and the assessment of spatial differences, 2) to evaluate the MP risks, and 3) to provide a science-based knowledge for the management of marine plastic litter in the Lebanese marine water, in the context of the UN Sustainable Development Goal 14.1 (SDG 14.1).

Case Study Sites
Surveys for marine litter were conducted onshore and offshore sites near Beirut and Tyre along the Lebanese coast. Beirut, the capital, is the most urbanized and industrialized area (population ≈2.4 million according to the World Population Review 2021). It is subject to high coastal pressures due to the presence of the Beirut port, three coastal landfills and numerous sewage outfalls. The capital suffers from solid waste problems, as there are many mismanaged coastal dumpsites polluting the coastal area. In addition, various effluents such as industrial discharges and wastewater from inhabited areas drain along by Beirut River and Antelias River, and are discharged directly into the sea. Tyre region does not suffer from the same anthropogenic pressure (population ≈ 200,000 according to the Lebanese Population -Central Administration of Statistics (CAS, 2014)), despite being the fourth densest coastal city in Lebanon. It is considered to be less polluted due to the lack of industrial activities in the region along with the presence of the Tyre Coast Nature Reserve (TCNR).

Sampling of Marine Litter
Sampling and analysis of macro and microplastics were carried out by the National Centre for Marine Sciences in Lebanon (CNRS-L). Samples were collected from 22 marine sites using a Manta net (W 60 cm x H 26 cm rectangular frame opening; 2.6 m length and net mesh size 200 μm) that was towed with a 7 m catamaran for 10 min at a speed of 3 knots and in low wind conditions (less than 0.4 miles/s as suggested by Kukulka et al., 2012). The number of items for each site is presented in Table 1.
Start-end position points were recorded on GPS. All tows were conducted from the ship's side and beyond the ship's wake. After completion of each tow, the net was washed thoroughly with seawater in order to collect all particles in the cod end (Zeri et al., 2018). The collected water in the cod end was transferred in glass jars containing 70% ethanol solution and were stored in a cooling box until further processing. In the laboratory, the samples were filtered on a stainless-steel sieve of an 850 μm mesh-size. Macroplastics were extracted and rinsed with filtered distilled water and the MPs present on the sieve were sorted visually under the stereomicroscope. In order to extract the MPs fraction below 850 μm, the recovered water was vacuum-filtered on GF/A filters and conserved inside a clean glass Petri dish. The solutions of hydrogen peroxide (H 2 O 2 ) and zinc chloride (ZnCl 2 ) were used to eliminate organic material (prior to microscopic observation and counting) and separate MPs, respectively (Guven et al., 2016;Kazour et al., 2019).
Microplastics were counted and photographed under a dissecting Stereomicroscope (Carl Zeiss, Stemi-305, Germany) equipped with AXIOCAM 105 color 5 MP. Photos were taken and processed with Zein platform (Image Analysis software). When observing particles, the filter inside the glass Petri dish was moved in a "zigzag" pattern from left to right. Microplastics particles were visually identified and assessed based on their shape, color and size (Cózar et al., 2015;Kazour et al., 2019). Two shapes of microplastics were observed: filament and fragment. Based on different colors, microplastics can be divided into ten categories: blue, black, yellow, transparent, white, grey, orange, purple, green and red.
Particles longer than 5 mm were considered as macroplastics and recorded only in fall. Several precautions were followed to minimize air-born contamination during sampling, sample preparation and analysis. Different procedures (filtration, measurements) were done under a laminar flow hood, wearing cotton lab coats. All surfaces and equipment were rinsed before usage with filtered distilled water. Only glass and stainless-steel laboratory supplies were used (Funnel, bottles, Petri dishes, Erlenmeyer flasks, filtration system) and were covered immediately after manipulation. All solutions (distilled water, H 2 O 2 and ZnCl 2 ) were filtered on glass fiber filters (GF/A filters, 1.6 μm porosity, Whatman, France) to eliminate unwanted particles or fibers. The contamination control procedure was followed according to the methodology mentioned by Kazour et al., 2019. During the lab work a quality control procedure was implemented using a Petri dish with distilled wather puted in the laminar flow hood were all procedure was done. After work the Petri dish were examined to make sure there is no macroplastic contamination in the working environment.

Microplastics Characterization and Identification
Visual observation often leads to an overestimation of the microplastics found in environmental and biological samples (Dekiff et al., 2014;Lenz et al., 2015). Micro-Raman Xplora Plus (HORIBA Scientific ® , France) was used to identify microplastics and macroplastics observed in this study. Due to the high time consumption, a subsample that represented 1% of random suspected plastics were characterized. Identification was carried using two lasers of a 532 and 785 nm wavelengths with a range of 200-3,400 cm −1 using two objectives (Olympus):×10 and ×100 (Kazour et al., 2019). The spectra is compared using the polymer identification software (KnowItAll, WILEY ® ) and a personal library made with specific polymers obtained from Goodfellow (France).

Microplastics Risk Assessment Method
To assess the degree of MPs pollution in sea surface of the two Lebanese coastal regions, an integrated Pollution Load Index (PLI) was calculated based on Tomlinson et al. (1980). The PLI is regarded as a standardized rule for monitoring the degree of pollution between different areas (Angulo, 1996). Therefore, this risk assessment model is robust and applicable to better understand and evaluate the ecological risks caused by microplastics pollution Xu et al., 2018;Liu et al., 2019). PLI at each site is related to MPs Concentration Factors (CFi) as given below: CFi of the MPs is the quotient of the MPs concentration at each station (Ci) and the minimum MPs concentration (Coi). The lowest concentration of MPs value detected in this study was considered as a Coi value. PLI area was calculated for both Beirut and Tyre regions.

Statistical Analyses and Maps
The statistical software R Studio (RStudio Team, 2020) was used for all statistical analyses. In addition, all maps were built using the Generic Mapping Tools (GMT), version GMT6 using the GMT coastline database GSHHG to show the Lebanese coast features and (GEBCO, 2020) Gridded Bathymetry Data for the bathymetry contour lines (GEBCO, 2020).
A Shapiro-Wilk normality test was first used to test for the normality of the sample using the shapiro. test function. A Kruskal-Wallis rank sum test was then used to examine the differences between total microplastics concentrations by sites according to several factors, such as, season, region, and sites using the kruskal. test function.
A Principal Component Analysis (PCA) was also performed using the FactoMineR package and the PCAmix function of the PCAmixdata package. PCA provides insight into the structure of the data by synthesizing the relationships between variables. PCA was implemented on the following numerical variables: Total MPs, total filaments and total fragments concentrations and the two categorical variables (season and region) where included in the analysis as additional variables with the following modalities: the 2 modalities of season (spring and fall), and the two modalities of region (Beirut and Tyre). Results are presented on the first two principal components representing the highest percentage of the total variance. Finally, a Pearson correlation test was performed using the cor. test function to test the relationship between total MPs concentration and PLI.

Macro Litter Abundance and Composition
All anthropogenic litter items were recorded but for the purpose of this study we focused only on plastics (Synthetic) which accounted for 99.2% of total collected items ( Figure 1). Floating natural materials (feather, corn, hair, sticks, watermelon seeds and fruits) accounted for 0.8% of total litter items recorded. Out of the 22 transects, 13 were macroplastics-free (59%). Floating macroplastics constituted 3.7% of the plastic litter. A total of 439 plastic items (hygienic pad cover, rope, beer tag, band-aid, fishing net, agricultural plastic fragments) were visually counted and correspond to a concentration of 0.45 ± 0.6 items/m 3 .
Macroplastics were observed at five sites in Tyre (TNR1, TNR2, Jamal 1, Jamal 2, and Jamal 3) with the highest number at Jamal 2 and 3 with 198 and 104 items, respectively. The highest concentrations were observed in Jamal 2 (1.95/m 3 ) followed by Jamal 3 and TNR2 (≈0.60/m 3 ) (Figure 3). Comparing the two regions, macroplastic items constituted 0.75 and 13.5% of the plastic litter in Beirut and Tyre, respectively.

Microplastics Concentrations and Seasonal Variation
A total of 23,023 MPs items were identified in the 22 sites and assessed according to color: black, blue, green, grey orange, purple, red, transparent, white and yellow. The number of MPs in spring and fall was 11,474 and 11,549, respectively (  Figures 4, 5, respectively. Except Costa 1 that showed a higher concentration in fall, all the other coastal sites in Beirut showed high concentrations during spring, with the highest in Manara 2 (109/m 3 ) and Dora 1 (86.7 item/m 3 ) ( Figure 4). As for Tyre, no exceptionally high concentrations were observed for all sites (all sites had a concentration less 25 item/m 3 during both seasons). Similar to Beirut, the highest concentrations were observed during the spring season for all sites, with the highest in TNR2 (24.8/m 3 ), NT1 (21.5/m 3 ), and TNR3 (21/m 3 ).
Microplastics detected in seawater samples from Beirut and Tyre were classified into two main shapes: fragments and filaments. The percentage of fragments was high in Beirut and Tyre during fall, representing 97 and 91% of the total MPs concentration, respectively. During spring, filaments dominated in Beirut with a percentage of 76.5% ( Figure 6). For both areas and seasons, fragments dominated over the filaments with 65.4 and 34.6%, respectively.
To affirm previous results, a PCA was performed using total fragment and filament concentrations as numerical variables, in addition to season and region as categorical variables. Results were presented on the first two principal components representing 67.97% of the total variance The PCA showed that the main variables contributing to axis 1 (41.71%) were total fragment and filament concentrations, as well as, seasons, whereas region contributed to axis 2 (26.26%). The PCA outcomes showed that total fragment and filament concentrations (summed up to become total MPs concentration) followed seasonal variability but not a spatial one (Figure 7).

Microplastics Colors
For all microplastics observed, blue colored items showed the highest concentrations, followed by white and black colored items (Figure 8). During spring, the blue color showed a high concentration in both regions followed by black in Beirut and grey in Tyre. The lowest concentration was recorded for the yellow color during spring in both areas. During fall, there was a predominance of white color in Beirut, and transparent one in Tyre (Figure 9).

Polymer Identification
Among the suspected plastics, 94% were identified as plastics. The other 6% were impossible to characterize and they were considered as "unidentified". In the subsample, four types of polymers were specified: Polypropylene (PP), Polyethylene (PE), Polystyrene (PS) and Polyamide (PA) (Figure 10). The identification of colored plastics demanded a precise analysis of the dye (colorant) spectrum in addition to that of the polymers (Supplementary Figure S1). PP was the most abundant polymer representing 48.5% of the samples. Whereas PA was only found in 1.5% of the sample. These polymers constituted the two observed shapes in all samples.

Risk Assessment of Microplastics
The MPs were indexed by their concentration using the Pollution Load Index (PLI) and the results are shown in Figure 11. The PLI values of sites from Tyre region were less than 10, which indicated Frontiers in Environmental Science | www.frontiersin.org July 2021 | Volume 9 | Article 699343 a minor risk category ( Table 2). There was a strong correlation (94%) between total MPs concentration and PLI (Pearson rank correlation, rp 0.94, df 40, p < 0.01). Moreover, significant differences in PLI were also detected according to season, whether taking the sample as a whole or taking each region separately (

Macroplastics Density
Floating macroplastics have been recognized as an ecological problem of global concern (Lebreton et al.   Frontiers in Environmental Science | www.frontiersin.org July 2021 | Volume 9 | Article 699343 reported for the Pelagos Sanctuary located in the Ligurian Sea (average 194 items/km 2 ; range 16-453 items/km 2 ) (Fossi et al., 2017) and the upper range of values found in Central and Western Mediterranean areas (maximum 195 items/km 2 ) (Suaria and Aliani, 2014). According to Table 3 and due to the differences in sampling methodologies and the calculation of concentrations, macroplastics concentrations cannot be considered comparable. The most important factors affecting the distribution of macroplastics in marine waters are the vicinity to marine litter sources and pathways with high uncertainties Frontiers in Environmental Science | www.frontiersin.org July 2021 | Volume 9 | Article 699343 7 regarding their fluxes, as well as, the specific oceanographic features prevailing in each sub-basin (Zeri et al., 2018).

Microplastics Concentration and Polymers Abundances
The MPs concentration recorded for the Lebanese marine environment in this study (12.13 items/m 3 ) was more than 2fold higher than the one recorded in 2018 (4.3 items/m 3 ) (Kazour et al., 2019). This concentration was also higher than those recorded in the surface waters of the South-East Levantine sub-basin (7.68 items/m 3 ) ( Van der Hal et al., 2017) and Iskenderun Bay (7.26 items/m 3 ) in the northeastern Mediterranean (Gündoğdu, 2017). The MPs concentration in Lebanese seawater is also higher than that recorded in the Gulf of Asinara-Sardinia (Western Mediterranean), Ligurian Sea (western Mediterranean), and the Turkish Mediterranean coast ( Table 4). The heterogeneity of the MPs concentrations and distribution levels may be linked to multiple factors, such as, sampling methodologies, digestion and density separation.
In this study, the high concentration of MPs during spring could be related to the seasonality effect. Rainfall during the wet  Frontiers in Environmental Science | www.frontiersin.org July 2021 | Volume 9 | Article 699343 8 season (winter) may explain these differences, as surface runoff can transport plastic litter from inland areas to streams and rivers (Shimizu et al., 2008;Ivar do Sul and Costa, 2013). Moreover, abundances and weights of plastic debris were demonstrated to be significantly higher in the wet season compared to the dry season (Cheung et al., 2016), similar to the obtained pattern in this study.
In addition, the influence of the Mediterranean Sea physical oceanographic features, such as, gyres, eddies, sea surface currents, and fronts might be the cause of heterogeneity in MPs density. In fact, the Eastern Mediterranean presented a high infestation of MPs (specifically small fragments) compared to the Western Mediterranean (Cózar et al., 2015). Grand fluxes of MPs (10-30 kg/km/day) were observed along the Levantine sub-basin coasts from Egypt to Lebanon, Syria, and Turkey (Constantino et al., 2019). The differences between the Eastern and Western Mediterranean basins are due to different circulation patterns: water currents in the Eastern Mediterranean play the role of attractors leading to the accumulation of MPs, while the currents in the Western Mediterranean promote circulation (Mansui et al., 2015). Furthermore, the complex configuration of the current circulation, with the presence of meso-scale and sub-basin scale eddies in the Levantine sub-basin could explain the distribution of MPs along the Lebanese coast (Menna et al., 2012).
Those factors can also influence the polymers distribution in our water samples. Polypropylene and Polyethylene were the most abundant types of plastics found. These two polymers are the most found in the Mediterranean Sea and are two of the greatest subsets of plastics manufactured worldwide (Geyer et al., 2017;De Haan et al., 2019). Their characteristics as flexible, lightweight and low density polymers allow them to float on the surface water (Kazour et al., 2019). They are highly used in packaging (plastic containers, bottles, films and bags), geotextiles and in fishing activities (i.e: fishing nets) (Suaria et al., 2016). The presence of Polyamide, a polymer with a higher denser of marine water (1.14 g/cm 3 ), in our samples might be connected to the turbulent mixing and currents that lead denser polymers to persist on the surface water. Previous studies have also described the presence of polymers with a higher density than seawater Suaria et al., 2016;Kazour et al., 2019). The presence of these polymers on the surface water and their potential degradation might affect the aquatic biota biological performances and induce oxidative stress (Sun et al., 2021). Due to the enormous time consuming process, the Frontiers in Environmental Science | www.frontiersin.org July 2021 | Volume 9 | Article 699343 percentage of identified particles represented 1% of the total sample. Such identification percentage is similar to various microplastics studies that have thousands of suspected MPs usually analyze up to 5% of the total items (Lares et al., 2018;Gies et al., 2018;Magnusson and Norèn, 2014;2-8%;Simon et al., 2018).

Influence of Plastic Colors on Fish Ingestion
In this study, the green and black colors were identified with the dominance of blue., Colored MPs are sometimes confused as preys and are potentially ingested by aquatic biota (Lopes et al., 2020). This ingestion occurs especially during the spring season (i.e. spawning period for most fish). and can affect the growth and  Frontiers in Environmental Science | www.frontiersin.org July 2021 | Volume 9 | Article 699343 development of larval fish (Lönnstedt and Eklöv, 2016). Moreover, microplastics was ingested by three anchovy species were fibers hade a highest concentration than fragments (Guntur et al., 2021). Microplastics are frequently detected in the gastrointestinal tracts of aquatic organisms (Lusher et al., 2017). A number of active and passive pathways have been suggested for fish, including the confusion of microplastic particles with prey, accidental uptake while foraging and transfer through the food chain (Roch et al., 2020). Study of sea birds indicates that plastic particles with algae living on the surface can release infochemicals to entice the sea birds to eat them directly (Cózar et al., 2017). Sea turtles which perceive floating plastics from below, preferentially ingest dark plastic fragments (Santos et al., 2016).
Microplastics that visually resemble the natural prey are actively ingested by several fish species such as Girella laevifrons, Decapterus muroadsi. A laboratory study of palm ruff Seriolella violacea, showed that particles with similar color to artificial food pellets were more often ingested (Mizraji et al., 2017;Ory et al., 2017;Ory et al., 2018).
When the density of food was generally lower, the number of microplastics in the gastrointestinal tract increased with MPs concentration in the water (Roch et al., 2020). Particles smaller than 5 μm can pass the gastrointestinal tract wall and bioaccumulation could arise when uptake exceeds release or when particles are assimilated in tissues or organs and their effects may be significant, especially in long-living species, with implications for food web transfer and fish as food items (Roch et al., 2020). The goldfish (Carassius auratus) ingests polyethylene (PE) microplastics and these processes are affected by color, and MPs shape. However, in the presence of food, more green and black MPs were ingested compared to red, blue, and white MPs while significantly higher amounts of microplastics films were ingested compared to fragments and filaments (Xiong et al., 2019). The blue color was the most frequent in water and shrimps Paratya australiensis (Nan et al., 2020). Only one study performed in Lebanon showed the presence of filaments in the stomachs of European anchovies (Engraulis encrasicolus) constituting 1.6% of the total prey content (Jemaa et al., 2016).

The Pollution Load Index of Microplastics
Based on the PLI results (Figure 9), a total of nine sites in Tyre and Beirut showed values less than 10 (minor pollution), and three sites in Beirut showed values between 10 and 20 (high pollution). The risk within the two regions was determined using PLI zone, and the results indicated a seasonal variability in each region. It was revealed that sites with a high MPs pollution load had a high concentration of MPs. Therefore, increased MPs    Ranjani et al. (2021) and Xu et al. (2018). pollution risk is dependent on high levels of accumulation in the environment and the existence of harmful microplastics polymers . The highest value in Dora can be attributed to the presence of Beirut River that diverse in north of Beirut facing the sampling site. According to the assessment result of PLI, the risk pollution of MPs in surface waters of Beirut and Tyre was minor and not serious. However, the risk level of MPs pollution in Beirut was greater than that of Tyre (mean PLI: 6.33 > 5.02). The risk of MPs pollution is not only driven by the accumulation of MPs caused by hydrological dynamics in the Lebanese marine area, but also through land-based sources like municipal wastes and coastal landfills, which increase the risk of MPs in the marine environment (Swati et al., 2017;Khalil et al., 2018;Kazour et al., 2019). The highest PLI of MPs was recorded at Dora 1 (PSI 17; high risk level), Manara 1 (PSI 12.7; high risk level) and Manara 2 (PSI 19.1; high risk level) in Beirut region.

Value of the pollution load index
The highest values of PLI can be correlated to the land-based source of MPs; for instance, Dora 1 was located near the economic harbor of Beirut (receiving ≈2,300 ships by year), the Beirut River and the Dora landfill. The highest level of PLI in Manara 2 can be related to the presence of a wastewater outflow at its proximity.
According to Plastics Europe (2020), plastic production is forecasted to increase by 4% every year. This current study highlights the presence of macroplastics and microplastics in the sea surface of Lebanon, all of which is related to the increased urbanization along the coast, currents dynamics, and other environmental factors. So far, there are no legislations that aim to manage plastic use in Lebanon, although several organizations/NGOs are working at a local scale to reduce the single use of plastics and encourage recycling. This study can help decision makers take proactive or reactive actions and face  marine litter pollution in the Lebanese coastal and marine environment. However, MPs in the marine environment constitute a worldwide problem that should be thoroughly tackled by all stakeholders and countries.

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

AUTHOR CONTRIBUTIONS
SJ: Conceptualization, methodology, investigation, writing original draft CM: Methodology, investigation review and editing ML: Conceptualization, software, statistic, review and editing AH: Review and editing MG: Software GK: Review and editing MF: Project administration, review and editing.

FUNDING
This work was conducted in the framework of the ELME project "Evaluation of the Lebanese Marine Environment: a multi-disciplinary study" and fully financed by the European Union Ref: ENI/2018/398-295.

ACKNOWLEDGMENTS
First, the authors would like to extend their gratitude to the European Union through funding the ELME project that aims to support the Lebanese marine environment. The authors would like to thank Anna Maria Kalaani and Ridge Antoun (Master Students (M2: BEMA) at the Lebanese University) for their help during sampling and treatment of samples. In addition, the authors' gratitude goes to IUCN ROWA through the Coastal Ecosystem Resilience (CER) project funded by the Royal Norwegian Embassy in Beirut. Finally, the authors would also like to thank the crew of CANA-CNRS R/V, Elie Tarek and Houssein Jaber from the National Center of Marine Sciences for their help during the sampling process.