Search

The so-called bi-spectral method retrieves cloud optical thickness (τ) and cloud droplet effective radius (re) simultaneously from a pair of cloud reflectance observations, one in a visible or near infrared (VIS/NIR) band and the other in a shortwave-infrared (SWIR) band. A cloudy pixel is usually assumed to be horizontally homogeneous in the retrieval. Ignoring sub-pixel variations of cloud reflectances can lead to a significant bias in the retrieved τ and re. In this study, we use the Taylor expansion of a two-variable function to understand and quantify the impacts of sub-pixel variances of VIS/NIR and SWIR cloud reflectances and their covariance on the τ and re retrievals. This framework takes into account the fact that the retrievals are determined by both VIS/NIR and SWIR band observations in a mutually dependent way. In comparison with previous studies, it provides a more comprehensive understanding of how sub-pixel cloud reflectance variations impact the τ and re retrievals based on the bi-spectral method. In particular, our framework provides a mathematical explanation of how the sub-pixel variation in VIS/NIR band influences the re retrieval and why it can sometimes outweigh the influence of variations in the SWIR band and dominate the error in re retrievals, leading to a potential contribution of positive bias to the re retrieval.

The bispectral method retrieves cloud optical thickness (τ) and cloud droplet effective radius (re) simultaneously from a pair of cloud reflectance observations, one in a visible or near-infrared (VIS/NIR) band and the other in a shortwave infrared (SWIR) band. A cloudy pixel is usually assumed to be horizontally homogeneous in the retrieval. Ignoring subpixel variations of cloud reflectances can lead to a significant bias in the retrieved τ and re. In the literature, the retrievals of τ and re are often assumed to be independent and considered separately when investigating the impact of subpixel cloud reflectance variations on the bispectral method. As a result, the impact on τ is contributed only by the subpixel variation of VIS/NIR band reflectance and the impact on re only by the subpixel variation of SWIR band reflectance. In our new framework, we use the Taylor expansion of a two-variable function to understand and quantify the impacts of subpixel variances of VIS/NIR and SWIR cloud reflectances and their covariance on the τ and re retrievals. This framework takes into account the fact that the retrievals are determined by both VIS/NIR and SWIR band observations in a mutually dependent way. In comparison with previous studies, it provides a more comprehensive understanding of how subpixel cloud reflectance variations impact the τ and re retrievals based on the bispectral method. In particular, our framework provides a mathematical explanation of how the subpixel variation in VIS/NIR band influences the re retrieval and why it can sometimes outweigh the influence of variations in the SWIR band and dominate the error in re retrievals, leading to a potential contribution of positive bias to the re retrieval. We test our framework using synthetic cloud fields from a large-eddy simulation and real observations from Moderate Resolution Imaging Spectroradiometer. The predicted results based on our framework agree very well with the numerical simulations. Our framework can be used to estimate the retrieval uncertainty from subpixel reflectance variations in operational satellite cloud products and to help understand the differences in τ and re retrievals between two instruments.

This study presents preliminary results on the effect of cirrus heterogeneities on top-of-atmosphere (TOA) simulated radiances or reflectances for MODIS channels centered at 0.86, 2.21, 8.56, 11.01 and 12.03 micrometers , and on cloud optical properties retrieved with a research-level optimal estimation method (OEM). Synthetic cirrus cloud fields are generated using a 3D cloud generator (3DCLOUD) and radiances/reflectances are simulated using a 3D radiative transfer code (3DMCPOL). We find significant differences between the heterogeneity effects on either visible and near-infrared (VNIR) or thermal infrared (TIR) radiances. However, when both wavelength ranges are combined, heterogeneity effects are dominated by the VNIR horizontal radiative transport effect. As a result, small optical thicknesses are overestimated and large ones are underestimated. Retrieved effective diameter are found to be slightly affected, contrarily to retrievals using TIR channels only.

Despite its importance as one of the key radiative properties that determines the impact of upper tropospheric clouds on the radiation balance, ice cloud optical thickness (IOT) has proven to be one of the more challenging properties to retrieve from space-based remote sensing measurements. In particular, optically thin upper tropospheric ice clouds (cirrus) have been especially challenging due to their tenuous nature, extensive spatial scales, and complex particle shapes and light scattering characteristics. The lack of independent validation motivates the investigation presented in this paper, wherein systematic biases between MODIS Collection 5 (C5) and CALIOP Version 3 (V3) unconstrained retrievals of tenuous IOT (< 3) are examined using a month of collocated A-Train observations. An initial comparison revealed a factor of two bias between the MODIS and CALIOP IOT retrievals. This bias is investigated using an infrared (IR) radiative closure approach that compares both products with MODIS IR cirrus retrievals developed for this assessment. The analysis finds that both the MODIS C5 and the unconstrained CALIOP V3 retrievals are biased (high and low, respectively) relative to the IR IOT retrievals. Based on this finding, the MODIS and CALIOP algorithms are investigated with the goal of explaining and minimizing the biases relative to the IR. For MODIS we find that the assumed ice single scattering properties used for the C5 retrievals are not consistent with the mean IR COT distribution. The C5 ice scattering database results in the asymmetry parameter (g) varying as a function of effective radius with mean values that are too large. The MODIS retrievals have been brought into agreement with the IR by adopting a new ice scattering model for Collection 6 (C6) consisting of a modified gamma distribution comprised of a single habit (severely roughened aggregated columns); the C6 ice cloud optical property models have a constant g approx. = 0.75 in the mid-visible spectrum, 5-15% smaller than C5. For CALIOP, the assumed lidar ratio for unconstrained retrievals is fixed at 25 sr for the V3 data products.This value is found to be inconsistent with the constrained (predominantly nighttime) CALIOP retrievals. An experimental data set was produced using a modified lidar ratio of 32 sr for the unconstrained retrievals (an increase of 28%), selected to provide consistency with the constrained V3 results. These modifications greatly improve the agreement with the IR and provide consistency between the MODIS and CALIOP products. Based on these results the recently released MODIS C6 optical products use the single habit distribution given above, while the upcoming CALIOP V4 unconstrained algorithm will use higher lidar ratios for unconstrained retrievals.

The Collection 6 (C6) MODIS (Moderate Resolution Imaging Spectroradiometer) land and atmosphere data sets are scheduled for release in 2014. C6 contains significant revisions of the calibration approach to account for sensor aging. This analysis documents the presence of systematic temporal trends in the visible and near-infrared (500 m) bands of the Collection 5 (C5) MODIS Terra and, to lesser extent, in MODIS Aqua geophysical data sets. Sensor degradation is largest in the blue band (B3) of the MODIS sensor on Terra and decreases with wavelength. Calibration degradation causes negative global trends in multiple MODIS C5 products including the dark target algorithm's aerosol optical depth over land and Ångstrom exponent over the ocean, global liquid water and ice cloud optical thickness, as well as surface reflectance and vegetation indices, including the normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI). As the C5 production will be maintained for another year in parallel with C6, one objective of this paper is to raise awareness of the calibration-related trends for the broad MODIS user community. The new C6 calibration approach removes major calibrations trends in the Level 1B (L1B) data. This paper also introduces an enhanced C6C calibration of the MODIS data set which includes an additional polarization correction (PC) to compensate for the increased polarization sensitivity of MODIS Terra since about 2007, as well as detrending and Terra- Aqua cross-calibration over quasi-stable desert calibration sites. The PC algorithm, developed by the MODIS ocean biology processing group (OBPG), removes residual scan angle, mirror side and seasonal biases from aerosol and surface reflectance (SR) records along with spectral distortions of SR. Using the multiangle implementation of atmospheric correction (MAIAC) algorithm over deserts, we have also developed a detrending and cross-calibration method which removes residual decadal trends on the order of several tenths of 1% of the top-of-atmosphere (TOA) reflectance in the visible and near-infrared MODIS bands B1-B4, and provides a good consistency between the two MODIS sensors. MAIAC analysis over the southern USA shows that the C6C approach removed an additional negative decadal trend of Terra (Delta)NDVI approx.0.01 as compared to Aqua data. This change is particularly important for analysis of vegetation dynamics and trends in the tropics, e.g., Amazon rainforest, where the morning orbit of Terra provides considerably more cloud-free observations compared to the afternoon Aqua measurements.

The role of clouds remains the largest uncertainty in climate projections. They influence solar and thermal radiative transfer and the earth's water cycle. Therefore, there is an urgent need for accurate cloud observations to validate climate models and to monitor climate change. Passive satellite imagers measuring radiation at visible to thermal infrared (IR) wavelengths provide a wealth of information on cloud properties. Among others, the cloud top height (CTH) - a crucial parameter to estimate the thermal cloud radiative forcing - can be retrieved. In this paper we investigate the skill of ten current retrieval algorithms to estimate the CTH using observations from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard Meteosat Second Generation (MSG). In the first part we compare ten SEVIRI cloud top pressure (CTP) data sets with each other. The SEVIRI algorithms catch the latitudinal variation of the CTP in a similar way. The agreement is better in the extratropics than in the tropics. In the tropics multi-layer clouds and thin cirrus layers complicate the CTP retrieval, whereas a good agreement among the algorithms is found for trade wind cumulus, marine stratocumulus and the optically thick cores of the deep convective system. In the second part of the paper the SEVIRI retrievals are compared to CTH observations from the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) and Cloud Profiling Radar (CPR) instruments. It is important to note that the different measurement techniques cause differences in the retrieved CTH data. SEVIRI measures a radiatively effective CTH, while the CTH of the active instruments is derived from the return time of the emitted radar or lidar signal. Therefore, some systematic differences are expected. On average the CTHs detected by the SEVIRI algorithms are 1.0 to 2.5 kilometers lower than CALIOP observations, and the correlation coefficients between the SEVIRI and the CALIOP data sets range between 0.77 and 0.90. The average CTHs derived by the SEVIRI algorithms are closer to the CPR measurements than to CALIOP measurements. The biases between SEVIRI and CPR retrievals range from −0.8 kilometers to 0.6 kilometers. The correlation coefficients of CPR and SEVIRI observations vary between 0.82 and 0.89. To discuss the origin of the CTH deviation, we investigate three cloud categories: optically thin and thick single layer as well as multi-layer clouds. For optically thick clouds the correlation coefficients between the SEVIRI and the reference data sets are usually above 0.95. For optically thin single layer clouds the correlation coefficients are still above 0.92. For this cloud category the SEVIRI algorithms yield CTHs that are lower than CALIOP and similar to CPR observations. Most challenging are the multi-layer clouds, where the correlation coefficients are for most algorithms between 0.6 and 0.8. Finally, we evaluate the performance of the SEVIRI retrievals for boundary layer clouds. While the CTH retrieval for this cloud type is relatively accurate, there are still considerable differences between the algorithms. These are related to the uncertainties and limited vertical resolution of the assumed temperature profiles in combination with the presence of temperature inversions, which lead to ambiguities in the CTH retrieval. Alternative approaches for the CTH retrieval of low clouds are discussed.

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.

In this paper we describe a general procedure for calculating synthetic sensor radiances from variable output from a global atmospheric forecast model. In order to take proper account of the discrepancies between model resolution and sensor footprint, the algorithm takes explicit account of the model subgrid variability, in particular its description of the probability density function of total water (vapor and cloud condensate.) The simulated sensor radiances are then substituted into an operational remote sensing algorithm processing chain to produce a variety of remote sensing products that would normally be produced from actual sensor output. This output can then be used for a wide variety of purposes such as model parameter verification, remote sensing algorithm validation, testing of new retrieval methods and future sensor studies.We show a specific implementation using the GEOS-5 model, the MODIS instrument and the MODIS Adaptive Processing System (MODAPS) Data Collection 5.1 operational remote sensing cloud algorithm processing chain (including the cloud mask, cloud top properties and cloud optical and microphysical properties products). We focus on clouds because they are very important to model development and improvement.

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.

The seasonal and spatial varability of Aerosol Above Cloud (AAC) properties are derived from passive satellite data for the year 2008. A significant amount of aerosols are transported above liquid water clouds on the global scale. For particles in the fine mode (i.e., radius smaller than 0.3 m), including both clear sky and AAC retrievals increases the global mean aerosol optical thickness by 25(+/- 6%). The two main regions with man-made AAC are the tropical Southeast Atlantic, for biomass burning aerosols, and the North Pacific, mainly for pollutants. Man-made AAC are also detected over the Arctic during the spring. Mineral dust particles are detected above clouds within the so-called dust belt region (5-40 N). AAC may cause a warming effect and bias the retrieval of the cloud properties. This study will then help to better quantify the impacts of aerosols on clouds and climate.

Clouds cover about 70% of the Earth's surface and play a dominant role in the energy and water cycle of our planet. Only satellite observations provide a continuous survey of the state of the atmosphere over the whole globe and across the wide range of spatial and temporal scales that comprise weather and climate variability. Satellite cloud data records now exceed more than 25 years in length. However, climatologies compiled from different satellite datasets can exhibit systematic biases. Questions therefore arise as to the accuracy and limitations of the various sensors. The Global Energy and Water cycle Experiment (GEWEX) Cloud Assessment, initiated in 2005 by the GEWEX Radiation Panel, provided the first coordinated intercomparison of publically available, standard global cloud products (gridded, monthly statistics) retrieved from measurements of multi-spectral imagers (some with multiangle view and polarization capabilities), IR sounders and lidar. Cloud properties under study include cloud amount, cloud height (in terms of pressure, temperature or altitude), cloud radiative properties (optical depth or emissivity), cloud thermodynamic phase and bulk microphysical properties (effective particle size and water path). Differences in average cloud properties, especially in the amount of high-level clouds, are mostly explained by the inherent instrument measurement capability for detecting and/or identifying optically thin cirrus, especially when overlying low-level clouds. The study of long-term variations with these datasets requires consideration of many factors. A monthly, gridded database, in common format, facilitates further assessments, climate studies and the evaluation of climate models.

Moderate Resolution Imaging Spectroradiometer (MODIS) retrievals of optical thickness and effective particle radius for liquid water and ice phase clouds employ a well-known VNIR/ SWIR solar reflectance technique. For this type of algorithm, we evaluate the quantitative uncertainty in simultaneous retrievals of these two cloud parameters to pixel-level radiometric calibration estimates and other fundamental (and tractable) error sources.

Satellite observations show that ice cloud effective radius (r(sub e)) increases with ice water content (IWC) but decreases with aerosol optical thickness (AOT). Using least-squares fitting to the observed data, we obtain an analytical formula to describe the variations of r(sub e) with IWC and AOT for several regions with distinct characteristics of r(sub e) -IWC-AOT relationships. As IWC directly relates to convective strength and AOT represents aerosol loading, our empirical formula provides a means to quantify the relative roles of dynamics and aerosols in controlling r(sub e) in different geographical regions, and to establish a framework for parameterization of aerosol effects on r(sub e) in climate models.

Large fine mode (sub-micron radius) dominated aerosols in size distributions retrieved from AERONET have been observed after fog or low-altitude cloud dissipation events. These column-integrated size distributions have been obtained at several sites in many regions of the world, typically after evaporation of low altitude cloud such as stratocumulus or fog. Retrievals with cloud processed aerosol are sometimes bimodal in the accumulation mode with the larger size mode often approx.0.4 - 0.5 microns radius (volume distribution); the smaller mode typically approx.0.12 to aprrox.0.20 microns may be interstitial aerosol that were not modified by incorporation in droplets and/or aerosol that are less hygroscopic in nature. Bimodal accumulation mode size distributions have often been observed from in situ measurements of aerosols that have interacted with clouds, and AERONET size distribution retrievals made after dissipation of cloud or fog are in good agreement with particle sizes measured by in situ techniques for cloud-processed aerosols. Aerosols of this type and large size range (in lower concentrations) may also be formed by cloud processing in partly cloudy conditions and may contribute to the shoulder of larger size particles in the accumulation mode retrievals, especially in regions where sulfate and other soluble aerosol are a significant component of the total aerosol composition. Observed trends of increasing aerosol optical depth (AOD) as fine mode radius increased suggests higher AOD in the near cloud environment and therefore greater aerosol direct radiative forcing than typically obtained from remote sensing, due to bias towards sampling at low cloud fraction.

Haywood et al. (2004) show that an aerosol layer above a cloud can cause a bias in the retrieved cloud optical thickness and effective radius. Monitoring for this potential bias is difficult because space ]based passive remote sensing cannot unambiguously detect or characterize aerosol above cloud. We show that cloud retrievals from aircraft measurements above cloud and below an overlying aerosol layer are a means to test this bias. The data were collected during the Intercontinental Chemical Transport Experiment (INTEX-A) study based out of Portsmouth, New Hampshire, United States, above extensive, marine stratus cloud banks affected by industrial outflow. Solar Spectral Flux Radiometer (SSFR) irradiance measurements taken along a lower level flight leg above cloud and below aerosol were unaffected by the overlying aerosol. Along upper level flight legs, the irradiance reflected from cloud top was transmitted through an aerosol layer. We compare SSFR cloud retrievals from below ]aerosol legs to satellite retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) in order to detect an aerosol ]induced bias. In regions of small variation in cloud properties, we find that SSFR and MODIS-retrieved cloud optical thickness compares within the uncertainty range for each instrument while SSFR effective radius tend to be smaller than MODIS values (by 1-2 microns) and at the low end of MODIS uncertainty estimates. In regions of large variation in cloud properties, differences in SSFR and MODIS ]retrieved cloud optical thickness and effective radius can reach values of 10 and 10 microns, respectively. We include aerosols in forward modeling to test the sensitivity of SSFR cloud retrievals to overlying aerosol layers. We find an overlying absorbing aerosol layer biases SSFR cloud retrievals to smaller effective radii and optical thickness while nonabsorbing aerosols had no impact.

The detection of mUltiple cloud layers using satellite observations is important for retrieval algorithms as well as climate applications. In this paper, we describe a relatively simple algorithm to detect multiple cloud layers and distinguish them from vertically-extended clouds. The algorithm can be applied to coincident passive sensors that derive both cloud-top pressure from the thermal infrared observations and an estimate of solar photon pathlength from UV, visible, or near-IR measurements. Here, we use data from the A-train afternoon constellation of satellites: cloud-top pressure, cloud optical thickness, the multi-layer flag from the Aqua MODerate-resolution Imaging Spectroradiometer (MODIS) and the optical centroid cloud pressure from the Aura Ozone Monitoring Instrument (OMI). For the first time, we use data from the CloudSat radar to evaluate the results of a multi-layer cloud detection scheme. The cloud classification algorithms applied with different passive sensor configurations compare well with each other as well as with data from CloudSat. We compute monthly mean fractions of pixels containing multi-layer and vertically-extended clouds for January and July 2007 at the OMI spatial resolution (l2kmx24km at nadir) and at the 5kmx5km MODIS resolution used for infrared cloud retrievals. There are seasonal variations in the spatial distribution of the different cloud types. The fraction of cloudy pixels containing distinct multi-layer cloud is a strong function of the pixel size. Globally averaged, these fractions are approximately 20% and 10% for OMI and MODIS, respectively. These fractions may be significantly higher or lower depending upon location. There is a much smaller resolution dependence for fractions of pixels containing vertically-extended clouds (approx.20% for OMI and slightly less for MODIS globally), suggesting larger spatial scales for these clouds. We also find higher fractions of vertically-extended clouds over land as compared with ocean, particularly in the tropics and summer hemisphere.

In situ measurements of ice crystal concentrations and sizes made with aircraft instrumentation over the past two decades have often indicated the presence of numerous relatively small (< 50 m diameter) crystals in cirrus clouds. Further, these measurements frequently indicate that small crystals account for a large fraction of the extinction in cirrus clouds. The fact that the instruments used to make these measurements, such as the Forward Scattering Spectrometer Probe (FSSP) and the Cloud Aerosol Spectrometer (CAS), ingest ice crystals into the sample volume through inlets has led to suspicion that the indications of numerous small ]crystals could be artifacts of large ]crystal shattering on the instrument inlets. We present new aircraft measurements in anvil cirrus sampled during the Tropical Composition, Cloud, and Climate Coupling (TC4) campaign with the 2 ] Dimensional Stereo (2D ]S) probe, which detects particles as small as 10 m. The 2D ]S has detector "arms" instead of an inlet tube. Since the 2D ]S probe surfaces are much further from the sample volume than is the case for the instruments with inlets, it is expected that 2D ]S will be less susceptible to shattering artifacts. In addition, particle inter ]arrival times are used to identify and remove shattering artifacts that occur even with the 2D ]S probe. The number of shattering artifacts identified by the 2D ]S interarrival time analysis ranges from a negligible contribution to an order of magnitude or more enhancement in apparent ice concentration over the natural ice concentration, depending on the abundance of large crystals and the natural small ]crystal concentration. The 2D ]S measurements in tropical anvil cirrus suggest that natural small ]crystal concentrations are typically one to two orders of magnitude lower than those inferred from CAS. The strong correlation between the CAS/2D ]S ratio of small ]crystal concentrations and large ]crystal concentration suggests that the discrepancy is likely caused by shattering of large crystals on the CAS inlet. We argue that past measurements with CAS in cirrus with large crystals present may contain errors due to crystal shattering, and past conclusions derived from these measurements may need to be revisited. Further, we present correlations between CAS spurious concentration and 2D ]S large ]crystal mass from spatially uniform anvil cirrus sampling periods as an approximate guide for estimating quantitative impact of large ]crystal shattering on CAS concentrations in previous datasets. We use radiative transfer calculations to demonstrate that in the maritime anvil cirrus sampled during TC4, small crystals indicated by 2D ]S contribute relatively little cloud extinction, radiative forcing, or radiative heating in the anvils, regardless of anvil age or vertical location in the clouds. While 2D ]S ice concentrations in fresh anvil cirrus may often exceed 1 cm.3, and are observed to exceed 10 cm.3 in turrets, they are typically ~0.1 cm.3 and rarely exceed 1 cm.3 (<1.4% of the time) in aged anvil cirrus. We hypothesize that isolated occurrences of higher ice concentrations in aged anvil cirrus may be caused by ice nucleation driven by either small ]scale convection or gravity waves. It appears that the numerous small crystals detrained from convective updrafts do not persist in the anvil cirrus sampled during TC ]4.

This study compares cirrus particle effective radius retrieved by a CALIPSO-like method with two similar methods using MODIS, MODI Airborne Simulator (MAS), and GOES imagery. The CALIPSO-like method uses lidar measurements coupled with the split-window technique that uses the infrared spectral information contained at the 8.65-micrometer, 11.15-micrometer and 12.05-micrometer bands to infer the microphysical properties of cirrus clouds. The two other methods, sing passive remote sensing at visible and infrared wavelengths, are the operational MODIS cloud products (referred to by its archival product identifier MOD06 for MODIS Terra) and MODIS retrievals performed by the CERES team at LaRC (Langley Research Center) in support of CERES algorithms; the two algorithms will be referred to as MOD06- and LaRC-method, respectively. The three techniques are compared at two different latitudes: (i) the mid-latitude ice clouds study uses 18 days of observations at the Palaiseau ground-based site in France (SIRTA: Site Instrumental de Recherche par Teledetection Atmospherique) including a ground-based 532 nm lidar and the Moderate Resolution Imaging Spectrometer (MODIS) overpasses on the Terra Platform, (ii) the tropical ice clouds study uses 14 different flight legs of observations collected in Florida, during the intensive field experiment CRYSTAL-FACE (Cirrus Regional Study of Tropical Anvils and cirrus Layers-Florida Area Cirrus Experiment), including the airborne Cloud Physics Lidar (CPL) and the MAS. The comparison of the three methods gives consistent results for the particle effective radius and the optical thickness, but discrepancies in cloud detection and altitudes. The study confirms the value of an active remote-sensing method (CALIPSO-like) for the study of sub-visible ice clouds, in both mid-latitudes and tropics. Nevertheless, this method is not reliable in optically very thick tropical ice clouds.

Moderate Resolution Imaging Spectroradiometer (MODIS) retrievals of cloud optical thickness and effective particle radius employ well-known solar reflectance techniques using pre-calculated reflectance look-up tables. We evaluate the quantitative uncertainty in simultaneous retrievals of cloud optical thickness and particle size for this type of algorithm. The technique uses sensitivity calculations derived from the reflectance look-up tables, coupled with estimates for the effect of various error terms on the uncertainty in inferring the actual cloud-top reflectance. The error terms include the effects of instrument calibration, surface spectral albedo, and atmospheric corrections on both water and ice cloud retrievals. Because particle shapes in ice clouds are highly variable, the effect of particle shape is analyzed separately with a more approximate method. Results will deal exclusively with pixel-level uncertainties associated with plane-parallel clouds; real-world radiative departures from a plane-parallel model are an additional consideration. While we demonstrate the uncertainty technique with operational 1 km MODIS retrievals from the Terra and Aqua satellite platforms, the technique is applicable to any reflectance-based satellite- or air-borne sensor retrieval using similar spectral channels.

The Moderate Resolution Imaging Spectroradiometer (MODE) was developed by NASA and launched onboard the Terra spacecraft on December 18,1999 and the Aqua spacecraft on April 26,2002. MODIS scans a swath width sufficient to provide nearly complete global coverage every two days from each polar-orbiting, sun-synchronous, platform at an altitude of 705 km, and provides images in 36 spectral bands between 0.415 and 14.235 pm with spatial resolutions of 250 m (2 bands), 500 m (5 bands) and 1000 m (29 bands). In this paper, we describe the radiative properties of clouds as currently determined from satellites (cloud fraction, optical thickness, cloud top pressure, and cloud effective radius), and highlight the global and regional cloud microphysical properties currently available for assessing climate variability and forcing. These include the latitudinal distribution of cloud optical and radiative properties of both liquid water and ice clouds, as well as joint histograms of cloud optical thickness and effective radius for selected geographical locations around the globe.

In order to perform satellite retrievals of cloud properties, it is important to account for the effect of the above-cloud atmosphere on the observations. The solar bands used in the operational MODIS Terra and Aqua cloud optical and microphysical algorithms (visible, NIR, and SWIR spectral windows) are primarily affected by water vapor, and to a lesser extent by well-mixed gases. For water vapor, the above-cloud column amount, or precipitable water, provides adequate information for an atmospheric correction; details of the vertical vapor distribution are not typically necessary for the level of correction required. Cloud-top pressure has a secondary effect due to pressure broadening influences. For well- mixed gases, cloud-top pressure is also required for estimates of above-cloud abundances. We present a method for obtaining above-cloud precipitable water over dark Ocean surfaces using the MODIS 0.94 pm vapor absorption band. The retrieval includes an iterative procedure for establishing cloud-top temperature and pressure, and is useful for both single layer water and ice clouds. Knowledge of cloud thermodynamic phase is fundamental in retrieving cloud optical and microphysical properties. However, in cases of optically thin cirrus overlapping lower water clouds, the concept of a single unique phase is ill- defined and depends, at least, on the spectral region of interest. We will present a method for multi-layer and multi-phase cloud detection which uses above-cloud precipitable water retrievals along with several existing MODIS operational cloud products (cloud-top pressure derived from a C02 slicing algorithm, IR and SWIR phase retrievals). Results are catagorized by whether the radiative signature in the MODIS solar bands is primarily that of a water cloud with ice cloud contamination, or visa-versa. Examples in polar and mid-latitude regions will be shown.

Spectral land surface albedo is an important parameter for describing the radiative properties of the Earth. Accordingly it reflects the consequences of natural and human interactions, such as anthropogenic, meteorological, and phenological effects, on global and local climatological trends. Consequently, albedos are integral parts in a variety of research areas, such as general circulation models (GCMs), energy balance studies, modeling of land use and land use change, and biophysical, oceanographic, and meteorological studies. The availability of global albedo data over a large range of spectral channels and at high spatial resolution has dramatically improved with the launch of the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument aboard NASA s Earth Observing System (EOS) Terra spacecraft in December 1999. However, lack of spatial and temporal coverage due to cloud and snow effects can preclude utilization of official products in production and research studies. We report on a technique used to fill incomplete MOD43 albedo data sets with the intention of providing complete value-added maps. The technique is influenced by the phenological concept that within a certain area, a pixel s ecosystem class should exhibit similar growth cycle events over the same time period. The shape of an area s phenological temporal curve can be imposed upon existing pixel-level data to fill missing temporal points. The methodology will be reviewed by showcasing 2001 global and regional results of complete albedo and NDVl data sets.

MODIS observations from the NASA EOS Terra spacecraft (launched in December 1999, 1030 local time equatorial crossing) have provided a unique data set of Earth observations. With the launch of the NASA Aqua spacecraft in May 2002 (1330 local time), two MODIS daytime (sunlit) and nighttime observations are now available in a 24 hour period, allowing for some measure of diurnal variability. We report on an initial analysis of several operational global (Level-3) cloud products from the two platforms. The MODIS atmosphere Level-3 products, which include clear-sky and aerosol products in addition to cloud products, are available as three separate files providing daily, eight-day, and monthly aggregations; each temporal aggregation is spatially aggregated to a 1 degree grid. The files contain approximately 600 statisitical datasets (from simple means and standard deviations to 1 - and 2-dimensional histograms). Operational cloud products include detection (cloud fraction), cloud-top properties, and daytimeonly cloud optical thickness and particle effective radius for both water and ice clouds. We will compare example global Terra and Aqua cloud fraction, optical thickness, and effective radius aggregations.

The Moderate Resolution Imaging Spectroradiometer (MODIS), developed as part of the Earth Observing System (EOS) and launched on Terra in December 1999 and Aqua in May 2002, is designed to meet the scientific needs for satellite remote sensing of clouds, aerosols, water vapor, and land and ocean surface properties. During the CRYSTAL-FACE experiment, numerous aircraft coordinated both in situ and remote sensing observations with the Terra and Aqua spacecraft. In this paper we will emphasize the optical, microphysical, and physical properties of both liquid water and ice clouds obtained from an analysis of the satellite observations over Florida and the Gulf of Mexico during July 2002. We will present the frequency distribution of liquid water and ice cloud microphysical properties throughout the region, separating the results over land and ocean. Probability distributions of effective radius and cloud optical thickness will also be shown.

The conventional Lorenz-Mie formalism is extended to the scattering process associated with a coated sphere embedded in an absorbing medium. It is shown that apparent and inherent scattering cross sections of a scattering particle, which are identical in the case of transparent host medium, are different if the host medium is absorptive. Here the inherent single-scattering properties are derived from the near-field information whereas the corresponding apparent counterparts are derived from the far-field asymptotic form of the scattered wave with scaling of host absorption that is assumed to be in an exponential form. The formality extinction and scattering efficiencies defined in the same manner as in the conventional sense can be unbounded. For a nonabsorptive particle embedded in an absorbing medium, the effect of host absorption on the phase matrix elements associated with polarization is significant. This effect, however, is largely reduced for strongly absorptive particles such as soot. For soot particles coated with water, the impurity can substantially reduce the single-scattering albedo of the particle if the size parameter is small. For water-coating soot and hollow ice spheres, it is shown that the phase matrix elements -P(sub 12)/P(sub 11) and P(sub 33)/P(sub 11) are unique if the shell is thin, as compared with the case for thick shell. Furthermore, the radiative transfer equation regarding a multidisperse particle system in an absorbing medium is discussed. It is illustrated that the conventional computation algorithms can be applied to solve the multiple scattering process if the scaled apparent single-scattering properties are applied.

The Moderate Resolution Imaging Spectroradiometer (MODIS) is one of five instruments aboard the Terra Earth Observing System platform launched in December 1999. With 36 spectral bands from the visible through the infrared, and spatial resolution from 250m to 1km, the instrument provides an unprecedented opportunity for global cloud studies. A comprehensive set of remote sensing algorithms for cloud masking and retrieval of cloud physical and optical properties have been developed by members of the MODIS atmosphere team. The archived products from these algorithms have applications in climate change studies, climate modeling, numerical weather prediction, as well as fundamental atmospheric research. In addition to an extensive cloud mask, products include cloud top physical parameters (temperature, pressure, emissivity), cloud phase, cloud optical parameters (optical depth, effective particle radius, water path), visible cirrus reflectance, a contrail flag, and other derived parameters. All products are archived into two categories: pixel-level retrievals at a 1 km or 5 km spatial resolution at nadir (referred to as Level-2 products) and 1 degree global gridded statistics (Level-3 products). An overview of the MODIS atmosphere algorithms and products, their status, validation activities, and early level-2 and -3 results will be presented.

The Moderate Resolution Imaging Spectroradiometer (MODIS) provides an unprecedented opportunity for global cloud studies with 36 spectral bands from the visible through the infrared, and spatial resolution from 250 m to 1 km at nadir. In particular, all solar window bands useful for simultaneous retrievals of cloud optical thickness and particle size (0.67, 0.86, 1.2, 1.6, 2.1, and 3.7 micron bands) are now available on a single satellite instrument/platform for the first time. An operational algorithm for the retrieval of these optical and cloud physical properties (including water path) have been developed for both liquid and ice phase clouds. The product is archived into two categories: pixel-level retrievals at 1 km spatial resolution (referred to as a Level-2 product) and global gridded statistics (Level-3 product). An overview of the MODIS cloud retrieval algorithm and early level-2 and -3 results will be presented. A number of MODIS cloud validation activities are being planned, including the recent Southern Africa Regional Science Initiative 2000 (SAFARI-2000) dry season campaign conducted in August/September 2000. The later part of the experiment concentrated on MODIS validation in the Namibian stratocumulus regime off the southwest coast of Africa. Early retrieval results from this regime will be discussed.

Photon transport in plane-parallel, vertically inhomogeneous clouds is investigated and applied to cloud remote sensing techniques that use solar reflectance or transmittance measurements for retrieving droplet effective radius. Transport is couched in terms of weighting functions which approximate the relative contribution of individual layers to the overall retrieval. Two vertical weightings are investigated, including one based on the average number of scatterings encountered by reflected and transmitted photons in any given layer. A simpler vertical weighting based on the maximum penetration of reflected photons proves useful for solar reflectance measurements. These weighting functions are highly dependent on droplet absorption and solar/viewing geometry. A superposition technique, using adding/doubling radiative transfer procedures, is derived to accurately determine both weightings, avoiding time consuming Monte Carlo methods. Superposition calculations are made for a variety of geometries and cloud models, and selected results are compared with Monte Carlo calculations. Effective radius retrievals from modeled vertically inhomogeneous liquid water clouds are then made using the standard near-infrared bands, and compared with size estimates based on the proposed weighting functions. Agreement between the two methods is generally within several tenths of a micrometer, much better than expected retrieval accuracy. Though the emphasis is on photon transport in clouds, the derived weightings can be applied to any multiple scattering plane-parallel radiative transfer problem, including arbitrary combinations of cloud, aerosol, and gas layers.

Cloud optical thickness and effective radius retrievals from solar reflectance measurements traditionally depend on a combination of spectral channels that are absorbing and non-absorbing for liquid water droplets. Reflectances in non-absorbing channels (e.g., 0.67, 0.86 micrometer bands) are largely dependent on cloud optical thickness, while longer wavelength absorbing channels (1.6, 2.1, and 3.7 micrometer window bands) provide cloud particle size information. Retrievals are complicated by the presence of an underlying ice/snow surface. At the shorter wavelengths, sea ice is both bright and highly variable, significantly increasing cloud retrieval uncertainty. However, reflectances at the longer wavelengths are relatively small and may be comparable to that of dark open water. Sea ice spectral albedos derived from Cloud Absorption Radiometer (CAR) measurements during April 1992 and June 1995 Arctic field deployments are used to illustrate these statements. A modification to the traditional retrieval technique is devised. The new algorithm uses a combination of absorbing spectral channels for which the snow/ice albedo is relatively small. Using this approach, preliminary retrievals have been made with the MODIS Airborne Simulator (MAS) imager flown aboard the NASA ER-2 during FIRE-ACE. Data from coordinated ER-2 and University of Washington CV-580 aircraft observations of liquid water stratus clouds on June 3 and June 6, 1998 have been examined. Size retrievals are compared with in situ cloud profile measurements of effective radius made with the CV-580 PMS FSSP probe, and optical thickness retrievals are compared with extinction profiles derived from the Gerber Scientific "g-meter" probe. MAS retrievals are shown to be in good agreement with the in situ measurements.

Cloud optical thickness and effective radius retrievals from solar reflectance measurements are traditionally implemented using a combination of spectral channels that are absorbing and non-absorbing for water particles. Reflectances in non-absorbing channels (e.g., 0.67, 0.86, 1.2 micron spectral window bands) are largely dependent on cloud optical thickness, while longer wavelength absorbing channels (1.6, 2. 1, and 3.7 micron window bands) provide cloud particle size information. Cloud retrievals over ice and snow surfaces present serious difficulties. At the shorter wavelengths, ice is bright and highly variable, both characteristics acting to significantly increase cloud retrieval uncertainty. In contrast, reflectances at the longer wavelengths are relatively small and may be comparable to that of dark open water. A modification to the traditional cloud retrieval technique is devised. The new algorithm uses only a combination of absorbing spectral channels for which the snow/ice albedo is relatively small. Using this approach, retrievals have been made with the MODIS Airborne Simulator (MAS) imager flown aboard the NASA ER-2 from May - June 1998 during the Arctic FIRE-ACE field deployment. Data from several coordinated ER-2 and University of Washington CV-580 in situ aircraft observations of liquid water stratus clouds are examined. MAS retrievals of optical thickness, droplet effective radius, and liquid water path are shown to be in good agreement with the in situ measurements. The initial success of the technique has implications for future operational satellite cloud retrieval algorithms in polar and wintertime regions.

Photon transport in a multiple scattering medium is critically dependent on scattering statistics, in particular the average number of scatterings. A superposition technique is derived to accurately determine the average number of scatterings encountered by reflected and transmitted photons within arbitrary layers in plane-parallel, vertically inhomogeneous clouds. As expected, the resulting scattering number profiles are highly dependent on cloud particle absorption and solar/viewing geometry. The technique uses efficient adding and doubling radiative transfer procedures, avoiding traditional time-intensive Monte Carlo methods. Derived superposition formulae are applied to a variety of geometries and cloud models, and selected results are compared with Monte Carlo calculations. Cloud remote sensing techniques that use solar reflectance or transmittance measurements generally assume a homogeneous plane-parallel cloud structure. The scales over which this assumption is relevant, in both the vertical and horizontal, can be obtained from the superposition calculations. Though the emphasis is on photon transport in clouds, the derived technique is applicable to any scattering plane-parallel radiative transfer problem, including arbitrary combinations of cloud, aerosol, and gas layers in the atmosphere.