What Does Nature Tell Us About Anthropogenic Aerosol Indirect Effects?
No abstract available
cloud regimes↗
Engineering topics
Oreopoulos, L.
Publications and source records attributed to Oreopoulos, L..
Modeling the Influences of Aerosols on Pre-Monsoon Circulation and Rainfall over Southeast Asia
We conduct several sets of simulations with a version of NASA's Goddard Earth Observing System, version 5, (GEOS-5) Atmospheric Global Climate Model (AGCM) equipped with a two-moment cloud microphysical scheme to understand the role of biomass burning aerosol (BBA) emissions in Southeast Asia (SEA) in the pre-monsoon period of February-May. Our experiments are designed so that both direct and indirect aerosol effects can be evaluated. For climatologically prescribed monthly sea surface temperatures, we conduct sets of model integrations with and without biomass burning emissions in the area of peak burning activity, and with direct aerosol radiative effects either active or inactive. Taking appropriate differences between AGCM experiment sets, we find that BBA affects liquid clouds in statistically significantly ways, increasing cloud droplet number concentrations, decreasing droplet effective radii (i.e., a classic aerosol indirect effect), and locally suppressing precipitation due to a deceleration of the autoconversion process, with the latter effect apparently also leading to cloud condensate increases. Geographical re-arrangements of precipitation patterns, with precipitation increases downwind of aerosol sources are also seen, most likely because of advection of weakly precipitating cloud fields. Somewhat unexpectedly, the change in cloud radiative effect (cloud forcing) at surface is in the direction of lesser cooling because of decreases in cloud fraction. Overall, however, because of direct radiative effect contributions, aerosols exert a net negative forcing at both the top of the atmosphere and, perhaps most importantly, the surface, where decreased evaporation triggers feedbacks that further reduce precipitation. Invoking the approximation that direct and indirect aerosol effects are additive, we estimate that the overall precipitation reduction is about 40% due to the direct effects of absorbing aerosols, which stabilize the atmosphere and reduce surface latent heat fluxes via cooler land surface temperatures. Further refinements of our two-moment cloud microphysics scheme are needed for a more complete examination of the role of aerosol-convection interactions in the seasonal development of the SEA monsoon.
SE Asia Monsoon↗
A Novel Method for Estimating Shortwave Direct Radiative Effect of Above-Cloud Aerosols Using CALIOP and MODIS Data
This paper describes an efficient and unique method for computing the shortwave direct radiative effect (DRE) of aerosol residing above low-level liquid-phase clouds using CALIOP and MODIS data. It accounts for the overlapping of aerosol and cloud rigorously by utilizing the joint histogram of cloud optical depth and cloud top pressure. Effects of sub-grid scale cloud and aerosol variations on DRE are accounted for. It is computationally efficient through using grid-level cloud and aerosol statistics, instead of pixel-level products, and a pre-computed look-up table in radiative transfer calculations. We verified that for smoke over the southeast Atlantic Ocean the method yields a seasonal mean instantaneous shortwave DRE that generally agrees with more rigorous pixel-level computation within 4. We have also computed the annual mean instantaneous shortwave DRE of light-absorbing aerosols (i.e., smoke and polluted dust) over global ocean based on 4 yr of CALIOP and MODIS data. We found that the variability of the annual mean shortwave DRE of above-cloud light-absorbing aerosol is mainly driven by the optical depth of the underlying clouds.
radiative forcing↗
A Novel Method for Estimating Shortwave Direct Radiative Effect of Above-cloud Aerosols over Ocean Using CALIOP and MODIS Data
This paper describes an efficient and unique method for computing the shortwave direct radiative effect (DRE) of aerosol residing above low-level liquid-phase clouds using CALIOP and MODIS data. It accounts for the overlapping of aerosol and cloud rigorously by utilizing the joint histogram of cloud optical depth and cloud top pressure. Effects of sub-grid scale cloud and aerosol variations on DRE are accounted for. It is computationally efficient through using grid-level cloud and aerosol statistics, instead of pixel-level products, and a pre-computed look-up table in radiative transfer calculations. We verified that for smoke over the southeast Atlantic Ocean the method yields a seasonal mean instantaneous shortwave DRE that generally agrees with more rigorous pixel-level computation within 4%. We have also computed the annual mean instantaneous shortwave DRE of light-absorbing aerosols (i.e., smoke and polluted dust) over global ocean based on 4 yr of CALIOP and MODIS data. We found that the variability of the annual mean shortwave DRE of above-cloud light-absorbing aerosol is mainly driven by the optical depth of the underlying clouds.
CALIOP and MODIS↗
Using ISCCP Weather States to Decompose Cloud Radiative Effects
The presentation will examine the shortwave (SW) and longwave (LW) cloud radiative effect CRE (aka "cloud radiative forcing") at the top-of-the-atmosphere and surface of ISCCP weather states (aka "cloud regimes") in three distinct geographical zones, one tropical and two mid-latitude. Our goal is to understand and quantify the contribution of the different cloud regimes to the planetary radiation budget. In the tropics we find that the three most convectively active states are the ones with largest SW, LW and net TOA CRE contributions to the overall daytime tropical CRE budget. They account for 59%, 71% and 55% of the total CRE, respectively. The boundary layer-dominated weather states account for only 34% of the total SW CRE and 41% of the total net CRE, so to focus only on them in cloud feedback studies may be imprudent. We also find that in both the northern and southern midlatitude zones only two weather states, the first and third most convectively active with large amounts of nimbostratus-type clouds, contribute ",40% to both the SW and net TOA CRE budgets, highlighting the fact that cloud regimes associated with frontal systems are not only important for weather (precipitation) but also for climate (radiation budget). While all cloud regimes in all geographical zones have a slightly larger SFC than TOA SW CRE, implying cooling of the surface and slight warming of the atmosphere, their LW radiative effects are more subtle: in the tropics the weather states with plentiful high clouds warm the atmosphere while those with copious amounts of low clouds cool the atmosphere. In both midlatitude zones only the weather states with peak cloud fractions at levels above 440 mbar warm the atmosphere while all the rest cool it. These results make the connection of the contrasting CRE effects to the atmospheric dynamics more explicit - "storms" tend to warm the atmosphere whereas fair weather clouds cool it, suggesting a positive feedback of clouds on weather systems. The breakdown of CRE by cloud regime are however not entirely similar between the two midlatitude zones. Despite the existence of an additional state in the nort!lern midlatitudes, only four weather states have net daytime CREs with absolute values above 100 Watts per square meter compared to six in the south. This reminds us that the environment where clouds occur also has a crucial role in determining their radiative effects. All the above make evident that reproducing grand averages of current CRE by climate models in only part of the challenge. If existing cloud regimes and shifts in their distributions and frequency of occurrence in a changed climate are not properly simulated, the radiative role of clouds will not be adequately predicted.
Oreopoulos, L.↗
The Performance of Current Atmospheric Radiation Codes in Phase I of CIRC
The Continual Intercomparison of Radiation Codes (CIRC) is intended as an evolving and regularly updated reference source for evaluation of radiative transfer (RT) codes used in Global Climate Models and other atmospheric applications. In our presentation we will discuss our evaluation of the performance of 13 shortwave and 11 longwave RT codes that participated in Phase I of CIRC. CIRC differs from previous intercomparisons in that it relies on an observationally validated catalogue of cases. The seven CIRC Phase I baseline cases, five cloud-free, and two with overcast liquid clouds, are built around observations by the Atmospheric Radiation Measurements (ARM) program that satisfy the goals .of Phase I, namely to examine RT model performance in realistic, yet not overly complex, atmospheric conditions. Besides the seven baseline cases, additional idealized "subcases" are also examined to facilitate interpretation of model errors. We will quantify individual model performance with respect to reference line-by-line calculations, and will also highlight RT code behavior for conditions of doubled CO2 , aspects of utilizing a spectral specification of surface albedo, and the impact of the inclusion of scattering in the thermal infrared. Our analysis suggests that RT codes should work towards improving their calculation of diffuse shortwave flux, shortwave absorption, treatment of spectral surface albedo, and shortwave CO2 forcing. Despite practical difficulties in comparing our results to previous results by the Intercomparison of Radiation Codes in Climate Models (ICRCCM) conducted about 20 years ago, it appears that the current generation of RT codes do indeed perform better than the codes of the ICRCCM era. By enhancing the range of conditions under which participating codes are tested, future CIRC phases will hopefully allow even more rigorous examination of RT code performance.
Oreopoulos, L.↗
Sensitivity of Cirrus and Mixed-phase Clouds to the Ice Nuclei Spectra in McRAS-AC: Single Column Model Simulations
The salient features of mixed-phase and ice clouds in a GCM cloud scheme are examined using the ice formation parameterizations of Liu and Penner (LP) and Barahona and Nenes (BN). The performance of LP and BN ice nucleation parameterizations were assessed in the GEOS-5 AGCM using the McRAS-AC cloud microphysics framework in single column mode. Four dimensional assimilated data from the intensive observation period of ARM TWP-ICE campaign was used to drive the fluxes and lateral forcing. Simulation experiments where established to test the impact of each parameterization in the resulting cloud fields. Three commonly used IN spectra were utilized in the BN parameterization to described the availability of IN for heterogeneous ice nucleation. The results show large similarities in the cirrus cloud regime between all the schemes tested, in which ice crystal concentrations were within a factor of 10 regardless of the parameterization used. In mixed-phase clouds there are some persistent differences in cloud particle number concentration and size, as well as in cloud fraction, ice water mixing ratio, and ice water path. Contact freezing in the simulated mixed-phase clouds contributed to transfer liquid to ice efficiently, so that on average, the clouds were fully glaciated at T approximately 260K, irrespective of the ice nucleation parameterization used. Comparison of simulated ice water path to available satellite derived observations were also performed, finding that all the schemes tested with the BN parameterization predicted 20 average values of IWP within plus or minus 15% of the observations.
Betancourt, R. Morales↗
Radiative Impacts of Cloud Heterogeneity and Overlap in an Atmospheric General Circulation Model
The radiative impacts of introducing horizontal heterogeneity of layer cloud condensate, and vertical overlap of condensate and cloud fraction are examined with the aid of a new radiation package operating in the GEOS-5 Atmospheric General Circulation Model. The impacts are examined in terms of diagnostic top-of-the-atmosphere shortwave (SW) and longwave (LW) cloud radiative effect (CRE) calculations for a range of assumptions and parameter specifications about the overlap. The investigation is conducted for two distinct cloud schemes, the one that comes with the standard GEOS-5 distribution, and another which has been recently used experimentally for its enhanced GEOS-5 distribution, and another which has been recently used experimentally for its enhanced cloud microphysical capabilities; both are coupled to a cloud generator allowing arbitrary cloud overlap specification. We find that cloud overlap radiative impacts are significantly stronger for the operational cloud scheme for which a change of cloud fraction overlap from maximum-random to generalized results to global changes of SW and LW CRE of approximately 4 Watts per square meter, and zonal changes of up to approximately 10 Watts per square meter. This is because of fewer occurrences compared to the other scheme of large layer cloud fractions and of multi-layer situations with large numbers of atmospheric being simultaneously cloudy, conditions that make overlap details more important. The impact on CRE of the details of condensate distribution overlap is much weaker. Once generalized overlap is adopted, both cloud schemes are only modestly sensitive to the exact values of the overlap parameters. We also find that if one of the CRE components is overestimated and the other underestimated, both cannot be driven towards observed values by adjustments to cloud condensate heterogeneity and overlap alone.
Oreopoulos, L.↗
Performance of McRAS-AC in the GEOS-5 AGCM: Part 1, Aerosol-Activated Cloud Microphysics, Precipitation, Radiative Effects, and Circulation
A revised version of the Microphysics of clouds with Relaxed Arakawa-Schubert and Aerosol-Cloud interaction (McRAS-AC), including, among others, the Barahona and Nenes ice nucleation parameterization, is implemented in the GEOS-5 AGCM. Various fields from a 10-year long integration of the AGCM with McRAS-AC were compared with their counterparts from an integration of the baseline GEOS-5 AGCM, and with satellite data as observations. Generally using McRAS-AC reduced biases in cloud fields and cloud radiative effects are much better over most of the regions of the Earth. Two weaknesses are identified in the McRAS-AC runs, namely, too few cloud particles around 40S-60S, and too high cloud water path during northern hemisphere summer over the Gulf Stream and North Pacific. Sensitivity analyses showed that these biases potentially originated from biases in the aerosol input. The first bias is largely eliminated in a sensitivity test using 50% smaller aerosol particles, while the second bias is much reduced when interactive aerosol chemistry was turned on. The main drawback of McRAS-AC is dearth of low-level marine stratus clouds, probably due to lack of dry-convection, not yet implemented into the cloud scheme. Despite these biases, McRAS-AC does simulate realistic clouds and their optical properties that can improve with better aerosol-input and thereby has the potential to be a valuable tool for climate modeling research because of its aerosol indirect effect simulation capabilities involving prediction of cloud particle number concentration and effective particle size for both convective and stratiform clouds is quite realistic.
Sud, Y. C.↗
Allowing for Horizontally Heterogeneous Clouds and Generalized Overlap in an Atmospheric GCM
While fully accounting for 3D effects in Global Climate Models (GCMs) appears not realistic at the present time for a variety of reasons such as computational cost and unavailability of 3D cloud structure in the models, incorporation in radiation schemes of subgrid cloud variability described by one-point statistics is now considered feasible and is being actively pursued. This development has gained momentum once it was demonstrated that CPU-intensive spectrally explicit Independent Column Approximation (lCA) can be substituted by stochastic Monte Carlo ICA (McICA) calculations where spectral integration is accomplished in a manner that produces relatively benign random noise. The McICA approach has been implemented in Goddard's GEOS-5 atmospheric GCM as part of the implementation of the RRTMG radiation package. GEOS-5 with McICA and RRTMG can handle horizontally variable clouds which can be set via a cloud generator to arbitrarily overlap within the full spectrum of maximum and random both in terms of cloud fraction and layer condensate distributions. In our presentation we will show radiative and other impacts of the combined horizontal and vertical cloud variability on multi-year simulations of an otherwise untuned GEOS-5 with fixed SSTs. Introducing cloud horizontal heterogeneity without changing the mean amounts of condensate reduces reflected solar and increases thermal radiation to space, but disproportionate changes may increase the radiative imbalance at TOA. The net radiation at TOA can be modulated by allowing the parameters of the generalized overlap and heterogeneity scheme to vary, a dependence whose behavior we will discuss. The sensitivity of the cloud radiative forcing to the parameters of cloud horizontal heterogeneity and comparisons of CERES-derived forcing will be shown.
Lee, D.↗
Precipitation Characteristics of ISCCP Cloud Regimes for Improving Model Hydrological Budgets
The key in unraveling relationships between precipitation and atmospheric circulations is their common linkage to clouds. Clouds can be described in a variety of ways and several approaches can be adopted to examine their connections to precipitation. We claim that when cloud regimes (aka weather states) from the International Satellite Cloud Climatology Project (ISCCP) are used to conditionally sample/sort and average precipitation data, useful insight and GCM-appropriate diagnostics on the origins and distribution of precipitation can be obtained. The ISCCP cloud regimes are mesoscale (2.5 ) cloud mixtures determined by cluster analysis on joint histograms of cloud optical thickness and cloud top pressure inferred from geostationary and polar orbiter satellite passive retrievals. The ISCCP cloud regime data are combined with GPCP IDD merged surface precipitation data and/or higher temporal and spatial resolution TRMM Multi-Satellite Precipitation Analysis (TMPA) data. The analysis is performed separately for three geographical zones, tropics, and northern/southern midlatitudes (for GPCP; only the tropics can be examined with TMPA data). Our presentation aspires to provide answers to the following questions: (l) What is the mean and variability of surface precipitation produced by each cloud regime at the time of regime occurrence? (2) What is the relative contribution of each cloud regime to the total precipitation within its geographical zone? (3) What is the geographical distribution of precipitation corresponding to particular cloud regime? (4) To what extent are the cloud regimes distinct in terms of their precipitation characteristics and is the regime ordering in terms of convective strength consistent with the observed precipitation intensity?
Lee, D.↗
Advancements in the Representation of Cloud-Aerosol Microphysics in the GEOS-5 AGCM
Despite numerous challenges, the physical parameterization of cloud-aerosol interactions in atmospheric GCMs has become a top priority for advancement because of our need to simulate and understand past, current, and future indirect effects of aerosols on clouds. The challenges stem from the involvement of wide range of cloud-scale dynamics and aerosol activation physical processes. Cloud dynamics modulate cloud areal extent and condensate, while aerosol activation depends on aerosol mass load, size distribution, internal mixing state, and nucleating properties, and ultimately determines cloud optical properties via particle sizes. Both macro- and micro-scale processes are obviously important for cloud-radiation interactions. We will present the main features of cloud microphysical properties in the GEOS- 5 Atmospheric GCM (AGCM) as simulated by the McRAS-AC (Microphysics of Clouds with Relaxed Arakawa-Schubert and Aerosol-Cloud interaction) scheme. McRAS-AC uses Fountoukis and Nenes (2005) aerosol activation for liquid clouds, and has an option for either Liu and Penner (2005) or Barahona and Nenes (2008, 2009) aerosol activation for ice clouds. Aerosol loading (on-line or climatological) comes from GOCART, with an assumed log-normal size distribution. Other features of McRAS-AC are level-by-level cloud-scale thermodynamics, and Seifert-Beheng (2001)-type precipitation microphysics, particularly from moist convection. Results from Single-Column Model simulations will be shown to demonstrate how cloud radiative properties, lifetimes, and precipitation are influenced by different parameterization assumptions. Corresponding fields from year-long simulations of the full AGCM will also be presented with geographical distributions of cloud effective particle sizes compared to satellite retrievals. While the primary emphasis will be on current climate, simulation results with perturbed aerosol loadings will also be shown to expose the radiative sensitivity of the microphysical parameterization.
Lee, D.↗
Performance of Goddard Earth Observing System GCM Column Radiation Models under Heterogeneous Cloud Conditions
We test the performance of the shortwave (SW) and longwave (LW) Column Radiation Models (CORAMs) of Chou and collaborators with heterogeneous cloud fields from a global single-day dataset produced by NCAR's Community Atmospheric Model with a 2-D CRM installed in each gridbox. The original SW version of the CORAM performs quite well compared to reference Independent Column Approximation (ICA) calculations for boundary fluxes, largely due to the success of a combined overlap and cloud scaling parameterization scheme. The absolute magnitude of errors relative to ICA are even smaller for the LW CORAM which applies similar overlap. The vertical distribution of heating and cooling within the atmosphere is also simulated quite well with daily-averaged zonal errors always below 0.3 K/d for SW heating rates and 0.6 K/d for LW cooling rates. The SW CORAM's performance improves by introducing a scheme that accounts for cloud inhomogeneity. These results suggest that previous studies demonstrating the inaccuracy of plane-parallel models may have unfairly focused on worst scenario cases, and that current radiative transfer algorithms of General Circulation Models (GCMs) may be more capable than previously thought in estimating realistic spatial and temporal averages of radiative fluxes, as long as they are provided with correct mean cloud profiles. However, even if the errors of the particular CORAMs are small, they seem to be systematic, and the impact of the biases can be fully assessed only with GCM climate simulations.
Oreopoulos, L.↗
Overlap Properties of Clouds Generated by a Cloud Resolving Model
In order for General Circulation Models (GCMs), one of our most important tools to predict future climate, to correctly describe the propagation of solar and thermal radiation through the cloudy atmosphere a realistic description of the vertical distribution of cloud amount is needed. Actually, one needs not only the cloud amounts at different levels of the atmosphere, but also how these cloud amounts are related, in other words, how they overlap. Currently GCMs make some idealized assumptions about cloud overlap, for example that contiguous cloud layers overlap maximally and non-contiguous cloud layers overlap in a random fashion. Since there are difficulties in obtaining the vertical profile of cloud amount from observations, the realism of the overlap assumptions made in GCMs has not been yet rigorously investigated. Recently however, cloud observations from a relatively new type of ground radar have been used to examine the vertical distribution of cloudiness. These observations suggest that the GCM overlap assumptions are dubious. Our study uses cloud fields from sophisticated models dedicated to simulate cloud formation, maintenance, and dissipation called Cloud Resolving Models . These models are generally considered capable of producing realistic three-dimensional representation of cloudiness. Using numerous cloud fields produced by such a CRM we show that the degree of overlap between cloud layers is a function of their separation distance, and is in general described by a combination of the maximum and random overlap assumption, with random overlap dominating as separation distances increase. We show that it is possible to parameterize this behavior in a way that can eventually be incorporated in GCMs. Our results seem to have a significant resemblance to the results from the radar observations despite the completely different nature of the datasets. This consistency is encouraging and will promote development of new radiative transfer codes that will estimate the radiation effects of multi-layer cloud fields more accurately.
Oreopoulos, L.↗
Consistency of ARESE II Cloud Absorption Estimates and Sampling Issues
Data from three cloudy days (March 3, 21, 29, 2000) of the ARM Enhanced Shortwave Experiment II (ARESE II) were analyzed. Grand averages of broadband absorptance among three sets of instruments were compared. Fractional solar absorptances were approx. 0.21-0.22 with the exception of March 3 when two sets of instruments gave values smaller by approx. 0.03-0.04. The robustness of these values was investigated by looking into possible sampling problems with the aid of 500 nm spectral fluxes. Grand averages of 500 nm apparent absorptance cover a wide range of values for these three days, namely from a large positive (approx. 0.011) average for March 3, to a small negative (approximately -0.03) for March 21, to near zero (approx. 0.01) for March 29. We present evidence suggesting that a large part of the discrepancies among the three days is due to the different nature of clouds and their non-uniform sampling. Hence, corrections to the grand average broadband absorptance values may be necessary. However, application of the known correction techniques may be precarious due to the sparsity of collocated flux measurements above and below the clouds. Our analysis leads to the conclusion that only March 29 fulfills all requirements for reliable estimates of cloud absorption, that is, the presence of thick, overcast, homogeneous clouds.
Oreopoulos, L.↗
Radiative Smoothing in Clouds at Transparent and Absorbing Wavelengths
For absorbing and transparent wavelengths, we discuss the effect of horizontal solar radiative fluxes in clouds on the accuracy of a conventional plane-parallel radiative transfer calculations for a single pixel, known as the Independent Pixel Approximation (IPA). We address the question of correlations between horizontal fluxes, IPA accuracies and radiative smoothing. By smoothing we understand a radiative transfer process whereby radiation does not follow the small-scale fluctuations of cloud structure, producing much smoother radiation fields. The scale eta that characterizes this process is called "radiative smoothing scale." We relate radiative smoothing to the photon's horizontal displacement that characterizes a "spot" of reflected light associated with a point source. We generalize the "spot-size" estimate derived for conservative scattering using the diffusion theory to the case of non-conservative scattering. For reflected light, theoretical results are confirmed with numerical simulations. The radiative smoothing scale eta is a critical value where IPA effectively breaks down; for scales smaller than TI, real radiation field are much smoother than their IPA counterparts for the same cloud structure. In addition to the estimate of il for absorbing wavelengths, we show that: (1) with more absorption, the scale break determined by eta in a log-log plot of wavenumber spectra moves towards smaller scales and (2) the smaller eta the flatter the small-scale slope which means less radiative smoothing, thus more accuracy in the IPA reflection.
Marshak, A.↗
Path Radiance Technique for Retrieving Aerosol Optical Thickness over Land
The key issue in retrieving aerosol optical thickness over land from short-wave satellite radiances is to identify and separate the signal due to scattering by a largely transparent aerosol layer from the noise due to reflection by the background surface, where the signal is relatively uniform compared to the highly inhomogeneous surface contribution. Sensitivity studies in aerosol optical thickness retrievals reveal that the apparent reflectance at the top of the atmosphere is very susceptible to the surface reflectance, especially when aerosol optical thickness is small. Uncertainties associated with surface reflectance estimation can greatly amplify the error of the aerosol optical thickness retrieval. To reduce these uncertainties, we have developed a "path radiance" method to retrieve aerosol optical thickness over land by extending the traditional technique that uses the "dark object" approach to extract the aerosol signal. This method uses the signature of the correlation of visible and mid-IR reflectance at the surface, and couples the correlation with the atmospheric effect. We have applied this method to a Landsat TM (Thematic Mapper) image acquired over the Oklahoma Southern Great Plains (SGP) site of DoE's ARM (Department of Energy's Atmospheric Radiation Measurement) program on September 27, 1997, a very clear day (aerosol optical thickness of 0.07 at 0.5 pm) during the first Landsat IOP (Intensive Observation Period). The retrieved mean aerosol optical thickness for TM band 1 at 0.49 pm and band 3 at 0.66 pm agree very well with the ground-based sun-photometer measurements at the ARM site. The ability to retrieve small aerosol optical thickness makes this path radiance technique promising. More importantly, the path radiance is relatively insensitive to surface inhomogeneity. The retrieved mean path radiances in reflectance units have very small standard deviations for both TM blue and red bands. This small variability of path radiance further supports the current aerosol retrieval method.
Tsay, S.-C.↗
Path Radiance Technique for Retrieving Aerosol Optical Thickness Over Land
The key issue in retrieving aerosol optical thickness over land from shortwave satellite radiances is to identify and separate the signal due to scattering by a largely transparent aerosol layer from the noise due to reflection by the background surface, where the signal is relatively uniform compared to the highly inhomogeneous surface contribution. Sensitivity studies in aerosol optical thickness retrievals reveal that the apparent reflectance at the top of the atmosphere is very susceptible to the surface reflectance, especially when aerosol optical thickness is small. Uncertainties associated with surface reflectance estimation can greatly amplify the error of the aerosol optical thickness retrieval. To reduce these uncertainties, we have developed a "path radiance" method to retrieve aerosol optical thickness over land by extending the traditional technique that uses the "dark object" approach to extract the aerosol signal. This method uses the signature of the correlation of visible and mid-IR reflectance at the surface, and couples the correlation with the atmospheric effect. We have applied this method to a Landsat TM (Thematic Mapper) image acquired over the Oklahoma Southern Great Plains (SGP) site of DoE's ARM (Department of Energy's Atmospheric Radiation Measurement) program on September 27,1997, a very clear day (aerosol optical thickness of 0.07 at 0.5 microns) during the first Landsat Intensive Observation Period (IOP). The retrieved mean aerosol optical thickness for TM band 1 at 0.49 micron and band 3 at 0.66 micron agree very well with the ground-based sun-photometer measurements at the ARM site. The ability to retrieve small aerosol optical thickness makes this path radiance technique promising. More importantly, the path radiance is relatively insensitive to surface inhomogeneity. The retrieved mean path radiances in reflectance units have very small standard deviations for both TM blue and red bands. This small variability of path radiance further supports the current aerosol retrieval method.
Tsay, S.-C.↗