Search NASA⌕ Search

Engineering topics

S Joseph Munchak

Publications and source records attributed to S Joseph Munchak.

A Comparative Evaluation of Snowflake Particle Shape Estimation Techniques Used By the Precipitation Imaging Package (PIP), Multi-Angle Snowflake Camera (MASC), and Two-Dimensional Video Disdrometer (2DVD)

Measurements of snowflake particle shape are important for studying snow microphysics. While a number of instruments exist that are capable of measuring particle shape, this study focuses on the measurement techniques of three digital video disdrometers: the Precipitation Imaging Package (PIP), the Multi-Angle Snowflake Camera (MASC), and the Two-Dimensional Video Disdrometer (2DVD). To gain a better understanding of the relative strengths and weaknesses of these instruments and to provide a foundation upon which comparisons can be made between studies using data from different instruments, we perform a comparative analysis of the shape measurement algorithms employed by each of the three instruments by applying the algorithms to snowflake images captured by PIP during the ICE-POP 2018 field campaign. Our analysis primarily focuses on the measurement of the aspect ratio of either the particle itself, in the case of PIP and MASC, or of the particle bounding box, in the case of PIP and 2DVD. Both PIP and MASC use shape-fitting algorithms to measure aspect ratio. While our analysis of the MASC aspect ratio suggests that the measurements are reliable, our findings indicate that both the ellipse and rectangle aspect ratios produced by PIP underperformed considerably due to the shortcomings of the PIP shape-fitting techniques. We also demonstrate that reliable measurements of aspect ratio can be retrieved from PIP by reprocessing the raw PIP images using either the MASC ellipse-fitting algorithm or a tensor-based ellipse-fitting algorithm. Because of differences in instrument design, 2DVD produces measurements of particle horizontal and vertical extent rather than length and width. Furthermore, the 2DVD measurements of particle horizontal extent can be contaminated by horizontal particle motion. Our findings indicate that, although the correction technique used to remove the horizontal motion contamination performs remarkably well with snowflakes despite being designed for use with raindrops, the 2DVD measurements of particle horizontal extent are less reliable than those measured by PIP.

Charles N Helms↗

Snow Microphysical Retrieval from the NASA D3R Radar During ICE-POP 2018

A method is developed to use both polarimetric and dual-frequency radar measurements to retrieve microphysical properties of falling snow. It is applied to the Ku- and Ka-band measurements of the NASA dual-polarization, dual-frequency Doppler radar (D3R) obtained during the International Collaborative Experiments for PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign and incorporates the Atmospheric Radiative Transfer Simulator (ARTS) microwave singles cattering property database for oriented particles. The retrieval uses optimal estimation to solve for several parameters that describe the particle size distribution (PSD), relative contribution of pristine, aggregate, and rimed ice species, and the orientation distribution along an entire radial simultaneously. Examination of Jacobian matrices and averaging kernels shows that the dual-wavelength ratio (DWR) measurements provide information regarding the characteristic particle size, and to a lesser extent, the rime fraction and shape parameter of the size distribution, whereas the polarimetric measurements provide information regarding the mass fraction of pristine particles and their characteristic size and orientation distribution. Thus, by combining the dual-frequency and polarimetric measurements, some ambiguities can be resolved that should allow a better determination of the PSD and bulk microphysical properties (e.g., snowfall rate) than can be retrieved from single-frequency polarimetric measurements or dual-frequency, single-polarization measurements. The D3R ICE-POP retrievals were validated using Precipitation Imaging Package (PIP) and Pluvio weighing gauge measurements taken nearby at the May Hills ground site. The PIP measures the snow PSD directly, and its measurements can be used to derived the snowfall rate (volumetric and water equivalent), mean volume-weighted particle size, and effective density, as well as particle aspect ratio and orientation. Four retrieval experiments were performed to evaluate the utility of different measurement combinations: Kuonly, DWR-only, Ku-pol, and All-obs. In terms of correlation, the volumetric snowfall rate (r D 0:95) and snow water equivalent rate (r D 0:92) were best retrieved by the Ku-pol method, while the DWR-only method had the lowest magnitude bias for these parameters (-31% and -8 %, respectively). The methods that incorporated DWR also had the best correlation to particle size (r D 0:74 and r D 0:71 for DWR-only and All-obs, respectively), although none of the methods retrieved density particularly well (r D 0:43 for Allobs). The ability of the measurements to retrieve mean aspect ratio was also inconclusive, although the polarimetric methods (Ku-pol and All-obs) had reduced biases and mean absolute error (MAE) relative to the Ku-only and DWR-only methods. The significant biases in particle size and snowfall rate appeared to be related to biases in the measured DWR, emphasizing the need for accurate DWR measurements and frequent calibration in future D3R deployments.

S Joseph Munchak↗

Linkage among Ice Crystal Microphysics, Mesoscale Dynamics and Cloud and Precipitation Structures Revealed By Collocated Microwave Radiometer and Multi-Frequency Radar Observations

Ice clouds and falling snow are ubiquitous globally and play important roles in the Earth’s radiation budget and precipitation processes. Ice particle microphysical properties (e.g., size, habit and orientation) are not only influenced by ambient environment’s dynamic and thermodynamic conditions, but also intimately connect to the cloud radiative effects and particle fall speeds, which therefore impact up to the future climate projection and down to the details of the surface precipitation (e.g., onset-time, location, type and strength). Our previous work revealed that high-frequency Polarimetric radiance Difference (PD) from passive microwave sensors is a good indicator of the bulk aspect ratio of horizontally oriented ice particles that are often occur inside anvil clouds and/or stratiform precipitations. In this current work, we further investigate the dynamic/thermodynamic mechanisms and cloud/precipitation structures associated with ice-phase microphysics corresponding to different PD signals. In order to do so, collocated CloudSat radar (W-band) and Global Precipitation Measurement Dual-frequency Precipitation Radar (GPM-DPR, Ku/Ka bands) observations as well as European Centre for Medium-Range Weather Forecasts (ECMWF) atmosphere background profiles are grouped according to the magnitude of PD for only stratiform precipitation and/or anvil cloud scenes. We found that horizontally-oriented snow aggregates or large snow particles are likely the major contributor to the high-PD signals at 166 GHz, while low-PD magnitudes can be attributed to small cloud ice, randomly oriented snow aggregates, riming snow or super-cooled water. Further, high (low) PD scenes are found to be associated with stronger (weaker) wind shear and higher (lower) ambient humidity, both of which help promote (prohibit) the growth of frozen particles and the organization of convective systems. An ensemble of squall line cases is studied at the end to demonstrate that the PD asymmetry in the leading and trailing edges of the deep convection line is closely tied to the anvil cloud and stratiform precipitation layers respectively, suggesting the potential usefulness of PD as a proxy of stratiform/convective precipitation flag, as well as a proxy of convection life stage.

Ice crystal microphysics↗

Raindrop Signature from Microwave Radiometer over Deserts

Rainfall estimates from spaceborne microwave radiometers form the foundation of global precipitation datasets. Since the beginning of the satellite microwave rainfall estimation era in the1980s, the primary signature leveraged over land for these estimates has been the brightness temperature (TB) depression due to ice particle scattering. Contrary to this practice, time series analyses based on observations from two spaceborne radars and two spaceborne radiometers reveal a TB increase due to raindrop emission as the primary cloud particle signature over desert terrain. Low surface emissivity supports the use of liquid raindrop emission as the primary signature over desert surfaces. In these regions, the surface rain rate better correlates with the liquid raindrop emission signal than with the scattering induced by ice further aloft, suggesting a new potential for improving rainfall estimation over deserts by exploiting the liquid raindrop emission signature.

Yalei You↗

Reconciling CloudSat and GPM Estimates of Falling Snow

Satellite-based estimates of falling snow have been provided by CloudSat (launched in 2006) and the Global Precipitation Measurement (GPM) core satellite (launched in 2014). The CloudSat estimates are derived from W-band radar measurements whereas the GPM estimates are derived from its scanning Ku- and Ka-band Dual-Frequency Precipitation Radar (DPR) and 13-channel microwave imager (GMI). Each platform has advantages and disadvantages: CloudSat has higher resolution (approximately 1.5 km) and much better sensitivity (-28 dBZ), but poorer sampling (nadir-only and daytime-only since 2011) and the reflectivity-snowfall (Z-S) relationship is poorly constrained with single-frequency measurements. Meanwhile, DPR suffers from relatively poor resolution (5 km) and sensitivity (approximately 13 dBZ), but has cross-track scanning capability to cover a 245-km swath. Additionally, where Ku and Ka measurements are available, the conversion of reflectivity to snowfall rate is better-constrained than with a single frequency.

S Joseph Munchak↗

Evaluation of Precipitation Detection over Various Surfaces from Passive Microwave Imagers and Sounders

During the middle part of this decade a wide variety of passive microwave imagers and sounders will be unified in the Global Precipitation Measurement (GPM) mission to provide a common basis for frequent (3 hr), global precipitation monitoring. The ability of these sensors to detect precipitation by discerning it from non-precipitating background depends upon the channels available and characteristics of the surface and atmosphere. This study quantifies the minimum detectable precipitation rate and fraction of precipitation detected for four representative instruments (TMI, GMI, AMSU-A, and AMSU-B) that will be part of the GPM constellation. Observations for these instruments were constructed from equivalent channels on the SSMIS instrument on DMSP satellites F16 and F17 and matched to precipitation data from NOAA's National Mosaic and QPE (NMQ) during 2009 over the continuous United States. A variational optimal estimation retrieval of non-precipitation surface and atmosphere parameters was used to determine the consistency between the observed brightness temperatures and these parameters, with high cost function values shown to be related to precipitation. The minimum detectable precipitation rate, defined as the lowest rate for which probability of detection exceeds 50%, and the detected fraction of precipitation, are reported for each sensor, surface type (ocean, coast, bare land, snow cover) and precipitation type (rain, mix, snow). The best sensors over ocean and bare land were GMI (0.22 mm/hr minimum threshold and 90% of precipitation detected) and AMSU (0.26 mm/hr minimum threshold and 81% of precipitation detected), respectively. Over coasts (0.74 mm/hr threshold and 12% detected) and snow-covered surfaces (0.44 mm/hr threshold and 23% detected), AMSU again performed best but with much lower detection skill, whereas TMI had no skill over these surfaces. The sounders (particularly over water) benefited from the use of re-analysis data (vs. climatology) to set the a-priori atmospheric state and all instruments benefit from the use of a conditional snow cover emissivity database over land. It is recommended that real-time sources of these data be used in the operational GPM precipitation algorithms.

GPM↗

Detection Thresholds of Falling Snow from Satellite-Borne Active and Passive Sensors

Precipitation, including rain and snow, is a critical part of the Earth's energy and hydrology cycles. Precipitation impacts latent heating profiles locally while global circulation patterns distribute precipitation and energy from the equator to the poles. For the hydrological cycle, falling snow is a primary contributor in northern latitudes during the winter seasons. Falling snow is the source of snow pack accumulations that provide fresh water resources for many communities in the world. Furthermore, falling snow impacts society by causing transportation disruptions during severe snow events. In order to collect information on the complete global precipitation cycle, both liquid and frozen precipitation must be collected. The challenges of estimating falling snow from space still exist though progress is being made. These challenges include weak falling snow signatures with respect to background (surface, water vapor) signatures for passive sensors over land surfaces, unknowns about the spherical and non-spherical shapes of the snowflakes, their particle size distributions (PSDs) and how the assumptions about the unknowns impact observed brightness temperatures or radar reflectivities, differences in near surface snowfall and total column snow amounts, and limited ground truth to validate against. While these challenges remain, knowledge of their impact on expected retrieval results is an important key for understanding falling snow retrieval estimations. Since falling snow from space is the next precipitation measurement challenge from space, information must be determined in order to guide retrieval algorithm development for these current and future missions. This information includes thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types. For example, can a lake effect snow system with low (approx 2.5 km) cloud tops having an ice water content (IWC) at the surface of 0.25 g / cubic m and dendrite snowflakes be detected? If this information is known, we can focus retrieval efforts on detectable storms and concentrate advances on achievable results. Here, the focus is to determine thresholds of detection for falling snow for various snow conditions over land and lake surfaces. The results rely on simulated Weather Research Forecasting (WRF) simulations of falling snow cases since simulations provide all the information to determine the measurements from space and the ground truth. Sensitivity analyses were performed to better ascertain the relationships between multifrequency microwave and millimeter-wave sensor observations and the falling snow/underlying field of view. In addition, thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types were studied. Results will be presented for active radar at Ku, Ka, and W-band and for passive radiometer channels from 10 to 183 GHz.

Gail Skofronick-Jackson↗

Reconfiguring COSSIR for the Next Generation of Cloud and Precipitation Science

NASA originally developed the Compact Scanning Sub-millimeter Imaging Radiometer (CoSSIR) in the early 2000s, where it was own in the CRYSTAL-FACE and TC4 field experiments. With renewed interest in sub-millimeter remote sensing motivated by the upcoming launch of the EUMETSAT-sponsored Ice Cloud Imager (ICI), along with its potential for inclusion as part of the NASA Decadal Survey Aerosols, CLouds, Convection, and Precipitation (ACCP) observing architecture, CoSSIR is being refurbished to prepare for routine suborbital submillimeter measurements this decade. Upgrades include hardware (procurement of shared parts so that CoSMIR and CoSSIR can fly concurrently) and a change to the channels, with a focus on differential scattering and polarization at eight frequencies ranging from 170 to 640 GHz. As part of the refurbishment, the instrument name has been changed to Configurable Scanning Sub-millimeter Instrument/Radiometer, retaining the same acronym, in order to emphasize the science objectives that are enabled by CoSSIR's unique scanning geometry and channel set. These anticipated capabilities include two-dimensional pro ling of clouds and precipitation and characterization of the geometry (shape and orientation) of large cloud and precipitating ice particles.

Radiometry↗

Polarized Observations of Falling Snow from IMPACTS

The polarization of electromagnetic waves (radio waves and light) provides useful information on the dynamics and structure of precipitating storms. Polarization results from aligned precipitation particles interacting with the electromagnetic waves. This orientation is related to the dynamics of the weather system, and the orientation effects the radiative properties of the clouds (which is related to climate sensivity and feedbacks). We use observations from the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Storms (IMPACTS), a NASA field campaign aiming to understand winter storms that regularly impact the East Coast of the United States, to understand the relationship between properties of snowstorms and polarization of millimeter waves (89 GHz and 166 GHz). We see that height and intensity are strong indicators of the polarization signatures observed in these systems.

Ian S Adams↗

Assessment of Retrieved GMI Emissivity Over Land, Snow and Sea Ice in the GEOS System

Measurements from microwave sounders and imagers provide a valuable source of information including atmospheric temperature and water vapor in Numerical Weather Prediction (NWP) systems that assimilate these observations directly over water surfaces (oceans and other large water bodies). In a recent decadal survey, targeted observables in the Planetary Boundary Layer (PBL) were cited as a key need for future observations (NASEM, 2018). Microwave observations which sense in the PBL are currently available, however, utilizing surface-sensitive microwave observations for atmospheric data assimilation remains a challenge over land, snow and sea ice. This is in part due to the inability of surface emissivity models used by NWP data assimilation systems to simulate observations with sufficient accuracy. The GEOS-ADAS (Todling and el Akkraoui, 2018) which utilizes the Community Radiative Transfer Model (CRTM) (Han, 2006; Chen 2009) is no exception. The ECMWF system has retrieved instantaneous surface emissivity from surface-sensitive channels for SSMI/S and MHS radiance observations, and apply these estimates to the closest channels higher in frequency (Baordo and Geer 2016) in the calculation of simulated radiances. This approach currently is also being tested in the GEOS-ADAS for AMSU-A and ATMS radiances (Zhu et al. 2021). No or minimal emissivity spectral variability has been assumed in the above-mentioned studies. Recently, work by Munchak et al., 2020 (hereby referred to as M2020) provided a new database for emissivity over land, snow and sea ice retrieved from the NASA Global Precipitation Mission (GPM). Compared with Tool to Estimate Land Surface Emissivities at Microwave (TELSEM2; Wang et al., 2017), M2020 provides emissivities for more frequencies(i.e., 10.7 GHz V/H). Moreover, this database is unique in that it utilizes both active and passive data to retrieve surface emissivity and normalized radar cross section. While the emissivity values may be useful for other sensors, they are most applicable to the GPM Microwave Imager (GMI). In this work the GEOS-ADAS is modified to utilize emissivity values from Munchak et al, 2020 in place of values used by CRTM. Presently, only GMI radiances over ocean are used in the operational GEOS-ADAS. This study will focus on the GMI radiances over land, snow, and ice, as a first attempt to evaluate GMI radiances over these non-water surface types. Two cases are then presented, one with one week of observation minus background departures using the modified GEOS-ADAS, and one utilizing the original GEOS-ADAS. It should be noted that the surface emissivity models in CRTM are not state of the art and are scheduled to be replaced by the Community Surface Emissivity Module (CSEM; Chen and Weng, 2016). Simulations using default CRTM emissivity values are used merely as reference comparing against M2020, and is not a thorough comparison against other more state of the art modules such as CSEM.

GMI↗