Impact of computational methods and spectral models on the retrieval of optical properties via spectral optimization

Spectral optimization algorithm (SOA) is a well-accepted scheme for the retrieval of water constituents from the measurement of ocean color radiometry. It defines an error function between the input and output remote sensing reflectance spectrum, with the latter modeled with a few variables that represent the optically active properties, while the variables are solved numerically by minimizing the error function. In this paper, with data from numerical simulations and field measurements as input, we evaluate four computational methods for minimization (optimization) for their efficiency and accuracy on solutions, and illustrate impact of bio-optical models on the retrievals. The four optimization routines are the Levenberg-Marquardt (LM), the Generalized Reduced Gradient (GRG), the Downhill Simplex Method (Amoeba), and the Simulated Annealing-Downhill Simplex (i.e. SA + Amoeba, hereafter abbreviated as SAA). The Garver-Siegel-Maritorena SOA model is used as a base to test these computational methods. It is observed that 1) LM is the fastest method, but SAA has the largest number of valid retrievals; 2) the quality of final solutions are strongly influenced by the forms of spectral models (or eigen functions); and 3) dynamically-varying eigen functions are necessary to obtain smaller errors for both reflectance spectrum and retrievals. Results of this study provide helpful guidance for the selection of a computational method and spectral models if an SOA scheme is to be used to process ocean color images.