Retrieval of aerosol optical thickness and size ... - Wiley Online Library

10 downloads 84 Views 1MB Size Report
Jan 27, 1999 - assuming an aerosol model representative of the local conditions ... MEDIS, only a limited number of the primary parameters of the aerosol .... allowed us to look at the same target at different view angles and also to ... overtly the Meridian ship on its approach to New York on July. 20. ...... 1816-1823, 1984.
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 104, NO. D2, PAGES 2261-2278, JANUARY 27, 1999

Retrieval of aerosol optical thickness and size distribution over ocean from the MODIS airborne simulator during TARFOX

D. Tanr6,1L. A. Remer, 2 Y. J. Kaufman, 2 S. Mattoo, 3 P.V. Hobbs, 4 J.M. Livingston, 5 P. B. Russell, 5 A. Smirnov 6 Abstract. Radiationand in-situmeasurements collectedduringthe TroposphericAerosol RadiativeForcingObservational Experiment(TARFOX) areusedto testthe methodfor remote sensingof aerosolpropertiesandloadingfrom theMODIS instrument. MODIS, a Moderate ResolutionImagingSpectroradiometer, will be launchedin 1999 aboardthe first EOS (Earth ObservingSystem).Followingthe MODIS procedure[Tanrd et al., 1997],the spectralradianceat the top of the atmosphere (TOA) measuredoverthe oceanin a wide spectralrange(0.55-2.13 !,tm) is usedto derivethe aerosolopticalthickness(proportionalto the aerosoltotalloading)andthe aerosolsizedistribution(integratedoverthe verticalcolumn)of the ambient(undisturbed)aerosol

bycomparing measured radiances withvalues in look-up table(LUT).TheLUT includes thegasphaseoxidationaccumulation mode,cloud-phase accumulation mode,anda coarsemodethat represents maritimeparticles(salt)anddust.In eachinversion,oneaccumulation andonecoarse modecanbe retrieved.The inversionretrievestheratioof thecontribution to the optical thicknesses of the two particlemodesandthe meanparticlesizethatbestfits the measurements. This algorithmis successfully appliedto the datasetsacquiredduringTARFOX. The MODIS _airborne simulator(MAS) aboardtheNASA ER-2 aircraftflew severaltimesduringthe experimentabovethe Universityof WashingtonC-131A researchaircrafton whichthe six-channel AmesAirborneTrackingSunPhotometer(AATS-6) wasmounted.It flew alsoabovesurfacebasedSunphotometers. Opticalthicknesses (at)• = 550 nm) aswell asthespectraldependence from the variousdatasetscomparevery well. 1. Introduction

assuming an aerosol model representative of the local

Interestin global aerosolmonitoringhas increasedrecently following the recognition that their contributions to atmosphericand Earth processesare very important. Since frequent and global-scale measurements of the Earth's atmosphereare achievableonly by observationsfrom space, efforts have been made to develop new satellite sensors and new strategiesfor retrieving key aerosolparameters[Kaufman et al., 1997]. Early studies [Griggs, 1975; Fraser, 1976; Mekler et al., 1977; Carlson, 1979; Koepke and Quenzel, 1979; Norton et

al., 1980] showedthe potentialof visible satellite imagery to

conditions [Jankowiak and Tanrd, 1992; Husar et al., 1997'

Moulin et al., 1997]. New techniquesusingthe long recordsof previous satellite data, like the UV channels of TOMS/Nimbus-7 instrument[Herman et al., 1997] and the two visible channelsof AVHRR/NOAA [Nakajima and Higurashi, 1997], are also being developed;they are expectedto better characterizethe aerosol componentthrough an estimate of their size (AngstrOmexponent)and absorptionproperties. With the new generation of satellite sensors, like MODIS/EOS [Salomonson, 1989], MISR/EOS [Diner et al., 1989], EOSP [Travis, 1993], and POLDER/ADEOS

[Deschampset al., 1994], more parameterscan be derived. Correspondinginversionschemesuse the spectraldependence channels on the initial instruments, one channel for GOES and of the multiple channelsatellite radiances[Tanrdet al., 1996; Meteosat and two channels for AVHRR/NOAA, these Stowe et al., 1997], the multiangle capability [Martonchick algorithms can derive only the total aerosol loading, and Diner, 1992], the polarized ratio of the reflected light [Mishchenko and Travis, 1997] or both [Deschampset al., •Laboratoire d'OptiqueAtmosph6rique,Villeneuved'Ascq, 1994]. Although more sophisticated,all of these new methods FRANCE. still require assumptionsand have their own limitations. 2Laboratory forAtmospheres, NASAGoddard Space FlightCenter, This paper focuseson the multi-wavelength approach to Greenbelt,Maryland. aerosol retrieval from remote sensing. The algorithm 3Science Application Corporation, Vienna, Virginia. developed for the Moderate Resolution Imaging 4University ofWashington, Seattle. Spectroradiometer(MODIS)is applied to the data acquiredby 5NASAAmesResearch Center,MoffettField,California. the M__ODISairborne simulator (MAS) [King et al., 1996]

detect aerosols over ocean. Because of the limited

number of

6Science Systems andApplications Inc.,Laboratory for Terrestrial aboard the NASA ER-2 aircraft during the Tropospheric

Physics,NASA GoddardSpaceFlight Center,Greenbelt,Maryland.

Copyright1999by the AmericanGeophysicalUnion.

Aerosol

Radiative

(TARFOX)

Forcing

Observational Experiment

[Russell et al., this issue (a)]. The TARFOX

Papernumber 1998JD200077.

campaign was conductedJuly 10-31, 1996, off the United States mid-Atlantic coast and was designed to reduce

0148-0227/99/1998JD200077509.00.

uncertainties 2261

in the aerosol

radiative

effects.

The retrieved

2262

TANR• ET AL.'RETRIEVALOFAEROSOL PROPERTIES OVERTHEOCEAN

aerosolpropertiesare testedagainst groundand airborne measurements providedby the Aerosol Robotic Network (AEReNE2)field networkof Sun-skyradiometers [Holbenet al., 1998; Rerneret al., this issue],andby the in situ and

the coarsemode(radius> 1.0 gm), the optical thickness, and

et al., this issue (b)].

measurements.

2. Descriptionof the Algorithm

combinationof two lognormal laws, which correspond to the accumulation(small) and coarse (large) modes.The physical

the mean particle size. A complete descriptionof the algorithmis given by Tanrdet al. [1997]. The algorithmis basedon a comparisonbetweenthe measurements and the precomputedradiances;the solution correspondsto the radiometric measurementsobtained aboard the University of Washington C-131Aaircraft[Hobbset al., this issue;Russell aerosol model within a look-up table that best fits the We assume that the aerosol size distribution

follows

a

TheMEDIS algorithm 'tested inthispaper uses thesolar- propertiesof the aerosolmodel(i.e., the particlerefractive reflected spectral radiances, normalized to a central index, the moderadius,and the width of the mode) are given in wavelength to derivethe aerosolsizedistribution; it thenuses Table 1. These characteristics have been derived from in situ the retrieved aerosol size distribution and the corresponding measurementsand are representative of several different phasefunctionfor an accurate derivation of the aerosoloptical aerosolregimes.Six valuesof the aerosolloading,expressed by theopticalthickness 'ca at •k= 550 nm are considered in the thickness.Even with the large spectralcoverageprovidedby MEDIS, only a limitednumberof the primaryparameters of look-uptable (LUT) ('ca= 0.0, 0.2, 0.5, 1.0, 2.0, and3.0). To the aerosol size distribution can be retrieved; the rest of the

describethe surface,we use the rough ocean model proposed

parameters haveto be assumed. A sensitivity study[Tanrdet al., 1997] showedthat the mainaerosolparameters, whichcanbe derived,arethe optical

by Coxand Munk [ 1954];for the Fresnelreflectionon the sea surface,threetypical valuesof wind speedareused:2.0, 6.0, and10.0 m s'•. The percentage of the seacoveredby foam followsKoepke's[1984] model,and the spectralreflectanceof

thicknessratio of the accumulationmode (radius< 0.5 gm) and

Table1. Description of Aerosol Models Usedin theLook-Up Table:Values of Real(mr)andImaginary (mi)Partsof the

Spectral Refractive Index, Median Radius (rg),Standard Deviation ({J),andEffective Particle Radius (reff) oftheLognormal Size Distributions

mr - i mi

)• = 470- > 860 nm

mr - i mi

)• = 1240 nm

mr - i mi

mr - i mi

)•=1640 nm

)•=2130 nm

rg

•J

reff

Comment

0.035

0.40

0.05

watersoluble

Smallparticles

1

1.45-0.0035i

1.45-0.0035i

1.43-0.01i

1.40-0.005i

type

2

1.45-0.0035i

1.45-0.0035i

1.43-0.01i

1.40-0.005i

0.07

0.40

0.10

water soluble type

3

1.45-0.0035i

1.45-0.0035i 1.43-0.01i

1.40-0.005i

0.06

0.60

0.15

water•oluble type

4

1.40-0.0020i

1.40-0.0020i

1.39-0.005i

1.36-0.003i

0.08

0.60

0.20

water soluble with humidity

5

1.40-0.0020i

1.40-0.0020i

1.39-0.005i

1.36-0.003i

0.10

0.60

0.25

water soluble

0.40

0.60

0.98

wet sea-salt

with humidity

Largeparticles

6

1.45-0.0035i

1.45-0.0035i

1.43-0.0035i

1.43-0.0035i

type

7

1.45-0.0035i

1.45-0.0035i

1.43-0.0035i

1.43-0.0035i

0.60

0.60

1.48

wet sea-salt type

8

1.45-0.0035i

1.45-0.0035i

1.43-0.0035i

1.43-0.0035i

0.80

0.60

1.98

wet sea-salt type

9

1.53-0.008i

1.46-0.008i

1.40-0.008i

1.22-0.009i

0.40

0.60

0.98

dust-liketype

10

1.53-0.008i

1.46-0.008i

1.40-0.008i

1.22-0.009i

0.50

0.80

2.50

dust-liketype

11

1.53-0.008i

1.46-0.008i

1.40-0.008i

1.22-0.009i

1.00

0.80

4.95

dust-liketype

Smallparticles correspond totheaccumulation (small) modeandlargeparticles tothecoarse (large) mode.

TANRI• ET AL.: RETRIEVALOF AEROSOLPROPERTIES OVERTHE OCEAN the foam is taken from Whitlock

et al. [1982].

Recent

measurements [Frouin et al., 1996] lead to a different spectral behavior

of the foam reflectance.

Since the use of these new

valueswouldrequireadjustments in the LUT, these adjustments are not done here, but we recognizethat this omission may result in additional uncertainties. Finally, the water-leaving radianceis 0.0 in all channelsexceptat X,= 555 nm where a constant value of 0.005 (expressedin reflectance units) is assumed[Gordon, 1997]. This assumptionis questionablein coastal regions. The inversion

scheme uses six channels between X, = 555

nm and 2130 nm. The blue MODIS channel (470 nm) is not consideredbecauseof the uncertainties of the water-leaving radiance which result from fluctuations in the chlorophyll content. We apply a version of the operational mask [Ackerman et al., 1998] which is based on threshold values in the visible and thermal infrared channels. In this study, we do not correct for ozone contamination

in channels centered at X,

= 555 and 659 nm but will do so in the satellite

version

of the

algorithm. This point is discussed in section 6. The water vapor(channelscenteredat X,= 865 and2130 nm) absorption

2263

of the track were separatedby 10 km. Figure 1 shows the flight patternon July 25, 1996. Looping aroundthe racetrack allowed us to look at the same target at different view angles and

also

to

monitor

the

time

evolution

of

the

derived

parameters.

The ER-2 was directedto fly either in the solar plane or in the plane perpendicularto the solar plane to avoid Sun glint. For low Sun elevations the ER-2 flew in the solar plane, scanningside to side since there was no probability of glint along the scan. For high Sun elevations the ER-2 flew in the perpendicularplane to obtain glint free measurementsin the backscatteredpart of the scan. The ER-2 was also directed to overtly the Meridian ship on its approachto New York on July 20. On July 26 it overflew groundinstrumentationon Bermuda Island as well as on the Meridian ship as it left this island. The ER-2 flew on nine daysin TARFOX (Table 2). Data from six flights are analyzedhere; they correspondto a large range of atmosphericconditions,from very clean on July 20 and 21 to very hazy on July 17 and 24 with intermediate conditions on July 25 and 26.

featuresare weak enoughto be neglected.Sincethe time period Validation Data Set of TARFOX was free of major volcanic effects, no 3.2. stratosphericcorrectionwas applied. To avoid measurements To validate the MAS optical thickness results, the sixcontaminatedby Sun glint, we rejecteddata acquiredwithin a channel Ames Airborne Tracking Sun Photometer (AATS-6) coneof 30ø half angle around the speculardirection. Since the aboardthe University of WashingtonC- 131A aircraft was used wind speedvalue was not available during the experiment, a [Matsumoto et al., 1987; Russell et al., this issue (b)]. AATS-

singlevalue(6.0 m s'l) is usedin thealgorithm. Theimpactof

this assumptionon the retrieval is also discussedin section6. The retrievedaerosolparameterswe examine are the optical

thickness('Ca)at )• = 550 nm, the effectiveradius(reft) of the size distribution, and the ratio rl between the two modes. The rl corresponds to ('cs/('cs+ 'c•))where'csand 'c• are respectively, the optical thicknessof the small (accumulation)and the large (coarse) modes at 550 nm. The general accuracy of the retrieved parameterswas estimatedin a previous paper to be

A'ca = _+0.05+_0.05'c a (at X, = 550 nm), Areff = 0.3reft, and Arl = _-/-0.25[Tanrd et al., 1997].

3. 3.1.

Data MAS/ER-2

Data

The MODIS airborne simulator (MAS) flew on the NASA ER-2 aircraft in TARFOX. The instrumentis describedby King et al. [1996]. It is a scanningspectrometerwith 50 channels

coveringthe rangeX,= 550 nm to 14,300 nm, but only the five channelscentered at X,= 549, 655, 867, 1643, and 2105 nm, equivalent to the MODIS/EOS channels [Salomonson et al., 1989], are used. The instrument scans in the across-track direction.The half scanangle is 42.96ø, which correspondsto

6 measured the solar direct-beam

transmission

in four aerosol

bands centered at X,= 380, 451, 525, and 1021 nm. Aerosol optical depthsand uncertaintieswere derivedfrom the AATS-6 measurements,as describedby Russell et al. [1993]. These data allow validation of the optical thickness derived at 550 nm from MAS and also the spectral dependencerelated to the aerosolsize distribution.To compare the MAS inversion with the C-131A measurements,the flight of the ER-2 was timed to be above the C-131A when it was at its lowest altitude, just above sea level. In this way the optical thicknessmeasuredby the AATS-6 looking upwardrepresentsthe sametotal column value retrievedfrom the MAS instrumentlooking downward. The retrievedoptical thicknesseswere also comparedwith surface-basedSun Photometer measurements[Remer et al., this issue]. Instruments were put on two ships cruising between New York City and the BermudaIslands. We also used a Sun Photometer of the AERONET program [Holben et al., 1998] locatedat the biological stationin the BermudaIslands. The MAS inversions represent the whole atmospheric column. The optical thicknessesmeasuredby Sun Photometers

on the ground or from a low-flying aircraft are also representative of the whole column. The in situ aerosol size distribution measurementsmade with the p_.assivecavity a_erosolspectrometer p_.robe(PCASP) [Hobbs et al., 1998] a swath width of 36 km at the nominal altitude of 20 km. The spatial resolutionis 50 m at nadir, and pixels are aggregated aboardthe C-131A aircraft are representativeof one or more into boxes of 40 x 40 pixels in order to derive the aerosol atmosphericlayers.The PCASP has 15 bins covering the size parameters on a scaleof 2 x 2 km2. Thecalibration of the range 0.10 to 3.00gm diameter.The aerosol size distribution shortwavechannels(from X,= 550 to 2500 nm) is performed derivedfrom the PCASP has not been adjustedto the ambient relative humidity, which may produce differences with the on the ground using an integrating sphere. There is no onboard calibration calibration

is tested

for the shortwave channels, but the

actual aerosol

in

distributionis performedusingthe effective radiusreft defined

a thermal

vacuum

chamber

and

a

temperaturecorrection applied during data processing; the accuracyof the calibration is approximately5%. The ER-2 usually flew along a "racetrack"patternover open water, approximately160 km in length in which the two legs

by

size distribution.

reft =

The

validation

I•0 r3n(r)dr 0•r2n(r)dr

of

the

size

(1)

2264

TANRI• ET AL.: RETRIEVAL OF AEROSOLPROPERTIESOVER THE OCEAN

July 17, 1996, flight 10 Tau

Radius

Ratio

July 17, 1996, flight 11 0.8

0.6

0.6

0.5

1.0

:.

-....... 0.9

....

0.4

0.8

0.3

0.7

0.2

0.6

0.1

0.5

0.4

0.2

0.0

July 17, 1996, flight 12

II

:-am

Plate 1. Maps of aerosol parametersderivedfrom the MAS measurementson July 17. Optical thickness (left), effective radius (middle) and ratio rl (right). Informationon tracks 10, 11, and 12 is given in Table 3.

TANR• ET AL.: RETRIEVAL OF AEROSOLPROPERTIESOVER THE OCEAN

I• I,.I,•,HT g!!. ] $3

25 ,,IIJLY ] EJq5

n/ [' '7.ca

PIFIS-5 O / ¾J:3 f LR5E

Figure 1. Exampleof an ER-2 flight patternduringthe TARFOX experiment.This pattern was performedon July 25, 1996.The distancebetweenthe two legsis around10 km, andthe legsare 160 km long. Table 2. Summaryof MAS Data CollectedDuring TARFOX Flight

Time,

Location

Comment

Atmospheric

Number

1996

UTC

Conditions

96-146

July14

1707-1825

mid-Atlantic states

fightaborted

/

96-147

July16

1709-1826

Delawarepeninsula/

testflight

/

Atlantic Ocean

96-148

July17

1712-2130

Delawarepeninsula/ Atlantic Ocean

96-149

July20

1723-2230

Delawarepeninsula/ Atlantic Ocean

96-150

July21

1720-2047

Delawarepeninsula/ Atlantic Ocean

96-151

July22

1759-1936

mid-Atlantic coast/ Atlantic

96-152

July24

1710-1949

July25

1508-1815

mid-Atlantic coast/ mid-Atlantic coast/ Atlantic Ocean

96-154

July26

1809-2305

AtlanticOcean/ Bermuda

coordination

C-131A/ship coordination

C-130/ship coordination

testflight

cloudy,somecloudless holes,hazy very clean conditions, low humidity very clean conditions, low humidity

/

Ocean

Atlantic Ocean

96-153

C-131A

C-131A

very hazy and cloudy

coordination

C-131A,C-130

cloud free and hazy

coordination

ship coordination

cloudy near the coast and clear later

2265

2266

TANRl•ET AL.' RETRIEVALOFAEROSOLPROPERTIES OVERTHEOCEAN

4. Results

4.1.

July

17, 1996

For each selectedflight, maps of the following aerosol The meteorologicalconditions on July 17 were not ideal productsare provided:(1) the optical thicknessXa at )• = 550 sincelow-level "fog-likeclouds"weremoving from the shore nm; (2) the effective radius reft of the complete size eastwardandhigh-level cirruswerealso moving eastwardand expanding in extent. Nevertheless a hole within the cloud cover developedandprovidedabout 1 hourof clear-skydata, coincident measurements were obtained from the two aircraft starting at around1830 UTC, which correspondedto ER-2 or from the ER-2 and the surface. Summaries of the ER-2 flight tracks 10, 11, and 12 (Plate 1). At that time, the C-131A flightsare givenin Table 3. We will comparethe MAS optical wasat its lowest altitude(100 feet) at 37.44øN/74.09øW,just thicknessderivedat )• = 550 nm and its spectraldependence below the ER-2 on track 10 when there was Sun glint in the computed from the retrieved aerosol model with the Sun middleof thescan.A few minuteslaterat 1900 UTC the image Photometermeasurements. The effectiveparticleradiusderived was almostglint free (track 11), but the C-131A was then at an from the MAS is comparedwith the values obtained from the altitudeof 3750 m. The ER-2 passed againoverthe first leg of in situ measurements aboard the C-131A aircraft over a vertical the track pattern(track 12). profile usingequation(1). In only two flights did the C-131A In Plate 1, we definethe glint mask to be +30ø aroundthe complete a vertical profile coincident with a MAS retrieval. Fresnelreflection.This thresholdvalue shouldbe increased up Thus comparisons of reft are only madetwice. to 35ø sincethe image was obviouslystill contaminated(i.e., distribution; (3) the ratio q. We analyzed those parts of the flight tracks where

Table 3. Characteristicsof ER-2 Flight Tracks That Were Selectedand Analyzed.

Flight Track

Time, Lat(ø)/Long(ø) Hour Minute Second

SolarZenithand AzimuthAngles(ø)

Heading (ø)

July 17, 1996 10

1833:27 1849:24

37.46øN/74.19øW 38.53øN/72.47øW

{25.6o;236.7 ø} {30.00;242.8 ø}

48.6 ø

11

1852:38 1906:59

38.64øN/72.53øW 37.55øN/74.17øW

{30.50;242.9 ø} {31.4ø;246.9 ø}

230.5 ø

12

1910:45 1926:33

37.43øN/74.24øW 38.50øN/72.52øW

{32.1ø;248.6 ø} {36.70;252.7 ø}

49.7 ø

July 20, 1996 04

1818:42 1834:06

36.88øN/73.32øW 38.44øN/74.33øW

{23.90;232.0 ø} {26.60;233.9 ø}

333.6 ø

15

2109:14 2122:36

37.18øN/72.54øW 37.93øN/70.91øW

{56.8ø;271.7 ø} {60.70;274.2 ø}

59.6 ø

July 24, 1996 08

1857:20 1908:08

38.39øN/73.65øW 38.98øN/75.02øW

{31.4ø;241.5 ø} { 32.60;242.3 ø}

296.7 ø

09

1911:30 1920:40

39.21 øN/74.89øW 38.70øN/73.73øW

{33.2ø;243.2 ø} {35.40;247.7 ø}

119.4 ø

03

1544:55

38.04øN/74.01

1555:03

37.07øN/74.76øW

{23.5ø;133.8ø}

38.13øN/73.95øW 37.08øN/74.75øW

{21.3ø;211.8 ø} { 21.4ø;217.8 ø}

210.3 ø

36.91 øN/74.68øW 37.88øN/73.93øW

{21.9ø;221.8 ø} {24.40;226.5ø}

49.7 ø

July 25, 1996

13

1806:00 1815:56

øW

July 26, 1996

04

2102:40 2116:16

33.13øN/65.93øW 32.26øN/64.52øW

{61.6ø;275.4ø} {65.70;278.2ø}

124.8ø

The universaltime, the latitudeandkongitude,andthe solarzenithandazimuthanglesare givenfor the centerof the scan.Valuesat the beginning(top) andat the end (bottom)of each trackare provided.The headingis the ER-2 azimuth.

TANRl• ET AL.: RETRIEVAL OF AEROSOLPROPERTIESOVER THE OCEAN 1.00

,

0.30 17 July 96

0.01

MAS/Track MAS/Track MAS/Track

ß

AATS-6

10 11 12

........

0.2



0.20



0.15



0.10

'-,-,

0.05

0.00 0.4

0.6

0.8

1.0



I

I

3.0

0.0

Figure 2. Aerosol optical thickness measuredby the Sun photometer(AATS-6) aboardthe C- 131A aircraft and derived from MAS data as a functionof wavelength.The comparisonis madeover a point located at 37.49øN/74.16øW around 1830 UTC. Information on tracks 10, 11, and 12 is given in Table 3.

the retrievedoptical thicknessdisplays an unexpectedpattern in the middle of the scan). A strong, shallow fog layer was identified by the AATS-6 aboard the C-131A [see Russell et al., this issue (a), Plate 3d] which can explain the larger optical thicknessobservedin the southernpart of the images. The _LidarAtmospheric Sensing _Experiment(LASE) airborne lidar Ilsmail et al., 1997] aboard the ER-2 [see Russell et al., this issue (a), Plate 3c] shows that the layer was not uniform,

but there is no informationon the spatialextent of the layer to our results.

Figure 2 showsthe optical thicknessesmeasuredby the Sun photometeron the C-131A and those derivedfrom the MAS data as a function of wavelength.The comparisonis madeover

a pointlocatedat 37.49øN/74.16øWwhenthe ER-2 wasflying

I

17 July 96 18:30-19:00 N37.49-W74.16

• 0.5

• 1.0

1.5

Wavelength(•tm)

confirm

I

• 0.25

0.10

o [] /x



2267

2.0

2.5

Altitude

(km)

3.0

3.5

4.0

Figure 3. Effective particle radiusmeasured by the Passive Cavity Aerosol SpectrometerProbe (PCASP) aboard the C131A aircraft as a function of altitude.

In Figure4, we showthe MAS opticalthicknessaveragedin the cloud free region correspondingto a geographical area defined by points located at 38.17øN/74.0øW and 38.27øN/74.10øW; the time of the acquisition was around 1830 UTC. The time correspondsto a portion of the ER-2 track 4. The comparisonis madewith the Sun photometer measurements aboard the C- 131A acquired over 38.20øN/73.98øWat 1842 UTC. The agreementis still good (A'C= 0.015); the retrievedvalueof 'ca at )• = 550 nm is 0.036 when the Sun photometer measured around 0.05. The agreementbetweenspectraldependence is not so good since the MAS gave an almost constant value, while the Sun photometer data show a slight decrease as a function of wavelength. Becauseof the extremely clean conditions, the relative accuracyof the inversion is decreasedsignificantly.

Nevertheless,the algorithm gives the optical thickness'ca within a +0.015 accuracy.

above the C-131A. The agreement is very good (i.e., the

optical thicknessat )• = 550 nm is very close to the Sun photometer value). The derived spectral dependencealso matches the measurementvery well, which means that the aerosol model derived from the ER-2 remote sensing measurementshas the right optical properties. The aerosoleffective radiusmeasuredby the PCASP during its profile is plotted as a function of altitude in Figure 3. Values are in the range 0.10 to 0.15 •tm. The median values computedfrom the ER-2 measurements over the uniform part of the aerosollayer for tracks10, 11, and 12 are 0.219, 0.184, and 0.183 gm, respectively. Mdszaros [1970] showedthat an ammoniumsulfate particle with a dry radiusof 0.14gm can increaseits size by a factor of 2 at a relative humidity of 80%. It is excludedto state that similar effects occurredduring the TARFOX experiment, but correcting the aerosol size distribution for relative humidity, which is not a trivial matter, would definitely producelarger sizes. The decreasein effective radius observed for tracks 11 and 12 may actually occur, but there are no in situ measurements to confirm it.

4.2.

July

20,

1996

July 20 was a very clean day with low humidity, so we expectedvery low aerosol optical thickness. The ER-2 flew perpendicularto the solar plane.

1.00

i

i

i

i

o

MAS/Track

ß

AATS-6

i

i

i

]

04

20 July 96

O.lO

:

:

I

O.Ol 0.2

0.4

i

i

0.6

i

i

0.8

i

I

1.0

i

3.0

Wavelength (•tm) Figure 4. Aerosol optical thickness measuredby the Sun photometer (AATS-6) aboardthe C- 131A aircraft and derived from MAS data as a function of the wavelength. The Sun

photometer data were acquired over a point at 38.20øN/73.98øW at 1842 UTC. The optical thicknesses derivedfrom the MAS are averagedover a geographical area defined by points located at 38.17øN/74.0øW and 38.27øN/74.1øW;the time of the acquisitionwas around1830 UTC. Informationon track 4 is given in Table 3.

2268

TANRI•ET AL.:RETRIEVALOFAEROSOLPROPERTIES OVERTHEOCEAN 1.00

i

]

i

!

i

[

[

cloudsto theobservation area,andbandsof cirruswerequickly movingover the oceanfrom the continent.We analyzedER-2 tracks8 and9, whichcorrespond to thebestpartsof the flight consideringthe occurrence of cloudsandthe presenceof Sun glint duringthe first loops. Track 8 corresponds to the time whenthe C-131Aaircraftmadeits verticalprofile, although the exactlocationwaswithin the part of the MAS image

[

20 July 96

0.10

affectedby the Sun glint.

Plate 3 showsmapsof the derivedaerosolpropertiesfor both tracks8 and9. The imageis affectedby Sunglint, but thereis sufficientunaffected oceansurfacein the scanaway

SHIP

MAS/Track 0.01

I

i

0.2

I

15

I

0.4

I

0.6

I

0.8

I

fromtheSun. Notethatourcloudmaskdoesnotrejectthetiny

[

i

cloudsin the top left part of the images,but this will not be

1.0

3.0

Wavelength(gm)

the casewith the operational MODIS cloudmask. An aerosol

plumecanbe seenin the centerof both ER-2 flight tracks,

whichmeansthat the plume persistedin the areafor at least a Figure 5. Aerosol optical thicknessmeasuredby the CIMEL half hour.The largeropticalthicknessis connectedto smaller instrumentaboardthe Meridian ship and derivedfrom the MAS data as a function of wavelength. The ship data are averaged particles(as seenin the ratio B). The aerosol optical thickness measuredby the Sun over a geographical area defined by points located at 37.44øN/70.28øW

and 38.02øN/70.50øW;

the time of the

acquisition was between 2011 and 2124 UTC. The optical thicknesses derived from the MAS are averaged over a geographical area defined by points located at 37.40øN/71.86øW and 37.50øN/71.56øW; the time of the acquisitionwas around 2120 UTC. Information on track 15 i s given in Table 3.

The ER-2 then flew over the cruiseshipMeridian (track 15). In Figure 5, the aerosol optical thicknesses measuredby the CIMEL instrumentaboardthe Meridian ship and derivedfrom MAS dataare plottedas a functionof the wavelength.The ship data are averagedover a geographical area defined by points located at 37.44øN/70.28øW

and 38.02øN/70.50øW;

the time

of the acquisition was between 2011 and 2124 UTC. The optical thicknessesderivedfrom MAS are averagedover a cloud free geographical area defined by points located at 37.40øN/71.86øW

and 37.50øN/71.56øW;

compared in Figure6. The Sunphotometerdatawereacquired over point 38.80øN/74.33øWat 1900 UTC. The optical thicknessesderived from MAS are averagedover two geographical areas defined by points located at 38.90øN/74.47øW for zone a and at 38.82øN/74.52øW for zone

b; zone a is within the aerosolplume, and zone b is on the

edgeof theplumecloseto the endof the scan.Theagreement is perfectin zonea but the spectraldependence is differentfor zoneb. AlthoughtheSunphotometer datawereacquiredwithin the glint, we can reasonablyassumethat the AATS-6 was samplingthe aerosolplume observedin zone a. 4.4.

July

to be nonuniform.

These

inversion

variations

scheme

and/or

could

result

from

in the calibration

artifacts

in

1.00



did not confirm

was around0.15gm.

this

of the instrument.

since the effective

radius

24 July 96

0.10

the

It can also be attributedto the higher humidity in the presence of cloudsat the beginningand at the end of the track.The large valuesof reft are consistent with the weak spectraldependence observed in the Sun photometer data but the PCASP measurements

1996

The conditions on this day were not perfect because althoughthe sky was clearat the beginning of the flight, the

Maps of the aerosol parametersderived from track 15 are shown in Plate 2. The optical thickness shows an angular dependenceand varies between 0.045 and 0.095 along the scan.

25,

the time of the

acquisition was around 2118 UTC. The agreement is better thanfor the case shown in Figure 4. In this case, the spectral dependenceis better retrieved, although the optical thickness derivedat 550 nm is again smaller than the Sun photometer measurements. Note that the optical thicknessis twice as great as in Figure 4, which can explain the better agreement.The aerosollayer was also observedby the Sun photometeron the C-131A

photometer on the C-131A and derived from MAS data are

ß

AATS-6

[] ¸

MAS/Track MAS/Track

I

0.01

0.2

I

0.4

I

I

0.6

08/b 08/a I

I

I

I

0.8 1.0

I

3.0

Wavelength(gm) Figure 6. Aerosol optical thicknessmeasuredby the Sun photometer (AATS-6) aboard the C- 131A aircraft and derived from MODIS airborne simulator (MAS) data as a function of

wavelength. The Sun photometer data were acquiredat 38.80øN/74.33øW at 1900 UTC. The optical thicknesses derivedfrom the MAS are averagedover two geographical On this day the meteorologicalsituationwas againnot ideal areasdefinedby pointslocatedat 38.90øN/74.47øWfor zone a for deriving aerosolparametersfrom the MAS. A stationary andat 38.82øN/74.52øW for zoneb; the timeof the acquisition front, orientatedeast-westand crossing the shoreline over was around1904 UTC. Information on track 8 is given in North Carolina, was slowly moving northward,bringing Table 3. 4.3.

July

24,

1996

TANRI•ET AL.:RETRIEVALOFAEROSOLPROPERTIES OVERTHE OCEAN

2269

July 20, 1996, flight 15 Tau

Radius

!

Ratio

g

0.20

1.0

1.0

0.9

0.8

0.8

0.7

0.6 0.6 0.10

0.5 0.4 0.4

0.05 0.3

0.2

0.2

o.oo .

:,

4•I' 11

r lB

sa

.

II

,

0.1

.""

gl

ß

!

I .d

II

Plate 2. Sameas Plate 1 but for July 20. Informationon track 15 is given in Table 3.

0.0

2270

TANRI• ET AL.:RETRIEVALOF AEROSOLPROPERTIES OVERTHE OCEAN

July 24, 1996, flight 8 Tau

Radius

0.8

Ratio

0.6

1.0

0.5

0.9

0.4

0.8

0.3

0.7

0.2

0.6

0.1

O.5

0.6

0.4

0.2

0.0

July 24, 1996, flight 9

0.8

0.6

1.0

0.5

0.9

0.4

0.8

0.3

0.7

0.2

0.6

0.1

0.5

0.6

0.4

0.2

0.0

Plate 3. Same as Plate 1 but for July 24. Informationon tracks8 and 9 is given in Table 3.

TANR!• ET AL.: RETRIEVAL OF AEROSOL PROPERTIESOVER THE OCEAN

2271

July 25, 1996, flight 3 Tau

Radius

0.8

Ratio

0.6

1.0

0.5

0.9

0.4

0.8

0.3

0.7

0.2

0.6

o. 1

0.5

0.6

0.4

0.2

o.o

July 25, 1996, flight 13

0.8

0.6

1.0

0.5

0.9

0.6

0.4

'

0.8

0.4

0.2

o.o

0.3

0.7

0.2

0.6

o. 1

0.5

Plate 4. Same as Plate 1 but for July 25. Informationon tracks3 and 13 is given in Table 3.

2272

TANRI• ET AL.: RETRIEVALOF AEROSOLPROPERTIES OVERTHE OCEAN

July 26, 1996, flight 4 Tau

Radius

Ratio

i

0.8

.

"

0.6

1.0

0.5

0.9

0.4

0.8

0.3

0.7

I m

i

i

i

0.6

0.4

0.2

. 0.2

I .

ß

ß

0.0

I I

I

I

.

I

0.6

.

0.]

0.5

I I ii

I

i

i

m

m

m

t

m"

I'

Plate 5. Same as Plate 1 but for July 26. Information on track 4 is given in Table 3.

TANRI•ET AL.:RETRIEVALOFAEROSOLPROPERTIES OVERTHEOCEAN

MASimages wereaffected bySunglint(track3). Laterin the afternoon,whenthe geometrical conditionswere more

2273

0.50

favorable for avoiding the Sun glint, there was more

cloudiness. Nevertheless, holesin the cloudcovercan be

25 July 96

0.40

foundon tracks 12 and 13 that correspondto coincident Sun

photometer measurements performed aboard theC-131Awhen

it was flying atitslowest altitude.

Oøo o

0.30

c• o oo o

oo

Plate 4 showsthe aerosolpropertiesfor both tracks3 and

13. Track3 wasstronglyaffected by Sunglint, andlarge

0.20

instabilities in the aerosolretrievalcanbe seenon the edgeof the glint, mainly in the valuesof the ratio fl. An aerosol

0.10

plume canbeseenin thebottom portion of theplate.Track

o

AATS-6

ß ß

MAS/Track MAS/Track

13 12

13 is moreinterestingbecausebothedgesof the scanare glint

free.Again,theaerosol layerisnotuniform asseen alsoin the

0.00

I 37.6

Sun photometermeasurements.

37.7





37.8

37.9

• 38.0

• 38.1

I 38.2

38.3

In Figure7, the aerosoloptical thicknessmeasured by the Latitude (ø) Sunphotometer aboardtheC- 131Aandderivedfrom MAS data are plotted as a functionof wavelength.The Sun photometer Figure 8. Aerosol optical thickness measuredby the Sun data were acquiredat latitude 37.84øN along the transect photometer(AATS-6) aboardthe C- 13IA aircraft and derived from the MAS data as a function of latitude. The Sun performedat 74.15øW around 1849 UTC. The optical photometer data were acquired along the transect performed at thicknessesderivedfrom the MAS were for a point locatedat 74.15øW between 1845 and 1857 UTC. Values were

37.85øN/74.15øW; the time of the acquisition was around 1815 UTC. The agreement is again very good, the small discrepanciescan be attributedto different acquisitiontimes. Aerosol optical thicknessesmeasured by the Sun photometer

interpolatedat )•= 550 nm. The optical thicknessesderived from the MAS were obtained over the same meridian 74.15øW at 1750 UTC for track 12 and at 1815 UTC for track 13.

along the transectare plotted as a functionof the latitude in Figure8. Values have been interpolatedat 3•=550 nm for comparisonwith the MAS retrievals. Two points resulting from the MAS inversionhave been plotted' one for a pixel of

140



track 12 located on the transect and the other for track 13. The 120

-

100

-

imageshowsvaluesbetween0.25 and 0.28gm (Figure9a),

80

-

which is twice larger than the PCASP measurementsalong its profile (Figure 9b). Again, the PCASP measurements are not

60

-

40

-

agreementis perfect. Finally, the histogram of the effective particle radius derivedfrom the MAS for track 13 over the left part of the

adjusted totheambient relative humidity.

20 0

1.00

0.00

25 July 96

0.05

0.10

0.15

0.20

750

0.25 I

0.30 I

0.35

0.40

I

800

25 July 96 0.10

850

0.01

'

0.2

o

MAS/Track

ß

AATS-6

i

i

0.4

i

i

0.6

900

13

i

i

i

[

0.8 1.0

'

-

1000

-

ß.

o

o

3.0 1050

Wavelength(gm) Figure

950

0.00

0.05

0.10

0.15

0.20

0.25

0.30

0.35

0.40

Effective Radius (gm)

7. Aerosol optical thickness measuredby the Sun

photometer(AATS-6) aboardthe C- 13IA aircraft and derived from the MAS data as a function of wavelength. The Sun

Figure 9. (a) Histogram of the effective particle radius derivedfrom the MAS for track 13 over the left part of the image (see Plate 4). (b) Effective particle radius measuredby

photometerdata were acquiredat 1849 UTC over a point at latitude 37.84øN on the transect performed at 74.15øW. The optical thicknessesderivedfrom the MAS were obtained at 37.85øN/74.15øW; the time of the acquisition was around

during a vertical profile performedbetween 1845 and 1930

1815 UTC.

UTC.

the PCASP aboard the C-131A

aircraft at 38,64øN/74.15øW

2274

TANR• ET AL.: RETRIEVAL OF AEROSOLPROPERTIESOVER THE OCEAN 1.00

,

,

'

'

'

'

'

I

i

:

0.50

. 0.40

-

26 July

ß _

Bermuda Island 0.30

ß 0.10

-

..........................

ß

.le 0.20

::

ß e e.e ee

ß

eel

-

MAS/Track 04 SHIP

o

0.10

26 July96

ß

(a)

0.01

i

i

0.2

i

0.4

i

i

i

0.6

i

J

0.8 1.0

3.0

0.00 8:00:00

12:00:00

16:00:00

20:00:00

1.00

Time (GMT)

• ••

Figure 10. Aerosol optical thickness measuredat the biologicalstationin BermudaIsland as a functionof the time.

26 July 96:

0.10

4.5.

July

26,

1996

This flight can be dividedinto two parts.The ER-2 first flew a "racetrack"pattern over the target area in order to obtain coincident

measurements with

the C-131A.

However,

cloudiness

was variable

with

broken

clouds.

The

aerosol

optical thickness varied between 0.15 and 0.50 along the transect;a similar increase was observedby the ground-based CIMEL Sun photometer during the day at the biological stationin Bermuda (Figure 10). In Figure 11 we comparethe MAS-retrieved values to the Sun photometer data from the Meridian ship as it left the BermudaIslands(Figure 1la) and to the CIMEL measurements from the biological stationnear the time of the ER-2 overpass time (Figure 1lb). Again, the agreementis very good.

MAS/Track

o /•

Biological Station 20:56 Biological Station 21:06

¸

the

flight did not go as plannedbecauseof cloudiness.The ER-2 then flew southeastto overtly the Bermuda Islands. Plate 5 illustrates the results. During the flight the

ß

04

Biological Station 21:17 i

0.01

i

0.2

i

0.4

i

i

i

0.6

0o)

i J

0.8

i

1.0

3.0

Wavelength(pm) Figure 11. (a) Aerosol optical thickness measuredby the CIMEL instrument aboardthe ship and derivedfrom the MAS instrument. The ship data correspondto 32.23øN/64.34øW; the acquisition time was around2145 [ITC. The MAS optical thicknesses are averagedover a zone from 32.34øN/64.83øW to 32.11 øN/64.78øW; the time was around 2115 UTC. (b) Aerosol optical thicknessesmeasuredat the biological station in Bermuda

Island at three different

times.

5. Summary of Results To illustrate the performancesof the algorithm, we will summarize the results for those days when coincident verification measurements were available. The optical

thicknesses Xa derivedat t = 550 nm from theMAS (Figure 12) comparewith the "ground-truth"measurementswith a better accuracythan expected.The error can be describedby Axa = +0.01+0.05Xa (instead of the theoretical value of Axa = +0.05+0.05xa [Tanr6 et al., 1997] for the aerosol type encounteredduring the experiment. Concerning the effective particle radius,comparisonson July 17 (reft = 0.18 to 0.22 pm from the MAS, and reff = 0.10 to 0.15 pm from the PCASP) andon July 25 (reft = 0.25 to 0.28 !.tmfrom the MAS, and reft = 0.10 to 0.15 pm from the PCASP)give an accuracy

lowerthanourpreviousestimate(Areff= 0.3reft),buta part of the discrepancycan clearly be attributedto humidityeffects. Becauseof the difficulty in makingan accuratevalidationof reft, we also comparedthe Angstr6m exponent which is representative of the aerosol size distribution (Table 4). Except for July 20 the agreementis good. Recall that the Angstr6m exponent is not computedfrom exactly the same wavelengthsfor the different instruments, which may explain

0.60 0.50

0.40

/ / / / /

0.30

/

/

/

0.20

/ / /

/

0.10

/

0.00 0.00

,

0.10

I

0.20

,

I

0.30

,

I

0.40

,

I

,

0.50

0.60

Measured Optical Thickness()•=550nm) Figure 12. Comparisonof the aerosoloptical thicknessat t = 550 nm measuredby the CIMEL instrument,aboardthe C131A aircraft and values derived from the MAS all coincident measurements.

instrument

for

TANRl• ET AL.: RETRIEVAL OF AEROSOLPROPERTIESOVER THE OCEAN

2275

Table 4. Angstr6mExponentDerived from the MAS Instrument(C•MAS) and from the SurfaceMeasurements:Sun photometer(c•sp) and AirborneSunphotometer(C•AATS_6)

Flight

XMAS

{XMAS

OtAATS_ 6

asp

Comment

Track

July 17 10 11 12

x = 0.57 x = 0.49 x = 0.55

1.60 1.65 1.68

1.58 1.58 1.58

-

goodmatches(time and position) between ER-2 (track 10) and C-131A; high spatial variability

0.55 (ship)

aerosolopticalthicknessis too small to derive the size

July 20 04 15

x < 0.05 x < 0.07

-0.15 -0.11

1.13 -

distribution

July 24 08

x = 0.49

1.62 a

x = 0.37

1.47 •

1.73

goodmatches(time and position) between ER-2 and C-131A;

presenceof an aerosol plume

July 25 13

x = 0.29

1.23

1.60

goodmatch(position)between ER-2 andC-131A, gap in time (45min); high spatialvariability

July 26 04

x = 0.28

0.97

1.24 (Bermuda Islands)

good match (time) between ER-2 and surface instrumentation;not the same location

t•MAS is computed at 3,=550and865 nm,t•Sp at 3,= 499 and865 nm,O•AATS_ 6 at 3,= 525 and1025nm.

a Corresponds to thetwozonesdefinedin section 4.3.

b Corresponds tothetwozones defined insection 4.3. part of the discrepancy. In addition, the larger discrepancy (Aot= 0.35) occurred when there was a delay of 45 min betweenmeasurements.When there is both a temporal and a geographicalmatch, the differenceis aroundA•z = 0.10.

wind speed. Larger values of x are obtained the smaller the wind speed because the corresponding foam contribution decreases and this decreasehas to be compensatedby a larger optical thickness. Also, the effective particle radius (Figure

13b) is hardly affected.Sincethe aerosolcontributionis larger in the visible when the wind speedis smaller, the resulting 6. Discussion stronger spectral dependence is then interpreted as a In this section we discussthe impactsof potential sources contribution of smaller particles. Similar computationswere performedfor the clearestday, July 20 (Figure 14). The above of errors, such as surfaceconditions, the glint mask, and the conclusions were still valid althoughthe impact was twice as ozone correction. large (Axa = +0.02). Since the relative contribution of the atmosphereis smallerin this case, resultsare slightly more 6.1. Impact of Wind Speed affectedby the changesin the surfacereflectance. To check the impact of wind speed on the retrieval, we performedthe inversion for the three values included in the

LUT (i.e., 2.0, 6.0, and 10.0 ms'i). In Figure13, results obtainedfor the two extremevalues(2.0 and 10.0 ms-1) are plotted as a function of the resultsobtained for the mean value

(6.0 ms'i). July 26 wasselected because therewasno glint within the image, and the opticalthicknesses displayeda large variability. The impact on 'ca (Figure 13a) is very limited, around0.01, which meansthat the coupling term betweenthe atmosphereand the oceansurfaceis almost independentof the

6.2.

Impact

of the Glint

Mask

We have noticed previously that the glint mask is very important since the reflection of the directSunbeam is very dependenton the wind speed.To checkthis, we selecteda day affected by the presence of the glint (July 17) and we performed the inversion for the three wind speeds. The differencesin the optical thickness between two inversions are shownin Figure 15 as a functionof the glint angle. Angles

2276

TANRI•ET AL.:RETRIEVALOF AEROSOLPROPERTIES OVERTHE OCEAN These results show that the Sun glint requiresa very large

0.50

' 0.40

I

'

I

'

I

'•

maskaroundthe specularreflection. it has to be testedagain on actual MODIS data since there may be contamination by stray light in the airborne MAS version, which increasesthe sensitivity.

ß v=2ms -1

6.3.

0.30

0.20

••,c• 0.10

July 26,Track 04

0.20

0.30

0.40

0.50

Retrieved Optical Thickness (atv=6ms 4) 0.70

I

'

I

of Ozone Correction

The impact of the ozone correction is estimated by assuminga variation of 80 Dobson units in the apparent ozone content; this correspondsto a variation of 27 Dobson

(a) '•• , I , I )•=550nm , I.... (a)

0.10

Impact

'

unit for an air mass of 3. Recall that the difference

between the

two extreme models (tropical and sub arctic winter) is around 25 Dobson

units.

The quality of the inversion is clearly not affected by variationsin the ozone content.For the clearestday (July 20) the differencebetween both optical thicknesses is less than

0.01 at )• = 550 nm (Figure 16). 0.60

o

v=10ms4

ß

v=2ms -1

7.

Conclusion

0.50

The potential of the MODIS data for retrieving aerosol parameters (optical thickness, the effective radius of the particle size distribution,and the ratio of the accumulationand coarsemodes),is demonstrated by the presentstudy.The data acquiredduringTARFOX confirm the expectedaccuracyof the retrievedsize distribution. The accuracyof the derivedvalues of the optical thicknessis even better than anticipated.

o o

0.40

o

0.30

July26, Track 04 (b)

0.20



0.10

0.10

0.20

0.30

0.40

0.50

I

0.60

0.70

0.10

RetrievedEffectiveRadius (at v=6ms-1) Figure

0.00

13. (a) Aerosol optical thickness derived from the

MAS instrumentassuminga wind speedof 2.0 ms-• and

10.0ms '•asa function oftheoptical thickness assuming a '•' =

-0.10

windspeedof 6.0 ms-•.Thecomparison is doneat 550 nm for July 26. (b) As in Figure 13a but for the effective particle radius.

m ii

July 17, Track12

-0.20

v=6.0ms4 and v=2.0ms4

! smaller

than

30 ø are not

shown

because this

was the

value

previouslyselectedfor the mask. It is clear that directionsare increasingly affected as specular reflection is approaching. For an angle of 35ø from the specularreflection the impact is aroundAxa = + 0.10, and the impact is completelynegligible for anglesgreater than 45 ø. 0.10

-0.30

-

ß

ß

(a) -0.40

I

I

20.0 25.0 30.0 0.10

I

0.00



35.0

•1

40.0

45.0

50.0

55.0

I

I

I

I

~ß ..........

...............

&

60.0



0.O9

-010

0.08 0.07

-0.20

0.06 -0.30

0.05 0.04 0.03

0.02

%

-

-

-

July17,Track12

v=6.0rns qand v=10.0ms 4

-



July 20,Track 15

-/ 0.02

-]

k=õ50nm (a)1 0.03

0.04

0.05

0.06

0.07

0.08

0.09

0.10

Retrieved Optical Thickness (atv=6ms q) Figure 14. As in Figure 13a but for July 20.

-0.40

I

I

I

I

I

(b) i

I

20.00 25.00 30.00 35.00 40.00 45.00 50.00 55.00 60.00

Glint Angle (ø) Figure 15. Differences in the retrieved optical thickness as a function of the glint angle assumingdifferent wind speed

valuesin the inversion.(a) Wind speedsof 2.0 ms-• and 6.0 ms-•. (b) Windspeeds of 6.0 ms4and10.0ms4.

TANRl• ET AL.: RETRIEVAL OF AEROSOLPROPERTIESOVER THE OCEAN

2277

Earth ObservingSystemera, J. Geophys.Res., 102, 17,081-17,107, 1997.

Griggs,M., Measurements of atmospheric aerosolopticalthicknessover waterusingERTS-1 data,J. Air Pollut. ControlAssoc.,25, 622-626,

0.11

1975.

0.10

Herman, J., P. K. Bhartia, O. Torres, C. Hsu, C. Seftor, and E. Celarier, Global distributionof UV-absorbing aerosolsfrom Nimbus 7/TOMS data,J. Geophys.Res., 102, 16,911-16,922, 1997. Hobbs, P. V., An overview of the University of Washingtonairborne

0.09

0.08

measurements and resultsfrom the TroposphericAerosol Radiative Forcing ObservationalExperiment (TARFOX), J. Geophys. Res.,

0.07

this issue.

0.06

-

0.05

_//if-

0.04

0.04

_•-

0.05

0.06

July20,Track15 -

Holben, B. N., et al., AERONET--A federated instrument network and data archive for aerosol characterization, Remote Sens. Environ., 66, 1-16, 1998.

L=550nm _ 0.07

0.08

0.09

0.10

0.11

AerosolOpticalThickness Figure 16. Impact of the ozone correctionon the aerosol opticalthickness.The comparison is doneat 550 nm for July 20.

Althoughthe presentresultsneed to be confirmedfor other aerosol types (e.g., dust and smoke), the MODIS aerosol productsover ocean can be usedwith some confidence.The problem of the glint mask is still unresolved,but actual MODIS data are requiredto solveit.

Husar, R. B., J. M. Prospero,and L. L. Stowe, Characterizationof troposphericaerosolsover the ocean usingthe NOAA advanced very high resolution radiometer optical thickness operational product,J. Geophys.Res., 102, 16,889-16,909, 1997.

Ismail, S., E. V. Browell, A. S. Moore, W. C. Edwards, K. Browns, S. A.

Kooi, V. G. Brackett, and M. B. Clayton, LASE measurementsof aerosol, cloud and water vapor profiles during TARFOX field experiment,Eos Trans.AGU, 78(17), SpringMeet. Suppl.,S82, 1997. Jankowiak, I., and D. Tanrr, Climatology of Saharan dust events observedfrom Meteosat imagery over Atlantic Ocean, Method and preliminaryresults,J. Clim., 5, 646-656, 1992. Kaufman, Y. J., D. Tanrr, H. R. Gordon, T. Nakajima, J. Lenoble, R. Frouin,H. Grassl,B. M. Herman,M. King, andP.M. Teillet, Passive remote sensingof troposphericaerosoland atmosphericcorrection for the aerosoleffects,J. Geophys.Res., 102, 16,815-16,830, 1997. King, M. D., et al., Airborne scanningspectrometerfor remote sensing of cloud, aerosol, water vapor and surface properties,J. Atmos. Oceanic Technol., 13, 777-794, 1996.

Koepke, P., Effective reflectanceof oceanic whitecaps,Appl. Opt., 23, 1816-1823, 1984.

Acknowledgments.This researchwas conductedas part of the Koepke, P., and H. Quenzel, Turbidityof the atmospheredetermined TroposphericAerosol Radiative Forcing ObservationalExperiment from satellitecalculationof optimumviewing geometry,J. Geophys. (TARFOX),

which is a contribution to the International Global

Res., 84, 7847-7855, 1979.

AtmosphericChemistry(IGAC) core project of the International Martonchik, J. V., and D. J. Diner, Retrieval of aerosol optical Geosphere-Biosphere Programme(IGBP). Financialsupportfor the propertiesfrom multi-angle satellite imagery, IEEE Trans. Geosci. measurementsand analyses was provided by the U.S. National Remote Sens., 30, 223-230, 1992. Aeronauticsand Space Administration,National Science Foundation, Matsumoto, T., P. B. Russell, C. Mina, and W. Van Ark, Airborne Office of Naval Research, and National Oceanographic and tracking sunphotometer,J. Atmos. Oceanic Technol., 4, 336-339, 1987. Atmospheric Administration, by theU.K. MeteorologicalOffice andby the French Science Foundation (Centre National de la Recherche Mekler, Y., H. Quenzel, G. Ohring, and I. Marcus, Relative atmospheric aerosolcontentfrom ERTS observations, J. Geophys. Scientifique(CNRS)) and the French Space Administration(Centre Res., 82,967-972, 1977. Nationald'EtudesSpatiales(CNES)). We wouldlike to thankthe NASA ER-2 crew, especiallypilotsJan Nystromand Ken Broda for their M6szaros,E., Seasonaland diurnal variations of the size distributionof excellentwork andextremeflexibility duringthe TARFOX deployment. atmospheric sulfateparticles,Tellus,22, 235-238, 1970. Mishchenko, M. I., and L. D. Travis, Satellite retrieval of aerosol Thanksalsoto Rong-RongLi for datapreparation andprocessing.

propertiesover the oceanusingpolarization as well as intensityof References Ackerman, S. A., K. I. Strabala, W. P. Menzel, R. A. Frey, C. C. Moeller, and L. E. Gumley, DiscriminatingClear-sky from Clouds with MODIS, J. Geophys.Res.,in press,1998.

reflectedsunlight,J. Geophys.Res.,102, 16,989-17,014,1997. Moulin, C., F. Guillard, F. Dulac, and C.E. Lambert, Long-term daily monitoringof Saharandustload over oceanusingMeteosatISCCPB2 data, Methodologyand preliminaryresultsfor 1983-1994in the Mediterranean,J. Geophys.Res.,102, 16,947-16,958,1997.

Nakajima,T., and A. Higurashi,AVHRR remotesensingof aerosol opticalpropertiesin the PersianGulf region, summer1991, J. Geophys. Res.,102, 16,935-16,946,1997.

Carlson,T. N., Atmosphericturbidityin Saharandust outbreaksas determinedby analysesof satellitebrightnessdata, Mon. Weather Norton,C. C., F. R. Mosher,B. Hinton, D. W. Martin, D. Santek,and W. Rev., 107, 322-335, 1979. Kuhlow, A model for calculatingdesert aerosol turbidity over Cox, C., and W. Munk, Statisticsof the sea surface derived from sun oceansfrom geostationary satellitedata,J. Appl.Meteorol.,19, 633glitter,J. Mar. Res.,13, 198-208,1954. 642, 1980. Deschamps, P. Y., F. M. Br6on,M. Leroy,A. Podaire,A. Bricaud,J. C. Remer, L. A., Y. J. Kaufman, and B. N. Holben, Interannual variation Buriez, and G. Sbze, The POLDER mission: Instrument of aerosolphysicaland opticalcharacteristics on the east coast of

the United States,J. Geophys.Res.this issue. Russell,P. B., et al., Pinatuboand pre-Pinatubooptical depth spectra: Mauna Loa measurements,comparisons,inferred particle size Diner, D. J., et al., MISR: A Multiangleimagingspectroradiometer for distributions,radiative effects, and relationshipto lidar data, J. geophysicaland climatologicalresearchfrom EOS, IEEE Trans. Geosci. Remote Sens., 27, 200-214, 1989. Geophys.Res., 98, 22,969-22,985, 1993. Fraser, R. S., Satellite measurementof mass of Saharan dust in the Russell,P. B., P. V. Hobbs, and L. L. Stowe, Aerosol propertiesand radiativeeffectsin the U.S. east coasthaze plume:An overview of atmosphere, Appl. Opt.,15, 2471-2479,1976. the Tropospheric Aerosol Radiative Forcing Observational Frouin,R., M. Schwindling,and P. Y Deschamps,Spectralreflectance of sea foam in the visible and near infrared: In situ measurements Experiment(TARFOX), J. Geophys.Res.,thisissue(a). and remote sensing implications, J. Geophys. Res., 101, Russell,P. B., J. M. Livingston,P. Hignett,S. Kinne, J. Wong, A. Chien, 14,361-14,371, 1996. R. Bergstrom., P. Durkee, and P. V. Hobbs, Aerosol-induced radiativeflux changesoff the U.S. mid-Atlantic coast: Comparison Gordon,H. R., Atmospheric corrections of oceancolorimageryin the characteristics and scientific objectives, IEEE RemoteSens., 32, 598-615, 1994.

Trans. Geosci.

2278

TANRI• ET AL.: RETRIEVALOF AEROSOLPROPERTIES OVERTHE OCEAN

of valuescalculatedfrom Sun photometerand in situ data with those measuredby airbornepyranometer,J. Geophys.Res., this issue(b). Salomonson,V. V., W. L. Barnes, P. W. Maymon, H. E. Montgomery, and H. Ostrow, MODIS: advancedfacility instrumentfor studiesof the earth as a system, IEEE Trans. Geosci.Remote Sens.,27, 145153, 1989.

Stowe,L. L., A.M. Ignatov, and R. R. Singh,Development,validation, and potential enhancementsto the second-generationoperational aerosolproductat the National EnvironmentalSatellite, Data, and Information Service of the National Oceanic and Atmospheric Administration,J. Geophys.Res., 102, 16,923-16,934,1997. Tanr6, D., M. Herman, and Y. J. Kaufman, Information on the aerosol

size distributioncontainedin the solar reflected spectralradiances, J. Geophys.Res., 101, 19,043-19,060, 1996. Tanr6, D., Y. J. Kaufman, M. Herman, and S. Mattoo, Remotesensing of aerosolpropertiesover oceans usingthe MODIS/EOS spectral radiances,J. Geophys.Res., 102, 16,971-16,988,1997. Travis, L. D., EOSP: Earth observingscanningpolarimeter,in EOS ReferenceHandbook1993, editedby G. Asrarand D.J. Dokken, pp. 74-75, NASA, 1993.

Whitlock, C. H., D. S. Bartlett, and E. A. Gurganus,Sea foam reflectance and influence on optimum wavelength for remote sensingof oceanaerosols,Geophys.Res.Lett., 9, 719-722, 1982. D. Tanr6, Laboratoired'OptiqueAtmosph6rique,UA CNRS 713, U.S.T. de Lille, Bat. P5, 59655 - Villeneuve d'Ascq,France (email: Didier.Tanre@univ-lille 1.fr.) Y. J. Kaufman and L. A. Remer, LaboratoryFor Atmospheres, NASA Goddard Space Flight Center, Code 913, Greenbelt, MD 20771. (email: [email protected].) S. Mattoo, ScienceApplicationCorporation,Vienna,VA 22180. P. V. Hobbs,Universityof Washington,Seattle,WA 98195. J. M. LivingstonandP. B. Russell,NASA Ames ResearchCenter, Moffett Field, CA 94035.

A. Smirnov,ScienceSystemsand Applications,Inc., Laboratory for Terrestrial Physics,NASA GoddardSpace Flight Center, Code 923, Greenbelt, MD 20771.

(ReceivedAugust28, 1998;revisedNovember2, 1998; acceptedNovember 3, 1998.)