Next Article in Journal
Monitoring Greenhouse Gases from Space
Next Article in Special Issue
Detection of Spatial and Temporal Patterns of Liana Infestation Using Satellite-Derived Imagery
Previous Article in Journal
Analysis of Agronomic Drought in a Highly Anthropogenic Context Based on Satellite Monitoring of Vegetation and Soil Moisture
Previous Article in Special Issue
A Remote Sensing Approach to Understanding Patterns of Secondary Succession in Tropical Forest
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water Mixing Conditions Influence Sentinel-2 Monitoring of Chlorophyll Content in Monomictic Lakes

1
Department of Spatial Sciences, Faculty of Environmental Sciences, Czech University of Life Sciences Prague, Kamýcká 129, 165 00 Praha Suchdol, Czech Republic
2
Department of Sciences, University Roma Tre, Viale Guglielmo Marconi 446, 00146 Rome, Italy
3
Department of Physical Geography and Geoecology, Faculty of Science, Charles University, Albertov 6, 128 00 Praha, Czech Republic
4
Italian National Agency for New Technologies, Energy and Sustainable Economic Development (ENEA) CR Casaccia, 00123 Rome, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(14), 2699; https://doi.org/10.3390/rs13142699
Submission received: 5 June 2021 / Revised: 2 July 2021 / Accepted: 6 July 2021 / Published: 8 July 2021
(This article belongs to the Special Issue Feature Paper Special Issue on Ecological Remote Sensing)

Abstract

:
Prompt estimation of phytoplankton biomass is critical in determining the ecological quality of freshwaters. Remote Sensing (RS) may provide new opportunities to integrate with situ traditional monitoring techniques. Nonetheless, wide regional and temporal variability in freshwater optical constituents makes it difficult to design universally applicable RS protocols. Here, we assessed the potential of two neural networks-based models, namely the Case 2 Regional CoastColour (C2RCC) processor and the Mixture Density Network (MDN), applied to MSI Sentinel-2 data for monitoring Chlorophyll (Chl) content in three monomictic volcanic lakes while accounting for the effect of their specific water circulation pattern on the remotely-sensed and in situ data relation. Linear mixed models were used to test the relationship between the remote sensing indices calculated through C2RCC (INN) and MDN (IMDN), and in situ Chl concentration. Both indices proved to explain a large portion of the variability in the field data and exhibited a positive and significant relationship between Chl concentration and satellite data, but only during the mixing phase. The significant effect of the water circulation period can be explained by the low responsiveness of the RS approaches applied here to the low phytoplankton biomass, typical of the stratification phase. Sentinel-2 data proved their valuable potential for the remote sensing of phytoplankton in small inland water bodies, otherwise challenging with previous sensors. However, caution should be taken, since the applicability of such an approach on certain water bodies may depend on hydrological and ecological parameters (e.g., thermal stratification and seasonal nutrient availability) potentially altering RS chlorophyll detection by neural networks-based models, despite their alleged global validity.

Graphical Abstract

1. Introduction

It is widely shared by scientists and public opinion that freshwater ecosystems are of high importance for human wellbeing, although diverse activities act as environmental stressors detrimentally impacting the aquatic habitat quality, functionality, and status [1]. Major pressures such as nutrient pollution, deforestation, and urbanization are strongly reducing water quality, leading to the degradation of freshwater habitats [2,3,4]. In the EU territory, preserving and monitoring water quality has raised growing concern, leading to the adoption of the Water Framework Directive (Directive 2000/60/EC or WFD). WFD encourages policymakers, environmental managers, and researchers to plan and adopt yearly monitoring programmes to improve the ecological quality of all waters within the territory. Phytoplankton biomass monitoring is a conventional method, also included in the WFD guidelines, for assessing water quality status in both inland [5] and marine waters [6]. Specifically, since chlorophyll-a (Chl-a) is the main pigment in photosynthetic organisms, it represents a suitable proxy for phytoplankton biomass estimation [7], thus relating to nutrient concentration and to freshwater trophic state.
However, the need for consistent, long-term assessment of freshwater ecological status urges integrative, frequent, and synoptic approaches. While in situ sampling of photosynthetic pigments is expensive and labour-intensive, as well as spatially and temporally limited, Remote Sensing (RS) provides new opportunities to integrate and improve traditional monitoring techniques [8,9,10]. Further, the increasing number of orbiting satellites makes remote sensing a cost-effective, energy- and time-saving method without equal. In particular, the recent launch of high spatial resolution satellites providing free and open data (e.g., Landsat-8 and Sentinel-2) has enhanced remote freshwater investigations, resulting in an ever-growing number of studies [11,12,13,14,15] with the specific features of the MultiSpectral Instrument (MSI) onboard Sentinel-2 satellites offering additional advantages for water quality measurement [16]. The MSI sensor measures Earth’s reflected radiance in 13 spectral bands. Given its high spatial resolution (10–60 m) in the visible region, and its 12-bit radiometric resolution, MSI allows fine-scale measurement of water properties, enabling the monitoring of small water bodies [13,17]. Furthermore, Sentinel-2 temporal resolution of 2–5 days provides frequent data acquisition, allowing constraints to be overcome (such as frequent cloud cover) and to detect rapid changes of phytoplankton biomass [18].
Remote sensing-based detection of Chl relies on the application of specific indices depending on the biogeochemical characteristics of the water body, which leads to the distinction between case-I and case-II waters [19,20]. The former are waters where phytoplankton and its degradation products dominate optical properties (e.g., open oceans), while inland waters are categorised as case-II waters, being subject to potentially large and independent variations of Chl and other optically active suspended substances. Given the wide regional and temporal variability in constituents of case-II waters, a large number of Chl-retrieving algorithms have been developed and tested so far (e.g., semi-empirical methods or semi-analytical approaches) [15,17,21,22,23,24,25,26,27,28,29], making the selection of the most suitable index sometimes intricate [27,30,31]. Indeed, specific algorithms are often optimized and validated for commonly understood but not always straightforward to define water types, e.g., turbid [32] or clear water [33]. Other statistical methods like neural networks (NN) based on spectral inversion algorithms for radiative transfer simulations [34] have successfully been applied, is intended for use under different bio-optical regimes and without a priori knowledge of the in-water optical conditions (being trained to varying parameter concentration and optical property ranges), sometimes outperforming other methods [27,31,35,36,37,38].
Among these, the Case-2 Regional CoastColour processor (C2RCC) [39] is a “user-friendly” method using MSI data that is readily available in the ESA’s SentiNel Application Platform (SNAP) [28]. C2RCC functioning relies on a set of neural networks deriving water inherent optical properties, which are then converted into Chl-a and total suspended matter (TSM) concentrations. The C2RCC processor is an improvement of the Case 2 Regional processor originally developed by Doerffer and Schiller [34], integrating an atmospheric correction part dedicated for water observation, and has been trained to cover even extreme ranges of scattering and absorption [39]. On the other hand, the Mixture Density Network (MDN) is a recently developed algorithm by Pahlevan et al. [27], which is a variation of Multilayer Perceptrons NN methods and relies on a class of neural networks that estimates multimodal Gaussian distributions over a range of solutions. The MDN algorithm has been trained with 1000 co-located in situ Rrs- Chl-a pairs and was initially tested with Sentinel-2 and Sentinel-3 data, outperforming previously published algorithms (e.g., OC and Blend) [27]. Still, both C2RCC and MDN are far from being considered universally applicable approaches for Chl retrieval, and their efficiency needs to be tested under different conditions.
Water quality analysis of volcanic lakes has rarely been addressed by applying RS techniques [40,41] and, as far as we know, none of these studies has investigated warm monomictic volcanic lakes using MSI data. In such lakes, temperatures never drop below 4 °C and thermal stratification occurs during summer, when air temperature increases, and surface water is heated more rapidly. Water circulation occurs only once a year when surface waters cool to a temperature equal to the bottom waters. In most temperate and sub-tropical warm monomictic lakes, water column mixing brings hypolimnetic nutrients to the surface, leading to annual maximum phytoplankton proliferation [42,43,44]. Such a marked variation in primary production between the thermally stratified season and the full water column mixing could define a sharp temporal divergence in the viability of RS observation of such water bodies.
Accordingly, the purpose of this paper is to assess the potential of two open-source algorithms (i.e., C2RCC and MDN) applied to MSI Sentinel-2 data to retrieve Chl content in three monomictic volcanic lakes (Central Italy) while accounting for their water circulation pattern. Given the appropriateness of both processors for a broad range of bio-optical conditions, we hypothesize that the water circulation pattern of such lakes (i.e., one mixing and one stratification phase) will not affect the strength of the relationship between remote sensing and field data. Due to growing anthropogenic activity and to an increasing frequency of algal blooms, the need for trophic level assessment of Italian volcanic lakes has become essential for their preservation [45].

2. Materials and Methods

2.1. Study Area

This study was carried out in three warm monomictic volcanic lakes located in the Latium region, Italy: Lake Albano, Lake Bracciano, and Lake Nemi (Figure 1).
These lakes are all located in craters originated from Quaternary volcanic activity [46], which determines their circular shape and great depth compared to the relatively low surface area [44] (further described in Appendix A). As their water input is nearly exclusively provided by underground springs and rainfall, these lakes are characterized by long water renewal times, which can represent an additional risk factor if deterioration phenomena occur [45,47].
Based on the temperature data of the water column collected by the Regional Agency for Environmental Protection of Lazio (ARPA Lazio), in all three lakes, the seasonal thermal stratification for the year 2019 started in spring, in agreement with findings from previous studies reporting similar trends [42,45]. According to the water circulation phases, chlorophyll concentration data were grouped as follows (Figure 2): mixing from November to April, stratification from May to October.

2.2. Digital Data Collection and Processing

Sentinel-2 MSI L1C products were downloaded from the Copernicus Open Access Hub (https://scihub.copernicus.eu (accessed on 28 October 2019)). L1C products contain orthorectified, georeferenced Top of Atmosphere (TOA) reflectance in UTM map projection with WGS84 datum. Matchups with in situ measurements were selected with a mean time difference between suitable Sentinel-2 and in situ data of ~2 days. For only one matchup a maximum time difference of 8 days with Sentinel-2 MSI overpass was allowed (Table 1), provided that weather conditions were stable between the image acquisition and field sampling and that no algal blooms occurred. This time window was set to allow sufficient matchups between ground data and Sentinel-2 imagery.
Satellite images were visually checked to make sure that the studied lakes were not covered by clouds or cloud shadows before application by exploiting the spectral band at 1375 nm (band 10), which enables cirrus detection. Due to cloud cover on images, no satellite data are available for two sampling dates for Lake Albano and Lake Nemi.
Two different indices were implemented, one calculated through the neural networks of the C2RCC processor (INN) and one calculated through the MDN algorithm (IMDN). The Case 2 Regional Coast Colour (C2RCC) processor (v1.0) is readily available on SNAP among the Thematic Water Processing tools. Getting as input the water leaving reflectances from the atmosphere part, the in-water part performs the inversion of the spectrum into inherent optical properties (IOPs), i.e., absorption and backscattering coefficients, from which Chl and total suspended matter concentrations are calculated. The simulations on which the C2RCC is built are based on a large set of optical data and concentration measurements, enabling the processor to cover extreme ranges of scattering and absorption and making it suitable for a large variety of study areas [34]. To better adapt the processing to the conditions of the study area, C2RCC allows flexibility in adjusting ancillary parameters. Among the processing parameters, we adjusted elevation, salinity, and temperature according to local values registered during each sampling campaign. The CHL exponent and the CHL factor required to better estimate Chl concentration remained unchanged from the default values, as it was not possible to retrieve such parameters from the lakes under study. To test the suitability of this approximation, we assessed the relationship between in situ measured Chl and Chl absorption coefficient (iop_apig) calculated by the C2RCC (the results are presented in Appendix B).
The MDN algorithm was recently developed by Pahlevan et al. [27] and is freely available on GitHub (https://github.com/STREAM-RS/STREAM-RS (accessed on 26 May 2021)). Such processor models conditional probability distributions of the target variables (i.e., Chl-a) as a mixture of multiple Gaussian functions from remote sensing reflectance (Rrs). The Gaussians are subsequently combined to form the final output estimation through maximum likelihood [48]. In the case of the MSI sensor, the processor requires as input Rrs values at the following wavelengths: 443 nm, 490 nm, 560 nm, 665 nm, 705 nm, 740 nm, and 783 nm (i.e., Band 1 through Band 7). In the present study, Rrs reflectance values were calculated via the C2RCC atmospheric correction.
The SNAP Pixel Extraction tool was later used to extract pixel values from the GeoTIFF products obtained from both types of processing. To minimize errors due to boat drift and GPS positional inaccuracies, the median value of a 3 × 3 pixel window centred at each sampling point was extracted and used for further statistical analysis.

2.3. In Situ Data

To validate remote sensing data, a total of 30 sampling campaigns were conducted over the period March-October 2019. Samples were collected at five sampling sites distributed uniformly across each lake. The sampling site coordinates were recorded, and geolocation uncertainties related to boat drift and GPS positional inaccuracies were taken into account in the satellite data analysis (3 × 3 pixel windows, see Section 2.2). Concurrently, epilimnetic water temperature, pH, salinity, and dissolved oxygen (DO) concentration were recorded at each sampling site using a Hach HQ40d Portable Multi-Parameter Meter.
To estimate water Chl content, 2 L water samples, one at each sampling site, were collected from the surface layer (max. 50 cm depth, according to previous literature) with a Van Dorn water sampler [13,21,26,27,28,29,49]. Samples were stored on ice in sterile, light-reflective containers (to prevent pigment degradation).
Within 24 h of collection, water samples underwent laboratory measurement of their Chl content. In this study, total Chl concentration measurements were carried out using the Phyto-PAM—Multiple Excitation Wavelength Phytoplankton & Photosynthesis Analyzer (Heinz Walz GmbH, Effeltrich, Germany). This instrument allows non-invasive and rapid biomass estimates since it does not require sample treatments before carrying out the analysis. The functioning of this fluorimeter relies on the selective amplification of a fluorescence signal, measured with the help of short (μsec) pulses of measuring light at a high repetition rate (1200 Hz). These light pulses are provided by an array of light-emitting diodes at four distinct light wavelengths (470, 520, 645, and 665 nm). A photomultiplier in conjunction with a pulse amplifier is used as a detector for fluorescence at wavelengths above 710 nm. The information consisting of four independent fluorescence signals is further processed by computer-controlled data analysis using the dedicated PhytoWin-software. Three measurements were carried out for each water sample. A Zero offset (Zoff) determination preceded fluorescence measurement to remove other fluorescing substances’ contribution to the measured signal (e.g., humic acids).
As remote sensing indices are developed to retrieve Chl-a concentration and the Phyto-PAM instrument measures total Chl concentration, a standard spectrophotometric method was applied to a representative subset of water samples to validate the use of the Phyto-PAM results (further description of the validation process can be found in Appendix C). We assumed that Phyto-PAM’s concentration values were consistent with Chl-a concentration. Water samples were filtered through Whatman GF/F glass microfiber filters (0.7 μm pore size) under vacuum pressure. Filters were subsequently put in 5 mL of 90% acetone solution and then homogenised, and acetone solution was added to the emulsion to reach the final volume of 10 mL. Sample extracts were stored at 4 °C in the dark for 24 h for pigment extraction. Extracts were clarified by centrifugation (3500 rpm, 12 min) and then transferred to a 15-mL vial to measure their final volumes. Triplicate absorbance measurements were carried out on each water sample using a dual-beam spectrophotometer. The resulting absorbance values at 630 nm, 647 nm, 664 nm, and 750 nm were entered into the equation by Jeffrey & Humphrey (1975) [50] to estimate Chl-a concentration.

2.4. Statistical Analysis

To test the relationship between remote sensing derived indices and in situ Chl concentration, we assumed a positive linear relationship between the two sets of data. Due to the spatial interdependence of the samples in our dataset, we used linear mixed models (LMM) as they allowed us to use the lake identity as a random intercept effect. We used R software v. 3.6.1 [51] and lme4 package [52] to fit the linear mixed-effect model, using the indices as the dependent variable and the in situ Chl content as a fixed-effect independent variable. Further, we considered the interaction of in situ Chl concentration with the water mixing period (i.e., mixing from November to April, stratification from May to October), deleting the intercept to see the effect of the interaction. Finally, to determine the amount of variance in our in situ data explained by the models, we considered the model’s conditional R2 (R2c), related to the variance explained by both fixed and random effects, and marginal R2 (R2m), related to the variance explained only by fixed effects [53]. No evident deviations from normality emerged after visual inspection of residual plots.

3. Results

3.1. In Situ Data

Descriptive statistical results of measured parameters are shown in Table 2. Over the sampling period, Chl and Chl-a concentrations exhibited comparable trends in all three lakes, with the highest values recorded in early spring and minimum values in July (Table 2). Epilimnetic water temperatures showed their lowest values in March and maximum values in July, followed by a slow decrease. Comparable values of DO and pH were observed among the three lakes, whereas salinity was significantly lower in lake Nemi compared to the other two lakes.

3.2. Indices Performance Assessment

When analysing the results for the INN and IMDN models, we found significant dependency on the water mixing period. In particular, the relationship between INN and the in situ data was significant when considering data in the mixing phase, while there was no significant relationship during the stratification period (Figure 3a). For the mixing phase, there was a positive relationship between the two sets of data, while for the stratification phase the relationship was negative, although not significant (Figure 3b).
With the IMDN index, similar outcomes were obtained in comparison to the model with INN. The relationship between Chl data and remote sensing data was positive and significant during the water mixing phase, whereas there was no significant relationship at low Chl concentrations, that is, during the stratification phase (Figure 4).
Furthermore, the variance explained by the two indices showed comparable values of R2m (INN = 0.331, IMDN = 0.316) and R2c (INN = 0.475, IMDN = 0.467), indicating that the neural networks-based approaches implemented in C2RCC and MDN are able to explain a satisfactory portion of the variability in the in situ data.

4. Discussion

The aim of this paper was to test the performance of the C2RCC and MDN processors applied to Sentinel-2 products for Chl concentration monitoring in three volcanic lakes of central Italy while assessing the influence of the water mixing conditions on the relationship between remote sensing and field data. We chose to apply (a) the C2RCC processor due to its ease of operation within the ESA’s SentiNel Application Platform (SNAP) and wide applicability, thus potentially making it a favoured choice for many users; (b) the MDN algorithm as it is a recently developed method that proved to outperform previous, well-established methods when applied to Sentinel-2 data. Given the purpose of the present study, our efforts were focused on a restricted (but representative) number of comparable lakes over one seasonal cycle, unlike other studies investigating one [12,29] or several water bodies [15] on a single or limited dates.
This study showed that MSI data has potential for Chl estimation in monomictic volcanic lakes, proving Sentinel-2 satellites to be powerful tools for water quality monitoring. Indeed, the two models based on INN and IMDN both exhibited a positive and significant relationship between Chl concentration and satellite data. Nonetheless, contrary to our expectations, such a significant relationship was observed only during the mixing phase, indicating that the water circulation pattern of such lakes (i.e., one mixing and one stratification phase) may affect the strength of the relationship between remote sensing and field data when using C2RCC and MDN processors.
Seemingly, the significant effect of thermal stratification on the performance of the RS-based index can be explained by the low responsiveness of the RS approach applied here to the low amount of phytoplankton biomass. Indeed, the large amount of nutrients made available by the full water circulation taking place during the mixing phase increases phytoplankton biomass, which can be more easily captured by RS techniques. Conversely, when the stratification phase is established, the water column is divided into three compartments: (1) an upper layer (epilimnion) of warm, oxygenated water, where the majority of primary production takes place; (2) a thin middle layer (metalimnion) where temperature changes more rapidly; and (3) a cold, low in dissolved oxygen layer (hypolimnion) where the decomposition of sedimenting organisms from the epilimnion occurs [43,44]. Hence, during thermal stratification, the organisms’ metabolism within the epilimnion causes a progressive reduction of nutrient availability with a consequent phytoplankton mortality increase in this layer [54,55,56,57,58]. Accordingly, when considering data related to the stratification phase, the absence of a significant relationship between Chl concentration surveyed by the Phyto-PAM and the one retrieved by RS can have a threefold explanation.
First, several studies proved the important (and sometimes dominant) contribution that inorganic and organic sediments have on the optical properties of case-II waters [59,60,61,62,63,64]. Their effects on total backscattering are magnified when phytoplankton biomass is low. In this case, the remotely-sensed signal can be noticeably impacted by non-phytoplankton water constituents (e.g., minerals and CDOM) which are not correlated with Chl-a concentration [34,62,65,66]. In the case of spectral inversion algorithms, as in C2RCC, the ability of radiative transfer simulations to simultaneously solve several IOPs (i.e., absorption and backscattering coefficients) is sensitive to errors due to the existence of ambiguous solutions resulting from the additive nature of the IOPs. That is to say, different combinations of IOPs of each water component can lead to an identical sum of IOPs and, thus, to similar reflectance values [67,68]. On the contrary, the MDN algorithm attempts to overcome this one-to-many issue by modeling a conditional probability distribution of the target variables given input Rrs data rather than directly modeling the Rrs- Chl-a relationship like standard NNs [27,48,69]. However, the results obtained in the present study show a comparable performance, indicating that this non-uniqueness problem is yet to be fully solved for the investigated lakes. Indeed, such ambiguity error on the total absorption coefficient is the highest for highly absorbing water bodies [67]. Therefore, even though the methods used in this study are suitable for phytoplankton monitoring, the difficulty of sensing low Chl-a signals limits their effectiveness only to the mixing phase of the yearly water circulation pattern of this type of lakes. Moreover, the performance of this kind of approach is dependent on the ranges of the training parameters. Indeed, when the bio-optical properties of the investigated area are beyond training ranges, this is expected to adversely impact the NN efficiency of retrieving Chl concentration [70,71].
Second, the absence of a significant correspondence between remotely-sensed data and field data could be concurrently ascribed to the MSI sensor specifications in terms of Signal-to-Noise Ratio (SNR). Sentinel-2 MSI SNR (generally below 150:1 for most of its spectral bands) is much lower than the minimum SNR required for Chl retrievals in marine waters, which is set at 400:1 for all visible bands and at 600:1 for NIR bands [72,73]. However, improving SNR values when designing sensors requires trade-offs in terms of reduced spatial and spectral resolution, as well as other radiometric characteristics [72,74]. Due to the relatively low Sentinel-2 SNR, it is reasonable to assume that the MSI sensor mainly measures noise when Chl concentration is low. Consequently, a higher relative error in the sensed signal is expected [75], which affects the significance of the indices during the stratification phase.

5. Conclusions

Sentinel-2 data proved their valuable potential for remote sensing of phytoplankton in small inland water bodies, otherwise challenging with previous sensors. This study confirmed the feasibility of C2RCC and MDN applied to Sentinel 2 data to monitor chlorophyll content in warm monomictic lakes. The approach presented has a good potential for operational monitoring of small and medium-sized water bodies required, for instance, by the European Union Water Framework Directive. Nonetheless, it is worth highlighting the limits intrinsically embedded in our dataset. Indeed, this study was carried out on a group of peculiar lakes and over just one year, without an optimal (although acceptable) time match-up between ground data and Sentinel-2 imagery. Additionally, the absence of simultaneous field radiometry data did not allow an accurate validation of the atmospheric correction.
Caution should be taken when relying on remote sensing for water quality monitoring (i.e., C2RCC and MDN) because the consistency of such a method may depend on the hydrological and ecological processes associated with the water mixing conditions, which have not been addressed in the algorithm yet. Especially the low chlorophyll concentrations during the stratification phase can prove challenging for the performance of the tested approach, despite its claimed broad applicability. From the present study, it seems evident that preliminary ecological knowledge of the water body under investigation is of paramount importance for a proper application of the Chl-retrieving algorithms, even on a temporal scale. Therefore, indiscriminate use of the neural networks-based approaches assessed in this study without focusing on the local biological and hydrodynamic processes could lead to inconsistent results, because water thermal stratification of lakes could significantly affect retrieval.

Author Contributions

Conceptualization, M.P., M.S. (Massimiliano Scalici) and M.M.; methodology, D.M., J.K., M.S. (Maria Sighicelli) and M.M.; software, D.M. and J.K.; validation, M.P., L.C. and M.M.; formal analysis, M.P., L.C. and M.M.; investigation, M.P., M.S. (Maria Sighicelli) and F.L.; resources, M.S. (Massimiliano Scalici) and M.M.; data curation, M.P., L.C. and M.M.; writing—original draft preparation, M.P.; writing—review and editing, M.P., M.S. (Massimiliano Scalici), L.C., D.M., J.K., M.S. (Maria Sighicelli) and M.M.; visualization, M.P., L.C. and M.M.; supervision, M.S. (Massimiliano Scalici) and M.M.; project administration, M.S. (Massimiliano Scalici) and M.M.; funding acquisition, M.S. (Massimiliano Scalici) and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the grant of Excellence Departments, MIUR-Italy (ARTICOLO 1 COMMI 314-337, LEGGE 232/2016).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to private restriction. Authors agree on considering data of own property, that can be shared with other research groups on request.

Acknowledgments

Authors would like to acknowledge the data support from the Regional Agency for Environmental Protection of Lazio (ARPA Lazio).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Lake Bracciano is included in Bracciano-Martignano Regional Park and designated as a Site of Community Importance (SCI-IT6030010). Lake Bracciano has recently been used as a drinking water supply in the event of a city water emergency. The lake water crisis of 2017 caused by dry weather and constant uptake led to a substantial water-level reduction (up to −198 cm) and to water quality deterioration, causing growing concern for the preservation of the ecosystem integrity.
Lake Albano is ranked as meso-eutrophic [47,76], with progressive deterioration over time due to an increase of sewage discharge [77] and over-extraction of groundwater within the watershed. Cyanotoxic algal blooms regularly occur in Lake Albano, when the water column mixing brings hypolimnetic nutrients to the surface [76,78,79]. Even though its quality is gradually deteriorating, Lake Albano is designated SCI (IT6030038).
Together with Lake Albano, Lake Nemi is part of Castelli Romani Regional Park. Recent data reveal an improvement in Lake Nemi’s water quality. After the implementation of domestic sewage diversion in the early 1990s, all biological and chemical data showed partial recovery, shifting Lake Nemi’s trophic status from eutro-hypereutrophic to meso-eutrophic [46,47,77].
Table A1. Main features of the volcanic lakes under investigation.
Table A1. Main features of the volcanic lakes under investigation.
Lake BraccianoLake AlbanoLake Nemi
Location (Lat., Lon.)42°07′16″N 12°13′55″E41°45′0″N 12°39′54″E41°42′44″N 12°42′09″E
Max. depth (m)16517527.5
Mean elevation (m a.s.l.)164293316
Surface area (km2)57.56.01.6
Volume (106 m3)505046426.5
Renewal time (yr)13747.615
OutflowsArrone river (currently dry in its first stretch), Paul aqueductNo natural outletsNo natural outlets

Appendix B

Table A2. Marginal R2 (R2m) and conditional R2 (R2c) of the linear mixed model based on the chlorophyll absorption coefficients retrieved from the C2RCC (iop_apig).
Table A2. Marginal R2 (R2m) and conditional R2 (R2c) of the linear mixed model based on the chlorophyll absorption coefficients retrieved from the C2RCC (iop_apig).
R2mR2c
0.3260.476
Figure A1. (a) Coefficient plot of the linear mixed model designed with in situ measured Chl and Chl absorption coefficient (apig) calculated by the C2RCC. Dots represent coefficient estimates; the outer error bars are 95% confidence intervals and the inner error bars are 50% confidence intervals. The 95% coefficient interval for stratification includes zero, which indicates a non-significant relationship between measured chlorophyll content and pigment absorption coefficient estimates. (b) Plot of the interaction effect between apig and water circulation phase.
Figure A1. (a) Coefficient plot of the linear mixed model designed with in situ measured Chl and Chl absorption coefficient (apig) calculated by the C2RCC. Dots represent coefficient estimates; the outer error bars are 95% confidence intervals and the inner error bars are 50% confidence intervals. The 95% coefficient interval for stratification includes zero, which indicates a non-significant relationship between measured chlorophyll content and pigment absorption coefficient estimates. (b) Plot of the interaction effect between apig and water circulation phase.
Remotesensing 13 02699 g0a1

Appendix C

A correlation test between Phyto-PAM chlorophyll values and the ones obtained spectrophotometrically was performed. Before carrying out this test, it has been necessary to assess the effective relationship between total chlorophyll and chlorophyll-a, as well as testing whether this latter pigment was the major contributor to total Chl concentration. To this end, a T-test was used to calculate the differences between Chl-a and total chlorophyll concentration. The results indicate that Chl-a does not differ significantly from total chlorophyll (p-value = 0.402). The association between total Chl and Chl-a values was assessed using a correlation test, whose result was significant (p-value < 2.22 × 10−16), with r = 0.99. Based on the aforementioned results, total chlorophyll and Chl-a concentrations measured spectrophotometrically showed comparable values and trends. The result of the correlation test exhibits a significant and high association (p-value < 2.22 × 10−16, r = 0.95) between total Chl estimated by Phyto-PAM and through spectrophotometric measurements.

References

  1. Vörösmarty, C.J.; McIntyre, P.B.; Gessner, M.O.; Dudgeon, D.; Prusevich, A.; Green, P.; Glidden, S.; Bunn, S.E.; Sullivan, C.A.; Liermann, C.R.; et al. Global threats to human water security and river biodiversity. Nature 2010, 467, 555–561. [Google Scholar] [CrossRef]
  2. Brönmark, C.; Hansson, L.A. Environmental issues in lakes and ponds: Current state and perspectives. Environ. Conserv. 2002, 29, 290–307. [Google Scholar] [CrossRef] [Green Version]
  3. Duker, L.; Borre, L. Biodiversity Conservation of the World’s Lakes: A Preliminary Framework for Identifying Priorities. Available online: www.worldlakes.org/uploads/report2.pdf (accessed on 5 June 2021).
  4. Dudgeon, D.; Arthington, A.H.; Gessner, M.O.; Kawabata, Z.I.; Knowler, D.J.; Lévêque, C.; Naiman, R.J.; Prieur-Richard, A.H.; Soto, D.; Stiassny, M.L.J.; et al. Freshwater biodiversity: Importance, threats, status and conservation challenges. Biol. Rev. Camb. Philos. Soc. 2006, 81, 163–182. [Google Scholar] [CrossRef]
  5. Padisák, J.; Borics, G.; Grigorszky, I.; Soróczki-Pintér, É. Use of phytoplankton assemblages for monitoring ecological status of lakes within the water framework directive: The assemblage index. Hydrobiologia 2006, 553, 1–14. [Google Scholar] [CrossRef]
  6. Garmendia, M.; Borja, Á.; Franco, J.; Revilla, M. Phytoplankton composition indicators for the assessment of eutrophication in marine waters: Present state and challenges within the European directives. Mar. Pollut. Bull. 2013, 66, 7–16. [Google Scholar] [CrossRef] [PubMed]
  7. Brönmark, C.; Hansson, L.-A. The Biology of Lakes and Ponds; Oxford University Press: Oxford, UK, 2017. [Google Scholar]
  8. Chen, Q.; Zhang, Y.; Ekroos, A.; Hallikainen, M. The role of remote sensing technology in the EU water framework directive (WFD). Environ. Sci. Policy 2004, 7, 267–276. [Google Scholar] [CrossRef]
  9. Williamson, C.E.; Saros, J.E.; Vincent, W.F.; Smol, J.P. Lakes and reservoirs as sentinels, integrators, and regulators of climate change. Limnol. Oceanogr. 2009, 54, 2273–2282. [Google Scholar] [CrossRef]
  10. Palmer, S.C.J.; Kutser, T.; Hunter, P.D. Remote sensing of inland waters: Challenges, progress and future directions. Remote Sens. Environ. 2015, 157, 1–8. [Google Scholar] [CrossRef] [Green Version]
  11. Bresciani, M.; Cazzaniga, I.; Austoni, M.; Sforzi, T.; Buzzi, F.; Morabito, G.; Giardino, C. Mapping phytoplankton blooms in deep subalpine lakes from Sentinel-2A and Landsat-8. Hydrobiologia 2018, 824, 197–214. [Google Scholar] [CrossRef] [Green Version]
  12. Dörnhöfer, K.; Göritz, A.; Gege, P.; Pflug, B.; Oppelt, N. Water constituents andwater depth retrieval from Sentinel-2A-A first evaluation in an oligotrophic lake. Remote Sens. 2016, 8, 941. [Google Scholar] [CrossRef] [Green Version]
  13. Grendaitė, D.; Stonevičius, E.; Karosienė, J.; Savadova, K.; Kasperovičienė, J. Chlorophyll-a concentration retrieval in eutrophic lakes in Lithuania from Sentinel-2 data. Geol. Geogr. 2018, 4. [Google Scholar] [CrossRef]
  14. Liu, H.; Li, Q.; Shi, T.; Hu, S.; Wu, G.; Zhou, Q. Application of sentinel 2 MSI images to retrieve suspended particulate matter concentrations in Poyang Lake. Remote Sens. 2017, 9, 761. [Google Scholar] [CrossRef] [Green Version]
  15. Toming, K.; Kutser, T.; Laas, A.; Sepp, M.; Paavel, B.; Nõges, T. First experiences in mapping lakewater quality parameters with sentinel-2 MSI imagery. Remote Sens. 2016, 8, 640. [Google Scholar] [CrossRef] [Green Version]
  16. D’Odorico, P.; Gonsamo, A.; Damm, A.; Schaepman, M.E. Experimental evaluation of sentinel-2 spectral response functions for NDVI time-series continuity. IEEE Trans. Geosci. Remote Sens. 2013, 51, 1336–1348. [Google Scholar] [CrossRef]
  17. Kutser, T.; Paavel, B.; Verpoorter, C.; Ligi, M.; Soomets, T.; Toming, K.; Casal, G. Remote sensing of black lakes and using 810 nm reflectance peak for retrieving water quality parameters of optically complex waters. Remote Sens. 2016, 8, 497. [Google Scholar] [CrossRef]
  18. Beck, R.; Zhan, S.; Liu, H.; Tong, S.; Yang, B.; Xu, M.; Ye, Z.; Huang, Y.; Shu, S.; Wu, Q.; et al. Comparison of satellite reflectance algorithms for estimating chlorophyll-a in a temperate reservoir using coincident hyperspectral aircraft imagery and dense coincident surface observations. Remote Sens. Environ. 2016, 178, 15–30. [Google Scholar] [CrossRef] [Green Version]
  19. Gordon, H.R.; Morel, A.Y. Remote assessment of ocean color for interpretation of satellite visible imagery a review. In Lecture Notes on Coastal and Estuarine Studies; Springer: New York, NY, USA, 1983; ISBN 0387909230. [Google Scholar]
  20. Morel, A. Optical modeling of the upper ocean in relation to its biogenous matter content (case I waters). J. Geophys. Res. 1988, 93, 10749. [Google Scholar] [CrossRef] [Green Version]
  21. Ha, N.T.T.; Thao, N.T.P.; Koike, K.; Nhuan, M.T. Selecting the best band ratio to estimate chlorophyll-a concentration in a tropical freshwater lake using sentinel 2A images from a case study of Lake Ba Be (Northern Vietnam). ISPRS Int. J. Geo-Inf. 2017, 6, 290. [Google Scholar] [CrossRef]
  22. Bonansea, M.; Rodriguez, M.C.; Pinotti, L.; Ferrero, S. Using multi-temporal Landsat imagery and linear mixed models for assessing water quality parameters in Río Tercero reservoir (Argentina). Remote Sens. Environ. 2015, 158, 28–41. [Google Scholar] [CrossRef]
  23. Bresciani, M.; Stroppiana, D.; Odermatt, D.; Morabito, G.; Giardino, C. Assessing remotely sensed chlorophyll-a for the implementation of the Water Framework Directive in European perialpine lakes. Sci. Total Environ. 2011, 409, 3083–3091. [Google Scholar] [CrossRef] [Green Version]
  24. Giardino, C.; Bresciani, M.; Villa, P.; Martinelli, A. Application of remote sensing in water resource management: The case study of Lake Trasimeno, Italy. Water Resour. Manag. 2010, 24, 3885–3899. [Google Scholar] [CrossRef]
  25. Mishra, S.; Mishra, D.R. Normalized difference chlorophyll index: A novel model for remote estimation of chlorophyll-a concentration in turbid productive waters. Remote Sens. Environ. 2012, 117, 394–406. [Google Scholar] [CrossRef]
  26. Hafeez, S.; Wong, M.; Ho, H.; Nazeer, M.; Nichol, J.; Abbas, S.; Tang, D.; Lee, K.; Pun, L. Comparison of machine learning algorithms for retrieval of water quality indicators in case-II waters: A case study of Hong Kong. Remote Sens. 2019, 11, 617. [Google Scholar] [CrossRef] [Green Version]
  27. Pahlevan, N.; Smith, B.; Schalles, J.; Binding, C.; Cao, Z.; Ma, R.; Alikas, K.; Kangro, K.; Gurlin, D.; Hà, N.; et al. Seamless retrievals of chlorophyll-a from Sentinel-2 (MSI) and Sentinel-3 (OLCI) in inland and coastal waters: A machine-learning approach. Remote Sens. Environ. 2020, 240, 111604. [Google Scholar] [CrossRef]
  28. Ansper, A.; Alikas, K. Retrieval of chlorophyll a from Sentinel-2 MSI data for the European Union water framework directive reporting purposes. Remote Sens. 2019, 11, 64. [Google Scholar] [CrossRef] [Green Version]
  29. Dörnhöfer, K.; Klinger, P.; Heege, T.; Oppelt, N. Multi-sensor satellite and in situ monitoring of phytoplankton development in a eutrophic-mesotrophic lake. Sci. Total Environ. 2018, 612, 1200–1214. [Google Scholar] [CrossRef] [PubMed]
  30. Odermatt, D.; Gitelson, A.; Brando, V.E.; Schaepman, M. Review of constituent retrieval in optically deep and complex waters from satellite imagery. Remote Sens. Environ. 2012, 118, 116–126. [Google Scholar] [CrossRef] [Green Version]
  31. Palmer, S.C.J.; Odermatt, D.; Hunter, P.D.; Brockmann, C.; Présing, M.; Balzter, H.; Tóth, V.R. Satellite remote sensing of phytoplankton phenology in Lake Balaton using 10years of MERIS observations. Remote Sens. Environ. 2015, 158, 441–452. [Google Scholar] [CrossRef] [Green Version]
  32. Gitelson, A.A.; Schalles, J.F.; Hladik, C.M. Remote chlorophyll-a retrieval in turbid, productive estuaries: Chesapeake Bay case study. Remote Sens. Environ. 2007, 109, 464–472. [Google Scholar] [CrossRef]
  33. Belzile, C.; Vincent, W.F.; Howard-Williams, C.; Hawes, I.; James, M.R.; Kumagai, M.; Roesler, C.S. Relationships between spectral optical properties and optically active substances in a clear oligotrophic lake. Water Resour. Res. 2004, 40, 1–12. [Google Scholar] [CrossRef] [Green Version]
  34. Doerffer, R.; Schiller, H. The MERIS case 2 water algorithm. Int. J. Remote Sens. 2007, 28, 517–535. [Google Scholar] [CrossRef]
  35. Schroeder, T.; Schaale, M.; Fischer, J. Retrieval of atmospheric and oceanic properties from MERIS measurements: A new Case-2 water processor for BEAM. Int. J. Remote Sens. 2007, 28, 5627–5632. [Google Scholar] [CrossRef]
  36. Odermatt, D.; Giardino, C.; Heege, T. Chlorophyll retrieval with MERIS Case-2-Regional in perialpine lakes. Remote Sens. Environ. 2010, 114, 607–617. [Google Scholar] [CrossRef] [Green Version]
  37. Cui, T.; Zhang, J.; Groom, S.; Sun, L.; Smyth, T.; Sathyendranath, S. Validation of MERIS ocean-color products in the Bohai Sea: A case study for turbid coastal waters. Remote Sens. Environ. 2010, 114, 2326–2336. [Google Scholar] [CrossRef]
  38. Minghelli-Roman, A.; Laugier, T.; Polidori, L.; Mathieu, S.; Loubersac, L.; Gouton, P. Satellite survey of seasonal trophic status and occasional anoxic “malaïgue” crises in the Thau lagoon using MERIS images. Int. J. Remote Sens. 2011, 32, 909–923. [Google Scholar] [CrossRef] [Green Version]
  39. Brockmann, C.; Doerffer, R.; Peters, M.; Kerstin, S.; Embacher, S.; Ruescas, A. Evolution of the C2RCC neural network for Sentinel 2 and 3 for the retrieval of ocean colour products in normal and extreme optically complex waters. In Proceedings of the ESASP, Prague, Czech Republic, 9–13 May 2016. [Google Scholar]
  40. Mancino, G.; Nolè, A.; Urbano, V.; Amato, M.; Ferrara, A. Assessing water quality by remote sensing in small lakes: The case study of Monticchio lakes in southern Italy. IForest 2009, 2, 154–161. [Google Scholar] [CrossRef] [Green Version]
  41. Oppenheimer, C. Remote sensing of the colour and temperature of volcanic lakes. Int. J. Remote Sens. 1997, 18, 5–37. [Google Scholar] [CrossRef]
  42. Ellwood, N.T.W.; Albertano, P.; Galvez, R.; Funiciello, R.; Mosello, R. Water chemistry and trophic evaluation of Lake Albano (Central Italy): A four year water monitoring study. J. Limnol. 2009, 68, 288–303. [Google Scholar] [CrossRef]
  43. Wetzel, R.G. Limnology: Lake and River Ecosystems; Elsevier Academic Press, 2001; ISBN 9780127447605. [Google Scholar]
  44. O’Sullivan, P.E.; Reynolds, C.S. The Lakes Handbook: Limnology and Limnetic Ecology; Blackwell Publishing: Oxford, UK, 2004; ISBN 0-632-04797-6. [Google Scholar]
  45. Bruno, M.; Marchiori, E.; Mecozzi, M.; Congestri, R.; Melchiorre, S.; Falleni, F.; Nusca, A. Risanamento Trofico Negli Ecosistemi Lacustri: Confronto Fra i Laghi di Bracciano e Martignano; Rapporti ISTISAN; Istituto Superiore di Sanità: Rome, Italy, 2006. [Google Scholar]
  46. Margaritora, F.G. Limnology in Latium: The volcanic lakes. Mem. Dell’Istituto Ital. Di Idrobiol. 1992, 50, 319–336. [Google Scholar]
  47. Margaritora, F.G.; Bazzanti, M.; Ferrara, O.; Mastrantuono, L.; Seminara, M.; Vagaggini, D. Classification of the ecological status of volcanic lakes in Central Italy. J. Limnol. 2003, 62, 49–59. [Google Scholar] [CrossRef] [Green Version]
  48. Pahlevan, N.; Smith, B.; Binding, C.; Gurlin, D.; Li, L.; Bresciani, M.; Giardino, C. Hyperspectral retrievals of phytoplankton absorption and chlorophyll-a in inland and nearshore coastal waters. Remote Sens. Environ. 2021, 253, 112200. [Google Scholar] [CrossRef]
  49. Binding, C.E.; Greenberg, T.A.; McCullough, G.; Watson, S.B.; Page, E. An analysis of satellite-derived chlorophyll and algal bloom indices on Lake Winnipeg. J. Great Lakes Res. 2018, 44, 436–446. [Google Scholar] [CrossRef]
  50. Jeffrey, S.W.; Humphrey, G.F. New spectrophotometric equations for determining chlorophylls a, b, c1 and c2 in higher plants, algae and natural phytoplankton. Biochem. Physiol. Pflanz. 1975, 167, 191–194. [Google Scholar] [CrossRef]
  51. R Core Team R. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  52. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  53. Nakagawa, S.; Schielzeth, H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 2013, 4, 133–142. [Google Scholar] [CrossRef]
  54. Vincent, W.F.; Gibbs, M.M.; Dryden, S.J. Accelerated eutrophication in a New Zealand lake: Lake rotoiti, central north island. N. Z. J. Mar. Freshw. Res. 1984, 18, 431–440. [Google Scholar] [CrossRef]
  55. Von Westernhagen, N.; Hamilton, D.P.; Pilditch, C.A. Temporal and spatial variations in phytoplankton productivity in surface waters of a warm-temperate, monomictic lake in New Zealand. Hydrobiologia 2010, 652, 57–70. [Google Scholar] [CrossRef] [Green Version]
  56. Vincent, W.F. Phytoplankton production and winter mixing: Contrasting effects in two oligotrophic lakes. J. Ecol. 1983, 71, 1–20. [Google Scholar] [CrossRef]
  57. Hamilton, D.P.; O’Brien, K.R.; Burford, M.A.; Brookes, J.D.; McBride, C.G. Vertical distributions of chlorophyll in deep, warm monomictic lakes. Aquat. Sci. 2010, 72, 295–307. [Google Scholar] [CrossRef] [Green Version]
  58. Flint, E.A. Phytoplankton in seven monomictic lakes near rotorua, New Zealand. N. Z. J. Bot. 1977, 15, 197–208. [Google Scholar] [CrossRef] [Green Version]
  59. Gons, H.J.; Auer, M.T.; Effler, S.W. MERIS satellite chlorophyll mapping of oligotrophic and eutrophic waters in the Laurentian Great Lakes. Remote Sens. Environ. 2008, 112, 4098–4106. [Google Scholar] [CrossRef]
  60. Gitelson, A.A.; Dall’Olmo, G.; Moses, W.; Rundquist, D.C.; Barrow, T.; Fisher, T.R.; Gurlin, D.; Holz, J. A simple semi-analytical model for remote estimation of chlorophyll-a in turbid waters: Validation. Remote Sens. Environ. 2008, 112, 3582–3593. [Google Scholar] [CrossRef]
  61. Moses, W.J.; Gitelson, A.A.; Berdnikov, S.; Povazhnyy, V. Satellite estimation of chlorophyll-a concentration using the red and NIR bands of MERISThe azov sea case study. IEEE Geosci. Remote Sens. Lett. 2009, 6, 845–849. [Google Scholar] [CrossRef]
  62. Binding, C.E.; Greenberg, T.A.; Bukata, R.P. The MERIS Maximum Chlorophyll Index; its merits and limitations for inland water algal bloom monitoring. J. Great Lakes Res. 2013, 39, 100–107. [Google Scholar] [CrossRef]
  63. Dall’Olmo, G.; Gitelson, A.A. Effect of bio-optical parameter variability and uncertainties in reflectance measurements on the remote estimation of chlorophyll-a concentration in turbid productive waters: Modeling results. Appl. Opt. 2006, 45, 3577–3592. [Google Scholar] [CrossRef] [Green Version]
  64. Dall’Olmo, G.; Gitelson, A.A.; Rundquist, D.C.; Leavitt, B.; Barrow, T.; Holz, J.C. Assessing the potential of SeaWiFS and MODIS for estimating chlorophyll concentration in turbid productive waters using red and near-infrared bands. Remote Sens. Environ. 2005, 96, 176–187. [Google Scholar] [CrossRef]
  65. Gitelson, A.A.; Yacobi, Y.Z. Reflectance in the red and near infra-red ranges of the spectrum as tool for remote chlorophyll estimation in inland waters—Lake kinneret case study. In Proceeding of the Eighteenth Convention of Electrical and Electronics Engineers in Israel, Tel Aviv, Israel, 7–8 March 1995; pp. 1–5. [Google Scholar] [CrossRef]
  66. Zheng, G.; DiGiacomo, P.M. Remote sensing of chlorophyll-a in coastal waters based on the light absorption coefficient of phytoplankton. Remote Sens. Environ. 2017, 201, 331–341. [Google Scholar] [CrossRef]
  67. Defoin-Platel, M.; Chami, M. How ambiguous is the inverse problem of ocean color in coastal waters? J. Geophys. Res. Ocean. 2007, 112, 1–16. [Google Scholar] [CrossRef] [Green Version]
  68. Matthews, M.W. A current review of empirical procedures of remote sensing in Inland and near-coastal transitional waters. Int. J. Remote Sens. 2011, 32, 6855–6899. [Google Scholar] [CrossRef]
  69. Smith, B.; Pahlevan, N.; Schalles, J.; Ruberg, S.; Errera, R.; Ma, R.; Giardino, C.; Bresciani, M.; Barbosa, C.; Moore, T.; et al. A Chlorophyll-a Algorithm for Landsat-8 Based on Mixture Density Networks. Front. Remote Sens. 2021, 1, 1–17. [Google Scholar] [CrossRef]
  70. Mishra, D.R.; Ogashawara, I.; Gitelson, A.A. Bio-Optical Modeling and Remote Sensing of Inland Waters; Elsevier: Amsterdam, The Netherlands, 2017; ISBN 9780128046548. [Google Scholar]
  71. Palmer, S.C.J.; Hunter, P.D.; Lankester, T.; Hubbard, S.; Spyrakos, E.; Tyler, A.N.; Présing, M.; Horváth, H.; Lamb, A.; Balzter, H.; et al. Validation of Envisat MERIS algorithms for chlorophyll retrieval in a large, turbid and optically-complex shallow lake. Remote Sens. Environ. 2015, 157, 158–169. [Google Scholar] [CrossRef] [Green Version]
  72. Qi, L.; Lee, Z.; Hu, C.; Wang, M. Requirement of minimal signal-to-noise ratios of ocean color sensors and uncertainties of ocean color products. J. Geophys. Res. Ocean. 2017, 122, 1–22. [Google Scholar] [CrossRef]
  73. Wang, M.; Gordon, H.R. Sensor performance requirements for atmospheric correction of satellite ocean color remote sensing. Opt. Express 2018, 26, 7390. [Google Scholar] [CrossRef] [PubMed]
  74. Hu, C.; Feng, L.; Lee, Z.; Davis, C.O.; Mannino, A.; McClain, C.R.; Franz, B.A. Dynamic range and sensitivity requirements of satellite ocean color sensors: Learning from the past. Appl. Opt. 2012, 51, 6045–6062. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Jorge, D.S.F.; Barbosa, C.C.F.; de Carvalho, L.A.S.; Affonso, A.G.; Lobo, F.D.L.; Novo, E.M.L.D.M. SNR (signal-to-noise ratio) impact on water constituent retrieval from simulated images of optically complex Amazon lakes. Remote Sens. 2017, 9, 644. [Google Scholar] [CrossRef] [Green Version]
  76. Chondrogianni, C.; Ariztegui, D.; Guilizzoni, P.; Lami, A. ALBANO e NEMI. Mem. Ist. Ital. Idrobiol. 1996, 55, 17–22. [Google Scholar]
  77. Medici, F. Laghi Albano e di Nemi: Carenza idrica e alterazione della qualità delle acque. Geol. Dell’ambiente Period. Trimest. Soc. Ital. Geol. Ambient. 2005, 1, 8–11. [Google Scholar]
  78. Carapezza, M.L.; Lelli, M.; Tarchini, L. Geochemistry of the Albano and Nemi crater lakes in the volcanic district of Alban Hills (Rome, Italy). J. Volcanol. Geotherm. Res. 2008, 178, 297–304. [Google Scholar] [CrossRef] [Green Version]
  79. Cioni, R.; Guidi, M.; Raco, B.; Marini, L.; Gambardella, B. Water chemistry of Lake Albano (Italy). J. Volcanol. Geotherm. Res. 2003, 120, 179–195. [Google Scholar] [CrossRef]
Figure 1. Location of the three investigated volcanic lakes within the study area. Datum: WGS 84, UTM 32N.
Figure 1. Location of the three investigated volcanic lakes within the study area. Datum: WGS 84, UTM 32N.
Remotesensing 13 02699 g001
Figure 2. Water thermal profiles in 2019 of: (a) Lake Bracciano, (b) Lake Albano, and (c) Lake Nemi (data from two sampling dates in the month of May are available: May = 6 May 2019 and May2 = 31 May 2019). Blue lines represent temperature profiles of the mixing phase, during which little or no differences between the surface and deep layers of the water column were observed. Red lines represent thermal stratification, where the water column is divided into epilimnion (top, warmer layer), metalimnion (with rapid temperature change), and hypolimnion (bottom, colder layer).
Figure 2. Water thermal profiles in 2019 of: (a) Lake Bracciano, (b) Lake Albano, and (c) Lake Nemi (data from two sampling dates in the month of May are available: May = 6 May 2019 and May2 = 31 May 2019). Blue lines represent temperature profiles of the mixing phase, during which little or no differences between the surface and deep layers of the water column were observed. Red lines represent thermal stratification, where the water column is divided into epilimnion (top, warmer layer), metalimnion (with rapid temperature change), and hypolimnion (bottom, colder layer).
Remotesensing 13 02699 g002
Figure 3. (a) Coefficient plot of the linear mixed model designed with the index based on the C2RCC (INN). Dots represent coefficient estimates; the outer error bars are 95% confidence intervals and the inner error bars are 50% confidence intervals. The 95% coefficient interval for stratification includes zero, which indicates a non-significant relationship between measured and predicted chlorophyll content. (b) Plot of the interaction effect between Chl concentration and water circulation phase.
Figure 3. (a) Coefficient plot of the linear mixed model designed with the index based on the C2RCC (INN). Dots represent coefficient estimates; the outer error bars are 95% confidence intervals and the inner error bars are 50% confidence intervals. The 95% coefficient interval for stratification includes zero, which indicates a non-significant relationship between measured and predicted chlorophyll content. (b) Plot of the interaction effect between Chl concentration and water circulation phase.
Remotesensing 13 02699 g003
Figure 4. (a) Coefficient plot of the linear mixed model designed with the index based on the MDN algorithm (IMDN). Dots represent coefficient estimates; the outer error bars are 95% confidence intervals and the inner error bars are 50% confidence intervals. The 95% coefficient interval for stratification includes zero, which indicates a non-significant relationship between measured and predicted chlorophyll content. (b) Plot of the interaction effect between Chl concentration and water circulation phase.
Figure 4. (a) Coefficient plot of the linear mixed model designed with the index based on the MDN algorithm (IMDN). Dots represent coefficient estimates; the outer error bars are 95% confidence intervals and the inner error bars are 50% confidence intervals. The 95% coefficient interval for stratification includes zero, which indicates a non-significant relationship between measured and predicted chlorophyll content. (b) Plot of the interaction effect between Chl concentration and water circulation phase.
Remotesensing 13 02699 g004
Table 1. Sentinel-2 image acquisition dates and sampling dates (“-“ indicates the absence of suitable satellite data).
Table 1. Sentinel-2 image acquisition dates and sampling dates (“-“ indicates the absence of suitable satellite data).
Lake AlbanoLake BraccianoLake Nemi
Sampling DateImage Acquisition DateTime Difference (d)Sampling DateImage Acquisition DateTime Difference (d)Sampling DateImage Acquisition DateTime Difference (d)
19 March 201922 March 2019318 March 201922 March 2019419 March 201922 March 20193
2 April 20191 April 201911 April 201930 March 201922 April 20191 April 20191
16 April 2019--17 April 201919 April 2019216 April 2019--
16 May 2019-- 16 May 2019--
28 May 20195 June 2019830 May 201931 May 2019129 May 20195 June 20197
12 June 201915 June 2019310 June 201913 June 2019312 June 201915 June 20193
27 June 201925 June 2019226 June 201925 June 2019127 June 201925 June 20192
9 July 20195 July 20194 9 July 20195 July 20194
23 July 201925 July 2019225 July 201925 July 2019023 July 201925 July 20192
5 September 20193 September 2019217 September 201918 September 201915 September 20193 September 20192
21 October 201923 October 20192 21 October 201923 October 20192
Table 2. Statistical description (mean, minimum, maximum, and standard deviation) of water properties: total chlorophyll concentration (Chl) measured spectrophotometrically, chlorophyll-a concentration (Chl-a) measured spectrophotometrically, epilimnetic water temperature (Temp), salinity (Sal), pH, dissolved oxygen concentration (DO) in the studied lakes.
Table 2. Statistical description (mean, minimum, maximum, and standard deviation) of water properties: total chlorophyll concentration (Chl) measured spectrophotometrically, chlorophyll-a concentration (Chl-a) measured spectrophotometrically, epilimnetic water temperature (Temp), salinity (Sal), pH, dissolved oxygen concentration (DO) in the studied lakes.
LakeChl (mg/L)Chl-a (mg/L)Temp (°C)
MeanMinMaxSDMeanMinMaxSDMeanMinMaxSD
Bracciano1.680.443.730.951.290.442.510.6018.3710.4726.476.03
Albano3.690.4220.24.563.170.3319.14.221.111.130.16.3
Nemi2.150.4210.862.811.870.338.552.520.410.328.76.1
LakeSal (‰)PHDO (mg/L)
MeanMinMaxSDMeanMinMaxSDMeanMinMaxSD
Bracciano0.270.260.270.018.317.148.920.379.578.310.380.79
Albano0.240.220.250.018.428.048.650.150.278.5710.688.42
Nemi0.160.150.160.008.307.368.690.278.967.5910.310.94
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Perrone, M.; Scalici, M.; Conti, L.; Moravec, D.; Kropáček, J.; Sighicelli, M.; Lecce, F.; Malavasi, M. Water Mixing Conditions Influence Sentinel-2 Monitoring of Chlorophyll Content in Monomictic Lakes. Remote Sens. 2021, 13, 2699. https://doi.org/10.3390/rs13142699

AMA Style

Perrone M, Scalici M, Conti L, Moravec D, Kropáček J, Sighicelli M, Lecce F, Malavasi M. Water Mixing Conditions Influence Sentinel-2 Monitoring of Chlorophyll Content in Monomictic Lakes. Remote Sensing. 2021; 13(14):2699. https://doi.org/10.3390/rs13142699

Chicago/Turabian Style

Perrone, Michela, Massimiliano Scalici, Luisa Conti, David Moravec, Jan Kropáček, Maria Sighicelli, Francesca Lecce, and Marco Malavasi. 2021. "Water Mixing Conditions Influence Sentinel-2 Monitoring of Chlorophyll Content in Monomictic Lakes" Remote Sensing 13, no. 14: 2699. https://doi.org/10.3390/rs13142699

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop