Search NASA⌕ Search

SEARCH · Search NASA

Results for “cloud modeling”

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 541 records · Page 30

Extended atmospheres of outer planet satellites and comets

Model analysis of the extended atmospheres of outer planet satellites and comets are discussed. Understanding the neutral hydrogen distribution in the Saturn system concentrated on assessing the spatial dependence of the lifetime of hydrogen atoms and on obtaining appropriately sorted Lyman ALPHA data from the Voyager 1 UVS instrument. Progress in the area of the extended cometary atmospheres included analysis of Pioneer Venus Layman alpha observations of Comet P/Encke with the fully refined hydrogen cloud model, development of the basic carbon and oxygen models, and planning for the Pioneer Venus UVS observations of Comets P/Giacobini-Zinner and P/Halley.

Smyth, W. H.↗

Two space scatterer formalism calculation of bulk parameters of thunderclouds

In a previous study, we used a modified two-space scatterer formalism of Twersky to establish for a cloud modeled as a statistically homogeneous distribution of spherical water droplets, the dispersion relations that determine its bulk propagation numbers and bulk indexes of refraction in terms of the vector equivalent scattering amplitude and the dyadic scattering amplitude of the single water droplet in isolation. The results were specialized to the forward direction of scattering while demanding that the scatterers preserve the incident polarization. We apply this approach to obtain specific numerical values for the macroscopic parameters of the cloud. We work with a cloud of density rho = 100 cm(exp -3), a wavelength lambda = 0.7774 microns, and with spherical water droplets of common radius alpha = 10 microns. In addition, the scattering medium is divided into three parts, the medium outside the cloud, moist air (the medium inside the cloud but outside the droplets), and the medium inside the spherical water droplets. The results of this report are applicable to a cloud of any geometry since the boundary does not interfere with the calculations. Also, it is important to notice the plane wave nature of the incidence wave in the moist atmosphere.

Phanord, Dieudonne D.↗

Supersaturation, buoyancy, and deep convection dynamics

Motivated by recent discussions concerning differences of convective dynamics in polluted and pristine environments, the so-called convective invigoration in particular, this paper provides an analysis of factors affecting convective updraft buoyancy, such as the in-cloud supersaturation, condensate and precipitation loading, and entrainment. We use the deep convective period from simulations of daytime convection development over land discussed in our previous publications. An entraining parcel framework is used in the theoretical analysis. We show that for the specific case considered here, finite (positive) supersaturation noticeably reduces pseudo-adiabatic parcel buoyancy and cumulative convective available potential energy (cCAPE) in the lower troposphere. This comes from keeping a small fraction of the water vapor in a supersaturated state and thus reducing the latent heating. Such a lower-tropospheric impact is comparable to the effects of condensate loading and entrainment in the idealized parcel framework. For the entire tropospheric depth, loading and entrainment have a much more significant impact on the total CAPE. For the cloud model results, we compare ensemble simulations applying either a bulk microphysics scheme with saturation adjustment or a more comprehensive double-moment scheme with supersaturation prediction. We compare deep convective updraft velocities, buoyancies, and supersaturations from all ensembles. In agreement with the parcel analysis, the saturation-adjustment scheme provides noticeably stronger updrafts in the lower troposphere. For the simulations predicting supersaturation, there are small differences between pristine and polluted conditions below the freezing level that are difficult to explain by standard analysis of the in-cloud buoyancy components. By applying the piggybacking technique, we show that the lower-tropospheric buoyancy differences between pristine and polluted simulations come from a combination of temperature (i.e., latent heating) and condensate loading differences that work together to make polluted buoyancies and updraft velocities slightly larger when compared to their pristine analogues. Overall, the effects are rather small and contradict previous claims of a significant invigoration of deep convection in polluted environments.

54 ENVIRONMENTAL SCIENCES↗

The influence of multiple groups of biological ice nucleating particles on microphysical properties of mixed-phase clouds observed during MC3E

Abstract. A new empirical parameterization (EP) for multiple groups of primary biological aerosol particles (PBAPs) is implemented in the aerosol–cloud model (AC) to investigate their roles as ice nucleating particles (INPs). The EP describes the heterogeneous ice nucleation by (1) fungal spores, (2) bacteria, (3) pollen, (4) detritus of plants, animals, and viruses, and (5) algae. Each group includes fragments from the originally emitted particles. A high-resolution simulation of a midlatitude mesoscale squall line by AC is validated against airborne and ground observations. Sensitivity tests are carried out by varying the initial vertical profiles of the loadings of individual PBAP groups. The resulting changes in warm and ice cloud microphysical parameters are investigated. The changes in warm microphysical parameters, including liquid water content and cloud droplet number concentration, are minimal (<10 %). Overall, PBAPs have little effect on the ice number concentration (<6 %) in the convective region. In the stratiform region, increasing the initial PBAP loadings by a factor of 1000 resulted in less than 40 % change in ice number concentrations. The total ice concentration is mostly controlled by various mechanisms of secondary ice production (SIP). However, when SIP is intentionally shut down in sensitivity tests, increasing the PBAP loading by a factor of 100 has an effect of less than 3 % on the ice phase. Further sensitivity tests revealed that PBAPs have little effect on surface precipitation and on the shortwave and longwave flux (<4 %) for a 100-fold perturbation in PBAPs.

54 ENVIRONMENTAL SCIENCES↗

Multiparameter radar study of a microburst - Comparison with model results

Radar observations and model results are used to investigate the microphysical evolution of an isolated, intense storm observed on July 20, 1988 during the Microburst and Severe Thunderstorm experiment. The storm grew to a height of 14 km and upon collapsing, produced heavy rain, pea-sized hail, and a microburst at the surface. The radar observations indicate that the initial precipitation development was by collision-coalescence. As the storm intensified, accretional growth became dominant leading to rapid precipitation development. Radar-derived rainfall rates peaked around 150 to 190 mm/h. Each morning during the experiment, a two-dimensional, time-dependent cloud model, initialized with the morning sounding, was run. The model results from the July 20 sounding are compared to the radar observations. Good agreement is shown in some aspects of the storm development, although the numerical simulation predicted a more vigorous storm than actually developed.

Tuttle, John D.↗

The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era

Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope (JWST) era, high-resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constraints on the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models. In this paper, we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide-separation giant planets with key molecular species such as CH 4 , H 2 O, CO, and NH 3 in disequilibrium. Our grid covers atmospheres with T eff ∈ [500 K, 1300 K], log g ∈ [3.0, 5.5] (cgs) and an eddy diffusion parameter of logK zz = 2,4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH 4 , H 2 O, CO, and NH 3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗

Breakups are complicated: an efficient representation of collisional breakup in the superdroplet method

Abstract. A key constraint of particle-based methods for modeling cloud microphysics is the conservation of total particle number, which is required for computational tractability. The process of collisional breakup poses a particular challenge to this framework, as breakup events often produce many droplet fragments of varying sizes, which would require creating new particles in the system. This work introduces a representation of collisional breakup in the so-called “superdroplet” method which conserves the total number of superdroplets in the system. This representation extends an existing stochastic collisional-coalescence scheme and samples from a fragment size distribution in an additional Monte Carlo step. This method is demonstrated in a set of idealized box model and single-column warm-rain simulations. We further discuss the effects of the breakup dynamic and fragment size distribution on the particle size distribution, hydrometeor population, and microphysical process rates. Box model experiments serve to characterize the impacts of properties such as coalescence efficiency and fragmentation function on the relative roles of collisional breakup and coalescence. The results demonstrate that this representation of collisional breakup can produce a stationary particle size distribution, in which breakup and coalescence rates are approximately equal, and that it recovers expected behavior such as a reduction in precipitate-sized particles in the column model. The breakup algorithm presented here contributes to an open-source pythonic implementation of the superdroplet method, PySDM, which will facilitate future research using particle-based microphysics.

de Jong, Emily (ORCID:0000000253104554)↗

Carbon Dioxide Convection in the Martian Polar Night and Its Implications for Polar Processes

Each Martian year nearly 30% of the atmosphere is exchanged with the polar ice caps. This exchange occurs through a combination of direct surface condensation and atmospheric precipitation of carbon dioxide. It has long been thought the amount of condensation within the polar night is maintained by a balance between diabatic processes such as radiative cooling and latent heating from condensing CO2. This assumption manifests itself in Mars General Circulation Models (GCM) in such a way as to never allow the atmospheric temperature to dip below the saturation temperature of CO2. However, observations from Mars Global Surveyor (MGS) Radio Science (RS) and the Thermal Emission Spectrometer (TES) have demonstrated this assumption to be, at best, approximate. Both RS and TES observations within the polar nights of both poles indicate substantial supersaturated regions with respect to CO2. The observed temperature profiles suggest conditionally unstable regions containing planetary significant amounts of potential convective energy. Presented here are estimates of the total planetary inventory of convective available potential energy (CAPE) and the potential convective energy flux (PCEF). The values for CAPE and PCEF are derived from RS temperature profiles and compared to Mars GCM results using a new convective CO2 cloud model that allows for the formation of CAPE.

Colaprete, A.↗

Graupel and Snow Determination in Convection Using Dual-Wavelength Radar

Convective storms contain significant amounts of snow and graupel that each have different scattering properties for radar and radiometer remote sensing. TRMM TMI algorithms have relied heavily on radiative transfer calculations based on cloud model predicted microphysics information. In general, without in situ aircraft measurements, it is very difficult to separate snow and graupel using remote measurements. At the 94 GHZ frequency, snow is weakly attenuating while graupel is strongly attenuating. With a non-attenuating frequency such as X-band, one can distinguish Mie scattering from larger hydrometeors as well as attenuation from graupel. The presence of supercooled water will also cause attenuation of the 94 GHz radar signal, but this is mainly in convective updrafts. We will present preliminary results aimed at distinguishing snow from graupel using dual-wavelength (9.6 GHz, 94 GHz) measurements from aircraft overflights of convective storms in Florida. Measurements were made with the 9.6 GHz ER-2 Doppler Radar (EDOP) and the 94 GHz Cloud Radar System (CRS) over a number of different precipitation systems. Analysis of the vertical profiles at the two wavelengths provides information on snow and graupel regions, as well as supercooled water. The presentation will include examples from different cases and discussion of the results.

Heymsfield, Gerald↗

The Updated GEO Population for ORDEM 3.1

The limited availability of data for satellite fragmentations and debris in the geosynchronous orbit (GEO) region creates challenges to building accurate models for the orbital debris environment at such altitudes. Updated methods to properly incorporate and extrapolate measurement data have become a cornerstone of the GEO component in the newest version of the NASA Orbital Debris Engineering Model (ORDEM), ORDEM 3.1. For the GEO region, the Space Surveillance Network (SSN) catalog provides coverage down to a limit of approximately 1 m. A more statistically complete representation of the GEO population for smaller objects, which can pose a high risk to operational spacecraft, is thus dependent on dedicated observations by instruments optimized to observe debris smaller than the SSN cataloging threshold. For ORDEM 3.1, optical data from the Michigan Orbital DEbris Survey Telescope (MODEST) provided the input for building the GEO population down to approximately 30 cm (converting absolute magnitude to size). For smaller sizes, the size distribution of debris in the MODEST dataset was extrapolated down to 10 cm, and orbital parameters were estimated based on the orbits of the larger objects. When compared to previous versions of the model, significant improvements were made to the process of building the GEO population in ORDEM 3.1, both in the assessment of fragmentation debris in the data and assignment of orbital elements within the model. A so-called “debris ring filter,” based on a range of angles between an orbit’s angular momentum vector and that of the stable Laplace plane, was applied to the data to reduce biases from non- GEO objects, such as objects in a GEO-transfer orbit. In addition, a new approach was implemented to assign noncircular mean motions and eccentricities to the fragmentation debris observed by MODEST because the short observation window (5 min) in GEO limits orbit resolution to a circular orbit assumption for assigning orbital parameters. For ORDEM 3.1, non-circular orbital elements were assigned using relationships that were identified between mean motion and the angle between the orbit plane and the stable Laplace plane, as well as between mean motion and eccentricity, based on breakup clouds modeled by the NASA Standard Breakup Model. This approach has yielded a high-fidelity GEO model that has been validated with data from more recent MODEST observation campaigns.

Manis, A.↗

Towards a program of record of inland water quality: Exploiting present and heritage multispectral sensors for maximum information extraction

Degradation of Earth’s inland water resources due to anthropogenic perturbations and climate anomalies at both local and global scales continues to place human health at substantial risk. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This research exploits recent advancements in bio-optical modeling, cloud computing, and machine learning to enhance our capacity to leverage present and heritage satellite data. Recent research suggests that sensors with low spectral resolution, such as Sentinel 2 and Landsat missions, contain enough hidden spectral variation which can be exploited using data-driven approaches. The availability of three decades of archival imagery will open doors to discover global trends of eutrophication and increased cyanobacteria dominance and provide valuable insight to the development of predictive methodologies. Preliminary efforts in synthetic emulation of global natural inland waters will be discussed and contextualized against satellite radiometric measurement uncertainty, satellite data product uncertainties and causal signal ambiguity over the visible wavelength range, supported by high quality field and image data for selected inland aquatic sites. Insights on water quality estimation via data-driven machine learning models versus matrix inversions will be discussed, and how we can exploit spectral-spatial relationships in high spatial resolution data. A cross-sensor synergistic approach with detailed uncertainty analysis based on optical water types, will allow for unprecedented global snapshots of fine scale ecological dynamics of inland waters.

Inland↗

A model of the galactic tidal interaction with the Oort comet cloud

A model of the in situ Oort cloud which is isotropic with a random distribution of perihelia directions and angular momenta is considered. The energy distribution adopted has a continuous range of values appropriate for long-period (greater than 200 yr) comets. The evolution of 1,000,000 comets are tracked as they are made observable by the galactic tidal torque. Detailed results are obtained for the predicted distribution of new comets. Predictions of observed distributions are presented and compared with the random in situ results as well as with the actual observed distributions of class I comets. The predictions are in reasonable agreement with actual observations and, in many cases, are significantly different from random when perihelia directions are separated into galactic northern and southern hemispheres.

Matese, John J.↗

Adaptive Grid Redistribution for a 1D Model of Turbulence and Clouds

In global atmospheric models, resolving stratocumulus (Sc) in the vertical is computationally expensive. However, Sc appear only under special meteorological conditions. Therefore, there is motivation to refine the vertical grid levels adaptively. In order to facilitate the possibility of parallelization on graphical processing units, our grid adaptation method prescribes the number of vertical levels a priori. Then grid levels are relocated toward altitude ranges in need of refinement. Because the method relocates existing grid levels, rather than adding extra levels, there is a risk of creating regions with overly coarse grid spacing, that is, voids in the grid mesh. To prevent such voids from forming, a simple method is developed to impose a maximum grid spacing. To decide where to place enhanced resolution, the authors develop an empirical mesh refinement criterion. It refines grid spacing near the ground, near strong temperature gradients, and within clouds. Our grid adaptation method is implemented in a single-column model and evaluated on four test cases: decaying stratocumulus, developing shallow cumulus, a quasi-stationary stratocumulus deck, and the diurnal cycle of a dry boundary layer. In the stratocumulus cases, mesh refinement leads to improvements in both the time evolution of fields and their time averages. The other two cases show smaller differences.

Carstensen, Steffen [Univ. of Wisconsin, Milwaukee↗

The Impact on Simulated Storm Structure and Intensity of Variations in the Lifted Condensation Level and the Level of Free Convection

The sensitivities of convective storm structure and intensity to changes in the altitudes of the prestorm environmental lifted condensation level and level of free convection axe studied using a full-physics three-dimensional cloud model. Matrices of simulations are conducted for a range of LCL=LFC altitudes, using a single moderately-sheared curved hodograph trace in conjunction with 1 convective available potential energy values of either 800 or 2000 J/kg, with the matrices consisting of all four combinations of two distinct choices of buoyancy and shear profile shape. For each value of CAPE, the LCL=LFC altitudes are also allowed to vary in a series of simulations based on the most highly compressed buoyancy and shear profiles for that CAPE, with the environmental buoyancy profile shape, subcloud equivalent potential temperature, subcloud lapse rates of temperature and moisture, and wind profile held fixed. For each CAPE, one final simulation is conducted using a near optimal LFC, but a lowered LCL, with a neutrally buoyant environmental thermal profile specified in between. Results show that, for the buoyancy-starved small-CAPE environments, the simulated storms are supercells and are generally largest and most intense when LCL=LFC altitudes lie in the approximate range 1.5-2.5 km above the surface. The simulations show similar trends for the shear-starved large-CAPE environments, except that conversion from supercell to multicell morphology frequently occurs when the LCL is high. For choices of LCL=LFC height within the optimal 1.5-2.5 km range, peak storm updraft overturning efficiency may approaches unity relative to parcel theory, while for lower LCL=LFC heights, overturning efficiency is reduced significantly. The enhancements of overturning efficiency and updraft diameter with increasing LFC height are shown to be the result of systematic increases in the mean equivalent potential temperature of the updraft at cloud base. For the shear-starved environments, the tendency for outflow dominance is eliminated, but a large overturning efficiency maintained, when a low LCL is used in conjunction with a high LFC. The result regarding outflow dominance at high LCL is consistent with expectations, but the beneficial effect of a high LFC on convective overturning efficiency has not previously been widely recognized. The simulation findings here also appear to be consistent with statistics from previous severe storm environment climatologies, but provide a new framework for interpreting those statistics.

McCaul, Eugene W., Jr.↗

Model of a rotating magnetic cloud

The possibility that magnetic clouds rotate while they propagate antisunward was investigated. Magnetic clouds are modeled as magnetic flux ropes which rotate rigidly about the axis of symmetry. An ideal magnetohydrodynamic model, in which the evolution of the magnetic structure is related to the time evolution of the angular frequency, is developed. A class of 'separable' magnetic fields is employed to reduce the problem to a nonlinear ordinary differential equation for the evolution function, and it is solved numerically. The corresponding effective potential gives rise to two modes of evolution--expansion and oscillation--depending on the energy and on the value of a dimensionless parameter, k. Parameter k depends on the gas pressure, the ratio of the magnetic field components, and the frequency of rotation. There is a critical value of k, k(sub c), above which the oscillatory regime disappears and the flux rope invariably expands, regardless of the energy. Below k(sub c) the energy determines whether the configuration is confined or unbounded. Rotation always helps expansion by lowering the potential barrier. A data example was studied and features which are interpreted as signatures of rotation are presented. The angular speed is comparable to the Alfven speed, and the core of the rotating cloud completes on average one full revolution every three days at 1 AU. The parameter k is calculated from observations, and it is found to be close to, but below, critical. Only three out of the nine clouds examined showed signatures of rotation. Theoretical analysis suggests that close to the Sun rotation effects may play a more important role in the evolution of magnetic clouds than 1 AU.

Farrugia, C. J.↗

Pragmatic protocols for working cleanly when measuring ice nucleating particles

We report that measuring ice nucleating particles (INPs) is critical for understanding, and modeling, cloud formation, reflectivity, and precipitation patterns. However, because INPs are very rare in the atmosphere, but abundant—sometimes alarmingly so—in the dust that covers all work surfaces, contamination during processing of collected aerosol samples can be a significant hindrance to obtaining accurate measures of INPs. In preparing this technical note, we questioned the cleanliness of every collection and processing step involved in making immersion freezing measurements of INPs re-suspended from filter samples. The aim was to identify, and then minimize, all potential sources of significant contamination, including containers and tools used to store filters and prepare liquid suspensions, gloves, work surfaces, and the polymerase chain reaction (PCR) trays commonly used for sample analysis. While plasticware released few INPs, most gloves readily shed them, but this can be mitigated by using washed cleanroom vinyl or polyethylene gloves. Work surfaces, even those that had been cleaned, were prodigious reservoirs of contaminating INPs, but simply covering them with aluminum foil will provide an INP-free surface. By applying these practices, we developed a method to reduce the background level of INPs on Nuclepore™ polycarbonate filters in our tests to 0 at -25 °C and < 20 at -27 °C, making them suitable for use in all sampling environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Pressure Drag on Rising Buoyant Thermals in a Neutrally Stable Environment

Abstract This study examines dynamic pressure drag on rising dry buoyant thermals. A theoretical expression for drag coefficient C d as a function of several other nondimensional parameters governing thermal dynamics is derived based on combining the thermal momentum budget with the similarity theory of Scorer. Using values for these nondimensional parameters from previous studies, the theory suggests drag on thermals is small relative to that on solid spheres in laminar or turbulent flow. Two sets of numerical simulations of thermals in an unstratified, neutrally stable environment using an LES configuration of the Cloud Model 1 (CM1) are analyzed. One set has a relatively low effective Reynolds number R e and the other has an order-of-magnitude-higher R e ; these produce laminar and turbulent resolved-scale flows, respectively. Consistent with the theoretical C d , the magnitude of drag is small in all simulations. However, whereas the laminar thermals have C d ≈ 0.01, the turbulent thermals have weakly negative drag ( C d ≈ −0.1). This difference is explained by the laminar thermals having near vertical symmetry but the turbulent thermals exhibiting considerable vertical asymmetry of their azimuthally averaged flows. In the laminar thermals, buoyancy rapidly becomes concentrated around the main centers of rotation located along the horizontal central axis, leading to expansion of thermals via baroclinic vorticity generation but doing little to break vertical symmetry of the flow. Vertical asymmetry of the azimuthally averaged flow of turbulent thermals is attributed mainly to small-scale resolved eddies that are concentrated in the upper part of the thermals.

Meteorology & Atmospheric Sciences↗

ARES I: WASP-76 b, A Tale of Two HST Spectra

We analyze the transmission and emission spectra of the ultra-hot Jupiter WASP-76 b, observed with the G141 grism of the Hubble Space Telescope’s (HST) Wide Field Camera 3 (WFC3). We reduce and fit the raw data for each observation using the open-source software Iraclis before performing a fully Bayesian retrieval using the publicly available analysis suite TauREx 3. Previous studies of the WFC3 transmission spectra of WASP-76 b found hints of titanium oxide (TiO) and vanadium oxide (VO) or non-gray clouds. Accounting for a fainter stellar companion to WASP-76, we reanalyze this data and show that removing the effects of this background star changes the slope of the spectrum, resulting in these visible absorbers no longer being detected, eliminating the need for a non-gray cloud model to adequately fit the data but maintaining the strong water feature previously seen. However, our analysis of the emission spectrum suggests the presence of TiO and an atmospheric thermal inversion, along with a significant amount of water. Given the brightness of the host star and the size of the atmospheric features, WASP-76 b is an excellent target for further characterization with HST, or with future facilities, to better understand the nature of its atmosphere, to confirm the presence of TiO and to search for other optical absorbers.

79 ASTRONOMY AND ASTROPHYSICS↗