Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice cloud models”

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 217 records · Page 12

Observations of Recent Arctic Sea Ice Volume Loss and Its Impact on Ocean-Atmosphere Energy Exchange and Ice Production

Using recently developed techniques we estimate snow and sea ice thickness distributions for the Arctic basin through the combination of freeboard data from the Ice, Cloud, and land Elevation Satellite (ICESat) and a snow depth model. These data are used with meteorological data and a thermodynamic sea ice model to calculate ocean-atmosphere heat exchange and ice volume production during the 2003-2008 fall and winter seasons. The calculated heat fluxes and ice growth rates are in agreement with previous observations over multiyear ice. In this study, we calculate heat fluxes and ice growth rates for the full distribution of ice thicknesses covering the Arctic basin and determine the impact of ice thickness change on the calculated values. Thinning of the sea ice is observed which greatly increases the 2005-2007 fall period ocean-atmosphere heat fluxes compared to those observed in 2003. Although there was also a decline in sea ice thickness for the winter periods, the winter time heat flux was found to be less impacted by the observed changes in ice thickness. A large increase in the net Arctic ocean-atmosphere heat output is also observed in the fall periods due to changes in the areal coverage of sea ice. The anomalously low sea ice coverage in 2007 led to a net ocean-atmosphere heat output approximately 3 times greater than was observed in previous years and suggests that sea ice losses are now playing a role in increasing surface air temperatures in the Arctic.

Kurtz, N. T.↗

Physical Models of Layered Polar Firn Brightness Temperatures from 0.5 to 2 GHz

We investigate physical effects influencing 0.5-2 GHz brightness temperatures of layered polar firn to support the Ultra Wide Band Software Defined Radiometer (UWBRAD) experiment to be conducted in Greenland and in Antarctica. We find that because ice particle grain sizes are very small compared to the 0.5-2 GHz wavelengths, volume scattering effects are small. Variations in firn density over cm- to m-length scales, however, cause significant effects. Both incoherent and coherent models are used to examine these effects. Incoherent models include a 'cloud model' that neglects any reflections internal to the ice sheet, and the DMRT-ML and MEMLS radiative transfer codes that are publicly available. The coherent model is based on the layered medium implementation of the fluctuation dissipation theorem for thermal microwave radiation from a medium having a nonuniform temperature. Density profiles are modeled using a stochastic approach, and model predictions are averaged over a large number of realizations to take into account an averaging over the radiometer footprint. Density profiles are described by combining a smooth average density profile with a spatially correlated random process to model density fluctuations. It is shown that coherent model results after ensemble averaging depend on the correlation lengths of the vertical density fluctuations. If the correlation length is moderate or long compared with the wavelength (approximately 0.6x longer or greater for Gaussian correlation function without regard for layer thinning due to compaction), coherent and incoherent model results are similar (within approximately 1 K). However, when the correlation length is short compared to the wavelength, coherent model results are significantly different from the incoherent model by several tens of kelvins. For a 10-cm correlation length, the differences are significant between 0.5 and 1.1 GHz, and less for 1.1-2 GHz. Model results are shown to be able to match the v-pol SMOS data closely and predict the h-pol data for small observation angles.

layered polar firn emission↗

Seasonal Variation in Martian Water Ice Cloud Particle Size

We employ the complete record of Mars limb-viewing observations by the Compact Reconnaissance Imaging Spectrometer, totaling 922 successful retrievals, to chart the seasonal variation in water ice cloud particle size. Low-latitude clouds exhibit particle size sorting with altitude at all seasons, with particles ranging from >3.0 m in effective radius to 1.01.5 m. However, at a given altitude, the ice particle size often varies through the year. Ice particle sizes in the tropics follow a different seasonal cycle than ice mixing ratio. North polar clouds have a complex seasonal cycle, with a sloped size profile with altitude as the polar hood cloud forms in early fall trending to a roughly uniform 1.5 m profile by the end of winter. These data are an important metric to validate and improve general circulation model microphysical routines and to better understand how martian water ice clouds influence the atmospheric radiative budget.

Scott D. Guzewich↗

Diabatic forcing and intialization with assimilation of cloud water and rainwater in a forecast model

In this study, diabatic forcing, and liquid water assimilation techniques are tested in a semi-implicit hydrostatic regional forecast model containing explicit representations of grid-scale cloud water and rainwater. Diabatic forcing, in conjunction with diabatic contributions in the initialization, is found to help the forecast retain the diabatic signal found in the liquid water or heating rate data, consequently reducing the spinup time associated with grid-scale precipitation processes. Both observational Special Sensor Microwave/Imager (SSM/I) and model-generated data are used. A physical retrieval method incorporating SSM/I radiance data is utilized to estimate the 3D distribution of precipitating storms. In the retrieval method the relationship between precipitation distributions and upwelling microwave radiances is parameterized, based upon cloud ensemble-radiative model simulations. Regression formulae relating vertically integrated liquid and ice-phase precipitation amounts to latent heating rates are also derived from the cloud ensemble simulations. Thus, retrieved SSM/I precipitation structures can be used in conjunction with the regression-formulas to infer the 3D distribution of latent heating rates. These heating rates are used directly in the forecast model to help initiate Tropical Storm Emily (21 September 1987). The 14-h forecast of Emily's development yields atmospheric precipitation water contents that compare favorably with coincident SSM/I estimates.

Raymond, William H.↗

The albedo of snow for partially cloudy skies

The input parameters of the model are atmospheric precipitable water, ozone content, turbidity, cloud optical thickness, size and shape of ice crystal of snow and surface pressure. The model outputs spectral and integrated solar flux snow reflectance as a function of solar elevation and fractional cloudcover. The model is illustrated using representative parameters for the Antarctic coastal regions. The albedo for a clear sky depends inversely on the solar elevation. At high elevation the albedo depends primarily upon the grain size; at low elevation this dependence is on grain size and shape. The gradient of the albedo-elevation curve increases as the grains get larger and faceted. The albedo for a dense overcast is a few percent higher than the clear sky albedo at high elevations. A simple relation between the grain size and the overcast albedo is obtained. For a set of grain size and shape, the albedo matrices (the albedo as a function of solar elevation and fractional cloudcover) are tabulated.

Choudhury, B. J.↗

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗

Microwave radiances from precipitating clouds containing aspherical ice, combined phase, and liquid hydrometeors

A numerical algorithm based on Eddington's second approximation to the equation of radiative transfer has been developed in order to compute the radiances with horizontal and vertical polarization that emerge from precipitating clouds. This algorithm yields a rapid solution to problems pertaining to clouds with vertically inhomogeneous structure. Precipitating clouds containing liquid, mixtures of phases, and ice hydrometeors are modelled. It is shown that the lower-frequency radiances are sensitive to liquid precipitation at low altitudes while the higher-frequency radiances are more sensitive to the ice hydrometeors at the cloud tops. The extinction coefficients of aspherical hydrometeors are presented as a function of rainfall rates. The fact that vertically polarized radiances are warmer than horizontally polarized ones at high rainfall rates, with the difference diminishing at lower frequencies, is due to aspherical ice hydrometeors in the upper regions of precipitating clouds.

Wu, R.↗

Modeling of Highly Instrumented Honeywell Turbofan Engine Tested with Ice Crystal Ingestion in the NASA Propulsion System Laboratory

The Propulsion Systems Laboratory (PSL), an altitude test facility at NASA Glenn Research Center, has been used to test a highly instrumented turbine engine at simulated altitude operating conditions. This is a continuation of the PSL testing that successfully duplicated the icing events that were experienced in a previous engine (serial LF01) during flight through ice crystal clouds, which was the first turbofan engine tested in PSL. This second model of the ALF502R-5A serial number LF11 is a highly instrumented version of the previous engine. The PSL facility provides a continuous cloud of ice crystals with controlled characteristics of size and concentration, which are ingested by the engine during operation at simulated altitudes. Several of the previous operating points tested in the LF01 engine were duplicated to confirm repeatability in LF11. The instrumentation included video cameras to visually illustrate the accretion of ice in the low pressure compressor (LPC) exit guide vane region in order to confirm the ice accretion, which was suspected during the testing of the LF01. Traditional instrumentation included static pressure taps in the low pressure compressor inner and outer flow path walls, as well as total pressure and temperature rakes in the low pressure compressor region. The test data was utilized to determine the losses and blockages due to accretion in the exit guide vane region of the LPC. Multiple data points were analyzed with the Honeywell Customer Deck. A full engine roll back point was modeled with the Numerical Propulsion System Simulation (NPSS) code. The mean line compressor flow analysis code with ice crystal modeling was utilized to estimate the parameters that indicate the risk of accretion, as well as to estimate the degree of blockage and losses caused by accretion during a full engine roll back point. The analysis provided additional validation of the icing risk parameters within the LPC, as well as the creation of models for estimating the rates of blockage growth and losses.

Simulation↗

The Cirrus Parcel Model Comparison Project

The cirrus Parcel Model Comparison Project involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 centimeters per second, respectively. Five models are participating in the project. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found amongst the models for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Systematic bias exists between results of a model based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies. This information is necessary, but not sufficient, to construct consistent formulae for the two approaches. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period.

Lin, Ruei-Fong↗

Case study of a low-reflectivity pulsating microburst - Numerical simulation of the Denver, 8 July 1989, storm

Results from a 3D simulation of a very strong microburst, detected by the Low-Level Windshear Alert System, are presented. The simulation was performed using the Terminal Area Simulation System (TASS) model (Proctor, 1987) which is a 3D nonhydrostatic cloud model that includes parameterizations for both liquid and ice-phase microphysics. Emphasis is placed om the pulsating characteristic and the very-low radar reflectivity of the event. The simulation revealed that the dissipating cloud induces a type-1 pulsating microburst event, with at least three distinct microburst pulses over a period of 20 min.

Proctor, Fred H.↗

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗

Anvil Glaciation in a Deep Cumulus Updraught over Florida Simulated with the Explicit Microphysics Model. I: Impact of Various Nucleation Processes

Simulations of a cumulonimbus cloud observed in the Cirrus regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) with an advanced version of the Explicit Microphysics Model (EMM) are presented. The EMM has size-resolved aerosols and predicts the time evolution of sizes, bulk densities and axial ratios of ice particles. Observations by multiple aircraft in the troposphere provide inputs to the model, including observations of the ice nuclei and of the entire size distribution of condensation nuclei. Homogeneous droplet freezing is found to be the source of almost all of the ice crystals in the anvil updraught of this particular model cloud. Most of the simulated droplets that freeze to form anvil crystals appear to be nucleated by activation of aerosols far above cloud base in the interior of the cloud ("secondary" or "in cloud" droplet nucleation). This is partly because primary droplets formed at cloud base are invariably depleted by accretion before they can reach the anvil base in the updraught, which promotes an increase with height of the average supersaturation in the updraught aloft. More than half of these aerosols, activated far above cloud base, are entrained into the updraught of this model cloud from the lateral environment above about 5 km above mean sea level. This confirms the importance of remote sources of atmospheric aerosol for anvil glaciation. Other nucleation processes impinge indirectly upon the anvil glaciation by modifying the concentration of supercooled droplets in the upper levels of the mixed-phase region. For instance, the warm-rain process produces a massive indirect impact on the anvil crystal concentration, because it determines the mass of precipitation forming in the updraught. It competes with homogeneous freezing as a sink for cloud droplets. The effects from turbulent enhancement of the warm-rain process and from the nucleation processes on the anvil ice properties are assessed.

Phillips, Vaughan T. J.↗

Impact of GMI All-Sky Radiance Assimilation in the NASA GEOS Forecast System

The assimilation of cloud- and precipitation-affected ("all-sky") radiances has become an important focus of development at most numerical weather prediction centers. Efforts at the Global Modeling and Assimilation Office (GMAO) have focused on all-sky assimilation of GPM Microwave Imager (GMI) radiances, which became operational in the GEOS real-time production system in July 2018. Implementation of the all-sky capability required several upgrades to the GEOS hybrid 4D-EnVar assimilation infrastructure including the addition of control variables for cloud liquid, cloud ice, rain and snow, enhancements to the radiative transfer model, new hybrid background and observational error models, and modified quality control and bias correction procedures. This talk describes the impact of GMI all-sky radiance assimilation on GEOS analyses and forecasts as determined from examination of various metrics including statistics of background departures and analysis increments, forecast skill scores, and forecast sensitivity observation impact (FSOI) calculations. It is shown that in addition to the hydrometeors themselves, the initial wind, temperature and pressure fields all undergo significant dynamic adjustment in response to the analyzed cloud and precipitation features. Assimilation of GMI radiances leads to improved forecasts of lower tropospheric wind, temperature and humidity, especially in the tropics. The largest forecast improvements occur during the first 48 hours, with diminishing impact thereafter. However, combining GMI all-sky assimilation with improvements to the GEOS model physics as in the recent implementation of the real-time production system, extends these forecast improvements well in to the medium range. FSOI results based on a 24-hr moist global energy norm show that GMI radiances provide nearly uniform beneficial impact throughout the tropics, with more mixed impacts in the subtropics. While the overall impact of GMI is smaller than that of other, much more numerous microwave and hyperspectral infrared radiance types, its impact is among the largest of all radiance types on a per-observation basis.

Gelaro, Ron↗

A Single-column Model Ensemble Approach Applied to the TWP-ICE Experiment

Single-column models (SCM) are useful test beds for investigating the parameterization schemes of numerical weather prediction and climate models. The usefulness of SCM simulations are limited, however, by the accuracy of the best estimate large-scale observations prescribed. Errors estimating the observations will result in uncertainty in modeled simulations. One method to address the modeled uncertainty is to simulate an ensemble where the ensemble members span observational uncertainty. This study first derives an ensemble of large-scale data for the Tropical Warm Pool International Cloud Experiment (TWP-ICE) based on an estimate of a possible source of error in the best estimate product. These data are then used to carry out simulations with 11 SCM and two cloud-resolving models (CRM). Best estimate simulations are also performed. All models show that moisture-related variables are close to observations and there are limited differences between the best estimate and ensemble mean values. The models, however, show different sensitivities to changes in the forcing particularly when weakly forced. The ensemble simulations highlight important differences in the surface evaporation term of the moisture budget between the SCM and CRM. Differences are also apparent between the models in the ensemble mean vertical structure of cloud variables, while for each model, cloud properties are relatively insensitive to forcing. The ensemble is further used to investigate cloud variables and precipitation and identifies differences between CRM and SCM particularly for relationships involving ice. This study highlights the additional analysis that can be performed using ensemble simulations and hence enables a more complete model investigation compared to using the more traditional single best estimate simulation only.

numerical weather forecasting↗

Inference of precipitation through thermal infrared measurements of soil moisture

The physics of microwave radiative transfer is well understood so that causal models can be assembled which relate the observed brightness temperatures to assumed distributions of hydrometeors (both liquid and ice), non-precipitating clouds, water vapor oxygen, and surface conditions. Present models assume a Marshall Palmer size distribution of liquid hydrometers from the surface to the freezing level (near the 0 C isotherm) and a variable thickness of frozen hydrometeors above that with various reasonable distribution of the other relevant constituents. The validity of such models is discussed. All uncertainties in the rain rate retrieval algorithms can be expressed in terms of specific model uncertainties which can be addressed through appropriate measurements. Those factors which must be known to achieve umambiguous results can be identified so that rainfall measuring algorithms can be developed and improved. The emissivity of the underlying surface significantly affects the contrast that may be measured between areas covered by rain and those which are dry. Sensing strategies for measuring rain over the ocean and rain over land are reviewed.

Wetzel, P. J.↗

Cirrus Clouds on Mars: Data Analysis and GCM Modeling

The goal of the investigation was to obtain a better retrieval of cirrus ice concentrations [to] gain a quantitative understanding of the cirrus clouds in the upper atmosphere by analyzing The Thermal Emission Spectrometer (TES) data obtained by the Mars Global Surveyor (MGS). We study the global distribution of dust and ice for a Martian year to address the fundamental question of the importance of ice in the Martian atmosphere.

Yung, Yuk L.↗

Cloud Influence on ERA5 and AMPS Surface Downwelling Longwave Radiation Biases in West Antarctica

The surface downwelling longwave radiation component (LW[down arrow]) is crucial for the determination of the surface energy budget and has significant implications for the resilience of ice surfaces in the polar regions. Accurate model evaluation of this radiation component requires knowledge about the phase, vertical distribution, and associated temperature of water in the atmosphere, all of which control the LW[down arrow] signal measured at the surface. In this study, we examine the LW[down arrow] model errors found in the Antarctic Mesoscale Prediction System (AMPS) operational forecast model and the ERA5 reanalysis model relative to observations from the AWARE campaign at McMurdo Station and the West Antarctic Ice Sheet (WAIS) Divide. The errors are calculated separately for observed clear-sky conditions, ice-cloud occurrences, and liquid-bearing cloud layer (LBCL) occurrences. The analysis results show a tendency in both models at each site to underestimate the LW[down arrow] during clear sky conditions, high error variability (standard deviations > 20 W/m[exp2]) during any type of cloud occurrence, and negative LW biases when LBCLs are observed (bias magnitudes > 15 W/m[exp2] in tenuous LBCL cases; > 43 W/m[exp2] in optically thick/opaque LBCLs instances). We suggest that a generally dry and liquid-deficient atmosphere responsible for the identified LW[down arrow] biases in both models is the result of excessive ice formation and growth, which could stem from model initial and lateral boundary conditions, microphysics scheme, aerosol representation, and/or limited vertical resolution.

Israel Silber↗

Improving the Representation of Snow Crystal Properties Within a Single-Moment Microphysics Scheme

As computational resources continue their expansion, weather forecast models are transitioning to the use of parameterizations that predict the evolution of hydrometeors and their microphysical processes, rather than estimating the bulk effects of clouds and precipitation that occur on a sub-grid scale. These parameterizations are referred to as single-moment, bulk water microphysics schemes, as they predict the total water mass among hydrometeors in a limited number of classes. Although the development of single moment microphysics schemes have often been driven by the need to predict the structure of convective storms, they may also provide value in predicting accumulations of snowfall. Predicting the accumulation of snowfall presents unique challenges to forecasters and microphysics schemes. In cases where surface temperatures are near freezing, accumulated depth often depends upon the snowfall rate and the ability to overcome an initial warm layer. Precipitation efficiency relates to the dominant ice crystal habit, as dendrites and plates have relatively large surface areas for the accretion of cloud water and ice, but are only favored within a narrow range of ice supersaturation and temperature. Forecast models and their parameterizations must accurately represent the characteristics of snow crystal populations, such as their size distribution, bulk density and fall speed. These properties relate to the vertical distribution of ice within simulated clouds, the temperature profile through latent heat release, and the eventual precipitation rate measured at the surface. The NASA Goddard, single-moment microphysics scheme is available to the operational forecast community as an option within the Weather Research and Forecasting (WRF) model. The NASA Goddard scheme predicts the occurrence of up to six classes of water mass: vapor, cloud ice, cloud water, rain, snow and either graupel or hail.

Molthan, Andrew L.↗