Search NASASearch

Engineering topics

Ali Tokay

Publications and source records attributed to Ali Tokay.

At least 19 records

Retrieval of Normalized Gamma Size Distribution Parameters Using Precipitation Imaging Package (Pip) Snowfall Observations During Ice-Pop 2018

Parameters of the normalized gamma particle size distribution (PSD) have been retrieved from the Precipitation Image Package (PIP) snowfall observations collected during the International Collaborative Experiment - PyeongChang Olympics and Paralympic (ICE-POP 2018). Two of the gamma PSD parameters, the mass weighted particle diameter (Dmass) and the normalized intercept parameter NW, have median values of 1.15-1.31 mm and 2.84-3.04 log(mm-1 m-3), respectively. This range arises from the choice of the relationship between the maximum versus equivalent diameter, Dmx-Deq, and the relationship between the Reynolds and Best numbers, Re-X. Normalization of snow water equivalent rate (SWER) and ice water content (W) by NW reduces the range in NW resulting in well fitted power law relationship, between SWER/NW and Dmass and between W/NW and Dmass. The bulk descriptors of snowfall are calculated from PIP observations and from the gamma PSD with values of the shape parameter (μ) ranging from -2 to 10. NASA’s Global Precipitation Measurement (GPM) mission, which adopted the normalized gamma PSD, assumes μ = 2 and μ = 3 in its two separate algorithms. The mean fractional bias (MFB) of the snowfall parameters changes with μ, where the functional dependence on μ depends on the specific snowfall parameter of interest. The MFB of the total concentration was underestimated by 0.23-0.34 when μ = 2 and by 0.29-0.40 when μ = 3, while the MFB of SWER had a much narrower range (-0.03 to 0.04) for the same μ values.

Snowfall

Evaluation of SWER(Ze) Relationships by Precipitation Imaging Package (PIP) during ICE-POP 2018

Improving estimation of snow water equivalent rate (SWER) from radar reflectivity (Ze), known as a SWER(Ze) relationship, is a priority for NASA’s Global Precipitation Measurement (GPM) mission ground validation program as it is needed to comprehensively validate spaceborne precipitation retrievals. This study investigates the performance of eight operational and four research-based SWER(Ze) relationships utilizing Precipitation Imaging Probe (PIP) observations from the International Collaborative Experiment for Pyeongchang 2018 Olympic and Paralympic Winter Games (ICE-POP 2018) field campaign. During ICE-POP 2018, there were 10 snow events that are classified by synoptic conditions as either cold low or warm low, and a SWER(Ze) relationship is derived for each event. Additionally, a SWER(Ze) relationship is derived for each synoptic classification by merging all events within each class. Two new types of SWER(Ze) relationships are derived from PIP measurements of bulk density and habit classification. These two physically based SWER(Ze) relationships provided superior estimates of SWER when compared to the operational, event-specific, and synoptic SWER(Ze) relationships. For estimates of the event snow water equivalent total, the event-specific, synoptic, and best-performing operational SWER(Ze) relationships outperformed the physically based SWER(Ze) relationship, although the physically based relationships still performed well. This study recommends using the density or habit-based SWER(Ze) relationships for microphysical studies, whereas the other SWER(Ze) relationships are better suited toward hydrologic application.

Ali Tokay

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

Retrieval of Snow Water Equivalent by the Precipitation Imaging Package (PIP) over Northern Great Lakes

Performance of the Precipitation Imaging Package (PIP) for estimating the snow water equivalent (SWE) is evaluated through a comparative study with the collocated National Oceanic and Atmospheric Administration National Weather Service snow stake field measurements. The PIP together with a vertically pointing radar, a weighing bucket gauge, and a laser-optical disdrometer was deployed at the NWS Marquette, Michigan office building for a long-term field study supported by the National Aeronautics and Space Administration’s Global Precipitation Measurement mission Ground Validation program. The site was also equipped with a weather station. During the 2017-18 winter, PIP functioned nearly uninterrupted at frigid temperatures accumulating 2345.8 mm of geometric snow depth over a total of 499 hours. This long record consists of 30 events, and the PIP-retrieved and snow stake field measured SWE differed less than 15% in every event. For the windy events, the wind corrected stake field record agreed better with the PIP-retrieved SWE. Two of the major events with the longest duration and the highest accumulation are examined in detail. The particle mass with a given diameter was much lower during a shallow, colder, uniform lake-effect event than in the deep, less cold, and variable synoptic event. This study demonstrated that PIP is a robust instrument for operational use, and is reliable for deriving the bulk properties of falling snow.

Ali Tokay

Snowfall in the Northern Great Lakes: Lessons Learned from a Multisensor Observatory

A multisensor snowfall observational suite has been deployed at the Marquette, Michigan, National Weather Service Weather Forecast Office (KMQT) since 2014. Micro Rain Radar (MRR; profiling radar), Precipitation Imaging Package (PIP; snow particle imager), and ancillary ground-based meteorological observations illustrate the unique capabilities of these combined instruments to document radar and concomitant microphysical properties associated with northern Great Lakes snowfall regimes. Lake-effect, lake-orographic, and transition event case studies are presented that illustrate the variety of snowfall events that occur at KMQT. Case studies and multiyear analyses reveal the ubiquity of snowfall produced by shallow events. These shallow snowfall features and their distinctive microphysical fingerprints are often difficult to discern with conventional remote sensing instruments, thus highlighting the scientific and potential operational value of MRR and PIP observations. The importance of near-surface lake-orographic snowfall enhancement processes in extreme snowfall events and regime-dependent snow particle microphysical variability controlled by regime and environmental factors are also highlighted.

Lake effects

Assessment of Ku- and Ka-band Dual-Frequency Radar for Snow Retrieval

Dual-frequency Ku/Ka-band radar retrievals of snow parameters such as liquid-equivalent snowfall rate (R) and mass-weighted diameter (Dm) have two principal errors, namely, the differences between the assumed particle size distribution (PSD) model from the actual PSD and inadequacies in characterizing the single-scattering properties of snowflakes. Regarding the first issue, this study, based on radar simulations from a large amount of observed PSD data, shows that there exist relatively high correlations between the estimated snow parameters and their true values derived directly from the measured PSD. For PSD data with R greater than 0.1 mm h−1, a gamma PSD model with a fixed shape factor (μ) equal to 0 (or exponential distribution) provides the best estimates of R and Dm. This is despite negative biases of up to −15 % in R and underestimates and overestimates in Dm for small and large Dm, respectively. The μ = 0 assumption, however, produces relatively poor estimates of normalized intercepts of a gamma PSD (Nw), whereas the best estimates are obtained when μ is considered either 3 or 6. However, the use of an inappropriate scattering table increases the errors in snow retrieval. Simple evaluations are made for cases where the scattering databases used for the algorithm input differ from that used for retrieval. The mismatched scattering databases alone could cause at least 30 %-50 % changes in the estimates of snow water content (SWC) and R and could affect the retrievals of Dm and Nw and their dependence on μ.

Liang Liao

Characteristics of DSD Bulk Parameters: Implication to Radar Rain Retrieval

With the use of 213,456 one-minute measured data of droplet-size distribution (DSD) of rain collected during several National Aeronautics and Space Administration (NASA)-sponsored field campaigns, the relationships between rainfall rate R, mass-weighted diameter D(m) and normalized intercept parameter N(w) of the gamma DSD are studied. It is found, based on the simulations of the gamma DSD model, that R, D(m) and N(w) are closely interrelated, and that the ratio of R to N(w) is solely a function of D(m), independent of the shape factor μ of the gamma distribution. Furthermore, the model-produced ratio agrees well with those from the DSD data. When a power-law equation is applied to fit the model data, we have: R=aN(w)D(sup b, sub m) , where a=1.588×10(exp −4) , b=4.706 . Analysis of two-parameter relationships such as R–D(m), N(w)–R and N(w)–D(m) reveals that R and D(m) are moderately correlated while N(w) and D(m) are negatively correlated. N(w) and R, however, are uncorrelated. The gamma DSD model also reveals that variation of R–D(m) relation is caused primarily by N(w). For the application of the Ku- and Ka-band dual-frequency radar for the retrieval of the DSD bulk parameters as well as the specific radar attenuations, the study is carried out to relate the dual-frequency radar reflectivity factors to the DSD and attenuation parameters.

Liang Liao

The GPM Ground Validation Program

We present a detailed overview of the structure and activities associated with the NASA-led ground-validation component of the NASA-JAXA Global Pre­cipitation Measurement (GPM) mission. The overarching philosophy and approaches for NASA's GV program are presented with primary focus placed on aspects of direct validation and a summary of physical validation campaigns and results. We describe a spectrum of key instruments, methods, field campaigns and data products developed and used by NASA's GV team to verify GPM level-2 precipitation products in rain and snow. We describe the tools and analysis framework used to confirm that NASA's Level-I science requirements for GPM are met by the GPM Core Observatory. Examples of routine validation activities related to verification of Integrated Multi­satellitE Retrievals for GPM (IMERG) products for two different regions of the globe (Korea and the US) are provided, and a brief analysis related to IMERG performance in the extreme rainfall event associated with Hurricane Florence is discussed.

precipitation

Comparison of Raindrop Size Distribution between NASA’s S-Band Polarimetric Radar and Two-Dimensional Video Disdrometers

The ground-based-radar-derived raindrop size distribution (DSD) parameters—mass-weighted drop diameter Dmass and normalized intercept parameter NW—are the sole resource for direct validation of the National Aeronautics and Space Administration (NASA) Global Precipitation Measurement (GPM) mission Core Observatory satellite-based retrieved DSD. Both Dmass and NW are obtained from radar-measured reflectivity ZH and differential reflectivity ZDR through empirical relationships. This study uses existing relationships that were determined for the GPM ground validation (GV) program and directly compares the NASA S-band polarimetric radar (NPOL) observables of ZH and ZDR and derivedDmass and NW with those calculated by two-dimensional video disdrometer (2DVD). The joint NPOL and 2DVD datasets were acquired during three GPM GV field campaigns conducted in eastern Iowa, southern Appalachia, and western Washington State. The comparative study quantifies the level of agreement for ZH, ZDR, Dmass, and log(NW) at an optimum distance (15–40 km) from the radar as well as at distances greater than 60 km from radar and over mountainous terrain. Interestingly, roughly 10%–15% of the NPOL ZH–ZDR pairs were well outside the envelope of 2DVD-estimated ZH–ZDR pairs. The exclusion of these pairs improved the comparisons noticeably.

Ali Tokay

Comparison of Raindrop Size Distribution between NASA's S-Band Polarimetric Radar and Two-Dimensional Video Disdrometers

The ground-based-radar-derived raindrop size distribution (DSD) parameters—mass-weighted drop diameter D mass and normalized intercept parameter NW—are the sole resource for direct validation of the National Aeronautics and Space Administration (NASA) Global Precipitation Measurement (GPM ) mission Core Observatory satellite-based retrieved DSD. Both D mass and NW are obtained from radar-measured reflectivity Z Hand differential reflectivity ZDR through empirical relationships. This study uses existing relationships that were determined for the GPM ground validation (GV) program and directly compares the NASA S-band polarimetric radar (NPOL) observables of Z Hand Z DR and derived D mass and NW with those calculated by two-dimensional video disdrometer (2DVD). The joint NPOL and 2DVD datasets were acquired during three GPM GV field campaigns conducted in eastern Iowa, southern Appalachia, and western Washington State. The comparative study quantifies the level of agreement for ZH, ZDR, D mass, and log (NW) at an optimum distance (15–40 km) from the radar as well as at distances greater than 60 km from radar and over mountainous terrain. Interestingly, roughly 10%–15% of the NPOLZH–Z DR pairs were well outside the envelope of 2DVD-estimatedZH–Z DR pairs. The exclusion of these pairs improved the comparisons noticeably.

Ali Tokay

Development and Evaluation of the Raindrop Size Distribution Parameters for the NASA Global Precipitation Measurement Mission Ground Validation Program

The National Aeronautics and Space Administration Global Precipitation Measurement (GPM) mission ground validation program uses dual-polarization radar moments to estimate raindrop size distribution (DSD) parameters, the mass-weighted mean drop diameter Dmass, and normalized intercept parameter NW, to validate the GPM Core Observatory–derived DSD parameters. The disdrometer-based Dmass and NW are derived through empirical relationships between Dmass and differential reflectivity ZDR, and between NW, reflectivity ZH, and Dmass. This study employs large datasets collected from two-dimensional video disdrometers (2DVD) during six different field studies to derive the requisite empirical relationships. The uncertainty of the derived Dmass(ZDR) relationship is evaluated through comparisons of 2DVD-calculated and ZDR-estimated Dmass, where ZDR is calculated directly from 2DVD observations. Similarly, the uncertainty of the NW(ZH, Dmass) relationship is evaluated through 2DVD-calculated and Dmass and ZH-estimated NW, where Dmass and ZH are directly calculated from 2DVD observations. This study also presents the sensitivity of Dmass(ZDR) relationships to climate regime and to disdrometer type after developing three additional Dmass(ZDR) relationships from second-generation Particle Size Velocity (PARSIVEL2) disdrometer (P2) observations collected in the Pacific Northwest, in Iowa, and at Kwajalein Atoll in the tropical Pacific Ocean. The application of P2-derived Dmass(ZDR) relationship based on precipitation in the northwestern United States to P2 observations collected over the tropical ocean resulted in the highest error among comparisons of the three datasets.

Ali Tokay

Overview of the New Version 3 NASA Micro-Pulse Lidar Network (MPLNET) Automatic Precipitation Detection Algorithm

Precipitation modifies atmospheric column thermodynamics through the process of evaporation and serves as a proxy for latent heat modulation. For this reason, a correct precipitation parameterization (especially for low-intensity precipitation) within global scale models is crucial. In addition to improving our modeling of the hydrological cycle, this will reduce the associated uncertainty of global climate models in correctly forecasting future scenarios, and will enable the application of mitigation strategies. In this manuscript we present a proof of concept algorithm to automatically detect precipitation from lidar measurements obtained from the National Aeronautics and Space Administration Micropulse lidar network (MPLNET). The algorithm, once tested and validated against other remote sensing instruments, will be operationally implemented into the network to deliver a near real time (latency <1.5 h) rain masking variable that will be publicly available on MPLNET website as part of the new Version 3 data products. The methodology, based on an image processing technique, detects only light precipitation events (defined by intensity and duration) such as light rain, drizzle, and virga. During heavy rain events, the lidar signal is completely extinguished after a few meters in the precipitation or it is unusable because of water accumulated on the receiver optics. Results from the algorithm, in addition to filling a gap in light rain, drizzle, and virga detection by radars, are of particular interest for the scientific community as they help to fully characterize the aerosol cycle, from emission to deposition, as precipitation is a crucial meteorological phenomenon accelerating atmospheric aerosol removal through the scavenging effect. Algorithm results will also help the understanding of long term aerosol–cloud interactions, exploiting the multi-year database from several MPLNET permanent observational sites across the globe. The algorithm is also applicable to other lidar and/or ceilometer network infrastructures in the framework of the Global Aerosol Watch (GAW) aerosol lidar observation network (GALION).

Simone Lolli

Assessing Dual-Polarization Radar Estimates of Extreme Rainfall During Hurricane Harvey

Hurricane Harvey hit the Texas Gulf Coast as a major hurricane on August 25, 2017 before exiting the state as a tropical storm on August 29, 2017. Left in its wake was historic flooding, with some locations measuring more than 60 inches of rain over a five-day period. The WSR-88D radar (KHGX) maintained operations for the entirety of the event. Rain gauge data from the Harris County Flood Warning System (HCFWS) was used for validation with the full radar data set to retrieve daily and event-total precipitation estimates for the period August 25-29, 2017. The KHGX precipitation estimates were then compared to the HCFWS gauges. Three different hybrid polarimetric rainfall retrievals were used, along with attenuation-based retrieval that employs the radar-observed differential propagation. An advantage of using a attenuation-based retrieval is its immunity to partial beam blockage and calibration errors in reflectivity and differential reflectivity. All of the retrievals are susceptible to changes in the observed Drop Size Distribution (DSD). No in situ DSD data were available over the study area, so changes in the DSD were interpreted by examining the observed radar data. We examined the parameter space of two key values in the attenuation retrieval to test the sensitivity of the rain retrieval. Selecting a value of α=0.015 and β=0.600 and β=0.625 provided the best overall results, relative to the gauges, but more work needs to be done to develop an automated technique to account for changes in the ambient DSD.

David B Wolff

Assessing Dual-Polarization Radar Estimates of Extreme Rainfall During Hurricane Harvey

Hurricane Harvey hit the Texas Gulf Coast as a major hurricane on 25 August 2017 before exiting the state as a tropical storm on 29 August 2017. Left in its wake was historic flooding, with some locations measuring more than 60 in. (150 cm) of rain over a 5-day period. The WSR-88D radar (KHGX) maintained operations for the entirety of the event. Rain gauge data from the Harris County Flood Warning System (HCFWS) was used for validation with the full radar dataset to retrieve daily and event-total precipitation estimates for the period 25–29 August2017. The KHGX precipitation estimates were then compared with the HCFWS gauges. Three different hybrid polarimetric rainfall retrievals were used, along with attenuation-based retrieval that employs the radar-observed differential propagation. An advantage of using a attenuation-based retrieval is its immunity to partial beam blockage and calibration errors in reflectivity and differential reflectivity. All of the retrievals are susceptible to changes in the observed drop size distribution (DSD). No in situ DSD data were available over the study area, so changes in the DSD were interpreted by examining the observed radar data. We examined the parameter space of two key values in the attenuation retrieval to test the sensitivity of the rain retrieval. Selecting a value ofa50.015 andb50.600 provided the best overall results, relative to the gauges, but more work needs to be done to develop an automated technique to account for changes in the ambient DSD.

David Billings Wolff

Assessment of Ku- and Ka-Band Dual-Frequency Radar for Retrieval of Snow Properties

Current scattering look-up tables for snow assume a constant mass density along with a Gamma particle size distribution (PSD). The first assumption is tested against scattering parameters from simulated particle models generated by Dr. Kwo-Sen Kuo at GSFC and Dr. Guosheng Liu at FSU. Good agreement of the scattering parameters is found with the FSU results if the mass is taken to be the same as the mass of the simulated particle and the mass density is taken to be 0.2 g/cm cu. For the GSFC data base, good agreement is found if the mass density is taken to be between 0.1 and 0.2 g/cm cu. The second assumption of a Gamma PSD is tested against measured PSD's along with a m-d (mass-dimension) relationship. The degree of agreement depends on the value of, the 'shape' parameter in the Gamma distribution but to a lesser degree on the m-d relationship (of the three that were examined). A shortcoming of the simulated snow particle data bases is the lack of large particles. As a consequence, larger values of DFR (dual-frequency ratio) that are commonly seen in airborne and spaceborne measurements cannot be reproduced from the tables. This situation is expected to improve as scattering parameters from larger particles are included in the databases.

Liang Liao