A Comparative Time-Series Analysis of Ocean Color Products from MODIS/Aqua and SeaWiFS

Bryan Franz (SAIC)
Ocean Discipline Processing Group
2 July 2004

Introduction

This document presents a comparison of MODIS/Aqua Level-3 products with equivalent SeaWiFS products. Both the SeaWiFS and MODIS/Aqua products are those currently processed and distributed by the Ocean Discipline Processing Group at NASA/GSFC (or the project formerly known as SeaWiFS/SIMBIOS). The analysis looks at average values of coincident retrievals on global and regional spatial scales, and presents the results as a comparative time-series over the common mission lifespan.

Description of Analysis

The purpose of this analysis is to produce a statistically rigorous comparison of equivalent MODIS and SeaWiFS ocean color products, to provide insight to the users of these data sets as to their relative agreement. The analysis also serves to provide feedback to the calibration and algorithm development process as to the relative impact of calibration and algorithm changes on global spatial scales and life-of-mission time scales. The ocean color products compared are the standard chlorophyll products derived from each mission data set, as well as the normalized water-leaving radiances in the four closest visible wavelengths. The equivalent wavelengths are listed in Table 1. The chlorophyll products compared are the standard "chlor_a" products produced for each sensor, which means the OC3M algorithm for MODIS and the OC4v4 algorithm for SeaWiFS. These are both emperical max-band-ratio algorithms developed by O'Reilly (O'Reilly et al. 2000).

Band SeaWiFS MODIS
1 412 412
2 443 443
3 490 488
5 555 551

Table 1: Band Correspondence (nm)



The analysis began with MODIS and SeaWiFS Level-3 equal-area binned products, which were composited over 8-day periods. The standard SeaWiFS products are 9-km resolution, but the standard MODIS products are 4.6-km resolution. To allow for a direct, bin-for-bin comparison, the MODIS products were rebinned to the SeaWiFS 9-km resolution using standard binning algorithms.

With 8-day composited SeaWiFS and MODIS data products in an equivalent form, the data sets were further divided into several geographic subsets. Three global subsets were defined, corresponding to clear water, deep water, and coastal water. The deep water subset consists of all bins where water depth is greater than 1000 meters. Clear water was defined as deep water where the retrieved chlorophyll is less than 0.15 mg/m^3. For the chlorophyll test, both SeaWiFS and MODIS retrievals were required to be below the clear-water threshold. Coastal water was defined as all bins where water depth is between 30 and 1000 meters, as defined by a shallow water mask and the deep water mask.

When comparing the clear-water subsetted data, it should be noted that anomalously high chlorophyll retrievals from either sensor can significantly alter the geographic distribution of selected bins. In contrast, the deep-water and coastal subsets are purely geographic in selection criteria. The coastal subset, however, is more likely to contain regions of significant variability in water structure and atmospheric conditions, as well as case 2 water types, all of which can lead to greater retrieval uncertainty and larger differences between the two sensors. The deep-water subset is, therefore, the most stable subset for cross-sensor comparison of retrieved oceanic optical properties. The geographic extent of all three global subsets will vary, however, with the seasonal change in earth illumination and thus sensor imaging duty cycle. Figure 1 presents a sample pair of MODIS and SeaWiFS deep water subsetted chlorophyll images for one 8-day period in May of 2003. The images show the geographic extent of the common-binned, deep-water subset, and they provide some insight into the qualitative agreement between the two sensors.

SeaWiFS

MODIS

Log10(Chla), 0.01 - 1.0 mg/m^3

Figure 1: Sample Chlorophyll Images, Deep-Water Subset
Days 137-144 of 2003



In addition to the global subsets, a set of zonal subsets was defined to provide a systematic means for investigating latitudinally-dependent differences between the two sensors. A longitudinal segment of the Pacific from 170W to 150W was divided into 10-deg latitude zones. These zonal subsets are summarized in Table 2.

Region
ID
Minimum
Latitude
Maximum
Latitude
Minimum
Longitude
Maximum
Longitude
PacN5040.050.0-170.0-150.0
PacN4030.040.0-170.0-150.0
PacN3020.030.0-170.0-150.0
PacN2010.020.0-170.0-150.0
PacN100.010.0-170.0-150.0
PacS10-10.00.0-170.0-150.0
PacS20-20.0-10.0-170.0-150.0
PacS30-30.0-20.0-170.0-150.0
PacS40-40.0-30.0-170.0-150.0
PacS50-50.0-40.0-170.0-150.0

Table 3: Zonal Subset Definitions



Trending Analysis

For each sensor, for each 8-day product, the filled bins associated with a particular subset were identified and used to compute the mean, standard deviation, and average observation time. Figure 2 shows an example of a typical trend plot derived from this analysis. For the plot on the left, the common MODIS and SeaWiFS bins for the deep-water subset were spatially averaged for each 8-day-binned water-leaving radiance product, and the resulting means were then plotted as a function of time. MODIS is shown as dashed lines. The colors indicate different bands, as summarized in Table 1. The plot on the right shows the same data as a ratio, with MODIS means normalized by SeaWiFS means. Similarly, Figure 3 shows the temporal trends in mean chlorophyll.

Figure 2: SeaWiFS and MODIS Water-Leaving Radiance Comparison, Deep-Water Subset



Figure 3: SeaWiFS and MODIS Chlorophyll Comparison, Deep-Water Subset



Links to additional trend plots are provided below, along with tabulated data and global images. The subset images show mapped chlorophyll for each 8-day period, for the common bins associated with each geographic subset. These images allow for a qualitative assessment of the agreement between the two sensors, and indicate the spatial extent of the subsetted, common bins. Full product images of the chlorophyll and radiance data are also provided, allowing comparison between the two sensors prior to subsetting or reduction to common bins. Finally, tabulated results of the mean and standard deviation for each product, for each 8-day subset are provided. The tabulated means are the values plotted in the trend plots.

Time-Series Trend Plots

  • Global Water-Leaving Radiance Mission Trend Plots
  • Global Chlorophyll Mission Trend Plots
  • Zonal Water-Leaving Radiance Mission Trend Plots
  • Zonal Chlorophyll Mission Trend Plots
  • Subset Images

  • Deep Water Chlorophyll Images
  • Clear Water Chlorophyll Images
  • Coastal Water Chlorophyll Images
  • PacN50 Chlorophyll Images
  • PacN40 Chlorophyll Images
  • PacN30 Chlorophyll Images
  • PacN20 Chlorophyll Images
  • PacN10 Chlorophyll Images
  • PacS10 Chlorophyll Images
  • PacS20 Chlorophyll Images
  • PacS30 Chlorophyll Images
  • PacS40 Chlorophyll Images
  • PacS50 Chlorophyll Images
  • Tabulated Results

  • Deep Water Statistics
  • Clear Water Statistics
  • Coastal Water Statistics
  • Discussion of Results

    Trend Statistics

    Chlor_aChlor_a nLw_412nLw_412 nLw_443nLw_443 nLw_490nLw_490 nLw_555nLw_555
    Sensor Subset meanstdev meanstdev meanstdev meanstdev meanstdev
    SeaWiFS Clear 0.077 0.0035 2.197 0.0719 1.856 0.0493 1.249 0.0179 0.296 0.0042
    MODIS Clear 0.077 0.0036 2.152 0.0926 1.833 0.0729 1.282 0.0377 0.315 0.0123
    SeaWiFS Deep 0.183 0.0125 1.736 0.0529 1.530 0.0372 1.134 0.0186 0.334 0.0068
    MODIS Deep 0.189 0.0154 1.689 0.0850 1.496 0.0599 1.148 0.0371 0.344 0.0159
    SeaWiFS Coastal 1.115 0.2043 0.778 0.0580 0.862 0.0475 0.888 0.0420 0.474 0.0285
    MODIS Coastal 1.340 0.2352 0.727 0.0837 0.848 0.0561 0.916 0.0465 0.501 0.0299

    Table 4: Global Trend Statistics



    Summary

    References

    Clark, D. K., M. E. Feinholz, M. A. Yarbrough, B. C. Johnson, S. W. Brown, Y. S. Kim, and R. A. Barnes, Overview of the radiometric calibration of MOBY, Proc. Spie, 4483, 64-76 (2001).

    Eplee, R.E., Jr., R.A. Barnes, and F.S. Patt, 2003a: "Changes to the on-orbit calibration of SeaWiFS." In: Patt, F.S., R.A. Barnes, R.E. Eplee, Jr., B.A. Franz, W.D. Robinson, G.C. Feldman, S.W. Bailey, P.J. Werdell, R. Frouin, R.P. Stumpf, R.A. Arnone, R.W. Gould, Jr., P.M. Martinolich, and V. Ransibrahmanakul, Algorithm Updates for the Fourth SeaWiFS Data Reprocessing, NASA Tech. Memo. 2003--206892, Vol. 22, S.B. Hooker and E.R. Firestone, Eds., NASA Goddard Space Flight Center, Greenbelt, Maryland, (in press).

    Eplee, R.E., Jr., R.A. Barnes, S.W. Bailey, and P.J. Werdell, 2003b: "Changes to the vicarious calibration of SeaWiFS." In: Patt, F.S., R.A. Barnes, R.E. Eplee, Jr., B.A. Franz, W.D. Robinson, G.C. Feldman, S.W. Bailey, P.J. Werdell, R. Frouin, R.P. Stumpf, R.A. Arnone, R.W. Gould, Jr., P.M. Martinolich, and V. Ransibrahmanakul, Algorithm Updates for the Fourth SeaWiFS Data Reprocessing, NASA Tech. Memo. 2003--206892, Vol. 22, S.B. Hooker and E.R. Firestone, Eds., NASA Goddard Space Flight Center, Greenbelt, Maryland, (in press).

    O'Reilly, J., S. Maritorena, M. O'Brien, D. Siegel, D. Toole, D. Menzies, R. Smith, J. Mueller, B. Mitchell, M. Kahru, F. CHavez, P. Strutton, G. Cota, S. Hooker, C. McClain, K. Carder, F. Muller-Karger, L. Harding, A. Magnuson, D. Phinney, G. Moore, J. Aiken, K. Arrigo, R. Letelier and M. Culver (2000). SeaWiFS postlaunch technical report series, Volume 11, SeaWiFS postlaunch calibration and validation analyses, Part 3, NASA Technical Memorandum.


    Bryan A. Franz
    Last modified: Thu Sep 12 14:33:29 EDT 2002