SEARCH · Search NASA

Results for “hydrometeor”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Aircraft microwave observations and simulations of deep convection at 18-183 GHz

A radiative transfer model is used to simulate the aircraft multichannel microwave observations and to study a storm system observed on June 29, 1986 during the Cohmex program. The model, which is a revision of the model used by Wilheit et al. (1982), is described. The microwave upwelling brightness temperatures at the top of the atmosphere are calculated from hydrometeor profiles derived from ground-based radar data. The model is used to study various vertical hydrometeor phase profiles and particle size distributions. The results are compared with aircraft data from the Multispectral Cloud Radiometer and the Microwave Precipitation Radiometer. The results show that the brightness temperatures are directly influenced by the complex nature of hydrometeors, including size distribution and the mixture of ice and water phases.

Adler, R. F.↗

NASA NTRS1990

Aircraft microwave observations and simulations of deep convection from 18 to 183 GHz. II - Model results

In this June 29, 1986 case study, a radiative transfer model is used to simulate the aircraft multichannel microwave brightness temperatures presented in the Adler et al. (1990) paper and to study the convective storm structure. Ground-based radar data are used to derive hydrometeor profiles of the storm, based on which the microwave upwelling brightness temperatures are calculated. Various vertical hydrometeor phase profiles and the Marshall and Palmer (M-P, 1948) and Sekhon and Srivastava (S-S, 1970) ice particle size distributions are experimented in the model. The results are compared with the aircraft radiometric data. The comparison reveals that the M-P distribution well represents the ice particle size distribution, especially in the upper tropospheric portion of the cloud; the S-S distribution appears to better simulate the ice particle size at the lower portion of the cloud, which has a greater effect on the low-frequency microwave upwelling brightness temperatures; and that, in deep convective regions, significant supercooled liquid water (about 0.5 g/cu m) may be present up to the -30 C layer, while in less convective areas, frozen hydrometeors are predominant above -10 C level.

Yeh, Hwa-Young M.↗

NASA NTRS2001

TRMM Precipitation Radar and Microwave Imager Observations of Convective and Stratiform Rain and Their Theoretical Implications

Theoretical models of radiative transfer are developed to simulate the 85 GHz brightness temperature, T85, observed by the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer as a function of rain rate. These simulations are performed separately over regions of convective and stratiform rain. For the purpose of constructing vertical profiles of hydrometeors in these regions, guidance is taken from TRMM Precipitation Radar (PR) observations, and from recent investigations of LDR and Z(sub DR) measurements made by land-based polarimetric radars. We find that for a given rain rate, the extinction in 85 GHz due to hydrometeors above the freezing level is relatively weak in convective regions compared to that in stratiform regions. The hydrometeor profile above the freezing level responsible for this weak extinction in convective regions is inferred, from LDR and Z(sub DR) observations and theoretical considerations, to contain two layers: 1) a layer of two km thickness with mixed phase particles, liquid drops and graupel, directly above the freezing level, and 2) a layer of graupel extending from the top of the mixed-phase layer to the cloud top. Strong extinction in the stratiform regions is inferred to result from slowly-falling, low-density ice aggregates (snow) above the freezing level. These theoretical results are consistent with T85 measured by TMI and the rain rate deduced from PR for convective and stratiform rain regions. On the basis of this study we infer that the accuracy of the rain rate sensed by TMI depends critically on the specification of the convective or stratiform nature of the rain.

Prabhakara, C.↗

NASA NTRS2004

Classification of Tropical Oceanic Precipitation Using High Altitude Aircraft Microwave and Electric Field Measurements

A physically intuitive and computationally simple precipitation mapping algorithm has been developed for use with the airborne Advanced Microwave Precipitation Radiometer (AMPR). The algorithm is based on microwave emission and scattering properties of precipitation. Specifically, emission by liquid water allows increasing brightness temperatures at low frequencies to be interpreted as increasing rain rates. Scattering by large hydrometeors (particularly graupel and hail) causes relative minima in the brightness temperatures, with progressively larger hydrometeors scattering progressively longer wavelengths. The vigor of convection is therefore ascertained according to which wavelengths are being significantly scattered. The combination of emission and scattering information from four microwave channels is used to assign a precipitation category, which is related to the liquid rain rate, the vertical extent of precipitation, and the vigor of convection. The qualitative precipitation categories output by the passive microwave algorithm have been verified using coincident radar (ER-2 Doppler Radar - EDOP) and electric field measurements (Lightning Instrument Package - LIP). These coincident measurements can subsequently be used to quantify rain rates, hydrometeor contents, and vertical profiles that are typical for each precipitation category. This algorithm has been developed using an airborne platform. Comparisons are being made with other airborne, satellite, and ground-based radar and radiometer data. This technique shows promise both as a research tool and potentially as a real-time analysis tool, which could be applied to either traditional or uninhabited aerial vehicles.

Cecil, Daniel J.↗

NASA NTRS2004

Classification of Tropical Oceanic Precipitation Using High Altitude Aircraft Microwave and Electric Field Measurements

Cecil, Daniel J.↗

NASA NTRS2005

Interannual Variations in Tropical Upper-Tropospheric Humidity: Understanding Tropical Convective and Dynamical Processes

Uncertainty remains as to what extent variability in mid to upper tropospheric moisture, especially over the tropics, behaves as constant relative humidity during interannual climate variations associated with ENSO. Systematic variations in HIRS 6.7 micron and MLS 205 GHz suggest that dry subtropical regions evolving during warm SST events depress relative humidity, but the interpretation of these events is still uncertain. Additional specific concerns have to do with regional signatures of convective processes: How does the origin of dry air in the eastern subtropical N. Pacific differ in ENSO warm versus cold years? The dynamics of Rossby wave forcing by convective heating, subtropical jet stream dynamics, and dynamics driven subsidence all come into play here. How variations in precipitating ice hydrometeors from tropical anvils relate to variations in UTH is also a subject of debate? Do variations in precipitating ice, cloud cover and water vapor behavior show any support for the Iris-hypothesis mechanism? Here we examine historical records of SSM/T-2 data to gain a better physical understanding of the effects of deep convective moisture sources and dynamically-induced vertical circulations on UTH. These high frequency microwave measurements (183.3 GHz) take advantage of far less sensitivity to cloud hydrometeors than the 6.7 micron data to yield a record of upper tropospheric relative humidity. Furthermore, signatures of precipitating ice from these channels facilitate comparisons to TRMM hydrometeors detected by radar. In analyzing these observations, we isolate water vapor and temperature change components that affect brightness temperatures and the inferred relative humidity. Trajectory modeling is also used to understand interannual humidity anomalies in terms of outflow fbm convective regions and history of diabatically-driven sinking which modifies relative humidity.

Robertson, Franklin R.↗

NASA NTRS2005

On Modeling Air/Space-Borne Radar Returns in the Melting Layer

The bright band is the enhanced radar echo associated with the melting of hydrometeors in stratiform rain where the melting process usually occurs below 0 C isotherm over a distance of about 500m. To simulate this radar signature, a scattering model of melting snow is proposed in which the fractional water content is prescribed as a function of the radius of a spherical mixed- phase particle consisting of air, ice and water. The model is based on the observation that melting starts at the surface of the particle and then gradually develops towards the center. To compute the scattering parameters of a non-uniform melting particle, the particle is modeled as a sphere represented by a collection of 64(exp 3) cubic cells of identical size where the probability of water at any cell is prescribed as a function of the radius. The internal field of the particle, used for deriving the effective dielectric constant, is computed by the Conjugate Gradient and Fast Fourier Transform (CGFFT) numerical methods. To make computations of the scattering parameters more efficient, a multi-layer stratified-sphere scattering model is introduced after demonstrating that the scattering parameters of the non-uniformly melting particle can be accurately reproduced by the stratified sphere. In conjunction with a melting layer model that describes the melting fractions and fall velocities of hydrometeors as a function of the distance from the 0 C isotherm, the stratified-sphere model is used to simulate the radar bright band profiles. These simulated profiles are shown to compare well with measurements from the Precipitation Radar (PR) aboard the Tropical Rainfall Measuring Mission (TRMM) satellite and a dual-wavelength airborne radar. The results suggest that the proposed model of a melting snow particle may be useful in studying the characteristics of the bright-band in particular and mixed- phase hydrometeors in general.

Liao, Liang↗

NASA NTRS2006

Percentage Contributions from Atmospheric and Surface Features to Computed Brightness Temperatures

Over the past few years, a few solid precipitation detection and retrieval algorithms have been developed and shown to be applicable for snowing clouds and blizzards. Current precipitating snow retrieval algorithms require the use of millimeter-wave radiometer observations. The millimeter-wave frequencies are especially sensitive to the scattering and emission properties of frozen particles due to the ice particle refractive index. These channels can also be used to discern information about the frozen particles above the melting layer. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. This work will show the percentage of the brightness temperature resulting from the liquid hydrometeor, frozen hydrometeor, relative humidity, and surface contributions. The focus will be on precipitating snow events and millimeter-wave frequencies however, other events and frequencies will be included in the analysis.

SkofronickJackson, Gail↗

NASA NTRS2006

A Climatology of Midlatitude Continental Clouds from the ARM SGP Site: Cloud Fraction and Surface Radiative Forcing - Part II

Data collected at the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) Central Facility (SCF) are analyzed to determine the monthly and hourly variations of cloud fraction and radiative forcing between January 1997 and December 2002. Cloud fractions are estimated for total cloud cover and for single-layered low (0-3 km), middle (3-6 km), and high clouds (more than 6 km) using ARM SCG ground-based paired lidar-radar measurements. Shortwave (SW) and longwave (LW) fluxes are derived from up- and down-looking standard precision spectral pyranometers and precision infrared radiometer measurements with uncertainties of approximately 10 Wm(exp -2). The annual averages of total, and single-layered low, middle and high cloud fractions are 0.49, 0.11, 0.03, and 0.17, respectively. Both total and low cloud amounts peak during January and February and reach a minimum during July and August, high clouds occur more frequently than other types of clouds with a peak in summer. The average annual downwelling surface SW fluxes for total and low clouds (151 and 138 Wm(exp-2), respectively) are less than those under middle and high clouds (188 and 201 Wm(exp -2), respectively), but the downwelling LW fluxes (349 and 356 Wm(exp -2)) underneath total and low clouds are greater than those from middle and high clouds (337 and 333 Wm(exp -2)). Low clouds produce the largest LW warming (55 Wm(exp -2) and SW cooling (-91 Wm(exp -2)) effects with maximum and minimum absolute values in spring and summer, respectively. High clouds have the smallest LW warming (17 Wm(exp -2)) and SW cooling (-37 Wm(exp -2)) effects at the surface. All-sky SW CRF decreases and LW CRF increases with increasing cloud fraction with mean slopes of -0.984 and 0.616 Wm(exp -2)%(exp -1), respectively. Over the entire diurnal cycle, clouds deplete the amount of surface insolation more than they add to the downwelling LW flux. The calculated CRFs do not appear to be significantly affected by uncertainties in data sampling and clear-sky screening. Traditionally, cloud radiative forcing includes, not only the radiative impact of the hydrometeors, but also the changes in the environment. Taken together over the ARM SCF, changes in humidity and surface albedo between clear and cloudy conditions offset approximately 20% of the NET radiative forcing caused by the cloud hydrometeors alone. Variations in water vapor, on average, account for 10% and 83% of the SW and LW CRFs, respectively, in total cloud cover conditions. The error analysis further reveals that the cloud hydrometeors dominate the SW CRF, while water vapor changes are most important for LW flux changes in cloudy skies. Similar studies over other locales are encouraged where water and surface albedo changes from clear to cloudy conditions may be much different than observed over the ARM SCF.

Xi, B.↗

NASA NTRS2009

Developing Lightning Prediction Tools for the CCAFS Dual-Polarimetric Radar

NASA Marshall Space Flight Center and the University of Alabama Huntsville are collaborating with the 45th Weather Squadron (45WS) to develop improved lightning prediction capabilities for the new C-band dual-polarimetric weather radar being acquired for use by 45WS and launch weather forecasters at Cape Canaveral Air Force Station (CCAFS). In particular, these algorithms will focus on lightning onset, cessation and combined lightning-radar applications for convective winds assessment. Research using radar reflectivity (Z) data for prediction of lightning onset has been extensively discussed in the literature and subsequently applied by launch weather forecasters as it pertains to lightning nowcasting. Currently the forecasters apply a relatively straight forward but effective temperature-Z threshold algorithm for assessing the likelihood of lightning onset in a given storm. In addition, a layered VIL above the freezing level product is used as automated guidance for the onset of lightning. Only limited research and field work has been conducted on lightning cessation using Z and vertically-integrated Z for determining cessation. Though not used operationally vertically-integrated Z (basis for VIL) has recently shown promise as a tool for use in nowcasting lightning cessation. The work discussed herein leverages and expands upon these and similar reflectivity-threshold approaches via the application/addition of over two decades of polarimetric radar research focused on distinct multi-parameter radar signatures of ice/mixed-phase initiation and ice-crystal orientation in highly electrified convective clouds. Specifically, our approach is based on numerous previous studies that have observed repeatable patterns in the behavior of the vertical hydrometeor column as it relates to the temporal evolution of differential reflectivity and depolarization (manifested in either LDR or p(sub hv)), development of in-situ mixed and ice phase microphysics, electric fields, and ensuing lightning in the sub-tropical/tropical convection typical of the southeastern U.S., Maritime Continent, and southwestern Amazon. The polarimetric signatures detected in this setting provide a basis for automated 3-D detection of hydrometeor types in fuzzy logic hydrometeor identification algorithms (HID). Our working hypothesis is that improvement in lightning onset warning lead time and specificity for a given storm, relative to application of a Z-threshold algorithm, should arise as a consequence of the ability of dual-polarimetric radar to unambiguously detect and identify (through HID algorithms) the updraft elevation of rain-water cores above the freezing level and subsequent onset of drop freezing, riming, and robust mixed phase processes leading to significant charge separation and lightning. This type of algorithm, though dependent on the quality of the polarimetric data should be less susceptible to variable Z-calibration that can impact a given Z-threshold approach. To facilitate development of the algorithm while the 45WS dual-pol radar is in its current test stages and to evaluate the impact of polarimetric data quality (e.g., modified scan parameters and sampling) on the ensuing algorithms, we are using the ARMOR C-band dual-pol radar in Huntsville combined with N. Alabama LMA data and ARMOR HID algorithms [NCAR algorithm modified for application at C-band] in a testbed fashion. For lightning cessation we are revisiting the application of differential propagation phase variables for the monitoring of ice crystal alignment driven by in-cloud electric fields combined with metrics of ice water path (i.e., vertically integrated reflectivity). Importantly it should be noted that this approach is still very much a research topic and as such, we will explore operational applications that involve radar frequencies other than C-Band by using the UAH MAX X-band dual-pol radar in slow staring modes.

Petersen, W. A.↗

NASA NTRS2010

Comparing Aircraft Observations of Snowfall to Forecasts Using Single or Two Moment Bulk Water Microphysics Schemes

High resolution weather forecast models with explicit prediction of hydrometeor type, size distribution, and fall speed may be useful in the development of precipitation retrievals, by providing representative characteristics of frozen hydrometeors. Several single or double-moment microphysics schemes are currently available within the Weather Research and Forecasting (WRF) model, allowing for the prediction of up to three ice species. Each scheme incorporates different assumptions regarding the characteristics of their ice classes, particularly in terms of size distribution, density, and fall speed. In addition to the prediction of hydrometeor content, these schemes must accurately represent the vertical profile of water vapor to account for possible attenuation, along with the size distribution, density, and shape characteristics of ice crystals that are relevant to microwave scattering. An evaluation of a particular scheme requires the availability of field campaign measurements. The Canadian CloudSat/CALIPSO Validation Project (C3VP) obtained measurements of ice crystal shapes, size distributions, fall speeds, and precipitation during several intensive observation periods. In this study, C3VP observations obtained during the 22 January 2007 synoptic-scale snowfall event are compared against WRF model output, based upon forecasts using four single-moment and two double-moment schemes available as of version 3.1. Schemes are compared against aircraft observations by examining differences in size distribution, density, and content. In addition to direct measurements from aircraft probes, simulated precipitation can also be converted to equivalent, remotely sensed characteristics through the use of the NASA Goddard Satellite Data Simulator Unit. Outputs from high resolution forecasts are compared against radar and satellite observations emphasizing differences in assumed crystal shape and size distribution characteristics.

Molthan, Andrew L.↗

NASA NTRS2010

Estimating Glaciation Temperature of Deep Convective Clouds with Remote Sensing Data

Major uncertainties exist for observing and modeling ice content inside deep convective clouds (DCC). One of the difficulties has been the lack of characterization of vertical profiles of cloud hydrometeor phase. Here we propose a technique to estimate the DCC glaciation temperature using passive remote sensing data. It is based on a conceptual model of vertical hydrometeor size profiles inside DCCs. Estimates from the technique agree well with our general understanding of the problem. Furthermore, the link between vertical profiles of cloud particle size and hydrometeor thermodynamic phase is confirmed by a 3-13 cloud retrieval technique. The technique is applied to aircraft measurements of cloud side reflectance and the result was compared favorably with an independent retrieval of thermodynamic phase based on different refractive indices at 2.13 micron and 2.25 micron. Possible applications of the technique are discussed.

Yuan, Tianle↗

NASA NTRS2009

Spatial and Temporal Extrapolation of Disdrometer Size Distributions Based on a Lagrangian Trajectory Model of Falling Rain

Methodologies to improve disdrometer processing, loosely based on mathematical techniques common to the field of particle flow and fluid mechanics, are examined and tested. The inclusion of advection and vertical wind field estimates appear to produce significantly improved results in a Lagrangian hydrometeor trajectory model, in spite of very strict assumptions of noninteracting hydrometeors, constant vertical air velocity, and time independent advection during the scan time interval. Wind field data can be extracted from each radar elevation scan by plotting and analyzing reflectivity contours over the disdrometer site and by collecting the radar radial velocity data to obtain estimates of advection. Specific regions of disdrometer spectra (drop size versus time) often exhibit strong gravitational sorting signatures, from which estimates of vertical velocity can be extracted. These independent wind field estimates become inputs and initial conditions to the Lagrangian trajectory simulation of falling hydrometeors.

Lane, John E.↗

NASA NTRS2011

Surface and Atmospheric Contributions to Passive Microwave Brightness Temperatures for Falling Snow Events

Physically based passive microwave precipitation retrieval algorithms require a set of relationships between satellite -observed brightness temperatures (TBs) and the physical state of the underlying atmosphere and surface. These relationships are nonlinear, such that inversions are ill ]posed especially over variable land surfaces. In order to elucidate these relationships, this work presents a theoretical analysis using TB weighting functions to quantify the percentage influence of the TB resulting from absorption, emission, and/or reflection from the surface, as well as from frozen hydrometeors in clouds, from atmospheric water vapor, and from other contributors. The percentage analysis was also compared to Jacobians. The results are presented for frequencies from 10 to 874 GHz, for individual snow profiles, and for averages over three cloud-resolving model simulations of falling snow. The bulk structure (e.g., ice water path and cloud depth) of the underlying cloud scene was found to affect the resultant TB and percentages, producing different values for blizzard, lake effect, and synoptic snow events. The slant path at a 53 viewing angle increases the hydrometeor contributions relative to nadir viewing channels. Jacobians provide the magnitude and direction of change in the TB values due to a change in the underlying scene; however, the percentage analysis provides detailed information on how that change affected contributions to the TB from the surface, hydrometeors, and water vapor. The TB percentage information presented in this paper provides information about the relative contributions to the TB and supplies key pieces of information required to develop and improve precipitation retrievals over land surfaces.

Skofronick-Jackson, Gail↗

NASA NTRS2016

Deployment and Performance of the NASA D3R During the GPM OLYMPEx Field Campaign

The NASA D3R was successfully deployed and operated throughout the NASA OLYMPEx field campaign. A differential phase based attenuation correction technique has been implemented for D3R observations. Hydrometeor classification has been demonstrated for five distinct classes using Ku-band observations of both convection and stratiform rain. The stratiform rain hydrometeor classification is compared against LDR observations and shows good agreement in identification of mixed-phase hydrometeors in the melting layer.

Sea measurements↗

NASA NTRS2018

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗

NASA NTRS2021

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗

NASA NTRS2022

Implementation of a Discrete Dipole Approximation Scattering Database into Community Radiative Transfer Model

The Community Radiative Transfer Model (CRTM) is a fast model that requires bulk optical properties of hydrometeors in the form of lookup tables to simulate all-sky satellite radiances. Current cloud scattering lookup tables of CRTM were generated using the Mie-Lorenz theory thus assuming spherical shapes for all frozen habits, while actual clouds contain frozen hydrometeors with different shapes. The Discrete Dipole Approximation (DDA) technique is an effective technique for simulating the optical properties of non-spherical hydrometeors in the microwave region. This paper discusses the implementation and validation of a comprehensive DDA cloud scattering database into CRTM for the microwave frequencies. The original DDA database assumes total random orientation in the calculation of single scattering properties. The mass scattering parameters required by CRTM were then computed from single scattering properties and water content dependent particle size distributions. The new lookup tables eliminate the requirement for providing the effective radius as input to CRTM by using the cloud water content for the mass dimension. A collocated dataset of short-term forecasts from Integrated Forecast System of the European Centre for Medium-Range Weather Forecasts and satellite microwave data was used for the evaluation of results. The results overall showed that the DDA lookup tables, in comparison with the Mie tables, greatly reduce the differences among simulated and observed values. The Mie lookup tables especially introduce excessive scattering for the channels operating below 90 GHz and low scattering for the channels above 90 GHz.

Discrete dipole↗