Remote Sensing Retrieval of Cloud Top Height Using Neural Networks and Data from Cloud-Aerosol Lidar with Orthogonal Polarization

In order to enhance the retrieval accuracy of cloud top height (CTH) from MODIS data, neural network models were employed based on Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data. Three types of methods were established using MODIS inputs: cloud parameters, calibrated radiance, and a combination of both. From a statistical standpoint, models with combination inputs demonstrated the best performance, followed by models with calibrated radiance inputs, while models relying solely on calibrated radiance had poorer applicability. This work found that cloud top pressure (CTP) and cloud top temperature played a crucial role in CTH retrieval from MODIS data. However, within the same type of models, there were slight differences in the retrieved results, and these differences were not dependent on the quantity of input parameters. Therefore, the model with fewer inputs using cloud parameters and calibrated radiance was recommended and employed for individual case studies. This model produced results closest to the actual cloud top structure of the typhoon and exhibited similar cloud distribution patterns when compared with the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) CTHs from a climatic statistical perspective. This suggests that the recommended model has good applicability and credibility in CTH retrieval from MODIS images. This work provides a method to improve accurate CTHs from MODIS data for better utilization.

Evaluation of CanESM Cloudiness, Cloud Type and Cloud Radiative Forcing Climatologies Using the CALIPSO-GOCCP and CERES Datasets

In this study, the annual and seasonal climatology of cloud fraction (CF) and cloud type simulated by the Canadian Environmental System Models (CanESMs) version 5 (CanESM5) and version 2 (CanESM2) at their fully coupled and AMIP configurations were validated against the CALIPSO-GOCCP-based CF. The CFs produced using the CALIPSO-COSP simulator based on the CanESMs data at their atmospheric (AMIP) configuration are also evaluated. The simulated shortwave, longwave, and net cloud radiative forcing using the AMIP version of the CanESM5 were also validated against satellite observations based on the recent CERES radiation satellite products. On average, all models have a negative bias in the total CF with global mean biases (MBs) of 2%, 2.4%, 3.9%, 6.4,%, 5.6%, and 7.1% for the coupled-CanESM5, AMIP-CanESM5, COSP-AMIP-CanESM5, coupled-CanESM2, AMIP-CanESM2, and COSP-AMIP-CanESM2, respectively, indicating that the CanESM5 has a smaller MB. There were no significant differences between AMIP and coupled versions of the model, but the COSP-based model-simulated data showed larger biases. Although the models captured well the climatological features of CF, they also exhibited a significant bias in CF reaching up to 40% over some geographical locations. This is particularly prevalent over the low level (LL) marine stratocumulus/cumulus, convectively active tropical latitudes that are normally dominated by high level (HL) clouds and at the polar regions where all models showed negative, positive, and positive bias corresponding to these locations, respectively. The AMIP-CanESM5 model performed reasonably well simulating the global mean cloud radiative forcing (CRF) with slight negative biases in the NetCRF at the TOA and surface that would be expected if the model has a positive bias in CF. This inconsistent result may be attributed to the parameterization of the optical properties in the model. The geographical distributions of the model bias in the NetCRF, however, can be significant reaching up to ±40 Wm−2 depending on the location and atmospheric level. The Pearson correlation showed that there is a strong correlation between the global distribution of model bias in NetCRF and CF and it is significantly influenced by the LL and HL clouds.

Comparison of the Spatial and Temporal Variability of Cloud Amounts over China Derived from Different Satellite Datasets

Various cloud cover products have been developed over the past few decades, but their uncertainties have not been sufficiently assessed, especially at a regional scale, which is vital for the application of satellite products to climate studies. In this study, we compare the spatial–temporal variability of the cloud amount over China from the 11 datasets provided by the Global Energy and Water Cycle Experiment (GEWEX) cloud assessment project at a horizontal resolution of 1° × 1° from the 1980s to 2000s, using the site data as a reference. The differences among these datasets are quantified in terms of the standard deviations and the correlation coefficients between different datasets. Most of the datasets show a similar spatial distribution of total cloud amounts (TCAs), but their magnitudes differ. The standard deviations of the annual, winter, and summer mean TCA are approximately 9–18% for the regional mean TCAs over the four typical regions of China, including the northwestern region (NW), northeastern region (NE), Tibetan Plateau region (TP), and southern China region (SC), with the largest standard deviations of 13–18% in the TP. By analyzing the factors that influence the satellite inversion data, such as the observation instrument, inversion algorithm, and observation time, we found that the difference caused by the observation instrument or algorithm is greater than the effect of the observation time, and the satellite cloud datasets with better recognition capability for cloud types show lower uncertainties when compared with the station observation. In terms of seasonal cycle, except HIRS and MODIS-ST, most satellite datasets can reproduce the observed seasonal cycle with the largest TCA in summer and the smallest TCA in autumn and winter. For the interannual variation, ISCCP-D1, MODIS-CE, and MODIS-ST are most consistent with the site data for the annual mean TCA, and two of the remaining datasets (PATMOSX and TOVSB) show more consistent temporal variations with the site observation in summer than in winter, especially over NW and NE regions. In general, MODIS-CE shows the best performance in reproducing the spatial pattern and interannual variation of TCA amongst the 11 satellite datasets, and PATMOSX, MODIS-ST, CALIPSO-GOCCP, and CALIPSO-ST also show relatively good performance.

Cloud-Top Height Comparison from Multi-Satellite Sensors and Ground-Based Cloud Radar over SACOL Site

Cloud-top heights (CTH), as one of the representative variables reflecting cloud macro-physical properties, affect the Earth–atmosphere system through radiation budget, water cycle, and atmospheric circulation. This study compares the CTH from passive- and active-spaceborne sensors with ground-based Ka-band zenith radar (KAZR) observations at the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) site for the period 2013–2019. A series of fundamental statistics on cloud probability in different limited time and areas at the SACOL site reveals that there is an optimal agreement for both cloud frequency and fraction derived from space and surface observations in a 0.5° × 0.5° box area and a 40-min time window. Based on the result, several facets of cloud fraction (CF), cloud overlapping, seasonal variation, and cloud geometrical depth (CGD) are investigated to evaluate the CTH retrieval accuracy of different observing sensors. Analysis shows that the CTH differences between multi-satellite sensors and KAZR decrease with increasing CF and CGD, significantly for passive satellite sensors in non-overlapping clouds. Regarding passive satellite sensors, e.g., Moderate Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua, the Multi-angle Imaging SpectroRadiometer (MISR) on Terra, and the Advanced Himawari Imager on Himawari-8 (HW8), a greater CTH frequency difference exists between the upper and lower altitude range, and they retrieve lower CTH than KAZR on average. The CTH accuracy of HW8 and MISR are susceptible to inhomogeneous clouds, which can be reduced by controlling the increase of CF. Besides, the CTH from active satellite sensors, e.g., Cloud Profiling Radar (CPR) on CloudSat, and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO), agree well with KAZR and are less affected by seasonal variation and inhomogeneous clouds. Only CALIPSO CTH is higher than KAZR CTH, mainly caused by the low-thin clouds, typically in overlapping clouds.

Analysis of Near-Cloud Changes in Atmospheric Aerosols Using Satellite Observations and Global Model Simulations

This paper examines cloud-related variations of atmospheric aerosols that occur in partly cloudy regions containing low-altitude clouds. The goal is to better understand aerosol behaviors and to help better represent the radiative effects of aerosols on climate. For this, the paper presents a statistical analysis of a multi-month global dataset that combines data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) satellite instruments with data from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) global reanalysis. Among other findings, the results reveal that near-cloud enhancements in lidar backscatter (closely related to aerosol optical depth) are larger (1) over land than ocean by 35%, (2) near optically thicker clouds by substantial amounts, (3) for sea salt than for other aerosol types, with the difference from dust reaching 50%. Finally, the study found that mean lidar backscatter is higher near clouds not because of large-scale variations in meteorological conditions, but because of local processes associated with individual clouds. The results help improve our understanding of aerosol-cloud-radiation interactions and our ability to represent them in climate models and other atmospheric models.

Measurements of Cloud Radiative Effect across the Southern Ocean (43° S–79° S, 63° E–158° W)

The surface radiation environment over the Southern Ocean within the region bound by 42.8° S to 78.7° S and 62.6° E to 157.7° W is summarised for three austral summers. This is done using ship-based measurements with the combination of downwelling radiation sensors and a cloud imager. We focus on characterising the cloud radiative effect (CRE) under a variety of conditions, comparing observations in the open ocean with those in the sea ice zone. For comparison with our observed data, we obtained surface data from the European Centre for Medium-Range Weather Forecasts fifth reanalysis (ERA5). We found that the daily average cloud fraction was slightly lower in ERA5 compared with the observations (0.71 and 0.75, respectively). ERA5 also showed positive biases in the shortwave radiation effect and a negative bias in the longwave radiation effect. The observed mean surface CRE of −164 ± 100 Wm−2 was more negative than the mean surface CRE for ERA5 of −101 W m−2.

Parameterization of Radiation Fog-Top Height and Methods Evaluation in Tianjin

Different methods have been developed to estimate the fog-top height of radiation fog and evaluated using the measurements obtained from a 255-m meteorological tower located in Tianjin in 2016. Different indicators of turbulence intensity, friction velocity (u*), turbulence kinetic energy (TKE), and variance of vertical velocity (σw2) were used to estimate the fog-top height, respectively. Positive correlations between the fog-top height and u*, TKE, and σw2 were observed, with empirical parameterization schemes H = 583.35 × u * 1.12 , H = 205.4   ×   ( T K E ) 0.68 , and H = 420.10 × ( σ w 2 ) 0.51 being obtained. Among them, σw2 is the most appropriate indicators of turbulence intensity to estimate the fog-top height. Compared with sensible flux and condensation rate, the new form of convective velocity scale (w*) was the most appropriate indicator of buoyancy induced by radiative cooling, and the relationship H = 328.33 × w * 1.34 was obtained. σw2 and with w*, which represents the intensity of turbulence and buoyancy, were used to estimate the fog-top height. The relationship H = 396.26 ×   (σw + 0.1 ×   w*) − 16 was obtained, which can be used to accurately estimate the fog-top height. Moreover, the temperature convergence (TC) method was used to estimate the fog-top height; however, the results strongly rely on the threshold value.

Cloud Products from the Earth Polychromatic Imaging Camera (EPIC): Algorithms and Initial Evaluation

This paper presents the physical basis of the EPIC cloud product algorithms and an initial evaluation of their performance. Since June 2015, EPIC has been providing observations of the sunlit side of the Earth with its 10 spectral channels ranging from the UV to the near-IR. A suite of algorithms has been developed to generate the standard EPIC Level 2 Cloud Products that include cloud mask, cloud effective pressure/height, and cloud optical thickness. The EPIC cloud mask adopts the threshold method and utilizes multichannel observations and ratios as tests. Cloud effective pressure/height is derived with observations from the O₂ A-band (780 nm and 764 nm), and B-band (680 nm and 688 nm) pairs. The EPIC cloud optical thickness retrieval adopts a single channel approach where the 780 nm and 680 nm channels are used for retrievals over ocean and over land, respectively. Comparison with co-located cloud retrievals from geosynchronous earth orbit (GEO) and low earth orbit (LEO) satellites shows that the EPIC cloud product algorithms are performing well and are consistent with theoretical expectations. These products are publicly available at the Atmospheric Science Data Center at the NASA Langley Research Center for climate studies and for generating other geophysical products that require cloud properties as input.

MISR-GOES 3D Winds: Implications for Future LEO-GEO and LEO-LEO Winds

Global wind observations are fundamental for studying weather and climate dynamics and for operational forecasting. Most wind measurements come from atmospheric motion vectors (AMVs) by tracking the displacement of cloud or water vapor features. These AMVs generally rely on thermal infrared (IR) techniques for their height assignments, which are subject to large uncertainties in the presence of weak or reversed vertical temperature gradients near the planetary boundary layer (PBL) and tropopause folds. Stereo imaging can overcome the height assignment problem using geometric parallax for feature height determination. In this study we develop a stereo 3D-Wind algorithm to simultaneously retrieve AMV and height from geostationary (GEO) and low Earth orbit (LEO) satellite imagery and apply it to collocated Geostationary Operational Environmental Satellite (GOES) and Multi-angle Imaging SpectroRadiometer (MISR) imagery. The new algorithm improves AMV and height relative to products from GOES or MISR alone, with an estimated accuracy of <0.5 m/s in AMV and <200 m in height with 2.2 km sampling. The algorithm can be generalized to other LEO-GEO or LEO-LEO combinations for greater spatiotemporal coverage. The technique demonstrated with MISR and GOES has important implications for future high-quality AMV observations, for which a low-cost constellation of CubeSats can play a vital role.

Observational Constraints on Cloud Feedbacks: The Role of Active Satellite Sensors

Cloud profiling from active lidar and radar in the A-train satellite constellation has significantly advanced our understanding of clouds and their role in the climate system. Nevertheless, the response of clouds to a warming climate remains one of the largest uncertainties in predicting climate change and for the development of adaptions to change. Both observation of long-term changes and observational constraints on the processes responsible for those changes are necessary. We review recent progress in our understanding of the cloud feedback problem. Capabilities and advantages of active sensors for observing clouds are discussed, along with the importance of active sensors for deriving constraints on cloud feedbacks as an essential component of a global climate observing system.

Overview of MPLNET Version 3 Cloud Detection

The National Aeronautics and Space Administration Micropulse Lidar Network Version 3 cloud detection algorithm is described and its differences relative to the previous version highlighted. Clouds are identified from normalized Level 1 signal profiles using two complementary methods. The first considers signal derivatives vertically for resolving low-level clouds. The second, which resolves high-level clouds like cirrus, is based on signal uncertainties given the relatively low signal-to-noise ratio exhibited in the upper troposphere by eye-safe network instruments, especially during daytime. Furthermore, a multi-temporal averaging scheme is used to improve cloud detection under conditions of weak signal-to-noise. Diurnal and seasonal cycles of cloud occurrence frequency based on one year of measurements at the Goddard Space Flight Center (Greenbelt, MD) site are compared for the new and previous versions. The largest differences, and perceived improvement, in detection occurs for high clouds (above 5-km, mean sea level) which increase in occurrence by nearly 6%. There is also an increase in the detection of multi-layered cloud profiles from 9% to 20%. Macrophysical properties and estimates of cloud optical depth are presented for a transparent cirrus dataset. However, the limit to which molecular signal can be reliably retrieved above cirrus clouds occurs between cloud optical depths of 0.5 and 0.8.

Clouds at Barbados are representative of clouds across the trade wind regions in observations and climate models

Trade wind regions cover most of the tropical oceans, and the prevailing cloud type is shallow cumulus. These small clouds are parameterized by climate models, and changes in their radiative effects strongly and directly contribute to the spread in estimates of climate sensitivity. This study investigates the structure and variability of these clouds in observations and climate models. The study builds upon recent detailed model evaluations using observations from the island of Barbados. Using a dynamical regimes framework, satellite and reanalysis products are used to compare the Barbados region and the broader tropics. It is shown that clouds in the Barbados region are similar to those across the trade wind regions, implying that observational findings from the Barbados Cloud Observatory are relevant to clouds across the tropics. The same methods are applied to climate models to evaluate the simulated clouds. The models generally capture the cloud radiative effect, but underestimate cloud cover and show an array of cloud vertical structures. Some models show strong biases in the environment of the Barbados region in summer, weakening the connection between the regional biases and those across the tropics. Even bearing that limitation in mind, it is shown that covariations of cloud and environmental properties in the models are inconsistent with observations. The models tend to misrepresent sensitivity to moisture variations and inversion characteristics. These model errors are likely connected to cloud feedback in climate projections, and highlight the importance of the representation of shallow cumulus convection.

Uncertainties in cloud phase and optical thickness retrievals from the Earth Polychromatic Imaging Camera (EPIC)

This paper presents an investigation of the expected uncertainties of a single channel cloud optical thickness (COT) retrieval technique, as well as a simple cloud temperature threshold based thermodynamic phase approach, in support of the Deep Space Climate Observatory (DSCOVR) mission. DSCOVR cloud products will be derived from Earth Polychromatic Imaging Camera (EPIC) observations in the ultraviolet and visible spectra. Since EPIC is not equipped with a spectral channel in the shortwave or mid-wave infrared that is sensitive to cloud effective radius (CER), COT will be inferred from a single visible channel with the assumption of appropriate CER values for liquid and ice phase clouds. One month of Aqua MODIS daytime granules from April 2005 is selected for investigating cloud phase sensitivity, and a subset of these granules that has similar EPIC sun-view geometry is selected for investigating COT uncertainties. EPIC COT retrievals are simulated with the same algorithm as the operational MODIS cloud products (MOD06), except using fixed phase-dependent CER values. Uncertainty estimates are derived by comparing the single channel COT retrievals with the baseline bi-spectral MODIS retrievals. Results show that a single channel COT retrieval is feasible for EPIC. For ice clouds, single channel retrieval errors are minimal (< 2%) due to the particle size insensitivity of the assumed ice crystal (i.e., severely roughened aggregate of hexagonal columns) scattering properties at visible wavelengths, while for liquid clouds the error is mostly limited to within 10%, although for thin clouds (COT < 2) the error can be higher. Potential uncertainties in EPIC cloud masking and cloud temperature retrievals are not considered in this study.