On the Retrieval of the Water Quality Parameters from Sentinel-3/2 and Landsat-8 OLI in the Nile Delta’s Coastal and Inland Waters

: Reduced water quality due to the eutrophication process causes large economic losses worldwide. Multi-source remotely-sensed water quality monitoring can help provide effective water resource management. The research evaluates the retrieval of the water quality parameters: chlorophyll- a (Chl- a ), total suspended matter (TSM), and chromophoric dissolved organic matter (CDOM), over optically different water types. Cross-sensor performance analysis of three satellite data sources: Sentinel-3 Ocean Land Color Imager (OLCI), Euglenophycaea, and Bacillariophycaea) and the biological contents (carbohydrates, lipids, and proteins), arranged in decreasing order. The largest normalized RMSE estimates marked the downstream areas where the inﬂow of polluted water persistently brings nutrient loads of nitrogen and phosphorous compounds as well as substantial amounts of detrital particles and sediments discharged from the agricultural and industrial drains and the land use changes related to agricultural practices, resulting in the increase of water turbidity giving rise to inaccurate Chl- a estimates.


Introduction
Coastal and inland water conditions, in particular, the chlorophyll-a (Chl-a) content and its accurate estimation controlling parameters, such as the total suspended matter no rigorous theoretical basis to support them, which indicates that the models here could not be directly used for other waters [62].
Further, the fast microalgae growth variation hinders accurate estimation of their Chl-a and the biomass content due to the dynamic light regime and intensity, which enhance and fluctuate the photosynthetic efficiency inducing the microalgae growth rate and hence productivity (Chl-a and biomass concentration), which can be 4 to 4.5-folds higher under continuous illumination conditions [63,64].
The Nile delta coast of Egypt is an important region where periodic monitoring of water quality is of interest. A systematic study with a long-term perspective appraising the synoptic view and revisit time of remote sensing, and its inherent ocean color products, is needed to characterize the spatiotemporal variability of the water conditions affecting its quality related to the algal abundance, and the harmful algal blooms.
Therefore, the overarching focus of this research is to evaluate the effectiveness of the C2RCC in optically varying water types (marine coastal water and inland lakes) in the Nile Delta by analyzing the performance of the Sentinel-3 OLCI, Sentinel-2 MSI, and Landsat-8 OLI acquired during a 45 min overpass on the Nile Delta coast on 22 March 2020 for retrieving the Chl-a content and the factors affecting their retrieval, i.e., the spectral scattering of the TSM, and the absorbing CDOM contents, through the inter-comparison of the retrieved estimates with the standard reference CMEMS, S3 OLCI, and MODIS data products. The inter-comparison is appraised to examine its validity as an alternative approach to overcome the lack of ground-truth data to derive average and range levels estimates in the coastal water. Further objectives pursued were to validate the products from the MSI data for the inland Burullus Lake water acquired in the 18-20 September 2020 period, as the only available cloud-free sensor data that is nearly contemporaneous with in situ Chl-a measurements carried out during the 22-25 September 2020 period, and to demarcate common factors hindering accurate retrieval of the water quality parameters in the complex inland lake.
In this research, the relevant data sources used and the preprocessing routines for Sentinel-3 OLCI, Sentinel-2 MSI, and Landsat-8 OLI imagery are firstly described. Next, the C2RCC used for atmospheric correction is explained. Finally, the results in overlapping region of the studied data footprints are discussed in light of the validity of the retrieved water parameters and the estimation error sources both for the marine and the inland waters based on inter-sensor comparison of the 45 min sensed time data and validation in the local inland lake.

Study Area
The region of interest is located on the Mediterranean coast of the Nile Delta and its fringes with a special focus on the Burullus Lake environment representing the inland lake waters (Figure 1). Higher values of biomass, and consequently of Chl-a, may be found in areas influenced by watercourses runoff into the study area, which was selected because of its complex water dynamics and its dense anthropogenic pressure. Burullus Lake is the second largest of the northern coastal lakes along the Mediterranean coast of Egypt, covering about 420 km 2 commonly used for fishing and recreation activities. The lake is covered by floating vegetation, isolated sand bars, urban areas, and agriculture lands.
The physical and biological characteristics of the lake has been recently documented by Masoud et al. [65]. The lake water showed a depth average of 130.5 cm with a range level of 0.90-2 m, alkaline pH (7.78-8.92), and average salinity of 4.33 gm/L with a range of 0.3-17.19 gm/l. The phytoplankton displayed a high level of variation in the number of taxa (145 species) and abundance with maximum counts of 1005 × 10 3 units.L −1 to a minimum of 340 × 10 3 units.L −1 . Six groups of freshwater forms dominated by the Chlorophyceae, Bacillariophyceae, Cyanophyceae, Dinophyceae, Euglenophyceae, and Silicoflagillate, were arranged in decreasing order of abundance. of 0.3-17.19 gm/l. The phytoplankton displayed a high level of variation in the number of taxa (145 species) and abundance with maximum counts of 1005 × 10 3 units.L −1 to a minimum of 340 × 10 3 units.L −1 . Six groups of freshwater forms dominated by the Chlorophyceae, Bacillariophyceae, Cyanophyceae, Dinophyceae, Euglenophyceae, and Silicoflagillate, were arranged in decreasing order of abundance.

Satellite Data and Image Processing
Different data sets of S3A/OLCI, S2A/MSI, and L8/OLI imagery were searched and allocated for atmospheric correction assessment ( Figure 2). The imagery file naming, processing levels and sensing time are shown on Table 1. The sensing time period was 45 min for the three data sets, starting from 7:51 to 8:36 AM acquired on 22 March 2020. Image processing was carried out for the Landsat-8 Operational Land Imager (OLI) and for two Sentinel sensors developed by the European Space Agency (ESA) as part of the Copernicus Earth Observation Programme for the monitoring of natural water bio-optical properties: Sentinel-2 Multi-Spectral Instrument (S2-MSI) and Sentinel-3 Ocean and Land Color Instrument (S3-OLCI) [66,67]. The sensors slightly vary in the bandwidths and their spectral locations, as indicated in Table 2.

Satellite Data and Image Processing
Different data sets of S3A/OLCI, S2A/MSI, and L8/OLI imagery were searched and allocated for atmospheric correction assessment ( Figure 2). The imagery file naming, processing levels and sensing time are shown on Table 1. The sensing time period was 45 min for the three data sets, starting from 7:51 to 8:36 AM acquired on 22 March 2020. Image processing was carried out for the Landsat-8 Operational Land Imager (OLI) and for two Sentinel sensors developed by the European Space Agency (ESA) as part of the Copernicus Earth Observation Programme for the monitoring of natural water bio-optical properties: Sentinel-2 Multi-Spectral Instrument (S2-MSI) and Sentinel-3 Ocean and Land Color Instrument (S3-OLCI) [66,67]. The sensors slightly vary in the bandwidths and their spectral locations, as indicated in Table 2.    The Landsat-8 signal-to-noise ratio (SNR) has been shown to increase significantly with narrower bandwidth compared to Landsat's previous missions in the red and near-infrared (NIR), among others, and attains radiation resolution, increasing to 16 bits with a 16-day repeat cycle. All advances have attributed to improved pigment discrimination ability, which have proven useful for estimating concentrations of Chl-a in water bodies [53,68].
The Sentinel-2 Multi-Spectral Instrument (S2MSI) attains a high 12-bit radiometric resolution, wide spatial and temporal resolution; the location and bandwidth are narrower compared to the OLI data and has a wide field of view (FoV) that generates spectral images composed of 13 bands at high SNR and varying wavelengths ( Table 2). The use of S2-MSI images has been common in monitoring lakes and rivers and in developing predictive models for Chl-a (e.g., [69][70][71][72][73][74]). Cloud-free imageries were only available from the S2MSI acquisitioned within the same week of in situ sampling from 55 stations, on 18 September 2020 for the western part, and on 20 September for the eastern part of the inland Burullus Lake, that was validated against ground-truth in situ data.
The Sentinel 3 OLCI has two sets of data: unprocessed earth full resolution (EFR) data and the processed for water full resolution products (WFR). The OLCI large number of bands improve the atmospheric correction over the optically-complex waters [75], O 2 gas absorption correction, and the water constituents' retrieval. The ready-to-use OLCI Chl-a data have two products derived by applying the Neural Network (NN) Inverse Radiative Transfer Model (e.g., [76][77][78]) and the Ocean Color 4 for MERIS (OC4Me) of the chlorophyll-a (Chl-a), a semi-analytical model that is based on a polynomial algorithm that uses a maximum band ratio approach of reflectances at 443, 490, and 510 nm over the 560 nm [79] based on the analysis of AOPs measured in situ over the past decades in various oceanic regions (e.g., [79,80]).
The surface Chl-a of the Copernicus Marine Environment Monitoring Service (CMEMS), S3 WFR, and the MODIS are used as reference standard for evaluation. The Chl-a concentration of the coastal water estimated on 22 March 2020 sampled daily at a spatial resolution of 0.042 degrees is analyzed from the CMEMS data available at (https://resources.marine. copernicus.eu/product-detail/MEDSEA_ANALYSISFORECAST_BGC_006_014, last accessed 1 January 2022). The remotely-sensed water geochemical parameters of the pH, alkalinity, dissolved inorganic carbon, zooplankton biomass, chlorophyll-a, phytoplankton biomass, nitrate, phosphate, silicate, and ammonia is used for correlation and regression analysis. The MODIS Aqua measurements Level 3 (4 km), downloaded from the ocean color website (https://oceancolor.gsfc.nasa.gov/) were acquired on 22 March 2020 from 11:10 to 11:14 a.m.
A flowchart explaining the adopted methods and associated steps is shown on Figure 3. A flowchart explaining the adopted methods and associated steps is shown on Figure  3. All sensor data were first masked for clouds and cloud shadows, and subsequently atmospherically corrected using the C2RCC atmospheric correction scheme. The algorithm runs on the visible-NIR bands ( Table 2) commonly used for water color retrieval, with 5 bands for the OLI data, the MSI has 9 bands out of the sensor 13 bands at 10-60 m grid spacing, while OLCI has the 21 spectral bands in this region at a resolution of 300 m [81]. The satellite imageries were geometrically corrected to the UTM Zone 36, resampled, to 10 m resolution for MSI, 30 m for OLI, and 300 m for OLCI EFR; all were then subset to the study area. The land and cloud covers were masked using the metadata associated with the imageries. The atmospheric correction was then applied by appraising the use of the C2RCC algorithm implemented in the Sentinel Application Platform (SNAP) version 5.0 with the Sentinel-3 Toolbox Kit Module (S3TBX) version 5.0.1 for the three sensor data sets. Simply, C2RCC is based on water radiative transfer modeling and the vector successive order of scattering (SOS) atmospheric model with aerosol optical properties derived All sensor data were first masked for clouds and cloud shadows, and subsequently atmospherically corrected using the C2RCC atmospheric correction scheme. The algorithm runs on the visible-NIR bands ( Table 2) commonly used for water color retrieval, with 5 bands for the OLI data, the MSI has 9 bands out of the sensor 13 bands at 10-60 m grid spacing, while OLCI has the 21 spectral bands in this region at a resolution of 300 m [81]. The satellite imageries were geometrically corrected to the UTM Zone 36, resampled, to 10 m resolution for MSI, 30 m for OLI, and 300 m for OLCI EFR; all were then subset to the study area. The land and cloud covers were masked using the metadata associated with the imageries. The atmospheric correction was then applied by appraising the use of the C2RCC algorithm implemented in the Sentinel Application Platform (SNAP) version 5.0 with the Sentinel-3 Toolbox Kit Module (S3TBX) version 5.0.1 for the three sensor data sets. Simply, C2RCC is based on water radiative transfer modeling and the vector successive order of scattering (SOS) atmospheric model with aerosol optical properties derived from NASA AERONET-OC measurements [82,83]. The artificial neural network technology [84,85] is used where for every sensor data, a subset of bands is set as neural nets; the main net derives water-leaving reflectance 'ρ w (λ)' and the top-of-atmosphere (TOA) ρ toa (λ) radiances after atmospheric correction with an elaborate bio-optical oceanatmosphere model using a large database of water-leaving reflectances and relating to the optical properties from various coastal areas around the world. Subsequently, the trained neural nets parametrize the inverse relationship between inherent optical properties (IOPs) and reflectances (water-leaving ρ w (λ) from the ρ toa (λ)), allowing the retrieval of certain water quality parameters, such as chlorophyll-a, TSM, and the CDOM [86].
For all combination of the IOP components, based on a large world database of radiative transfer simulations in natural waters as "truth", and the trained IOP inversion net as "estimate", the difference between truth and estimate gives the uncertainty per IOP per pixel in the image [82], defined here as the associated error of the water color constituent. C2RCC proved applicable to a wide range of satellite sensors, such as SeaWiFS, MERIS, MODIS, OLI, MSI, and OLCI [85], and efficient in many world case studies for the retrieval of Ocean Color estimates; detailed information on processing is given in Brockmann et al. [85].
For the C2RCC implemented in SNAP, the local relationship between IOPs and concentrations of optical substances in the water that are regionally variable can be adjusted by adapting the ancillary parameters to local waterbody specific inherent optical properties (sIOPs), such as salinity, temperature, ozone, air pressure as well as the specific IOPs, namely the Chl-specific absorption coefficient and the specific scatter of TSM at 442 nm. For this study, the coastal water salinity was set to 39 PSU and the water temperature was set to 21.5 • C as common values for the study area in the Mediterranean Sea stated in several studies, which is different from the default setting (e.g., [87]) and the sea surface data of the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) for temperature; [88] for salinity). EUMETSAT data is available at https://www.eumetsat.int/S3b-sst-processing-baseline, last accessed 1 January 2022. The values of temperature and salinity of the studied period and the examined area investigated using the Copernicus products sampled hourly at a spatial resolution of 0.042 degrees (https://resources.marine.copernicus.eu/?option=com_csw&view=details& product_id=MEDSEA_ANALYSISFORECAST_PHY_006_013, last accessed 1 January 2022) were comparable to the values set for the analysis. The average values of temperature and salinity of Burullus Lake were set from in situ measurements [65], sampled in the period 22-25 March 2020 as 28.8 • C and 4.32 PSU, respectively.
After the AC and prior to pixel extraction, MSI and OLI pixels were resampled to 300 m and registered to match the OLCI georeferenced image to allow inter-comparison.

Inter-Sensor Comparison and Validation in Burullus Lake
Due to a lack of ground-truth data for the March 2020 sensor data, metrics of the Pearson correlation coefficients (r), least-square coefficient of determination (R 2 ) of multivariate regression analyses, along with the mean absolute error (MAE), and the root mean square error (RMSE) were employed and evaluated for sampling location points among the three sensor-based empirical models and compared with the standard measurements of the CMEMS, OLCI WFR, and MODIS data in order to examine the accuracy and cross-method consistency of water quality products.
The MAE and the RMSE metrics indicate the relative error of estimates and account for the proportionality of the errors with the concentration of the constituents.
Equations (1) and (2) are used for estimating the MAE and RMSE, respectively.
where N is the number of sample points, y i is the estimated water color parameter from the three data sets, and x i is the value of water color value of the standard products. Sampling of 275 locations in the coastal water and 30 points in the inland Burullus Lake were extracted from the overlapping region of the data footprints through on-screen digitizing selected to cover varying optical contents, water depths, and varying distances from the downstream areas of the inland watercourses. For validation in the inland lake water, in situ data from 55 sampling points across Burullus Lake, collected during 22-25 September 2020, were used [65] for evaluating the water products retrieved from the S2 MSI image acquired on 18-20 September 2020, where r, R 2 , RMSE have been quantified.

Distribution of Retrieved Water Color Parameters
Retrieved water color parameters of Chl-a, TSM, and CDOM along with their associated errors are shown for the S2MSI (Figure 4), selected for display to attain the finest resolution of 10 m among all data. The distribution of the retrieved parameters for the S3 OCLI (WFR and EFR), OLI, and MODIS are shown in the Supplementary file to this article. The distribution of water content was uneven with the maximum concentrations marked by the inland lake water and the downstream areas located in front of the outlet of the northern lakes into the coastal water, close to the shoreline. With focus on the large scale of 10 m S2MSI data (Figure 4), as was clear on all sensor data, bloom areas of Chl-a, TSM, and CDOM characterize the Rosetta promontory and Boughaz El-Burullus which discharge large amounts of water of intense nutrient loads into the sea, in particular, nitrogen and phosphorous compounds of vital importance to phytoplankton growth and productivity. A recent study by [65] indicated that Burullus Lake is demarcated with wide Chl-a blooms (av. 83.4 mg.m −3 and range levels of 53.76-129 mg.m −3 ) intensified by the inflow from agricultural drains in the east, and a fresh Nile water canal in the west. Chl-a concentrations of 2.144 mg m −3 indicated a phytoplankton bloom in the Southern Ocean [89].
Based on this value, marked blooms can be easily mapped and their controlling factors can be spatially better understood. These blooms, as indicated in our recent study, can be used in efficient monitoring of the harmful algal blooms and their biological contents of carbohydrates, lipids, and protein can be utilized for the phytoplankton bioenergy potential mapping practices.

Inter-Sensor Data Comparison
Inter-comparison of retrieved parameters applying C2RCC over optically different water types proved successful in the open coastal water, while revealing problems for the The Chl-a content showed a gradual decrease with a depth pattern parallel to the coastline, where local highs are located in front of the Rosetta and Damietta promontories and Boughaz El-Burullus ( Figure 4). The patterns of concentrations in all investigated sensors and reference data were congruent and showed similar trends in the coastal water and varied widely in the inland lake and fish farm waters.
Assessing the retrieved parameters for the 275 sampling points in the coastal water clarified that the Chl-a content averaged 3.14 mg m −3 in the range level of 0.39-4.81 mg m −3 for all data. TSM clarified a mean of 7.66 g m −3 ranging from 6.32 g m −3 to 10.18 g m −3 . CDOM averaged 0.18 m −1 with a range level of 0.13-0.30 m −1 (Table 3). With focus on the large scale of 10 m S2MSI data (Figure 4), as was clear on all sensor data, bloom areas of Chl-a, TSM, and CDOM characterize the Rosetta promontory and Boughaz El-Burullus which discharge large amounts of water of intense nutrient loads into the sea, in particular, nitrogen and phosphorous compounds of vital importance to phytoplankton growth and productivity. A recent study by [65] indicated that Burullus Lake is demarcated with wide Chl-a blooms (av. 83.4 mg.m −3 and range levels of 53.76-129 mg.m −3 ) intensified by the inflow from agricultural drains in the east, and a fresh Nile water canal in the west. Chl-a concentrations of 2.144 mg m −3 indicated a phytoplankton bloom in the Southern Ocean [89].
Based on this value, marked blooms can be easily mapped and their controlling factors can be spatially better understood. These blooms, as indicated in our recent study, can be used in efficient monitoring of the harmful algal blooms and their biological contents of carbohydrates, lipids, and protein can be utilized for the phytoplankton bioenergy potential mapping practices.

Inter-Sensor Data Comparison
Inter-comparison of retrieved parameters applying C2RCC over optically different water types proved successful in the open coastal water, while revealing problems for the inland lake water with estimating accurate variables induced mostly by the adjacency effect and the optical dominance of the CDOM absorption, which cause lower reflectances.

Correlation Analysis
Correlation analysis for 275 sampling points visually selected to cover varying depths along the coastal waters clarified that the Chl-a contents attained the most inferior coefficients for the reference CMEMS and the estimated OLI contents among all data (Table 4).  MSI proved second in rank where it showed strong positive correlation against WFR (r = 0.84) and OLI (r = 0.87). MSI proved superior over all sensors for estimating the CDOM content against WFR (NN) (r = 0.93), EFR (r = 0.80), and OLI (r = 0.80). Correlating the Chl-a contents among estimated and reference data (Table 5 and Figure 5) against the coastal water geochemistry (CMEMS) clarified a decreasing trend with depth, with strongest correlation marked by the WFR (NN) and MSI data (r => 0.80). MSI estimates are by far the strongest significant positive coefficients against all parameters, following in rank that against the reference CMEMS Chl-a contents (Table 5).

Regression Analysis
The satellite sensors also clarified overall good performance in the retrieval of Chl-a, TSM, and CDOM in the visible and NIR bands, which proved optimal in deriving the water color parameters. Regression analysis of the Chl-a contents of the EFR, MSI, and OLI against the reference standard CMEMS, WFR, and MODIS clarified that the MSI estimates are advantageous over all data followed by EFR, where OLI came last in rank ( Figure 6). MSI showed best fits against WFR (NN) (R 2 = 0.94), MODIS (R 2 = 0.77), WFR (OC4Me) (R 2 = 0.71), and CMEMS (R 2 = 0.43), in decreasing order.

Regression Analysis
The satellite sensors also clarified overall good performance in the retrieval of Chl-a, TSM, and CDOM in the visible and NIR bands, which proved optimal in deriving the water color parameters. Regression analysis of the Chl-a contents of the EFR, MSI, and OLI against the reference standard CMEMS, WFR, and MODIS clarified that the MSI estimates are advantageous over all data followed by EFR, where OLI came last in rank ( Figure 6). MSI showed best fits against WFR (NN) (R 2 = 0.94), MODIS (R 2 = 0.77), WFR (OC4Me) (R 2 = 0.71), and CMEMS (R 2 = 0.43), in decreasing order. Regression analysis of the TSM and CDOM (Table 6) clarified that the TSM estimates of MSI proved advantageous over all sensors with best fits (R 2 ) in the range level of 0.82-0.91, being maximum against OLI (R 2 = 0.91), followed by WFR (R 2 = 0.86), and EFR (R 2 = 0.82). CDOM estimates best fitted with WFR (R 2 = 0.73), OLI (R 2 = 0.70), and EFR (R 2 = 0.42), in a decreasing order. The R 2 proved low for the retrieved parameters in Lake Burullus compared to that in the coastal marine water (Table 7).  Regression analysis of the TSM and CDOM (Table 6) clarified that the TSM estimates of MSI proved advantageous over all sensors with best fits (R 2 ) in the range level of 0.82-0.91, being maximum against OLI (R 2 = 0.91), followed by WFR (R 2 = 0.86), and EFR (R 2 = 0.82). CDOM estimates best fitted with WFR (R 2 = 0.73), OLI (R 2 = 0.70), and EFR (R 2 = 0.42), in a decreasing order. The R 2 proved low for the retrieved parameters in Lake Burullus compared to that in the coastal marine water (Table 7).

Retrieval Performance Analysis Metrics
Over all the investigated sensors with local water salinity and temperature inputs in the coastal water, MSI ranked first achieving the smallest MAE and RMSE for the Chl-a contents, while the EFR proved superior with lowest MAE and RMSE for TSM and CDOM estimates ( Table 8).

Burullus Lake Water Quality Parametrs
Spatial distribution of the water quality parameters for the Chl-a, TSM, and CDOM content in Burullus Lake is shown on Figure 7.
The in-situ data indicated that the Burullus Lake water is hypertrophic with a wide variation of Chl-a blooms (av. 83 (Table 9). This confirmed a clear underestimation in the level range of 0.1-59% and average of 35.35% of the Chl-a, which becomes large at larger in situ Chl-a concentrations. Correlation analysis shown on Table 10 indicated that the MSI-derived Chl-a contents are much affected by the TSM (r = 0.71) and CDOM (r = 0.68) contents. Fair to weak, either positive or negative coefficients are clarified between the estimated water quality parameters against the in situ measured lake water variables (e.g., TDS, T • C, algal composition and abundance, and biological contents of Chl-a, proteins, carbohydrates, and lipids).

Burullus Lake Water Quality Parametrs
Spatial distribution of the water quality parameters for the Chl-a, TSM, and CDOM content in Burullus Lake is shown on Figure 7.  (Table 9). This confirmed a clear underestimation in the level range of 0.1-59% and average of 35.35% of the Chl-a, which becomes large at larger in situ Chl-a concentrations. Correlation analysis shown on Table 10 indicated that the MSI-derived Chl-a contents are much affected by the TSM (r = 0.71) and CDOM (r = 0.68) contents. Fair to weak, either positive or negative coefficients are clarified between the estimated water quality parameters against the in situ measured lake water variables (e.g., TDS, T°C, algal composition and abundance, and biological contents of Chl-a, proteins, carbohydrates, and lipids). Table 9. Summary statistics of the MSI Chl-a, TSM, and CDOM contents and the in-situ measured chlorophyll-a in Burullus Lake.   Table 9. Summary statistics of the MSI Chl-a, TSM, and CDOM contents and the in-situ measured chlorophyll-a in Burullus Lake.

Validation and RMSE Characterization
The RMSE between measured and estimated Chl-a content was then estimated, normalized between minimum and maximum values, and regressed against some lake water variables to characterize possible factors enhancing the disagreement magnitudes at the 55 sampling stations taking the spatial distribution of the in situ measured Chl-a contents as reference. The normalized RMSE was then mapped and lake water factors inducing retrieval of water quality parameters were identified and characterized.
The largest Chl-a contents associated with largest normalized RMSE estimates marked the downstream areas where the inflow of polluted water persistently brings nutrients loads of nitrogen and phosphorous compounds as well as substantial amounts of detrital particles and sediments flowing from the agricultural and industrial drains into the lake shallow zones from the south and east causing re-suspension of bottom sediments in the water column, resulting in the increase of water turbidity, and giving rise to inaccurate Chl-a estimates (Figure 9). These areas attained the largest contents of algal composition and abundance of total algal counts dominated by Chlorophyacea, Euglenophycaea, and Bacillariophycaea, which clarify largest biological contents of carbohydrates, lipids, and proteins. The intense environmental pollution marked by intense eutrophication dominates large error-prone areas mostly linked to climatic, hydrologic fluctuation of the drain water flows most affected by wind direction and turbulent water intermixing regimes, and associated nutrient load discharged from land use changes to fish farms and the related agricultural practices south of the lake. The dynamicity in shallow (<2 m depth) water spatio-temporal turbidity and associated pollution loads hindered accurate estimates of the lake water quality parameters from the S2 MSI imagery applying the C2RCC algorithm. The largest Chl-a contents associated with largest normalized RMSE estimates marked the downstream areas where the inflow of polluted water persistently brings nutrients loads of nitrogen and phosphorous compounds as well as substantial amounts of detrital particles and sediments flowing from the agricultural and industrial drains into the lake shallow zones from the south and east causing re-suspension of bottom sediments in the water column, resulting in the increase of water turbidity, and giving rise to inaccurate Chl-a estimates ( Figure 9). These areas attained the largest contents of algal composition and abundance of total algal counts dominated by Chlorophyacea, Euglenophycaea, and Bacillariophycaea, which clarify largest biological contents of carbohydrates, lipids, and proteins. The intense environmental pollution marked by intense eutrophication dominates large error-prone areas mostly linked to climatic, hydrologic fluctuation of the drain water flows most affected by wind direction and turbulent water intermixing regimes, and associated nutrient load discharged from land use changes to fish farms and the related agricultural practices south of the lake. The dynamicity in shallow (<2 m depth) water spatio-temporal turbidity and associated pollution loads hindered accurate estimates of the lake water quality parameters from the S2 MSI imagery applying the C2RCC algorithm.

Discussion
The used multi-sensor data proved promising to a large extent of success for the retrieval assessment of the complex case 2 coastal water quality parameters, including Chla, TSM, and CDOM. Mean of 3.14 mg m −3 , 7.66 g m −3 , and 0.18 m −1 with range levels of 0.39-4.81 mg m −3 , 6.32-10.18 g m −3 , and 0.13-0.30 m −1 were estimated for the Chl-a, TSM, and CDOM, respectively. In situ Chl-a data available from recent local reports and litera- Figure 9. Spatial distribution of the in situ measured Chl-a content and the associated normalized RMSE against the MSI derived Chl-a content.

Discussion
The used multi-sensor data proved promising to a large extent of success for the retrieval assessment of the complex case 2 coastal water quality parameters, including Chl-a, TSM, and CDOM. Mean of 3.14 mg m −3 , 7.66 g m −3 , and 0.18 m −1 with range levels of 0.39-4.81 mg m −3 , 6.32-10.18 g m −3 , and 0.13-0.30 m −1 were estimated for the Chl-a, TSM, and CDOM, respectively. In situ Chl-a data available from recent local reports and literature clarified good agreement with the mean and the range level estimated in the present research. The Egyptian Environmental Affairs Agency (EEAA) report in 2019 for the Mediterranean coast of Egypt (https://www.eeaa.gov.eg/en-us/mediacenter/reports/ projectstudies/eimp.aspx, last accessed on 28 January 2022) clarified a yearly Chl-a average of about 2 mg m −3 and range level of 0.42-5.33 mg m −3 . In addition, in situ monthly Chl-a from March 2019 to February 2020 at the beaches in Alexandria [90], located on the western part of the study area's coastal water clarified an average of 2.47 with a range level of 1.37-3.19 mg m −3 . Further, in situ Chl-a measurements for 692 points distributed in the Mediterranean Sea showed a range level of 0.01-4 mg m −3 , which were comparable to Chl-a estimated from the S3 OLCI and the CMEMS data acquired between the years 2016-2018 [91].
Sensor differences in spatial and spectral samplings, induced and disclosed how well the sensors are suited for the seamless generation of the water quality parameters. The C2RCC provided the pixel-wise certainty level of the estimated parameter concentration and also the associated uncertainty. The certainty maps provide a support to understand the challenges in ocean color monitoring by the studied sensors where these certainty maps can be used as a mask, to disregard areas with relative high uncertainties, and keep the estimates, where the computed statistical measures are valid. C2RCC was originally developed for marine waters using bands at the blue and green regions of the electromagnetic spectrum [92].
In the coastal water, close and parallel to the shoreline, C2RCC proved promising where at shallow depths, nutrient loads (e.g., nitrogen and phosphorous compounds of vital importance to the phytoplankton growth) flow into the sea from the Nile branches and outlets of the coastal inland lakes, the increasing trends of which can indicate the coastal aquatic ecosystems' eutrophication.
C2RCC failed to give accurate results, with an underestimation of the Chl-a concentrations in an irregular shallow (<2 m depth) hypertrophic lake with toxic phytoplankton blooms due to high CDOM, which strongly absorbs in the blue region, and masks the Chl-a absorption [93,94]. This is affected by turbidity loaded with TSM and CDOM in the downstream areas of drains flowing into the lake as well as by the coastline adjacency posing invalid mixed pixels caused by molecular and aerosol scattering where the target pixels are affected by an increase of signal due to the proximity of a neighboring bright land surface, and mostly induced by the shallow bottom sediment reflectance. Similar results and conclusions were drawn from eutrophic lakes in Lithuania [95] and in CDOMdominated waters in Estonia and Sweden [96], and in similar turbid Italian lakes [97], which consolidate the validity of the C2RCC in varied regions. In summary, TSM is retrieved well from all data at a content below 10 gm −3 , while Chl-a seems to be retrieved rather well at concentrations below 5 µg m −3 . The clear underestimation of Chl-a at larger concentrations in Burullus Lake has also been observed using MSI and OLCI data above optically-complex worldwide lake waters (e.g., [98,99]. The ever-enhancing spatial, spectral, and radiometric resolutions, and the signal-tonoise ratio of the Landsat-8 and Sentinel-2/3 over existing ocean color capable missions, such as the MODIS, validated their joint evaluation in studying heterogeneous coastal and inland waters, where the typical 1 km coarse resolution of existing global sensors cannot resolve the fine spatial and the spectral dynamics of the water constituents. This further consolidates their environmental monitoring capabilities for aquatic ecosystems, especially in coastal environments. Further, the incorporation of the new red edge of increased spectral reflectances, results from the presence of partly submersed vegetation (e.g., [100,101]) or algal bloom surface expressions (e.g., [102,103]), with various edge spectral bands at varying spatial resolutions; three for the MSI sensor: B5 (705 nm) of 10 m, B6 (740 nm) of 20 m, and B7 (783 nm) of 60 m; and eight bands for the OLCI sensor (673.5-778.75 nm), compared to the OLI sensor, improve the accuracy of estimating various water bio-optical variables [104,105].
The spectral enhancement reduces artifacts and biases in Chl-a retrieval due to residual glint, stray light, atmospheric correction errors, and white or spectrally-linear bias errors. Moreover, with the potential for higher temporal resolution of about 2-3 days of the three sensors, it is now possible to integrate and fuse the products from these satellites' higher temporal monitoring of aquatic systems. Currently in-orbit Sentinel-2A/B twin satellites provide a five-day revisit time at the equator and even less time at higher latitudes, enabling continuous synoptic monitoring for water quality in small areas where the Sentinel-3 OLCI 300 m spatial resolution is not sufficient. In addition, a multi-sensor approach using both Sentinel3/2 is also suitable for the operational complete monitoring over broader coastal regions. The distinct improved performance of the S2 MSI and S3 OLCI (EFR) over Landsat-8 OLI is expected due to the lack of necessary optical feature along the red edge spectral region (700-750 nm). Furthermore, the phytoplankton usually forms very elongated and thin slicks of varying widths, a few to tens of meters, so S3 OCLI's spatial resolution (300 m) is not sufficient to resolve the patchiness which has been achieved by S2 MSI, due to its 10 m spatial resolution. The advantage of using next-generation optical sensors such as S2 to supplement the information gathered from in situ observations of water quality dynamics is of key importance (e.g., [106,107]).

Conclusions
In this study, the efficiency of multi-sensor data was tested against reference standard CMEMS, OLCI WFR, and the MODIS measurements. Despite unavoidably owing to the lack of ground-truth data, the inter-comparison of the multi-sensor and reference data proved promising in spatially quantifying the water quality estimates, where the average values and range levels of water quality parameters of the coastal water were very close for all sensor and reference data. For all data, Chl-a, TSM, and CDOM averaged 3.14 mg m −3 , 7.66 g m −3 , and 0.18 m −1 with ranges of 0.39-4.81 mg m −3 , 6.32-10.18 g m −3 and 0.13-0.30 m −1 , respectively. Spatial patterns of retrieved water content proved congruent with similar trends in the coastal water and varied widely in the inland lake and fish farm waters. The uneven distribution of water content showed local highs marked the inland lakes and their downstream areas in the coastal water, evidently in front of the Rosetta and Damietta promontories and Boughaz El-Burullus.
Intensive eutrophication is clarified by comparing the coastal water quality estimates to that of the inland Burullus Lake water. Additionally, MSI-retrieved Chl-a when validated against in situ data highlighted the common factors leading to the C2RCC underestimation of the retrieved parameters in the inland Burullus Lake water. Overall, Chl-a is retrieved better than other water constituents for all investigated data followed by TSM and CDOM. Local water salinity and temperature inputs into the C2RCC clarified the suitability of the MSI products, which came first in rank achieving the lowest MAE and RMSE for the Chl-a retrieval in the coastal water that was followed by the EFR and then by the OLI. EFR retrieved the TSM and CDOM much more effectively compared to the MSI and OLI.
In Burullus Lake, validation of MSI Chl-a retrieval against ground truth data confirmed a clear underestimation in the level range of 0.1-59% and average of 35.35%. The error increases at larger Chl-a concentrations, which is more affected by increasing TSM and CDOM contents. The RMSE estimated in the inland lake showed a more varied and strong effect on the error magnitude of the lake water's physical characteristics, such as turbidity loaded by nutrients along with composition and abundance of the Chlorophyacea, Euglenophycaea, and Bacillariophycaea, which clarify the largest biological contents of carbohydrates, lipids, and proteins. This research demonstrated the strong advantage of S2 MSI (10 m) and the S3 OLCI (300 m) data for the retrieval of Chl-a and the TSM/CDOM, respectively, as well as their better spectral, spatial, and temporal resolution than that of OLI. Therefore, in terms of their frequency and synoptic observations, the value added by the Copernicus products is of paramount importance for water quality monitoring plans and for ecological and management purposes at local, regional, and national scales.
Ocean observation technology is ever evolving; our analysis consolidates a rationale for water color estimates indicative of the pattern distribution of phytoplankton biomass Chl-a content that remains the best proxy for studies of primary productivity with successful application on the Nile Delta coastal waters where in situ measurements are lacking. Hot spots of maximum chlorophyll-a are related to the suitable growth conditions (nutrient input during column mixing periods), and reflect an increase in algal biomass and the primary productivity in front of watercourses discharging large nutrient loads into the sea.
The research keeps abreast of the advanced atmospheric correction algorithm, C2RCC, for retrieval of the coastal water quality parameters from three operational ocean and land color images within a 45 min overpass time on the Nile Delta within 4 h on 22 March 2020. The data inter-comparison enabled the evaluation of the algorithm with success in the coastal waters and highlighted the reasons for its pitfall in the inland lake waters that can be used for monitoring long-term changes. Moreover, results of the remotely-sensed water quality products prove indispensable for further community-wide validations and ensure performance under different environmental conditions.
In conclusion and compared to previous studies, the current study further confirms the suitability of the visible and NIR bands in the prediction of chlorophyll-a in complex case 2 waters. The concentration of chlorophyll-a presenting the active optical properties in the visible and near-infrared wavelengths, estimated using the C2RCC, indicated the abundance of algae in the area.
Future research will be devoted to the detailed estimation of Burullus Lake's water color parameters using the same three sensor data sets employed in this research at different seasonal times but with in situ spectroscopic measurements of the water physico-chemical characteristics and the algal biomass types, abundances, and compositional biological contents. Furthermore, various atmospheric correction algorithms will be evaluated. This is to find out key relationships among these parameters and investigate their coherence with those data of possible determination from the remote sensing data in this region, which is undergoing fast ecological changes in the context of a changing climate leading to environmental, social, and economic impacts. It has a promising potential to be periodically employed for cost-effective water-quality monitoring and leads to a reliable and continuously updated database for better water management plans in a GIS environment.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/w14040593/s1, Figure S1: Retrieved water color parameters of (a) Chl-a, TSM, and CDOM from S3 OCLI (WFR), and (b) their associated errors, Figure S2: Retrieved water color parameters of (a) Chl-a, TSM, and CDOM from S3 OCLI (EFR), and (b) associated errors, Figure S3: Retrieved water color parameters of (a) Chl-a, TSM, and CDOM from OLI, and (b) and their associated errors, Figure S4:  Data Availability Statement: The data sets generated and analyzed during the current study are available from the author upon request.