Remote sensing of the ocean contributions from ultraviolet to near-infrared using the shortwave infrared bands: simulations

In the remote sensing of the ocean near-surface properties, it is essential to derive accurate water-leaving radiance spectra through the process of the atmospheric correction. The atmospheric correction algorithm for Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) and Moderate Resolution Imaging Spectroradiometer (MODIS) uses two near-infrared (NIR) bands at 765 and 865 nm (748 and 869 nm for MODIS) for retrieval of aerosol properties with assumption of the black ocean at the NIR wavelengths. Modifications are implemented to account for some of the NIR ocean contributions for the productive but not very turbid waters. For turbid waters in the coastal regions, however, the ocean could have significant contributions in the NIR, leading to significant errors in the satellite-derived ocean water-leaving radiances. For the shortwave infrared (SWIR) wavelengths (approximately > 1000 nm), water has significantly larger absorption than those for the NIR bands. Thus the black ocean assumption at the SWIR bands is generally valid for turbid waters. In addition, for future sensors, it is also useful to include the UV bands to better quantify the ocean organic and inorganic materials, as well as for help in atmospheric correction. Simulations are carried out to evaluate the performance of atmospheric correction for nonabsorbing and weakly absorbing aerosols using the NIR bands and various combinations of the SWIR bands for deriving the water-leaving radiances at the UV (340 nm) and visible wavelengths. Simulations show that atmospheric correction using the SWIR bands can generally produce results comparable to atmospheric correction using the NIR bands. In particular, the water-leaving radiance at the UV band (340 nm) can also be derived accurately. The results from a sensitivity study for the required sensor noise equivalent reflectance, (NE Delta rho), [or the signal-to-noise ratio (SNR)] for the NIR and SWIR bands are provided and discussed.     

Improving atmospheric correction for highly productive coastal waters using the short wave infrared retrieval algorithm with water-leaving reflectance constraints at 412 nm

The recently developed short wave infrared (SWIR) atmospheric correction algorithm for ocean color retrieval uses long wavelength channels to retrieve atmospheric parameters to avoid bright pixel contamination. However, this retrieval is highly sensitive to errors in the aerosol model, which is magnified by the higher variability of aerosols observed over urban coastal areas. While adding extra regional aerosol models into the retrieval lookup tables would tend to increase retrieval error since these models are hard to distinguish in the IR, we explore the possibility that for highly productive waters with high colored dissolved organic matter, an estimate of the 412 nm channel water-leaving reflectance can be used to constrain the aerosol model retrieval and improve the water-leaving reflectance retrieval. Simulations show that this constraint is particularly useful where aerosol diversity is significant. To assess this algorithm we compare our retrievals with the operational SeaWiFS Data Analysis System (SeaDAS) SWIR and near infrared retrievals using in situ validation data in the Chesapeake Bay and show that, especially for absorbing aerosols, significant improvement is obtained. Further insight is also obtained by the intercomparison of retrieved remote sensing reflectance images at 443 and 551 nm, which demonstrates the removal of anomalous artifacts in the operational SeaDAS retrieval.     

Assessment of a bidirectional reflectance distribution correction of above-water and satellite water-leaving radiance in coastal waters

Water-leaving radiances, retrieved from in situ or satellite measurements, need to be corrected for the bidirectional properties of the measured light in order to standardize the data and make them comparable with each other. The current operational algorithm for the correction of bidirectional effects from the satellite ocean color data is optimized for typical oceanic waters. However, versions of bidirectional reflectance correction algorithms specifically tuned for typical coastal waters and other case 2 conditions are particularly needed to improve the overall quality of those data. In order to analyze the bidirectional reflectance distribution function (BRDF) of case 2 waters, a dataset of typical remote sensing reflectances was generated through radiative transfer simulations for a large range of viewing and illumination geometries. Based on this simulated dataset, a case 2 water focused remote sensing reflectance model is proposed to correct above-water and satellite water-leaving radiance data for bidirectional effects. The proposed model is first validated with a one year time series of in situ above-water measurements acquired by collocated multispectral and hyperspectral radiometers, which have different viewing geometries installed at the Long Island Sound Coastal Observatory (LISCO). Match-ups and intercomparisons performed on these concurrent measurements show that the proposed algorithm outperforms the algorithm currently in use at all wavelengths, with average improvement of 2.4% over the spectral range. LISCO's time series data have also been used to evaluate improvements in match-up comparisons of Moderate Resolution Imaging Spectroradiometer satellite data when the proposed BRDF correction is used in lieu of the current algorithm. It is shown that the discrepancies between coincident in-situ sea-based and satellite data decreased by 3.15% with the use of the proposed algorithm. This confirms the advantages of the proposed model over the current one, demonstrating the need for a specific case 2 water BRDF correction algorithm as well as the feasibility of enhancing performance of current and future satellite ocean color remote sensing missions for monitoring of typical coastal waters.     

Remote sensing of ocean color: a methodology for dealing with broad spectral bands and significant out-of-band response

A methodology for delineating the influence of finite spectral bandwidths and significant out-of-band response of sensors for remote sensing of ocean color is developed and applied to the Sea-viewing Wide-Field-of-view Sensor (SeaWiFS). The basis of the method is the application of the sensor's spectral-response functions to the individual components of the top-of-the-atmosphere (TOA) radiance rather than the TOA radiance itself. For engineering purposes, this approach allows one to assess easily (and quantitatively) the potential of a particular sensor design for meeting the system-sensor plus algorithms-performance requirements. In the case of the SeaWiFS, two significant conclusions are reached. First, it is found that the out-of-band effects on the water-leaving radiance component of the TOA radiance are of the order of a few percent compared with a sensor with narrow spectral response. This implies that verification that the SeaWiFS system-sensor plus algorithms-meets the goal of providing the water-leaving radiance in the blue in clear ocean water to within 5% will require measurements of the water-leaving radiance over the entire visible spectrum as opposed to just narrow-band (10-20-nm) measurements in the blue. Second, it is found that the atmospheric correction of the SeaWiFS can be degraded by the influence of water-vapor absorption in the shoulders of the atmospheric-correction bands in the near infrared. This absorption causes an apparent spectral variation of the aerosol component between these two bands that will be uncharacteristic of the actual aerosol present, leading to an error in correction. This effect is dependent on the water-vapor content of the atmosphere. At typical water-vapor concentrations the error is larger for aerosols with a weak spectral variation in reflectance than for those that display a strong spectral variation. If the water-vapor content is known, a simple procedure is provided to remove the degradation of the atmospheric correction. Uncertainty in the water-vapor content will limit the accuracy of the SeaWiFS correction algorithm.     

Atmospheric correction over case 2 waters with an iterative fitting algorithm

A modular atmospheric correction algorithm is proposed that uses atmospheric and water contents models to predict the visible and near-infrared reflectances observed by a satellite over water. These predicted values are compared with the satellite reflectances at each pixel, and the model parameters changed iteratively with an error minimization algorithm. The default atmospheric model uses single-scattering theory with a correction for multiple scattering based on lookup tables. With this model we used parameters of the proportions of three tropospheric aerosol types. For the default water content model we need the parameters of the concentrations of chlorophyll, inorganic sediment, and gelbstoff. The diffuse attenuation and backscatter coefficients attributed to these constituents are calculated and used to derive the water-leaving reflectance. Products include water-leaving reflectance, concentrations of water constituents, and aerosol optical depth and type. We demonstrate the application of the method to sea-viewing wide field-of-view sensor by using model data.     

Calibration of SeaWiFS. II. Vicarious techniques

We present an overview of the vicarious calibration of the Sea-Viewing Wide Field-of-view Sensor (SeaWiFS). This program has three components: the calibration of the near-infrared bands so that the atmospheric correction algorithm retrieves the optical properties of maritime aerosols in the open ocean; the calibration of the visible bands against in-water measurements from the Marine Optical Buoy (MOBY); and a calibration-verification program that uses comparisons between SeaWiFS retrievals and globally distributed in situ measurements of water-leaving radiances. This paper describes the procedures as implemented for the third reprocessing of the SeaWiFS global mission data set. The uncertainty in the near-infrared vicarious gain is 0.9%. The uncertainties in the visible-band vicarious gains are 0.3%, corresponding to uncertainties in the water-leaving radiances of approximately 3%. The means of the SeaWiFS/in situ matchup ratios for water-leaving radiances are typically within 5% of unity in Case 1 waters, while chlorophyll a ratios are within 1% of unity. SeaWiFS is the first ocean-color mission to use an extensive and ongoing prelaunch and postlaunch calibration program, and the matchup results demonstrate the benefits of a comprehensive approach.     

Aerosol polarization effects on atmospheric correction and aerosol retrievals in ocean color remote sensing

The current ocean color data processing system for the Sea-viewing Wide Field-of-View Sensor (SeaWiFS) and the moderate resolution imaging spectroradiometer (MODIS) uses the Rayleigh lookup tables that were generated using the vector radiative transfer theory with inclusion of the polarization effects. The polarization effects, however, are not accounted for in the aerosol lookup tables for the ocean color data processing. I describe a study of the aerosol polarization effects on the atmospheric correction and aerosol retrieval algorithms in the ocean color remote sensing. Using an efficient method for the multiple vector radiative transfer computations, aerosol lookup tables that include polarization effects are generated. Simulations have been carried out to evaluate the aerosol polarization effects on the derived ocean color and aerosol products for all possible solar-sensor geometries and the various aerosol optical properties. Furthermore, the new aerosol lookup tables have been implemented in the SeaWiFS data processing system and extensively tested and evaluated with SeaWiFS regional and global measurements. Results show that in open oceans (maritime environment), the aerosol polarization effects on the ocean color and aerosol products are usually negligible, while there are some noticeable effects on the derived products in the coastal regions with nonmaritime aerosols.     

Contribution of water-leaving radiances to multiangle, multispectral polarimetric observations over the open ocean: bio-optical model results for case 1 waters

Multiangle, multispectral photopolarimetry of atmosphere-ocean systems provides the fullest set of remote sensing information possible on the scattering properties of aerosols and on the color of the ocean. Recent studies have shown that inverting such data allows for the potential of separating the retrieval of aerosol properties from ocean color monitoring in the visible part of the spectrum. However, the data in these studies were limited to those principal plane observations where the polarization of water-leaving radiances could be ignored. Examining similar potentials for off-principal plane observations requires the ability to assess realistic variations in both the reflectance for and bidirectionality of polarized water-leaving radiances for such viewing geometries. We provide hydrosol models for use in underwater light scattering computations to study such variations. The model consists of two components whose refractive indices resemble those of detritus-minerallike and planktonlike particles, whose size distributions are constrained by underwater light linear polarization signatures, and whose mixing ratios change as a function of particulate backscattering efficiency. Multiple scattering computations show that these models are capable of reproducing realistic underwater light albedos for wavelengths ranging from 400 to 600 nm, and for chlorophyll a concentrations ranging from 0.03 to 3.0 mg/m(3). Numerical results for spaceborne observations of the reflectance for total and polarized water-leaving radiances are provided as a function of polar angles, and the change in these reflectances with wavelength, chlorophyll a concentration, and hydrosol model are discussed in detail for case 1 (open ocean) waters.     

Optical remote sensing of marine constituents in coastal waters: a feasibility study

Optical remote sensing of ocean color is a well-established technique for inferring ocean properties. However, most retrieval algorithms are based on the assumption that the radiance received by satellite instruments is affected only by the phytoplankton pigment concentration and correlated substances. This assumption works well for open ocean water but becomes questionable for coastal waters. To reduce uncertainties associated with this assumption, we developed a new algorithm for the retrieval of marine constituents in a coastal environment. We assumed that ocean color can be adequately described by a three-component model made up of chlorophyll a, suspended matter, and yellow substance. The simultaneous retrieval of these three marine constituents and of the atmospheric aerosol content was accomplished through an inverse-modeling scheme in which the difference between simulated radiances exiting the atmosphere and radiances measured with a satellite sensor was minimized. Simulated radiances were generated with a comprehensive radiative transfer model that is applicable to the coupled atmosphere-ocean system. The method of simulated annealing was used to minimize the difference between measured and simulated radiances. To evaluate the retrieval algorithm, we used simulated (instead of measured) satellite-received radiances that were generated for specified concentrations of aerosols and marine constituents, and we tested the ability of the algorithm to retrieve assumed concentrations. Our results require experimental validation but show that the retrieval of marine constituents in coastal waters is possible.     

SeaWiFS on-orbit gain and detector calibrations: effect on ocean products

The NASA Ocean Biology Processing Group's Calibration and Validation Team has analyzed the mission-long Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) on-orbit gain and detector calibration time series to verify that lunar calibrations, obtained at nonstandard gains and radiance ranges, are valid for Earth data collected at standard gains and typical ocean, cloud, and land radiances. For gain calibrations, a constant voltage injected into the postdetector electronics allows gain ratios to be computed for all four detectors in each band. The on-orbit lunar gain ratio time series show small drifts for the near infrared bands. These drifts are propagated into the ocean color data through the atmospheric correction parameter epsilon, which uses the 765/865 nm band ratio. An anomaly analysis of global mean normalized water-leaving radiances at 510 nm shows a small decrease over the mission, while an analysis of epsilon shows a corresponding increase. The drifts in the lunar time series for the 765 and 865 nm bands were corrected. An analysis of the revised water-leaving radiances at 510 nm shows the drift has been eliminated, while an analysis of epsilon shows a reduced drift. For detector calibrations, solar diffuser observations made by the individual detectors in each band allows the response of the detectors to be monitored separately. The mission-long time series of detector calibration data show that the variations in the response of the individual detectors are less than 0.5% over the mission for all bands except the 865 nm band, where the variations are less than 1%.     

Role of oceanic air bubbles in atmospheric correction of ocean color imagery

Ocean color is the radiance that emanates from the ocean because of scattering by chlorophyll pigments and particles of organic and inorganic origin. Air bubbles in the ocean also scatter light and thus contribute to the water-leaving radiance. This additional water-leaving radiance that is due to oceanic air bubbles could violate the black pixel assumption at near-infrared wavelengths and be attributed to chlorophyll in the visible. Hence, the accuracy of the atmospheric correction required for the retrieval of ocean color from satellite measurements is impaired. A comprehensive radiative transfer code for the coupled atmosphere--ocean system is employed to assess the effect of oceanic air bubbles on atmospheric correction of ocean color imagery. This effect is found to depend on the wavelength-dependent optical properties of oceanic air bubbles as well as atmospheric aerosols.     

Simultaneous determination of water-leaving reflectance and aerosol optical thickness from Coastal Zone Color Scanner measurements

We present an iterative algorithm for the simultaneous determination of water-leaving reflectances and aerosol optical thickness from the total outgoing radiances measured in the Coastal Zone Color Scanner (CZCS) visible and near-infrared channels. Numerical experiments were carried out to investigate the feasibility of the algorithm. The results show that the errors in determined water-leaving reflectance at 0.443 mum are approximately 10% in most cases. The errors in determined water-leaving reflectance at 0.55 mum are in the 4.05-7.2% range. The errors in the simultaneously determined aerosol optical thickness for all CZCS visible and near-infrared channels are less than approximately 10%.     

Atmospheric correction of SeaWiFS imagery for turbid coastal and inland waters

The standard SeaWiFS atmospheric correction algorithm, designed for open ocean water, has been extended for use over turbid coastal and inland waters. Failure of the standard algorithm over turbid waters can be attributed to invalid assumptions of zero water-leaving radiance for the near-infrared bands at 765 and 865 nm. In the present study these assumptions are replaced by the assumptions of spatial homogeneity of the 765:865-nm ratios for aerosol reflectance and for water-leaving reflectance. These two ratios are imposed as calibration parameters after inspection of the Rayleigh-corrected reflectance scatterplot. The performance of the new algorithm is demonstrated for imagery of Belgian coastal waters and yields physically realistic water-leaving radiance spectra. A preliminary comparison with in situ radiance spectra for the Dutch Lake Markermeer shows significant improvement over the standard atmospheric correction algorithm. An analysis is made of the sensitivity of results to the choice of calibration parameters, and perspectives for application of the method to other sensors are briefly discussed.     

Correction of Sun glint Contamination on the SeaWiFS Ocean and Atmosphere Products

For remote sensing of the ocean and atmosphere optical properties, the measurement of radiances affected by sun glint has to be avoided and/or masked out. There are usually no meaningful retrievals in regions significantly contaminated by sun glint. The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) is capable of operationally tilting the sensor 20 degrees away from nadir to minimize sun glint contamination. The sun glint mask is computed from the Cox and Munk model [J. Opt. Soc. Am. 44, 838-850 (1954)] and applied to the SeaWiFS data. However, sun glint is still a factor near the subsolar point. We present results that demonstrate the effect of sun glint contamination on retrievals of ocean bio-optical and atmospheric products. We show that, although sun glint contamination has a minor effect on retrieved ocean color products, the effect on retrieved atmospheric products, e.g., aerosol optical thickness, is significant. We describe a sun glint correction scheme implemented in the SeaWiFS data-processing system and compare the results with and without sun glint correction. With sun glint correction the derived ocean and atmospheric products are improved. Also, the sun glint masked area can be reduced and therefore can increase significantly the coverage area near the subsolar point.     

Moderate-resolution imaging spectroradiometer ocean color polarization correction

The polarization correction for the Moderate-Resolution Imaging Spectroradiometer (MODIS) instruments on the Terra and Aqua satellites is described. The focus is on the prelaunch polarization characterization and on the derivation of polarization correction coefficients for the processing of ocean color data. The effect of the polarization correction is demonstrated. The radiances at the top of the atmosphere need to be corrected by as much as 3.2% in the 412 nm band. The effect on the water-leaving radiances can exceed 50%. The polarization correction produces good agreement of the MODIS Aqua water-leaving radiance time series with data from another, independent satellite-based ocean color sensor, the Sea-Viewing Wide Field-of-View Sensor (SeaWiFS).     

Atmospheric correction of ocean color imagery: use of the junge power-law aerosol size distribution with variable refractive index to handle aerosol absorption

When strongly absorbing aerosols are present in the atmosphere, the usual two-step procedure of processing ocean color data-(1) atmospheric correction to provide the water-leaving reflectance (rho(w)), followed by (2) relating rho(w) to the water constituents-fails and simultaneous estimation of the ocean and aerosol optical properties is necessary. We explore the efficacy of using a simple model of the aerosol-a Junge power-law size distribution consisting of homogeneous spheres with arbitrary refractive index-in a nonlinear optimization procedure for estimating the relevant oceanic and atmospheric parameters for case 1 waters. Using simulated test data generated from more realistic aerosol size distributions (sums of log-normally distributed components with different compositions), we show that the ocean's pigment concentration (C) can be retrieved with good accuracy in the presence of weakly or strongly absorbing aerosols. However, because of significant differences in the scattering phase functions for the test and power-law distributions, large error is possible in the estimate of the aerosol optical thickness. The positive result for C suggests that the detailed shape of the aerosol-scattering phase function is not relevant to the atmospheric correction of ocean color sensors. The relevant parameters are the aerosol single-scattering albedo and the spectral variation of the aerosol optical depth. We argue that the assumption of aerosol sphericity should not restrict the validity of the algorithm and suggest an avenue for including colored aerosols, e.g., wind-blown dust, in the procedure. A significant advantage of the new approach is that realistic multicomponent aerosol models are not required for the retrieval of C.     

Atmospheric Correction of Ocean Color Imagery: Test of the Spectral Optimization Algorithm with the Sea-Viewing Wide Field-of-View Sensor

We implemented the spectral optimization algorithm [SOA; Appl. Opt. 37, 5560 (1998)] in an image-processing environment and tested it with Sea-viewing Wide Field-of-View Sensor (SeaWiFS) imagery from the Middle Atlantic Bight and the Sargasso Sea. We compared the SOA and the standard SeaWiFS algorithm on two days that had significantly different atmospheric turbidities but, because of the location and time of the year, nearly the same water properties. The SOA-derived pigment concentration showed excellent continuity over the two days, with the relative difference in pigments exceeding 10% only in regions that are characteristic of high advection. The continuity in the derived water-leaving radiances at 443 and 555 nm was also within ~10%. There was no obvious correlation between the relative differences in pigments and the aerosol concentration. In contrast, standard processing showed poor continuity in derived pigments over the two days, with the relative differences correlating strongly with atmospheric turbidity. SOA-derived atmospheric parameters suggested that the retrieved ocean and atmospheric reflectances were decoupled on the more turbid day. On the clearer day, for which the aerosol concentration was so low that relatively large changes in aerosol properties resulted in only small changes in aerosol reflectance, water patterns were evident in the aerosol properties. This result implies that SOA-derived atmospheric parameters cannot be accurate in extremely clear atmospheres.     

Columnar aerosol properties over oceans by combining surface and aircraft measurements: simulations

We report an algorithm that can be used to invert the radiance exiting the top and bottom of the atmosphere to yield the columnar optical properties of atmospheric aerosol under clear sky conditions over the oceans. The method is an augmentation of a similar algorithm presented by Wang and Gordon [Appl. Opt. 32, 4598 (1993)] that used only sky radiance, and therefore was incapable of retrieving the aerosol phase function at the large scattering angles that are of critical importance in remote sensing of oceanic and atmospheric properties with satellites. Well-known aerosol models were combined with radiative transfer theory to simulate pseudodata for testing of the algorithm. The tests suggest that it should be possible to retrieve the aerosol phase function and the aerosol single-scattering albedo accurately over the visible spectrum at aerosol optical thicknesses as large as 2.0. The algorithm is capable of retrievals with such large optical thicknesses because all significant orders of multiple scattering are included. We believe that combining an algorithm of this type with surface-based and high-altitude aircraft-based radiance measurements could be useful for studying aerosol columnar optical properties over oceans and large lakes. The use of the retrieval method is possible over the ocean because, unlike the land surface, the albedo of the ocean is low and spatially uniform.     

Atmospheric correction over case 2 waters with an iterative fitting algorithm: relative humidity effects

In algorithms for the atmospheric correction of visible and near-IR satellite observations of the Earth's surface, it is generally assumed that the spectral variation of aerosol optical depth is characterized by an Angstr?m power law or similar dependence. In an iterative fitting algorithm for atmospheric correction of ocean color imagery over case 2 waters, this assumption leads to an inability to retrieve the aerosol type and to the attribution to aerosol spectral variations of spectral effects actually caused by the water contents. An improvement to this algorithm is described in which the spectral variation of optical depth is calculated as a function of aerosol type and relative humidity, and an attempt is made to retrieve the relative humidity in addition to aerosol type. The aerosol is treated as a mixture of aerosol components (e.g., soot), rather than of aerosol types (e.g., urban). We demonstrate the improvement over the previous method by using simulated case 1 and case 2 sea-viewing wide field-of-view sensor data, although the retrieval of relative humidity was not successful.     

Atmospheric correction of satellite ocean color imagery: the black pixel assumption

The assumption that values of water-leaving radiance in the near-infrared (NIR) are negligible enable aerosol radiative properties to be easily determined in the correction of satellite ocean color imagery. This is referred to as the black pixel assumption. We examine the implications of the black pixel assumption using a simple bio-optical model for the NIR water-leaving reflectance [rho(w)(lambda(NIR))](N). In productive waters [chlorophyll (Chl) concentration >2 mg m(-3)], estimates of [rho(w)(lambda(NIR))](N) are several orders of magnitude larger than those expected for pure seawater. These large values of [rho(w)(lambda(NIR))](N) result in an overcorrection of atmospheric effects for retrievals of water-leaving reflectance that are most pronounced in the violet and blue spectral region. The overcorrection increases dramatically with Chl, reducing the true water-leaving radiance by roughly 75% when Chl is equal to 5 mg m(-3). Relaxing the black pixel assumption in the correction of Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) satellite ocean color imagery provides significant improvements in Chl and water-leaving reflectance retrievals when Chl values are greater than 2 mg m(-3). Improvements in the present modeling of [rho(w)(lambda(NIR))](N) are considered, particularly for turbid coastal waters. However, this research shows that the effects of nonzero NIR reflectance must be included in the correction of satellite ocean color imagery.