Predicting minimum uncertainties in the inversion of ocean color geophysical parameters based on Cramer-Rao bounds

We present an analytical approach based on Cramer-Rao Bounds (CRBs) to investigate the uncertainties in estimated ocean color parameters resulting from the propagation of uncertainties in the bio-optical reflectance modeling through the inversion process. Based on given bio-optical and noise probabilistic models, CRBs can be computed efficiently for any set of ocean color parameters and any sensor configuration, directly providing the minimum estimation variance that can be possibly attained by any unbiased estimator of any targeted parameter. Here, CRBs are explicitly developed using (1) two water reflectance models corresponding to deep and shallow waters, resp., and (2) four probabilistic models describing the environmental noises observed within four Sentinel-2 MSI, HICO, Sentinel-3 OLCI and MODIS images, resp. For both deep and shallow waters, CRBs are shown to be consistent with the experimental estimation variances obtained using two published remote-sensing methods, while not requiring one to perform any inversion. CRBs are also used to investigate to what extent perfect a priori knowledge on one or several geophysical parameters can improve the estimation of remaining unknown parameters. For example, using pre-existing knowledge of bathymetry (e.g., derived from LiDAR) within the inversion is shown to greatly improve the retrieval of bottom cover for shallow waters. Finally, CRBs are shown to provide valuable information on the best estimation performances that may be achieved with the MSI, HICO, OLCI and MODIS configurations for a variety of oceanic, coastal and inland waters. CRBs are thus demonstrated to be an informative and efficient tool to characterize minimum uncertainties in inverted ocean color geophysical parameters. © 2017 Optical Society of America OCIS codes: (010.4450) Oceanic optics; (100.3190) Inverse problems; (110.4234) Multispectral and hyperspectral imaging; (280.0280) Remote sensing and sensors. References and links 1. P. J. Werdell, B. A. Franz, S. W. Bailey, G. C. Feldman, E. Boss, V. E. Brando, M. Dowell, T. Hirata, S. J. Lavender, Z. Lee, H. Loisel, S. Maritorena, F. Mélin, T. S. Moore, T. J. Smyth, D. Antoine, E. Devred, O. H. F. D’Andon, and A. Mangin, “Generalized ocean color inversion model for retrieving marine inherent optical properties,” Appl. Opt. 52, 2019–2037 (2013). 2. J. D. Hedley, C. M. Roelfsema, I. Chollett, A. R. Harborne, S. F. Heron, S. J. Weeks, W. J. Skirving, A. E. Strong, C. M. Eakin, T. R. Christensen, V. Ticzon, S. Bejarano, and P. J. Mumby, “Remote sensing of coral reefs for monitoring and management: A review,” Remote Sens. 8, 118 (2016). 3. R. A. Garcia, P. R. Fearns, and L. I. McKinna, “Detecting trend and seasonal changes in bathymetry derived from HICO imagery: A case study of Shark Bay, Western Australia,” Remote Sens. Environ. 147, 186–205 (2014). Vol. 26, No. 2 | 22 Jan 2018 | OPTICS EXPRESS A1 #302324 https://doi.org/10.1364/OE.26.0000A1 Journal © 2018 Received 13 Jul 2017; revised 5 Oct 2017; accepted 8 Oct 2017; published 27 Nov 2017 4. T. S. Moore, J. W. Campbell, and M. D. Dowell, “A class-based approach to characterizing and mapping the uncertainty of the MODIS ocean chlorophyll product,” Remote Sens. Environ. 113, 2424–2430 (2009). 5. M. S. Salama and A. Stein, “Error decomposition and estimation of inherent optical properties,” Appl. Opt. 48, 4947–4962 (2009). 6. J. D. Hedley, C. M. Roelfsema, S. R. Phinn, and P. J. Mumby, “Environmental and sensor limitations in optical remote sensing of coral reefs: Implications for monitoring and sensor design,” Remote Sens. 4, 271–302 (2012). 7. W. W. Gregg, “Assimilation of SeaWiFS ocean chlorophyll data into a three-dimensional global ocean model,” J. Mar. Syst. 69, 205–225 (2008). 8. Z. Lee, R. Arnone, C. Hu, P. J. Werdell, and B. Lubac, “Uncertainties of optical parameters and their propagations in an analytical ocean color inversion algorithm,” Appl. Opt. 49, 369–381 (2010). 9. N. Pahlevan, S. Sarkar, and B. A. Franz, “Uncertainties in coastal ocean color products: Impacts of spatial sampling,” Remote Sens. Environ. 181, 14–26 (2016). 10. Z. Lee and K. L. Carder, “Effect of spectral band numbers on the retrieval of water column and bottom properties from ocean color data,” Appl. Opt. 41, 2191–2201 (2002). 11. P. Wang, E. S. Boss, and C. Roesler, “Uncertainties of inherent optical properties obtained from semianalytical inversions of ocean color,” Appl. Opt. 44, 4074–4085 (2005). 12. Z. Lee, A. Weidemann, and R. Arnone, “Combined Effect of reduced band number and increased bandwidth on shallow water remote sensing: The case of Worldview 2,” IEEE Trans. Geosci. Remote Sens. 51, 2577–2586 (2013). 13. D. R. Lyzenga, “Passive remote sensing techniques for mapping water depth and bottom features,” Appl. Opt. 17, 379–383 (1978). 14. Z. Lee, K. L. Carder, C. D. Mobley, R. G. Steward, and J. S. Patch, “Hyperspectral remote sensing for shallow waters: 2. deriving bottom depths and water properties by optimization,” Appl. Opt. 38, 3831–3843 (1999). 15. V. Brando, J. Anstee, M. Wettle, A. Dekker, S. Phinn, and C. Roelfsema, “A physics based retrieval and quality assessment of bathymetry from suboptimal hyperspectral data,” Remote Sens. Environ. 113, 755–770 (2009). 16. J. Hedley, C. Roelfsema, and S. R. Phinn, “Efficient radiative transfer model inversion for remote sensing applications,” Remote Sens. Environ. 113, 2527–2532 (2009). 17. A. G. Dekker, S. R. Phinn, J. Anstee, P. Bissett, V. E. Brando, B. Casey, P. Fearns, J. Hedley, W. Klonowski, Z. P. Lee, M. Lynch, M. Lyons, C. Mobley, and C. Roelfsema, “Intercomparison of shallow water bathymetry, hydro-optics, and benthos mapping techniques in Australian and Caribbean coastal environments,” Limnol. Oceanogr. Methods 9, 396–425 (2011). 18. S. Jay and M. Guillaume, “A novel maximum likelihood based method for mapping depth and water quality from hyperspectral remote-sensing data,” Remote Sens. Environ. 147, 121–132 (2014). 19. L. I. McKinna, P. R. Fearns, S. J. Weeks, P. J. Werdell, M. Reichstetter, B. A. Franz, D. M. Shea, and G. C. Feldman, “A semianalytical ocean color inversion algorithm with explicit water column depth and substrate reflectance parameterization,” J. Geophys. Res. Oceans 120, 1741–1770 (2015). 20. S. Jay and M. Guillaume, “Regularized estimation of bathymetry and water quality using hyperspectral remote sensing,” Int. J. Remote Sens. 37, 263–289 (2016). 21. S. Jay, M. Guillaume, A. Minghelli, Y. Deville, M. Chami, B. Lafrance, and V. Serfaty, “Hyperspectral remote sensing of shallow waters: Considering environmental noise and bottom intra-class variability for modeling and inversion of water reflectance,” Remote Sens. Environ. 200, 352–367 (2017). 22. J. E. O’Reilly, S. Maritorena, B. G. Mitchell, D. A. Siegel, K. L. Carder, S. A. Garver, M. Kahru, and C. McClain, “Ocean color chlorophyll algorithms for SeaWiFS,” J. Geophys. Res. Oceans 103, 24937–24953 (1998). 23. Z. Lee, K. L. Carder, and R. A. Arnone, “Deriving inherent optical properties from water color: a multiband quasi-analytical algorithm for optically deep waters,” Appl. Opt. 41, 5755–5772 (2002). 24. D. B. Gillis, J. H. Bowles, and W. J. Moses, “Improving the retrieval of water inherent optical properties in noisy hyperspectral data through statistical modeling,” Opt. Express 21, 21306–21316 (2013). 25. P. H. Garthwaite, I. T. Jolliffe, and B. Jones, Statistical inference (Oxford University, 2002). 26. A. Roueff, A. Arnaubec, P. C. Dubois-Fernandez, and P. Refregier, “Cramer–rao lower bound analysis of vegetation height estimation with random volume over ground model and polarimetric sar interferometry,” IEEE Geosci. Remote Sens. Lett 8, 1115–1119 (2011). 27. X. Liu, S. Bourennane, and C. Fossati, “Denoising of hyperspectral images using the PARAFAC model and statistical performance analysis,” IEEE Trans. Geosci. Remote Sens. 50, 3717–3724 (2012). 28. S. Jay and M. Guillaume, “Estimation of water column parameters with a maximum likelihood approach,” in “3rd Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing,” (IEEE, 2011), pp. 1–4. 29. J. Hedley, C. Roelfsema, B. Koetz, and S. Phinn, “Capability of the Sentinel-2 mission for tropical coral reef mapping and coral bleaching detection,” Remote Sens. Environ. 120, 145–155 (2012). 30. R. A. Garcia, L. I. McKinna, J. D. Hedley, and P. R. Fearns, “Improving the optimization solution for a semi-analytical shallow water inversion model in the presence of spectrally correlated noise,” Limnol. Oceanogr. Methods 12, 651–669 (2014). 31. V. E. Brando and A. G. Dekker, “Satellite hyperspectral remote sensing for estimating estuarine and coastal water quality,” IEEE Trans. Geosci. Remote Sens. 41, 1378–1387 (2003). 32. Z. Lee, K. L. Carder, C. D. Mobley, R. G. Steward, and J. S. Patch, “Hyperspectral remote sensing for shallow waters. i. a semianalytical model,” Appl. Opt. 37, 6329–6338 (1998). Vol. 26, No. 2 | 22 Jan 2018 | OPTICS EXPRESS A2


Introduction
Optical remote sensing is of major importance to monitor optically deep and shallow waters from air-and satellite-borne sensors.In the case of optically deep waters, it enables the estimation of the inherent optical properties (IOPs) of the upper water layer, i.e., the spectral absorption and backscattering coefficients.In addition to pure water, the IOPs can be related to three optically active water constituents, namely, phytoplankton, colored dissolved organic matter (CDOM) and non-algal suspended particles.Spatially-and temporally-resolved knowledge of these variables is critical to understand biogeochemical processes such as carbon exchanges, phytoplankton biodiversity shifts, and responses to climatic disturbances [1].In the case of optically shallow waters, optical remote sensing also enables the estimation of depth and bottom cover (both potentially greatly affecting the water-leaving radiance), which has important implications, e.g., for monitoring coral reefs [2] or detecting seasonal changes in bathymetry [3].Providing reliable uncertainties (e.g., expressed in terms of mean absolute difference [4] or root mean square deviation [5]) to remote-sensing products is an important task since such knowledge informs one of the confidence level of these products, which is critical for many applications.For example, knowledge of chlorophyll-a uncertainty is needed to characterize the uncertainty budget of primary production algorithms [4].Uncertainties are also required to guide the design of future sensors [6], or to improve the predictions of ocean models through data assimilation [7].Accurate knowledge of uncertainties is especially critical for remote sensing of aquatic environments, since water absorption often leads to low reflectances compared to other Earth surfaces such as soil or vegetation, and therefore to high uncertainties.
First of all, it is worth reminding that the uncertainty in a given ocean color product has several components, each of which arises from a particular stage of the whole calibration, measurement and processing chain.When assessing uncertainties by comparing the actual product values measured in-situ with their satellite-derived estimates, one actually takes into account the uncertainties related to, e.g., in-situ measurements, correction schemes (e.g., atmospheric correction), bio-optical modeling, sensor noise and radiometric specifications, as well as those induced by the inversion method [4,5,8,9].Note that the uncertainty in the inversion of a particular ocean color parameter usually depends on the actual parameter value itself: for example, in the case of shallow water remote sensing, the uncertainties in the inversion of bathymetry and bottom cover increase with depth since the bottom influence becomes negligible as depth increases.Instead of providing a single uncertainty value for the whole range of possible parameter values, a more informative approach consists in subdividing the parameter space and computing the uncertainty for each subdivision [4,8,10].For example, in [4], the authors characterized the chlorophyll-a uncertainty of the Aqua MODerate resolution Imaging Spectroradiometer (MODIS) for eight optical water types, each of which is defined by a mean reflectance spectrum and a spectral covariance matrix.The uncertainty can then be computed for each image pixel based on the memberships to these water types.Still for optically deep waters, other approaches have also been proposed to provide an uncertainty value for each estimate.In [11], the uncertainty is derived from the results of more than 1,000 inversions performed for various spectral shapes of phytoplankton, CDOM and suspended particle IOPs.In [8], it is obtained by propagating uncertainties through the Quasi-Analytical Algorithm.In [1], it is based on the shape of the cost function (that is used for model inversion) around the optimum through the use of the Jacobian matrix.
Simulations have also widely been used to evaluate for some of the uncertainty components.The two main advantages are that (1) one can generate as many simulations as needed, and (2) there is no error in ground truth measurements, correction schemes and modeling, such that only the uncertainties related to the sensor radiometric specifications, the noise model and the inverse method, are evaluated.Such a separation is critical to quantify the influence of each source of uncertainty [5,6,8,9].The use of simulations is especially interesting for shallow water remote sensing because of the higher number of model parameters compared to deep waters [6].Simulated data have been used to compare sensors [12] or to identify the main factors confounding the discrimination of coral reefs [6].However, the number of parameters is such that, for computational convenience, the analysis is often simplified by sampling the parameter space and by using a few bottom covers, sensor configurations or signal-to-noise ratios (SNRs).
Most of the above mentioned approaches thus require remote-sensing data (either real or simulated) to evaluate uncertainties, and the implementation of these approaches may be computationally expensive.In addition, all of them depend on the inversion method, although a wide variety of statistically-and physically-based methods have been developed for remote sensing of optically shallow [13][14][15][16][17][18][19][20][21] and deep waters [1,[22][23][24].Therefore, one cannot know to what extent the uncertainties obtained using a given algorithm could be improved using another one.As emphasized in [6], theoretical studies remain required to characterize the remote-sensing uncertainties, but also to help the design of appropriate sensors.
In this context, the development of Cramer-Rao Bounds (CRBs) is relevant since CRBs provide theoretical lower bounds for the variances of unbiased estimators [25].For a given water column and noise covariance matrix, the CRB analytical expressions thus directly provide the minimum achievable variances in unbiased estimators for all targeted ocean color geophysical parameters.Since CRBs only depend on the Fisher information related to the modeled reflectance data [25], they are algorithm-independent and do not require one to invert remote-sensing data.Applications of CRBs in remote sensing include polarimetric synthetic aperture radar interferometry [26] or denoising of hyperspectral images [27], but, to our knowledge, they have not been fully developed in the context of ocean color remote sensing.They have only been introduced in [28] for shallow waters in a simple yet not realistic case (i.e., one unknown parameter, known bottom cover, diagonal spectral covariance matrix).However, a more comprehensive study is still required to evaluate the interest of CRBs for aquatic remote sensing, as knowledge of such minimum inversion uncertainties may be very useful, e.g., to design satellite-borne sensors or to provide lower bounds for remote-sensing data assimilation within oceanic models.
In this paper, we use the reflectance model presented in Section 2 to develop the CRB expressions for both optically deep and shallow waters (Section 3 and Appendix A).Three applications are then presented in Section 4. First, the minimum uncertainties in the inversion of ocean color parameters, as provided by the CRBs, are compared to the experimental estimation standard deviations obtained using the LS (Least-Square) [1,14] and MILE (MaxImum Likelihood estimation including Environmental noise) [21] methods (Section 4.2).Then, the CRBS are used to investigate how perfect a priori knowledge on some model parameters can improve the estimation of the remaining unknown parameters (Section 4.3).Finally, the CRBs are used to predict the minimum uncertainties in the inversion of ocean color parameters for four sensor configurations, namely, the HICO (Hyperspectral Imager for the Coastal Ocean), MODIS, Sentinel-2 MSI (MultiSpectral Instrument) and Sentinel-3 OLCI (Ocean and Land Colour Instrument) configurations (Section 4.4).Note that, unless otherwise stated, the term "inversion uncertainty" does not hereafter correspond to the total uncertainty in a given ocean color product.Instead, it refers to the uncertainty in the inversion of a given ocean color parameter, and results from the propagation of uncertainties in the bio-optical modeling of reflectance (such uncertainties being accounted for by the noise model presented in Section 2.1) through the inversion process.The term "inversion uncertainty" is here seen as the estimator standard deviation, provided this estimator is unbiased.

Noise probabilistic modeling
Similarly as in [21,29,30], the spectral variability of subsurface remote-sensing reflectance r r s = [r r s (λ 1 ), ..., r r s (λ L )] t (where L is the number of wavebands) is described using a multivariate Gaussian distribution, i.e., r r s = µ + n s where µ and n s are L × 1 vectors.The mean subsurface remote-sensing reflectance spectrum µ is parameterized using either a shallow (Section 2.2) or deep (Section 2.3) water reflectance model.n s is a Gaussian vector with zero mean and covariance matrix Γ s that accounts for the environmental noise [31] (note that all the per-band noise variances are given by the diagonal of Γ s ).The environmental noise includes all the sources of above-water reflectance variability that are not accounted for by the bio-optical model µ (e.g., sensor noise, variable residuals from atmospheric correction or effects related to the rough water surface) and may be therefore scene-specific.In practice, it is generally estimated over a homogeneous area of optically deep water.

Bio-optical modeling for optically shallow waters
For optically shallow waters, the mean subsurface remote-sensing reflectance µ(λ) at a given wavelength λ is given by the semi-analytical model developed in [14,15,32]: where r r s,∞ (λ) is the deep water reflectance, H is the depth, and the bottom albedo is given by a linear mixture of two substratum albedos ρ b,1 (λ) and ρ b,2 (λ) with fractional covers B and (2) (5) where θ s and θ v are respectively the solar and viewing zenith angles.Finally, the absorption and backscattering coefficients are given by the sum of the contributions of optically active water constituents, i.e., pure water, phytoplankton, CDOM and suspended particles [17]: where a w (λ) and b bw (λ) are the absorption and backscattering coefficients of pure water [33,34], a 0 (λ) and a 1 (λ) are empirical spectra [32], a φ (440) is the phytoplankton absorption coefficient at 440 nm, a g (440) is the CDOM absorption coefficient at 440 nm, a * g (λ) = exp [−0.015(λ − 440)], b bp (550) is the particle backscattering coefficient at 550 nm and b * bp (λ) = (550/λ) 0.5 .Note that, in the following, the known a w (λ), a 0 (λ), a 1 (λ), a * g (λ), b bw (λ), b * bp (λ), ρ b,1 (λ) and ρ b,2 (λ) spectra were oversampled to 1 nm step before being spectrally-degraded according to each tested sensor radiometric configuration (Section 4.1), assuming a rectangular spectral response for each waveband.For a given sun-sensor geometry, the shallow water reflectance is thus fully determined by the five geophysical parameters H, a φ (440), a g (440), b bp (550) and B.

Development of Cramer-Rao lower bounds
Let us define ∆ as the vector of unknown ocean color geophysical parameters (i.e., H, a φ (440), a g (440), b bp (550) and/or B) to be estimated from the remote-sensing observation r r s (either for optically deep or shallow waters).According to [25], the variance of any unbiased estimator ∆ i (r r s ) of any element ∆ i of ∆ is bounded as where CRB(∆) is a square matrix with dimensions equal to the length of ∆.The CRB of ∆ i estimation thus corresponds to the i th element of the diagonal of CRB(∆).Further, the inverse of CRB(∆) is the Fisher information matrix I F (∆) [25], whose element i, j is where P(r r s |∆) is the likelihood of observing r r s given ∆.In the multivariate Gaussian case with mean µ and covariance matrix Γ s (Section 2.1), Eq. ( 10) can be simplified [35]: where tr is the trace operator.As Γ s does not depend on unknown parameters, Eq. ( 11) becomes When using the models described in Section 2, the determination of I F (∆) thus only requires the calculation of derivatives ∂µ/∂∆ i , as well as the estimation of the environmental noise Γ s .The components of ∆ are H, a φ (440), a g (440), b bp (550) and/or B in the case of shallow waters (Section 2.2), while they are a φ (440), a g (440) and/or b bp (550) in the case of deep waters (Section 2.3).The analytical expressions of the five derivatives are developed in Appendix A. I F (∆) is then inverted numerically using a standard matrix inversion algorithm to obtain the CRB matrix CRB(∆).Based on such a modeling, CRB(∆) only includes the uncertainties related to the considered Fisher information (which depends on the bio-optical model, environmental noise and sensor radiometric specifications), but not those related to, e.g., ground truth measurements, systematic errors in correction schemes (e.g., atmospheric correction), or inverse methods.
In other words, CRB(∆) are the minimum uncertainties in ∆(r r s ) due to the propagation of uncertainties in the bio-optical model µ through the inversion process.Uncertainties in µ refer to all the factors causing pixel-to-pixel reflectance variations that are not accounted for by µ, such factors being instead described by the noise covariance matrix (in this paper, Γ s ).Note that the CRBs provided in Section 4 are actually the square roots of CRBs, denoted [CRB(∆)] i,i 0.5 for parameter ∆ i (see Eq. ( 9)).As shown in Eq. ( 9), [CRB(∆)] i,i 0.5 has the same nature as a standard deviation and it is thus expressed in the same unit as ∆ i .From a practical point of view, given a set of ocean color parameters ∆ and a noise covariance matrix Γ s , [CRB(∆)] i,i 0.5 directly provides the minimum estimation standard deviation that can be possibly (but not necessarily) attained by any unbiased estimator ∆ i (r r s ).

Estimation of environmental noise
A reliable estimate of the environmental noise matrix Γ s based on real data is necessary to obtain consistent CRBs.In this paper, four estimates of Γ s were derived from images acquired by four satellite-borne multi-, super-and hyperspectral sensors, namely Sentinel-2 MSI (MultiSpectral Instrument), HICO (Hyperspectral Imager for the Coastal Ocean), Sentinel-3 OLCI (Ocean and Land Colour Instrument) and MODIS (Table 1).Besides having different radiometric configurations, these sensors cover a wide range of spatial resolutions, namely, 60 m for MSI (after downsampling the 10 and 20 m bands by averaging over square cells), 90 m for HICO, 300 m for OLCI full resolution and 1,000 m for MODIS reduced resolution.While the spatial resolutions and the radiometric configurations of MODIS and OLCI have been designed to monitor oceanic and coastal waters from regional to global scales, the high spatial resolution provided by MSI is convenient to monitor coastal and inland waters from local to regional scales.Incidentally, note that the inversion performances obtained using the original 10 m and/or 20 m resolution MSI bands (which are obviously noisier than the 60 m resolution bands used in this paper) cannot be better than the CRBs of MSI shown in Section 4.4.Finally, though being no longer operational, HICO was also included in this study because it provided both high spatial and spectral resolutions as well as a high SNR, all of which are critical to study shallow areas.Fig. 1.Environmental noise estimation based on the MSI, HICO, OLCI and MODIS (columns 1 to 4, resp.)images (row 1), where unreliable pixels (e.g., land and clouds) are masked.For each image, the white arrow in row 2 indicates the minimum of the modified ALCL-generated image (corresponding to the most homogeneous area), where the Γ s matrix (row 3) is estimated.In row 3, the size of each uniformly colored area denotes the sensor bandwidths at the corresponding wavebands (no spectral data measured in the black areas).Note that the colorbar is the same for each row.
The satellite images used for Γ s estimation were acquired over the Western Mediterranean region.They were directly obtained at Level-2A (Fig. 1).The MSI bottom-of-atmosphere reflectance image (dimensionless) was downloaded from the Copernicus Open Access Hub (scihub.copernicus.eu/dhus/#/home)and converted into remote-sensing reflectance by dividing it by a factor of π.Note that, although the Sentinel-2 Sen2cor toolbox used for atmospheric correction was not originally designed for aquatic scenes, several recent studies showed that it could allow satisfactory estimations of water reflectance [36,37].The HICO remote-sensing reflectance image was provided by the Oregon State University, College of Earth, Ocean, and Atmospheric Sciences (hico.coas.oregonstate.edu).The OLCI marine product (dimensionless) was downloaded from the Copernicus Online Data Access (coda.eumetsat.int)and converted into remote-sensing reflectance by dividing it by a factor of π.Finally, the MODIS remote-sensing reflectance image was obtained from the NASA Ocean Color data archive (oceancolor.gsfc.nasa.gov).
While sunglint was found to have a negligible influence on the considered OLCI and MODIS coarse spatial resolution data (especially due to the tilted viewing configuration of OLCI), it was corrected for MSI and HICO data by using the algorithm proposed in [38].For this purpose, note that, for both MSI and HICO sensors, the near-infrared reflectance used in the deglint algorithm was the band value at 842 nm for MSI, and the value averaged over the ten bands located between 800 and 850 nm for HICO.Estimating the sunglint contribution outside the spectral range where Γ s was derived, ensures that the environmental noise is not underestimated at all the considered wavebands (because the sunglint contribution is basically subtracted from all these wavebands).Finally, subsurface remote-sensing reflectance data were obtained after correcting for the air/water interface [23].
For each image, the Automated Local Convergence Locator (ALCL) method [39] was implemented to find a homogeneous area of optically-deep water in the image, where Γ s could then be objectively and reliably estimated (Fig. 1).ALCL consists in finding an area where the spectrally-averaged standard deviation of reflectance computed over a square cell does not change when increasing the cell size.The ALCL criterion to be minimized is the slope of the linear regression between the reflectance standard deviation and the cell size.In this paper, the ALCL criterion was multiplied by the spectrally-averaged standard deviation of reflectance at the maximum cell size in order to favor areas of low environmental noise.As shown in Fig. 1, such a modified ALCL criterion takes higher values around coastal areas, where the water column may be highly spatially variable.Note also the red squares on the modified ALCL-generated maps of MSI and HICO that indicate the presence of boats and boat-induced waves (Fig. 1).On the other hand, the modified ALCL criterion takes lower values over homogeneous deep water areas (Fig. 1).Implementing an ALCL-like algorithm is critical as it ensures that, in the selected area, the spectral variability is only due to the environmental noise and not to possible changes in bathymetry, bottom type and/or water clarity.
The HICO configuration is used in Section 4.2, Section 4.3 and Section 4.4, because, as a hyperspectral sensor, HICO can accurately characterize both deep and shallow waters due to its large number of narrow wavebands.The other sensor configurations are used in Section 4.4.

Ocean color remote-sensing methods
In this study, the minimum inversion uncertainties provided by CRBs are compared to the experimental estimation uncertainties obtained using two published methods [1,14,21].Both methods consist in iteratively optimizing a cost function that relates the measured reflectance r r s and the model µ = µ(∆) (Eq.( 1) and Eq. ( 8)) in order to derive the ocean color geophysical parameters ∆.The widely used Least-Square (LS) estimator [1,14] is given by The MILE (MaxImum Likelihood estimation including Environmental noise) estimator [21] utilizes the covariance matrix Γ s and is given by the minimum of the Mahalanobis distance: The bottom was defined as a linear mixture of 50% sand and 50% seagrasses, both reflectance spectra being extracted from the data set presented in [21].Finally, the solar and viewing zenith angles were set to 35 o and 0 o , respectively.For each depth value, the "mvnrnd" MATLAB function was used to randomly generate 2,000 noise-perturbed spectra using the model in Eq. ( 1) as the mean vector and the HICO Γ s matrix (Fig. 1) as the spectral covariance matrix.Each spectrum was then inverted using LS and MILE to estimate the unknown vector of parameters ∆ = [H, a φ (440), a g (440), b bp (550), B].Note that initialization was performed using a Latin Hypercube Sampling scheme as detailed in [21].For each depth value and each parameter ∆ i , the experimental standard deviations of LS and MILE estimates (hereafter referred to as LS and MILE inversion uncertainties) were calculated over the 2,000 samples and compared to [CRB(∆)] i,i 0.5 (Eq.( 9)).Note that the optimization domain was not bounded since setting lower and upper bounds amounts to introducing a priori knowledge within the inversion, which may make the comparison with CRBs misleading when the inversion uncertainty is high.Preliminary tests actually showed that bounding the optimization domain had a negligible effect on the retrieval, the only non-negligible influence being observed for H and B estimations over quasi optically deep waters.The experimental estimation biases (i.e., estimates of systematic errors) are not presented here as they were found to be mostly negligible for every parameter and both methods when H is lower than 10-15 m.Note that, for higher depths, compensations between H and B may lead to non-negligible biases (that increase with H) for both parameters [21], thus implying that the inversion uncertainty cannot be represented only by CRBs but also by biases.The results show that the CRBs obtained without a priori knowledge vary similarly to LS and MILE inversion uncertainties (Fig. 2).For example, due to the decreasing bottom detectability, the CRBs of H and B increase with depth and reach 17 and 85%, respectively, for H = 15 m.Also, the estimation of water clarity parameters (i.e., a φ (440) and b bp (550) in Fig. 2) becomes unreliable for very shallow waters as the bottom influence on water-leaving radiance overshadows those of water constituents [12,20].Note that this would be even more critical if the bottom intra-class variability (which is one of the most confounding factors for shallow water remote sensing [29]) could be taken into account in the modeling.The decrease in performance observed for such shallow waters is markedly more visible for b bp (550) estimation (whose CRB reaches 60% for H = 0.50 m) because the spectral signatures of sand and particle backscattering are similar, i.e., both increase the reflectance at all wavebands.
The results also show that the MILE inversion uncertainties almost equal the CRBs for every parameter and H ≤ 10 m, i.e., the MILE estimator is efficient [25].This result theoretically demonstrates that no method can perform significantly better than MILE if the bio-optical and probabilistic models are perfect.Note that such a conclusion may change when using actual spectra for which the water column may be not accurately described by the models.For H ≥ 10 m, progressively increasing biases may lead LS and MILE inversion uncertainties for H and B to increase with depth faster than the CRBs.

Optically deep waters
LS and MILE are also compared based on a synthetic HICO-like hyperspectral data set corresponding to optically deep waters.Unlike optically shallow waters for which depth and bottom cover can vary independently from water clarity, optically deep waters generally show correlations between phytoplankton, CDOM and suspended particle concentrations (especially in the open ocean).In this paper, these correlations are taken into account by assessing LS and MILE inversion uncertainties for the eight water optical types defined in [4].For each type, the mean subsurface remote-sensing reflectance spectrum provided in [4] was inverted using LS and the model in Eq. ( 8).This resulted in eight sets of a φ (440), a g (440) and b bp (550) values corresponding to an increasing turbidity, ranging from low-chlorophyll blue waters to turbid sediment-dominated waters (Table 2).For each of these eight sets, the same methodology as for shallow waters was then applied to generate 2,000 noise-perturbed spectra based on the HICO Γ s matrix, to invert the model (Eq.( 8)) and to compute the experimental standard deviations of a φ (440), a g (440) and b bp (550) estimates.Note that the experimental estimation bias was negligible (i.e., well below 10%) for the three parameters and the two methods, thus making CRBs appropriate indicators of minimum achievable inversion uncertainties.
Overall, for every water type, the CRBs obtained without a priori knowledge vary similarly to LS and MILE inversion uncertainties (Fig. 3).In the case of absorbing components (Figs.3(a) and 3(b)), the CRBs of a φ (440) and a g (440) depend on the a g (440)/a φ (440) ratio as noticed in [10].For example, for types 1 to 3, the total water absorption in the blue-green domain is dominated by CDOM for these simulations (Table 2), meaning that the CDOM influence overshadows that of phytoplankton [10].a g (440) (whose CRB ranges between 6 and 13%) can therefore be estimated more accurately than a φ (440) (whose CRB ranges between 16 and 30%).Also, the retrieval accuracies of a φ (440) and a g (440) increase with backscattering (see types 5 to 8 in Figs.3(a) and  3(b)).As also demonstrated, e.g., in [41], higher backscattering indeed leads to higher reflectance, and therefore, to higher SNR and improved estimation accuracies.
Particle backscattering (whose CRB is lower than 7% for every water type) can generally be retrieved more accurately than a φ (440) and a g (440) (Fig. 3(c)), which is in agreement with previous findings [10,11].For a given value of b bp (550), the CRB of b bp (550) increases with total absorption (see types 1 to 3 in Fig. 3(c)) due to the corresponding decreases in reflectance and SNR.On the contrary, b bp (550) can be estimated more accurately for higher b bp (550) values for the same reasons as for a φ (440) and a g (440), i.e., due to the increase in SNR (see types 5 to 8 in Fig. 3(c)).
Similarly to shallow waters, for the three parameters, the MILE inversion uncertainties are closer to the CRBs than are the LS ones, demonstrating the interest of taking spectral covariance into account when inverting the model based on HICO data [24].Note that this is more evident for b bp (550) estimation as also observed in [24].This may be due to the fact that the water-leaving radiance is affected more similarly over all wavebands by particle backscattering Table 2. a φ (440), a g (440) and b bp (550) estimated values for the water types defined in [4].
The results obtained without a priori knowledge for both optically shallow (Fig. 2) and deep waters (Fig. 3) show that, overall, the CRBs accurately describe the experimental estimation standard deviations obtained from simulated reflectance data.Therefore, no simulations and inversions will be performed in the following when investigating the influence of perfect a priori knowledge on model inversion (Section 4.3) and when presenting the CRBs obtained for MSI, OLCI and MODIS configurations (Section 4.4).

Using CRBs to assess the influence of a priori knowledge on model inversion
CRBs can be used to study how perfect a priori knowledge on some of the ocean color parameter(s) improves the estimation of remaining unknown parameters [26].For example, for shallow waters, depth may be known from a LiDAR-or acoustic-derived bathymetric map, and used as an input to the inversion algorithm.In the latter case, ∆ = [a φ (440), a g (440), b bp (550), B], so I F (∆) and CRB(∆) become 4 × 4 matrices.Note that such a study can also provide useful information on the interest of strictly bounding the optimization domain for one or several parameter(s) during model inversion (if the water column is sufficiently well known a priori), so as to significantly improve the estimation of remaining parameters.shows that, for the considered simulation conditions, using perfect a priori knowledge on H has little effect on a φ (440) retrieval.It has slightly more influence on b bp (550) estimation, e.g., at 5 m, the CRB is about two-fold lower than that obtained without a priori knowledge (Fig. 2(c)).However, the greatest gain in performance is obtained for B retrieval (Fig. 2(d)), whose CRB is seven-fold lower at 15 m when using perfect a priori knowledge on H.Such knowledge indeed prevents the compensations between H and B as mentioned in Section 4.2.2.Combining optical data with acoustic or LiDAR data thus appears to be relevant to improve the remote sensing of benthic habitats, confirming the results obtained in [42] based on real data.
Using perfect a priori knowledge on both H and B significantly improves the a φ (440) retrieval for very shallow waters only (e.g., by a factor 1.6 for H = 1 m), for which the bottom has the strongest influence (Fig. 2(b)).It has a much stronger influence on b bp (550) estimation, whose CRB at 1 m depth is seven-fold lower than that obtained without a priori knowledge (Fig. 2(c)).This improvement is mainly due to the similar spectral signatures of sandy bottoms and particle backscattering (Section 4.2.2):knowing the former significantly enhances the retrieval of the latter.Interestingly, this result theoretically proves the relevance of the Shallow Water Inversion Model [19] that uses known benthic albedo and bathymetric maps as inputs to the inversion algorithm to improve the remote sensing of water composition in shallow waters.
Ultimately, as expected, the highest increase in performance is obtained when all the parameters except the one to be estimated are known.

Optically deep waters
Figures 3(a) and 3(b) show that, for the eight water types, the CRB of a φ (440) (resp., a g (440)) is about two-fold lower when a g (440) (resp., a φ (440)) is known, and is close to the CRB obtained with known a g (440) and b bp (550) (resp., a φ (440) and b bp (550)).In other words, the accuracy of a φ (440) retrieval mainly depends on that of a g (440) retrieval and reciprocally.In particular, such a result indicates that the use of long ultraviolet bands (characterized by low chlorophyll absorption and high CDOM absorption) as suggested in [1] should be helpful in improving the decoupling and retrieval of both absorbing components.Figures 3(b) and 3(c) also show that the a g (440) and b bp (550) retrievals are uncorrelated, since knowledge on b bp (550) does not improve a g (440) estimation and reciprocally.CDOM absorption indeed affects the water-leaving radiance mainly in the blue-green domain while particle backscattering affects the radiance at all wavebands, which allows both effects to be accurately identified.

Using CRBs to predict minimum inversion uncertainties for MSI, HICO, OLCI and MODIS radiometric configurations
Beforehand, it is worth emphasizing that the CRBs of the four investigated sensors depend on various factors, including not only their radiometric configurations (e.g., number and position of wavebands, as well as their bandwiths) but also their spatial resolutions for example.Therefore, the effects of these factors on CRBs cannot be distinguished from each other based on the following results only.Furthermore, the environmental noise covariance matrices presented in Fig. 1 may be scene-specific, leading in turn to somewhat scene-specific CRBs as well.These important aspects must be kept in mind in the following when comparing the CRBs.

Optically shallow waters
In Fig. 4, CRBs are represented for the same depth range, water composition and bottom cover as in Fig.  areas.Similar conclusions can be drawn about the CRBs of OLCI, although its higher spatial resolution (300 m) may allow OLCI to characterize some large coastal areas.Actually, only the spatial resolutions of HICO (90 m) and MSI (60 m) are adapted for shallow water remote sensing.As a hyperspectral sensor, which has been shown to better identify the influences of all parameters than multispectral sensors [10,12], and despite its generally higher noise variances per band for the considered scenes (Fig. 1), HICO provides significantly lower CRBs than MSI (bearing in mind that HICO has a coarser spatial resolution).Such a result emphasizes the great potential of HICO-like sensors for accurately characterizing coastal areas.
Though not being specifically developed for shallow water remote sensing, MSI provides accurate inversion performances for H, B and, to a lesser extent, a φ (440), for H ≤ 5 m with CRBs generally lower than 25%.The high spatial and temporal resolutions offered by Sentinel-2 MSI thus make it an outstanding opportunity to monitor environments such as shallow coral reefs (e.g., to detect coral bleaching) that are an important subject of on-going research within the Sen2Coral ESA project, e.g., [43].Except for b bp (550), MSI generally shows high CRBs when H > 5 m.In this case, water absorption is such that the bottom is nearly invisible at the four wavebands beyond 650 nm (see the MSI wavebands in Fig. 1).For H > 5 m and the considered water and substrate compositions, the influences of all the parameters except b bp (550) can be mainly observed in the first three remaining wavebands.This means that, without a priori knowledge on H, a φ (440), a g (440) and/or B, the inversion is ill-posed and thus leads to poor estimation results for these four parameters.Meanwhile, b bp (550) can still be accurately retrieved because it also affects the reflectance at the last four wavebands.Beyond a given limiting optical depth, further constraining the inversion (e.g., by setting strict lower and upper bounds for some of the parameters to be retrieved) is thus required to improve the overall retrieval (note that the relevance of strictly bounding each parameter can be investigated as proposed in Section 4.3).It is worth emphasizing that such a limiting optical depth (and more generally, all the CRBs shown in Fig. 2 and Fig. 4) not only depends on depth, but also on water and substrate compositions.For example, the CRB of H could be shown to decrease when the water clarity and substrate brightness increase.

Optically deep waters
In Fig. 5 and Fig. 6, the CRBs of a φ (440) and b bp (550) are plotted for a wide variety of oceanic, coastal and inland waters.Note that this representation does not consider correlations between parameters, so some of the presented a φ (440) and b bp (550) combinations are unlikely to occur in natural waters, e.g., high a φ (440) and low b bp (550) [44].Yet, Fig. 5 and Fig. 6 provide a more comprehensive overview of the waters possibly encountered on Earth than the water types  defined in [4].Note that the ranges of CRB values obtained for Case-1 clear oceanic and Case-2 coastal and inland deep waters [45] (Fig. 5, Fig. 6) are further emphasized in Table 3.
The CRBs of a φ (440) and b bp (550) show similar patterns for the four sensors.For a given b bp (550) value, the CRB of a φ (440) (row 1 in Fig. 5 and Fig. 6) generally shows a bowl-shaped pattern with respect to a φ (440).It is minimum at the value a φ,min (440), with a φ,min (440) slightly greater than a g (440) (a φ,min (440) ≈ 0.03 m −1 for Fig. 5 and a φ,min (440) ≈ 0.15 m −1 for Fig. 6).For a φ (440) < a φ,min (440), the CRB of a φ (440) increases due to the overshadowing influence of CDOM (Section 4.2.3), while for a φ (440) > a φ,min (440), the increase is due to the decreasing SNR and progressive reflectance saturation.For a given b bp (550) value, the CRB of b bp (550) generally increases with a φ (440) and is less dependent on a g (440) than the CRB of a φ (440) (row 2 in Fig. 5 and Fig. 6).Finally, for a given a φ (440) value, the CRBs of a φ (440) and b bp (550) decrease as b bp (550) increases due to the associated increase in SNR (Section 4.2.3).
Because of the convenient choice of wavebands as well as its coarser spatial resolution leading to higher SNRs as compared to the other tested sensors (Fig. 1), MODIS mostly obtains the best CRBs for a φ (440) and b bp (550), i.e., generally below 25%, for the two a g (440) values (Fig. 5 and Fig. 6).For Case-1 waters, the CRBs are even lower than 8% (Fig. 5 and Table 3), thus demonstrating the strong potential of MODIS for mapping water clarity at the global scale.
Overall, for both a g (440) values, the CRBs of HICO and OLCI show rather similar patterns.For Case-1 waters, the CRBs obtained for these two sensors range between 5 and 23% for a φ (440), and between 2 and 11% for b bp (550) (Fig. 5 and Table 3).Their higher spatial resolutions as compared to MODIS also allows them to conveniently monitor Case-2 coastal and inland waters.For such waters, the CRBs of HICO and OLCI are, respectively, better than 22 and 10% for a φ (440), while they are, respectively, better than 17 and 10% for b bp (550) (Fig. 6 and Table 3).
These results thus emphasize the potential of these two sensors for characterizing a wide range of waters, ranging from oceanic clear waters to coastal and inland turbid waters.The comparable CRBs obtained for OLCI and HICO also suggest that a lower SNR (e.g., due to a finer spatial resolution) may be somewhat compensated for by a higher number of bands (potentially leading to a higher amount of information).
Overall, MSI provides higher CRBs compared to HICO, OLCI and MODIS.For example, the CRB of a φ (440) cannot be better than 27% for Case-1 waters, and even dramatically increases up to more than 200% for very low values of a φ (440) and b bp (550) (Fig. 5 and Table 3).These poorer performances are not only due to the finer spatial resolution of MSI (that limits the SNR), but also to its lower number of bands and coarser spectral resolution in the chlorophyll-sensitive 400-700 nm domain (Fig. 1).However, MSI enables a more accurate a φ (440) retrieval over Case-2 waters, the associated CRB ranging from 3 to 45% (Fig. 6 and Table 3).For such turbid waters, the red-edge MSI band at 705 nm is critical to achieve these performances as it enables the amplitude of the reflectance peak near 700 nm (that correlates with phytoplankton concentration [44]) to be captured.Note that the good performances obtained for Case-2 waters are in agreement with the recent results obtained in [46] using actual MSI data.Finally, it is worth noting that the MSI CRB for b bp (550) is better than 27% for both Case-1 and Case-2 waters (Fig. 5, Fig. 6 and Table 3).These results thus suggest that MSI could be used for detecting and tracking sediment plumes and algal blooms at fine spatial and temporal scales.This is especially important for monitoring coastal and inland waters that have a smaller geographic footprint compared to the open ocean [47].

Conclusions and perspectives
This study focuses on the development and application of Cramer-Rao bounds (CRBs) for ocean color remote sensing of deep and shallow waters.In both cases, the CRBs are computed using standard probabilistic and bio-optical models that represent the variability of water subsurface remote-sensing reflectance.Overall, CRBs are shown to provide valuable information about the best inversion performances that may be expected for a given sensor configuration and a given application.In this paper, CRBs are used to evaluate the efficiency of the LS and MILE inversion methods and to investigate to what extent the use of perfect a priori knowledge on one or several ocean color geophysical parameters can improve the estimation of remaining unknown parameters.In addition, CRBs are used to predict the minimum uncertainties in the inversion of ocean color calculated similarly for these three parameters.In the case of P, it is given by Equations ( 18)-( 21) also require us to calculate the derivatives of κ and u: Since b b does not depend on P (Eq.( 7)), we have ∂b b /∂P = 0.The last derivative is given by ∂a ∂P = a 0 + a 1 (1 + ln P).

4. 2 . 2 .
Optically shallow watersLS and MILE are compared based on a synthetic hyperspectral data set generated similarly as in[21,29,30].The wavebands considered here were those of HICO, as introduced in Section 4.1.The water column influence was studied for various depth values, ranging from 0.5 to 15 m by 1 m step.All the remaining model parameters were fixed.In particular, we considered an intermediate water clarity[40], i.e., a φ (440) = 0.05 m −1 , a g (440) = 0.1 m −1 , and b bp (550) = 0.01 m −1 .

Fig. 2 .
Fig. 2. Square roots of CRBs of (a) H, (b) a φ (440), (c) b bp (550) and (d) B (all in [%]) versus depth for optically shallow water and HICO configuration.Other model parameters are fixed: a φ (440) = 0.05 m −1 , a g (440) = 0.1 m −1 , b bp (550) = 0.01 m −1 , and the bottom is a mix of sand and seagrasses with B = 0.5.Various situations are considered: no a priori knowledge (Section 4.2.2), and perfect a priori knowledge on H, on both H and B, or on every parameter except the targeted one (Section 4.3.1).The LS and MILE experimental standard deviations obtained without a priori knowledge are also displayed.

Fig. 3 .
Fig. 3. Square roots of CRBs of (a) a φ (440), (b) a g (440) and (c) b bp (550) (all in [%]) for optically deep water and HICO configuration, and for the eight water types defined in [4] and corresponding to the parameter values shown in Table 2. Various situations are considered: no a priori knowledge (Section 4.2.3), and perfect a priori knowledge on a φ (440), on a g (440), on b bp (550), or on every parameter except the targeted one (Section 4.3.2).The LS and MILE experimental standard deviations obtained without a priori knowledge are also displayed.

4. 3 . 1 .
Figure 2(b)  shows that, for the considered simulation conditions, using perfect a priori knowledge on H has little effect on a φ (440) retrieval.It has slightly more influence on b bp (550) estimation, e.g., at 5 m, the CRB is about two-fold lower than that obtained without a priori knowledge (Fig.2(c)).However, the greatest gain in performance is obtained for B retrieval (Fig.2(d)), whose CRB is seven-fold lower at 15 m when using perfect a priori knowledge on H.Such knowledge indeed prevents the compensations between H and B as mentioned in Section 4.2.2.Combining optical data with acoustic or LiDAR data thus appears to be relevant to improve the remote sensing of benthic habitats, confirming the results obtained in[42] based on real data.Using perfect a priori knowledge on both H and B significantly improves the a φ (440) retrieval for very shallow waters only (e.g., by a factor 1.6 for H = 1 m), for which the bottom has the strongest influence (Fig.2(b)).It has a much stronger influence on b bp (550) estimation, whose CRB at 1 m depth is seven-fold lower than that obtained without a priori knowledge (Fig.2(c)).This improvement is mainly due to the similar spectral signatures of sandy bottoms and particle backscattering (Section 4.2.2):knowing the former significantly enhances the retrieval of the latter.Interestingly, this result theoretically proves the relevance of the Shallow Water Inversion Model[19] that uses known benthic albedo and bathymetric maps as inputs to the inversion algorithm to improve the remote sensing of water composition in shallow waters.Ultimately, as expected, the highest increase in performance is obtained when all the parameters except the one to be estimated are known.

Table 1 .
Features of the remote-sensing images used to estimate the environmental noise.
2. Overall, MODIS obtains lower CRBs than the other sensors for H, a φ (440), b bp (550) and B. While this result indicates that MODIS can theoretically remotely sense shallow waters with a satisfactory accuracy, the low spatial resolution considered here (1 km) generally makes this sensor inappropriate for such applications due to the difficulty in properly focusing on coastal