Search NASA⌕ Search

SEARCH · Search NASA

Results for “cloud inhomogeneity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Liquid Water Cloud Measurements Using the Raman Lidar Technique: Current Understanding and Future Research Needs

This paper describes recent work in the Raman lidar liquid water cloud measurement technique. The range-resolved spectral measurements at the National Aeronautics and Space Administration Goddard Space Flight Center indicate that the Raman backscattering spectra measured in and below low clouds agree well with theoretical spectra for vapor and liquid water. The calibration coefficients of the liquid water measurement for the Raman lidar at the Atmospheric Radiation Measurement Program Southern Great Plains site of the U.S. Department of Energy were determined by comparison with the liquid water path (LWP) obtained with Atmospheric Emitted Radiance Interferometer (AERI) and the liquid water content (LWC) obtained with the millimeter wavelength cloud radar and water vapor radiometer (MMCR-WVR) together. These comparisons were used to estimate the Raman liquid water cross-sectional value. The results indicate a bias consistent with an effective liquid water Raman cross-sectional value that is 28%-46% lower than published, which may be explained by the fact that the difference in the detectors' sensitivity has not been accounted for. The LWP of a thin altostratus cloud showed good qualitative agreement between lidar retrievals and AERI. However, the overall ensemble of comparisons of LWP showed considerable scatter, possibly because of the different fields of view of the instruments, the 350-m distance between the instruments, and the horizontal inhomogeneity of the clouds. The LWC profiles for a thick stratus cloud showed agreement between lidar retrievals andMMCR-WVR between the cloud base and 150m above that where the optical depth was less than 3. Areas requiring further research in this technique are discussed.

vapor and liquid water↗
DOE OSTI2021

Observations of Aerosol, Cloud, Turbulence, and Radiation Properties at the Top of the Marine Boundary Layer over the Eastern North Atlantic Ocean: The ACORES Campaign

We report on the Azores Stratocumulus Measurements of Radiation, Turbulence and Aerosols (ACORES) campaign, which took place around Graciosa and Pico Islands/Azores in July 2017. The main objective was to investigate the vertical distribution of aerosol particles, stratocumulus microphysical and radiative properties, and turbulence parameters in the eastern North Atlantic. The vertical exchange of mass, momentum, and energy between the free troposphere (FT) and the cloudy marine boundary layer (MBL) was explored over a range of scales from submeters to kilometers. To cover these spatial scales with appropriate measurements, helicopter-borne observations with unprecedented high resolution were realized using the Airborne Cloud Turbulence Observation System (ACTOS) and Spectral Modular Airborne Radiation Measurement System–Helicopter-Borne Observations (SMART-HELIOS) instrumental payloads. Here, the helicopter-borne observations were combined with ground-based aerosol measurements collected at two continuously running field stations on Pico Mountain (2,225 m above sea level, in the FT), and at the Atmospheric Radiation Measurement (ARM) station on Graciosa (at sea level). First findings from the ACORES observations we are discussing in the paper are as follows: (i) we have observed a high variability of the turbulent cloud-top structure on horizontal scales below 100 m with local temperature gradients of up to 4 K over less than 1 m vertical distance, (ii) we have collected strictly collocated radiation measurements supporting the relevance of small-scale processes by revealing significant inhomogeneities in cloud-top brightness temperature to scales well below 100 m, and (iii) we have concluded that aerosol properties are completely different in the MBL and FT with often-complex stratification and frequently observed burst-like new particle formation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Low-level liquid cloud properties during ORACLES retrieved using airborne polarimetric measurements and a neural network algorithm

In this study we developed a neural network (NN) that can be used to retrieve cloud microphysical properties from multiangular and multispectral polarimetric remote sensing observations. This effort builds upon our previous work, which explored the sensitivity of neural network input, architecture, and other design requirements for this type of remote sensing problem. In particular this work introduces a framework for appropriately weighting total and polarized reflectances, which have vastly different magnitudes and measurement uncertainties. The NN is trained using an artificial training set and applied to research scanning polarimeter (RSP) data obtained during the ORACLES field campaign (ObseRvations of Aerosols above CLouds and their intEractionS). The polarimetric RSP observations are unique in that they observe the same cloud from a very large number of angles within a variety of spectral bands, resulting in a large dataset that can be explored rapidly with a NN approach. The usefulness of applying a NN to a dataset such as this one stems from the possibility of rapidly obtaining a retrieval that could be subsequently applied as a first guess for slower but more rigorous physical-based retrieval algorithms. This approach could be particularly advantageous for more complicated atmospheric retrievals – such as when an aerosol layer lies above clouds like in ORACLES. For RSP observations obtained during ORACLES 2016, comparisons between the NN and standard parametric polarimetric (PP) cloud retrieval give reasonable results for droplet effective radius (r(e): R=0.756, RMSE=1.74 µm) and cloud optical thickness (τ: R=0.950, RMSE=1.82). This level of statistical agreement is shown to be similar to comparisons between the two most well-established cloud retrievals, namely, the polarimetric and the bispectral total reflectance cloud retrievals. The NN retrievals from the ORACLES 2017 dataset result in retrievals of r(e) (R=0.54, RMSE=4.77 µm) and τ (R=0.785, RMSE=5.61) that behave much more poorly. In particular we found that our NN retrieval approach does not perform well for thin (τ<3), inhomogeneous, or broken clouds. We also found that correction for above-cloud atmospheric absorption improved the NN retrievals moderately – but retrievals without this correction still behaved similarly to existing cloud retrievals with a slight systematic offset.

SPECTRAL-RESOLUTION LIDAR↗
NASA NTRS2016

Observation of the Spectrally Invariant Properties of Clouds in Cloudy-to-Clear Transition Zones During the MAGIC Field Campaign

We use the spectrally invariant method to study the variability of cloud optical thickness τau and droplet effective radius r(sub eff) in transition zones (between the cloudy and clear sky columns) observed from Solar Spectral Flux Radiometer (SSFR) and Shortwave Array Spectroradiometer-Zenith (SASZe) during the Marine ARM GPCI Investigation of Clouds (MAGIC) field campaign. The measurements from the SSFR and the SASZe are different, however inter-instrument differences of self-normalized measurements (divided by their own spectra at a fixed time) are small. The spectrally invariant method approximates the spectra in the cloud transition zone as a linear combination of definitely clear and cloudy spectra, where the coefficients, slope and intercept, characterize the spectrally invariant properties of the transition zone. Simulation results from the SBDART (Santa Barbara DISORT Atmospheric Radiative Transfer) model demonstrate that (1) the slope of the visible band is positively correlated with the cloud optical thickness τ while the intercept of the near-infrared band has high negative correlation with the cloud drop effective radius r(sub eff)even without the exact knowledge of τau; (2) the above relations hold for all Solar Zenith Angle (SZA) and for cloud-contaminated skies. In observations using redundant measurements from SSFR and SASZe, we find that during cloudy-to-clear transitions, (a) the slopes of the visible band decrease, and (b) the intercepts of the near-infrared band remain almost constant near cloud edges. The findings in simulations and observations suggest that, while the optical thickness decreases during the cloudy-to-clear transition, the cloud drop effective radius does not change when cloud edges are approached. These results support the hypothesis that inhomogeneous mixing dominates near cloud edges in the studied cases.

spectral invariance↗
NASA NTRS2013

CALIOP and AERONET Aerosol Optical Depth Comparisons: One Size Fits None

We compare the aerosol optical depths (AOD) retrieved from backscatter measurements of the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aboard the Cloud Aerosol Lidar Infrared Pathfinder Satellite Observations (CALIPSO) satellite with coincident Aerosol Robotic Network (AERONET) measurements. Overpass coincidence criteria of +/- 2 h and within a 40 km radius are satisfied at least once at 149 globally distributed AERONET sites from 2006 to 2010. Most data pairs (>80%) use AERONET measurements acquired +/- 30 min of the overpass. We examine the differences in AOD estimates between CALIOP and AERONET for various aerosol, environmental, and geographic conditions. Results show CALIOP AOD are lower than AERONET AOD especially at low optical depths as measured by AERONET (500 nm AOD<0.1). Furthermore, the median relative AOD difference between the two measurements is 25% of the AERONET AOD for AOD>0.1. Differences in AOD between CALIOP and AERONET are possibly due to cloud contamination, scene inhomogeneity, instrument view angle differences, CALIOP retrieval errors, and detection limits. Comparison of daytime to nighttime number of 5 km 60m (60m in the vertical) features detected by CALIOP show that there are 20% more aerosol features at night. We find that CALIPSO and AERONET do not agree on the cloudiness of scenes. Of the scenes that meet the above coincidence criteria, CALIPSO finds clouds in more than 45% of the coincident atmospheric columns AERONET classifies as clear.

Omar, A. H.↗
DOE OSTI2021

Remote Sensing of Electric Fields Observed Within Winter Precipitation During the 2020 Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) Field Campaign

Aircraft electric fields from NASA's Lightning Instrument Package (LIP) were coupled with other airborne and ground-based, and in situ measurements to understand electrification in winter clouds that did not produce lightning. The measurements were made during seven research flights by a NASA ER-2 during the 2020 Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. Observed total electric field magnitudes were as high as 80 V m -1 and variability in the electric field was observed along the flight path of the ER-2, indicating horizontal and/or vertical inhomogeneity in the cloud's electrical structure. X-band airborne radar data indicated 20-dBZ echo tops above 5 km in regions where electrification exceeded 10 V m -1 . In these regions, 85-GHz brightness temperatures (TB) from an airborne radiometer were lower than 265 K, with the lowest T B (~210 K) associated with ice scattering collocated with the strongest electric fields. In situ microphysical measurements from the NASA P-3 aircraft on February 7 indicated that regions near strong electric field contained supercooled water, rimed ice hydrometeors, ice water p- content as high as 1 g m -3 , liquid water content as high as 0.15 g m -3 , and supersaturation as high as 3.5%. Furthermore, these observations support the role of mixed phase microphysics in the generation of electric fields in clouds. In three case studies, ground based S-band polarimetric radar observed depolarization streaks in differential reflectivity near areas where the strongest electrification was observed. This observation reinforces the utility of depolarization streaks to identify areas of electrification prior to lightning occurrence.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Multiple scattering depolarization in marine stratus clouds: Lidar experiments

The depolarization of ruby lidar backscattering caused by multiple scattering in marine stratus clouds was examined systematically from a field site on the southern California coast. Investigated were the effects on the linear depolarization (delta) of lidar receiver field of view (FOV), elevation angle and laser beam pointing errors. An approximately linear increase in maximum delta values was observed with increasing receiver FOV, and the importance of accurate transmitter/receiver beam alignment was demonstrated during experiments in which the laser axis was deliberately misaligned. An elevation angle dependence to the delta values was observed as a consequence of the natural vertical inhomogeneity of water cloud content above the cloud base. Time histories of the depolarization characteristics of dissipating stratus clouds revealed significant spatial and temporal variability in delta values attributed to cloud composition variations. Employing a 1 mrad transmitter FOV, maximum delta values of 0.21 nd 0.33 were observed with 1 and 3 mrad receiver FOVs, respectively, from the low stratus clouds. The fundamental causes and effects on the lidar equation of multiple scattering are also discussed.

Sassen, K.↗
NASA NTRS1998

Monte Carlo Calculations of Polarized Microwave Radiation Emerging from Cloud Structures

The last decade has seen tremendous growth in cloud dynamical and microphysical models that are able to simulate storms and storm systems with very high spatial resolution, typically of the order of a few kilometers. The fairly realistic distributions of cloud and hydrometeor properties that these models generate has in turn led to a renewed interest in the three-dimensional microwave radiative transfer modeling needed to understand the effect of cloud and rainfall inhomogeneities upon microwave observations. Monte Carlo methods, and particularly backwards Monte Carlo methods have shown themselves to be very desirable due to the quick convergence of the solutions. Unfortunately, backwards Monte Carlo methods are not well suited to treat polarized radiation. This study reviews the existing Monte Carlo methods and presents a new polarized Monte Carlo radiative transfer code. The code is based on a forward scheme but uses aliasing techniques to keep the computational requirements equivalent to the backwards solution. Radiative transfer computations have been performed using a microphysical-dynamical cloud model and the results are presented together with the algorithm description.

Kummerow, Christian↗
NASA NTRS1989

Fractal statistics of cloud fields

Landsat Multispectral Scanner (MSS) and Thematic Mapper (TM) data, with 80 and 30 m spatial resolution, respectively, have been employed to study the spatial structure of boundary-layer and intertropical convergence zone (ITCZ) clouds. The probability distributions of cloud areas and cloud perimeters are found to approximately follow a power-law, with a different power (i.e., fractal dimension) for each cloud type. They are better approximated by a double power-law behavior, indicating a change in the fractal dimension at a characteristic size which depends upon cloud type. The fractal dimension also changes with threshold. The more intense cloud areas are found to have a higher perimeter fractal dimension, perhaps indicative of the increased turbulence at cloud top. A detailed picture of the inhomogeneous spatial structure of various cloud types will contribute to a better understanding of basic cloud processes, and also has implications for the remote sensing of clouds, for their effects on remote sensing of other parameters, and for the parameterization of clouds in general circulation models, all of which rely upon plane-parallel radiative transfer algorithms.

Cahalan, Robert F.↗
NASA NTRS1994

The giant molecular cloud Monoceros R2. 1: Shell structure

We have obtained a 45 sec resolution, Nyquist-sampled map in CO J = 1-0 covering approximately a 3 deg x 3 deg region of the giant molecular cloud Monoceros R2. The map consists of 167,000 spectra observed with the 15 element focal-plane array system on the FCRAO 14 m telescope. The data reveal that the large-scale structure of Mon R2 is dominated by a is approximately 30 pc diameter largely hemispherical shell containing approximately 4 x 10(exp 4) solar mass of molecular material and expanding at approximately 3-4 km s(exp -1) with symmetric axis roughly along the line of sight. The dynamical timescale of the shell is estimated to be approximately 4 x 10(exp 6) yr, which is consistent with the age of main-sequence stars powering the clusters of reflection nebulea in this region. There is no evidence for a redshifted shell on the far side of the interior 'bubble,' which is largely devoid of molecular material. Distortions of the shell are obvious, suggesting inhomogeneity of the cloud and possible presence of a magnetic field prior to its formation. Dense clumps in Mon R2, including the main core and the GGD 12-15 core, appear to be condensations located on the large shell. The reflection nebulea with their illuminating stars as well as embedded IRAS sources suggest that triggered star formation has taken place over a large part of the Mon R2 shell.

Xie, Taoling↗
NASA NTRS2010

Effects of Ice Particle Size Vertical Inhomogeneity on the Passive Remote Sensing of Ice Clouds

The solar reflectance bi-spectral (SRBS) and infrared split-window (IRSpW) methods are two of the most popular techniques for passive ice cloud property retrievals from multispectral imagers. Ice clouds are usually assumed to be vertically homogeneous in global operational algorithms based on these methods, although significant vertical variations of ice particle size are typically observed in ice clouds. In this Study we investigate uncertainties in retrieved optical thickness, effective particle size, and ice water path introduced by a homogeneous cloud assumption in both the SRBS and IRSpW methods, and focus on whether the assumption can lead to significant discrepancies between the two methods. The study simulates the upwelling spectral radiance associated with vertically structured clouds and passes the results through representative SRBS and IRSpW retrieval algorithms. Cloud optical thickness is limited to values for which IRSpW retrievals are possible (optical thickness less than about 7). When the ice cloud is optically thin and yet has a significant ice particle size vertical variation, it is found that both methods tend to underestimate the effective radius and ice water path. The reason for the underestimation is the nonlinear dependence of ice particle scattering properties (extinction and single scattering albedo) on the effective radius. Because the nonlinearity effect is stronger in the IRSpW than the SRBS method, the IRSpW-based IWP tends to be smaller than the SRBS counterpart. When the ice cloud is moderately optically thick, the IRSpW method is relatively insensitive to cloud vertical structure and effective radius retrieval is weighted toward smaller ice particle size, while the weighting function makes the SRBS method more sensitive to the ice particle size in the upper portion of the cloud. As a result, when ice particle size increases monotonically toward cloud base, the two methods are in qualitative agreement; in the event that ice particle size decreases toward cloud base, the effective radius and ice water path retrievals based on the SRBS method are substantially larger than those from the IRSpW. The main findings of this Study Suggest that the homogenous cloud assumption can affect the SRBS and IRSpW methods to different extents and, consequently, can lead to significantly different retrievals. Therefore caution should be taken when comparing and combining the ice cloud property retrievals from these two methods.

Zhang, Zhibo↗
NASA NTRS2004

The Effect of Surface Heterogeneity on Cloud Absorption Estimates

This study presents a systematic and quantitative analysis of the effect of inhomogeneous surface albedo on shortwave cloud absorption estimates. We use 3D radiative transfer modeling with gradually complex clouds over a simplified surface to calculate cloud absorption. We find that averaging surface albedo always underestimates cloud absorption, and thus accounting for surface heterogeneity always enhances cloud absorption. However, the impact on cloud absorption estimates is not enough to explain the discrepancy between measured and model calculated shortwave cloud absorptions.

Chiu, Jui-Yuan C.↗
NASA NTRS1983

The influence of rain and clouds on a satellite dual frequency radar altimeter system operating at 13 and 35 GHz

The effects of inhomogeneous spatial attenuation resulting from clouds and rain on the altimeter estimate of the range to mean sea level are modelled. It is demonstrated that typical cloud and rain attenuation variability at commonly expected spatial scales can significantly degrade altimeter range precision. Rain cell and cloud scale sizes and attenuations are considered as factors. The model simulation of altimeter signature distortion is described, and the distortion of individual radar pulse waveforms by different spatial scales of attenuation is considered. Examples of range errors found for models of a single cloud, a rain cell, and cloud streets are discussed.

Walsh, E. J.↗
DOE OSTI2026

Generating a stratocumulus-like cloud top in a convection-cloud chamber

Stratocumulus-topped boundary layers play a crucial role in influencing daily weather and earth energy balance. Entrainment at the stratocumulus cloud top affects the cloud’s lifetime, precipitation, and radiative properties, but our understanding remains limited due to the lack of resolution in both field observations and numerical simulations. A recently proposed convection-cloud chamber with detailed control of sidewall temperatures can provide a unique opportunity to explore this mechanism in a laboratory setting. In this work, we use numerical simulations to demonstrate that this design can produce a cloud top that mimics the entrainment interfacial layer in a stratocumulus cloud. Our results show that a steady-state cloud can be formed by cooling the lower portions of the sidewalls and warming the bottom surface, while a temperature inversion at the cloud top can be generated by keeping the upper sidewalls and top surface warmer than the bottom. The turbulent kinetic energy profile and budget are similar to those found in a convective boundary layer, and inhomogeneous mixing near the cloud top can be observed. These findings significantly enhance the scientific value of constructing the tall convection-cloud chamber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High-resolution lidar observations of sedimentation-induced size sorting of droplets near a laboratory cloud top

Cloud optical properties and precipitation, which are crucial to weather and climate, are strongly influenced by cloud microphysical properties that are still poorly understood. Here, we develop a high-resolution time-correlated single-photon-counting lidar and apply it to observe cloud microphysical properties at one-centimeter range resolution in a convection chamber under well-controlled conditions. Together with concurrent in-situ measurements and theoretical analysis, our lidar observations indicate that although turbulent mixing tends to homogenize the cloud in the bulk region, entrainment and sedimentation cause inhomogeneities in droplet concentrations near the cloud top. Specifically, the topmost region is directly affected by entrainment, and lidar profiles show clear evidence of entrained air and detrained cloud filament. The transition region below exhibits vertical size sorting of cloud droplets caused by sedimentation. Our results suggest that using a single sedimentation velocity for all cloud droplets, as is done in many atmospheric models, overlooks key physics relevant to the microphysical structure near the cloud top. In conclusion, our conceptual model used to describe these measurements can serve as a step toward improving the current modeling of processes in the cloud top region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Polarimetric Retrievals of Surface and Cirrus Clouds Properties in the Region Affected by the Deepwater Horizon Oil Spill

In 2010, the Goddard Institute for Space Studies (GISS) Research Scanning Polarimeter (RSP) performed several aerial surveys over the region affected by the oil spill caused by the explosion of the Deepwater Horizon offshore platform. The instrument was deployed on the NASA Langley B200 aircraft together with the High Spectral Resolution Lidar (HSRL), which provides information on the distribution of the aerosol layers beneath the aircraft, including an accurate estimate of aerosol optical depth. This work illustrates the merits of polarization measurements in detecting variations of ocean surface properties linked to the presence of an oil slick. In particular, we make use of the degree of linear polarization in the glint region, which is severely affected by variations in the refractive index but insensitive to the waviness of the water surface. Alterations in the surface optical properties are therefore expected to directly affect the polarization response of the RSP channel at 2264 nm, where both molecular and aerosol scattering are negligible and virtually all of the observed signal is generated via Fresnel reflection at the surface. The glint profile at this wavelength is fitted with a model which can optimally estimate refractive index, wind speed and direction, together with aircraft attitude variations affecting the viewing geometry. The retrieved refractive index markedly increases over oil-contaminated waters, while the apparent wind speed is significantly lower than in adjacent uncontaminated areas, suggesting that the slick dampens high-frequency components of the ocean wave spectrum. The constraint on surface reflectance provided by the short-wave infrared channels is a cornerstone of established procedures to retrieve atmospheric aerosol microphysical parameters based on the inversion of the RSP multispectral measurements. This retrieval, which benefits from the ancillary information provided by the HSRL, was in this specific case hampered by prohibitive variability in atmospheric conditions (very inhomogeneous aerosol distribution and cloud cover). Although the results presented for the surface are essentially unaffected, we discuss the results obtained by typing algorithms in sorting the complex mix of aerosol types, and show evidence of oriented ice in cirrus clouds present in the area. In this context, polarization measurements at 1880 nm were used to infer ice habit and cirrus optical depth, which was found in the subvisual/threshold-visible regime, confirming the utility of the aforementioned RSP channel for the remote sensing of even thin cold clouds.

Polarization↗
NASA NTRS2001

Atomic Oxygen Abundance in Molecular Clouds: Absorption Toward Sagittarius B2

We have obtained high-resolution (approximately 35 km/s) spectra toward the molecular cloud Sgr B2 at 63 micrometers, the wavelength of the ground-state fine-structure line of atomic oxygen (O(I)), using the ISO-LWS instrument. Four separate velocity components are seen in the deconvolved spectrum, in absorption against the dust continuum emission of Sgr B2. Three of these components, corresponding to foreground clouds, are used to study the O(I) content of the cool molecular gas along the line of sight. In principle, the atomic oxygen that produces a particular velocity component could exist in any, or all, of three physically distinct regions: inside a dense molecular cloud, in the UV illuminated surface layer (PDR) of a cloud, and in an atomic (H(I)) gas halo. For each of the three foreground clouds, we estimate, and subtract from the observed O(I) column density, the oxygen content of the H(I) halo gas, by scaling from a published high-resolution 21 cm spectrum. We find that the remaining O(I) column density is correlated with the observed (13)CO column density. From the slope of this correlation, an average [O(I)]/[(13)CO] ratio of 270 +/- 120 (3-sigma) is derived, which corresponds to [O(I)]/[(13)CO] = 9 for a CO to (13)CO abundance ratio of 30. Assuming a (13)CO abundance of 1x10(exp -6) with respect to H nuclei, we derive an atomic oxygen abundance of 2.7x10(exp -4) in the dense gas phase, corresponding to a 15% oxygen depletion compared to the diffuse ISM in our Galactic neighborhood. The presence of multiple, spectrally resolved velocity components in the Sgr B2 absorption spectrum allows, for the first time, a direct determination of the PDR contribution to the O(I) column density. The PDR regions should contain O(I) but not (13)CO, and would thus be expected to produce an offset in the O(I)-(13)CO correlation. Our data do not show such an offset, suggesting that within our beam O(I) is spatially coexistent with the molecular gas, as traced by (13)CO. This may be a result of the inhomogeneous nature of the clouds.

Lis, D. C.↗