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.
Analysis of one decade of radar‐lidar and Geostationary Operational Environmental Satellite (GOES) observations at the Department of Energy (DOE) Atmospheric Radiation Measurement Program (ARM) ...Southern Great Plains (SGP) site reveals that there is excellent agreement in the long‐term mean cloud fractions (CFs) derived from the surface and GOES data, and the CF is independent of temporal resolution and spatial scales for grid boxes of size 0.5° to 2.5°. When computed over a a 0.5 h (4 h) period, cloud frequency of occurrence (FREQ) and amount when present (AWP) derived from the point surface data agree very well with the same quantities determined from GOES for a 0.5° (2.5°) region centered on the DOE ARM SGP site. The values of FREQ (AWP) derived from the radar‐lidar observations at a given altitude increase (decrease) as the averaging period increases from 5 min to 6 h. Similarly, CF at a given altitude increases as the vertical resolution increases from 90 to 1000 m. The profiles of CF have distinct bimodal vertical distributions, with a lower peak between 1 and 2 km and a higher one between 8 and 11 km. The 10 year mean total CF, 46.9%, varies seasonally from a summer minimum of 39.8% to a maximum of 54.6% during the winter. The annual mean CF is 1%–2% less than that from previous studies, ∼48%–49%, because fewer clouds occurred during 2005 and 2006, especially during winter. The differences in single‐ and multilayered CFs between this study and an earlier analysis can be explained by the different temporal resolutions used in the two studies, where single‐layered CFs decrease but multilayered CFs increase from a 5 min resolution to a 1 h resolution. The vertical distribution of nighttime GOES high cloud tops agrees well with surface observations, but during the daytime, fewer high clouds are retrieved by the GOES analysis than seen from the surface observations. The FREQs for both daytime and nighttime GOES low cloud tops are significantly higher than surface observations, but the CFs are in good agreement.
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
The mean structure and diurnal cycle of southeast (SE) Atlantic boundary layer clouds are described with satellite observations and multiscale modeling framework (MMF) simulations during austral ...spring (September–November). Hourly resolution cloud fraction (CF) and cloud-top height (HT
) are retrieved fromMeteosat-9radiances using modified Clouds and the Earth’s Radiant Energy System (CERES) Moderate Resolution Imaging Spectroradiometer (MODIS) algorithms, whereas liquid water path (LWP) is from the University of Wisconsin microwave satellite climatology. The MMF simulations use a 2D cloud-resolving model (CRM) that contains an advanced third-order turbulence closure to explicitly simulate cloud physical processes in every grid column of a general circulation model. The model accurately reproduces themarine stratocumulus spatial extent and cloud cover. The mean cloud cover spatial variability in the model is primarily explained by the boundary layer decoupling strength, whereas a boundary layer shoaling accounts for a coastal decrease in CF. Moreover, the core of the stratocumulus cloud deck is concomitant with the location of the strongest temperature inversion. Although the model reproduces the observed westward boundary layer deepening and the spatial variability of LWP, it overestimates LWP by 50%. Diurnal cycles ofHT
, CF, and LWP from satellites and the model have the same phase, with maxima during the early morning and minima near 1500 local solar time, which suggests that the diurnal cycle is driven primarily by solar heating. Comparisons with the SE Pacific cloud deck indicate that the observed amplitude of the diurnal cycle is modest over the SE Atlantic, with a shallower boundary layer as well. The model qualitatively reproduces these interregime differences.
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.
Analysis of the cloud components in numerical weather prediction models using advanced data assimilation techniques has been a prime topic in recent years. In this research, the variational data ...assimilation (DA) system for the Weather Research and Forecasting (WRF) Model (WRFDA) is further developed to assimilate satellite cloud products that will produce the cloud liquid water and ice water analysis. Observation operators for the cloud liquid water path and cloud ice water path are developed and incorporated into the WRFDA system. The updated system is tested by assimilating cloud liquid water path and cloud ice water path observations from Global Geostationary Gridded Cloud Products at NASA. To assess the impact of cloud liquid/ice water path data assimilation on short-term regional numerical weather prediction (NWP), 3-hourly cycling data assimilation and forecast experiments with and without the use of the cloud liquid/ice water paths are conducted. It is shown that assimilating cloud liquid/ice water paths increases the accuracy of temperature, humidity, and wind analyses at model levels between 300 and 150 hPa after 5 cycles (15 h). It is also shown that assimilating cloud liquid/ice water paths significantly reduces forecast errors in temperature and wind at model levels between 300 and 150 hPa. The precipitation forecast skills are improved as well. One reason that leads to the improved analysis and forecast is that the 3-hourly rapid update cycle carries over the impact of cloud information from the previous cycles spun up by the WRF Model.
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.
Cloud properties are essential for the Clouds and the Earth’s Radiant Energy System (CERES) Project, enabling accurate interpretation of measured broadband radiances, providing a means to understand ...global cloud-radiation interactions, and constituting an important climate record. Producing consistent cloud retrievals across multiple platforms is critical for generating a multidecadal cloud and radiation record. Techniques used by CERES for retrievals from measurements by the MODerate-Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua platforms are adapted for the application to radiances from the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership to continue the CERES record beyond the MODIS era. The algorithm adjustments account for spectral and channel differences, use revised reflectance models, and set new thresholds for detecting thin cirrus clouds at night. Cloud amounts from VIIRS are less than their MODIS counterparts by 0.016 during the day and 0.026 at night, but trend consistently over the 2012–2020 period. The VIIRS mean liquid water cloud fraction differs by ~0.01 from the MODIS amount. The average cloud heights from VIIRS differ from the MODIS heights by less than 0.2 km, except the VIIRS daytime ice cloud heights, which are 0.4 km higher. The mean VIIRS nonpolar optical depths are 17% (1%) larger (smaller) than those from MODIS for liquid (ice) clouds. The VIIRS cloud hydrometeor sizes are generally smaller than their MODIS counterparts. Discrepancies between the MODIS and VIIRS properties stem from spectral and spatial resolution differences, new tests at night, calibration inconsistencies, and new reflectance models. Many of those differences will be addressed in future editions.
The decades-long Clouds and Earth’s Radiant Energy System (CERES) Project includes both cloud and radiation measurements from instruments on the Aqua, Terra, and Suomi National Polar-orbiting ...Partnership (SNPP) satellites. To build a reliable long-term climate data record, it is important to determine the accuracies of the parameters retrieved from the sensors on each satellite. Cloud amount, phase, and top height derived from radiances taken by the Visible Infrared Imaging Radiometer Suite (VIIRS) on the SNPP are evaluated relative to the same quantities determined from measurements by the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on the Cloud Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) spacecraft. The accuracies of the VIIRS cloud fractions are found to be as good as or better than those for the CERES amounts determined from Aqua MODerate-resolution Imaging Spectroradiometer (MODIS) data and for cloud fractions estimated by two other operational algorithms. Sensitivities of cloud fraction bias to CALIOP resolution, matching time window, and viewing zenith angle are examined. VIIRS cloud phase biases are slightly greater than their CERES MODIS counterparts. A majority of cloud phase errors are due to multilayer clouds during the daytime and supercooled liquid water clouds at night. CERES VIIRS cloud-top height biases are similar to those from CERES MODIS, except for ice clouds, which are smaller than those from CERES MODIS. CERES VIIRS cloud phase and top height uncertainties overall are very similar to or better than several operational algorithms, but fail to match the accuracies of experimental machine learning techniques. The greatest errors occur for multilayered clouds and clouds with phase misclassification. Cloud top heights can be improved by relaxing tropopause constraints, improving lapse-rate to model temperature profiles, and accounting for multilayer clouds. Other suggestions for improving the retrievals are also discussed.
As part of the Southeast United States-based Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC4RS), and collinear with part of the Southeast ...Atmosphere Study, the University of Wisconsin High Spectral Resolution Lidar system was deployed to the University of Alabama from 19 June to 4 November 2013. With a collocated Aerosol Robotic Network (AERONET) sun photometer, a nearby Chemical Speciation Network (PM2.5) measurement station, and near daily ozonesonde releases for the August-September SEAC4RS campaign, the site allowed the regions first comprehensive diurnal monitoring of aerosol particle vertical structure. A 532nm lidar ratio of 55 sr provided good closure between aerosol backscatter and AERONET (aerosol optical thickness, AOT). A principle component analysis was performed to identify key modes of variability in aerosol backscatter. ''Fair weather'' days exhibited classic planetary boundary layer structure of a mixed layer accounting for approx. 50% of AOT and an entrainment zone providing another 25%. An additional 5-15% of variance is gained from the lower free troposphere from either convective detrainment or frequent intrusions of western United States biomass burning smoke. Generally, aerosol particles were contained below the 0 C level, a common level of stability in convective regimes. However, occasional strong injections of smoke to the upper troposphere were also observed, accounting for the remaining 10-15% variability in AOT. Examples of these common modes of variability in frontal and convective regimes are presented, demonstrating why AOT often has only a weak relationship to surface PM2.5 concentration.