Impacts of Aerosol Scavenging and Processing on the Transition of a Stratocumulus Cloud System to Open Cells: A Comparison of Lagrangian and Bin Microphysics Schemes in LES
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.
Abstract. Extensive regions of marine boundary layer cloud impact the radiative balance through their significant shortwave albedo while having little impact on outgoing longwave radiation. Despite this importance, these cloud systems remain poorly represented in large-scale models due to difficulty in representing the processes that drive their life cycle and coverage. In particular, the mesoscale organization and cellular structure of marine boundary clouds have important implications for the subsequent cloud feedbacks. In this study, we use long-term (2013–2018) observations from the Atmospheric Radiation Measurement (ARM) Facility's Eastern North Atlantic (ENA) site on Graciosa Island, Azores, Portugal, to identify cloud cases with open- or closed-cellular organization. More than 500 h of each organization type are identified. The ARM observations are combined with reanalysis and satellite products to quantify the cloud, precipitation, aerosol, thermodynamic, and large-scale synoptic characteristics associated with these cloud types. Our analysis shows that both cloud organization populations occur during similar sea surface temperature conditions, but the open-cell cases are distinguished by stronger cold-air advection and large-scale subsidence compared to the closed-cell cases, consistent with their formation during cold-air outbreaks. We also find that the open-cell cases were associated with deeper boundary layers, stronger low-level winds, and higher rain rates compared to their closed-cell counterparts. Finally, raindrops with diameters larger than 1 mm were routinely recorded at the surface during both populations, with a higher number of large drops during the open-cellular cases. The similarities and differences noted herein provide important insights into the environmental and cloud characteristics during varying marine boundary layer cloud mesoscale organization and will be useful for the evaluation of model simulations for ENA marine clouds.
Abstract. A large part of the uncertainty in climate projections comes from uncertain aerosol properties and aerosol–cloud interactions as well as the difficulty in remotely sensing them. The southeastern Atlantic functions as a natural laboratory to study biomass-burning smoke and to constrain this uncertainty. We address these gaps by comparing the Weather Research and Forecasting with Chemistry Community Atmosphere Model (WRF-CAM5) to the multi-campaign observations ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS), CLARIFY (CLoud–Aerosol–Radiation Interaction and Forcing), and LASIC (Layered Atlantic Smoke Interactions with Clouds) in the southeastern Atlantic in August 2017 to evaluate a large range of the model's aerosol chemical properties, size distributions, processes, and transport, as well as aerosol–cloud interactions. Overall, while WRF-CAM5 is able to represent smoke properties and transport, some key discrepancies highlight the need for further analysis. Observations of smoke composition show an overall decrease in aerosol mean diameter as smoke ages over 4–12 d, while the model lacks this trend. A decrease in the mass ratio of organic aerosol (OA) to black carbon (BC), OA:BC, and the OA mass to carbon monoxide (CO) mixing ratio, OA:CO, suggests that the model is missing processes that selectively remove OA from the particle phase, such as photolysis and heterogeneous aerosol chemistry. A large (factor of ∼2.5) enhancement in sulfate from the free troposphere (FT) to the boundary layer (BL) in observations is not present in the model, pointing to the importance of properly representing secondary sulfate aerosol formation from marine dimethyl sulfide and gaseous SO2 smoke emissions. The model shows a persistent overprediction of aerosols in the marine boundary layer (MBL), especially for clean conditions, which multiple pieces of evidence link to weaker aerosol removal in the modeled MBL than reality. This evidence includes several model features, such as not representing observed shifts towards smaller aerosol diameters, inaccurate concentration ratios of carbon monoxide and black carbon, underprediction of heavy rain events, and little evidence of persistent biases in modeled entrainment. The average below-cloud aerosol activation fraction (NCLD/NAER) remains relatively constant in WRF-CAM5 between field campaigns (∼0.65), while it decreases substantially in observations from ORACLES (∼0.78) to CLARIFY (∼0.5), which could be due to the model misrepresentation of clean aerosol conditions. WRF-CAM5 also overshoots an observed upper limit on liquid cloud droplet concentration around NCLD= 400–500 cm−3 and overpredicts the spread in NCLD. This could be related to the model often drastically overestimating the strength of boundary layer vertical turbulence by up to a factor of 10. We expect these results to motivate similar evaluations of other modeling systems and promote model development to reduce critical uncertainties in climate simulations.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Equilibrium solutions are constructed using an idealized model for mixed and partially mixed boundary layers to show the dependence of the equilibrium cloud layer thickness on different physical processes. The condition when a mixed layer becomes cloud-free is reduced to a comparison between three pressure scales, associated with the surface wind and divergence, the radiative cooling and the subsaturation of the entrained air. The evaporation of layer cloud by cloud-top mixing is also formulated in terms of simple pressure scales.
Examination of cloud radiance fields derived from satellite observations sometimes indicates the existence of a range of scales over which the statistics of the field are scale invariant. Many methods were developed to quantify this scaling behavior in geophysics. The usefulness of such techniques depends both on the physics of the process being robust over a wide range of scales and on the availability of high resolution, low noise observations over these scales. These techniques (area perimeter relation, distribution of areas, estimation of the capacity, d0, through box counting, correlation exponent) are applied to the high resolution satellite data taken during the FIRE experiment and provides initial estimates of the quality of data required by analyzing simple sets. The results of the observed fields are contrasted with those of images of objects with known characteristics (e.g., dimension) where the details of the constructed image simulate current observational limits. Throughout when cloud elements and cloud boundaries are mentioned; it should be clearly understood that by this structures in the radiance field are meant: all the boundaries considered are defined by simple threshold arguments.
A bulk planetary boundary layer (PBL) model was developed with a simple internal vertical structure and a simple second-order closure, designed for use as a PBL parameterization in a large-scale model. The model allows the mean fields to vary with height within the PBL, and so must address the vertical profiles of the turbulent fluxes, going beyond the usual mixed-layer assumption that the fluxes of conservative variables are linear with height. This is accomplished using the same convective mass flux approach that has also been used in cumulus parameterizations. The purpose is to show that such a mass flux model can include, in a single framework, the compensating subsidence concept, downgradient mixing, and well-mixed layers.
Explore the source record for details and available documents.
Two tasks were completed during this period. In the first, we examined the polarization of millimeter-wavelength radar beams scattered by ice crystals. Because of their non-spherical shape and size, ice crystals depolarize the incident polarized radar beam. In principle, this depolarization can be used to identify ice from liquid water, as well as provide some information on size. However, the amount of de-polarization is small, producing only a weak signal at the receiver. Our task was to determine the magnitude of such a signal and decide if our radar would be capable of measuring it under typical cirrus conditions. The theoretical study was carried out by Henrietta Lemke, a visiting graduate student from Germany. She had prior experience using a discrete dipole code to compute scattering depolarization. Dr. Kultegin Aydin of the Penn State Electrical Engineering Department, who is also expert in this area, consulted with us on this project at no cost to the project. Our conclusion was that the depolarization signal is too weak to be usefully measured by our system. Therefore we proceeded no further in this study. The second task involved the study of the effect of stratus microphysics on surface cloud forcing. Manajit Sengupta, a graduate student, and the project PI jointly carried out this task. The study used data culled from over a year of continuous radar and radiometer observations at the Atmospheric Radiation Measurement (ARM) site in Oklahoma. The study compared solar radiation calculations made using constant microphysics with calculations made using a retrieved mean particle size. The results showed that on average the constant microphysics produced the correct forcing when compared with the observed forcing. We conclude, therefore, that there is little impetus on radiation grounds alone to include explicit microphysics in climate models. The question of pollutant particle emission impacts on microphysics remains to be resolved. A manuscript is in preparation and will be submitted this year.
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.
The 2004 Multi-Phase Arctic Cloud Experiment (M-PACE) IOP at the ARM NSA site focused on measuring the properties of autumn transition-season arctic stratus and the environmental conditions controlling them, including concentrations of heterogeneous ice nuclei. Our work aims to use a large-eddy simulation (LES) code with embedded size-resolved aerosol and cloud microphysics to identify factors controlling multi-phase arctic stratus. Our preliminary simulations of autumn transition-season clouds observed during the 1994 Beaufort and Arctic Seas Experiment (BASE) indicated that low concentrations of ice nuclei, which were not measured, may have significantly lowered liquid water content and thereby stabilized cloud evolution. However, cloud drop concentrations appeared to be virtually immune to changes in liquid water content, indicating an active Bergeron process with little effect of collection on drop number concentration. We will compare these results with preliminary simulations from October 8-13 during MPACE. The sensitivity of cloud properties to uncertainty in other factors, such as large-scale forcings and aerosol profiles, will also be investigated. Based on the LES simulations with M-PACE data, preliminary results from the NASA GlSS single-column model (SCM) will be used to examine the sensitivity of predicted cloud properties to changing cloud drop number concentrations for multi-phase arctic clouds. Present parametrizations assumed fixed cloud droplet number concentrations and these will be modified using M-PACE data.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
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.
Stratocumulus clouds, with their low cloud base and top, affect the atmospheric boundary layer wind and turbulence profile, thereby modulating wind energy resources. GOES satellite data reveal an abundance of stratocumulus clouds in the late spring and summer months off the coast of northern and central California, where there are active plans to deploy floating offshore wind farms at two lease areas (near Morro Bay and Humboldt). Since the fall of 2020, two buoys equipped with multiple instruments, including Doppler lidar, have been deployed for about 1 year in these wind farm lease areas to assess the rotor layer wind conditions in these locations. The objective of this study is to evaluate how well the High-Resolution Rapid Refresh (HRRR) model represents stratocumulus cloud characteristics and turbine-relevant rotor layer winds (surface to 300 m) by comparing HRRR simulations with buoy and satellite observations. We first find that the HRRR model reproduces the seasonal cycle of cloud top height reasonably well in these regions. However, during the warm season – especially at Morro Bay – the HRRR-simulated stratocumulus clouds tend to have lower tops by about 150 m and exhibit weaker diurnal cycles than satellite observations. Our analysis also shows that rotor layer wind speeds and vertical shear are stronger at Humboldt than at Morro Bay, and both are generally stronger under clear-sky conditions. Finally, the HRRR model bias in rotor layer wind speed is small under cloudy conditions but larger and dependent on observed wind speed under clear skies. Specifically, HRRR underestimates wind speeds at Morro Bay and overestimates them at Humboldt under clear-sky conditions.
Cloud liquid water path (L) adjusts to perturbations in cloud droplet number concentration (N) over time. We explore the magnitude and timescale of this adjustment in nocturnal non-precipitating marine stratocumuli using large eddy simulations of baseline conditions and aerosol seeding experiments for 22 meteorological conditions. The results confirm that the L adjustment (δL) slope (k) is more negative for simulation pairs with relatively low N and less negative for high N. Overall, k is unlikely to be lower than −0.4 within 24 h of seeding start, meaning the L adjustment is unlikely to fully offset the brightening due to the Twomey effect. After seeding, the δL becomes increasingly negative, which can be characterized by an exponential convergence. This evolution is governed by a short timescale around 5 h and lasts for around 8–12 h. It is driven by the feedback between entrainment, L, and boundary layer (BL) turbulence. Other processes, including radiation, surface fluxes, and subsidence, respond to the seeding weakly. This short timescale is insensitive to the amount of seeding, making the evolution of δL and some other deviations similar for different seeding amounts after appropriate scaling. The timescale of k evolution is closely related to the δL timescale and hence also short, while it could also be affected by the δN evolution. The results are most relevant to conditions where seeding is applied to a large area of marine stratocumulus in well-mixed and overcast BL where shear is not a primary source of turbulence.