Cloud droplet size and liquid water path retrievals from zenith radiance measurements: examples from the Atmospheric Radiation Measurement Program and the Aerosol Robotic Network

The ground-based Atmospheric Radiation Measurement Program (ARM) and NASA Aerosol Robotic Net- work (AERONET) routinely monitor clouds using zenith ra- diances at visible and near-infrared wavelengths. Using the transmittance calculated from such measurements, we have developed a new retrieval method for cloud effective droplet size and conducted extensive tests for non-precipitating liquid water clouds. The underlying principle is to combine a liquid-water-absorbing wavelength (i.e., 1640 nm) with a non-water-absorbing wavelength for acquiring information on cloud droplet size and optical depth. For simulated stratocumulus clouds with liquid water path less than 300 g m−2 and horizontal resolution of 201 m, the retrieval method underestimates the mean effective radius by 0.8μm, with a root-mean-squared error of 1.7 μm and a relative deviation of 13%. For actual observations with a liquid water path less than 450 g m−2 at the ARM Oklahoma site during 2007– 2008, our 1.5-min-averaged retrievals are generally larger by around 1 μm than those from combined ground-based cloud radar and microwave radiometer at a 5-min temporal resolution. We also compared our retrievals to those from combined shortwave flux and microwave observations for relatively homogeneous clouds, showing that the bias between these two retrieval sets is negligible, but the error of 2.6 μm and the relative deviation of 22 % are larger than those found in our simulation case. Finally, the transmittance-based cloud effective droplet radii agree to better than 11 % with satellite observations and have a negative bias of 1 μm. Overall, the retrieval method provides reasonable cloud effective radius estimates, which can enhance the cloud products of both ARM and AERONET.

Diurnal cycle of land surface temperature in a desert encroachment zone as observed from satellites

Climate variability in the African Soudano-Sahel savanna zone has attracted much attention because of the persistence of anomalously low rainfall. Past efforts to monitor the climate of this region have focused on rainfall and vegetation conditions, while land surface temperature (LST) has received less attention. Remote sensing of LST is feasible and possible at global scale. Most remotely sensed estimates of LST are based on the National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR) that are limited in their ability to capture the full diurnal cycle. Although more frequent observations are available from past geostationary satellites, their spatial resolution is coarser than that of polar orbiting satellites. In this study, the improved capabilities of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on the METEOSAT Second Generation (MSG) instrument are used to remotely sense the LST in the African Soudano-Sahel savanna zone at a resolution of 3 km and 15 minutes. In support of the Radiative Atmospheric Divergence using the ARM Mobile Facility (AMF), GERB and AMMA Stations (RADAGAST) project, African Monsoon Multidisciplinary Analyses (AMMA) project and the Department of Energy's Atmospheric Radiation Measurement (ARM) program, the ARM Mobile Facility was deployed during 2006 in this climatically sensitive region, thereby providing a unique opportunity to evaluate remotely sensed algorithms for deriving LST.

Overview of observations from the RADAGAST experiment in Niamey, Niger: 2. Radiative fluxes and divergences

Broadband shortwave and longwave radiative fluxes observed both at the surface and from space during the Radiative Atmospheric Divergence using ARM Mobile Facility, GERB data and AMMA Stations (RADAGAST) experiment in Niamey, Niger, in 2006 are presented. The surface fluxes were measured by the Atmospheric Radiation Measurement (ARM) Program Mobile Facility (AMF) at Niamey airport, while the fluxes at the top of the atmosphere (TOA) are from the Geostationary Earth Radiation Budget (GERB) instrument on the Meteosat-8 satellite. The data are analyzed as daily averages, in order to minimize sampling differences between the surface and top of atmosphere instruments, while retaining the synoptic and seasonal changes that are the main focus of this study. A cloud mask is used to identify days with cloud versus those with predominantly clear skies. The influence of temperature, water vapor, aerosols, and clouds is investigated. Aerosols are ubiquitous throughout the year and have a significant impact on both the shortwave and longwave fluxes. The large and systematic seasonal changes in temperature and column integrated water vapor (CWV) through the dry and wet seasons are found to exert strong influences on the longwave fluxes. These influences are often in opposition to each other, because the highest temperatures occur at the end of the dry season when the CWV is lowest, while in the wet season the lowest temperatures are associated with the highest values of CWV. Apart from aerosols, the shortwave fluxes are also affected by clouds and by the seasonal changes in CWV. The fluxes are combined to provide estimates of the divergence of radiation across the atmosphere throughout 2006. The longwave divergence shows a relatively small variation through the year, because of a partial compensation between the seasonal variations in the outgoing longwave radiation (OLR) and surface net longwave radiation. A simple model of the greenhouse effect is used to interpret this result in terms of the dependence of the normalized greenhouse effect at the TOA and of the effective emissivity of the atmosphere at the surface on the CWV. It is shown that, as the CWV increases, the atmosphere loses longwave energy to the surface with about the same increasing efficiency with which it traps the OLR. When combined with the changes in temperature, this maintains the atmospheric longwave divergence within the narrow range that is observed. The shortwave divergence is mainly determined by the CWV and aerosol loadings and the effect of clouds is much smaller than on the component fluxes.

Cloud optical depth retrievals from the Aerosol Robotic Network (AERONET) cloud mode observations

Cloud optical depth is one of the most poorly observed climate variables. The new “cloud mode” capability in the Aerosol Robotic Network (AERONET) will inexpensively yet dramatically increase cloud optical depth observations in both number and accuracy. Cloud mode optical depth retrievals from AERONET were evaluated at the Atmospheric Radiation Measurement program’s Oklahoma site in sky conditions ranging from broken clouds to overcast. For overcast cases, the 1.5 min average AERONET cloud mode optical depths agreed to within 15% of those from a standard ground‐based flux method. For broken cloud cases, AERONET retrievals also captured rapid variations detected by the microwave radiometer. For 3 year climatology derived from all nonprecipitating clouds, AERONET monthly mean cloud optical depths are generally larger than cloud radar retrievals because of the current cloud mode observation strategy that is biased toward measurements of optically thick clouds. This study has demonstrated a new way to enhance the existing AERONET infrastructure to observe cloud optical properties on a global scale.

Spectrally-invariant behavior of zenith radiance around cloud edges simulated by radiative transfer

In a previous paper, we discovered a surprising spectrally-invariant relationship in shortwave spectrometer observations taken by the Atmospheric Radiation Measurement (ARM) program. The relationship suggests that the shortwave spectrum near cloud edges can be determined by a linear combination of zenith radiance spectra of the cloudy and clear regions. Here, using radiative transfer simulations, we study the sensitivity of this relationship to the properties of aerosols and clouds, to the underlying surface type, and to the finite field-of-view (FOV) of the spectrometer. Overall, the relationship is mostly sensitive to cloud properties and has little sensitivity to other factors. At visible wavelengths, the relationship primarily depends on cloud optical depth regardless of cloud phase function, thermodynamic phase and drop size. At water-absorbing wavelengths, the slope of the relationship depends primarily on cloud optical depth; the intercept, by contrast, depends primarily on cloud absorbing and scattering properties, suggesting a new retrieval method for cloud drop effective radius. These results suggest that the spectrally-invariant relationship can be used to infer cloud properties near cloud edges even with insufficient or no knowledge about spectral surface albedo and aerosol properties.

An assessment of aerosol‐cloud interactions in marine stratus clouds based on surface remote sensing

An assessment of aerosol-cloud interactions (ACI) from ground-based remote sensing under coastal stratiform clouds is presented. The assessment utilizes a long-term, high temporal resolution data set from the Atmospheric Radiation Measurement (ARM) Program deployment at Pt. Reyes, California, United States, in 2005 to provide statistically robust measures of ACI and to characterize the variability of the measures based on variability in environmental conditions and observational approaches. The average ACIN (= dlnNd/dlna, the change in cloud drop number concentration with aerosol concentration) is 0.48, within a physically plausible range of 0–1.0. Values vary between 0.18 and 0.69 with dependence on (1) the assumption of constant cloud liquid water path (LWP), (2) the relative value of cloud LWP, (3) methods for retrieving Nd, (4) aerosol size distribution, (5) updraft velocity, and (6) the scale and resolution of observations. The sensitivity of the local, diurnally averaged radiative forcing to this variability in ACIN values, assuming an aerosol perturbation of 500 c-3 relative to a background concentration of 100 cm-3, ranges betwee-4 and -9 W -2. Further characterization of ACI and its variability is required to reduce uncertainties in global radiative forcing estimates.

Remote sensing of cloud properties using ground-based measurements of zenith radiance

We have conducted the first extensive field test of two new methods to retrieve optical properties for overhead clouds that range from patchy to overcast. The methods use measurements of zenith radiance at 673 and 870 nm wavelengths and require the presence of green vegetation in the surrounding area. The test was conducted at the Atmospheric Radiation Measurement Program Oklahoma site during September–November 2004. These methods work because at 673 nm (red) and 870 nm (near infrared (NIR)), clouds have nearly identical optical properties, while vegetated surfaces reflect quite differently. The first method, dubbed REDvsNIR, retrieves not only cloud optical depth τ but also radiative cloud fraction. Because of the 1-s time resolution of our radiance measurements, we are able for the first time to capture changes in cloud optical properties at the natural timescale of cloud evolution. We compared values of τ retrieved by REDvsNIR to those retrieved from downward shortwave fluxes and from microwave brightness temperatures. The flux method generally underestimates τ relative to the REDvsNIR method. Even for overcast but inhomogeneous clouds, differences between REDvsNIR and the flux method can be as large as 50%. In addition, REDvsNIR agreed to better than 15% with the microwave method for both overcast and broken clouds. The second method, dubbed COUPLED, retrieves τ by combining zenith radiances with fluxes. While extra information from fluxes was expected to improve retrievals, this is not always the case. In general, however, the COUPLED and REDvsNIR methods retrieve τ to within 15% of each other.

Moisture vertical structure, column water vapor, and tropical deep convection

The vertical structure of the relationship between water vapor and precipitation is analyzed in 5 yr of radiosonde and precipitation gauge data from the Nauru Atmospheric Radiation Measurement (ARM) site. The first vertical principal component of specific humidity is very highly correlated with column water vapor (CWV) and has a maximum of both total and fractional variance captured in the lower free troposphere (around 800 hPa). Moisture profiles conditionally averaged on precipitation show a strong association between rainfall and moisture variability in the free troposphere and little boundary layer variability. A sharp pickup in precipitation occurs near a critical value of CWV, confirming satellite-based studies. A lag–lead analysis suggests it is unlikely that the increase in water vapor is just a result of the falling precipitation. To investigate mechanisms for the CWV–precipitation relationship, entraining plume buoyancy is examined in sonde data and simplified cases. For several different mixing schemes, higher CWV results in progressively greater plume buoyancies, particularly in the upper troposphere, indicating conditions favorable for deep convection. All other things being equal, higher values of lower-tropospheric humidity, via entrainment, play a major role in this buoyancy increase. A small but significant increase in subcloud layer moisture with increasing CWV also contributes to buoyancy. Entrainment coefficients inversely proportional to distance from the surface, associated with mass flux increase through a deep lower-tropospheric layer, appear promising. These yield a relatively even weighting through the lower troposphere for the contribution of environmental water vapor to midtropospheric buoyancy, explaining the association of CWV and buoyancy available for deep convection.

Cloud statistics and cloud radiative effect for a low-mountain site

In 2007, the Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) was operated for a nine-month period in the Murg Valley, Black Forest, Germany, in support of the Convective and Orographically-induced Precipitation Study (COPS). The synergy of AMF and COPS partner instrumentation was exploited to derive a set of high-quality thermodynamic and cloud property profiles with 30 s resolution. In total, clouds were present 72% of the time, with multi-layer mixed phase (28.4%) and single-layer water clouds (11.3%) occurring most frequently. A comparison with the Cloudnet sites Chilbolton and Lindenberg for the same time period revealed that the Murg Valley exhibits lower liquid water paths (LWPs; median = 37.5 g m−2) compared to the two sites located in flat terrain. In order to evaluate the derived thermodynamic and cloud property profiles, a radiative closure study was performed with independent surface radiation measurements. In clear sky, average differences between calculated and observed surface fluxes are less than 2% and 4% for the short wave and long wave part, respectively. In cloudy situations, differences between simulated and observed fluxes, particularly in the short wave part, are much larger, but most of these can be related to broken cloud situations. The daytime cloud radiative effect (CRE), i.e. the difference of cloudy and clear-sky net fluxes, has been analysed for the whole nine-month period. For overcast, single-layer water clouds, sensitivity studies revealed that the CRE uncertainty is likewise determined by uncertainties in liquid water content and effective radius. For low LWP clouds, CRE uncertainty is dominated by LWP uncertainty; therefore refined retrievals, such as using infrared and/or higher microwave frequencies, are needed.

The evaluation of CloudSat and CALIPSO ice microphysical products using ground-based cloud radar and lidar observations

In this paper, the statistical properties of tropical ice clouds (ice water content, visible extinction, effective radius, and total number concentration) derived from 3 yr of ground-based radar–lidar retrievals from the U.S. Department of Energy Atmospheric Radiation Measurement Climate Research Facility in Darwin, Australia, are compared with the same properties derived using the official CloudSat microphysical retrieval methods and from a simpler statistical method using radar reflectivity and air temperature. It is shown that the two official CloudSat microphysical products (2B-CWC-RO and 2B-CWC-RVOD) are statistically virtually identical. The comparison with the ground-based radar–lidar retrievals shows that all satellite methods produce ice water contents and extinctions in a much narrower range than the ground-based method and overestimate the mean vertical profiles of microphysical parameters below 10-km height by over a factor of 2. Better agreements are obtained above 10-km height. Ways to improve these estimates are suggested in this study. Effective radii retrievals from the standard CloudSat algorithms are characterized by a large positive bias of 8–12 μm. A sensitivity test shows that in response to such a bias the cloud longwave forcing is increased from 44.6 to 46.9 W m−2 (implying an error of about 5%), whereas the negative cloud shortwave forcing is increased from −81.6 to −82.8 W m−2. Further analysis reveals that these modest effects (although not insignificant) can be much larger for optically thick clouds. The statistical method using CloudSat reflectivities and air temperature was found to produce inaccurate mean vertical profiles and probability distribution functions of effective radius. This study also shows that the retrieval of the total number concentration needs to be improved in the official CloudSat microphysical methods prior to a quantitative use for the characterization of tropical ice clouds. Finally, the statistical relationship used to produce ice water content from extinction and air temperature obtained by the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite is evaluated for tropical ice clouds. It is suggested that the CALIPSO ice water content retrieval is robust for tropical ice clouds, but that the temperature dependence of the statistical relationship used should be slightly refined to better reproduce the radar–lidar retrievals.

The statistical properties of tropical ice clouds generated by the West-African and Australian monsoons from ground-based radar-lidar observations

In the present paper we characterize the statistical properties of non-precipitating tropical ice clouds (deep ice anvils resulting from deep convection and cirrus clouds) over Niamey, Niger, West Africa, and Darwin, northern Australia, using ground-based radar–lidar observations from the Atmospheric Radiation Measurement (ARM) programme. The ice cloud properties analysed in this paper are the frequency of ice cloud occurrence, cloud fraction, the morphological properties (cloud-top height, base height, and thickness), the microphysical and radiative properties (ice water content, visible extinction, effective radius, terminal fall speed, and concentration), and the internal cloud dynamics (in-cloud vertical air velocity). The main highlight of the paper is that it characterizes for the first time the probability density functions of the tropical ice cloud properties, their vertical variability and their diurnal variability at the same time. This is particularly important over West Africa, since the ARM deployment in Niamey provides the first vertically resolved observations of non-precipitating ice clouds in this crucial area in terms of redistribution of water and energy in the troposphere. The comparison between the two sites also provides an additional observational basis for the evaluation of the parametrization of clouds in large-scale models, which should be able to reproduce both the statistical properties at each site and the differences between the two sites. The frequency of ice cloud occurrence is found to be much larger over Darwin when compared to Niamey, and with a much larger diurnal variability, which is well correlated with the diurnal cycle of deep convective activity. The diurnal cycle of the ice cloud occurrence over Niamey is also much less correlated with that of deep convective activity than over Darwin, probably owing to the fact that Niamey is further away from the deep convective sources of the region. The frequency distributions of cloud fraction are strongly bimodal and broadly similar over the two sites, with a predominance of clouds characterized either by a very small cloud fraction (less than 0.3) or a very large cloud fraction (larger than 0.9). The ice clouds over Darwin are also much thicker (by 1 km or more statistically) and are characterized by a much larger diurnal variability than ice clouds over Niamey. Ice clouds over Niamey are also characterized by smaller particle sizes and fall speeds but in much larger concentrations, thereby carrying more ice water and producing more visible extinction than the ice clouds over Darwin. It is also found that there is a much larger occurrence of downward in-cloud air motions less than 1 m s−1 over Darwin, which together with the larger fall speeds retrieved over Darwin indicates that the life cycle of ice clouds is probably shorter over Darwin than over Niamey.

CloudSat as a Global Radar Calibrator

The calibration of the CloudSat spaceborne cloud radar has been thoroughly assessed using very accurate internal link budgets before launch, comparisons with predicted ocean surface backscatter at 94 GHz, direct comparisons with airborne cloud radars, and statistical comparisons with ground-based cloud radars at different locations of the world. It is believed that the calibration of CloudSat is accurate to within 0.5–1 dB. In the present paper it is shown that an approach similar to that used for the statistical comparisons with ground-based radars can now be adopted the other way around to calibrate other ground-based or airborne radars against CloudSat and/or to detect anomalies in long time series of ground-based radar measurements, provided that the calibration of CloudSat is followed up closely (which is the case). The power of using CloudSat as a global radar calibrator is demonstrated using the Atmospheric Radiation Measurement cloud radar data taken at Barrow, Alaska, the cloud radar data from the Cabauw site, Netherlands, and airborne Doppler cloud radar measurements taken along the CloudSat track in the Arctic by the Radar System Airborne (RASTA) cloud radar installed in the French ATR-42 aircraft for the first time. It is found that the Barrow radar data in 2008 are calibrated too high by 9.8 dB, while the Cabauw radar data in 2008 are calibrated too low by 8.0 dB. The calibration of the RASTA airborne cloud radar using direct comparisons with CloudSat agrees well with the expected gains and losses resulting from the change in configuration that required verification of the RASTA calibration.

The variability of tropical ice cloud properties as a function of the large-scale context from ground-based radar-lidar observations over Darwin, Australia

The high complexity of cloud parameterizations now held in models puts more pressure on observational studies to provide useful means to evaluate them. One approach to the problem put forth in the modelling community is to evaluate under what atmospheric conditions the parameterizations fail to simulate the cloud properties and under what conditions they do a good job. It is the ambition of this paper to characterize the variability of the statistical properties of tropical ice clouds in different tropical "regimes" recently identified in the literature to aid the development of better process-oriented parameterizations in models. For this purpose, the statistical properties of non-precipitating tropical ice clouds over Darwin, Australia are characterized using ground-based radar-lidar observations from the Atmospheric Radiation Measurement (ARM) Program. The ice cloud properties analysed are the frequency of ice cloud occurrence, the morphological properties (cloud top height and thickness), and the microphysical and radiative properties (ice water content, visible extinction, effective radius, and total concentration). The variability of these tropical ice cloud properties is then studied as a function of the large-scale cloud regimes derived from the International Satellite Cloud Climatology Project (ISCCP), the amplitude and phase of the Madden-Julian Oscillation (MJO), and the large-scale atmospheric regime as derived from a long-term record of radiosonde observations over Darwin. The vertical variability of ice cloud occurrence and microphysical properties is largest in all regimes (1.5 order of magnitude for ice water content and extinction, a factor 3 in effective radius, and three orders of magnitude in concentration, typically). 98 % of ice clouds in our dataset are characterized by either a small cloud fraction (smaller than 0.3) or a very large cloud fraction (larger than 0.9). In the ice part of the troposphere three distinct layers characterized by different statistically-dominant microphysical processes are identified. The variability of the ice cloud properties as a function of the large-scale atmospheric regime, cloud regime, and MJO phase is large, producing mean differences of up to a factor 8 in the frequency of ice cloud occurrence between large-scale atmospheric regimes and mean differences of a factor 2 typically in all microphysical properties. Finally, the diurnal cycle of the frequency of occurrence of ice clouds is also very different between regimes and MJO phases, with diurnal amplitudes of the vertically-integrated frequency of ice cloud occurrence ranging from as low as 0.2 (weak diurnal amplitude) to values in excess of 2.0 (very large diurnal amplitude). Modellers should now use these results to check if their model cloud parameterizations are capable of translating a given atmospheric forcing into the correct statistical ice cloud properties.

RACORO extended-term, aircraft observations of boundary-layer clouds

A first-of-a-kind, extended-term cloud aircraft campaign was conducted to obtain an in-situ statistical characterization of continental boundary-layer clouds needed to investigate cloud processes and refine retrieval algorithms. Coordinated by the Atmospheric Radiation Measurement (ARM) Aerial Facility (AAF), the Routine AAF Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations (RACORO) field campaign operated over the ARM Southern Great Plains (SGP) site from 22 January to 30 June 2009, collecting 260 h of data during 59 research flights. A comprehensive payload aboard the Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft measured cloud microphysics, solar and thermal radiation, physical aerosol properties, and atmospheric state parameters. Proximity to the SGP's extensive complement of surface measurements provides ancillary data that supports modeling studies and facilitates evaluation of a variety of surface retrieval algorithms. The five-month duration enabled sampling a range of conditions associated with the seasonal transition from winter to summer. Although about two-thirds of the cloud flights occurred in May and June, boundary-layer cloud fields were sampled under a variety of environmental and aerosol conditions, with about 77% of the flights occurring in cumulus and stratocumulus. Preliminary analyses illustrate use of these data to analyze cloud-aerosol relationships, characterize the horizontal variability of cloud radiative impacts, and evaluate surface-based retrievals. We discuss how an extended-term campaign requires a simplified operating paradigm that is different from that used for typical, short-term, intensive aircraft field programs.

Physical interpretation of the spectral radiative signature in the transition zone between cloud-free and cloudy regions

One-second-resolution zenith radiance measure- ments from the Atmospheric Radiation Measurement pro- gram’s new shortwave spectrometer (SWS) provide a unique opportunity to analyze the transition zone between cloudy and cloud-free air, which has considerable bearing on the aerosol indirect effect. In the transition zone, we find a re- markable linear relationship between the sum and difference of radiances at 870 and 1640 nm wavelengths. The intercept of the relationship is determined primarily by aerosol prop- erties, and the slope by cloud properties. We then show that this linearity can be predicted from simple theoretical con- siderations and furthermore that it supports the hypothesis of inhomogeneous mixing, whereby optical depth increases as a cloud is approached but the effective drop size remains un- changed.