Comparison of three SeaWiFS atmospheric correction algorithms for turbid waters using AERONET-OC measurements

The use of satellites to monitor the color of the ocean requires effective removal of the atmospheric signal. This can be performed by extrapolating the aerosol optical properties in the visible from the near-infrared (NIR) spectral region assuming that the seawater is totally absorbant in this latter part of the spectrum. However, the non-negligible water-leaving radiance in the NIR which is characteristic of turbid waters may lead to an overestimate of the atmospheric radiance in the whole visible spectrum with increasing severity at shorter wavelengths. This may result in significant errors, if not complete failure, of various algorithms for the retrieval of chlorophyll-a concentration, inherent optical properties and biogeochemical parameters of surface waters.This paper presents results of an inter-comparison study of three methods that compensate for NIR water-leaving radiances and that are based on very different hypothesis: 1) the standard SeaWiFS algorithm (Stumpf et al., 2003; Bailey et al., 2010) based on a bio-optical model and an iterative process; 2) the algorithm developed by Ruddick et al. (2000) based on the spatial homogeneity of the NIR ratios of the aerosol and water-leaving radiances; and 3) the algorithm of Kuchinke et al. (2009) based on a fully coupled atmosphere-ocean spectral optimization inversion. They are compared using normalized water-leaving radiance nL_w in the visible. The reference source for comparison is ground-based measurements from three AERONET-Ocean Color sites, one in the Adriatic Sea and two in the East Coast of USA.Based on the matchup exercise, the best overall estimates of the nL_w are obtained with the latest SeaWiFS standard algorithm version with relative error varying from 14.97% to 35.27% for λ = 490 nm and λ = 670 nm respectively. The least accurate estimates are given by the algorithm of Ruddick, the relative errors being between 16.36% and 42.92% for λ = 490 nm and λ = 412 nm, respectively. The algorithm of Kuchinke appears to be the most accurate algorithm at 412 nm (30.02%), 510 (15.54%) and 670 nm (32.32%) using its default optimization and bio-optical model coefficient settings.Similar conclusions are obtained for the aerosol optical properties (aerosol optical thickness τ(865) and the Ångström exponent, α(510, 865)). Those parameters are retrieved more accurately with the SeaWiFS standard algorithm (relative error of 33% and 54.15% for τ(865) and α(510, 865)).A detailed analysis of the hypotheses of the methods is given for explaining the differences between the algorithms. The determination of the aerosol parameters is critical for the algorithm of Ruddick et al. (2000) while the bio-optical model is critical for the algorithm of Stumpf et al. (2003) utilized in the standard SeaWiFS atmospheric correction and both aerosol and bio-optical model for the coupled atmospheric-ocean algorithm of Kuchinke. The Kuchinke algorithm presents model aerosol-size distributions that differ from real aerosol-size distribution pertaining to the measurements. In conclusion, the results show that for the given atmospheric and oceanic conditions of this study, the SeaWiFS atmospheric correction algorithm is most appropriate for estimating the marine and aerosol parameters in the given turbid waters regions.     

A three component classification of phytoplankton absorption spectra: Application to ocean-color data

A method is presented to identify absorption characteristics of three optically-distinct phytoplankton classes from a suite of measurements of total phytoplankton absorption coefficient and chlorophyll-a concentration by successive application of the two-population absorption model of Sathyendranath et al. (2001) and Devred et al. (2006a). The total phytoplankton absorption coefficient at multiple wavelengths is expressed as the weighted sum of the absorption coefficients of each class at those wavelengths. The resultant system of equations is solved under some constraints to derive the fraction of each class present in any given sample of seawater, given the spectrum of total phytoplankton absorption coefficient. When applied to a large database, the results compare well with phytoplankton size-classes derived from pigment composition, so that we can assume that the three phytoplankton classes derived from absorption coefficients are representative of the pico-, nano- and microphytoplankton size classes. A modification is proposed to the pigment-based phytoplankton size classification of Uitz et al. (2006) to account for the effect of fucoxanthin associated with nanophytoplankton. Comparison between satellite and in situ data demonstrates the potential of satellite ocean-color data to yield the distribution of phytoplankton size classes from space. The algorithm is applied to phytoplankton absorption coefficients derived from remotely-sensed reflectance values collected by SeaWiFS over the Northwest Atlantic in 2007. Monthly composites for April, August and November, representative of Spring, Summer and Fall, give synoptic views of the phytoplankton community structure: a Spring bloom dominated by microphytoplankton is followed by a second, less intense, bloom in the Fall dominated by nanophytoplankton. Picophytoplankton are dominant in the study area in Summer.     

Cross-site consistent in situ measurements for satellite ocean color applications: The BiOMaP radiometric dataset

The development and assessment of satellite ocean color products require quality assured in situ data representative of the variety of bio-optical regimes encountered in the different seas. The measurement program named Bio-Optical mapping of Marine Properties (BiOMaP) fulfills this requirement by using identical instruments and applying cross-site consistent methods for the characterization of seawater inherent and apparent optical properties in the various European seas. This work introduces the BiOMaP radiometric data and describes their application to the validation of primary ocean color products. Within this framework, the radiometric data are discussed through the spectral shape and amplitude of normalized water-leaving radiances (L_WN). Specifically, the spectral shape is expressed through the Principal Component Analysis of L_WN(λ)/L_WN(555) while the amplitude is represented by L_WN(555). The resulting distribution of BiOMaP data in a three dimensional feature space demonstrates a continuity of cases across the investigated marine regions confirming a wide representativity of bio-optical regimes. The application of BiOMaP data to the validation of remote sensing reflectance from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) and the Moderate Resolution Imaging Spectroradiometer (MODIS), indicates improved performance of the SeaWiFS Data Analysis System (SeaDAS, version 6.1) atmospheric correction. In particular, the comparison of satellite and in situ matchups in the blue spectral region shows biases of a few percent with respect to the much larger reported in studies relying on earlier SeaDAS versions. Matchup analyses, restricted to the Eastern Mediterranean, Black and Baltic Seas, indicate marked regional differences likely explained by the diversity of water and aerosol types.     

Ocean sand ridge signatures in the Bohai Sea observed by satellite ocean color and synthetic aperture radar measurements

Satellite measurements from Synthetic Aperture Radar (SAR) and the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Aqua platform are used to study the ocean sand ridges in the eastern Bohai Sea in China. Even though the imaging mechanisms for SAR and MODIS-Aqua remote sensing are different, the sand ridges are shown to have exactly the same patterns in images from both sensors. Therefore, the location, extension and coverage of the ocean sand ridges can be detected and cross-examined by both SAR and MODIS-Aqua observations. Satellite images show quite different sand ridge distribution pattern from the published bathymetry map (based on in situ data) that shows six sand ridges in the area. 10 finger-shaped sand ridges are identified from satellite observations. The tidal-current/sand-ridge interaction driven physical and optical changes are assessed and evaluated. The existence of sand ridges causes enhanced water diffuse attenuation coefficient K_d(490) and elevated normalized water-leaving radiance at the red and near-infrared (NIR) wavelengths. The sea surface over the sand ridges experiences significant seasonal variability of water turbidity and shows remarkable differences from nearby ocean regions. During winter, K_d(490) values are about 2-3 m^− 1 over the ridges, while the maximum K_d(490) in the neighboring oceans is approximately 1.5 m^− 1. In summer, the enhancement of the sea surface turbidity is less significant than that which occurs in winter.     

Modeling spectral reflectance of optically complex waters using bio-optical measurements from Tokyo Bay

This study presents an approach for optimally parameterizing a reflectance model. A parameterization scheme is realized based on a comprehensive bio-optical data set, including subsurface downwelling and upwelling irradiance spectra, absorption spectra of particle and dissolved substances, as well as chlorophyll and total suspended matter concentrations at 45 stations near Tokyo Bay between 1982 and 1984. The irradiance reflectance model is implemented with three-component inherent optical property submodels.In this parameterization scheme, an unsupervised classification was applied in the hyper-spectral space of reflectance, leading to three spectrally distinct optical water types. The reflectance model was parameterized for the entire data set, and then parameterized for each of the water types. The three sets of type-specific model parameters, which define corresponding IOP submodels, are believed to accommodate differences in the optical properties of the in-water constituents. The parameterized reflectance model was evaluated by both reconstructing measured reflectance spectra and solving for the nonlinear inverse problem to retrieve in-water constituent concentrations. The model accuracy was significantly improved in the forward direction for classified waters over that of non-classified waters, but no significant improvement was achieved in the retrieval accuracy (inverse direction). A larger data set with greater resolution of constituent inherent optical properties would likely improve the modeling results.     

Evaluation of MODIS SWIR and NIR-SWIR atmospheric correction algorithms using SeaBASS data

Using the NASA maintained ocean optical and biological in situ data that were collected during 2002-2005, we have evaluated the performance of atmospheric correction algorithms for the ocean color products from the Moderate Resolution Imaging Spectroradiometer (MODIS) on Aqua. Specifically, algorithms using the MODIS shortwave infrared (SWIR) bands and an approach using the near-infrared (NIR) and SWIR combined method are evaluated, compared to the match-up results from the NASA standard algorithm (using the NIR bands). The in situ data for the match-up analyses were collected mostly from non-turbid ocean waters. It is critical to assess and understand the algorithm performance for deriving MODIS ocean color products, providing science and user communities with the important data quality information. Results show that, although the SWIR method for data processing has generally reduced the bias errors, the noise errors are increased due mainly to significantly lower sensor signal-noise ratio (SNR) values for the MODIS SWIR bands, as well as the increased uncertainties using the SWIR method for the atmospheric correction. This has further demonstrated that future ocean color satellite sensors will require significantly improved sensor SNR performance for the SWIR bands. The NIR-SWIR combined method, for which the non-turbid and turbid ocean waters are processed using the NIR and SWIR method, respectively, has been shown to produce improved ocean color products.     

Evaluation of shortwave infrared atmospheric correction for ocean color remote sensing of Chesapeake Bay

The NASA Moderate Resolution Imaging Spectroradiometer onboard the Aqua platform (MODIS-Aqua) provides a viable data stream for operational water quality monitoring of Chesapeake Bay. Marine geophysical products from MODIS-Aqua depend on the efficacy of the atmospheric correction process, which can be problematic in coastal environments. The operational atmospheric correction algorithm for MODIS-Aqua requires an assumption of negligible near-infrared water-leaving radiance, nL_w(NIR). This assumption progressively degrades with increasing turbidity and, as such, methods exist to account for non-negligible nL_w(NIR) within the atmospheric correction process or to use alternate radiometric bands where the assumption is satisfied, such as those positioned within shortwave infrared (SWIR) region of the spectrum. We evaluated a decade-long time-series of nL_w(λ) from MODIS-Aqua in Chesapeake Bay derived using NIR and SWIR bands for atmospheric correction. Low signal-to-noise ratios (SNR) for the SWIR bands of MODIS-Aqua added noise errors to the derived radiances, which produced broad, flat frequency distributions of nL_w(λ) relative to those produced using the NIR bands. The SWIR approach produced an increased number of negative nL_w(λ) and decreased sample size relative to the NIR approach. Revised vicarious calibration and regional tuning of the scheme to switch between the NIR and SWIR approaches may improve retrievals in Chesapeake Bay, however, poor SNR values for the MODIS-Aqua SWIR bands remain the primary deficiency of the SWIR-based atmospheric correction approach.     

Software for geologic processing of landsat imagery

The effective delineation of rock boundaries and other geologic features from LANDSAT imagery usually requires the use of scenes that have been processed digitally to optimize information content and presentation. Geological processing programs correct for atmospheric distortions and standardize reflectances to uniform ranges for within and between scene continuity. Enhancement procedures then are applied to increase total image resolution. Once the images have been enhanced suitably, secondary procedures such as ratioing, density slicing, and edge detection can be used to highlight specific image properties. The ideal geological software system for LANDSAT processing is a combination of programs that can be ordered to produce the desired resolution of reflectance zones while retaining the structural and texural information that is vital to accurate interpretation.     

Development of a coupled land-atmosphere satellite data assimilation system for improved local atmospheric simulations

This study developed a coupled land-atmosphere satellite data assimilation system as a new physical downscaling approach, by coupling a mesoscale atmospheric model with a land data assimilation system (LDAS). The LDAS consists of a land surface scheme as the model operator, a radiative transfer model as the observation operator, and the simulated annealing method for minimizing the difference between the observed and simulated microwave brightness temperature. The atmospheric model produces forcing data for the LDAS, and the LDAS produces better initial surface conditions for the modelling system. This coupled system can take into account land surface heterogeneities through assimilating satellite data for a better precipitation prediction. To assess the effectiveness of the new system, 3-dimensional numerical experiments were carried out in a mesoscale area of the Tibetan Plateau during the wet monsoon season. The results show significant improvement compared with a no assimilation regional atmospheric model simply nested from the global model. The surface soil moisture content and its distribution from the assimilation system were more consistent to in situ observations. These better surface conditions affect the land-atmosphere interactions through convection systems and lead to better atmospheric predictability as confirmed by satellite-based cloud observations and in situ sounding observations. Through the use of satellite brightness temperature, the developed coupled land-atmosphere assimilation system has shown potential ability to provide better initial surface conditions and its inputs to the atmosphere and to improve physical downscaling through regional models.     

An improved in-situ bio-optical data set for ocean color algorithm development and satellite data product validation

Global satellite ocean color instruments provide the scientific community a high-resolution means of studying the marine biosphere. Satellite data product validation and algorithm development activities both require the substantial accumulation of high-quality in-situ observations. The NASA Ocean Biology Processing Group maintains a local repository of in-situ marine bio-optical data, the SeaWiFS Bio-optical Archive and Storage System (SeaBASS), to facilitate their ocean color satellite validation analyses. Data were acquired from SeaBASS and used to compile a large set of coincident radiometric observations and phytoplankton pigment concentrations for use in bio-optical algorithm development. This new data set, the NASA bio-Optical Marine Algorithm Data set (NOMAD), includes over 3400 stations of spectral water-leaving radiances, surface irradiances, and diffuse downwelling attenuation coefficients, encompassing chlorophyll a concentrations ranging from 0.012 to 72.12 mg m^− 3. Metadata, such as the date, time, and location of data collection, and ancillary data, including sea surface temperatures and water depths, accompany each record. This paper describes the assembly and evaluation of NOMAD, and further illustrates the broad geophysical range of stations incorporated into NOMAD.     

Multiobjective PSO based adaption of neural network topology for pixel classification in satellite imagery

The proposed work involves the multiobjective PSO based adaption of optimal neural network topology for the classification of multispectral satellite images. It is per pixel supervised classification using spectral bands (original feature space). This paper also presents a thorough experimental analysis to investigate the behavior of neural network classifier for given problem. Based on 1050 number of experiments, we conclude that following two critical issues needs to be addressed: (1) selection of most discriminative spectral bands and (2) determination of optimal number of nodes in hidden layer. We propose new methodology based on multiobjective particle swarm optimization (MOPSO) technique to determine discriminative spectral bands and the number of hidden layer node simultaneously. The accuracy with neural network structure thus obtained is compared with that of traditional classifiers like MLC and Euclidean classifier. The performance of proposed classifier is evaluated quantitatively using Xie-Beni and β indexes. The result shows the superiority of the proposed method to the conventional one.     

Modified differential evolution based fuzzy clustering for pixel classification in remote sensing imagery

The problem of classifying an image into different homogeneous regions is viewed as the task of clustering the pixels in the intensity space. In particular, satellite images contain landcover types some of which cover significantly large areas, while some (e.g., bridges and roads) occupy relatively much smaller regions. Detecting regions or clusters of such widely varying sizes presents a challenging task. A modified differential evolution based fuzzy clustering technique, is proposed in this article. Real-coded encoding of the cluster centres is used for this purpose. Results demonstrating the effectiveness of the proposed technique are provided for several synthetic and real life data sets as well as for some benchmark functions. Different landcover regions in remote sensing imagery have also been classified using the proposed technique to establish its efficiency. Statistical significance tests have been performed to establish the superiority of the proposed algorithm.     

A sun glint correction method for hyperspectral imagery containing areas with non-negligible water leaving NIR signal

Present sun glint removal methods overcorrect data collected in very shallow (less than 2 m) waters if the sensors used do not have bands in far infrared part of the spectrum. The reason is assuming of zero water leaving signal at near infrared (NIR) wavelengths. This assumption is not valid in very shallow waters, but also in areas where aquatic vegetation reaches water surface and in case of phytoplankton blooms that reach very high biomass or form surface scum. We propose an alternative method that can be used for successful glint removal even if the sensor does not have spectral bands beyond 800 nm. The proposed method utilises the presence and depth of the oxygen absorption feature near 760 nm as an indicator of glint contamination. This method allows removing sun glint from hyperspectral imagery preserving shape and magnitude of reflectance spectra in the cases where the negligible water leaving NIR signal is not valid.     

Validation of SeaWiFS chlorophyll a concentrations in the Southern Ocean: A revisit

Surface chlorophyll a concentrations (C_a, mg m^− 3) in the Southern Ocean estimated from SeaWiFS satellite data have been reported in the literature to be significantly lower than those measured from in situ water samples using fluorometric methods. However, we found that high-resolution (∼ 1 km²/pixel) daily SeaWiFS C_a (C_a^SWF) data (SeaDAS4.8, OC4v4 algorithm) was an accurate measure of in situ C_a during January-February of 1998-2002 if concurrent in situ data measured by HPLC (C_a^HPLC) instead of fluorometric (C_a^Fluor) measurements were used as ground truth. Our analyses indicate that C_a^Fluor is 2.48 ± 2.23 (n = 647) times greater than C_a^HPLC between 0.05 and 1.5 mg m^− 3 and that the percentage overestimation of in situ C_a by fluorometric measurements increases with decreasing concentrations. The ratio of C_a^SWF/C_a^HPLC is 1.12 ± 0.91 (n = 96), whereas the ratio of C_a^SWF/C_a^Fluor is 0.55 ± 0.63 (n = 307). Furthermore, there is no significant bias in C_a^SWF (12% and − 0.07 in linear and log-transformed C_a, respectively) when C_a^HPLC is used as ground truth instead of C_a^Fluor. The high C_a^Fluor/C_a^HPLC ratio may be attributed to the relatively low concentrations of chlorophyll b (C_b/C_a = 0.023 ± 0.034, n = 482) and relatively high concentrations of chlorophyll c (C_c/C_a = 0.25 ± 0.59, n = 482) in the phytoplankton pigment composition when compared to values from other regions. Because more than 90% of the waters in the study area, as well as in the entire Southern Ocean (south of 60° S), have C_a^SWF between 0.05 and 1.5 mg m^− 3, we consider that the SeaWiFS performance of C_a retrieval is satisfactory and for this C_a range there is no need to further develop a “regional” bio-optical algorithm to account for the previous SeaWiFS “underestimation”.     

An intercomparison of bio-optical techniques for detecting dominant phytoplankton size class from satellite remote sensing

Satellite remote sensing of ocean colour is the only method currently available for synoptically measuring wide-area properties of ocean ecosystems, such as phytoplankton chlorophyll biomass. Recently, a variety of bio-optical and ecological methods have been established that use satellite data to identify and differentiate between either phytoplankton functional types (PFTs) or phytoplankton size classes (PSCs). In this study, several of these techniques were evaluated against in situ observations to determine their ability to detect dominant phytoplankton size classes (micro-, nano- and picoplankton). The techniques are applied to a 10-year ocean-colour data series from the SeaWiFS satellite sensor and compared with in situ data (6504 samples) from a variety of locations in the global ocean. Results show that spectral-response, ecological and abundance-based approaches can all perform with similar accuracy. Detection of microplankton and picoplankton were generally better than detection of nanoplankton. Abundance-based approaches were shown to provide better spatial retrieval of PSCs. Individual model performance varied according to PSC, input satellite data sources and in situ validation data types. Uncertainty in the comparison procedure and data sources was considered. Improved availability of in situ observations would aid ongoing research in this field.     

Survival analysis of a neotropical rainforest using multitemporal satellite imagery

We present results of an analysis of deforestation at a UNESCO Biosphere Reserve, the Parque National Yasuní, located in the rainforests of eastern Ecuador using multitemporal Landsat TM and ETM+ satellite imagery. Using survival analysis, we assessed both current and future trends in deforestation rates, and investigated the impact of spatial, cultural, and economic factors on deforestation. These factors included the distance from roads, rivers, research facilities, oil facilities, markets and towns, and land ownership by colonists, native inhabitants, and an oil company. We found the annual rate of deforestation is currently only 0.11%, but that this rate is increasing with time and, assuming that the trend of increasing rate of forest loss continues, we would predict that by 2063, 50% of the forest within 2 km of an oil access road will be lost to unhindered colonization and anthropogenic conversion. The Quechua colonists are associated with areas of the highest rate of deforestation, followed by the native Huaorani and the lowest region of deforestation was in areas occupied by a local oil company. By far, the strongest predictor of where deforestation is predicted to occur was proximity to the road. Proximity to research sites, oil facilities, market, and rivers significantly decreases deforestation rates, and proximity to towns significantly increases deforestation rates.     

Monitoring turbidity in Tampa Bay using MODIS/Aqua 250-m imagery

We developed an approach to map turbidity in estuaries using a time series (May 2003 to April 2006) of 250-m resolution images from the Moderate Resolution Imaging Spectroradiometer (MODIS) on NASA's Aqua satellite, using Tampa Bay as a case study. Cross-calibration of the MODIS 250-m data (originally designed for land use) with the well-calibrated MODIS 1-km ocean data showed that the pre-launch radiometric calibration of the 250-m bands was adequate. A simple single scattering atmospheric correction provided reliable retrievals of remote sensing reflectance at 645 nm (0.002 < R_rs(645) < 0.015 sr^− 1, median bias = − 7%, slope = 0.95, intercept = 0.00, r² = 0.97, n = 15). A more rigorous approach, using a multiple scattering atmospheric correction of the cross-calibrated at-sensor radiances, retrieved similar R_rs(645). R_rs(645) estimates, after stringent data quality control, showed a close correlation with in situ turbidity (turbidity = 1203.9 × R_rs(645)^1.087, 0.9 < turbidity < 8.0 NTU, r² = 0.73, n = 43). MODIS turbidity imagery derived using the developed approach showed that turbidity in Hillsborough Bay (HB) was consistently higher than that in other sub-regions except in August and September, when higher concentrations of colored dissolved organic matter seem to have caused underestimates of turbidity. In comparison, turbidity in Middle Tampa Bay (MTB) was generally lowest among the Bay throughout the year. Both Old Tampa Bay (OTB) and Low Tampa Bay (LTB) showed marked seasonal variations with higher turbidity in LTB during the dry season and in OTB during the wet season, respectively. This seasonality is linked to wind-driven bottom resuspension events in lower portion of the Bay and river inputs of sediments in the upper portion of the Bay. The Bay also experiences significant interannual variation in turbidity, which was attributed primarily to changes in wind forcing. Compared with the once-per-month, non-synoptic in situ surveys, synoptic and frequent sampling facilitated by satellite remote sensing provides improved assessments of turbidity patterns and thus a valuable tool for operational monitoring of water quality of estuarine and coastal waters such as in Tampa Bay.     

Estimating the concentration of optically active constituents of sea water by Takagi-Sugeno models with quadratic rule consequents

Determining the concentrations of chlorophyll, suspended particulate matter and coloured dissolved organic matter in the sea water is basic to support the monitoring of upwelling phenomena, algae blooms, and changes in the marine ecosystem. Since these concentrations affect the spectral distribution of the solar light back-scattered by the water body, their estimation can be computed by using a set of remotely sensed multispectral measurements of the reflected sunlight. In this paper, the relation between the concentrations of interest and the average subsurface reflectances is modelled by means of a set of second-order Takagi-Sugeno (TS) fuzzy rules. Unlike first-order TS rules, which adopt linear functions as consequent, second-order TS rules exploit quadratic functions, thus improving the modelling capability of the rule in the subspace determined by the antecedent. First, we show how we can build a second-order TS model through a simple transformation, which allows estimating the consequent parameters using standard linear least-squares algorithms, and by adopting one of the most used methods proposed in the literature to generate first-order TS models. Then, we compare first-order and second-order TS models against mean square error and interpretability of rules. We highlight how second-order TS models allow us to achieve better approximation than first-order TS models though maintaining interpretability of the rules. Finally, we show how second-order TS models perform considerably better (the mean square error is lower by two orders of magnitude) than the specific implementations of radial basis function networks and multi-layer perceptron networks used in previous papers for the same application domain.