Search NASA⌕ Search

SEARCH · Search NASA

Results for “cloud microphysics parameterization”

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 55 records · Page 3

A unified parameterization of clouds and turbulence using CLUBB and subcolumns in the Community Atmosphere Model

Most global climate models parameterize separate cloud types using separate parameterizations. This approach has several disadvantages, including obscure interactions between parameterizations and inaccurate triggering of cumulus parameterizations. Alternatively, a unified cloud parameterization uses one equation set to represent all cloud types. Such cloud types include stratiform liquid and ice cloud, shallow convective cloud, and deep convective cloud. Vital to the success of a unified parameterization is a general interface between clouds and microphysics. One such interface involves drawing Monte Carlo samples of subgrid variability of temperature, water vapor, cloud liquid, and cloud ice, and feeding the sample points into a microphysics scheme. This study evaluates a unified cloud parameterization and a Monte Carlo microphysics interface that has been implemented in the Community Atmosphere Model (CAM) version 5.3. Model computational expense is estimated, and sensitivity to the number of subcolumns is investigated. Results describing the mean climate and tropical variability from global simulations are presented. The new model shows a degradation in precipitation skill but improvements in shortwave cloud forcing, liquid water path, long-wave cloud forcing, precipitable water, and tropical wave simulation.

54 ENVIRONMENTAL SCIENCES↗

The Shape and Density Evolution of Snow Aggregates

Abstract Snow aggregates evolve into a variety of observed shapes and densities. Despite this diversity, models and observational studies employ fractal or Euclidean geometric measures that are assumed universal for all aggregates. This work therefore seeks to improve understanding and representation of snow aggregate geometry and its evolution by characterizing distributions of both observed and Monte Carlo–generated aggregates. Two separate datasets of best-fit ellipsoid estimates derived from Multi-Angle Snowflake Camera (MASC) observations suggest the use of a bivariate beta distribution model for capturing aggregate shapes. Product moments of this model capture shape effects to within 4% of observations. This mathematical model is used along with Monte Carlo simulated aggregates to study how combinations of monomer properties affect aggregate shape evolution. Plate aggregates of any aspect ratio produce a consistent ellipsoid shape evolution whereas thin column aggregates evolve to become more spherical. Thin column aggregates yield fractal dimensions much less than the often-assumed value of 2.0. Ellipsoid densities and fractal analogs of density (lacunarity) are much more variable depending on combinations of monomer size and shape. Simple mathematical scaling relationships can explain the persistent triaxial ellipsoid shapes that appear in both observed and modeled aggregates. Overall, both simulations and observations prove aggregates are rarely oblate. Therefore, the use of the proposed bivariate ellipsoid distribution in models will allow for similar-sized aggregates to exhibit a realistic dispersion of masses and fall speeds.

54 ENVIRONMENTAL SCIENCES↗

How Snow Aggregate Ellipsoid Shape and Orientation Variability Affects Fall Speed and Self-Aggregation Rates

Snow aggregate shapes and orientations have long been known to exhibit substantial variability. Despite this observed variability, most weather and climate prediction models use fixed power-law functions that deterministically map particle size to mass and fall speed. As such, integrated quantities like precipitation and self-aggregation rates currently ignore nonlinear effects resulting from variation in shape and orientation for aggregates of the same size. This study therefore develops an analytic framework that couples an empirically based bivariate distribution of ellipsoid shapes to classical hydrodynamic theory so as to capture an appropriate dispersion of masses, projected areas, and fall speeds for an assumed size distribution. For a fixed aggregate size, shape variations produce approximately ±0.13 m s -1 standard deviation of fall speed which increases the mass flux fall speed dispersion by more than 100% over traditional microphysics models. This increased fall speed dispersion results predominantly from shape-induced mass dispersion whereas orientation and drag dispersion play a lesser role. Shape variations can increase mass- and reflectivity-weighted fall speeds by up to 60% of traditional models whereas self-aggregation rates can increase by a factor of 100 for very small slope parameters. This implies that aggregate shape variations effectively forestall the theorized onset of fall speed distribution narrowing and subsequent quenching of the aggregation process. As a result, it is likely that secondary ice formation is necessary to prevent an ever decreasing slope parameter. The mathematical theory presented in this study is used to develop simple correction factors for snow forecast and climate models.

54 ENVIRONMENTAL SCIENCES↗

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. Furthermore, the combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

54 ENVIRONMENTAL SCIENCES↗

Fragmentation in Collisions of Snow with Graupel/Hail: New Formulation from Field Observations

Abstract Secondary ice production (SIP) has been attributed to the generation of most ice particles observed in precipitating clouds with cloud tops warmer than −36°C, from various aircraft- and ground-based field observations across the globe. One of the known SIP mechanisms is fragmentation during collisions among ice particles. It has been studied with our theoretical formulation, which has been applied in microphysical schemes of atmospheric models in a few studies. These have predicted an extensive impact on cloud glaciation and radiative properties. However, there has been a lack of experimental field studies, especially involving naturally falling snowflakes, to better understand this particular mechanism of SIP. This study reports the first field measurements with modern technology for fragmentation during collisions between naturally falling snowflakes and graupel/hail particles. This was observed with an innovatively designed portable chamber that was deployed outdoors in northern Sweden. Applying the observations from this field-based study, we optimized the existing formulation for predicting numbers of fragments from collisions of snow with graupel/hail. The observations show the average numbers of fragments per collision for dendritic (3–12 mm) and nondendritic (1–3 mm) snow were about 12 and 1, respectively. This represents a boost of predicted fragment numbers relative to our original formulation published in 2017. The updated formulation for breakup in ice–ice collisions can be implemented in the microphysical schemes of atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

Advanced Two-Moment Bulk Microphysics for Global Models. Part II: Global Model Solutions and Aerosol–Cloud Interactions*

A modified microphysics scheme is implemented in the Community Atmosphere Model, version 5 (CAM5). The new scheme features prognostic precipitation. The coupling between the microphysics and the rest of the model is modified to make it more flexible. Single-column tests show the new microphysics can simulate a constrained drizzling stratocumulus case. Substepping the cloud condensation (macrophysics) within a time step improves single-column results. Simulations of mixed-phase cases are strongly sensitive to ice nucleation. The new microphysics alters process rates in both single-column and global simulations, even at low (200 km) horizontal resolution. Thus, prognostic precipitation can be important, even in low-resolution simulations where advection of precipitation is not important. Accretion dominates as liquid water path increases in agreement with cloud-resolving model simulations and estimates from observations. The new microphysics with prognostic precipitation increases the ratio of accretion over autoconversion. Additionally, the change in process rates appears to significantly reduce aerosol–cloud interactions and indirect radiative effects of anthropogenic aerosols by up to 33% (depending on substepping) to below 1 W m-2 of cooling between simulations with preindustrial (1850) and present-day (2000) aerosol emissions.

58 GEOSCIENCES↗

Improved Convective Ice Microphysics Parameterization in the NCAR CAM Model

Partitioning deep convective cloud condensates into components that sediment and detrain, known to be a challenge for global climate models, is important for cloud vertical distribution and anvil cloud formation. In this study, we address this issue by improving the convective microphysics scheme in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 (CAM5.3). The improvements include: (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. Results from the default and modified CAM5.3 models were evaluated against observations from the U.S. Department of Energy Tropical Warm Pool-International Cloud Experiment (TWP-ICE) field campaign. The default model overestimates ice amount, which is largely attributed to the underestimation of convective ice particle sedimentation. By considering cloud ice sedimentation and rimed ice particles and applying a new convective snow terminal velocity parameterization, the vertical distribution of ice amount is much improved in the midtroposphere and upper troposphere when compared to observations. The vertical distribution of ice condensate also agrees well with observational best estimates upon considering snow detrainment. Further, comparison with observed convective updrafts reveals that current bulk model fails to reproduce the observed updraft magnitude and occurrence frequency, suggesting spectral distributions be required to simulate the subgrid updraft heterogeneity.

54 ENVIRONMENTAL SCIENCES↗

Investigation of liquid cloud microphysical properties of deep convective systems: 1. Parameterization raindrop size distribution and its application for stratiform rain estimation

Abstract To investigate liquid‐phase ( T > 3°C) cloud and precipitation microphysical properties within Deep Convective Systems (DCSs), eight DCS cases sampled by the University of North Dakota Citation II research aircraft during Midlatitude Continental Convective Clouds Experiment were selected. A full spectrum of raindrop size distribution (DSD) was constructed from 120 µm to 4000 µm through a combination of two‐dimensional cloud probe (120 to 900 µm) and High Volume Precipitation Spectrometer (900 to 4000 µm) data sets. A total of 1126 five second DSDs have been used to fit to Gamma and Exponential functions within the stratiform rain (SR) regions of DCSs. The Gamma shape μ Γ and slope λ Γ parameters are then compared with those derived from surface disdrometer measurements. The similar μ Γ ‐ λ Γ relationships but different μ Γ and λ Γ value ranges from two independent platforms at different elevations may represent the real nature of DSD shape information in clouds and at the surface. To apply the exponentially fitted DSD parameters to precipitation estimation using Next Generation Weather Radar (NEXRAD) radar reflectivity factor Z e , the terms N 0 E and λ E have been parameterized as a function of Z e using an empirical N 0 E ‐ λ E relationship. The averaged SR rain rate retrieved from this study is almost identical to the surface measurements, while the NEXRAD Q2 precipitation is twice as large. The comparisons indicate that the new DSD parameterization scheme is robust, while the Q2 SR precipitation estimation based on Marshall‐Palmer Z‐R relationship, where a constant DSD intercept parameter ( N 0 E ) was assumed, needs to be improved for heavy precipitation cases.

Wang, Jingyu↗

Coupling Spectral-bin Cloud Microphysics with the MOSAIC Aerosol Model in WRF-Chem: Methodology and Results for Marine Stratocumulus Clouds

Aerosol-cloud interaction processes can be represented more physically with bin cloud microphysics relative to bulk microphysical parameterizations. However, due to computational power limitations in the past, bin cloud microphysics was often run with very simple aerosol treatments. The purpose of this study is to represent better aerosol-cloud interaction processes in the Chemistry version of Weather Research and Forecast model (WRF-Chem) at convection-permitting scales by coupling spectral-bin cloud microphysics (SBM) with the MOSAIC sectional aerosol model. A flexible interface is built that exchanges cloud and aerosol information between them. The interface contains a new bin aerosol activation approach, which replaces the treatments in the original SBM. It also includes the modified aerosol resuspension and in-cloud wet removal processes with the droplet loss tendencies and precipitation fluxes from SBM. The newly coupled system is evaluated for two marine stratocumulus cases over the Southeast Pacific Ocean with either a simplified aerosol setup or full-chemistry. We compare the aerosol activation process in the newly-coupled SBM-MOSAIC against the SBM simulation without chemistry using a simplified aerosol setup, and the results show consistent activation rates. A longer time simulation reinforces that aerosol resuspension through cloud drop evaporation plays an important role in replenishing aerosols and impacts cloud and precipitation in marine stratocumulus clouds. Evaluation of the coupled SBM-MOSAIC with full-chemistry using aircraft measurements suggests that the new model works realistically for the marine stratocumulus clouds, and improves the simulation of cloud microphysical properties compared to a simulation using MOSAIC coupled with the Morrison two-moment microphysics.

Gao, Wenhua↗

Final Technical Report for "Reducing tropical precipitation biases in CESM"

In state-of-the-art climate models, each cloud type is treated using its own separate cloud parameterization and its own separate microphysics parameterization. This use of separate schemes for separate cloud regimes is undesirable because it is theoretically unfounded, it hampers interpretation of results, and it leads to the temptation to overtune parameters. In this grant, we have created a climate model that contains a unified cloud parameterization (“CLUBB”) and a unified microphysics parameterization (“MG2”). In this model, all cloud types --- including marine stratocumulus, shallow cumulus, and deep cumulus --- are represented with a single equation set. This model improves the representation of convection in the Tropics. The model has been compared with ARM observations. The chief benefit of the project is to provide a climate model that is based on a more theoretically rigorous formulation.

54 ENVIRONMENTAL SCIENCES↗

Cold-Season Precipitation Sensitivity to Microphysical Parameterizations: Hydrologic Evaluations Leveraging Snow Lidar Datasets

Abstract Cloud microphysical processes are an important facet of atmospheric modeling, as they can control the initiation and rates of snowfall. Thus, parameterizations of these processes have important implications for modeling seasonal snow accumulation. We conduct experiments with the Weather Research and Forecasting (WRF V4.3.3) Model using three different microphysics parameterizations, including a sophisticated new scheme (ISHMAEL). Simulations are conducted for two cold seasons (2018 and 2019) centered on the Colorado Rockies’ ∼750-km 2 East River watershed. Precipitation efficiencies are quantified using a drying-ratio mass budget approach and point evaluations are performed against three NRCS SNOTEL stations. Precipitation and meteorological outputs from each are used to force a land surface model (Noah-MP) so that peak snow accumulation can be compared against airborne snow lidar products. We find that microphysical parameterization choice alone has a modest impact on total precipitation on the order of ±3% watershed-wide, and as high as 15% for certain regions, similar to other studies comparing the same parameterizations. Precipitation biases evaluated against SNOTEL are 15% ± 13%. WRF Noah-MP configurations produced snow water equivalents with good correlations with airborne lidar products at a 1-km spatial resolution: Pearson’s r values of 0.9, RMSEs between 8 and 17 cm, and percent biases of 3%–15%. Noah-MP with precipitation from the PRISM geostatistical precipitation product leads to a peak SWE underestimation of 32% in both years examined, and a weaker spatial correlation than the WRF configurations. We fall short of identifying a clearly superior microphysical parameterization but conclude that snow lidar is a valuable nontraditional indicator of model performance.

54 ENVIRONMENTAL SCIENCES↗

Quantifying and attributing time step sensitivities in present-day climate simulations conducted with EAMv1

Abstract. This study assesses the relative importance of time integration error in present-day climate simulations conducted with the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) at 1∘ horizontal resolution. We show that a factor-of-6 reduction of time step size in all major parts of the model leads to significant changes in the long-term mean climate. Examples of changes in 10-year mean zonal averages include the following: up to 0.5 K of warming in the lower troposphere and cooling in the tropical and subtropical upper troposphere, 1 %–10 % decreases in relative humidity throughout the troposphere, and 10 %–20 % decreases in cloud fraction in the upper troposphere and decreases exceeding 20 % in the subtropical lower troposphere. In terms of the 10-year mean geographical distribution, systematic decreases of 20 %–50 % are seen in total cloud cover and cloud radiative effects in the subtropics. These changes imply that the reduction of temporal truncation errors leads to a notable although unsurprising degradation of agreement between the simulated and observed present-day climate; to regain optimal climate fidelity in the absence of those truncation errors, the model would require retuning. A coarse-grained attribution of the time step sensitivities is carried out by shortening time steps used in various components of EAM or by revising the numerical coupling between some processes. Our analysis leads to the finding that the marked decreases in the subtropical low-cloud fraction and total cloud radiative effect are caused not by the step size used for the collectively subcycled turbulence, shallow convection, and stratiform cloud macrophysics and microphysics parameterizations but rather by the step sizes used outside those subcycles. Further analysis suggests that the coupling frequency between the subcycles and the rest of EAM significantly affects the subtropical marine stratocumulus decks, while deep convection has significant impacts on trade cumulus. The step size of the cloud macrophysics and microphysics subcycle itself appears to have a primary impact on cloud fraction in the upper troposphere and also in the midlatitude near-surface layers. Impacts of step sizes used by the dynamical core and the radiation parameterization appear to be relatively small. These results provide useful clues for future studies aiming at understanding and addressing the root causes of sensitivities to time step sizes and process coupling frequencies in EAM. While this study focuses on EAMv1 and the conclusions are likely model-specific, the presented experimentation strategy has general value for weather and climate model development, as the methodology can help researchers identify and understand sources of time integration error in sophisticated multi-component models.

54 ENVIRONMENTAL SCIENCES↗

Collaborative Research: Quantifying the Uncertainties of Aerosol Indirect Effects and Impacts on Decadal-Scale Climate Variability in NCAR CAM5 and CESM1

The goal of this proposed project is to assess the climatic importance and sensitivity of aerosol indirect effect (AIE) to cloud and aerosol processes and feedbacks, which include organic aerosol hygroscopicity, cloud condensation nuclei (CCN) activation kinetics, Giant CCN, cloud-scale entrainment, ice nucleation in mixed-phase and cirrus clouds, and treatment of subgrid variability of vertical velocity. A key objective was to link aerosol, cloud microphysics and dynamics feedbacks in CAM5 with a suite of internally consistent and integrated parameterizations that provide the appropriate degrees of freedom to capture the various aspects of the aerosol indirect effect. The proposal integrated new parameterization elements into the cloud microphysics, moist turbulence and aerosol modules used by the NCAR Community Atmospheric Model version 5 (CAM5). The CAM5 model was then used to systematically quantify the uncertainties of aerosol indirect effects through a series of sensitivity tests with present-day and preindustrial aerosol emissions. New parameterization elements were developed as a result of these efforts, and new diagnostic tools & methodologies were also developed to quantify the impacts of aerosols on clouds and climate within fully coupled models. Observations were used to constrain key uncertainties in the aerosol-cloud links. Advanced sensitivity tools were developed and implements to probe the drivers of cloud microphysical variability with unprecedented temporal and spatial scale. All these results have been published in top and high impact journals (or are in the final stages of publication). This proposal has also supported a number of outstanding graduate students.

54 ENVIRONMENTAL SCIENCES↗

Collaborative Research: Reducing tropical precipitation biases in CESM — Tests of unified parameterizations with ARM observations

In state-of-the-art climate models, each cloud type is treated using its own separate cloud parameterization and its own separate microphysics parameterization. This use of separate schemes for separate cloud regimes is undesirable because it is theoretically unfounded, it hampers interpretation of results, and it leads to the temptation to overtune parameters. In this grant, we are creating a climate model that contains a unified cloud parameterization and a unified microphysics parameterization. This model will be used to address the problems of excessive frequency of drizzle in climate models and excessively early onset of deep convection in the Tropics over land. The resulting model will be compared with ARM observations.

54 ENVIRONMENTAL SCIENCES↗

Understanding Changes in Cloud Simulations from E3SM Version 1 to Version 2

This study documents clouds simulated by the Energy Exascale Earth System Model (E3SM) version 2 (E3SMv2) and attempts to understand what causes the model behavior change in clouds relative to E3SMv1. This is done by analyzing the last 30-year (1985–2014) data from the 165-year historical simulations using E3SMv1 and v2 and four sensitivity tests to isolate the impact of changes in model parameter choices in its turbulence, shallow convection, and cloud macrophysics parameterization (CLUBB), microphysical parameterization (MG2), and deep convection scheme (ZM), as well as model physics changes in convective triggering. It is shown that E3SMv2 significantly improves the simulation of subtropical coastal stratocumulus clouds (Sc) and clouds with optical depth larger than 3.6 over the stratocumulus to cumulus transition regimes, where Shortwave Cloud Radiative Effect (SWCRE) is also improved, and the Southern Ocean (SO) while seeing an overall slight degradation in low clouds over other tropical and subtropical oceans. The better performance in E3SMv1 over those regions is partially due to error compensation between its simulated optically thin and intermediate low clouds for which E3SMv2 actually improves simulation of optically intermediate low clouds. Sensitivity tests indicate that the changes in low clouds are primarily due to the tuning made in CLUBB. The impact of the ZM tuning is mainly on optically intermediate and thick high clouds, contributing to improved SWCRE and Longwave Cloud Radiative Effect (LWCRE). The impact of the MG2 tuning and the new convective trigger is primarily on the high latitudes and the Southern Ocean (SO). They have a relatively smaller impact on clouds than do the CLUBB and ZM tunings. This study offers additional insights about clouds simulated in E3SMv2 by utilizing multiple data sets and the COSP diagnostic tool as well as through sensitivity tests. The improved understanding will benefit the future E3SM developments.

54 ENVIRONMENTAL SCIENCES↗

Understanding changes in cloud simulations from E3SM version 1 to version 2

Abstract. This study documents clouds simulated by the Energy Exascale Earth System Model (E3SM) version 2 (E3SMv2) and attempts to understand what causes the model behavior change in clouds relative to E3SMv1. This is done by analyzing the last 30-year (1985–2014) data from the 165-year historical simulations using E3SMv1 and v2 and four sensitivity tests to isolate the impact of changes in model parameter choices in its turbulence, shallow convection, and cloud macrophysics parameterization (Cloud Layers Unified By Binormals, CLUBB); microphysical parameterization (MG2); and deep-convection scheme (ZM), as well as model physics changes in convective triggering. It is shown that E3SMv2 significantly improves the simulation of subtropical coastal stratocumulus clouds and clouds with optical depth larger than 3.6 over the stratocumulus-to-cumulus transition regimes, where the shortwave cloud radiative effect (SWCRE) is also improved, and the Southern Ocean (SO) while seeing an overall slight degradation in low clouds over other tropical and subtropical oceans. The better performance in E3SMv1 over those regions is partially due to error compensation between its simulated optically thin and intermediate low clouds for which E3SMv2 actually improves simulation of optically intermediate low clouds. Sensitivity tests indicate that the changes in low clouds are primarily due to the tuning done in CLUBB. The impact of the ZM tuning is mainly on optically intermediate and thick high clouds, contributing to an improved SWCRE and longwave cloud radiative effect (LWCRE). The impact of the MG2 tuning and the new convective trigger is primarily on the high latitudes and the SO. They have a relatively smaller impact on clouds than CLUBB tuning and ZM tuning do. This study offers additional insights into clouds simulated in E3SMv2 by utilizing multiple data sets and the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) diagnostic tool as well as sensitivity tests. The improved understanding will benefit future E3SM developments.

58 GEOSCIENCES↗