Sensitivity analysis of the dark spectrum fitting atmospheric correction for metre- and decametre-scale satellite...

The performance of the dark spectrum fitting (DSF) atmospheric correction algorithm is evaluated using matchups between metre- and decametre-scale satellite imagery as processed with ACOLITE and measurements from autonomous PANTHYR hyperspectral radiometer systems deployed in the Adriatic and North Sea. Imagery from the operational land imager (OLI) on Landsat 8, the multispectral instrument (MSI) on Sentinel-2 A and B, and the PlanetScope CubeSat constellation was processed for both sites using a fixed atmospheric path reflectance in a small region of interest around the system's deployment location, using a number of processing settings, including a new sky reflectance correction. The mean absolute relative differences (MARD) between in situ and satellite measured reflectances reach <20% in the Blue and 11% in the Green bands around 490 and 560 nm for the best performing configuration for MSI and OLI. Higher relative errors are found for the shortest Blue bands around 440 nm (30-100% MARD), and in the Red-Edge and near-infrared bands (35-100% MARD), largely influenced by the lower absolute data range in the observations. Root mean squared differences (RMSD) increase from 0.005 in the NIR to about 0.015-0.020 in the Blue band, consistent with increasing atmospheric path reflectance. Validation of the Red-Edge and NIR bands on Sentinel-2 is presented, as well as for the first time, the Panchromatic band (17-26% MARD) on Landsat 8, and the derived Orange contra-band (8-33% MARD for waters in the algorithm domain, and around 40-80% MARD overall). For Sentinel-2, excluding the SWIR bands from the DSF gave better performances, likely due to calibration issues of MSI at longer wavelengths. Excluding the SWIR on Landsat 8 gave good performance as well, indicating robustness of the DSF to the available band set. The DSF performance was found to be rather insensitive to (1) the wavelength spacing in the lookup tables used for the atmospheric correction, (2) the use of default or ancillary information on gas concentration and atmospheric pressure, and (3) the size of the ROI over which the path reflectance is estimated. The performance of the PlanetScope constellation is found to be similar to previously published results, with the standard DSF giving the best results in the visible bands in terms of MARD (24-40% overall, and 18-29% for the turbid site). The new sky reflectance correction gave mixed results, although it reduced the mean biases for certain configurations and improved results for the processing excluding the SWIR bands, giving lower RMSD and MARD especially at longer wavelengths (>600 nm). The results presented in this article should serve as guidelines for general use of ACOLITE and the DSF.