muri 24 vandelst

JCSDA Infrared Sea Surface Emissivity Model Status: 

JCSDA Infrared Sea Surface Emissivity Model Status Paul van Delst 2nd MURI Workshop 27-28 April 2004 Madison WI

Introduction: 

Introduction Global Data Assimilation System (GDAS) at NCEP/EMC previously used IRSSE model based on Masuda. Doesn’t include effect of enhanced emission due to reflection from sea surface. Only an issue for larger view angles. Coarse frequency resolution. Upgraded the model Use Wu-Smith methodology to compute sea surface emissivity spectra. Reflectivity is average of horizontal and vertical components. Assume that IR sensors are not sensitive to the different polarisations. Refractive index data used: Hale & Querry for real part (pure water) Segelstein for imaginary part (pure water) Friedman for salinity/chlorinity correction Instrument SRFs used to produce sensor channel emissivities. These are the predicted quantities.

IRSSE Model (1): 

IRSSE Model (1) Started with model used in ISEM-6 (Sherlock,1999). where and N1, N2 are integers. The coefficients c0, c1, and c2 for a set of N1 and N2 are determined by regression with a maximum residual cutoff of =0.0002. Only wind speeds of 0.0ms-1 were fit in ISEM-6. The variation of emissivity with wind speed (for HIRS Ch8) was found to be much more than 0.0002.

Wind Speed Dependence of Emissivity: 

Wind Speed Dependence of Emissivity Larger 

IRSSE Model (2): 

IRSSE Model (2) Since the variation with wind speed was greater than 0.0002, the exponents, N1 and N2, of the emissivity model were also allowed to vary. For integral values of N1 and N2 their variation with wind speed suggested inverse relationships for both. The exponents were changed to floating point values, and the fitting exercise was repeated. The result shows a smooth relationship.

Wind Speed Dependence of Integral Exponents: 

Wind Speed Dependence of Integral Exponents

Wind Speed Dependence of Real Exponents: 

Wind Speed Dependence of Real Exponents

IRSSE Model (3): 

IRSSE Model (3) The model was slightly changed to, where v is the wind speed in ms-1. Generating the coefficients For a series of wind speeds, the coefficients ci were obtained. Interpolating coefficients for each ci as a function of wind speed were determined. These are stored in the model datafiles. Using the model For a given wind speed, the ci are computed. These coefficients are then used to compute the view angle dependent emissivity

Emissivity Coefficient Variation By Channel for NOAA-17 HIRS/3: 

Emissivity Coefficient Variation By Channel for NOAA-17 HIRS/3

Emissivity Coefficient Variation By Channel for AIRS M8 (~850-900cm-1): 

Emissivity Coefficient Variation By Channel for AIRS M8 (~850-900cm-1)

TOA TB Residuals for NOAA-17 HIRS.RMS for all wind speeds: 

TOA TB Residuals for NOAA-17 HIRS. RMS for all wind speeds

TOA TB Residuals for AIRS 281 subset.RMS for all wind speeds: 

TOA TB Residuals for AIRS 281 subset. RMS for all wind speeds

TOA TB Residuals for NOAA-17 HIRS.RMS for all wind speeds; only 0ms-1  predicted: 

TOA TB Residuals for NOAA-17 HIRS. RMS for all wind speeds; only 0ms-1  predicted

TOA TB Residuals for AIRS 281 subset.RMS for all wind speeds; only 0ms-1  predicted: 

TOA TB Residuals for AIRS 281 subset. RMS for all wind speeds; only 0ms-1  predicted

TOA TB Residuals: 

TOA TB Residuals When wind speed is taken into account: Residuals are relatively independent of view angle and channel. Magnitudes (Ave., RMS, and Max) are ~10-4–10-3K. When only 0.0ms-1 emissivities are predicted: Residuals peak for largest view angles. Shortwave channels appear to be more sensitive. Magnitudes can be > 0.1K for high view angles. For angles < 40-45, residuals are typically <0.02K

Code Availability: 

Code Availability Three parts of the code Code to compute spectral emissivities (Fortran90) and refractive index netCDF datafiles Code to fit model and produce coefficients (IDL) IRSSE model code (Fortran90) and coefficient datafiles. (Operational code used in the GDAS.) IRSSE model code and datafiles available at http://cimss.ssec.wisc.edu/~paulv Follow the “Infrared Sea Surface Emissivity (IRSSE) Model” link.

Code Availability: 

Code Availability

Issues: 

Issues Use of Cox-Munk probability distribution function (PDF) for slopes of wind driven waves. Experimental data obtained for slopes <0.36. Extrapolations for larger slopes. PDF can have (unphysical) negative probabilities for these larger slopes. Ebuchi and Kizu (2002) PDF derived slope statistics may be more applicable to satellite-based remote sensing. Much larger data sample using GMS-5 visible images and NSCAT, ERS-1, and ERS-2 scatterometer data products. Narrower PDF and less asymmetry relative to wind direction compared with Cox-Munk. Effect of spatial resolution (“smearing” of wind fields) and wave growth dependency explored (shape of waves change with age; younger wind waves are steeper and more asymmetric, older waves are more symmetric, sinusoidal). Refractive index data still an issue, as well as the salinity/chlorinity corrections to fresh water from Friedman (1969).

Further work: 

Further work Investigate impact of JCSDA IRSSE model in the GDAS. Initial tests with the new model show more data is making it past quality control. Further validation of the model with measurements. AERI measurements from 1995 field experiment show that the new model is better at larger angles. More AERI measurements from the CSP tropical western Pacific cruise (1996) will be used for further validation. Investigation of using bicubic spline interpolation to extract IRSSE data from wind speed/view angle database. Surface of emissivities as a function of wind speed and view angle is very smooth, so fit equation may be overkill. Investigation of integration accuracy issue. A very few frequency/wind speed/view angle combinations in the emissivity spectra calculations have shown sensitivity to the integration accuracy over azimuth angle. Solved by higher integration accuracy, but at a computational cost.

Extra Stuff: 

Extra Stuff

TOA TB Residuals for NOAA-17 HIRS.MAX for all wind speeds: 

TOA TB Residuals for NOAA-17 HIRS. MAX for all wind speeds

TOA TB Residuals for AIRS 281 subset.MAX for all wind speeds: 

TOA TB Residuals for AIRS 281 subset. MAX for all wind speeds

TOA TB Residuals for NOAA-17 HIRS.MAX for all wind speeds; only 0ms-1  predicted: 

TOA TB Residuals for NOAA-17 HIRS. MAX for all wind speeds; only 0ms-1  predicted

TOA TB Residuals for AIRS 281 subset.MAX for all wind speeds; only 0ms-1  predicted: 

TOA TB Residuals for AIRS 281 subset. MAX for all wind speeds; only 0ms-1  predicted

Integration accuracy (1): 

Integration accuracy (1) It was noticed that anomalous “bumps” appeared in some coefficients. AIRS module 8 (M8) was affected most. Caused by integration accuracy in code that produces the emissivity spectra. Lower limit of integration over azimuth angle is determined by the accuracy, . In most cases  = 10-5 was sufficient.  = 10-6 was used for all computation except for frequencies around 880cm-1 where  = 10-7 was needed. Lower accuracy == Faster computation For the affected frequencies/wind speeds at a single angle, computation time increased from 6m30s to 4h03m18s!

Integration accuracy (2): 

Integration accuracy (2) Note 6ms-1 results AIRS M8 (~850-900cm-1) coefficients

Integration accuracy (3): 

Integration accuracy (3) E.g.: AIRS M8 ch700 (880.409cm-1) Note anomalous values at 6ms-1. For all affected channels, it’s caused by one “bad” point in the emissivity spectra.

Integration accuracy (4): 

Integration accuracy (4)

Integration accuracy (5): 

Integration accuracy (5) It is not clear why computed emissivities at certain frequencies/wind speeds/angles are sensitive to the integration accuracy. May be due in part to limited precision of the refractive index and salinity/chlorinity correction data – these are functions of frequency only. So, one would think this should affect results at more than a few isolated wind speeds and view angles. Effect of anomalous model coefficients produces an emissivity error of ~0.0003. This is small (effect on TB is also small), but is about 2x the typical RMS emissivity residual.