Simulation of Mesoscale Cellular Convection in Marine Stratocumulus. Part I: Drizzling Conditions
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Drizzle initiation through the collision and coalescence of cloud droplets plays a crucial role in warm cloud precipitation. Recent theoretical studies suggest that the influence of collisional growth on the droplet size distribution can be quantified by a non-dimensional drizzle number (Dz). Here, large-eddy simulations with Lagrangian microphysics are employed to evaluate the theory by simulating a tall convection-cloud chamber under various conditions. Results show that the smaller the Dz, the larger the impact of collisions on the right tail of the droplet size distribution, consistent with the theory. The simulations confirm that the collision rate can be estimated from the droplet size distribution interacting only with cloud droplets of the same size at the mode radius. This suggests that the idealized theory can be a useful tool to design a cloud chamber for drizzle investigation, as well as to represent drizzle formation in models of real atmospheric clouds.
Drizzle, a common feature of marine boundary layer clouds formed through collision coalescence, plays a key role in cloud microphysics and evolution. Yet, simultaneously retrieving cloud and drizzle properties from remote-sensing observations remains challenging because drizzle droplets often dominate radar signals, masking cloud contributions. The goal of the proposed research is to provide constraints for the process of autoconversion and accretion using ARM cloud measurements. Specifically, we provide concurrent retrievals of cloud and drizzle that allows users to derive corresponding autoconversion and accretion rates.
In this study, aircraft observations from the Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) Twin Otter are used to characterize the variability in drizzle, cloud, and aerosol properties associated with cloud rifts and the surrounding solid clouds observed off the coast of California. A flight made on 16 July 1999 provided measurements directly across an interface between solid and rift cloud conditions. Aircraft instrumentation allowed for measurements of aerosol, cloud droplet, and drizzle spectra. CCN concentrations were measured in addition to standard thermodynamic variables and the winds. A Forward Scatter Spectrometer Probe (FSSP) measured size distribution of cloud-sized droplets. A Cloud Imaging Probe (CIP) was used to measure distributions of drizzle-sized droplets. Aerosol distributions were obtained from a Cloud Aerosol Scatterprobe (CAS). The CAS probe measured aerosols, cloud droplets and drizzle-sized drops; for this study. The CAS probe was used to measure aerosols in the size range of 0.5 micron - 1 micron. Smaller aerosols were characterized using an Ultrafine Condensation Particle Counter (CPC) sensor. The CPC was used to measure particles with diameters greater than 0.003 micron. By subtracting different count concentrations measured with the CPC, this probe was capable of identifying ultrafine particles those falling in the size range of 3 nanometers - 7 nanometers that are believed to be associated with new particle production.
Drizzle is a common feature of warm stratiform clouds and it influences their radiative effects by modulating their physical properties and lifecycle. An important component of drizzle formation are processes that lead to a broadening of the droplet size distribution (DSD). Here, we examine observations of cloud and drizzle properties retrieved using colocated airborne measurements from the Research Scanning Polarimeter and the Third Generation Airborne Precipitation Radar. We observe a bimodal DSD as the aircraft transects drizzling open-cells whereby the larger mode reaches a maximum size near cloud center and the smaller mode remains relatively constant in size. We review similarities between our observations with droplet growth processes and their connections with precipitation onset. We estimate droplet sedimentation using the cloud top DSD and find a correlation with rain water path of 0.82. We also examine how changes in liquid water paths and droplet concentrations may act to enhance or suppress precipitation.
Vertically resolved thermodynamic cloud-phase classifications are essential for studies of atmospheric cloud and precipitation processes. The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Thermodynamic Cloud Phase (THERMOCLDPHASE) value-added product (VAP) uses a multi-sensor approach to classify the thermodynamic cloud phase by combining lidar backscatter and depolarization, radar reflectivity, Doppler velocity, spectral width, microwave-radiometer-derived liquid water path, and radiosonde temperature measurements. The measured pixels are classified as ice, snow, mixed phase, liquid (cloud water), drizzle, rain, and liq_driz (liquid+drizzle). We use this product as the ground truth to train three machine learning (ML) models to predict the thermodynamic cloud phase from multi-sensor remote sensing measurements taken at the ARM North Slope of Alaska (NSA) observatory: a random forest (RF), a multi-layer perceptron (MLP), and a convolutional neural network (CNN) with a U-Net architecture. Evaluations against the outputs of the THERMOCLDPHASE VAP with 1 year of data show that the CNN outperforms the other two models, achieving the highest test accuracy, F1 score, and mean intersection over union (IOU). Analysis of ML confidence scores shows that ice, rain, and snow have higher confidence scores, followed by liquid, while mixed, drizzle, and liq_driz have lower scores. Feature importance analysis reveals that the mean Doppler velocity and vertically resolved temperature are the most influential data streams for ML thermodynamic cloud-phase predictions. Lidar measurements exhibit lower feature importance due to rapid signal attenuation caused by the frequent presence of persistent low-level clouds at the NSA site. The ML models' generalization capacity is further evaluated by applying them at another Arctic ARM site in Norway using data taken during the ARM Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) field campaign. The models demonstrated similar performance to that observed at the NSA site. Finally, we evaluate the ML models' response to simulated instrument outages and signal degradation and show that a CNN U-Net model trained with input channel dropouts performs better when input fields are missing.
Vertically resolved thermodynamic cloud phase classifications are essential for studies of atmospheric cloud and precipitation processes. The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) THERMOCLDPHASE Value-Added Product (VAP) uses a multi-sensor approach to classify thermodynamic cloud phase by combining lidar backscatter and depolarization, radar reflectivity, Doppler velocity, spectral width, microwave radiometer-derived liquid water path, and radiosonde temperature measurements. The measured voxels are classified as ice, snow, mixed-phase, liquid (cloud water), drizzle, rain, and liq_driz (liquid+drizzle). We use this product as the ground truth to train three machine learning (ML) models to predict the thermodynamic cloud phase from multi-sensor remote sensing measurements taken at the ARM North Slope of Alaska (NSA) observatory: a random forest (RF), a multilayer perceptron (MLP), and a convolutional neural network (CNN) with a U-Net architecture. Evaluations against the outputs of the THERMOCLDPHASE VAP with one year of data show that the CNN outperforms the other two models, achieving the highest test accuracy, F1-score, and mean Intersection over Union (IOU). Analysis of ML confidence scores shows ice, rain, and snow have higher confidence scores, followed by liquid, while mixed, drizzle, and liq_driz have lower scores. Feature importance analysis reveals that the mean Doppler velocity and vertically resolved temperature are the most influential datastreams for ML thermodynamic cloud phase predictions. The ML models’ generalization capacity is further evaluated by applying them at another Arctic ARM site in Norway using data taken during the ARM Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) field campaign. Finally, we evaluate the ML models’ response to simulated instrument outages and signal degradation.
A series of test flights was conducted by the U. S. Navy over a 3- year period to evaluate the effects of icing on the operation of the ZPG-2 airship. In supercooled. clouds, ice formed only on the forward edges of small protuberances and wires and presented no serious hazard to operation. Ice accretions of the glaze type which occurred in conditions described as freezing drizzle adversely affected various components to a somewhat greater extent. The results indicated, a need for protection of certain components such as antennas, propellers, and certain parts of the control system. The tests showed that icing of the large surface of the envelope occurred only in freezing rain or drizzle. Because of the infrequent occurrence of these conditions, the potential maximum severity could not be estimated from the test results. The increases in heaviness caused by icing in freezing rain and drizzle were substantial, but well within the operational capabilities of the airship. In order to estimate the potential operational significance of icing in freezing rain, theoretical calculations were used to estimate: (1) the rate of icing as a function of temperature and rainfall intensity, (2) the climatological probability of occurrence of various combinations of these variables, and (3) the significance of the warming influence of the ocean in alleviating freezing-rain conditions. The results of these calculations suggest that, although very heavy icing rates are possible in combinations of low temperature and high rainfall rate, the occurrence of such conditions is very infrequent in coastal areas and virtually impossible 200 or 300 miles offshore.
A two-layer model of the marine boundary layer is described. The model is used to simulate both stratocumulus and shallow cumulus clouds in downstream simulations. Over cold sea surfaces, the model predicts a relatively uniform structure in the boundary layer with 90%-100% cloud fraction. Over warm sea surfaces, the model predicts a relatively strong decoupled and conditionally unstable structure with a cloud fraction between 30% and 60%. A strong large-scale divergence considerably limits the height of the boundary layer and decreases relative humidity in the upper part of the cloud layer; thus, a low cloud fraction results. The efffects of drizzle on the boundary-layer structure and cloud fraction are also studied with downstream simulations. It is found that drizzle dries and stabilizes the cloud layer and tends to decouple the cloud from the subcloud layer. Consequently, solid stratocumulus clouds may break up and the cloud fraction may decrease because of drizzle.
The indirect radiative forcing of aerosols is poorly constrained by the observational data underlying the simple cloud parameterizations in GCMs. signal of cloud response to increased aerosol concentrations from meteorological noise. Recent satellite observations indicate a significant decrease of cloud water in ship tracks, in contrast to an ensemble of in situ measurements showing no average change in cloud water relative to the surrounding clouds. Both results contradict the expectation of cloud water increasing in polluted clouds. We find through large-eddy simulations of stratocumulus that the trend in the satellite data is likely an artifact of sampling only overcast clouds. The simulations instead show cloud cover increasing with droplet concentrations. The simulations also show that increases in cloud water from suppressing drizzle by increased droplet concentrations are favored at night or at extremely low droplet concentrations. At typical droplet concentrations we find that the Twomey effect on cloud albedo is amplified very little by the secondary indirect effect of drizzle suppression, largely because the absorption of solar radiation by cloud water reduces boundary-layer mixing in the daytime and thereby restricts any possible increase in cloud water from drizzle suppression. The cloud and boundary layer respond to radiative heating variations on a time scale of hours, and on longer time scales respond to imbalances between large-scale horizontal advection and the entrainment of inversion air. We analyze the co-varying response of cloud water, cloud thickness, width of droplet size distributions, and dispersion of the optical depth, as well as the overall response of cloud albedo, to changes in droplet concentrations. We also dissect the underlying physical mechanisms through sensitivity studies. Ship tracks represent an ideal natural laboratory to extricate the
A cloud rift is characterized as a large-scale, persistent area of broken, low reflectivity stratocumulus clouds usually surrounded by a solid deck of stratocumulus. A rift observed off the coast of Monterey Bay, California on 16 July 1999 was studied to compare the aerosol and cloud microphysical properties in the rift with those of the surrounding solid stratus deck. Variables measured from an instrumented aircraft included temperature, water vapor, and cloud liquid water. These measurements characterized the thermodynamic properties of the solid deck and rift areas. Microphysical measurements made included aerosol, cloud drop and drizzle drop concentrations and cloud condensation nuclei (CCN) concentrations. The microphysical characteristics in a solid stratus deck differ substantially from those of a broken, cellular rift where cloud droplet concentrations are a factor of 2 lower than those in the solid cloud. Further, CCN concentrations were found to be about 3 times greater in the solid cloud area compared with those in the rift and aerosol concentrations showed a similar difference as well. Although drizzle was observed near cloud top in parts of the solid stratus cloud, the largest drizzle rates were associated with the broken clouds within the rift area. In addition to marked differences in particle concentrations, evidence of a mesoscale circulation near the solid cloud rift boundary is presented. This mesoscale circulation provides a mechanism for maintaining a rift, but further study is required to understand the initiation of a rift and the conditions that may cause it to fill.
A large convection cloud chamber has been proposed for exploring aerosol–cloud–drizzle interactions under well‐controlled turbulent conditions. Recent theoretical and numerical studies suggest that a convection cloud chamber with two heated and two cooled sidewalls can significantly enhance the liquid water content and thus benefit drizzle initiation. However, a chamber with such a sidewall configuration develops stable stratification and extremely weak turbulence therein. In this study, we conduct large‐eddy simulations of a tall convection chamber with five different sidewall configurations consisting of alternating warm and cold patches. For each configuration, the total surface area of warm patches equals that of cold patches, resulting in the same expected cloud‐free supersaturation based on a flux budget model. Results show that changing the sidewall configuration, while keeping all other factors constant, can substantially enhance turbulent mixing and improve the uniformity of thermodynamic and cloud microphysical properties in the bulk region of the chamber. In addition, turbulence strength is positively correlated with liquid water content and negatively correlated with cloud droplet number concentration, consistent with theoretical predictions. Our results highlight the advantage of building a large cloud chamber using modular patches with individually controllable temperature and humidity to achieve well‐mixed conditions.
In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.
The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s Thermodynamic Cloud Phase (THERMOCLDPHASE) Value-Added Product (VAP) provides vertically resolved thermodynamic cloud phase classifications (Zhang and Levin 2024). This VAP applies the multi-sensor methodology introduced by Shupe (2007) to identify cloud phase at the pixel level as liquid, drizzle, liquid + drizzle, rain, ice, snow, or mixed-phase. In addition, the VAP determines the overall cloud layer phase—classified as liquid, mixed-phase, or ice—based on the fraction of ice-containing pixels within the entire layer.
The G-Band Radar Demonstration for Microphysics (GRDM) campaign took place at the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) from March 15 to April 30, 2024. This was a deployment of two of NASA’s Jet Propulsion Laboratory (JPL) radars and one radar from Brookhaven National Laboratory to demonstrate the utility of high-frequency millimeter-wave radars for remote sensing of stratocumulus microphysical properties. The radars were deployed on the Ellen Browning Scripps Memorial Pier alongside the AMF instruments (Figure 1). The radars include a Ka-band (35 GHz), W-band (94 GHz), and four G-band (158, 165, 174, and 240 GHz) channels. The 240 GHz and W-band channels provide complete Doppler spectra, which are useful for advanced analysis. The 158-175 GHz channels are sensitive to the water vapor profile and are useful for attenuation correction. These radars complement the high-sensitivity ARM KAZR. The goal of the deployment was to observe drizzling stratocumulus and demonstrate the capabilities of the multifrequency radar data set to constrain profiles of liquid water content and drizzle drop characteristic size. The data are still being analyzed. The methodology to derive the cloud and precipitation parameters will exploit differential attenuation and differential reflectivity between low-frequency (Ka-band) and high-frequency (G-band) channels. The method will also exploit the capability of the G-band observations to constrain the attenuation due to water vapor. These observations will quantify the capabilities and limitations of the emerging technology of G-band radars for constraint stratocumulus cloud microphysics, which are key to constraining aerosol-cloud-precipitation interactions and low-cloud climate feedback.
The Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) was an aircraft-based campaign (15 January–26 February 2018) that deployed in situ probes and remote sensors to investigate low-level clouds over the Southern Ocean (SO). A novel methodology was developed to identify cloud boundaries and classify cloud phases in single-layer, low-level marine boundary layer (MBL) clouds below 3 km using the HIAPER Cloud Radar (HCR) and in situ measurements. The cloud base and top heights derived from HCR reflectivity, Doppler velocity, and spectrum width measurements agreed well with corresponding lidar-based and in situ estimates of cloud boundaries, with mean differences below 100 m. A liquid water content–reflectivity (LWC-Z) relationship, LWC = 0.70Z0.29, was derived to retrieve the LWC and liquid water path (LWP) from HCR profiles. The cloud phase was classified using HCR measurements, temperature, and LWP, yielding 40.6% liquid, 18.3% mixed-phase, and 5.1% ice samples, along with drizzle (29.1%), rain (3.2%), and snow (3.7%) for drizzling cloud cases. The classification algorithm demonstrates good consistency with established methods. This study provides a framework for the boundary and phase detection of MBL clouds, offering insights into SO cloud microphysics and supporting future efforts in satellite retrievals and climate model evaluation.
This study investigates marine and continental stratocumulus (Sc) cloud properties obtained from an automated implementation of a multispectral photometer retrieval. Photometer methods simultaneously retrieve cloud optical depth (τ) and cloud droplet effective radius (r e ), with estimates for liquid water path (LWP) calculated on the availability of those quantities. These applied methods evaluate retrieved cloud properties for Sc identified during a recent 6 year period over the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) program sites in Oklahoma, USA (SGP) and in the Azores, Portugal (ENA). Modest agreement in key quantity retrievals is found between the routine photometer products and multisensor collocated profiling references. Cumulative breakdowns contingent on cloud thickness indicate increases in all retrieved quantities in thicker clouds, with larger discrepancies in the relative performance between the retrievals collected in the presence of drizzle. Under continental cloud conditions, the clouds of a similar thickness and r e to those sampled under marine conditions report a factor of 1.5 larger τ and LWP. An r 2 ≅0.65 is found between photometer τ retrievals and shadowband radiometer measurements, with photometer retrievals reporting a high (relative) bias. The τ intercomparisons indicate that variability between retrievals is a factor of three larger than errors reported from individual retrieval input perturbation tests. Photometer r e retrievals suggest a low r 2 (< 0.1) having a standard deviation ≅ 3 µm when compared to ARM baseline multi-sensor radar/radiometer references (accounting for offsets in the cloud droplet number concentration assumptions of the latter). However, photometer LWP calculations remain relatively unbiased in non-drizzling conditions, with errors O (50 g m −2 ) and r 2 ≅0.5 to collocated radiometer and interferometer references. Additional sensitivity tests for island influences on marine Sc properties suggest that while island-influenced winds may promote larger cloud LWP or thickness, the influence could be within retrieval method uncertainty and/or collocated instrument variability.