Clouds cover about 70% of 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 entire globe and across the wide range of spatial and temporal scales that compose weather and climate variability. Satellite cloud data records now exceed more than 25 years; however, climate data records must be compiled from different satellite datasets and can exhibit systematic biases. Questions therefore arise as to the accuracy and limitations of the various sensors and retrieval methods. The Global Energy and Water Cycle Experiment (GEWEX) Cloud Assessment, initiated in 2005 by the GEWEX Radiation Panel (GEWEX Data and Assessment Panel since 2011), provides the first coordinated intercomparison of publicly available, standard global cloud products (gridded monthly statistics) retrieved from measurements of multispectral 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 thermodynamic phase, and cloud radiative and bulk microphysical properties (optical depth or emissivity, effective particle radius, 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. The monthly gridded database presented here facilitates further assessments, climate studies, and the evaluation of climate models.
The effects of dust storms on cloud properties and Radiative Forcing (RF) are analyzed over Northwestern China from April 2001 to June 2004 using data collected by the MODerate Resolution Imaging ...Spectroradiometer (MODIS) and Clouds and the Earth's Radiant Energy System (CERES) instruments on the Aqua and Terra satellites. On average, ice cloud effective particle diameter, optical depth and ice water path of cirrus clouds under dust polluted conditions are 11%, 32.8%, and 42% less, respectively, than those derived from ice clouds in dust‐free atmospheric environments. Due to changes in cloud microphysics, the instantaneous net RF is increased from −161.6 W/m2 for dust‐free clouds to −118.6 W/m2 for dust‐contaminated clouds.
Retrieval of ice cloud properties using IR measurements has a distinct advantage over the visible and near‐IR techniques by providing consistent monitoring regardless of solar illumination ...conditions. Historically, the IR bands at 3.7, 6.7, 11.0, and 12.0 µm have been used to infer ice cloud parameters by various methods, but the reliable retrieval of ice cloud optical depth τ is limited to nonopaque cirrus with τ < 8. The Ice Cloud Optical Depth from Infrared using a Neural network (ICODIN) method is developed in this paper by training Moderate Resolution Imaging Spectroradiometer (MODIS) radiances at 3.7, 6.7, 11.0, and 12.0 µm against CloudSat‐estimated τ during the nighttime using 2 months of matched global data from 2007. An independent data set comprising observations from the same 2 months of 2008 was used to validate the ICODIN. One 4‐channel and three 3‐channel versions of the ICODIN were tested. The training and validation results show that IR channels can be used to estimate ice cloud τ up to 150 with correlations above 78% and 69% for all clouds and only opaque ice clouds, respectively. However, τ for the deepest clouds is still underestimated in many instances. The corresponding RMS differences relative to CloudSat are ~100 and ~72%. If the opaque clouds are properly identified with the IR methods, the RMS differences in the retrieved optical depths are ~62%. The 3.7 µm channel appears to be most sensitive to optical depth changes but is constrained by poor precision at low temperatures. A method for estimating total optical depth is explored for estimation of cloud water path in the future. Factors affecting the uncertainties and potential improvements are discussed. With improved techniques for discriminating between opaque and semitransparent ice clouds, the method can ultimately improve cloud property monitoring over the entire diurnal cycle.
Key Points
Information about nocturnal opaque ice cloud optical depth available in multispectral infrared data
Neural network algorithm developed to retrieve ice cloud optical depth from nighttime imager data
Total cloud optical depth can be estimated from neural net ice COD using a parameterization
Present‐day shortcomings in the representation of upper tropospheric ice clouds in general circulation models (GCMs) lead to errors in weather and climate forecasts as well as account for a source of ...uncertainty in climate change projections. An ongoing challenge in rectifying these shortcomings has been the availability of adequate, high‐quality, global observations targeting ice clouds and related precipitating hydrometeors. In addition, the inadequacy of the modeled physics and the often disjointed nature between model representation and the characteristics of the retrieved/observed values have hampered GCM development and validation efforts from making effective use of the measurements that have been available. Thus, even though parameterizations in GCMs accounting for cloud ice processes have, in some cases, become more sophisticated in recent years, this development has largely occurred independently of the global‐scale measurements. With the relatively recent addition of satellite‐derived products from Aura/Microwave Limb Sounder (MLS) and CloudSat, there are now considerably more resources with new and unique capabilities to evaluate GCMs. In this article, we illustrate the shortcomings evident in model representations of cloud ice through a comparison of the simulations assessed in the Intergovernmental Panel on Climate Change Fourth Assessment Report, briefly discuss the range of global observational resources that are available, and describe the essential components of the model parameterizations that characterize their “cloud” ice and related fields. Using this information as background, we (1) discuss some of the main considerations and cautions that must be taken into account in making model‐data comparisons related to cloud ice, (2) illustrate present progress and uncertainties in applying satellite cloud ice (namely from MLS and CloudSat) to model diagnosis, (3) show some indications of model improvements, and finally (4) discuss a number of remaining questions and suggestions for pathways forward.
A cloud frequency of occurrence matrix is generated using merged cloud vertical profiles derived from the satellite‐borne Cloud‐Aerosol Lidar with Orthogonal Polarization (CALIOP) and cloud profiling ...radar. The matrix contains vertical profiles of cloud occurrence frequency as a function of the uppermost cloud top. It is shown that the cloud fraction and uppermost cloud top vertical profiles can be related by a cloud overlap matrix when the correlation length of cloud occurrence, which is interpreted as an effective cloud thickness, is introduced. The underlying assumption in establishing the above relation is that cloud overlap approaches random overlap with increasing distance separating cloud layers and that the probability of deviating from random overlap decreases exponentially with distance. One month of Cloud‐Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) and CloudSat data (July 2006) support these assumptions, although the correlation length sometimes increases with separation distance when the cloud top height is large. The data also show that the correlation length depends on cloud top hight and the maximum occurs when the cloud top height is 8 to 10 km. The cloud correlation length is equivalent to the decorrelation distance introduced by Hogan and Illingworth (2000) when cloud fractions of both layers in a two‐cloud layer system are the same. The simple relationships derived in this study can be used to estimate the top‐of‐atmosphere irradiance difference caused by cloud fraction, uppermost cloud top, and cloud thickness vertical profile differences.
Retrievals of cloud microphysical properties based on passive satellite imagery are especially difficult over snow‐covered surfaces because of the bright and cold surface. To help quantify their ...uncertainties, single‐layered overcast liquid‐phase Arctic stratus cloud microphysical properties retrieved by using the Clouds and the Earth's Radiant Energy System Edition 2 and Edition 4 (CERES Ed2 and Ed4) algorithms are compared with ground‐based retrievals at the Atmospheric Radiation Measurement North Slope of Alaska (ARM NSA) site at Barrow, AK, during the period from March 2000 to December 2006. A total of 206 and 140 snow‐free cases (Rsfc ≤ 0.3), and 108 and 106 snow cases (Rsfc > 0.3), respectively, were selected from Terra and Aqua satellite passes over the ARM NSA site. The CERES Ed4 and Ed2 optical depth (τ) and liquid water path (LWP) retrievals from both Terra and Aqua are almost identical and have excellent agreement with ARM retrievals under snow‐free and snow conditions. In order to reach a radiation closure study for both the surface and top of atmosphere (TOA) radiation budgets, the ARM precision spectral pyranometer‐measured surface albedos were adjusted (63.6% and 80% of the ARM surface albedos for snow‐free and snow cases, respectively) to account for the water and land components of the domain of 30 km × 30 km. Most of the radiative transfer model calculated SW↓sfc and SW↑TOA fluxes by using ARM and CERES cloud retrievals and the domain mean albedos as input agree with the ARM and CERES flux observations within 10 W m−2 for both snow‐free and snow conditions. Sensitivity studies show that the ARM LWP and re retrievals are less dependent on solar zenith angle (SZA), but all retrieved optical depths increase with SZA.
Key Points
The CERES‐MODIS retrieved cloud microphysical properties agree well with ARM retrievals under both snow‐free and snow conditions
A radiation closure has been reached at both surface and TOA for both snow‐free and snow conditions
A domain mean albedo is used to make apples‐to‐apples comparisons in cloud properties and to reach a radiation closure study
The impact of horizontal heterogeneities, liquid water path (LWP from AMSR-E), and cloud fraction (CF) on MODIS cloud effective radius (re), retrieved from the 2.1 μm (re2.1) and 3.8 μm (re3.8) ...channels, is investigated for warm clouds over the southeast Pacific. Values of re retrieved using the CERES algorithms are averaged at the CERES footprint resolution (∼20 km), while heterogeneities (Hσ) are calculated as the ratio between the standard deviation and mean 0.64 μm reflectance. The value of re2.1 strongly depends on CF, with magnitudes up to 5 μm larger than those for overcast scenes, whereas re3.8 remains insensitive to CF. For cloudy scenes, both re2.1 and re3.8 increase with Hσ for any given AMSR-E LWP, but re2.1 changes more than for re3.8. Additionally, re3.8–re2.1 differences are positive (<1 μm) for homogeneous scenes (Hσ < 0.2) and LWP > 45 gm−2, and negative (up to −4 μm) for larger Hσ. While re3.8–re2.1 differences in homogeneous scenes are qualitatively consistent with in situ microphysical observations over the region of study, negative differences – particularly evinced in mean regional maps – are more likely to reflect the dominant bias associated with cloud heterogeneities rather than information about the cloud vertical structure. The consequences for MODIS LWP are also discussed.
The National Aeronautics and Space Administration's Clouds and the Earth's Radiant Energy System (CERES) Project was designed to improve our understanding of the relationship between clouds and solar ...and longwave radiation. This is achieved using satellite broad-band instruments to map the top-of-atmosphere radiation fields with coincident data from satellite narrow-band imagers employed to retrieve the properties of clouds associated with those fields. This paper documents the CERES Edition-2 cloud property retrieval system used to analyze data from the Tropical Rainfall Measuring Mission Visible and Infrared Scanner and by the MODerate-resolution Imaging Spectrometer instruments on board the Terra and Aqua satellites covering the period 1998 through 2007. Two daytime retrieval methods are explained: the Visible Infrared Shortwave-infrared Split-window Technique for snow-free surfaces and the Shortwave-infrared Infrared Near-infrared Technique for snow or ice-covered surfaces. The Shortwave-infrared Infrared Split-window Technique is used for all surfaces at night. These methods, along with the ancillary data and empirical parameterizations of cloud thickness, are used to derive cloud boundaries, phase, optical depth, effective particle size, and condensed/frozen water path at both pixel and CERES footprint levels. Additional information is presented, detailing the potential effects of satellite calibration differences, highlighting methods to compensate for spectral differences and correct for atmospheric absorption and emissivity, and discussing known errors in the code. Because a consistent set of algorithms, auxiliary input, and calibrations across platforms are used, instrument and algorithm-induced changes in the data record are minimized. This facilitates the use of the CERES data products for studying climate-scale trends.
One year of instantaneous top‐of‐atmosphere (TOA) and surface shortwave and longwave irradiances are computed using cloud and aerosol properties derived from instruments on the A‐Train Constellation: ...the Cloud‐Aerosol Lidar with Orthogonal Polarization (CALIOP) on the Cloud‐Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite, the CloudSat Cloud Profiling Radar (CPR), and the Aqua Moderate Resolution Imaging Spectrometer (MODIS). When modeled irradiances are compared with those computed with cloud properties derived from MODIS radiances by a Clouds and the Earth's Radiant Energy System (CERES) cloud algorithm, the global and annual mean of modeled instantaneous TOA irradiances decreases by 12.5 W m−2 (5.0%) for reflected shortwave and 2.5 W m−2 (1.1%) for longwave irradiances. As a result, the global annual mean of instantaneous TOA irradiances agrees better with CERES‐derived irradiances to within 0.5W m−2 (out of 237.8 W m−2) for reflected shortwave and 2.6W m−2 (out of 240.1 W m−2) for longwave irradiances. In addition, the global annual mean of instantaneous surface downward longwave irradiances increases by 3.6 W m−2 (1.0%) when CALIOP‐ and CPR‐derived cloud properties are used. The global annual mean of instantaneous surface downward shortwave irradiances also increases by 8.6 W m−2 (1.6%), indicating that the net surface irradiance increases when CALIOP‐ and CPR‐derived cloud properties are used. Increasing the surface downward longwave irradiance is caused by larger cloud fractions (the global annual mean by 0.11, 0.04 excluding clouds with optical thickness less than 0.3) and lower cloud base heights (the global annual mean by 1.6 km). The increase of the surface downward longwave irradiance in the Arctic exceeds 10 W m−2 (∼4%) in winter because CALIOP and CPR detect more clouds in comparison with the cloud detection by the CERES cloud algorithm during polar night. The global annual mean surface downward longwave irradiance of 345.4 W m−2 is estimated by combining the modeled instantaneous surface longwave irradiance computed with CALIOP and CPR cloud profiles with the global annual mean longwave irradiance from the CERES product (AVG), which includes the diurnal variation of the irradiance. The estimated bias error is −1.5 W m−2 and the uncertainty is 6.9 W m−2. The uncertainty is predominately caused by the near‐surface temperature and column water vapor amount uncertainties.
Key Points
CALIPSO CloudSat increases surface downward LW to 345Wm‐2
Multi‐layer clouds provide better agreement with CERES
Surface radiation and other surface flux uncertainties overlap
Marine boundary layer (MBL) cloud properties derived from the NASA Clouds and the Earth's Radiant Energy System (CERES) project using Terra and Aqua Moderate Resolution Imaging Spectroradiometer ...(MODIS) data are compared with observations taken at the Department of Energy Atmospheric Radiation Measurement (ARM) Mobile Facility at the Azores (AMF‐Azores) site from June 2009 through December 2010. Cloud properties derived from ARM ground‐based observations were averaged over a 1 h interval centered at the satellite overpass time, while the CERES‐MODIS (CM) results were averaged within a 30 km × 30 km grid box centered over the Azores site. A total of 63 daytime and 92 nighttime single‐layered overcast MBL cloud cases were selected from 19 months of ARM radar‐lidar and satellite observations. The CM cloud top/base heights (Htop/Hbase) were determined from cloud top/base temperatures (Ttop/Tbase) using a regional boundary layer lapse rate method. For daytime comparisons, the CM‐derived Htop (Hbase), on average, is 0.063 km (0.068 km) higher (lower) than its ARM radar‐lidar‐observed counterpart, and the CM‐derived Ttop and Tbase are 0.9 K less and 2.5 K greater than the surface values with high correlations (R2 = 0.82 and 0.84, respectively). In general, the cloud top comparisons agree better than the cloud base comparisons, because the CM cloud base temperatures and heights are secondary products determined from cloud top temperatures and heights. No significant day‐night difference was found in the analyses. The comparisons of MBL cloud microphysical properties reveal that when averaged over a 30 km × 30 km area, the CM‐retrieved cloud droplet effective radius (re) at 3.7 µm is 1.3 µm larger than that from the ARM retrievals (12.8 µm), while the CM‐retrieved cloud liquid water path (LWP) is 13.5 gm−2 less than its ARM counterpart (114.2 gm−2) due to its small optical depth (9.6 versus 13.7). The differences are reduced by 50% when the CM averages are computed only using the MODIS pixel nearest the AMF site. Using the effective radius retrieved using 2.1 µm channel to calculate LWP can reduce the difference between the CM and ARM microwave radiometer retrievals from −13.7 to 2.1 gm−2. The 10% differences between the ARM and CERES‐MODIS LWP and re retrievals are within the uncertainties of the ARM LWP (~20 gm−2) and re (~10%) retrievals; however, the 30% difference in optical depth is significant. Possible reasons contributing to this discrepancy are increased sensitivities in optical depth from both surface retrievals when τ~10 and topography. The τ differences vary with wind direction and are consistent with the island orography. Much better agreement in τ is obtained when using only those data taken when the wind is from the northeast, where topographical effects on the sampled clouds are minimal.
Key Points
CERES‐MODIS‐retrieved cloud height/temperature agrees well to the ARM measurementUsing re retrieved at 2.1 µm channel to calculate LWP can reduce the differenceThe discrepancy in tau from both surface retrievals and possible island effects