Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cloud parameterizations”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Demonstrating Hierarchical System Development With the Common Community Physics Package Single‐Column Model: A Case Study Over the Southern Great Plains

This study demonstrates a specific application of the hierarchical system development (HSD) approach to investigate, analyze, and attribute model issues within the Unified Forecast System (UFS), with a focus on process isolation. By evaluating a non‐precipitating, shallow cumulus case at the Atmospheric Radiation Measurement Southern Great Plains site in the UFS global forecast against the observation, the investigation identifies a warmer and deeper daytime convective planetary boundary layer (PBL) and misrepresented nocturnal PBL transition. Hypothesis testing, which employs the Common Community Physics Package (CCPP) single‐column model (SCM) and uses the same physics as the UFS global model, confirms that these issues are attributed to the model physics and initialization. Specifically, misrepresented PBL processes are linked to problematic surface condition and a lack of cloud formation, which may stem from deficiencies in PBL and cloud microphysics parameterizations and their interactions. The UFS initial condition contributes to an earlier, excessively collapsed daytime convective boundary layer and a lack of decoupling between the stable boundary layer and residual layer late in the afternoon. This work introduces an avenue for the community to engage with the application of HSD, along with the CCPP and CCPP SCM, to understand the interplay of model physics, disentangle the roles of model components, as well as facilitate model and forecast improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulated Precipitation Diurnal Variation With a Deep Convective Closure Subject to Shallow Convection in Community Atmosphere Model Version 5 Coupled With CLUBB

In order to improve the physical consistency between shallow and deep convection, we modify the deep convective closure in the Community Atmosphere Model version 5 (CAM5) coupled with a third-order turbulence closure parameterization (i.e., Cloud Layers Unified by Binormals; CLUBB). The revised closure reserves a portion of the total convective available potential energy for shallow convection via utilizing the heating and moistening profiles from CLUBB to distribute moisture and energy between shallow and deep convection. Simulations at two resolutions (i.e. 2° and 0.5°, respectively) are conducted to investigate the impacts of convective closure on the simulated precipitation diurnal variations. Results from low-resolution simulations show that the revised closure suppresses deep convection until the lower troposphere is sufficiently moistened by shallow convection, which improves the precipitation diurnal variation simulations compared to the default closure, with the precipitation diurnal peak over tropical lands delayed from 12LST to 19LST. The revised closure better simulates the diurnal variations for precipitation over the Asian monsoon region, such as the delayed precipitation onset, but still fails to well capture the nighttime peak for precipitation there. This deficiency is alleviated to some extent when applying the revised closure in highresolution simulations, but nighttime precipitation is still underestimated probably because key processes responsible for nighttime convection are missing. Altogether, our results indicate that establishing the consistency between shallow and deep convection is critical for the precipitation diurnal cycle simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Model Intercomparison Study of Aerosol‐Cloud‐Turbulence Interactions in a Cloud Chamber: 1. Model Results

This study presents the first model intercomparison of aerosol‐cloud‐turbulence interactions in a controlled cloudy Rayleigh‐Bénard Convection chamber environment, utilizing the Pi Chamber at Michigan Technological University. We analyzed simulated cloud chamber‐averaged statistics of microphysics and thermodynamics in a warm‐phase, cloudy environment under steady‐state conditions at varying aerosol injection rates. Simulation results from seven distinct models (DNS, LES, and a 1D turbulence model) were compared. Our findings demonstrate that while all models qualitatively capture observed trends in droplet number concentration, mean radius, and droplet size distributions at both high and low aerosol injection rates, significant quantitative differences were observed. Notably, droplet number concentrations varied by over two orders of magnitude between models for the same injection rates, indicating sensitivities to the model treatments in droplet activation and removal and wall fluxes. Furthermore, inconsistencies in vertical relative humidity profiles and in achieving steady‐state liquid water content suggest the need for further investigation into the mechanisms driving these variations. Despite these discrepancies, the models generally reproduced consistent power‐law relationships between the microphysical variables. This model intercomparison underscores the importance of controlled cloud chamber experiments for validating and improving cloud microphysical parameterizations. Recommendations for future modeling studies are also highlighted, including constraining wall conditions and processes, investigating droplet/aerosol removal (including sidewall losses), and conducting simplified experiments to isolate specific processes contributing to model divergence and reduce model uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Recent Advancements in the Numerical Simulation of Surface Irradiance for Solar Energy Applications

This paper briefly reviews the National Renewable Energy Laboratory's recent efforts to develop all-sky solar irradiance models for solar energy applications. The Fast All-sky Radiation Model for Solar applications (FARMS) uses the simulation of clear-sky transmittance and reflectance and a parameterization of cloud transmittance and reflectance to rapidly compute broadband irradiances on horizontal surfaces. The accuracy of FARMS is comparable to that of two-stream approximation, but it is approximately 1,000 times faster. A FARMS for Narrowband Irradiance over Tilted surfaces (FARMS-NIT) has been developed to compute spectral irradiances on photovoltaic (PV) panels in 2,002 wavelength bands. FARMS-NIT has been extended to bifacial PV panels by accounting for solar radiation reaching the backside of a PV panel.

14 SOLAR ENERGY↗
DOE OSTI2019

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↗
DOE OSTI2021

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↗
DOE OSTI2021

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↗
DOE OSTI2024

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↗
DOE OSTI2024

Perturbing Parameters to Understand Cloud Contributions to Climate Change

Abstract The sensitivity of cloud feedbacks to atmospheric model parameters is evaluated using a CAM6 perturbed parameter ensemble (PPE). The CAM6 PPE perturbs 45 parameters across 262 simulations, 206 of which are used here. The spread in the total cloud feedback and its six components across the CAM6 PPE are comparable to the spread across the CMIP6 and AMIP ensembles, indicating that parametric uncertainty mirrors structural uncertainty. However, the high-cloud altitude feedback is generally larger in the CAM6 PPE than WCRP assessment, CMIP6, and AMIP values. We evaluate the influence of each of the 45 parameters on the total cloud feedback and each of the six cloud feedback components. We also explore whether the CAM6 PPE can be used to constrain the total cloud feedback, with inconclusive results. Further, we find that despite the large parametric sensitivity of cloud feedbacks in CAM6, a substantial increase in cloud feedbacks from CAM5 to CAM6 is not a result of changes in parameter values. Notably, the CAM6 PPE is run with a more recent version of CAM6 (CAM6.3) than was used for AMIP (CAM6.0) and has a smaller total cloud feedback (0.56 W m −2 K −1 ) as compared to CAM6.0 (0.81 W m −2 K −1 ) owing primarily to reductions in low clouds over the tropics and midlatitudes. The work highlights the large sensitivity of cloud feedbacks to both parameter values and structural details in CAM6.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Evaluation of autoconversion and accretion enhancement factors in general circulation model warm-rain parameterizations using ground-based measurements over the Azores

A great challenge in climate modeling is how to parameterize subgrid cloud processes, such as autoconversion and accretion in warm-rain formation. In this study, we use ground-based observations and retrievals over the Azores to investigate the so-called enhancement factors, E auto and E accr , which are often used in climate models to account for the influence of subgrid variance of cloud and precipitation water on the autoconversion and accretion processes. E auto and E accr are computed for different equivalent model grid sizes. The calculated E auto values increase from 1.96 (30 km) to 3.2 (180 km), and the calculated E accr values increase from 1.53 (30 km) to 1.76 (180 km). Comparing the prescribed enhancement factors in Morrison and Gettleman (2008, MG08) to the observed ones, we found that a higher E auto (3.2) at small grids and lower E accr (1.07) are used in MG08, which might explain why most of the general circulation models (GCMs) produce too-frequent precipitation events but with too-light precipitation intensity. The ratios of the rain to cloud water mixing ratio ( q r / q c ) at E accr =1.07 and E accr =2.0 are 0.063 and 0.142, respectively, from observations, further suggesting that the prescribed value of E accr =1.07 used in MG08 is too small to simulate precipitation intensity correctly. Both E auto and E accr increase when the boundary layer becomes less stable, and the values are larger in precipitating clouds ( CLWP>75 gm -2 ) than those in non-precipitating clouds ( CLWP<75 gm -2 ). Therefore, the selection of E auto and E accr values in GCMs should be regime- and resolution-dependent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package, version 2: assessing aerosols, clouds, and aerosol–cloud interactions via field campaign and long-term observations

Abstract. Poor representations of aerosols, clouds, and aerosol–cloud interactions (ACIs) in Earth system models (ESMs) have long been the largest uncertainties in predicting global climate change. Huge efforts have been made to improve the representation of these processes in ESMs, and the key to these efforts is the evaluation of ESM simulations with observations. Most well-established ESM diagnostics packages focus on the climatological features; however, they lack process-level understanding and representations of aerosols, clouds, and ACIs. In this study, we developed the Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package to facilitate the routine evaluation of aerosols, clouds, and ACIs simulated the Energy Exascale Earth System Model (E3SM) from the US Department of Energy (DOE). This paper documents its version 2 functionality (ESMAC Diags v2), which has substantial updates compared with version 1 (Tang et al., 2022a). The simulated aerosol and cloud properties have been extensively compared with in situ and remote-sensing measurements from aircraft, ship, surface, and satellite platforms in ESMAC Diags v2. It currently includes six field campaigns and two permanent sites covering four geographical regions: the eastern North Atlantic, the central US, the northeastern Pacific, and the Southern Ocean. These regions produce frequent liquid- or mixed-phase clouds, with extensive measurements available from the DOE Atmospheric Radiation Measurement user facility and other agencies. ESMAC Diags v2 generates various types of single-variable and multivariable diagnostics, including percentiles, histograms, joint histograms, and heatmaps, to evaluate the model representation of aerosols, clouds, and ACIs. Select examples highlighting the capabilities of ESMAC Diags are shown using E3SM version 2 (E3SMv2). In general, E3SMv2 can reasonably reproduce many observed aerosol and cloud properties, with biases in some variables such as aerosol particle and cloud droplet sizes and number concentrations. The coupling of aerosol and cloud number concentrations may be too strong in E3SMv2, possibly indicating a bias in processes that control aerosol activation. Furthermore, the liquid water path response to a perturbed cloud droplet number concentration behaves differently in E3SMv2 and observations, which warrants further study to improve the cloud microphysics parameterizations in E3SMv2.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

The GISS global climate-middle atmosphere model. I - Model structure and climatology

A coupled global climate/middle atmosphere model (GCMAM) is developed by extending the Hansen et al. (1983) GISS global climate model to include the middle atmosphere up to 85 km. The model includes numerical solutions of the primitive equations, calculations of the radiative and surface fluxes, and a complete hydrologic cycle with convective and cloud-cover parameterizations. In addition, a parameterized gravity wave drag cycle is incorporated, in which gravity-wave momentum fluxes due to flow over topography, wind shear, and convection are calculated at each grid box, using theoretical relationships between the grid-scale variables and expected source strengths. The results of the GCMAM demonstrate that the model produces a reasonable simulation of the stratosphere. It is shown that the improvements over the previous version were achieved largely through the incorporation of the parameterized gravity wave drag.

Rind, D.↗
NASA NTRS1988

Use of the Terminal Area Simulation System (TASS) to study microburst wind shears

Ground-based and airborne Doppler radar and LIDAR systems are being designed to alert pilots when a hazardous windshear is present. A key element in this design effort is understanding the microburst itself. This is accomplished by means of the TASS model which was developed for NASA by Proctor (1987a,b). The time-dependent TASS model has two vesions: a two-demensional high resolution axisymmetric model, and a three-demensional model. The model includes a sophisticated parameterization of cloud microphysics and a friction layer, both of which are essential to a realistic simulation of the microburst phenomenon. The TASS model has been successfully tested on well-observed convective events.

Byrd, Gregory P.↗
NASA NTRS1992

Lidar Development for the Atmospheric Radiation Measurement (ARM) Program

The Atmospheric Radiation Measurement (ARM) Program was initiated by the Department of Energy's Atmospheric and Climate Research Division (DOE/ACRD) to improve parameterizations of cloud radiative effects in general circulation models of the Earth's climate. The Instrument Development Program (IDP), a major component of ARM, is directed toward the improvement of radiometric and remote sensing instrumentation pertinent to near real-time characterizations of atmospheric parameters. The ultimate goal of this program is to provide state-of-the-art instruments suitable for installation and extended use in ARM's field measurement component, the Cloud and Radiation Testbed (CART).

Griffin, J.↗
NASA NTRS2012

Comparing Physics Scheme Performance for a Lake Effect Snowfall Event in Northern Lower Michigan

High resolution forecast models, such as those used to predict severe convective storms, can also be applied to predictions of lake effect snowfall. A high resolution WRF model forecast model is provided to support operations at NWS WFO Gaylord, Michigan, using a 12 ]km and 4 ]km nested configuration. This is comparable to the simulations performed by other NWS WFOs adjacent to the Great Lakes, including offices in the NWS Eastern Region who participate in regional ensemble efforts. Ensemble efforts require diversity in initial conditions and physics configurations to emulate the plausible range of events in order to ascertain the likelihood of different forecast scenarios. In addition to providing probabilistic guidance, individual members can be evaluated to determine whether they appear to be biased in some way, or to better understand how certain physics configurations may impact the resulting forecast. On January 20 ]21, 2011, a lake effect snow event occurred in Northern Lower Michigan, with cooperative observing and CoCoRaHS stations reporting new snow accumulations between 2 and 8 inches and liquid equivalents of 0.1 ]0.25 h. The event of January 21, 2011 was particularly well observed, with numerous surface reports available. It was also well represented by the WRF configuration operated at NWS Gaylord. Given that the default configuration produced a reasonable prediction, it is used here to evaluate the impacts of other physics configurations on the resulting prediction of the primary lake effect band and resulting QPF. Emphasis here is on differences in planetary boundary layer and cloud microphysics parameterizations, given their likely role in determining the evolution of shallow convection and precipitation processes. Results from an ensemble of seven microphysics schemes and three planetary boundary layer schemes are presented to demonstrate variability in forecast evolution, with results used in an attempt to improve the forecasts in the 2011 ]2012 lake effect season.

Molthan, Andrew↗
NASA NTRS2013

Using Field and Satellite Measurements to Improve Snow and Riming Processes in Cloud Resolving Models

The representation of clouds in climate and weather models is a driver in forecast uncertainty. Cloud microphysics parameterizations are challenged by having to represent a diverse range of ice species. Key characteristics of predicted ice species include habit and fall speed, and complex interactions that result from mixed-phased processes like riming. Our proposed activity leverages Global Precipitation Measurement (GPM) Mission ground validation studies to improve parameterizations

Colle, Brian A.↗
NASA NTRS2014

Assessment of NASA GISS CMIP5 and Post-CMIP5 Simulated Clouds and TOA Radiation Budgets Using Satellite Observations: TOA Radiation Budget and CREs - Part 2

In Part I of this study, the NASA GISS Coupled Model Intercomparison Project (CMIP5) and post-CMIP5 (herein called C5 and P5, respectively) simulated cloud properties were assessed utilizing multiple satellite observations, with a particular focus on the southern midlatitudes (SMLs). This study applies the knowledge gained from Part I of this series to evaluate the modeled TOA radiation budgets and cloud radiative effects (CREs) globally using CERES EBAF (CE) satellite observations and the impact of regional cloud properties and water vapor on the TOA radiation budgets. Comparisons revealed that the P5- and C5-simulated global means of clear-sky and all-sky outgoing longwave radiation (OLR) match well with CE observations, while biases are observed regionally. Negative biases are found in both P5- and C5-simulated clear-sky OLR. P5-simulated all-sky albedo slightly increased over the SMLs due to the increase in low-level cloud fraction from the new planetary boundary layer (PBL) scheme. Shortwave, longwave, and net CRE are quantitatively analyzed as well. Regions of strong large-scale atmospheric upwelling/downwelling motion are also defined to compare regional differences across multiple cloud and radiative variables. In general, the P5 and C5 simulations agree with the observations better over the downwelling regime than over the upwelling regime. Comparing the results herein with the cloud property comparisons presented in Part I, the modeled TOA radiation budgets and CREs agree well with the CE observations. These results, combined with results in Part I, have quantitatively estimated how much improvement is found in the P5-simulated cloud and radiative properties, particularly over the SMLs and tropics, due to the implementation of the new PBL and convection schemes.

Model comparison↗