1. Introduction

The growing population and the industrially based lifestyle has led to an increase in the anthropogenic impact in the world and has brought with it new water challenges including increasing water supply, inacceptable water quality in many areas, increasing number of aquatic ecosystems that are in a danger of collapse and an increasing number of pollutants. The types and sources of pollutants caused by human’s activities, such as the pollution coming from industries, coastal runoff and wastewater effluents, are diverse and have a potential effect and fate on the environment and are threatening, alarmingly, the coastal areas [1] [2] .

The Mediterranean region links 21 countries and attracts millions of tourists, and this has caused the environ-mental degradation in the Mediterranean sea and this degradation is increasing [2] . Human’s activities have transformed the Mediterranean into one of the most polluted seas in the world [3] . An estimated 300 million people live along the shores on the Mediterranean Sea [4] . Adding to that, it occupies an area of 3 million square kilometers [4] . Both factors explain the cause of the high amount of pollution on the coastal area. Moreover, the pollution in the Mediterranean sea is long lasting because it is almost closed [4] .

In the Lebanese Mediterranean Sea, besides being affected by the pollution projected in the Mediterranean, a recent article in Beirut’s English language newspaper, the Lebanese Star noted that pollution caused by both the public and private sectors is reaching an alarming level. Sadly, the Lebanese public is ignoring this problem and it is estimated that at least 200,000 cubic meters of untreated sewage water are being dumped into the sea per day [5] . The situation is so urgent because Lebanon is one of few countries where nearly all sewage goes into the sea [5] . Lebanon’s sea pollution multiplied considerably during the 2006 war, when more than 15,000 of fuel oil spilled into the sea following the bombing of the main power station in Jieh by Israel [6] . This incident has been considered to be the country’s worse environmental nightmare. Beside this, not only sewage and industrial waste were dumped into the sea, but also thousands of tons on untreated solid wastes from a number of dumping sites, on the country’s coastline [5] . The oil spill along with decades of unsustainable coastal urbanization and unregulated dumping of untreated domestic and industrial waste have magnified the pollution in the sea water during the previous years [7] . The need for stopping the evolution of the pollutants nowadays in Lebanon and the need for reducing and limiting the pollution in the Lebanese sea are critical. For this purpose, the Lebanese coastal area needs to be monitored in order to detect and be aware of this evolution. Several studies have been done on detecting the marine pollution by analyzing the sea water chemically and physically in Lebanon [8] and by observing the numbers of the coastal landfills such as Borj Hammoud landfill [9] . In this project, the detection and measurement of pollution in the Lebanese marine environment are realized using Remote sensing techniques. Remote sensors exhibit potentially good pollution detection capabilities because of its ability to observe a wide area at a time [10] .

The overall objective of this study is to monitor and detect the pollution on the whole area of the Lebanese coast using satellite data acquired during 2014 and 2015 to detect: the frequently polluted sites by different types of pollutants; the pollution dimension to know the amount of pollutants dumped into the sea each month; the most polluted month in 2014 and 2015; and the most polluted year between 2014 and 2015.

For this purpose, this study will highlight 2 main parts that will be discussed in detail:

1) Part one, the use of remote sensing and its importance for the detection of marine pollution (2014-2015).

2) Part two, interpret and compare the Land sat images in order to obtain a clear idea about the polluted sites and the polluted months and the changes of pollution between 2014-2015.

2. Materials and Methods

2.1. Study Area

The coastal zone and the Mount Lebanon chain are characterized by their closeness, except in the north and the south of the country. The study area (Figure 1) includes the Lebanese coastline which extends about 230 kilometers in length from the northeast border at Aarida to the southwest border at Naqoura [11] . Based on a study done by the CDR in 2005, the coastal zones represent 8% of the total Lebanese surface area which is approximately 840 km² of the Lebanese territories. The study area includes 15 rivers [11] . These rivers transport the manmade wastes to the sea. The annual mean temperature in the coastal zone varies between 13.5˚C and 27˚C with an average annual rainfall of 600 mm [12] [13] . The coastal zone has a very high population density; it was estimated that more than 1.5 million inhabitants in the whole coastal area, which repre- sents 55 percent of the total population [11] . The major six cities that are highly popu-

lated and located directly on the coast line are from North to South: Tripoli, Jbeil, Jounieh, Beirut, Saida, and Tyre (Figure 1) [13] .

2.2. Landsat Satellite Images

Landsat satellite images are characterized by medium spatial resolution with good capability to discriminate between different seawater aspects: normal versus polluted sea water [13] [14] . Two consecutive years were studied in 2014 and 2015 using Land sat OLI 8 satellite images. This approach aims to identify the pollution sources and types. The analyzed images are taken monthly for each year. This study focuses on monitoring change in pollution sources in space and time as well as the evaluation of the persistence of the identified sources. The goal of making monthly maps is to: compare the pollution changes and propagation monthly and yearly in order to monitor the changes in the types of pollutants; if the pollutants dumped into the sea have increased or decreased; and to detect new polluted sites, permanent polluted sites and to identify the critical, at risk areas where pollution needs to be stopped. Data can be downloaded (at no charge) from Glo Vis, Earth Explorer, or via the Land sat Look Viewer within 24 hours after the acquisition. In this study, Earth Explorer was used for download. All data images selected for analysis require several processing steps. The goal of these processing steps is to increase the accuracy, the clarity and interpretability of the digital data during image analysis. In this study, images were processed using ENVI5.2 and ArcGIS 10.3.1 software. After image processing, identification of marine pollutants was performed usual visual interpretation: color, tone, texture, flow pattern, and temperature differences [13] [14] .

2.3. Image Processing

2.3.1. Multiband Selection

Multiband selection and it is done by interrelating three bands B₄ (which is specialized for the red color), B₂ (specialized for the green color) and B₁ (specialized for the blue color) as one set. A 421 band combination of the coastal waters is also essential in order to identify potential pollutants and other coastal phenomena [10] - [15] .

2.3.2. Image Equalization

Image equalization is used for image enhancement since it clearly discriminates different water colors and tex-tures, indicative of water quality by stretching and maximizing the features of a specific area [10] - [15] .

2.3.3. Thermal Image Processing

Density slice coloring was performed to classify temperatures into ranges assigning each range a different color. First, a selection of a single band (B₁₀) that indicates the thermal temperature, is used, in order to detect temper-ature anomalies. Then, a density slice is applied to elaborate classes that represent the temperatures. The selec-tion temperature intervals depend on the properties that need to be emphasized [10] - [13] and facilitates the identification of temperature anomalies. The plumes indicate the presence of discharge that have a temperature different from the temperature of the neighboring area. This approach compliments the analysis of spectral anomalies using the RGB image and 421 band combination [16] .

3. Results

Processed images from Land sat 8 OLI data showed a variety of anomalous features of the coast of Lebanon between 2014 and 2015. The different types of pollution have been identified directly from satellite images based on the anomalies in 421 band combination along with density sliced thermal bands. The location of pollution is made by detecting colors differentiation (visible contrast among the most distinguishable elements of identification) correlated with associated land features, such as the presence of landfills next to the coast (Borj Hammoud landfill per example); the thermal differentiation of warm water such as warm water coming from power plant.

In Table 1, it was found that Akkar, Nahr el Jaouz, Nahr Ibrahim and Jounieh areas presented some unique anomalies in the maps that needs to be explained:

3.1. Akkar Area

Based on the dimension of the anomaly in Akkar, it is found that during the dry season such as June 2014, May 2015 and August 2015, the anomaly had a bigger width than during the wet season. This might be caused by the lack of water movement which allows the anomaly to be seen clearly. In Akkar area, an anomaly is visible every month. The anomalous could be explained as follow [19] :

Based on a research done by Shaban (2010) titled as “Geothermal Water in Lebanon: An Alternative Energy Source”, it was mentioned that in Akkar region along the Syrian border, there is a large basaltic rock plateau, and several observations for hot or warm water from springs and seepages at this area have been found. Field tests show that the average temperature of the springs ranges between 50˚C - 65˚C over various time periods. Furthermore, the hydrologic phenomenon of hot and warm water is not limited to the terrestrial environment, but extends into the marine environment. It showed the presence of geothermal water in Lebanon in four major geothermal domains especially in Akkar region. The presence of geothermal water in Akkar might explain the appearance of anomalous features on the maps. Another explanation might cause the presence of these anomalies such as the pollution from al Oustowan river that flows into this area [19] .

3.2. Nahr El Jaouz Area

During February and November 2014, the amount of pollution was high compared to other months. These also corresponded with wet months, thus, river flow was greater and may have dispersed the plume. This cannot be said for 2015 when Nahr el Jaouz did not demonstrate an anomaly during the wet season but a large polluted area was visible in May 2015 and bigger than that of 2014. A feasible explanation might be the excessive dumping of wastewater.

3.3. Nahr Ibrahim Area

During the wet season of 2014, Nahr Ibrahim has presented the appearance of anomaly and not during the dry season 2014. This might be due to the intensive water discharge during the wet season and its dryness during the summer. During the summer 2015, Nahr Ibrahim presented some anomalies on the maps. Similarly to Nahr el Jaouz, this might be due to the release of wastewater in the dry river.

3.4. Jounieh Area

An anomaly is visible in the thermal maps, meaning there is a differentiation in the water temperature. This might be caused by the submarine groundwater discharge present in the Jounieh area or there might be presence of warm water coming from ships that is contributing to the temperature differentiation. This anomaly only appeared in few months in 2014, but it was present during all months in 2015. This means that during 2015, Jounieh was more frequently polluted than 2014.

4. Discussion

Based on the results, it is obvious that in 2015 the pollution dimension in most of the sites such as: Tripoli, Chekka, Nahr el Jaouz, Nahr Ibrahim, Jounieh, Nahr Beirut, Ouzai, Nahr el Awali and Jieh and Nahr el Litani was bigger than in 2014. This indicates that in 2015 there was in increase in pollution. This might be due to the lack of monitoring and management of coastal sites leading to an increase of waste dumping along the coast. Moreover, 2015 faced lot of rainy days that led to an increase of river discharges in the sea causing pollution to become more widely dispersed via increase flow at the river mouth.

4.1. Ranking of Anomalies in the Sites Based on Their Frequencies

The following Table 2 presents the ranking of anomalies based on their frequencies, going from the presumed most polluted site to the least polluted site. The determinations are based on the appearance and dimension of anomalies on a monthly and annual basis.

4.2. Influencing Factors

Overall sea current dynamics in the area are directed from the south to the north, which control the spreading regime of pollution and divert it to the north shoreline [20] . This explains why the pollution direction is directed to the north in the maps. During April 2014 and June 2014, the current was directed to the south in northern parts (Figure 2). This might be explained in a number of ways: 1) Due to the diversion of the water current caused by the river flow; 2) springs in the sea water; 3) by pollutant discharge and dispersion; 4) By short term wind.

Further considerations include the dry and wet seasons. Anomalies are most clear during the wet season associated with intensive river discharge creating turbidity. Some focus, therefore, was placed on when an anomaly appears during the dry season. These

cases suggest this might be caused by a certain type of pollution as opposed to turbidity and discharge [21] . Examples are provided in Figures 3 through 8 below.

4.3. Identification of Pollution Anomalies

Figure 3 shows the pollution in red and yellow. Areas of interest are labeled (1), (2), (3), (4).

(1) and (2) are the areas of potential pollution; (3) marks possible oil slicks; (4) is the passage of a ship.

In Figure 4 it can be observed that the areas (1) and (2) of potential pollution extracted from Figure 3 continue to exist however the ship and oil slicks are no longer present, possibly explaining the color differences.

A Principal Component Analysis (PCA) was performed to enhance all uncorrelated bands and maximize any differences in ocean waters. The purpose is to emphasize potential pollutants. The same areas of interest are la-beled and clearly identifiable. This map gave the same information as the 421 band combination. The color contrasts are more defined and no new sources of pollution are apparent.

Again it is clear that the areas of potential pollution from Figure 5 (1, 2) still exist and although they are represented in different shades between the 2 dates in Figure 6.

It is feasible that different pollutants are being dumped between the two locations or that ocean circulation is causing them to be redistributed.

The spectral curves of all Landsat bands were examined for a profile straddling the area of pollution in site 1 in order to determine areas with strong absorption. The premise was that the band with the most absorption could be divided by the 2 adjacent bands using the equation:

(1)

to emphasize unique properties pertaining to pollution. It is possible that a zone with concentrated pollutants is emphasized by the light blue plume (1).

The plots in Figure 7(b) were examined for absorption features for the polluted areas. B₁ and B₂ was identified as having absorption centered over the feature to be detected (red lines), thus is used to create the simple ratio image Figure 7(a).

This time the blue plume in Figure 8 has dispersed and there is no longer a spike for these pollutants in site 1. This would be a good site to monitor on an operational basis for periodic dumping of polluting material. For this date, signs of pollution are reduced.

5. Conclusions

In order to monitor and detect the marine pollution on the Lebanese coast during the year 2014 and 2015, satellite images facilitated observation of the entire coastal region and provided valuable information concerning the detection of anomalies in this wide area of study. This detailed reconnaissance mapping helps identify those areas susceptible to pollution requiring continued monitoring. These data can be used on an operational basis enabling changes to be studied every 16 days. Once data are obtained, permanent hotspots and anomalies on the whole Lebanese coast can be depicted. During the dry season, there is low movement of water causing a least extension of the anoma lies. In contrary, during the wet season, rivers had an intense flow into the sea which caused an intense water movement and wide extension of anomalies on the coast. Moreover, it could be determined that the permanent polluted coastal sites in 2014 and 2015 are Tripoli, Kalamoun, Chekka, Batroun and Amchit, Jbeil, Jounieh, Nahr Beirut and Ouzai. The other sites are not considered permanent because the pollution appeared intermittently.

Data that showed most of anomalies were during April and November 2014 and April 2015. Data also showed that the least appearance of anomalies were during July 2014 and July 2015 which also correlates with the dry season and the low extension of anomalies. Spatially, the length and width of each anomaly at each site was documented in order to obtain a clear idea about the dimension of each anomaly and its importance based on its size. It is clear that pollution anomalies have increased in size over time from 2014 to 2015.

In sum, utilizing satellite data to implement a national monitoring system on water quality for the Lebanese coastal environment is critical. The goal being to reduce the number of pollutants dumped into the marine environment.

Acknowledgements

Professor Nadine Nassif would like to extend her appreciation to the Fulbright scholarship obtained at Northeastern University, College of Professional Studies. The collaboration between the Lebanese University and Northeastern University is appreciated and ongoing.

*Detection of marine pollution by using remote sensing all along the Lebanese coast during the years of 2014-2015 qualitatively.