# Presentation: 861102

## Atmospheric Temperature, Moisture, Ozone, 

**and Motion – Infrared**
**(MOD-07)**
**Jun Li, Liam Gumley, Suzanne Wetzel-Seemann,**
**Chris Moeller, Steve Ackerman, Richard Frey,**
**Dave Santek, Jeff Key, Chris Velden,**
**and Paul Menzel**
**18 Dec 2001**

## VIIRS, MODIS, FY-1C, AVHRR

- VIIRS, MODIS, FY-1C, AVHRR

- H2O

- H2O

- H2O

- H2O

- O2

- CO2

- CO2

- H2O

- O2

- O3

- H2O

- O2

## Earth emitted spectra overlaid on Planck function envelopes

- H20

- CO2

- O3

- CO2

**Notes:**
This slide shows an observed infrared spectrum of the earth thermal emission of radiance to space.  The earth surface Planck blackbody - like radiation at 295 K is severely attenuated in some spectral regions.  Around the absorbing bands of the constituent gases of the atmosphere (CO2 at 4.3 and 15.0 um, H20 at 6.3 um, and O3 at 9.7 um), vertical profiles of atmospheric parameters can be derived.  Sampling in the spectral region at the center of the absorption band yields radiation from the upper levels of the atmosphere (e.g. radiation from below has already been absorbed by the atmospheric gas); sampling in spectral regions away from the center of the absorption band yields radiation from successively lower levels of the atmosphere.  Away from the absorption band are the windows to the bottom of the atmosphere. Surface temperatures of 296 K are evident in the 11 micron window region of the spectrum and tropopause emissions of 220 K in the 15 micron absorption band.  As the spectral region moves toward the center of the CO2 absorption band, the radiation temperature decreases due to the decrease of temperature with altitude in the lower atmosphere.
IR remote sensing (e.g. HIRS and GOES Sounder) currently covers the portion of the spectrum that extends from around 3 microns out to about 15 microns. Each measurement from a given field of view (spatial element) has a continuous spectrum that may be used to analyze the earth surface and atmosphere.  Until recently, we have used “chunks” of the spectrum (channels over selected wavelengths) for our analysis. In the near future, we will be able to take advantage of the very high spectral resolution information contained within the 3-15 micron portion of the spectrum.  From the polar orbiting satellites, horizontal resolutions on the order of  10 kilometers will be available, and depending on the year, we may see views over the same area as frequently as once every 4 hours (assuming 3 polar satellites with interferometers).  With future geostationary interferometers, it may be possible to view at 4 kilometer resolution with a repeat frequency of once every 5 minutes to once an hour, depending on the area scanned and spectral resolution and signal to noise required for given applications. 

## 	Atmospheric Profile Retrieval from MODIS Radiances

- 	                                            ps
- 	I  =  sfc B(T(ps)) (ps) -     B(T(p)) [ d(p) / dp ]  dp .
- 	                                           o
- 	I1, I2, I3, .... , In are measured with MODIS
- 	P(sfc) and T(sfc) come from ground based conventional observations
- 	(p) are calculated with physics models
**	**
**	**Regression relationship is inferred from (1) global set of in situ radiosonde
- 	   reports, (2) calculation of expected radiances, and (3) statistical
- 	   regression of observed raob profiles and calculated MODIS radiances
**	**
**	**Need RT model, estimate of sfc, and MODIS radiances

## MODIS bands 20-29

- MODIS bands 30-36

## Slide 6

## Slide 7

## Slide 8

## Slide 9

## MODIS Emissive Band 

- Cal/Val from ER-2 Platform

-  Transfer S-HIS cal to MAS
-  Co-locate MODIS FOV on MAS
-  Remove spectral, geometric dependence
-  WISC-T2000, SAFARI-2000, TX-2001

**Notes:**
We’ve collected ER-2 and Teraa data during:
1) WISC-T2000  (Spring 2000)
2) SAFARI-2000 (Fall 2000)
3) TX-2001  (Spring 2001)
Since MAS is basically a 1 K accuracy instrument, the SHIS plays a very important role by allowing the 0.5 K or better accuracy of SHIS to be transferred to MAS.

The MAS 50 m resolution data allows the MODIS footprint to be accurately mapped over MAS (SHIS 2 km footprint is bigger than MODIS 1 km)

A forward model (with radiosonde atmospheric characterization) is used to generate spectral (different bandpass) and viewing geometry (altitude difference) corrections.

This is the only airborne effort to validate EOS MODIS thermal band radiances. The approach will be applied to other new sensors (e.g. MODIS on Aqua, VIIRS).

## SHIS calibration transfer to MAS on March 22, 2001

- MAS BT

- MAS BT

- SHIS BT

- SHIS BT

- 3.90um

- 8.6um

- 11.0um

- 13.90um

## Slide 12

## Slide 13

**Notes:**
MODIS B36 and MAS B50 peak at different levels.  This creates sensitivity to atmospheric correction.

There is important atmospheric contribution above the ER-2 level in MODIS Ozone (30) and CO2 bands (35 and 36).

## O3

- Influence of Altitude Difference between MODIS and MAS

- Atmospheric absorption above
- the ER-2 altitude (20 km) is
- important for O3 and CO2
- sensitive bands.

- O3

- CO2

## Switch to 

- Switch to
- Side B

- S/MWIR bias
- adjustment to
- 79/110

- Side A

- MODIS
- first light
- Side A

- WISC-T2000

- SAFARI-2000

- TX-2001

- S/MWIR bias
- at 79/190

**Notes:**
MAS and SHIS missions under Terra have been conducted during 3 field activities since 1st light.
MODIS has been in a stable configuration since Nov 1, 2000.
I’ll speak primarily today about data collected during TX-2001, the most recent activity.

## Calibration

- MODIS - SHIS
- > 0 MODIS too warm
- < 0 MODIS too cold

- SpecNoise

- Residual

## MODIS NEdR Estimate

- Based on Earth Scene Data Day 01153, 20:10 UTC Clear scenes of the Pacific Ocean
- Note: Some SG present in MWIR Used 150 x 28 box (420 data points per detector)

- Band 20		3.7 um		.007 mW/m2/ster/cm-1
- Band 21 	3.9		.02
- Band 22 	3.9		.04
- Band 23 	4.0		.025
- Band 24 	4.45		.03
- Band 25 	4.5		.045
- Band 27		6.7		.08
- Band 28		7.3		.07
- Band 29		8.6		.25
- Band 30		9.7		.2
- Band 31		11.0		.3
- Band 32 	12.0		.3
- Band 33 	13.3		.4
- Band 34 	13.6		.6
- Band 35		13.9		.4
- Band 36		14.2		.5

## Band 34

- Noisy
- Detectors

## Detector to detector calibration

## Baja

***A known optical leak at 11um caused the image of the Baja peninsula to be present in MODIS 14.3um data.  Through testing, the pre-launch correction coefficients were revised, removing the contamination. ***

- 11um

- 14.3um

- 14.3um

- Pre launch correction

- Post launch  correction

**Correcting Crosstalk**

## 11 micron

- GOES 3 by 3 FOVs (30 km)

- MODIS 5 by 5 FOVs (5 km)

## 	Atmospheric Profile Retrieval from MODIS Radiances

- 	                                            ps
- 	I  =  sfc B(T(ps)) (ps) -     B(T(p)) [ d(p) / dp ]  dp .
- 	                                           o
- 	I1, I2, I3, .... , In are measured with MODIS
- 	P(sfc) and T(sfc) come from ground based conventional observations
- 	(p) are calculated with physics models
**	**
**	**Regression relationship is inferred from (1) global set of in situ radiosonde
- 	   reports, (2) calculation of expected radiances, and (3) statistical
- 	   regression of observed raob profiles and calculated MODIS radiances
**	**
**	**Need RT model, estimate of sfc, and MODIS radiances

## The MODIS AP algorithms are based on a regression procedure, and makes use of the NOAA-88 data set containing more than 7500 global profiles of temperature and moisture to determine the regression coefficients.  The radiative transfer calculation of the MODIS spectral band radiances is performed for each training profile using the Pressure layer Fast Algorithm for Atmospheric Transmittances (PFAAST) transmittance model.  This model has 101 pressure layer vertical coordinates from 0.1 to 1050 hPa and takes into account the satellite zenith angle, absorption by well-mixed gases (including nitrogen, oxygen, and carbon dioxide), water vapor, and ozone. The MODIS instrument noise is added into the calculated spectral band radiances, and these radiative transfer calculations provide a temperature-moisture-ozone profile and MODIS radiance pair for use in the statistical regression analysis.  

## 30km resolution GOES

- 500hPa T

- 5km resolution MODIS

## MODIS 2000/09/05-08 

**MODIS 2000/09/05-08**** **
- Daytime 850 hPa Temperature (K) for 4 days

## NOAA-15 AMSU-A 2000/09/05

**NOAA-15 AMSU-A 2000/09/05**
**Daytime 850 hPa Temperature (K) for one day**

## MODIS

**MODIS**
**AMSU**

**850 mb**
**24 Aug 00**

## MODIS

**MODIS**
**AMSU**

**500 mb**
**24 Aug 00**

## MODIS

**MODIS**
**AMSU**

**50 mb**
**24 Aug 00**

## 4m Surface Reflection must be considered 

- 4m Surface Reflection must be considered

- December 10, 2001               						    	     MODIS Meeting

**Notes:**
This figure compares MODIS, GOES, and Radiosonde TPW to the microwave radiometer TPW at the SGP-CART site for 25 cases between 15 March - 15 November, 2001.  These cases were for clear-sky conditions when the MODIS viewing zenith angle was < 32 degrees over the SGP-CART site.  

The left column shows a scatter-plot comparison of MODIS, GOES, and Radiosonde TPW (y-axis) vs. Microwave radiometer (x-axis).  The right column shows the same data, but shown as a difference, microwave radiometer minus (MODIS, GOES, or radiosonde). 

##            Mitigating Problems in MODIS AP Algorithm

**           ****Mitigating Problems in MODIS AP Algorithm**

- * Changed predictors band 24 and 25 (4.4 and 4.5 um) brightness temperatures to their BT difference to remove surface effects:	                    	Old algorithm had 12 predictors:
- 		individual bands 24, 25, and 27 through 36.
- 	New algorithm has 11 predictors:
- 		band 25 - 24 BT difference and
- 		individual bands 27 through 36.
- * Implemented *preliminary* radiance bias corrections.  Currently based only on CART site data, global biases are in progress.
- * Applied post-launch NEdT in place of pre-launch.

- December 18, 2001               						         MODIS Science Team Meeting

**Notes:**
1) Changed the predictors used for the statistical regression: The operational algorithm uses 12 individual infrared 
bands, with wavelengths 4.4-14.2m (MODIS bands 24,25,27-36).  The new algorithm uses these same 12     
bands but uses the difference between the 4.4 and 4.5m bands (MODIS 24, 25) instead of the individual bands.
2) Separated the statistical training data set into 5 brightness temperature classes: < 254, 254-275, 275-286, 286-
297, > 297 oK.  For example, a retrieval for a scene with MODIS brightness temperature < 254 oK uses only the 
subset of the training data for which a forward model calculation from the profile predicts a brightness 
temperature < 254 oK.
3) Updated the estimates of NEdT: Operational NEdT was based on pre-launch values.  New, post-launch estimates were computed and are used in the new algorithm.
4) Implemented a radiance bias correction: The operational MODIS TPW algorithm does not include any radiance 
bias correction.  The new radiance bias adjustments account for errors introduced in both the measurements and the algorithm, including the radiative transfer calculations.  

## New

- Operational

- 850

- 500

- 300

**Notes:**
This figure compares MODIS, GOES, and Radiosonde TPW to the microwave radiometer TPW at the SGP-CART site for 25 cases between 15 March - 15 November, 2001.  These cases were for clear-sky conditions when the MODIS viewing zenith angle was < 32 degrees over the SGP-CART site.  

The left column shows a scatter-plot comparison of MODIS, GOES, and Radiosonde TPW (y-axis) vs. Microwave radiometer (x-axis).  The right column shows the same data, but shown as a difference, microwave radiometer minus (MODIS, GOES, or radiosonde). 

## GOES 30 km resolution

- MODIS 5 km resolution

- TPW

- TPW

**GOES vs. MODIS 2000/06/30 1600 UTC**
**Total Precipitable Water (mm)**

## MODIS total precipitable water vapor shows a wet bias wrt GOES;

**MODIS total precipitable water vapor shows a wet bias wrt GOES;**
**bias 1.5 mm and rms of 3 mm; bias will be removed after more validation**

## MODIS 2000/09/05-08 

**MODIS 2000/09/05-08**** **
**Daytime Total Precipitable Water (cm)**
**values over land not shown to facilitate comparison with AMSU**

## Other MOD07 Products: Water Vapor High (mm)

- Operational

- New

**Notes:**
This figure compares MODIS, GOES, and Radiosonde TPW to the microwave radiometer TPW at the SGP-CART site for 25 cases between 15 March - 15 November, 2001.  These cases were for clear-sky conditions when the MODIS viewing zenith angle was < 32 degrees over the SGP-CART site.  

The left column shows a scatter-plot comparison of MODIS, GOES, and Radiosonde TPW (y-axis) vs. Microwave radiometer (x-axis).  The right column shows the same data, but shown as a difference, microwave radiometer minus (MODIS, GOES, or radiosonde). 

## GOES 30 km resolution

- MODIS 5 km resolution

- Ozone

- Ozone

## MODIS ozone is very close to the GOES ozone (over North America); 

**MODIS ozone is very close to the GOES ozone (over North America); **
**rms of about 10 Dobsons; polar extreme ozone values will be improved**

## MODIS 2000/09/05-08

**MODIS 2000/09/05-08**
**Daytime Total Ozone (Dobsons)**

## Earth Probe TOMS 2000/09/05

**Earth Probe TOMS 2000/09/05**
**Total Ozone (Dobsons)**

## MODIS Ozone

- TOMS Ozone

## Algorithm Change effect on Total Column Ozone (dob)

- Operational

- “New”

**Notes:**
This figure compares MODIS, GOES, and Radiosonde TPW to the microwave radiometer TPW at the SGP-CART site for 25 cases between 15 March - 15 November, 2001.  These cases were for clear-sky conditions when the MODIS viewing zenith angle was < 32 degrees over the SGP-CART site.  

The left column shows a scatter-plot comparison of MODIS, GOES, and Radiosonde TPW (y-axis) vs. Microwave radiometer (x-axis).  The right column shows the same data, but shown as a difference, microwave radiometer minus (MODIS, GOES, or radiosonde). 

## A Preview of ABI High Resolution Water Vapor Imagery1 km MODIS & 4x8 km GOES

## Four Panel Zoom of Cloud-Free Orographic Waves revealed in Water Vapor Imagery

**Notes:**
These events (mountain waves and gravity waves were detected in clear skies.

## Every

**Every**
**100 min**
**MODIS**
**covers**
**polar**
** ****regions**

- x

**Notes:**
Polar orbiters revisit the poles every 100 minutes.  Loops of these images can reveal polar
motions undetected by other systems. 

## 

- x

## Winds from MODIS: An Arctic Example

- Cloud-track winds (left) and water vapor winds (right) from MODIS for a case in the western Arctic.  The wind vectors were derived from a sequence of three images, each separated by 100 minutes.  They are plotted on the first 11 m (left) and 6.7 m (right) images in the sequence.

## Orbital Issues

- How often can wind vectors be obtained from a polar-orbiting satellites? The figure below shows the time of successive overpasses at a given latitude-longitude point on a single day with only the Terra satellite.  The figure at the upper right shows the frequency of "looks" by two satellites: Terra and (the future) Aqua.  The figure at the lower right shows the temporal sampling with five satellites.

**Notes:**
Example of reading the plots: In the first plot (lower left), at 60o latitude there are only two views in a 24-hour period, and they are separated by more than 10 hours.  At 80o there are many views separated by approximately 100 minutes (the orbital period), but there is still a period without any coverage for 13-14 hours  each day. 

## Early Estimates of UW MODIS Atmospheric Products Quality

- MODIS IR radiances agree to within 1.5 C with GOES and ER-2 MAS/SHIS
- MODIS layer tropospheric temperatures compare well with AMSU;
- 	rms better than 1 C, both within 2 C of radiosonde observations;
- 	validated with know characteristics
- 	CART site validation ongoing.
- MODIS layer dewpoint temperatures depict gradients very well;
- 	are within 3-4 C of radiosonde observations.
- MODIS ozone is very close to the GOES ozone (over North America);
- 	rms of about 10 Dobsons;
- 	polar extreme ozone values will be improved with more training data.
- 	validated with known characteristics.
- MODIS polar winds represent coherent atmospheric motion;
- 	model assimilation underway;
- 	geo-like quality observed;
- 	within 7 – 10 m/s of the few raobs available for validation.

## Total Precipitable Water Vapor – Infrared(MOD-05)

**Suzanne Wetzel-Seemann, Jun Li, Liam Gumley, **
**and Paul Menzel**
**18 Dec 2001**

## VIIRS, MODIS, FY-1C, AVHRR

- VIIRS, MODIS, FY-1C, AVHRR

- H2O

- H2O

- H2O

- H2O

- O2

- CO2

- CO2

- H2O

- O2

- O3

- H2O

- O2

## Earth emitted spectra overlaid on Planck function envelopes

- H20

- CO2

- O3

- CO2

**Notes:**
This slide shows an observed infrared spectrum of the earth thermal emission of radiance to space.  The earth surface Planck blackbody - like radiation at 295 K is severely attenuated in some spectral regions.  Around the absorbing bands of the constituent gases of the atmosphere (CO2 at 4.3 and 15.0 um, H20 at 6.3 um, and O3 at 9.7 um), vertical profiles of atmospheric parameters can be derived.  Sampling in the spectral region at the center of the absorption band yields radiation from the upper levels of the atmosphere (e.g. radiation from below has already been absorbed by the atmospheric gas); sampling in spectral regions away from the center of the absorption band yields radiation from successively lower levels of the atmosphere.  Away from the absorption band are the windows to the bottom of the atmosphere. Surface temperatures of 296 K are evident in the 11 micron window region of the spectrum and tropopause emissions of 220 K in the 15 micron absorption band.  As the spectral region moves toward the center of the CO2 absorption band, the radiation temperature decreases due to the decrease of temperature with altitude in the lower atmosphere.
IR remote sensing (e.g. HIRS and GOES Sounder) currently covers the portion of the spectrum that extends from around 3 microns out to about 15 microns. Each measurement from a given field of view (spatial element) has a continuous spectrum that may be used to analyze the earth surface and atmosphere.  Until recently, we have used “chunks” of the spectrum (channels over selected wavelengths) for our analysis. In the near future, we will be able to take advantage of the very high spectral resolution information contained within the 3-15 micron portion of the spectrum.  From the polar orbiting satellites, horizontal resolutions on the order of  10 kilometers will be available, and depending on the year, we may see views over the same area as frequently as once every 4 hours (assuming 3 polar satellites with interferometers).  With future geostationary interferometers, it may be possible to view at 4 kilometer resolution with a repeat frequency of once every 5 minutes to once an hour, depending on the area scanned and spectral resolution and signal to noise required for given applications. 

## 	Total Water Vapor Retrieval from MODIS IR Radiances

- 	                                            ps
- 	I  =  sfc B(T(ps)) (ps) -     B(T(p)) [ d(p) / dp ]  dp .
- 	                                           o
- 	I1, I2, I3, .... , In are measured with MODIS
- 	P(sfc) and T(sfc) come from ground based conventional observations
- 	(p) are calculated with physics models
**	**
**	**Regression relationship is inferred from (1) global set of in situ radiosonde
- 	    reports, (2) calculation of expected radiances, and (3) statistical non-linear
- 	    regression of observed Raob TPW and calculated MODIS radiances
- 	   (brightness temperatures)
**	**
**	**Need RT model, estimate of sfc, and MODIS radiances

## The MODIS TPW algorithm is based on a regression procedure, and makes use of the NOAA-88 data set containing more than 7500 global profiles of temperature and moisture to determine the regression coefficients.  The radiative transfer calculation of the MODIS spectral band radiances is performed for each training profile using the Pressure layer Fast Algorithm for Atmospheric Transmittances (PFAAST) transmittance model.  This model has 101 pressure layer vertical coordinates from 0.1 to 1050 hPa and takes into account the satellite zenith angle, absorption by well-mixed gases (including nitrogen, oxygen, and carbon dioxide), water vapor, and ozone. The MODIS instrument noise is added into the calculated spectral band radiances, and these radiative transfer calculations provide a temperature-moisture-ozone profile and MODIS radiance pair for use in the statistical regression analysis.  

## Total Water Vapour can be evaluated in multiple infrared window channels where absorption is weak, so that 

- 		w  =  exp[- kwu] ~  1 - kwu where w denotes window channel
- and
- 		dw  =  - kwdu
- What little absorption exists is due to water vapour, therefore, u is a measure of precipitable water vapour.  RTE in window region
- 				         us
- 		Iw  =  Bsw  (1-kwus) +  kw        Bwdu
-                                 		       o
- us represents total atmospheric column absorption path length due to water vapour, and s denotes surface.  Defining an atmospheric mean Planck radiance, then
- 			        _ 	        _        us	  us
- 	Iw  =  Bsw (1-kwus) + kwusBw  with  Bw  =    Bwdu /   du
- 		         			  o	  o
- Since Bsw is close to both Iw and Bw, first order Taylor expansion about the surface temperature Ts allows us to linearize the RTE with respect to temperature, so
-                                                                 _
- 	     Tbw  =  Ts (1-kwus)  +  kwusTw , where Tw is mean atmospheric temperature corresponding to Bw.

## For two window channels (11 and 12um) the following ratio can be determined.

- 					   _
- 		  Ts - Tbw1      	 kw1us(Ts -  Tw1)	    kw1
- 		_________ =         ______________   =   ___       _                                         _
- 	 	   	  		  _
- 	  	  Ts - Tbw2	                 kw1us(Ts  - Tw2)          kw2
- where the mean atmospheric temperature measured in the one window region is assumed to be comparable to that measured in the other, Tw1 ~ Tw2,
- Thus  it follows that
- 			              kw1
- 		Ts  =  Tbw1 +     _                       _	[Tbw1 - Tbw2]
- 		      	          kw2 - kw1
- and
- 			  Tbw - Ts
- 		us     =    _                     _     .
-    	      	       	      _
- 			kw (Tw - Ts)
- Obviously, the accuracy of the determination of the total water vapour concentration depends upon the contrast between the surface temperature, Ts, and
- 				       _
- the effective temperature of the atmosphere Tw

## GOES vs. MODIS 2000/06/30 1600 UTC

- Total Precipitable Water (mm)

## MODIS total precipitable water vapor shows a wet bias wrt GOES;

**MODIS total precipitable water vapor shows a wet bias wrt GOES;**
**bias 1.5 mm and rms of 3 mm; bias will be removed after more validation**

##                MODIS TPW June 2, 2001: Operational Algorithm

**               ****MODIS TPW June 2, 2001: Operational Algorithm**

- December 10, 2001               							         MODIS Meeting

**SWIR daytime reflection causing wet atm over African desert**

**Notes:**
Operational MODIS TPW product:
Notice the large area of very high TPW over desert areas.   

##            Mitigating Problems in MODIS TPW Algorithm

**           ****Mitigating Problems in MODIS TPW Algorithm**

- * Changed predictors band 24 and 25 (4.4 and 4.5 um) brightness temperatures to their BT difference to remove surface effects:	                    	Old algorithm had 12 predictors:
- 		individual bands 24, 25, and 27 through 36.
- 	New algorithm has 11 predictors:
- 		band 25 - 24 BT difference and
- 		individual bands 27 through 36.
- *  Separated training into five regression BT zones to include a broader range of moisture regimes
- * Implemented *preliminary* radiance bias corrections.  Currently based only on CART site data, global biases are in progress.
- * Applied post-launch NEdT in place of pre-launch.

- December 18, 2001               						         MODIS Science Team Meeting

**Notes:**
1) Changed the predictors used for the statistical regression: The operational algorithm uses 12 individual infrared 
bands, with wavelengths 4.4-14.2m (MODIS bands 24,25,27-36).  The new algorithm uses these same 12     
bands but uses the difference between the 4.4 and 4.5m bands (MODIS 24, 25) instead of the individual bands.
2) Separated the statistical training data set into 5 brightness temperature classes: < 254, 254-275, 275-286, 286-
297, > 297 oK.  For example, a retrieval for a scene with MODIS brightness temperature < 254 oK uses only the 
subset of the training data for which a forward model calculation from the profile predicts a brightness 
temperature < 254 oK.
3) Updated the estimates of NEdT: Operational NEdT was based on pre-launch values.  New, post-launch estimates were computed and are used in the new algorithm.
4) Implemented a radiance bias correction: The operational MODIS TPW algorithm does not include any radiance 
bias correction.  The new radiance bias adjustments account for errors introduced in both the measurements and the algorithm, including the radiative transfer calculations.  

## MODIS Radiance (BT) Bias wrt CART Site

**MODIS Radiance (BT) Bias wrt CART Site**

## MODIS NEdR Estimate

- Based on Earth Scene Data Day 01153, 20:10 UTC Clear scenes of the Pacific Ocean
- Note: Some SG present in MWIR Used 150 x 28 box (420 data points per detector)

- Band 20		3.7 um		.007 mW/m2/ster/cm-1
- Band 21 	3.9		.02
- Band 22 	3.9		.04
- Band 23 	4.0		.025
- Band 24 	4.45		.03
- Band 25 	4.5		.045
- Band 27		6.7		.08
- Band 28		7.3		.07
- Band 29		8.6		.25
- Band 30		9.7		.2
- Band 31		11.0		.3
- Band 32 	12.0		.3
- Band 33 	13.3		.4
- Band 34 	13.6		.6
- Band 35		13.9		.4
- Band 36		14.2		.5

##                 MODIS TPW June 2, 2001: 

-                 **MODIS TPW June 2, 2001: **
**             ****Operational Algorithm (IR ch 24 –> 36)**

**Notes:**
Operational MODIS TPW product:
Notice the large area of very high TPW over desert areas.   

##              MODIS TPW June 2, 2001: 

-              **MODIS TPW June 2, 2001: **
**           ****New Algorithm 1 (Ch 24 -25)**

**Notes:**
New Algorithm:
By using the difference between the 4.4 and 4.5 micron bands (24 & 25) instead of the individual bands, we have fixed the desert problem by removing the surface contribution of these shortwave infrared bands.  The bias correction and post-launch NeDT also changed the TPW but to a far lesser degree.

##              MODIS TPW June 2, 2001: 

-              **MODIS TPW June 2, 2001: **
**           ****New Algorithm 2 (5 zones)**

**Notes:**
New Algorithm #2:
This algorithm uses the same brightness temperature difference as on the previous slide but it also divides the training data set into 5 11micron brightness temperature groups.  This version is superior to the version on the previous page in most areas, with the exception of the noise introduced over the North African desert.  The noise is not too apparent in this global image.  We will probably be able to overcome the noise problem by adding more training data in the very high brightness temperature zone.

##                    MODIS TPW June 2, 2001: 

-                    **MODIS TPW June 2, 2001: **
**              ****New Algorithm 3 (5 overlapping zones)**

**Notes:**
New Algorithm #3:
This algorithm uses the brightness temperature difference in bands 24-25 for *land only* and uses the operational algorithm over ocean.  The oceans are moister here than in the operational image because I used a constant factor = 1.0 for both land and ocean.  Operational has a factor = 3.0.  The “factor” is the amount that the input NEdT is increased.  If we decide we like the approach of separating out the land/ocean then I could run a version with land factor = 1.0 and ocean factor = 3.0 to bring down the moisture over the tropics.

## 1

- 2
- 3

## MODIS TPW and GDAS TPW

- MODIS TPW and GDAS TPW

## New MODIS TPW Algorithm: Comparison with NOAA-15 Advanced Microwave Sounding Unit (AMSU) for June 2, 2001

**New MODIS TPW Algorithm: Comparison with NOAA-15 Advanced Microwave Sounding Unit (AMSU) for June 2, 2001**

- TPW (cm)

- MODIS new TPW (cm)
- June 2, 2001

- AMSU TPW (cm)
- June 2, 2001

**Notes:**
New MODIS algorithm (version #2 shown here) agrees quite well with the AMSU TPW.  AMSU can measure TPW even in the presence of cloud, MODIS measures only for clear sky.  MODIS image on this slide made by Liam Gumley.  

## Improvement in Desert TPW: North Africa

- TPW (mm)

**Improvement in Desert TPW: North Africa**
**June 2, 2001: 0830, 0835, 1010, 1015, 1150, 1155 UTC **

- Operational MOD07 TPW

- New MOD07 TPW

**Notes:**
“New” MODIS shown here is version #1

## Slide 71

## Comparison with MOD05: North Africa

- TPW (mm)

**Comparison with MOD05: North Africa**
**June 2, 2001: 0830, 0835, 1010, 1015, 1150, 1155 UTC**

- New MOD07 TPW

- MOD05 TPW

## GOES 8 & 10

- TPW (mm)

- June 2, 2001:1800 UTC

- MOD 07 June 2, 2001:
- 1640, 1645, 1820, 1825 UTC

- MOD 05  June 2, 2001:
- 1640, 1645, 1820, 1825 UTC

## CART Site TPW Comparison: Sample of One Case

- December 18, 2001               						 MODIS Science Team Meeting

**Notes:**
This is an example of the CART site TPW comparison for one case.  Results from 25 cases combined are shown on the next slide.  This compares the MODIS, GOES, and radiosonde TPW at the CART site to the microwave radiometer (MWR) TPW co-located in time with each observation.  For example, using the case on this slide: radiosonde TPW was compared with MWR at about 17:39 UTC, MODIS at 17:43, and GOES at 18:00 UTC. 

## CART Site TPW Comparison: 27 Cases 

- CART Site TPW Comparison: 27 Cases
- March 15 - December 2, 2001

- Operational
-   Algorithm

- New
- Algorithm

**MWR - MODIS = 0.576 mm**
***for TPW < 15mm:***
**MWR - MODIS = - 2mm**
***for TPW > 15mm:***
**MWR - MODIS = + 3.4mm**

**MWR – MODIS = - 4.4mm**
***for TPW < 15mm:***
**MWR - MODIS = - 4.5mm**
***for TPW > 15mm:***
**MWR – MODIS = - 4.4mm**

**Notes:**
This figure compares MODIS, GOES, and Radiosonde TPW to the microwave radiometer TPW at the SGP-CART site for 25 cases between 15 March - 15 November, 2001.  These cases were for clear-sky conditions when the MODIS viewing zenith angle was < 32 degrees over the SGP-CART site.  

The left column shows a scatter-plot comparison of MODIS, GOES, and Radiosonde TPW (y-axis) vs. Microwave radiometer (x-axis).  The right column shows the same data, but shown as a difference, microwave radiometer minus (MODIS, GOES, or radiosonde). 

## MODIS: Product Image Artifacts

** ****Product striping introduced because:**
**     ****- two mirror sides not identical and not characterized perfectly,**
**     ****- each detector calibrated independently,**
**     ****- some detectors "out of family" compared to majority for a given band,**
**     ****- out of band response influencing measured radiance.**
**Destriping approaches being tested:**
**     ****- generate granule overall mean and standard deviation **
**	****for each detector each mirror side;**
**      ****- in regions of low standard deviation, select reference detector **
**	****and compute radiance ratios of different detectors;**
**      ****- scale other detectors using computed radiance ratio.**
**	**

## Slide 77

## Original L1B (V003)

*Original L1B (V003)*

*Destriped*

*MODIS Band 27 (6.7 ****m), 2001-06-04 16:45 UTC*

## Early Estimates of UW MODIS IR TPW Product Quality

**MODIS IR total precipitable water vapor **
**	***** captures TPW gradients very well over oceans**
**	***** has been improved over daytime non vegetated land**
**	***** is very sensitive to multi-spectral striping**
**	***** shows ****	****a wet bias wrt microwave for TPW < 15mm, **
**	****	****a dry bias for TPW > 15mm**
**	***** bias will be removed after more validation. **
**MODIS and GOES TPW agree well with rms difference of 3 mm**
**MODIS TPW**
**	***** best IR TPW is five zone split window regression**
**	****	****on de-striped data****	**
**	***** vis-NIR TPW also available over daytime land**
**STATUS – Validated with known characteristics**

## Cloud Top Properties – Infrared

**(MOD-06)**
**Richard Frey, Bryan Baum, Kathy Strabala,**
**Shaima Nasiri, Hong Zhang, Jun Li, **
**and Paul Menzel**
**18 Dec 2001**

## Cloud Top Properties

*Objective*: Retrieve cloud top properties for every 5 x 5 box of 1000 m FOVs where at least 5 FOVs are cloudy.
*Method*: Longwave infrared CO2 slicing; tri-spectral for cloud phase.
*Products*: Cloud top pressure and temperature; Cloud fraction and effective emissivity; Cloud phase (water, ice, mixed).
- Frey, R. A. et al, 1999: A comparison of cloud top heights computed from airborne lidar and MAS radiance data. *J. Geophys. Res.*, 104, 24547-24555.

## VIIRS, MODIS, FY-1C, AVHRR

- VIIRS, MODIS, FY-1C, AVHRR

- H2O

- H2O

- H2O

- H2O

- O2

- CO2

- CO2

- H2O

- O2

- O3

- H2O

- O2

## Earth emitted spectra overlaid on Planck function envelopes

- H20

- CO2

- O3

- CO2

**Notes:**
This slide shows an observed infrared spectrum of the earth thermal emission of radiance to space.  The earth surface Planck blackbody - like radiation at 295 K is severely attenuated in some spectral regions.  Around the absorbing bands of the constituent gases of the atmosphere (CO2 at 4.3 and 15.0 um, H20 at 6.3 um, and O3 at 9.7 um), vertical profiles of atmospheric parameters can be derived.  Sampling in the spectral region at the center of the absorption band yields radiation from the upper levels of the atmosphere (e.g. radiation from below has already been absorbed by the atmospheric gas); sampling in spectral regions away from the center of the absorption band yields radiation from successively lower levels of the atmosphere.  Away from the absorption band are the windows to the bottom of the atmosphere. Surface temperatures of 296 K are evident in the 11 micron window region of the spectrum and tropopause emissions of 220 K in the 15 micron absorption band.  As the spectral region moves toward the center of the CO2 absorption band, the radiation temperature decreases due to the decrease of temperature with altitude in the lower atmosphere.
IR remote sensing (e.g. HIRS and GOES Sounder) currently covers the portion of the spectrum that extends from around 3 microns out to about 15 microns. Each measurement from a given field of view (spatial element) has a continuous spectrum that may be used to analyze the earth surface and atmosphere.  Until recently, we have used “chunks” of the spectrum (channels over selected wavelengths) for our analysis. In the near future, we will be able to take advantage of the very high spectral resolution information contained within the 3-15 micron portion of the spectrum.  From the polar orbiting satellites, horizontal resolutions on the order of  10 kilometers will be available, and depending on the year, we may see views over the same area as frequently as once every 4 hours (assuming 3 polar satellites with interferometers).  With future geostationary interferometers, it may be possible to view at 4 kilometer resolution with a repeat frequency of once every 5 minutes to once an hour, depending on the area scanned and spectral resolution and signal to noise required for given applications. 

## MODIS 

**MODIS **
**identifies**
**cloud**
**classes**

**Clear**
**Snow**
**Low Cloud**
**Mid-low **
**   ****Cloud**
**High Cloud**

## Clouds separate into classes 

**Clouds separate into classes **
**when multispectral radiance information is viewed**

## MODIS Emissive Band 

- Cal/Val from ER-2 Platform

-  Transfer S-HIS cal to MAS
-  Co-locate MODIS FOV on MAS
-  Remove spectral, geometric dependence
-  WISC-T2000, SAFARI-2000, TX-2001

**Notes:**
We’ve collected ER-2 and Teraa data during:
1) WISC-T2000  (Spring 2000)
2) SAFARI-2000 (Fall 2000)
3) TX-2001  (Spring 2001)
Since MAS is basically a 1 K accuracy instrument, the SHIS plays a very important role by allowing the 0.5 K or better accuracy of SHIS to be transferred to MAS.

The MAS 50 m resolution data allows the MODIS footprint to be accurately mapped over MAS (SHIS 2 km footprint is bigger than MODIS 1 km)

A forward model (with radiosonde atmospheric characterization) is used to generate spectral (different bandpass) and viewing geometry (altitude difference) corrections.

This is the only airborne effort to validate EOS MODIS thermal band radiances. The approach will be applied to other new sensors (e.g. MODIS on Aqua, VIIRS).

## RTE in Cloudy Conditions

- 	Iλ  =  η Icd + (1 - η) Ic     where cd = cloud, c = clear, η = cloud fraction
-  	             λ                           λ
- 			     o
- 	Ic  =  Bλ(Ts) λ(ps) +     Bλ(T(p)) dλ .
- 	 λ                               ps
- 	                                                       pc
- 	Icd  =  (1-ελ) Bλ(Ts) λ(ps) + (1-ελ)      Bλ(T(p)) dλ
-      	 λ			           ps
- 						         o
-    				 + ελ Bλ(T(pc)) λ(pc) +     Bλ(T(p)) dλ
-                                    				        pc
- ελ is emittance of cloud.  First two terms are from below cloud, third term is cloud contribution, and fourth term is from above cloud.  After rearranging
- 			          pc         dBλ
- 		Iλ - Iλc  =  ηελ      (p)  _        _   dp  .
-     	  		         ps           dp
- Techniques for dealing with clouds fall into three categories: (a) searching for cloudless fields of view, (b) specifying cloud top pressure and sounding down to cloud level as in the cloudless case, and (c) employing adjacent fields of view to determine clear sky signal from partly cloudy observations.

## Cloud Properties

- RTE for cloudy conditions indicates dependence of cloud forcing (observed minus clear sky radiance) on cloud amount () and cloud top pressure (pc)
- 	                              pc
- 	(I - Iclr)  =        dB .
- 	                             ps
- Higher colder cloud or greater cloud amount produces greater cloud forcing; dense low cloud can be confused for high thin cloud.  Two unknowns require two equations.
- pc can be inferred from radiance measurements in two spectral bands where cloud emissivity is the same.   is derived from the infrared window, once pc is known.  This is the essence of the CO2 slicing technique.

## MODIS bands 20-29

- MODIS bands 30-36

## Slide 90

## MAS CO2 heights validated with 4700 Lidar measurements

**MAS CO2 heights validated with 4700 Lidar measurements**

## Slide 92

## MODIS 2000/03/12 1730 UTC Cloud Top Pressure

## VIS     CM

- CTP       N

- MODIS
- Cloud
- Mask

- MODIS
- Cloud
- Properties

- r 3-4
- g 4-5
- b 5-6

- r>95
- b>75
- g>50
- y>25

## Comparison of CLS (nadir view), HIRS (3 hrs later), RAOB, & MODIS Cloud Properties

**Comparison of CLS (nadir view), HIRS (3 hrs later), RAOB, & MODIS Cloud Properties**

## Cloudtop Pressure

- Cloudtop Temperature

- Cloud Fraction

- Cloud Effective Emissivity

- MODIS 2000/09/05-08

## MODIS CO2 heights compared with GOES CO2 heights

**MODIS CO2 heights compared with GOES CO2 heights**

## MODIS CO2 heights compared with GOES CO2 heights

**MODIS CO2 heights compared with GOES CO2 heights**
**High clouds < 400 mb**

- MODIS sees high clouds lower in atmosphere

## MODIS CO2 heights compared with GOES CO2 heights

**MODIS CO2 heights compared with GOES CO2 heights**
**mid level clouds between 400 and 700 mb**

## MODIS CO2 heights compared with GOES CO2 heights

**MODIS CO2 heights compared with GOES CO2 heights**
**Low clouds > 700 mb**

- MODIS sees low clouds higher in atmosphere

## Slide 101

## Tri-spectral IR thermodynamic phase algorithm

- 8.6-11 vs 11-12
- when slope > 1 then ice
- when slope < 1 then water

**ice cloud**
**April 1996**
**Success**

**water cloud**
**Jan 1993**
**TOGA/**
**COARE**

- Strabala, Menzel, and Ackerman, 1994, *JAM*, **2**, 212-229.
- Baum et al, 2000, *JGR*, **105**, 11781-11792.

**Notes:**
Example from MAS where water and ice clouds are separated according to both the absolute BT 
differences as well as the  slope.

Solar refl. technique more problematic because of a greater dependence on particle shape, as 
well as solar geometry

## Multispectral data reveals improved information about ice / water clouds

**Notes:**
A well known application of IR window multispectral is discrimination of water cloud from ice cloud.  Because liquid phase water absorbs radiation differently than solid phase water (ice) at different wavelengths, it will be possible to discriminate between the two, as well as to detect clouds with mixed phase. The tri-spectral technique utilizes radiance measurements near 10.0, 11.5, and 13.0 microns; water is indicated when T10~T11 and T11<<T13 and ice is indicated when T10<<T11 and T11~T13.  

With geostationary tri-spectral observations it will be possible to watch clouds as they change phase, which is often indicative of where they are in their life cycle.

Knowing cloud type  and height (from CO2 slicing) is important for numerous applications - not the least of which is feeding advanced numerical weather prediction models that will soon be available.

## Cloud Thermodynamic Phase

- 21 April, 2001 at 1745Z
- ARM Southern Great Plains Site

- Kansas

- Oklahoma

**Cloud Thermodynamic Phase**

##  Simulations of Ice and Water Phase Clouds

***Ice clouds:***
-  BTD[8.5-11] > 0 over a large range
- of optical thicknesses
- 1.6-m reflectances for ice clouds
- tend to be lower than for water clouds

***Midlevel clouds:***
- BTD[8.5-11] < 0 for both water
- and ice clouds

***Low-level, warm clouds:***
- BTD[8.5-11] always negative
- 1.6-m reflectances for water clouds
- tend to be higher than for ice clouds

- Radiative transfer simulations for a cirrus model derived from in-situ data collected
- during the FIRE-I field campaign (upper two panels) and a water cloud having an
- effective radius of 10 m. Viewing angles are o = 45o,  = 20o, and  = 40o and
- midlatitude summer temperature and relative humidity profiles are assumed.

**Notes:**
Ice cloud calculations: 
Mixture of bullet rosettes, plates, solid columns, and aggregates. The exact mixture was determined
 from in-situ measurements provided by Andy Heymsfield and were from the FIRE I field campaign. 
So, we're simulating fairly thin midlatitude cirrus. The scattering properties for the individual crystals, 
for each size and wavelength, were provided by Ping Yang. Shaima and I worked out how to put all 
these pieces together (in-situ measurements and individual crystal scattering properties) to form 
new cirrus scattering models for a given wavelength. This work has been written up and submitted to 
JAM for publication (written by Shaima with Baum, Andy Heymsfield, Ping Yang, Dave Kratz as 
co-authors) this past December. Details are on her web site. 

Water cloud simulations: 
Gamma distribution of water droplets, with effective radius of 10 microns and a 
variance of 0.1. 

Overall: 
We used a midlatitude summer climatological profile. Left out aerosols altogether. Used Kratz's 
correlated-K distribution routines to account for gaseous absorption. Solar zenith = 45 degrees, 
viewing zenith = 20 degrees, and relative azimuth = 40 degrees. 

## 200 by 200 pixels of  MODIS Data from 15 Oct. 2000 at 1725Z

**200 by 200 pixels of  MODIS Data from 15 Oct. 2000 at 1725Z**

**Low Cloud ****(from MODIS Phase)**

**High Cloud ****(from MODIS Phase)**

**Clear ****(from MODIS Cloud Mask)**

**Other ****(to be determined)**

**11 ********m BT (K)**

**1.6 ********m Reflectance**

**210         230           250            270          290 **

## Cloud Overlap Detection

- Potentially Overlapped Pixels in Red

- 21 April, 2001 at 1745Z
- ARM Southern Great Plains Site

**Cloud Overlap Detection**

## MODIS Frequency of Co-occurrence Water Phase with 253 K < Tcld < 268 K  05 Nov. 2000 - Daytime Only

- frequency of occurrence in percent (%)

**MODIS Frequency of Co-occurrence ****Water Phase with 253 K < T****cld**** < 268 K **** 05 Nov. 2000 - Daytime Only**

## MODIS Frequency of Co-occurrence Water Phase with 253 < Tcld < 268 K  05 Nov. 2000 - Nighttime Only

- frequency of occurrence in percent (%)

**MODIS Frequency of Co-occurrence ****Water Phase with 253 < T****cld**** < 268 K **** 05 Nov. 2000 - Nighttime Only**

##  Water phase clouds with 238K < Tc < 253K

** ****Water phase clouds with 238K < T****c**** < 253K**

## Early Estimates of UW MODIS Cloud Products Quality

- Cloud mask has demonstrated advantages of new multispectral approach;
- 	sun glint, desert, and polar problems diminished.
- MODIS cloud top pressures compare well with GOES;
- 	aircraft validation is better than 50 mb;
- 	best product from de-striped data;
- 	validated with known characteristics.
- MODIS cloud phase determinations are revealing interesting patterns;
- 	first global day/night ice/water cloud determinations;
- 	validations pending (e.g. CLEX).

## Validation committee incognito

## Early Estimates of MODIS Cloud and Atmospheric Products Quality

**IR radiances agree to within 1.5 C with GOES and ER-2 MAS/SHIS**
**Cloud mask has demonstrated advantages of new multispectral approach;**
**	****sun glint, desert, and polar problems diminished.**
**Layer tropospheric temperatures compare well with AMSU**
**	****rms better than 1 C, both within 2 C of radiosonde observations;**
**Layer dewpoint temperatures depict gradients very well**
**	****are within 2-3 C rms of radiosonde observations. **
**IR total precipitable water vapor within 3 mm rms **
**	****captures TPW gradients very well over oceans**
**	****has been improved over daytime non vegetated land**
**Ozone is very close to the GOES ozone (over North America) **
**	****rms of about 10 Dobsons **
**	****polar extreme ozone values will be improved with more training data.**
**Polar winds represent coherent atmospheric motion;**
**	****geo-like quality observed**
**	****within 7 – 10 m/s of the few raobs available for validation**
**Cloud top pressures compare well with GOES**
**	****aircraft validation is better than 50 mb**
**Cloud phase determinations are revealing interesting patterns**
**	****first global day/night ice/water cloud determinations****	**

##  Provisional Products 

**• ****Product quality may not be optimal**
**• ****Incremental product improvements are still occurring**
**• ****General research community is encouraged to participate in the QA and validation of the product, but need to be aware that product validation and QA are ongoing**
**• ****Users are urged to contact science team representatives prior to use of the data in publications**
**• ****May be replaced in archive when the validated product becomes available**
** **
_**Validated Products**_
**• ****Formally validated product, although validation is still ongoing**
**• ****Uncertainties are well defined**
**• ****Ready for use in scientific publications, and by other agencies**
**• ****There may be later improved versions**
**• ****Earlier validated versions will be deleted from the archive after a 6 month overlap period, but code for earlier versions will be maintained indefinitely**