We are currently reprocessing the AMSR2 data.  When we begin processing again, we will deliver a new version, version 7.2

OISST Product Description

Introduction

The through-cloud capabilities of microwave radiometers provide a valuable picture of global sea surface temperature (SST). To utilize this, scientists at Remote Sensing Systems have calculated a daily, Optimally Interpolated (OI) SST product at 25 km and 9 km resolutions. This product is ideal for research activities in which a complete, daily SST map is more desirable than one with missing data due to orbital gaps or environmental conditions precluding SST retrieval.  The 9 km MW+IR OI SST product combines the through-cloud capabilities of the MW data with the high spatical resolution of the IR SST data. 

Improved global daily NRT SSTs are useful for a wide range of scientific and operational activities.

MW OI SST Products

Several optimally interpolated (OI) SST products are created from the microwave and infrared SSTs.

OI SST

Observations

Coverage

Time Span

tmi

TMI

40N to 40S

1998-Jan to present

amsre

AMSR-E

global

2002-Jun to Oct 4th

tmi_amsre

TMI, AMSR-E, WindSAT

global

2002-Jun to present

mw_ir

TMI, AMSR-E, WindSAT, Terra MODIS, Aqua MODIS

global

2006-Jan to present

The tmi_amsre and mw_ir SST data both use WindSAT data from 10/2011 onward, due to the failure of AMSR-E. We are currently working in reprocess all the datasets. In the future there will simply be a microwave OI SST at 25 km and a microwave+infrared SST at 9 km.

One analysis contains only TMI SSTs, another contains only AMSR-E SSTs, and the final analysis blends the two SSTs. The OI TMI SSTs are available (latitudes 40S-40N) from January 1998 to the present, OI AMSR-E SSTs are available globally from June 2002 to the present, and the combined global OI TMI+AMSR-E SSTs are available from June 2002 to the present. All products are updated several times daily in NRT at 25 km resolution. These fields are intended as research for the Multi-sensor Improved SST (MISST) project, which is a US contribution to the Global Ocean Data Assimilation Experiment (GODAE) High-Resolution SST Pilot Project (GHRSST-PP).

Sensor-specific Error Calculation

The following corrections and analyses of errors are necessary first steps towards producing the MW OI SST products. Each step is further described below.

Correcting for TMI’s Emissive Antenna

The antenna coating of the TMI sensor was oxidized in orbit soon after launch, causing errors in the TMI observations. A correction was developed (Wentz, 2001), but proved to be incomplete in removing the error. A bias still exists in TMI data, which is a function of local observation time (Gentemann, accepted JGR). To account for this, an additional correction is applied before TMI data are included in the OI analysis.

Estimation and Removal of Diurnal Warming

Before blending the satellite data, we consider the data sampling of each instrument. For example, the sun synchronous orbit of MODIS and AMSR-E on Aqua yields retrievals at a local time of approximately 1:30 AM and 1:30 PM. During the daytime over-pass, solar heating of the ocean surface can cause warming of up to 3° C (Price et al, 1986, Yokohama, 1996). Currently, many OI SST algorithms either ignore daytime retrievals or assign them a higher error than nighttime retrievals. While simply removing the daytime retrievals from the objective analysis does prevent warm retrievals from 'contaminating' the final product, the number of samples can be reduced by half. In well-sampled regions this may not impact the final product, but the infrared SSTs used in most analyses have large regions where few retrievals exist each month due to persistent cloud cover, making the daytime retrievals extremely valuable. Assigning the daytime retrievals a higher error (and therefore a smaller weight in the objective analysis) reduces diurnal 'contamination' of the dataset, but at the risk of still including some component of diurnal warming. The AMSR-E and TMI OI SSTs include day and night observations. To optimally utilize MW daytime retrievals, a simple empirical model of diurnal warming was developed that depends on solar insolation, wind speed, and local time of observation (Gentemann, 2003). Solar insolation is calculated as a function of latitude and day of year; wind speed is simultaneously retrieved with SST from radiometer observations. Using this diurnal model, all MW SSTs are 'normalized' to a foundation SST.  For more information, see the GHRSST definition.  

Sensor Errors for OI Analysis

Microwave SST retrieval errors are mainly a function of wind speed and SST. These errors are added in a root-sum-squared sense to the daily standard deviation (STD) derived from buoy collocations to obtain a total retrieval error.

Additional Quality Control

Some rain contaminated SSTs exist in the microwave data. At the edges of rain cells, there is often undetected rain that causes a biased SST retrieval. Two tests attempt to remove rain contaminated SSTs. First, at each SST retrieval the STD is calculated using all data within one day and 100 km of the cell. SSTs falling outside of 3 STDs are flagged and removed from the data set. This process is then further repeated to remove outliers. Next, the SST is compared to the previous day's OI SST value. Any SSTs within 100 km of a rain pixel that are more than 0.6 C warmer than the previous day's OI SST value are removed.

Some cloud contaminated SSTs exist in the infrared data. At the edges of cloud cells, there is often undetected cloud that causes a biased SST retrieval. We use a similar test as described above, to remove spurious cloud contaminated retreivals from the infrared SST.

Optimum Interpolation (OI)

After characterizing the errors listed above, the SSTs are blended together using the OI scheme described in Reynolds and Smith (1994). OI is a widely utilized method in oceanography and meteorology that makes use of the statistical properties of irregularly spaced data (in time and space) to interpolate the data onto a regularly sampled grid. For each dataset included in the analysis, error characteristics must be understood or at least estimated.

A first-guess field, the previous day's OI SST, is employed to calculate data increments, which are all nearby data minus the first-guess field. The new SST estimate is formed by a weighted sum of increments, with the weights calculated by the OI method, added to the first guess data. Correlation scales of 4 days and 100 km are used in determining the weights used in our methodology.

Known Problems

Undetected Sea Ice

Undetected sea ice is causing some unrealistically warm SST values to appear in these products. The problem is most apparent near ice edges, especially as the ice edge advances or retreats.

The first set of images (below) illustrates the problem occuring in the Beaufort Sea, Artic Ocean, over a seven day period. In the images on the left, the ~4° C (light blue, circled) SSTs are probably artifacts of a thin layer of sea ice or slush. As the sea ice solidifies, it becomes more accurately identified as the images progress in time towards the right.

undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice
color bar
undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice undetected sea ice

The second set of images (above) tracks retreating Antarctic sea ice over ten days. Here we see probable ice causing up to ~5° C warm artifacts.

Missing Data

Satellite instruments are occasionally unavailable. Near real time OI SST products will be created for the current day, even if no new observations exist. The OI method utilizes a first guess field, which in this analysis is the previous day's OI SST. If there are no new observations, the new SST estimate is the previous day's SST.

Instrument observations are missing for the following dates:

Instrument   

Missing Data

# Days

TMI

1999.01.04 - 1999.01.05 
2000.09.18 
2001.08.14 - 2001.08.16



3

AMSR-E

2002.06.28 
2002.07.30 - 2002.08.07 
2002.09.13 - 2002.09.19 
2003.10.30 - 2003.11.05 
2004.11.19 
2006.11.18 
2007.11.27 - 2007.11.28 
2010.02.02 - 2010.02.05








3

For example, AMSR-E was unavailable September 13-19, 2002. For these dates, the combined TMI + AMSR-E product accurately represents detailed daily SSTs in TMI range (±40°), but at latitudes greater than 40° the OI SST values change little because no AMSR-E observations exist. As more satellites are added to this analysis, the chance of this 'frozen' data will diminish. All OI analyses suffer from this problem.

Daily browse imagery for the TMI and AMSR-E instrument products can show the observations available on any given day.

Gridded Binary Data File Format

Each binary data file available at ftp.discover-earth.org/sst consists of three 0.25 x 0.25 degree grids (1440 x 720 array) of single byte values representing a given day's SSTs, interpolation ERROR estimates, and data MASKing information. The MASK has bit values: leftmost bit (bit 0) = land if set to 1, bit 1 = ice, bit 2 shows whether IR data were used, bit 3 shows whether MW data were used, and bit 4 is set to 1 for bad data. Interim products ("rt") are updated several times per day until the data become final ("v03").

File names follow these conventions:

MW OI SST

directory path

file name

TMI daily/tmi/ tmi.fusion.yyyy.doy.ver.gz
AMSR-E daily/amsre/ amsre.fusion.yyyy.doy.ver.gz
TMI + AMSR-E daily/tmi_amsre/ tmi_amsre.fusion.yyyy.doy.ver.gz

Where "yyyy", "doy", and "ver" stand for:

yyyy year "2002", "2003", etc.
doy day of year "001" (Jan-1), "002" (Jan-2), etc.
ver version
"rt" = near real time (interim product)
"v03"  = version 3 (final product)

The center of the first cell of the 1440 column and 720 row map is at 0.125 E longitude and -89.875 latitude. The center of the second cell is 0.375 E longitude, -89.875 latitude.

Byte values range from 0 to 255.   The read routines supplied insert the following specific values into the SST and ERROR data grids as below:

0 to 250 = valid SST data
251 = missing data
252 = sea ice
253, 254 = missing data
255 = land mass

Byte values 0 - 250 need to be scaled to obtain standard units:

SST: (byte value * 0.15) - 3.0 yields temperature between -3.0 and 34.5 °C

Thus, to convert SST byte values (0 - 250, inclusive) to degrees Centigrade, multiply by .15, then subtract 3.

All binary data files have gzip compression to reduce size and decrease transfer time.

Read Routines:

Read routines are available in IDL, Matlab and Fortran at:ftp.discover-earth.org/sst/daily/read_routines. These read routines (dated August 2010 or later) read the TMI, AMSRE and TMI-AMSRE v03 OISST files.

Acknowledgements

Microwave OI SST data are produced by Remote Sensing Systems and sponsored by National Oceanographic Partnership Program (NOPP), the NASA Earth Science Physical Oceanography Program, and the NASA MEaSUREs DISCOVER Project. Data are available at www.remss.com.

Research into SST blending, diurnal warming, observation errors, and near real-time validation of TMI and AMSR-E SST is supported by the NASA ESE Physical Oceanography Program (Dr. Eric Lindstrom) and the NASA ESE AMSR-E Science Team (Dr. Ramesh Kakar).

The distribution, web-interface, and visualization tools for these data sets are supported by the NASA ESE MEaSUREs Project (Dr. Lucia Tsaoussi).