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

Bryan A. Franz
SAIC, NASA SIMBIOS Project
June 2003

Abstract

Introduction

Data Sources



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

Table 1: Band Correspondence (nm)



Subset Definitions

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 (OC4v4 algorithm, O'Reilly et al., 2000) and MODIS (OC3M algorithm, O'Reilly et al., 2000) 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. Some caution should be exercized when comparing the clear-water subsetted data, as 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, six basin-scale subsets were analyzed. These included regions in the northern Pacific (PacN), north western Pacific (PacNW), south eastern Pacific (PacSE), northern Atlantic (AtlN), southern Atlantic (AtlS), and the southern Indian Ocean (IndS). To these, a smaller region near Hawaii was also added. All of these subset regions were adopted from Fougnie et al. 2002. Their locations are listed in Table 2. Based on the results of the regional analysis, yet another group of 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 3.

Region
ID
Minimum
Latitude
Maximum
Latitude
Minimum
Longitude
Maximum
Longitude
Hawaii18.019.9-158.5-156.5
PacN15.023.0-180.0-159.4
PacNW10.022.7139.5165.6
PacSE-44.9-20.7-130.2-89.0
AtlN17.027.0-62.5-44.2
AtlS-19.9-9.9-32.3-11.0
IndS-29.9-21.289.5100.1

Table 2: Regional Subset Definitions



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, where the SeaWiFS chlorophylls correspond to the standard chlor_a product (OC4v4 algorithm) and the MODIS chlorophyll corresponds with the standard chlor_a_2 product (OC3M algorithm). The solid vertical lines in the temporal trend plots are provided as a reference to indicate the transitions between MODIS Oceans calibration epochs. These epochs are the independent periods over which MODIS calibration corrections were derived and implemented by the MODIS Oceans group at the University of Miami (RSMAS). In most cases, these periods correspond with the calibration epochs used by the MODIS Calibration Support Team (MCST) for the adjustment of the Level-1B radiances, and they usually correspond with spacecraft safe-hold events or significant instrument state changes.

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
  • Regional Water-Leaving Radiance Mission Trend Plots
  • Regional 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
  • PacN Chlorophyll Images
  • PacNW Chlorophyll Images
  • PacSE Chlorophyll Images
  • AtlN Chlorophyll Images
  • AtlS Chlorophyll Images
  • IndS Chlorophyll Images
  • Hawaii 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
  • Full Product Images

  • Complete Chlor-a Images
  • Complete nLw(412) Images
  • Complete nLw(443) Images
  • Complete nLw(490) Images
  • Complete nLw(555) 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.073 0.0018 2.213 0.0490 1.862 0.0309 1.235 0.0101 0.287 0.0019
    MODIS Clear 0.087 0.0036 2.397 0.0816 2.110 0.0262 1.480 0.0284 0.397 0.0136
    SeaWiFS Deep 0.179 0.0138 1.736 0.0324 1.530 0.0208 1.123 0.0128 0.324 0.0082
    MODIS Deep 0.176 0.0097 1.890 0.0553 1.741 0.0249 1.343 0.0244 0.430 0.0130
    SeaWiFS Coastal 0.947 0.2891 0.820 0.0756 0.881 0.0649 0.861 0.0556 0.411 0.0280
    MODIS Coastal 0.665 0.0710 0.943 0.0514 1.036 0.0447 1.042 0.0271 0.510 0.0186

    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).

    Fougnie, B., P. Henry, A. Morel, D. Antoine, and F. Montagner, 2002: Identification and Characterization of Stable Homogeneous Oceanic Zones: Climatology and Impact on In-Flight Calibration of Space Sensors over Rayleigh Scattering. Ocean Optics XVI, Santa Fe, NM, November 18-22, 2002.

    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