Multiplatform Earth Observation Systems for Monitoring Water Quality in Vulnerable Inland Ecosystems: Maspalomas Water Lagoon

: The accurate monitoring of water quality indicators, bathymetry and distribution of benthic habitats in vulnerable ecosystems is key to assessing the e ﬀ ects of climate change, the quality of natural areas and to guide appropriate biodiversity, tourism or ﬁsheries policies. Coastal and inland water ecosystems are very complex but crucial due to their richness and primary production. In this context, remote sensing can be a reliable way to monitor these areas, mainly thanks to satellite sensors’ improved spatial and spectral capabilities and airborne or drone instruments. In general, mapping bodies of water is challenging due to low signal-to-noise (SNR) at sensor level, due to the very low reﬂectance of water surfaces as well as atmospheric e ﬀ ects. Therefore, the main objective of this work is to provide a robust processing framework to estimate water quality parameters in inland shallow waters using multiplatform data. More speciﬁcally, we measured chlorophyll concentrations ( Chl-a ) from multispectral and hyperspectral sensors on board satellites, aircrafts and drones. The Natural Reserve of Maspalomas, Canary Island (Spain), was chosen for the study because of its complexity as well as being an inner lagoon with considerable organic and inorganic matter and chlorophyll concentration. This area can also be considered a well-known coastal-dune ecosystem attracting a large amount of tourists. The water quality parameter estimated by the remote sensing platforms has been validated using co-temporal in situ measurements collected during ﬁeld campaigns, and quite satisfactory results have been achieved for this complex ecosystem. In particular, for the drone hyperspectral instrument, the root mean square error, computed to quantify the di ﬀ erences between the estimated and in situ chlorophyll-a concentrations, was 3.45 with a bias of 2.96.


Introduction
The protection of natural ecosystems is key to preserving biodiversity [1]. Unfortunately, bodies of water are very sensitive to natural phenomena, invasive species, human activities and other factors [2,3]. However, to reach their full potential, periodic monitoring of coastal and inner water ecosystems, and the analysis of multi-temporal changes are important for the sustainable management of natural resources as well as to understand our environment.
Monitoring the water quality is of great importance, since different concentrations of chlorophyll-a (Chl-a), suspended matter or CDOM (Colored Dissolved Organic Matter) can have a major effect on the biological activity in aquatic ecosystems [4]. For instance, very high CDOM concentrations decreases important touristic pressure and, specifically, the lagoon is characterized by significant concentrations of inorganic and organic matter derived from the seabed and heavy rains. High concentrations of chlorophyll and materials directly affect nutrient dynamics, phytoplankton productivity and the presence of aquatic vegetation [17].
In summary, contributions are two-fold. On one hand, regarding the most accurate correction methods applied to VHRMS satellite data and, on the other hand, a comparative assessment of multiplatform (satellite, aircraft and drone) remote sensing imagery to obtain accurate, high resolution water quality parameters, specifically Chl-a, in complex inner water areas.
This paper is structured as follows: first, we present the study area, as well as the satellite, airborne, drone and in situ data used; next, we detail the processing methodology for the mapping of water quality parameters using multiplatform multispectral and hyperspectral data; then, the main results are included; finally, we highlight the central conclusions.

Study Area
The Natural Reserve of Maspalomas, located in the south of Gran Canaria, is 403.9 hectares and includes a system of mobile dunes and a small lagoon known as "la Charca de Maspalomas" (Figure 1). The Maspalomas lagoon was formed when the flow from the Fataga ravine decreased due to climate change. This small lagoon combines salt water with that from the underground streams of the ravine and occasional runoffs, forming a unique ecosystem in the Canary Islands. The breakdown of the ecological balance of this area has triggered the disappearance or decrease of both animal and plant species, caused by the surrounding buildings, wastewater effluents, canyon canalization and transit of tourists and vehicles [18]. Therefore, the huge amount of tourists (2 million visitors each year) exerts a very high anthropic effect on this vulnerable ecosystem.
Remote Sens. 2020, 12, x FOR PEER REVIEW 3 of 18 concentrations of inorganic and organic matter derived from the seabed and heavy rains. High concentrations of chlorophyll and materials directly affect nutrient dynamics, phytoplankton productivity and the presence of aquatic vegetation [17]. In summary, contributions are two-fold. On one hand, regarding the most accurate correction methods applied to VHRMS satellite data and, on the other hand, a comparative assessment of multiplatform (satellite, aircraft and drone) remote sensing imagery to obtain accurate, high resolution water quality parameters, specifically Chl-a, in complex inner water areas.
This paper is structured as follows: first, we present the study area, as well as the satellite, airborne, drone and in situ data used; next, we detail the processing methodology for the mapping of water quality parameters using multiplatform multispectral and hyperspectral data; then, the main results are included; finally, we highlight the central conclusions.

Study Area
The Natural Reserve of Maspalomas, located in the south of Gran Canaria, is 403.9 hectares and includes a system of mobile dunes and a small lagoon known as "la Charca de Maspalomas" ( Figure  1). The Maspalomas lagoon was formed when the flow from the Fataga ravine decreased due to climate change. This small lagoon combines salt water with that from the underground streams of the ravine and occasional runoffs, forming a unique ecosystem in the Canary Islands. The breakdown of the ecological balance of this area has triggered the disappearance or decrease of both animal and plant species, caused by the surrounding buildings, wastewater effluents, canyon canalization and transit of tourists and vehicles [18]. Therefore, the huge amount of tourists (2 million visitors each year) exerts a very high anthropic effect on this vulnerable ecosystem. As show in Figure 1, "La Charca de Maspalomas" is considered a breathtaking oasis with an outstanding floristic and ornithological richness. One of the main characteristics of this ecosystem is As show in Figure 1, "La Charca de Maspalomas" is considered a breathtaking oasis with an outstanding floristic and ornithological richness. One of the main characteristics of this ecosystem is its seasonality, with large variations in the salinity and volume of water, as well as in its levels of phytoplankton, dissolved and suspended matter. For this reason, within the framework of this research, a detailed study of water quality has been carried out in the water lagoon´s natural ecosystem with the latest technology using high-resolution images gathered by multiple sensors on board satellites, airplanes and drones, as well as field measurements.

Multisensor Remotely Sensed and In Situ Field Data
Satellite data obtained from the WorldView multispectral sensor were planned to be acquired at the same time as the airborne and drone overflights. Thus, VHR WorldView-3 (WV-3) satellite imagery was considered in the analysis. The level 2 ortho-ready product has been used, with a radiometric resolution of 11 bits and a spatial resolution of 1.8/1.6 m, at the nadir, for the eight multispectral bands (0.400-1.040 µm), as provided in Table 1 [19].
The flight campaign of 2 June 2017 used the INTA C-212 aircraft, which was equipped with the Airborne Hyperspectral Scanner (AHS). This radiometer records data with 12-bits resolution and covers the 0.43 to 12.8 µm range with 80 channels. In our study, which focuses on water areas, only the first 20 channels were used, covering the visible and near infrared (VIS-NIR) bands from 0.43 to 1.015 µm, with 2.5 m spatial resolution [20].
Finally, drone flights were performed on 4 June 2018 and 23 July 2019. In this case, the PIKA-L sensor was configured to have a spatial resolution of 10 cm × 10 cm, using 150 channels in the range 400-1000 nm, with a bandwidth of 4 nm. Table 1 includes a summary of the spectral and spatial characteristics of the multiplatform sensors.  Figure 2 shows the remote sensing platforms and the corresponding imagery of the Maspalomas lagoon water ecosystem. Specifically, Figure 2a highlights the multiplatform Earth observation systems used in this study, while Figure 2b shows the Worldview, AHS and drone true color composite imagery of the study area processed in this work, respectively. In situ field data was acquired simultaneously to the AHS and drone campaigns. Due to access restrictions to the area of interest, during the first campaigns, only locations near the shore could be sampled. Specifically, the following campaigns were performed (see Figure 3 and Table 2):

•
In June 2017, in situ data was acquired simultaneously to the AHS campaign. As shown in Figure  3, a total of three points were sampled at three different depths to measure water quality parameters (red points) using an IDRONAUT multiparameter probe. Precise information about the latitude, longitude and time were provided by a Trimble DSM132 Global Navigation Satellite System(GNSS) receiver. Additional sites monitored during 2015 were also used in the analysis. During the 2017 campaign, reflectance measurements were collected using the ASD FieldSpec3 precision radiometer. The procedure described in [21] was followed, measuring the radiance emanating from the water surface at a zenith angle between 30° and 50° and with an azimuth angle between 90° and 180° relative to the sun's azimuth. The procedure was replicated five times and the mean and standard deviation value were derived. A MICROTOPS II solar light meter was also used and a subsequent processing was performed to get the reflectance [22]. In situ field data was acquired simultaneously to the AHS and drone campaigns. Due to access restrictions to the area of interest, during the first campaigns, only locations near the shore could be sampled. Specifically, the following campaigns were performed (see Figure 3 and Table 2):

•
In June 2017, in situ data was acquired simultaneously to the AHS campaign. As shown in Figure 3, a total of three points were sampled at three different depths to measure water quality parameters (red points) using an IDRONAUT multiparameter probe. Precise information about the latitude, longitude and time were provided by a Trimble DSM132 Global Navigation Satellite System(GNSS) receiver. Additional sites monitored during 2015 were also used in the analysis. During the 2017 campaign, reflectance measurements were collected using the ASD FieldSpec3 precision radiometer. The procedure described in [21] was followed, measuring the radiance emanating from the water surface at a zenith angle between 30 • and 50 • and with an azimuth angle between 90 • and 180 • relative to the sun's azimuth. The procedure was replicated five times and the mean and standard deviation value were derived. A MICROTOPS II solar light meter was also used and a subsequent processing was performed to get the reflectance [22].

•
In June, 2018, another field campaign was performed along with a drone flight. As shown in Figure 3, a total of three points were sampled (blue points) near to the lagoon shore measuring turbidity, chlorophyll-a/b and carothenoids using a Hydrolab HL4 multiparameter probe. A HYPER-V GNSS was used to get the time and precise coordinates. • During July, 2019, in situ data was also collected simultaneously to the drone flight. A total of seven points were sampled (black in Figure 3); however, for the two locations in the center of the lagoon, samples at two depths were collected. In this campaign, more specific water parameters were measured; specifically, Chl-a, Chl-b, turbidity, CDOM, carotenoids, phycobiliproteins and functional groups of planktonic microalgae. In addition, bathymetry and the geographic coordinates and time of each sampled were obtained.
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 18 • In June, 2018, another field campaign was performed along with a drone flight. As shown in Figure 3, a total of three points were sampled (blue points) near to the lagoon shore measuring turbidity, chlorophyll-a/b and carothenoids using a Hydrolab HL4 multiparameter probe. A HYPER-V GNSS was used to get the time and precise coordinates. • During July, 2019, in situ data was also collected simultaneously to the drone flight. A total of seven points were sampled (black in Figure 3); however, for the two locations in the center of the lagoon, samples at two depths were collected. In this campaign, more specific water parameters were measured; specifically, Chl-a, Chl-b, turbidity, CDOM, carotenoids, phycobiliproteins and functional groups of planktonic microalgae. In addition, bathymetry and the geographic coordinates and time of each sampled were obtained.

Multiplatform Water Quality Monitoring
The overall processing protocol to generate high resolution water quality parameters for multiplatform imagery is presented in Figure 4. This work deals exclusively with the water quality monitoring; thus, the retrieval and evaluation of the high resolution bathymetry has not been addressed. Next, the different steps involved are described in more detail.

Multiplatform Water Quality Monitoring
The overall processing protocol to generate high resolution water quality parameters for multiplatform imagery is presented in Figure 4. This work deals exclusively with the water quality monitoring; thus, the retrieval and evaluation of the high resolution bathymetry has not been addressed. Next, the different steps involved are described in more detail.

Radiometric and Atmospheric Corrections
Advanced Earth observation platforms and sensors, with enhanced capabilities, can contribute relevant and accurate information for the sustainability of vulnerable water ecosystems [23]. In particular, Worldview-2/3 platforms offer VHR imagery with additional VIS-NIR spectrum channels not included in former high resolution satellites (i.e., Geoeye, Quickbird, Ikonos, KOMPSAT, Pleiades etc.). Unfortunately, the energy collected by space sensors is influenced by atmospheric disturbances that hinders the accuracy at surface level. On the other hand, the new hyperspectral high-resolution instruments in airborne platforms, with new spectral bands, as well as the ability of some platforms to modify their sensing angle has highlighted the need to address atmospheric effects [20].
As the mapping of water areas is challenging, due to the low signal-to-noise received at sensor level, as a consequence of the minimum reflectance of water surfaces and atmospheric effects, precise preprocessing algorithms have to be applied in order to correct limitations in the sensor calibration, solar illumination geometry and viewing effects, as well as the atmospheric and sunglint disturbances.
Furthermore, as previously indicated, in bodies of water, radiometric and atmospheric corrections are fundamental steps in the VHR satellite data preprocessing protocol due to the low signal reaching the sensor from water surfaces. However, in some environments such as very shallow inner-lakes, radiance leaving the water surface may be considerably higher as a result of the dissolved and suspended matter back-scattering and the bottom albedo. Consequently, applying models developed for open ocean or coastal waters to remote sensing imagery often provide erroneous estimates and, therefore, complex atmospheric correction algorithms are required for inner shallow waters [24,25].
Specifically, in Eugenio et al. [10] for the Maspalomas lagoon ecosystem, a comparison was performed between the estimated reflectivity, provided by the application of different atmospheric correction algorithms, and the real reference spectral signature measured with a spectroradiometer, acquired at the time of the satellite overflight. It was demonstrated that model-based algorithms properly correct the atmospheric disturbances.
After a comprehensive assessment (see results section), for the WorldView-2/3 imagery we have selected the FLAASH (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes) atmospheric model properly adapted to the multispectral satellite channels. On the other hand, for the airborne AHS data, ATCOR-4 (Atmospheric CORrection) was selected to correct atmospheric and illumination disturbances. Finally, atmospheric corrections were not applied to the drone Pika-L

Radiometric and Atmospheric Corrections
Advanced Earth observation platforms and sensors, with enhanced capabilities, can contribute relevant and accurate information for the sustainability of vulnerable water ecosystems [23]. In particular, Worldview-2/3 platforms offer VHR imagery with additional VIS-NIR spectrum channels not included in former high resolution satellites (i.e., Geoeye, Quickbird, Ikonos, KOMPSAT, Pleiades etc.). Unfortunately, the energy collected by space sensors is influenced by atmospheric disturbances that hinders the accuracy at surface level. On the other hand, the new hyperspectral high-resolution instruments in airborne platforms, with new spectral bands, as well as the ability of some platforms to modify their sensing angle has highlighted the need to address atmospheric effects [20].
As the mapping of water areas is challenging, due to the low signal-to-noise received at sensor level, as a consequence of the minimum reflectance of water surfaces and atmospheric effects, precise preprocessing algorithms have to be applied in order to correct limitations in the sensor calibration, solar illumination geometry and viewing effects, as well as the atmospheric and sunglint disturbances.
Furthermore, as previously indicated, in bodies of water, radiometric and atmospheric corrections are fundamental steps in the VHR satellite data preprocessing protocol due to the low signal reaching the sensor from water surfaces. However, in some environments such as very shallow inner-lakes, radiance leaving the water surface may be considerably higher as a result of the dissolved and suspended matter back-scattering and the bottom albedo. Consequently, applying models developed for open ocean or coastal waters to remote sensing imagery often provide erroneous estimates and, therefore, complex atmospheric correction algorithms are required for inner shallow waters [24,25].
Specifically, in Eugenio et al. [10] for the Maspalomas lagoon ecosystem, a comparison was performed between the estimated reflectivity, provided by the application of different atmospheric correction algorithms, and the real reference spectral signature measured with a spectroradiometer, acquired at the time of the satellite overflight. It was demonstrated that model-based algorithms properly correct the atmospheric disturbances.
After a comprehensive assessment (see results section), for the WorldView-2/3 imagery we have selected the FLAASH (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes) atmospheric model properly adapted to the multispectral satellite channels. On the other hand, for the airborne AHS data, ATCOR-4 (Atmospheric CORrection) was selected to correct atmospheric and illumination disturbances. Finally, atmospheric corrections were not applied to the drone Pika-L hyperspectral imagery. The following is a brief description of the atmospheric models used in the Maspalomas lagoon.
where A, B, S and L o are determined from the MODTRAN-4 simulations for the selected model of atmosphere (maritime, in our case). S is the spherical albedo, L o is the radiance backscattered by the atmosphere, and A, and B are coefficients related to the geometric and atmospheric conditions. ρ SU is the surface reflectance, ρ e is the mean surface reflectance for the pixel and adjacent regions. The first term in (1) refers to the direct radiation reaching the sensor from the surface, while the second term relates to the radiance received by the sensor coming from the surface that is scattered by the atmosphere. ATCOR-4 is used for atmospheric correction of sensors flying at altitudes up to 20 km and the database was compiled with the MODTRAN-5 radiative transfer code, covering a wide range of weather conditions and sun angles. The spectral sampling distance is 0.4 nm, equidistant over the range from 350 nm to 2550 nm. Regarding the accuracy of the method, in the solar region and assuming the terrain is flat, the estimated surface reflectance errors range between ±2% (reflectance < 10%) and ±4% (reflectance > 40%). For rugged terrain, accuracy depends on different parameters such as the spatial resolution of the digital elevation model or the ortho-rectification precision [20,27].
For the ATCOR algorithm, also based on MODTRAN-4, the surface reflectance, without taking into account the adjacency effect, is found by: where a 0 and a 1 coefficients are obtained from the estimation of the main atmospheric parameters: water vapor column, aerosol type and optical thickness. As indicated, for the drone-carried hyperspectral system, the flight altitude is around 100 m, so as to simplify the atmospheric correction process using the hypothesis that the irradiance at the height of the drone is equal to the irradiance reaching the Earth's surface. In order to obtain the irradiance measurement, the drone carries a radiometer that points to the zenith with a cosine integrator, calculating the total irradiance in the 400-1000 nm range. With the information of the radiance and irradiance sensed simultaneously, we can get the reflectivity easily and with great accuracy (see Figure 5).
Remote Sens. 2020, 12, x FOR PEER REVIEW 8 of 18 hyperspectral imagery. The following is a brief description of the atmospheric models used in the Maspalomas lagoon. FLAASH [26] derives its physics-based model from the MODTRAN-4 radiative code. FLAASH estimates the reflectance at the surface by removing the absorption and scattering disturbances by: where , , and are determined from the MODTRAN-4 simulations for the selected model of atmosphere (maritime, in our case). is the spherical albedo, is the radiance backscattered by the atmosphere, and , and are coefficients related to the geometric and atmospheric conditions. is the surface reflectance, is the mean surface reflectance for the pixel and adjacent regions. The first term in (1) refers to the direct radiation reaching the sensor from the surface, while the second term relates to the radiance received by the sensor coming from the surface that is scattered by the atmosphere.
ATCOR-4 is used for atmospheric correction of sensors flying at altitudes up to 20 km and the database was compiled with the MODTRAN-5 radiative transfer code, covering a wide range of weather conditions and sun angles. The spectral sampling distance is 0.4 nm, equidistant over the range from 350 nm to 2550 nm. Regarding the accuracy of the method, in the solar region and assuming the terrain is flat, the estimated surface reflectance errors range between ±2% (reflectance < 10%) and ±4% (reflectance > 40%). For rugged terrain, accuracy depends on different parameters such as the spatial resolution of the digital elevation model or the ortho-rectification precision [20,27].
For the ATCOR algorithm, also based on MODTRAN-4, the surface reflectance, without taking into account the adjacency effect, is found by: where and coefficients are obtained from the estimation of the main atmospheric parameters: water vapor column, aerosol type and optical thickness.
As indicated, for the drone-carried hyperspectral system, the flight altitude is around 100 m, so as to simplify the atmospheric correction process using the hypothesis that the irradiance at the height of the drone is equal to the irradiance reaching the Earth's surface. In order to obtain the irradiance measurement, the drone carries a radiometer that points to the zenith with a cosine integrator, calculating the total irradiance in the 400-1000 nm range. With the information of the radiance and irradiance sensed simultaneously, we can get the reflectivity easily and with great accuracy (see Figure 5).

Water Column Modeling and Chlorophyll-A Estimation
For bodies of water, reflectivity recorded by the sensor, once the atmospheric effects are corrected, is gathered by two sources: the specular reflection of the surface of the water and the light backscattered contributed by the water column. Only the upward reflectivity of the water column provides information about the physic-optical phenomena produced within the aquatic environment. In coastal areas these phenomena can be modeled by the inherent optical properties (IOPs) of the water and modeling the bottom albedo taking into account its depth [28][29][30]. The complete expression of the radiative model is: where r m rs (λ) is the modeled reflectance below the sea surface, r m rs,∞ (λ) is the reflectance below the sea surface for deep waters, ρ alb (λ) is the seafloor albedo considering the Lambertian hypothesis, z represents the bathymetry, k d is the diffuse attenuation coefficient of the water, µ sw s is the cosine of the solar inclination angle and µ sw v is the cosine of the sensor vision angle and, finally, D c u is the light diffusion ascending factor due to the water backscattering, while D b u is due to the the bottom albedo. The model variables related to the coastal bottom are both the bottom albedo and the bathymetry, where the albedo is generated by a linear mixing of representative benthic area classes [31]: where abu i is the abundance of the benthic class i, bento i is the reflectivity of the benthic class i and n em is the number of benthic classes considered in the model. On the other hand, the parameters related to the inherent properties of water and, therefore, the elements in solution and suspension found in the water are r m rs,∞ and k d . The diffuse attenuation is the sum of the water absorption (a) and backscatter (b b ): r m rs,∞ is computed from the Gordon parameter u [32]: Albert and Mobley [33] proposed the following equation for the calculation of the r m rs parameter, where the sun and sensor angles are also taken into account: Finally, absorption and backscattering results from contributions from the water plus the sum of the different elements in solution and suspension are: where a w is the absorption of water molecules, a p the phytoplankton absorption, a gd the joint absorption of dissolved matter and detritus, b bw is the molecular backscattering of water (Rayleigh scattering) and b bp is the backscattering of suspended matter (Mie scattering). While a w and b bw are known values, the other values depend on the concentration of the different elements (phytoplankton concentration P, dissolved matter concentration and detritus G and particulate matter concentration M). In this way, we can propose a non-linear system of equations with as many equations as available spectral bands in the image and with as many variables as the modeled parameters (P, G, M, z, Nabu). Each of these equations will match the sensor bands by integrating the result of the equation across the entire bandwidth: where BandPass(λ) is the normalized bandpass filter of each sensor channel. Integration is performed with 5 nm deltas for multispectral bands, while monochromatic is assumed for hyperspectral sensor bands. To solve the model variables, it is necessary to compare the modeled reflectivity r m rs with the measured reflectivity r rs . For this, we have to convert the water surface into r rs . First we must eliminate the specular brightness of the water using the following equation, as studied in Martin et al. [9]: where ρ glint (band i ) is the surface reflectivity of band i affected by the specular brightness, Idir(band i ) is the direct irradiance of band i, Idir(NIR) is the irradiance of the NIR band and ρ m glint (NIR) is the surface reflectivity of the modeled NIR band.
The removal of sunglint is carried out by subtracting the solar brightness detected in the NIR band, where the absorption of water means that, practically, all the reflectivity detected comes from the sun´s specular brightness. However, due to the very high turbidity and low depth of the study area, we cannot neglect the inherent reflectivity contribution of the NIR band. To take this fact into account, the contribution obtained by the radiative model is subtracted from the NIR band. Therefore, it is necessary to introduce the specular solar brightness removal equation in the radiative transfer model resolution. Values in one-for-one of direct irradiance Idir(band i ) are obtained from the atmospheric model for the atmospheric conditions of the selected day. An example of the direct irradiance value for the 400-900 nm range is shown in Figure 6.
where ( ) is the normalized bandpass filter of each sensor channel. Integration is performed with 5 nm deltas for multispectral bands, while monochromatic is assumed for hyperspectral sensor bands. To solve the model variables, it is necessary to compare the modeled reflectivity with the measured reflectivity . For this, we have to convert the water surface into . First we must eliminate the specular brightness of the water using the following equation, as studied in Martin et al. [9]: where ( ) is the surface reflectivity of band i affected by the specular brightness, ( ) is the direct irradiance of band i, ( ) is the irradiance of the NIR band and ( ) is the surface reflectivity of the modeled NIR band.
The removal of sunglint is carried out by subtracting the solar brightness detected in the NIR band, where the absorption of water means that, practically, all the reflectivity detected comes from the sun´s specular brightness. However, due to the very high turbidity and low depth of the study area, we cannot neglect the inherent reflectivity contribution of the NIR band. To take this fact into account, the contribution obtained by the radiative model is subtracted from the NIR band. Therefore, it is necessary to introduce the specular solar brightness removal equation in the radiative transfer model resolution. Values in one-for-one of direct irradiance ( ) are obtained from the atmospheric model for the atmospheric conditions of the selected day. An example of the direct irradiance value for the 400-900 nm range is shown in Figure 6.
Taking all this into account, the problem can be represented regardless of the type of sensor used, via a system of equations for each band: Finally, ρ(band i ) can be converted to r rs by: ρ(band i ) = π * 0.5374r rs 1 − 0.3124r rs (12) Remote Sens. 2020, 12, 284

of 18
Taking all this into account, the problem can be represented regardless of the type of sensor used, via a system of equations for each band: . . .
The system of equations´solution is not easy because we have more equations than unknowns (especially in hyperspectral images). Furthermore, as it is a non-linear problem, there are multiple solutions for each unknown, so it is impossible to solve the problem directly. To solve it, a gradient optimizer (Levenberg-Marquardt) has been used [34]. This algorithm is based on the minimization of a cost function, which is the mean square error of the difference between the measured and modeled reflectivity: where ∆ min is the cost function to be minimized by the optimizer and n bands is the number of bands used.
Finally, chlorophyll-a concentration can be estimated from remotely sensed data by relating optical changes observed in the reflected light at specific wavelengths. Multispectral data from low resolution satellite sensors, such as MODIS, MERIS, SENTINEL-3 etc., can be used to estimate the water quality of inner waters [35]. However, to monitor a complex ecosystem, such as the Maspalomas lagoon, VHR multispectral or hyperspectral imagery from different platforms is required. As better sensors with greater number of channels and narrower bandwidths are used, it is possible to more clearly distinguish the absorption and backscattering contributions of the water quality parameters. In particular, chlorophyll-a concentration can be estimated from remotely sensed data by relating optical changes observed in the reflected light at specific wavelengths. As presented in Figure 7, photosynthetic pigments can be distinguished thanks to the two absorption peaks of Chl-a.
The system of equations´ solution is not easy because we have more equations than unknowns (especially in hyperspectral images). Furthermore, as it is a non-linear problem, there are multiple solutions for each unknown, so it is impossible to solve the problem directly. To solve it, a gradient optimizer (Levenberg-Marquardt) has been used [34]. This algorithm is based on the minimization of a cost function, which is the mean square error of the difference between the measured and modeled reflectivity: where Δ is the cost function to be minimized by the optimizer and is the number of bands used.
Finally, chlorophyll-a concentration can be estimated from remotely sensed data by relating optical changes observed in the reflected light at specific wavelengths. Multispectral data from low resolution satellite sensors, such as MODIS, MERIS, SENTINEL-3 etc., can be used to estimate the water quality of inner waters [35]. However, to monitor a complex ecosystem, such as the Maspalomas lagoon, VHR multispectral or hyperspectral imagery from different platforms is required. As better sensors with greater number of channels and narrower bandwidths are used, it is possible to more clearly distinguish the absorption and backscattering contributions of the water quality parameters. In particular, chlorophyll-a concentration can be estimated from remotely sensed data by relating optical changes observed in the reflected light at specific wavelengths. As presented in Figure 7, photosynthetic pigments can be distinguished thanks to the two absorption peaks of Chla.

Satellite/Airborne Atmospheric Algorithms Assessment
Next, we presented a summary of the most relevant results, in the atmospheric and deglinting modeling context, for the Maspalomas inner-lake, using imagery from Worldview and airborne platforms.

Satellite/Airborne Atmospheric Algorithms Assessment
Next, we presented a summary of the most relevant results, in the atmospheric and deglinting modeling context, for the Maspalomas inner-lake, using imagery from Worldview and airborne platforms. Figure 8 shows the water reflectance, for the visible and infrared satellite and airborne channels, in the lagoon, compared to in situ spectral data (Figure 3a left) measured simultaneously with remote sensors overflight. Both graphs show the mean reflectance averaged for the three sampled points measured during the 2017 field campaign.
Following this comparison, the most accurate results from the average measurements for the considered inner-lake points are included in Table 3 (see Table 2), in which we present the Root Mean Square Error(RMSE) and BIAS between the WorldView and AHS sensors corrected reflectance and in situ measurement. The equations used to calculate the error criteria are as follows: where ρ in−situ is the measured in situ surface reflectance and ρ sat/airbone is the surface reflectance estimated from the satellite/airborne data.
Remote Sens. 2020, 12, x FOR PEER REVIEW 12 of 18 Figure 8 shows the water reflectance, for the visible and infrared satellite and airborne channels, in the lagoon, compared to in situ spectral data (Figure 3a left) measured simultaneously with remote sensors overflight. Both graphs show the mean reflectance averaged for the three sampled points measured during the 2017 field campaign.
Following this comparison, the most accurate results from the average measurements for the considered inner-lake points are included in Table 3 (see Table 2), in which we present the Root Mean Square Error(RMSE) and BIAS between the WorldView and AHS sensors corrected reflectance and in situ measurement. The equations used to calculate the error criteria are as follows: where is the measured in situ surface reflectance and / is the surface reflectance estimated from the satellite/airborne data.    Finally, as previously indicated, removal of sunglint is necessary for the reliable retrieval of water quality parameters in inner-water environments. Generally speaking, sunglint techniques are based on the assumption that water reflectivity in NIR channels is always negligible. However, as studied by Martin et al. [9] that is not the case in very shallow environments or in waters with high concentrations of turbidity. As these are common situations in coastal and inland waters, there can be significant errors. Thus, we proposed to integrate the sunglint removal algorithm in the radiative transfer model to estimate the contribution of the coastal waters NIR reflectance, which will allow us to eliminate specular NIR reflectance contribution. Next, Figure 9 shows an example of a hyperspectral spectrum of the Maspalomas lagoon, where the reference measured at the laboratory is depicted in blue, the modeled values generated using the radiative transfer model with the airborne data are included in green and the reflectivity after removing the sunglint is displayed in orange.
Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 18 Finally, as previously indicated, removal of sunglint is necessary for the reliable retrieval of water quality parameters in inner-water environments. Generally speaking, sunglint techniques are based on the assumption that water reflectivity in NIR channels is always negligible. However, as studied by Martin et al. [9] that is not the case in very shallow environments or in waters with high concentrations of turbidity. As these are common situations in coastal and inland waters, there can be significant errors. Thus, we proposed to integrate the sunglint removal algorithm in the radiative transfer model to estimate the contribution of the coastal waters NIR reflectance, which will allow us to eliminate specular NIR reflectance contribution. Next, Figure 9 shows an example of a hyperspectral spectrum of the Maspalomas lagoon, where the reference measured at the laboratory is depicted in blue, the modeled values generated using the radiative transfer model with the airborne data are included in green and the reflectivity after removing the sunglint is displayed in orange. Figure 9. Spectral fitting of the model reflectance and the measured-deglinted reflectance (reference in blue and modeled reflectivity before and after the sunglint correction are presented in green and orange, respectively).
As an example, Figure 10 shows multiplatform color composite image (drone and satellite) before (left) and after (right) the sunglint contamination removal.

Inner Lake Monitoring: Chlorophyll-A Concentration
The ability to retrieve water quality parameters using different sensors has been evaluated in Maspalomas (Gran Canaria Island), not only a UNESCO Biosphere Reserve, but where monitoring of water quality is of special interest because of its exuberance, marine life and tourist activities. In particular, a corrected WV-3 multispectral imagery from the lagoon has been processed to generate Figure 9. Spectral fitting of the model reflectance and the measured-deglinted reflectance (reference in blue and modeled reflectivity before and after the sunglint correction are presented in green and orange, respectively).
As an example, Figure 10 shows multiplatform color composite image (drone and satellite) before (left) and after (right) the sunglint contamination removal.
Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 18 Finally, as previously indicated, removal of sunglint is necessary for the reliable retrieval of water quality parameters in inner-water environments. Generally speaking, sunglint techniques are based on the assumption that water reflectivity in NIR channels is always negligible. However, as studied by Martin et al. [9] that is not the case in very shallow environments or in waters with high concentrations of turbidity. As these are common situations in coastal and inland waters, there can be significant errors. Thus, we proposed to integrate the sunglint removal algorithm in the radiative transfer model to estimate the contribution of the coastal waters NIR reflectance, which will allow us to eliminate specular NIR reflectance contribution.
Next, Figure 9 shows an example of a hyperspectral spectrum of the Maspalomas lagoon, where the reference measured at the laboratory is depicted in blue, the modeled values generated using the radiative transfer model with the airborne data are included in green and the reflectivity after removing the sunglint is displayed in orange. Figure 9. Spectral fitting of the model reflectance and the measured-deglinted reflectance (reference in blue and modeled reflectivity before and after the sunglint correction are presented in green and orange, respectively).
As an example, Figure 10 shows multiplatform color composite image (drone and satellite) before (left) and after (right) the sunglint contamination removal.

Inner Lake Monitoring: Chlorophyll-A Concentration
The ability to retrieve water quality parameters using different sensors has been evaluated in Maspalomas (Gran Canaria Island), not only a UNESCO Biosphere Reserve, but where monitoring of water quality is of special interest because of its exuberance, marine life and tourist activities. In particular, a corrected WV-3 multispectral imagery from the lagoon has been processed to generate

Inner Lake Monitoring: Chlorophyll-A Concentration
The ability to retrieve water quality parameters using different sensors has been evaluated in Maspalomas (Gran Canaria Island), not only a UNESCO Biosphere Reserve, but where monitoring of water quality is of special interest because of its exuberance, marine life and tourist activities. In particular, a corrected WV-3 multispectral imagery from the lagoon has been processed to generate a high-resolution chlorophyll-a concentration map. Furthermore, hyperspectral imagery from airborne Remote Sens. 2020, 12, 284 14 of 18 and drone platforms have also been processed to assess the benefits of each technology in this challenging scenario. The true color composites of the data used in the study are included in Figure 2b. Figure 11 presents the mapping results for the chlorophyll-a concentration (µg/L) in the small lake using the multisensor imagery. Figure 11a shows the chlorophyll-a concentration map obtained by the WorldView-3 satellite, Figure 11b presents the result for the AHS radiometer and Figure 11c the map for the Pika-L instrument.
a high-resolution chlorophyll-a concentration map. Furthermore, hyperspectral imagery from airborne and drone platforms have also been processed to assess the benefits of each technology in this challenging scenario. The true color composites of the data used in the study are included in Figure 2b. Figure 11 presents the mapping results for the chlorophyll-a concentration (µg/L) in the small lake using the multisensor imagery. Figure 11a shows the chlorophyll-a concentration map obtained by the WorldView-3 satellite, Figure 11b presents the result for the AHS radiometer and Figure 11c the map for the Pika-L instrument.  Table 4 offers a summary of the statistical analysis, providing quantitative information between estimated and real chlorophyll concentrations. Finally, it is important to highlight that during the 2019 field campaign, after requesting the corresponding authorization to the conservation authorities, on this occasion a small boat was permitted to sail inside the lagoon and, therefore, samples could be collected at different points far from the border (black dots in Figure 3b). In general, chlorophyll concentration values were in the 25 to 35 µg/L range. The 2019 campaign was carried out during summer and, as a result, water quality conditions for the lagoon could be similar to the dates the remote sensing data was collected. Sometimes, water quality in the lagoon can change considerably after heavy rain but the Maspalomas area usually has an arid climate with over 250 sunny days per year and with low rainfall (between 10 and 15 days with rain per year), mainly during the winter.

Discussion
Concerning to the Satellite/Airborne atmospheric algorithms assessment, as analyzed in Section 2.3.1 and shown, qualitatively, in Figure 8 and, quantitatively, in Table 3, a detailed assessment of physical model-based atmospheric correction algorithms has been performed, applied to the  Table 4 offers a summary of the statistical analysis, providing quantitative information between estimated and real chlorophyll concentrations. Finally, it is important to highlight that during the 2019 field campaign, after requesting the corresponding authorization to the conservation authorities, on this occasion a small boat was permitted to sail inside the lagoon and, therefore, samples could be collected at different points far from the border (black dots in Figure 3b). In general, chlorophyll concentration values were in the 25 to 35 µg/L range. The 2019 campaign was carried out during summer and, as a result, water quality conditions for the lagoon could be similar to the dates the remote sensing data was collected. Sometimes, water quality in the lagoon can change considerably after heavy rain but the Maspalomas area usually has an arid climate with over 250 sunny days per year and with low rainfall (between 10 and 15 days with rain per year), mainly during the winter.

Discussion
Concerning to the Satellite/Airborne atmospheric algorithms assessment, as analyzed in Section 2.3.1 and shown, qualitatively, in Figure 8 and, quantitatively, in Table 3, a detailed assessment of physical model-based atmospheric correction algorithms has been performed, applied to the multispectral bands of the high-resolution satellite for inner-lake. In this work, this study has been extended to airborne hyperspectral data.
An acceptable matching in the spectral signatures can be seen ( Figure 8) between measured data and those estimated by the atmospheric correction algorithms, showing better performance in the visible wavelengths. Following this comparison, atmospheric models with the best statistical parameters are FLAASH, for the satellite imagery (RMSE: 0.1014 and BIAS: −0.0348), and ATCOR for the airborne data (RMSE: 0.0318 and BIAS: −0.0251).
On the other hand, regarding relation to inner lake monitoring, as shown in Figure 11, a corrected WV-3 multispectral imagery and hyperspectral imagery from airborne and drone platforms from the lagoon have been processed to assess the benefits of each technology in this challenging scenario. In Figure 11, we can appreciate better detail in the drone map as a result of its higher spatial and spectral resolution. However, it is more affected by the bottom disturbance near the shore. The satellite multispectral data used in the analysis provides an average result as it cannot properly capture the Chl-a spatial variation.
As provided in Table 4, comparing the in situ chlorophyll-a concentrations measured with those estimated by the remote sensing platforms, we can see that it was fairly accurate (RMSE values ranging between 3.49 and 6.75). RMSE and BIAS values for the hyperspectral drone sensor are lower than for the airborne and satellite estimates. Compared to the quantitative results, such statistical values are taken into account for worst case scenarios, as the samplings performed during the 2017 and 2018 campaigns were taken very close to the shoreline and, therefore, the bottom reflectivity affected, to a greater extent, the chlorophyll concentrations collected.
In this line, an on-going research effort is currently being undertaken to estimate very high-resolution bathymetry in this shallow and complex water ecosystem to better estimate quality parameters employing accurate bathymetric information.

Conclusions
An efficient monitoring of water quality parameters is important in coastal and inland water ecosystems to guide conservation managers regarding environmental protection, fisheries, tourism etc. In this context, remote sensing can be a fundamental tool to analyze vulnerable or complex water ecosystems. In particular, high-resolution hyperspectral or multispectral imaging systems, with spectral channels in the visible spectrum, can provide information on water quality parameters at finer temporal and spatial scales than field observations alone.
In any case, estimating water quality with remote sensing data is very challenging because a very small amount of energy leaving the water's surface reaches the sensor and it is mainly affected by noise from different sources. Therefore, precise calibrations and corrections are required to take into account the sensors spectral response, viewing angles, solar illumination geometries, atmospheric effects and sunglint disturbances.
In this work, we chose the natural reserve of Maspalomas, Canary Islands (Spain). More specifically, a complex natural lagoon (La Charca de Maspalomas) of brackish waters, located at the edge of the dunes and the sea, which attracts many species of bird.
First, accurate preprocessing tasks have been carried out to retrieve the water leaving reflectance. In particular, three advanced atmospheric correction algorithms based on radiative transfer models (ATCOR, FLAASH and 6S) were applied and statistically evaluated with respect to real data. Indeed, for the Maspalomas lagoon, the reflectance estimated by the model-based atmospheric correction techniques was compared with the reference water leaving reflectance measured by a field radiometer signature measured with a spectroradiometer, simultaneously to the satellite overflight.
It was confirmed that the advanced correction algorithms assessed correctly model the atmospheric absorption and scattering disturbances. In particular, the ATCOR and FLAASH models performed properly in the Maspalomas lagoon and, accordingly, they have been implemented and adapted to the particularities of the specific area and sensor. Next, an enhanced deglinting algorithm was applied.
To gather water quality, specifically chlorophyll-a concentration, a multichannel physics-based model was implemented. The complex model used in this work expands a previous methodology based on a combination of water column properties, bathymetry retrieval and seafloor abundances.
Data provided by the multispectral (Worldview satellite) and hyperspectral (drone and airborne) imagery, together with the new methodology presented, performs reasonably when finding chlorophyll concentrations in a complex shallow-water ecosystem. Regarding the multiplatform comparison, the drone provides the highest spectral and spatial information and, as a result, the most accurate results.
In summary, the multiplatform processing protocol developed has provided a synoptic and systematic framework to support inland shallow water ecosystems. It has been validated using a database of real measurements sampled in field campaigns. Satisfactory results have been achieved, but to further increase the accuracy, new models are considered to improve gathering of bathymetric and seafloor information.
Author Contributions: All of the authors contributed to design the research. F.E. were responsible for the data preprocessing. J.M. (Javier Marcello) and J.M. (Javier Martín) designed and performed the in situ measurements (field campaign) and water quality parameter processing. Finally, F.E. wrote the overall paper with contributions from the rest of authors, who were also involved in the manuscript's discussion and revision. All authors have read and agreed to the published version of the manuscript.