On the Adequacy of Representing Water Reflectance by Semi-Analytical Models in Ocean Color Remote Sensing

Deterministic or statistical inversion schemes to retrieve ocean color from space often use a simplified water reflectance model that may introduce unrealistic constraints on the solution, a disadvantage compared with standard, two-step algorithms that make minimal assumptions about the water signal. In view of this, the semi-analytical models of Morel and Maritorena (2001), MM01, and Park and Ruddick (2005), PR05, used in the spectral matching POLYMER algorithm (Steinmetz et al., 2011), are examined in terms of their ability to restitute properly, i.e., with sufficient accuracy, water reflectance. The approach is to infer water reflectance at MODIS wavelengths, as in POLYMER, from theoretical simulations (using Hydrolight with fluorescence and Raman scattering) and, separately, from measurements (AERONET-OC network). A wide range of Case 1 and Case 2 waters, except extremely turbid waters, are included in the simulations and sampled in the measurements. The reflectance model parameters that give the best fit with the simulated data or the measurements are determined. The accuracy of the reconstructed water reflectance and its effect on the retrieval of inherent optical properties (IOPs) is quantified. The impact of cloud and aerosol transmittance, fixed to unity in the POLYMER scheme, on model performance is also evaluated. Agreement is generally good between model results and Hydrolight simulations or AERONET-OC values, even in optically complex waters, with discrepancies much smaller than typical atmospheric correction errors. Significant differences exist in some cases, but having a more intricate model (i.e., using more parameters) makes convergence more difficult. The trade-off is between efficiency/robustness and accuracy. Notable errors are obtained when using the model estimates to retrieve IOPs. Importantly, the model parameters that best fit the input data, in particular chlorophyll-a concentration, do not represent adequately actual values. The reconstructed water reflectance should be used in bio-optical algorithms. While neglecting cloud and aerosol transmittances degrades the accuracy of the reconstructed water reflectance and the retrieved IOPs, it negligibly affects water reflectance ratios and, therefore, any variable derived from such ratios.


Introduction
Ocean color depends on the absorption and scattering properties of the constituents of the water body.Observation of spectral water reflectance from space is a major tool to gather information about water constituents and associated biogeochemical processes such as primary production.A number of atmospheric correction algorithms have been proposed to remove the influence of the atmosphere and surface in the satellite imagery and, therefore, retrieve the water signal (see e.g., [1]).One approach is to determine simultaneously the key properties of aerosols and water constituents of the coupled system by obtaining a best fit to the measured top-of-atmosphere (TOA) reflectance [2][3][4][5][6].The advantage of this approach, compared with the standard two-step algorithm first suggested by [7], resides in its ability to manage situations of both Case 1 and Case 2 waters, i.e., waters whose inherent optical properties (IOPs) are either dominated by phytoplankton or markedly influenced by other constituents (colored dissolved organic matter and inorganic mineral particles), respectively, but this requires that the selected water reflectance model should be able to represent well the bio-optical variability across diverse aquatic ecosystems.Such atmospheric correction algorithms, whether deterministic or statistical, often use a simplified water reflectance model, such as the three-component model of Sathyendranath et al. [8], so that the absorption and scattering coefficients of the ocean water can be specified by a few parameters (e.g., chlorophyll concentration, backscattering and absorption coefficients at a given wavelength).The model, however, may not be representative of worldwide water conditions [9-12], since many variables affecting reflectance are fixed at some average values.It is generally preferable to make no assumption about the variable to retrieve.Such model may affect the retrieval of water reflectance and the accuracy of the derived ocean color products.Steinmetz et al. [13] developed an algorithm called POLYMER to perform atmospheric correction in regions affected by sun glint and retrieve the spectrum of water reflectance.The POLYMER algorithm relies on two models: one is a simple polynomial atmospheric model fitting the scattering of the atmosphere and sun glint signal, and the other is a bio-optical water reflectance model.A spectral matching method is used to obtain the best fit of the TOA reflectance, and the water reflectance model parameters are retrieved in an iterative process.The POLYMER algorithm works effectively in the presence of Sun glint and thin clouds, increasing the useful coverage of satellite measurements.The European Space Agency (ESA) has selected POLYMER for operational generation of OC-CCI (Ocean Colour Climate Change Initiative) products [14,15].One can choose from two water reflectance models in the POLYMER code: A modified version of the semi-analytical model of Morel and Maritorena [16], denoted MM01, and the semi-analytical model of Park and Ruddick [17], denoted PR05.The PR05 model is used as default in the latest version of the POLYMER code and the MM01 model is provided for optional usage.The MM01 model, slightly modified, depends on chlorophyll-a concentration and a backscattering coefficient for non-algal particles.The PR05 model depends on chlorophyll-a concentration, a parameter specifying the contribution of algal and non-algal particles to the backscattering coefficient, and a parameter allowing different absorption coefficients for dissolved organic matter.Accurate representation of the water reflectance is crucial to the success of the POLYMER algorithm.The adequacy of those models, therefore, needs to be assessed.
The objective of this study is to examine MM01 and PR05 in terms of their ability to represent properly water reflectance in a POLYMER-type inversion scheme applied to typical ocean-color sensors measuring in the spectral range 400-900 nm.The models are first evaluated theoretically, using Hydrolight simulations with IOPs specified from an IOCCG synthesized dataset [18], and then experimentally, using AERONET-OC measurements [19].For this, POLYMER is applied to the simulated data or the measurements assuming no interference of the atmosphere and/or surface.The accuracy of the reconstructed water reflectance is quantified, as well as its influence on the derived IOPs using the Quasi-Analytical Algorithm (QAA) [20].The impact of cloud and aerosol transmittance on the retrieval of water reflectance, neglected in POLYMER, is also evaluated.

Reflectance Models
The original MM01 model [16] is slightly modified by taking into account the absorption and backscattering of non-algal particles.The MM01 model uses two parameters: the chlorophyll-a concentration ( [Chla]) and the backscattering coefficient of non-algal particles (bbs), which varies spectrally in λ −1 (λ is the wavelength).This allows applicability to Case 2 water situations.The model input also includes sun and sensor viewing geometry and wind speed.Specifically, the input geometry was set to match that of the specific dataset described later and the default value of wind speed was used, i.e., 5 m/s, as done in POLYMER.The output of this model is the spectral irradiance reflectance just below the surface, from which remote sensing reflectance, R rs , is deduced.In the transformation, the f/Q correction is performed according to Morel and Gentili [21].The similarity spectrum of Ruddick et al. [22] is used to extend the model from 700 to 900 nm.
The PR05 model is built from radiative transfer simulations for a wide range of IOPs that cover both Case 1 and Case 2 waters [17].It is based on a generic parameterization of reflectance as a function of the total backscattering coefficients over the sum of the total absorption coefficients and the total backscattering coefficients.The PR05 model depends on three parameters: [Chla], a phase function parameter defined by the contribution of suspended particles to the backscattering coefficient (fb), and a parameter allowing different absorption coefficients for dissolved organic matter (fa).In POLYMER the default setting of the PR05 model is using [Chla] and fb as free parameters and fixing fa to 1.In this study, the PR05 model with both two and three free parameters was evaluated.Raman scattering is taken into account in the variant used, by applying Raman correction coefficients from Westberry et al. [23].The model input also includes sun and sensor viewing geometry, which were set to match the geometry of the specific dataset described later.The default wind speed is 5 m/s.The output of this model is the spectral bidirectional water reflectance just above the surface (R w ), defined as πL w /E d where L w is water-leaving radiance and E d the downward solar irradiance just above the surface, from which R rs = L w /E d is deduced by dividing by π.

Inversion Scheme
The model parameters, i.e., [Chla] and bbs in MM01, and [Chla], fa, and fb in PR05, are obtained from the prescribed R rs at a set of wavelengths in the 400-900 nm spectral range using, as in POLYMER, the Nelder-Mead optimization scheme [24].A simplex method is used to minimize the cost function f, which is defined as: where Řrs (λ i ) is the modeled remote sensing reflectance at wavelength λ i , and R rs (λ i ) is the Hydrolight simulated or field measured remote sensing reflectance at wavelength λ i .For the MM01 model, the first iteration simplex is defined by initial values log 10 ([Chla]) = −1 and bbs = 0, and initial steps 0.05 and 0.0005 m −1 .For the PR05 model, the first iteration simplex is defined by initial values log 10 ([Chla]) = −1, log 10 (fb) = 0 and initial steps 0.2 and 0.2 when using two parameters, or by initial values log 10 ([Chla]) = -1, log 10 (fb) = 0, log 10 (fa) = 0 and initial steps 0.2, 0.2, and 0.2 when using three parameters.The optimization scheme stops when f reaches a threshold value of 0.005 or the iteration reaches a maximum number.

Data Source
Two datasets were used in this study to evaluate the performance of the MM01 and PR05 water reflectance models, i.e., their ability to properly reconstruct a prescribed water reflectance spectrum.One dataset was obtained from radiative transfer simulations, the other from measurements at a variety of coastal sites.
The first dataset contained 500 cases generated using Hydrolight.These cases cover a wide range of Case 1 and Case 2 waters, but do not include extremely turbid water situations, such as those encountered in estuaries and inland water bodies.The simulations were made from 300 to 900 nm at 10 nm intervals.The input data including the IOPs and concentrations of phytoplankton, colored dissolved organic matter (CDOM), and detritus/minerals were taken from the synthesized dataset from IOCCG Report 5 [18].The IOCCG synthesized dataset only provides IOPs at the wavelength range of 400-800 nm at 10 nm intervals.To extend the IOP values to the wavelengths below 400 nm (necessary to account for Raman scattering) and above 800 nm the following procedures were applied:

•
Phytoplankton absorption.The phytoplankton absorption, a ph , is expressed as a ph = [Chla] a * ph , where a * ph is the chlorophyll-specific absorption coefficient, both provided in the IOCCG dataset.The values of a * ph at 350-400 nm were adopted from Morel [25] and extrapolated to 300 nm, which was then normalized to the IOCCG a * ph value at 400 nm to ensure continuity.The a * ph values in the 800-900 nm range were assumed constant and fixed at the value at 800 nm.

•
Backscattering of phytoplankton.The attenuation of phytoplankton, c ph, was modeled as c ph = c ph (550) (550/λ) n1 in the IOCCG dataset.The value of n1 was determined using c ph at 400-800 nm, which was then used to extend c ph to 300-1000 nm.The backscattering of phytoplankton, bb ph , was thus obtained using bb ph = b ph (c ph − a ph ) with b ph equal to 0.01.

•
Backscattering of detritus, mineral, and other particles.The backscattering of detritus, mineral, and others, bb dm , was modeled as bb dm = bb dm (550) (550/ λ) n2 .The value of n2 was determined using bb dm at 400-800 nm, which was then used to extend bb dm to 300-900 nm.
All simulations were done with a solar zenith angle of 30 • and a wind speed of 5 m/s.Pure seawater properties were specified by Pope and Fry [26] and Smith and Baker [27] for absorption and scattering, respectively.Clear sky condition and infinitely deep homogeneous water were assumed.The direct and diffuse solar irradiance were simulated using a semi-empirical sky model [28] with the annual average sun-earth distance and ozone content of 300 DU as input.Raman scattering and chlorophyll and CDOM fluorescence were also included in all simulations.The simulated nadir-viewed remote sensing reflectance, R rs , was further interpolated into MODIS wavelengths (412, 443, 488, 531, 547, 667, 678, 748, and 869 nm).Note that the IOCCG synthesized IOPs do not cover all possible natural waters.However, as the models and parameters of describing input IOPs are based on extensive field measurements of oligotrophic and eutrophic waters [8, [29][30][31][32][33][34][35], the Hydrolight simulated data should be consistent with a wide range of field observations.
The second dataset originated from the ocean color component of the Aerosol Robotic Network (AERONET-OC), which provides consistent and accurate long-term measurements collected by globally distributed autonomous radiometer systems deployed on offshore fixed platforms [19].In the AERONET-OC dataset, the normalized water-leaving radiance L WN (λ) at various wavelengths in the visible and near-infrared spectral regions is the primary ocean color radiometric product.Data at quality-assurance level 2.0, normalized to the geometry of the Sun at zenith and the view at nadir, was used.The level 2.0 data is fully quality controlled including pre-and post-field calibration with differences smaller than 5%, automatic cloud removal, and manual inspection.Conversion to R rs was accomplished by multiplying L WN (λ) by extraterrestrial solar irradiance.Note that uncertainties in the AERONET-OC in-situ measurements are about 5% in the blue to green bands, and about 8% in the red band [19].The AERONET-OC measurements were collected using different radiometric instruments and thus the wavelengths used were not exactly the same.Since in this study R rs at MODIS wavelengths (see above) are considered, and since the simulated data set corresponds to those wavelengths, measurements at wavelengths located within the MODIS center wavelengths ± 3 nm were selected.At some of the AERONET-OC sites, in particular the Venise site (45 • N,12 • E), measurements at 547 ± 3 nm were missing.In the end measurements from a total of fourteen sites were used, including twelve coastal sites, i.e., COVE_SEAPRISM (36 and Zeebrugge-MOW1 (51 • N,2 • E), and two lake sites, i.e., Lake_Erie (41 • N,83 • W) and Palgrunden (58 • N, 13 • E).Additional criteria based on the similarity spectrum [22] were applied to check and remove outliers, i.e., only the data for which (i) R rs (869 nm) > −0.0001 sr −1 and (ii) the ratio R rs (667 nm)/R rs (869 nm) > 3 if R rs (869 nm) >0.0008 sr −1 were selected.This resulted in a total number of 9824 R rs measurements at wavelengths close to the MODIS wavelengths.For simplification, the AERONET-OC wavelengths will be referred to as 412, 443, 488, 531, 547, 667, and 869 nm, i.e., the MODIS wavelengths.The AERONET-OC measurements were made only under clear sky conditions.They were collected when wind speed is less than 15 m/s.Chlorophyll-a concentrations are provided in the AERONET-OC dataset, but not from field measurements; they are estimated through an iterative procedure making use of reflectance band ratios [19].AERONET-OC sites for which band-ratio algorithms exist naturally benefit from those regional algorithms.For several AERONET-OC sites, the lack of regional algorithms imposes the application of the standard OC2V4 band-ratio.Because of this, the expected accuracy of AERONET-OC [Chla] varies from site to site.

Model Performance with Hydrolight Simulations
Model performance (MM01, PR05) was first evaluated theoretically using the Hydrolight simulations.Prior to evaluation, Case 1 and Case 2 waters were separated.Distinguishing Case 1 and Case 2 waters was based on the criteria proposed by Morel and Bélanger [36] and Robinson et al. [37].According to Morel and Bélanger [36], the particulate scattering coefficient at 560 nm, b p (560), increases with increasing [Chla].The upper limit of b p (560) can be expressed as Such a relationship was used as input to the bio-optical model of Morel and Maritorena [13] to generate the threshold water irradiance reflectance just beneath the water surface, R w (560), as a function of [Chla].The effect of solar zenith angle was taken into consideration.For a given [Chla], if the observed R w (560) exceeded the threshold R w (560), the water was determined as Case 2, otherwise as Case 1.Note that the threshold of Morel and Bélanger [36] only distinguishes Case 2 turbid water with excessive sediments and detritus and does not separate yellow substance-dominated Case 2 waters.Therefore, the determined Case 1 water situations still included yellow substance-dominated Case 2 water situations.An additional threshold of remote sensing reflectance at 670 nm, i.e., R rs (670) >0.0012 sr -1 , now adopted as turbid water flag by the NASA Ocean Biology Processing Group (starting with the fourth processing), was also applied [37].This second threshold may still not favor the identification of CDOM-dominated waters.With these two criteria, there are a total of 55 spectra classified as Case 1 and 445 spectra as Case 2 (Figure 1).Statistics including number of points (N), bias, root-mean-squared error (RMSE), and coefficient of determination (R 2 ) were used to evaluate the model performance.
Figure 2 shows some examples of R rs spectra reconstructed using the MM01 and PR05 models.Only 2 variable parameters, i.e., [Chla] and fb, are used in PR05, with fa fixed to 1.The MM01 and PR05 spectra generally agree with the Hydrolight simulated spectra, especially for Case 1 water (Figure 2a-c).As [Chla] is increased (Figure 2b,c), the difference between reconstructed and simulated R rs around 680-700 nm becomes larger.This is expected since chlorophyll fluorescence effects are not considered in both models.The ability of the model to reproduce properly R rs is degraded for Case 2 water (Figure 2d-f).In the complex Case 2 situation of Figure 2f, one noticeable feature of the MM01 reconstructed R rs is the large discrepancy at 400-500 nm and 600-650 nm, which is attributed to neglecting CDOM absorption and crudely parameterizing the backscattering by non-algal particles.The disagreement around 550-650 nm between PR05 reconstructed and Hydrolight simulated R rs is due to the fact that the slope of CDOM absorption adopted in the PR05 model, modeled as a function of [Chla], does not reflect well the actual (i.e., prescribed) slope, which ultimately affects the red wavelengths since the spectral-matching scheme minimizes the overall difference.
with excessive sediments and detritus and does not separate yellow substance-dominated Case 2 waters.Therefore, the determined Case 1 water situations still included yellow substance-dominated Case 2 water situations.An additional threshold of remote sensing reflectance at 670 nm, i.e., Rrs(670) >0.0012 sr -1 , now adopted as turbid water flag by the NASA Ocean Biology Processing Group (starting with the fourth processing), was also applied [37].This second threshold may still not favor the identification of CDOM-dominated waters.With these two criteria, there are a total of 55 spectra classified as Case 1 and 445 spectra as Case 2 (Figure 1).Statistics including number of points (N), bias, root-mean-squared error (RMSE), and coefficient of determination (R 2 ) were used to evaluate the model performance.Figure 2 shows some examples of Rrs spectra reconstructed using the MM01 and PR05 models.Only 2 variable parameters, i.e., [Chla] and fb, are used in PR05, with fa fixed to 1.The MM01 and PR05 spectra generally agree with the Hydrolight simulated spectra, especially for Case 1 water (Figure 2a-c).As [Chla] is increased (Figure 2b,c), the difference between reconstructed and simulated Rrs around 680-700 nm becomes larger.This is expected since chlorophyll fluorescence effects are not considered in both models.The ability of the model to reproduce properly Rrs is degraded for Case 2 water (Figure 2d-f).In the complex Case 2 situation of Figure 2f, one noticeable feature of the MM01 reconstructed Rrs is the large discrepancy at 400-500 nm and 600-650 nm, which is attributed to neglecting CDOM absorption and crudely parameterizing the backscattering by nonalgal particles.The disagreement around 550-650 nm between PR05 reconstructed and Hydrolight simulated Rrs is due to the fact that the slope of CDOM absorption adopted in the PR05 model, modeled as a function of [Chla], does not reflect well the actual (i.e., prescribed) slope, which ultimately affects the red wavelengths since the spectral-matching scheme minimizes the overall difference.3, blue) being better than that for Case 2 water (Figs. 3, red), with overall bias and Root-Mean-Squared Error (RMSE) of 0.00038 (7.7%) and 0.00082 sr −1 (16.5%) at 443 nm, −0.00016 (−2.7%) and 0.00070 sr −1 (11.9%) at 488 nm, 0.00013 (2.1%) and 0.00075 sr −1 (12.1%) at 547 nm, and −0.00039 (−19.3) and 0.00124 sr −1 (61.4%) at 667 nm, respectively.The percent values in parentheses are biases or RMSEs normalized by the mean of the simulated data.The relatively high bias and RMSE for Case 2 water at 443 nm is due to the lack of CDOM absorption in the MM01 model and is consistent with Figure 2. The MM01 model performance is degraded at 667 nm.At this wavelength, large bias and RMSE are Figure 3 displays scatter plots of MM01 and PR05 reconstructed R rs versus Hydrolight simulated R rs at selected wavelengths, i.e., 443, 488, 547, and 667 nm, for all cases.The MM01 model performs generally well in reproducing R rs at 443, 488, and 547 nm, with the performance for Case 1 water (Figure 3, blue) being better than that for Case 2 water (Figure 3, red), with overall bias and Root-Mean-Squared Error (RMSE) of 0.00038 (7.7%) and 0.00082 sr −1 (16.5%) at 443 nm, −0.00016 (−2.7%) and 0.00070 sr −1 (11.9%) at 488 nm, 0.00013 (2.1%) and 0.00075 sr −1 (12.1%) at 547 nm, and −0.00039 (−19.3) and 0.00124 sr −1 (61.4%) at 667 nm, respectively.The percent values in parentheses are biases or RMSEs normalized by the mean of the simulated data.The relatively high bias and RMSE for Case 2 water at 443 nm is due to the lack of CDOM absorption in the MM01 model and is consistent with Figure 2. The MM01 model performance is degraded at 667 nm.At this wavelength, large bias and RMSE are found for Case 2 water, although the reconstruction of R rs is still good for Case 1 water.The degradation at 667 nm is due primarily to the backscattering of non-algal particles, not being sufficiently well treated by introducing bbs, and secondarily to fluorescence and CDOM absorption not considered.Since the spectral-matching scheme utilizes information at all the input wavelengths, neglecting CDOM absorption, influential in the blue, may also impact results in the red.The PR05 model shows better performance than the MM01 model in reconstructing R rs at 443, 488, 547, and 667 nm, with bias and RMSE generally reduced, i.e., 0.00024 (4.8%) and 0.00038 sr −1 (7.7%) at 443 nm, −0.00001 (−0.1%) and 0.00039 sr −1 (6.8%) at 488 nm, −0.00030 (−4.9%) and 0.00061 sr −1 (9.8%) at 547 nm, and 0.00015 (7.6%) and 0.00027 sr −1 (13.4%) at 667 nm, respectively.The relatively large bias and RMSE at 547 nm, as compared to 443 and 488 nm, is due to the CDOM absorption slope, which is not adequately represented, as well as the spectral-matching scheme that minimizes the overall difference.The reconstructed R rs at 667 nm is in good agreement with the Hydrolight values, which is consistent with the results shown in Figure 2. The improvement is dramatic at this wavelength when using PR05 instead of MM01 (e.g., RMSE reduced from 61.4% to 13.4%).
Remote Sens. 2019, 11 FOR PEER REVIEW 7 model shows better performance than the MM01 model in reconstructing Rrs at 443, 488, 547, and 667 nm, with bias and RMSE generally reduced, i.e., 0.00024 (4.8%) and 0.00038 sr −1 (7.7%) at 443 nm, −0.00001 (−0.1%) and 0.00039 sr −1 (6.8%) at 488 nm, −0.00030 (−4.9%) and 0.00061 sr −1 (9.8%) at 547 nm, and 0.00015 (7.6%) and 0.00027 sr −1 (13.4%) at 667 nm, respectively.The relatively large bias and RMSE at 547 nm, as compared to 443 and 488 nm, is due to the CDOM absorption slope, which is not adequately represented, as well as the spectral-matching scheme that minimizes the overall difference.The reconstructed Rrs at 667 nm is in good agreement with the Hydrolight values, which is consistent with the results shown in Figure 2. The improvement is dramatic at this wavelength when using PR05 instead of MM01 (e.g., RMSE reduced from 61.4% to 13.4%).Performance statistics were computed separately for Case 1 and Case 2 water situations at all MODIS wavelengths (Table 1).The model results generally agree with the Hydrolight values, even in optically complex waters, with PR05 outperforming MM01 in almost all situations.Based on Table 1, the MM01 model generates Rrs more accurately for Case 1 water than for Case 2 water from 412 nm to 547 nm as the bias and RMSE values are lower, which is consistent with Figures 2 and 3.The bias and RMSE for Case 1 water are less than 4.0% and 6.0%, respectively, while the bias can increase to 17.0% and the RMSE to 26.1% for Case 2 water.As wavelength becomes longer, i.e., from 667 nm to 869 nm, the bias and RMSE error increase for both Case 1 and Case 2 water situations.The reconstruction of Rrs at near infrared wavelengths, i.e., 748 and 869 nm, show large bias and RMSE (up to 45.6% and 78.3%, respectively).This is due mostly to the simplified treatment of non-algal particles via the addition of a backscattering coefficient with spectral dependence in  -1 (and no absorption).Extending Rrs from 700 to 900 nm using the average similarity spectrum for turbid waters may also introduce errors.Similarly, the PR05 model generates more accurate Rrs estimation for Case 1 water than for Case 2 water from 412 nm to 547 nm.At these wavelengths, the largest bias observed is 2.1% and the largest RMSE 4.8% for Case 1 water, and the bias and RMSE increase up to 5.2% and 9.7%, respectively, for Case 2 water.The bias and RMSE errors increase (up to 13.7% and 22.2%, respectively) at 667 and 678 nm, indicating significant degradation in Rrs reconstruction at these wavelengths.Much larger relative errors are observed at the near infrared wavelengths, for both Case Performance statistics were computed separately for Case 1 and Case 2 water situations at all MODIS wavelengths (Table 1).The model results generally agree with the Hydrolight values, even in optically complex waters, with PR05 outperforming MM01 in almost all situations.Based on Table 1, the MM01 model generates R rs more accurately for Case 1 water than for Case 2 water from 412 nm to 547 nm as the bias and RMSE values are lower, which is consistent with Figures 2 and 3.The bias and RMSE for Case 1 water are less than 4.0% and 6.0%, respectively, while the bias can increase to 17.0% and the RMSE to 26.1% for Case 2 water.As wavelength becomes longer, i.e., from 667 nm to 869 nm, the bias and RMSE error increase for both Case 1 and Case 2 water situations.The reconstruction of R rs at near infrared wavelengths, i.e., 748 and 869 nm, show large bias and RMSE (up to 45.6% and 78.3%, respectively).This is due mostly to the simplified treatment of non-algal particles via the addition of a backscattering coefficient with spectral dependence in λ −1 (and no absorption).Extending R rs from 700 to 900 nm using the average similarity spectrum for turbid waters may also introduce errors.Similarly, the PR05 model generates more accurate R rs estimation for Case 1 water than for Case 2 water from 412 nm to 547 nm.At these wavelengths, the largest bias observed is 2.1% and the largest RMSE 4.8% for Case 1 water, and the bias and RMSE increase up to 5.2% and 9.7%, respectively, for Case 2 water.The bias and RMSE errors increase (up to 13.7% and 22.2%, respectively) at 667 and 678 nm, indicating significant degradation in R rs reconstruction at these wavelengths.Much larger relative errors are observed at the near infrared wavelengths, for both Case 1 and Case 2 situations, i.e., an RMSE of 47.3 and 21.7%, respectively, at 748 nm and 102.0 and 72.6%, respectively, at 869 nm, which reflects the generally lower R rs , especially at 869 nm.This could also be partly due to the fact that the input absorption and scattering coefficients at these wavelengths are extrapolated values based on the values at 700 nm instead of based on field measurements.Compared with the POLYMER performance against in situ data obtained for the Medium Resolution Imaging Spectroradiometer (MERIS) in Case 1 water (Müller et al. [14]), the error in reconstructing R rs from MM01 or PR05 (Table 1) is much smaller than the algorithm error, with RMSE values of 0.00008 (MM01) and 0.00012 (PR05) instead of 0.0011 sr −1 at 443 nm, and 0.00008 (MM01) and 0.00008 (PR05) at 547 nm instead of 0.00032 sr −1 at 560 nm.When data collected in both Case 1 and 2 waters are considered [14,38], the POLYMER RMSE values for MERIS are increased to 0.0037 sr −1 and 0.0033 sr −1 at 443 and 560 nm, respectively, which is much larger than the R rs reconstruction errors, i.e., 0.00082 and 0.00038 at 443 nm and 547, respectively (Table 1).For MODIS-Aqua, these authors reported similar POLYMER performance, i.e., RMSE of 0.0036 at 443 nm and 0.0031 sr −1 at 547 nm.Thus, the error introduced by the semi-analytical modeling of water reflectance is unlikely to dominate the POLYMER performance, indicating that using MM01 and, especially, PR05 in the spectral-matching algorithm is adequate.Figure 4 displays scatter plots of the model reconstructed and Hydrolight simulated R rs ratios.Three ratios commonly used in ocean color algorithms [39,40] were selected, i.e., R rs (412)/R rs (547), R rs (443)/R rs (547), and R rs (488)/R rs (547).It is found that the MM01 and PR05 modeled band ratios R rs (412)/R rs (547), R rs (443)/R rs (547), and R rs (488)/R rs (547) are all in very good agreement with the Hydrolight simulations.The overall statistics indicate that the two models perform similarly in restituting R rs ratios: the bias is less than 6% and the RMSE less than 9% for the MM01 model, and the bias is less than 5% and the RMSE less than 10% for the PR05 model.Also, the performance of Case 1 water situations is similar to that of Case 2 water situations (Table 2).This indicates that both the MM01 and PR05 models are capable of provide accurate estimation of R rs ratios, which is important for ocean-color products using band ratio algorithms.If the MM01 and PR05 reconstructed R rs is used in the OCI algorithm [41], the resulted bias and RMSE for estimated [Chla] are −0.000(−0.2%) and 0.006 mgm −3 (4.3%), and −0.001 (−0.7%) and 0.008 mgm −3 (6.6%), respectively, i.e., much less than the errors caused by the OCI algorithm itself.This is due to the fact that the errors on R rs at the various wavelengths considered in the CI equation, i.e., CI ≈ R rs (555) − 0.5(R rs (443) + R rs (667)) where R rs (555) is specified as 0.93R rs (547) [41], are small (on the order of 10 −4 sr −1 ) in the applicability domain of the algorithm (clear, Case 1 waters).Note for both models the Nelder-Mead retrieved [Chla] differs substantially from the prescribed [Chla] (Figure 4).The bias and RMSE of log 10 ([Chla]) are 0.52 and 0.60, and 0.21 and 0.37, respectively, for MM01 and PR05 model, and the errors are similar for both Case 1 and Case 2 waters (Table 2).These values, especially biases, are much larger when compared with the MODIS-Aqua estimated [Chla] performance against in situ data which shows absolute bias and mean absolute error (MAE, always less or equal to RMSE) of 0.07 and 0.22 on log 10 scale, respectively [42].This is probably due to the variations in the relation between the absorption/ backscattering coefficients and the model parameters, which are not accounted for in the models.It is thus concluded that the [Chla] parameter in the models does not represent well the actual [Chla] and should not be used as [Chla] estimate.Instead [Chla] should be derived from the reconstructed R rs , and using R rs in band-ratio or OCI algorithms is appropriate.
restituting Rrs ratios: the bias is less than 6% and the RMSE less than 9% for the MM01 model, and the bias is less than 5% and the RMSE less than 10% for the PR05 model.Also, the performance of Case 1 water situations is similar to that of Case 2 water situations (Table 2).This indicates that both the MM01 and PR05 models are capable of provide accurate estimation of Rrs ratios, which is important for ocean-color products using band ratio algorithms.If the MM01 and PR05 reconstructed Rrs is used in the OCI algorithm [41], the resulted bias and RMSE for estimated [Chla] are −0.000(−0.2%) and 0.006 mgm -3 (4.3%), and −0.001 (−0.7%) and 0.008 mgm -3 (6.6%), respectively, i.e., much less than the errors caused by the OCI algorithm itself.This is due to the fact that the errors on Rrs at the various wavelengths considered in the CI equation, i.e., CI ≈ Rrs(555) − 0.5(Rrs(443) + Rrs(667)) where Rrs( 555) is specified as 0.93Rrs(547) [41], are small (on the order of 10 -4 sr -1 ) in the applicability domain of the algorithm (clear, Case 1 waters).Note for both models the Nelder-Mead retrieved [Chla] differs substantially from the prescribed [Chla] (Fig. 4).The bias and RMSE of log10( [Chla]) are 0.52 and 0.60, and 0.21 and 0.37, respectively, for MM01 and PR05 model, and the errors are similar for both Case 1 and Case 2 waters (Table 2).These values, especially biases, are much larger when compared with the MODIS-Aqua estimated [Chla] performance against in situ data which shows absolute bias and mean absolute error (MAE, always less or equal to RMSE) of 0.07 and 0.22 on log10 scale, respectively [42].This is probably due to the variations in the relation between the absorption/ backscattering coefficients and the model parameters, which are not accounted for in the models.It is thus concluded that the [Chla] parameter in the models does not represent well the actual [Chla] and should not be used as [Chla] estimate.Instead [Chla] should be derived from the reconstructed Rrs, and using Rrs in band-ratio or OCI algorithms is appropriate.547), as well as Nelder-Mead retrieved and prescribed [Chla].Blue circles represent Case 1 water and red represent Case 2 water.The bias and RMSE of [Chla] were calculated on log10 scale.
As seen in Figure 5, the errors generally increase when [Chla] increases for both MM01 and PR05 models.The bias and RMSE obtained for the PR05 model are generally lower than those for the MM01 model through the entire [Chla] range, confirming that the PR05 model provides better reconstruction of Rrs than the MM01 model.Specifically, the large bias and RMSE of the MM01 model, noticeably at 412 nm, 667 nm and 678 nm, are mainly related to the simplified treatment of backscattering of non-algal particles, as well as lack of fluorescence and absorption of CDOM and non-algal particles in the model.The large bias and RMSE at 547 nm of the PR05 model is due to the non-representative CDOM absorption slope, as discussed previously.The results also suggest Table 2. Performance statistics of MM01 and PR05 reconstructed R rs ratios and [Chla] for Case 1 and Case 2 water using Hydrolight simulations.The bias and RMSE of log10( [Chla]) are in the unit of mgm -3 .As seen in Figure 5, the errors generally increase when [Chla] increases for both MM01 and PR05 models.The bias and RMSE obtained for the PR05 model are generally lower than those for the MM01 model through the entire [Chla] range, confirming that the PR05 model provides better reconstruction of R rs than the MM01 model.Specifically, the large bias and RMSE of the MM01 model, noticeably at 412 nm, 667 nm and 678 nm, are mainly related to the simplified treatment of backscattering of non-algal particles, as well as lack of fluorescence and absorption of CDOM and non-algal particles in the model.The large bias and RMSE at 547 nm of the PR05 model is due to the non-representative CDOM absorption slope, as discussed previously.The results also suggest degraded performance for Case 2 waters since these waters have generally higher [Chla] than Case 1 waters (bias and RMSE are larger for productive waters).

Parameter
Case  The MM01 and PR05 reconstructed Rrs were also used as input to the QAA algorithm [20] to retrieve IOPs including absorption of detritus and CDOM at 443 nm (adg(443)), phytoplankton absorption at 443 nm (aph(443)), and backscattering of particles at 547 nm (bbp(547)).When using the Hydrolight simulated Rrs as input to the QAA, the bias and RMSE for retrieved adg(443), aph(443), and bbp(547) are −3.3% and 29.3%, −3.4% and 34.4%, and −1.8% and 55.6%, respectively (Figure 6a-c).These values quantify algorithm errors.When using the reconstructed Rrs as input, the retrieved adg(443), aph(443), and bbp(443) are largely biased, i.e., bias being −71.0%, 57.8%, and −23.7% for MM01, and −31.3%, −11.4%, and −14.6% for PR05, when compared with the prescribed values (Figure 6d-i).The The MM01 and PR05 reconstructed R rs were also used as input to the QAA algorithm [20] to retrieve IOPs including absorption of detritus and CDOM at 443 nm (a dg (443)), phytoplankton absorption at 443 nm (a ph (443)), and backscattering of particles at 547 nm (b bp (547)).When using the Hydrolight simulated R rs as input to the QAA, the bias and RMSE for retrieved a dg (443), a ph (443), and b bp (547) are −3.3% and 29.3%, −3.4% and 34.4%, and −1.8% and 55.6%, respectively (Figure 6a-c).These values quantify algorithm errors.When using the reconstructed R rs as input, the retrieved a dg (443), a ph (443), and b bp (443) are largely biased, i.e., bias being −71.0%, 57.8%, and −23.7% for MM01, and −31.3%, −11.4%, and −14.6% for PR05, when compared with the prescribed values (Figure 6d-i).The errors of using the MM01 and PR05 reconstructed R rs are much higher than those of using Hydrolight simulated R rs .Specifically, the performance of using the PR05 reconstructed R rs is better than that of using the MM01 reconstructed R rs with relatively lower bias and RMSE.The QAA algorithm involves two important parameters: u and χ, the first one depending on R rs values and the second one depending on R rs ratios.The value of b bp (547) depends on u at 547 nm, u(547), and the total absorption coefficient at 547 nm, a(547).Since u(547) is related to R rs (547) only and a(547) is related to the ratio R rs (443)/R rs (547), R rs (488)/R rs (547), and R rs (667)/R rs (547), it is expected that the retrieved b bp (547) values show larger bias and RMSE compared with the retrieval using the Hydrolight simulated R rs , and the errors are larger for the MM01 model than for the PR05 model based on the model performance in reconstructing R rs (667) (Figure 3).The value of a dg (443) is expressed as a function of b bp (547), R rs (412), R rs (443), and R rs (443)/R rs (547).Therefore, the error in the b bp (547) propagates further and results in much larger error in the retrieved a dg (443), as shown in Figure 6d,g.Also, the retrieval error in a dg (443) is less when using the PR05 model reconstructed R rs due to the smaller error in b bp (547).The value of a ph (443) is calculated as the total absorption minus a dg and the absorption of pure sea water, in which way the error in retrieving a ph (443) is compensated to some extent as the total absorption is derived using b bp (547) and R rs (443).It is concluded that using successfully the QAA algorithm requires accurate input of R rs and R rs ratios.Other IOP algorithms such as GSM [43] and GIOP [44] also use R rs as input and may have similar problems.Therefore, the MM01 and PR05 reconstructed R rs should be used with caution to retrieve IOPs.One may envision, however, applying these inversion schemes to reflectance ratios, as long as sensitivity to IOPs remains adequate.
errors of using the MM01 and PR05 reconstructed Rrs are much higher than those of using simulated Rrs.Specifically, the performance of using the PR05 reconstructed Rrs is better than that of using the MM01 reconstructed Rrs with relatively lower bias and RMSE.The QAA algorithm involves two important parameters: u and χ, the first one depending on Rrs values and the second one depending on Rrs ratios.The value of bbp(547) depends on u at 547 nm, u(547), and the total absorption coefficient at 547 nm, a(547).Since u(547) is related to Rrs(547) only and a(547) is related to the ratio Rrs(443)/Rrs(547), Rrs(488)/Rrs(547), and Rrs(667)/Rrs(547), it is expected that the retrieved bbp(547) values show larger bias and RMSE compared with the retrieval using the Hydrolight simulated Rrs, and the errors are larger for the MM01 model than for the PR05 model based on the model performance in reconstructing Rrs(667) (Figure 3).The value of adg( 443) is expressed as a function of bbp(547), Rrs(412), Rrs(443), and Rrs(443)/Rrs(547).Therefore, the error in the bbp(547) propagates further and results in much larger error in the retrieved adg(443), as shown in Figure 6d,g.Also, the retrieval error in adg( 443) is less when using the PR05 model reconstructed Rrs due to the smaller error in bbp(547).The value of aph( 443) is calculated as the total absorption minus adg and the absorption of pure sea water, in which way the error in retrieving aph( 443) is compensated to some extent as the total absorption is derived using bbp(547) and Rrs(443).It is concluded that using successfully the QAA algorithm requires accurate input of Rrs and Rrs ratios.Other IOP algorithms such as GSM [43] and GIOP [44] also use Rrs as input and may have similar problems.Therefore, the MM01 and PR05 reconstructed Rrs should be used with caution to retrieve IOPs.One may envision, however, applying these inversion schemes to reflectance ratios, as long as sensitivity to IOPs remains adequate.The three-parameter PR05 model was evaluated to figure out whether the R rs reconstruction would be improved by introducing more parameters (Table 3).Results show that by introducing the third parameter fa (see Section 2.1), the PR05 model does improve the estimation of R rs at all wavelengths except 531 nm, wavelength for which the bias and RMSE are slightly higher than using a two-parameter model.However, the improvement is not significant for 412, 443, 488, 547, and 667 nm, with bias and RMSE reduced by about 1% or even less.Much more improvement is seen at longer wavelengths (i.e., 678, 748, and 869 nm), but the differences between Hydrolight simulated and PR05 modeled R rs are still large, for example the RMSE at 869 nm is 0.00010 sr −1 (56.1%).Similar results of R rs are observed when separating Case 1 and Case 2 waters.Therefore, having a more intricate model may not be warranted.The trade-off is between efficiency/robustness and accuracy.It also should be noted that having more input parameters in the spectral-matching scheme might not guarantee convergence.It is concluded that the PR05 model with two parameters is preferred for R rs estimation than the three-parameter model (allows faster and easier convergence while keeping similar accuracy).

Sensivity of R rs Reconstruction to Atmospheric Transmittance
In the POLYMER algorithm, the TOA reflectance is decomposed as follows [13]: where t oz (λ) is the ozone transmittance, ρ mol (λ) is the Rayleigh scattering reflectance, ρ gli (λ) is the (non-spectral) sun glint reflectance affected by the direct transmission factor T(λ), ρ aer (λ) is the aerosol reflectance, ρ gam (λ) accounts for the coupling between Fresnel reflectance and scattering by molecules and aerosols, t(λ) is the total (i.e., direct + diffuse) transmittance for atmospheric scattering, and ρ w (λ,0+) is the water reflectance above the water-air interface (ρ w = πR rs ).In this decomposition, the atmosphere is assumed to be clear (i.e., no clouds).In the presence of clouds, Equation ( 2) should be modified to account for scattering by cloud droplets and its interaction with molecule/aerosol scattering and surface reflection, and t(λ) should include cloud transmittance, i.e., t(λ) = t mol (λ) t aer (λ) t cloud (λ), where subscripts mol, aer, and cloud refer to molecules, aerosols, and clouds, The transmittance t(λ) is pre-calculated with a successive orders of scattering (SOS) code [45] and stored in look-up tables, but only molecules are considered in the calculations, even though the algorithm is designed to work in the presence of aerosols and, at least, optically thin clouds.Neglecting the transmittance of aerosols and clouds is expected to affect the R rs reconstruction in the spectral-matching scheme.Therefore, it is important to evaluate the sensitivity of the R rs reconstruction using the MM01 and PR05 models to cloud and aerosol transmittance.
To perform such evaluation, the aerosol total transmittance t aer was modeled following Tanre et al. [46]: where θ s and θ v are the Sun and view zenith angles, respectively, and τ aer is the aerosol optical thickness, τ aer (λ) = 0.1(λ/865) −1 .The cloud transmittance t cloud was parameterized according to Fitzpatrick et al. [47]: With the calculated t aer and t cloud , the Hydrolight simulated R rs were transformed to R rs t aer t cloud , and the MM01 and PR05 model parameters were retrieved using R rs t aer t cloud instead of R rs as input.
Figures 7 and 8 display the resulting scatter plots of MM01 and PR05 reconstructed R rs , [Chla], and R rs ratios versus the Hydrolight simulations.The reconstructed R rs is systematically biased, i.e., by 20-40%, after taking into account the cloud and aerosol transmittance.However, the R rs ratios are still in good agreement with the Hydrolight simulated ratios, i.e., they exhibit bias and RMSE comparable to those obtained without considering the effects of transmittance.When using the MM01 and PR05 reconstructed R rs to estimate [Chla] with the OCI algorithm, the resulting bias and RMSE increase to 0.039 (32.5%) and 0.040 mgm -3 (33.1%), and 0.036 (29.9%) and 0.036 mgm -3 (30.5%), respectively.This is because the cloud and aerosol transmittances introduce a large bias, i.e., more than 20%, to R rs , which is further propagated to the CI index.The cloud and aerosol transmittances do not have much effect on the model-retrieved [Chla]; statistics are similar to those in Figure 4.However, again for both models, the Nelder-Mead [Chla] does not represent well the actual values.Note that most of the impact on the R rs reconstruction performance is due cloud transmittance (0.78 at 443 nm), as the aerosol transmittance (0.95 at 443 nm) is close to 1.When using the reconstructed Rrs as input to the QAA algorithm to retrieve adg(443), aph(443), and bbp(547), large errors are obtained compared with the results using the Hydrolight simulated reflectance spectra (Figure 9).As discussed in Section 3.1, the retrieved value of bbp( 547) is dependent on Rrs(547), Rrs(443)/Rrs(547), Rrs(488)/Rrs(547), and Rrs(667)/Rrs(547).Since the reconstructed Rrs( 547) is systematically biased when the effects of aerosol/cloud transmittance are taken into account (Figure 7), and even though the values of Rrs(443)/Rrs(547), Rrs(488)/Rrs(547), and Rrs(667)/Rrs(547) (Figure 8) are similar to those with no transmittance, the error in the retrieved bbp(547) is increased (compare Figure 9 with Figure 6).However, since adg(443) and aph(443) are not only determined by bbp(547) but also by Rrs(412), Rrs(443), and Rrs(443)/Rrs(547), in the end the errors somewhat compensate, yielding bias and RMSE for adg(443) and aph(443) that are similar to those before accounting for aerosol/cloud When using the reconstructed R rs as input to the QAA algorithm to retrieve a dg (443), a ph (443), and b bp (547), large errors are obtained compared with the results using the Hydrolight simulated reflectance spectra (Figure 9).As discussed in Section 3.1, the retrieved value of b bp (547) is dependent on R rs (547), R rs (443)/R rs (547), R rs (488)/R rs (547), and R rs (667)/R rs (547).Since the reconstructed R rs (547) is systematically biased when the effects of aerosol/cloud transmittance are taken into account (Figure 7), and even though the values of R rs (443)/R rs (547), R rs (488)/R rs (547), and R rs (667)/R rs (547) (Figure 8) are similar to those with no transmittance, the error in the retrieved b bp (547) is increased (compare Figure 9 with Figure 6).However, since a dg (443) and a ph (443) are not only determined by b bp (547) but also by R rs (412), R rs (443), and R rs (443)/R rs (547), in the end the errors somewhat compensate, yielding bias and RMSE for a dg (443) and a ph (443) that are similar to those before accounting for aerosol/cloud transmittance.Again, the errors in the retrieved IOPs are less when using the PR05 model reconstructed R rs , compared with the MM01 results.
Remote Sens. 2019, 11 FOR PEER REVIEW 15 transmittance.Again, the errors in the retrieved IOPs are less when using the PR05 model reconstructed Rrs, compared with the MM01 results.

Model Performance Using AERONET-OC Dataset
Model (MM01, PR05) performance was further evaluated using the AERONET-OC dataset.The AERONET-OC dataset does not include Rrs at all wavelengths, i.e., only Rrs at 412, 443, 488, 531, 547, 667, and 869 nm are available.To determine if missing Rrs at 678 and 748 nm affects the model

Model Performance Using AERONET-OC Dataset
Model (MM01, PR05) performance was further evaluated using the AERONET-OC dataset.The AERONET-OC dataset does not include R rs at all wavelengths, i.e., only R rs at 412, 443, 488, 531, 547, 667, and 869 nm are available.To determine if missing R rs at 678 and 748 nm affects the model performance, tests were made using the 500 Hydrolight simulated cases (Section 2.3) with full MODIS wavelengths and with partial MODIS wavelengths (excluding 678 and 748 nm).No significant difference (i.e., less than 1%) in bias and RMSE was observed between R rs , [Chla], and R rs ratios retrieved using the full or partial set of wavelengths.Therefore, the errors caused by missing wavelengths are regarded as negligible in the following study.Prior to evaluation, Case 1 and Case 2 water situations were separated by using the same criteria that were applied to the Hydrolight simulations, taking the solar zenith angle into account.This resulted in 6094 Case 1 water and 3730 Case 2 water situations (Figure 10).

Model Performance Using AERONET-OC Dataset
Model (MM01, PR05) performance was further evaluated using the AERONET-OC dataset.The AERONET-OC dataset does not include Rrs at all wavelengths, i.e., only Rrs at 412, 443, 488, 531, 547, 667, and 869 nm are available.To determine if missing Rrs at 678 and 748 nm affects the model performance, tests were made using the 500 Hydrolight simulated cases (section 2.3) with full MODIS wavelengths and with partial MODIS wavelengths (excluding 678 and 748 nm).No significant difference (i.e., less than 1%) in bias and RMSE was observed between Rrs, [Chla], and Rrs ratios retrieved using the full or partial set of wavelengths.Therefore, the errors caused by missing wavelengths are regarded as negligible in the following study.Prior to evaluation, Case 1 and Case 2 water situations were separated by using the same criteria that were applied to the Hydrolight simulations, taking the solar zenith angle into account.This resulted in 6094 Case 1 water and 3730 Case 2 water situations (Figure 10).As shown in Figure 11, top and Table 4, the MM01 model performs generally well in reconstructing R rs at 443, 488, and 547 nm.The bias and RMSE values are −0.00001(−0.2%) and 0.00032 sr −1 (9.2%), −0.00055 (−11.7%) and 0.00085 sr −1 (18.0%), 0.00025 (5.5%) and 0.00056 sr −1 (12.2%) respectively, i.e., higher than those obtained using Hydrolight simulations.Much larger bias and RMSE are found at 667 nm, which are −0.00015(−12.1%) and 0.0010 sr −1 (84.5%), respectively.Nevertheless, the errors in R rs are still smaller than typical POLYMER atmospheric correction errors.The PR05 model performs better in reconstructing R rs with generally lower bias and RMSE at each wavelength (Figure 11, bottom, and Table 4).Unlike for the MM01 model, the PR05 reconstructed R rs values at 667 nm are in good agreement with the AERONET-OC measurements, even though fluorescence is not considered in the PR05 model.Similar to the MM01 results, the bias and RMSE values of the PR05 reconstruction are higher than those obtained using Hydrolight simulations.As discussed in the previous section, the simplified representation of the backscattering of non-algal particles, together with the lack of consideration of fluorescence and CDOM and non-algal absorption in the MM01 model are the main contributors to the errors in reconstructed R rs .The PR05 model takes into account the absorption and backscattering of non-algal particles by using the [Chla] and fb parameters, but the slope used for CDOM absorption might not be sufficiently well representative, causing substantial errors at 547 nm for Case 2 water situations when applying the spectral-matching scheme for inversion.Note that the Hydrolight simulated R rs is not subjected to radiometric errors and other, e.g., environmental, perturbations and this may partly explain why the bias and RMSE in reconstructed R rs using the AERONET-OC dataset are larger than those obtained using Hydrolight simulations.
into account the absorption and backscattering of non-algal particles by using the [Chla] and fb parameters, but the slope used for CDOM absorption might not be sufficiently well representative, causing substantial errors at 547 nm for Case 2 water situations when applying the spectral-matching scheme for inversion.Note that the Hydrolight simulated Rrs is not subjected to radiometric errors and other, e.g., environmental, perturbations and this may partly explain why the bias and RMSE in reconstructed Rrs using the AERONET-OC dataset are larger than those obtained using Hydrolight simulations.547) are all in good agreement with the AERONET-OC measurements.The overall bias is less than 15% and the overall RMSE less than 18% for the MM01 model.The PR05 reconstruction gives slightly better results, with the overall bias less than 11% and the overall RMSE less than 16%.As expected, the errors in the reconstructed Rrs ratios are larger than those obtained using Hydrolight simulations.Table 5 shows that the model reconstructed Rrs ratios agree well with the prescribed values for both Case 1 and Case 2 waters, with errors slightly lower for Case 1 situations than for Case 2 situations.When applying the MM01 and PR05 reconstructed Rrs to estimate [Chla] using the OCI algorithm, the resulting bias and RMSE for [Chla] are 0.010 (4.8%) and 0.014 mgm -3 (6.8%), and 0.007 (3.6%) and 0.011 mgm -3 (5.4%), respectively.This confirms that both the MM01 and PR05 reconstructed Rrs can be used to estimate [Chla] using band ratio and OCI algorithms.Figure 12 also   Figure 12 displays scatter plots of the model reconstructed and AERONET-OC measured R rs ratios.Results show that the MM01 and PR05 band ratios R rs (412)/R rs (547), R rs (443)/R rs (547), and R rs (488)/ R rs (547) are all in good agreement with the AERONET-OC measurements.The overall bias is less than 15% and the overall RMSE less than 18% for the MM01 model.The PR05 reconstruction gives slightly better results, with the overall bias less than 11% and the overall RMSE less than 16%.As expected, the errors in the reconstructed R rs ratios are larger than those obtained using Hydrolight simulations.Table 5 shows that the model reconstructed R rs ratios agree well with the prescribed values for both Case 1 and Case 2 waters, with errors slightly lower for Case 1 situations than for Case 2 situations.When applying the MM01 and PR05 reconstructed R rs to estimate [Chla] using the OCI algorithm, the resulting bias and RMSE for [Chla] are 0.010 (4.8%) and 0.014 mgm -3 (6.8%), and 0.007 (3.6%) and 0.011 mgm -3 (5.4%), respectively.This confirms that both the MM01 and PR05 reconstructed R rs can be used to estimate [Chla] using band ratio and OCI algorithms.Figure 12 also shows that the model-retrieved [Chla] differs substantially from the AERONET-OC [Chla], both models tending to overestimate [Chla].However, for both models the overall agreement is better than when using Hydrolight simulations: the bias and RMSE on log10 scale are 0.33 and 0.42 instead of 0.52 and 0.60, respectively, for the MM01 model, and 0.12 and 0.23 instead of 0.21 and 0.37, respectively, for the PR05 model.The statistics are similar when Case 1 and Case 2 waters are used separately (Table 5).The comparisons suggest that the bio-optical relationships used in both models, based on field measurements, are more representative of actual situations than those in the Hydrolight simulations.Despite the improvement in [Chla] estimate, the model-retrieved [Chla] should still not be considered as representative of the actual [Chla].The bias and RMSE of [Chla] were calculated on log10 scale.When using the MM01 model for Rrs reconstruction, relatively large bias and RMSE are observed at 488 nm (when [Chla] is within 2-4 mg m −3 ) and 667 nm (when [Chla] is larger than 2 mg m −3 ) (Figure 13).Such behavior could be due to the fact that the absorption and backscattering of nonalgal particles, fluorescence, and CDOM absorption are not taken into account in the MM01 model.For other wavelengths the bias and RMSE are relatively small and do not quite change with [Chla], except when [Chla] is within 2-4 mg m −3 (small increase).This is also the case when using PR05 model.This increased bias and RMS when [Chla] is in the range 2-4 mg m −3 could be due to the fact that the AERONET-OC [Chla] are not measured but estimated from remote sensing reflectance [19], and that higher uncertainty exists in this [Chla] range.It is also possible the variability of chlorophylla absorption is not sufficiently well represented in the MM01 and PR05 model.The bias and RMSE are smaller over the entire [Chla] range when using the PR05 model, indicating the overall better performance of that model in reconstructing water reflectance, again confirming the results from simulations.When using the MM01 model for R rs reconstruction, relatively large bias and RMSE are observed at 488 nm (when [Chla] is within 2-4 mg m −3 ) and 667 nm (when [Chla] is larger than 2 mg m −3 ) (Figure 13).Such behavior could be due to the fact that the absorption and backscattering of non-algal particles, fluorescence, and CDOM absorption are not taken into account in the MM01 model.For other wavelengths the bias and RMSE are relatively small and do not quite change with [Chla], except when [Chla] is within 2-4 mg m −3 (small increase).This is also the case when using PR05 model.This increased bias and RMS when [Chla] is in the range 2-4 mg m −3 could be due to the fact that the AERONET-OC [Chla] are not measured but estimated from remote sensing reflectance [19], and that higher uncertainty exists in this [Chla] range.It is also possible the variability of chlorophyll-a absorption is not sufficiently well represented in the MM01 and PR05 model.The bias and RMSE are smaller over the entire Figure 13.Bias and RMSE as a function of [Chla], based on MM01 and PR05 reconstructed results using AERONET-OC dataset.Eight [Chla] bins were used and the number of points was identical for each bin.

Summary and Conclusions
The MM01 and PR05 models used in the POLYMER spectral-matching algorithm have been examined in terms of their ability to properly represent water reflectance over the wide range of situations/regimes expected in oceans.The MM01 model, in the modified version considered, depends on chlorophyll-a concentration and backscattering coefficient for non-algal particles.The PR05 model depends on chlorophyll-a concentration, a parameter specifying the contribution of algal and non-algal particles to the backscattering coefficient, and a parameter allowing different absorption coefficients for dissolved organic matter.The two models have been evaluated separately for Case 1 and Case 2 waters at MODIS wavelengths (412, 443, 488, 531, 547, 667, 678, 748, and 869 nm).
A theoretical evaluation was first performed using Hydrolight simulations with IOPs from an IOCCG synthesized dataset (500 cases) representing a wide range of field observations.Results show that the MM01 model generally performs well in retrieving Rrs at 412, 443, 488, 531, and 547 nm, even in optically complex waters.Neglecting CDOM absorption in the MM01 model is responsible for the relatively large bias and RMS found at 412 nm.The discrepancies found at 667 and 869 nm for Case 2 water are mainly due to the backscattering of inorganic particles not well represented in the MM01 model.Other factors, including compensation by the spectral-matching scheme of differences at shorter wavelengths caused by the lack of CDOM and non-algal absorption and Rrs extrapolation above 700 nm using the similarity spectrum, may also contribute to degradation of performance at 667 and 869 nm.Despite the errors in Rrs, the MM01 retrieved Rrs ratios are in very good agreement with the Hydrolight Rrs ratios.The PR05 model, on the other hand, shows better performance than the MM01 model in reconstructing Rrs.The PR05 retrieved Rrs values, obtained by fixing to unity the parameter allowing variable absorption by dissolved organic matter, agree with the Hydrolight values for both Case 1 and Case 2 situations at all visible wavelengths including 667 and 678 nm, even though fluorescence is not considered in the model.The relatively large bias and RMSE at 547  [Chla], based on MM01 and PR05 reconstructed results using AERONET-OC dataset.Eight [Chla] bins were used and the number of points was identical for each bin.

Summary and Conclusions
The MM01 and PR05 models used in the POLYMER spectral-matching algorithm have been examined in terms of their ability to properly represent water reflectance over the wide range of situations/regimes expected in oceans.The MM01 model, in the modified version considered, depends on chlorophyll-a concentration and backscattering coefficient for non-algal particles.The PR05 model depends on chlorophyll-a concentration, a parameter specifying the contribution of algal and non-algal particles to the backscattering coefficient, and a parameter allowing different absorption coefficients for dissolved organic matter.The two models have been evaluated separately for Case 1 and Case 2 waters at MODIS wavelengths (412, 443, 488, 531, 547, 667, 678, 748, and 869 nm).
A theoretical evaluation was first performed using Hydrolight simulations with IOPs from an IOCCG synthesized dataset (500 cases) representing a wide range of field observations.Results show that the MM01 model generally performs well in retrieving R rs at 412, 443, 488, 531, and 547 nm, even in optically complex waters.Neglecting CDOM absorption in the MM01 model is responsible for the relatively large bias and RMS found at 412 nm.The discrepancies found at 667 and 869 nm for Case 2 water are mainly due to the backscattering of inorganic particles not well represented in the MM01 model.Other factors, including compensation by the spectral-matching scheme of differences at shorter wavelengths caused by the lack of CDOM and non-algal absorption and R rs extrapolation above 700 nm using the similarity spectrum, may also contribute to degradation of performance at 667 and 869 nm.Despite the errors in R rs , the MM01 retrieved R rs ratios are in very good agreement with the Hydrolight R rs ratios.The PR05 model, on the other hand, shows better performance than the MM01 model in reconstructing R rs .The PR05 retrieved R rs values, obtained by fixing to unity the parameter allowing variable absorption by dissolved organic matter, agree with the Hydrolight values for both Case 1 and Case 2 situations at all visible wavelengths including 667 and 678 nm, even though fluorescence is not considered in the model.The relatively large bias and RMSE at 547 nm is due to the non-representative slope of CDOM absorption and the spectral-matching scheme that minimizes the overall difference.Large discrepancy with the Hydrolight R rs values is found at 748 and 869 nm, which reflects the generally lower R rs , especially at 869 nm, and could be partly due to the fact that the input absorption and scattering coefficients at these wavelengths are extrapolated from values at 700 nm instead of based on field measurements.Overall, the performance for Case 1 waters is better than that of for Case 2 waters.The PR05 reconstructed R rs ratios also agree well with the Hydrolight ratios.The reconstructed R rs from both models, when used in the OCI algorithm for oligotrophic waters, yields an accurate estimation of chlorophyll-a concentration.However, the MM01 and PR05 retrieved R rs , when used as input to the QAA algorithm, gives significant errors in the estimation of a ph (443), a dg (443), and b bp (547).This is due to the fact that the QAA algorithm uses R rs at all visible wavelengths and not only R rs ratios.Therefore, caution should be exercised when using the model-retrieved R rs in algorithms for IOP estimations such as QAA.Importantly, the model parameters that best fit the input data, in particular chlorophyll-a concentration, may not represent adequately actual values.The impact of cloud and aerosol transmittance (assumed to be 1 in POLYMER) on R rs reconstruction was also evaluated theoretically using Hydrolight simulations.While the cloud and aerosol transmittances do not affect significantly the reconstruction of R rs ratios, they do affect the accuracy of reconstructed R rs , the OCI estimated chlorophyll-a concentration, and the derived IOPs when using QAA and, presumably, other IOP algorithms.
The MM01 and PR05 models were further evaluated using the AERONET-OC dataset, which consists of 9824 samples from a total of 15 sites.Only R rs measurements at 412, 443, 488, 531, 547, 667, and 869 nm were available, but the missing data at 678 and 748 nm does not affect significantly the R rs retrievals as tested using the Hydrolight simulations.Similar conclusions to those obtained with simulated R rs are drawn for the model-retrieved R rs and R rs ratios, confirming the theoretical results, but with relatively higher bias and RMSE probably due to the uncertainties existing in field measurements, and the IOP variability not captured by both models.In any case, the uncertainty in reconstructing R rs is much smaller than typical POLYMER atmospheric correction errors obtained against in-situ measurements.Again, the model-retrieved chlorophyll-a concentration does not agree well with the AERONET-OC values, which may partly result from the application of non-site-specific bio-optical algorithms.Therefore, it is recommended that the reconstructed water reflectance, not the retrieved model parameters, should be used in bio-optical algorithms.
The three-parameter PR05 model was also evaluated using the Hydrolight simulations, to examine whether using more parameters would improve the reconstruction of the water reflectance.Reconstruction was improved with three parameters, but slightly.Using more parameters, however, increases computation time and is less likely to guarantee convergence of the spectral matching scheme.The trade-off is between efficiency/robustness and accuracy.It is concluded that the two-parameter PR05 model is preferred and appropriate to represent the reflectance of a wide range of Case 1 and case 2 waters in spectral-matching atmospheric correction schemes such as POLYMER, and probably other (e.g., statistical) one-step inversion schemes, at least until progress in atmospheric correction necessitates otherwise.

Figure 1 .
Figure 1.Spectra of remote sensing reflectance Rrs for Case 1(a) and Case 2(b) water using Hydrolight simulations.To distinguish Case 1 and 2 water situations, two thresholds were used: the threshold of Morel and Bélanger [32], i.e., the limiting Rw(560) values produced by the bio-optical model of Morel and Maritorena [13] with the upper limit of bp(560) = 0.69[Chla] 0.766 as input, and the threshold of Robinson et al. [33], i.e.Rrs(670) ≤0.0012 sr -1 for Case 1 water, otherwise Case 2 water.

Figure 1 .
Figure 1.Spectra of remote sensing reflectance R rs for Case 1 (a) and Case 2 (b) water using Hydrolight simulations.To distinguish Case 1 and 2 water situations, two thresholds were used: the threshold of Morel and Bélanger [32], i.e., the limiting R w (560) values produced by the bio-optical model of Morel and Maritorena [13] with the upper limit of b p (560) = 0.69[Chla] 0.766 as input, and the threshold of Robinson et al. [33], i.e.R rs (670) ≤ 0.0012 sr −1 for Case 1 water, otherwise Case 2 water.

Figure 3
Figure3displays scatter plots of MM01 and PR05 reconstructed Rrs versus Hydrolight simulated Rrs at selected wavelengths, i.e., 443, 488, 547, and 667 nm, for all cases.The MM01 model performs generally well in reproducing Rrs at 443, 488, and 547 nm, with the performance for Case 1 water (Figs. 3, blue) being better than that for Case 2 water (Figs.3, red), with overall bias and Root-Mean-Squared Error (RMSE) of 0.00038 (7.7%) and 0.00082 sr −1 (16.5%) at 443 nm, −0.00016 (−2.7%) and 0.00070 sr −1 (11.9%) at 488 nm, 0.00013 (2.1%) and 0.00075 sr −1 (12.1%) at 547 nm, and −0.00039 (−19.3) and 0.00124 sr −1 (61.4%) at 667 nm, respectively.The percent values in parentheses are biases or RMSEs normalized by the mean of the simulated data.The relatively high bias and RMSE for Case 2 water at 443 nm is due to the lack of CDOM absorption in the MM01 model and is consistent with Figure2.The MM01 model performance is degraded at 667 nm.At this wavelength, large bias and RMSE are

Figure 2 .
Figure 2. Examples of MM01 and PR05 reconstructed R rs spectra versus Hydrolight simulated R rs spectra, for Case 1 (a-c) and Case 2 (d-f) water situations.The spectral matching procedure was performed at MODIS wavelengths (412, 443, 488, 531, 547, 667, 678, 748, and 869 nm), and the retrieved model parameters were then used to produce modeled R rs spectra from 400 nm to 900 nm at 10 nm interval.Black lines represent Hydrolight simulations, red lines represent MM01 reconstructed R rs spectra, and blue lines PR05 reconstructed R rs spectra.The dots show model reconstructed R rs at MODIS wavelengths.

Figure 3 .
Figure 3.Comparison of MM01 (a-d) and PR05 (e-h) reconstructed remote sensing reflectance R rs with Hydrolight simulated R rs at selected MODIS wavelengths: 443 nm, 488 nm, 547 nm, and 667 nm.Blue circles represent Case 1 water and red represent Case 2 water.

Figure 5 .
Figure 5. Bias and RMSE as a function of[Chla], based on MM01 and PR05 reconstructed results using Hydrolight simulations.Eight[Chla]  bins were used and the number of points was identical for each bin.

Figure 5 .
Figure 5. Bias and RMSE as a function of[Chla], based on MM01 and PR05 reconstructed results using Hydrolight simulations.Eight[Chla]  bins were used and the number of points was identical for each bin.

Figure 6 .
Figure 6.Comparison of QAA derived a dg (443), a ph (443), and b bp (547) using Hydrolight simulated R rs (a-c), MM01 reconstructed R rs (d-f), and PR05 reconstructed R rs (g-i) versus prescribed values.Blue circles represent Case 1 water and red represent Case 2 water.
) where τ cloud is geometric-optics cloud optical thickness, α is surface albedo.Cloud optical depth τ cloud is set to 3 (thin cloud), and surface albedo α to 0.06.Sun and view zenith angles are those of the IOCCG dataset, i.e., 30 • and 0 • , respectively.The coefficients a 1 , b 1 , b 2 , b 3 , k 1 , k 2 , k 3 , c, and d are tabulated in Fitzpatrick et al. [47] as a function of wavelength.The resulting t cloud values were interpolated to the MODIS wavelengths.

Figure 7 .
Figure 7.Comparison of the MM01 (a-d) and PR05 (e-h) reconstructed remote sensing reflectance R rs with Hydrolight simulated R rs at selected MODIS wavelengths (443, 488, 547, and 667 nm), after applying cloud and aerosol transmittance.Blue circles represent Case 1 water and red represent Case 2 water.

Figure 9 .
Figure 9.Comparison of QAA derived a dg (443), a ph (443), and b bp (547) using MM01 reconstructed R rs (a-c) and PR05 reconstructed R rs (d-f) versus prescribed values, after applying cloud and aerosol transmittance.Blue circles represent Case 1 water and red represent Case 2 water.

Figure 10 .
Figure 10.AERONET-OC spectra of remote sensing reflectance Rrs of Case 1(a) and Case 2(b) waters.To distinguish Case 1 and 2 water situations, two thresholds were used: the threshold of Morel and Bélanger [32], i.e., the limiting Rw(560) values produced by the bio-optical model of Morel and Maritorena [13] with the upper limit of bp(560) = 0.69[Chla] 0.766 as input, and the threshold of Robinson et al. [33], i.e.Rrs(670) ≤ 0.0012 sr -1 for Case 1 water, otherwise Case 2 water.

Figure 10 .
Figure 10.AERONET-OC spectra of remote sensing reflectance R rs of Case 1 (a) and Case 2 (b) waters.To distinguish Case 1 and 2 water situations, two thresholds were used: the threshold of Morel and Bélanger [32], i.e., the limiting R w (560) values produced by the bio-optical model of Morel and Maritorena [13] with the upper limit of b p (560) = 0.69[Chla] 0.766 as input, and the threshold of Robinson et al. [33], i.e.R rs (670) ≤ 0.0012 sr −1 for Case 1 water, otherwise Case 2 water.

Figure 12
Figure 12 displays scatter plots of the model reconstructed and AERONET-OC measured Rrs ratios.Results show that the MM01 and PR05 band ratios Rrs(412)/Rrs(547), Rrs(443)/Rrs(547), and Rrs(488)/ Rrs(547) are all in good agreement with the AERONET-OC measurements.The overall bias is less than 15% and the overall RMSE less than 18% for the MM01 model.The PR05 reconstruction gives slightly better results, with the overall bias less than 11% and the overall RMSE less than 16%.As expected, the errors in the reconstructed Rrs ratios are larger than those obtained using Hydrolight simulations.Table5shows that the model reconstructed Rrs ratios agree well with the prescribed values for both Case 1 and Case 2 waters, with errors slightly lower for Case 1 situations than for Case 2 situations.When applying the MM01 and PR05 reconstructed Rrs to estimate[Chla]  using the OCI algorithm, the resulting bias and RMSE for[Chla]  are 0.010 (4.8%) and 0.014 mgm -3 (6.8%), and 0.007 (3.6%) and 0.011 mgm -3 (5.4%), respectively.This confirms that both the MM01 and PR05 reconstructed Rrs can be used to estimate[Chla]  using band ratio and OCI algorithms.Figure12 also

Figure 11 .
Figure 11.Comparison of MM01 (a-d) and PR05 (e-h) reconstructed remote sensing reflectance R rs with AERNOET-OC R rs at selected MODIS wavelengths: 443 nm, 488 nm, 547 nm, and 667 nm.Blue circles represent Case 1 water and red represent Case 2 water.

Table 4 .
Performance statistics of MM01 and PR05 reconstructed R rs for Case 1 and Case 2 water using AERONET-OC dataset.

Figure 13 .
Figure13.Bias and RMSE as a function of[Chla], based on MM01 and PR05 reconstructed results using AERONET-OC dataset.Eight[Chla]  bins were used and the number of points was identical for each bin.

Table 1 .
Performance statistics of MM01 and PR05 reconstructed R rs for Case 1 and Case 2 water using Hydrolight simulations.

Table 3 .
Performance statistics of three-parameter PR05 reconstructed R rs for Case 1 and Case 2 water using Hydrolight simulations.

Table 5 .
Performance statistics of MM01 and PR05 reconstructed Rrs ratios and log10([Chla]) for Case 1 and Case 2 water using AERONET-OC dataset.The bias and RMSE of log10([Chla]) are in the unit of mg m −3 .