Search NASA⌕ Search

Engineering topics

Liang Liao

Publications and source records attributed to Liang Liao.

13 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

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↗
NASA NTRS2022

A Generalized Dual-Frequency Ratio (DFR) Approach for Rain Retrievals

The dual-frequency ratio of radar reflectivity factors (DFR) has been shown to be a useful quantity as it is independent of the number concentration of the particle size distribution and primarily a function of the mass-weighted particle diameter, Dm. A drawback of DFR-related methods for rain estimation, however, is the non-unique relationship between Dm and DFR. At Ku- and Ka-band frequencies, two solutions for Dm exist when DFR is less than zero. This ambiguity generates multiple solutions for the range profiles of the particle size parameters. We investigate characteristics of these solutions for both the initial-value (forward) and final-value (backward) forms of the equations. To choose one of among many possible range profiles of Dm, number concentration and rain rate, R, independently measured path attenuations are used. For the backward approach the possibility exists of dispensing with externally measured path attenuations by achieving consistency between the input and output path attenuations. The methods are tested by means of a simulation based on disdrometer-measured raindrop size distributions and the results are compared with a simplified version of the operational R-Dm method.

Robert Meneghini↗
NASA NTRS2022

A Comprehensive Machine Learning Study to Classify Precipitation Type over Land from Global Precipitation Measurement Microwave Imager (GPM-GMI) Measurements

Precipitation type is a key parameter used for better retrieval of precipitation characteristics as well as to understand the cloud–convection–precipitation coupling processes. Ice crystals and water droplets inherently exhibit different characteristics in different precipitation regimes (e.g., convection, stratiform), which reflect on satellite remote sensing measurements that help us distinguish them. The Global Precipitation Measurement (GPM) Core Observatory’s microwave imager (GMI) and dual-frequency precipitation radar (DPR) together provide ample information on global precipitation characteristics. As an active sensor, the DPR provides an accurate precipitation type assignment, while passive sensors such as the GMI are traditionally only used for empirical understanding of precipitation regimes. Using collocated precipitation type flags from the DPR as the “truth”, this paper employs machine learning (ML) models to train and test the predictability and accuracy of using passive GMI-only observations together with ancillary information from a reanalysis and GMI surface emissivity retrieval products. Out of six ML models, four simple ones (support vector machine, neural network, random forest, and gradient boosting) and the 1-D convolutional neural network (CNN) model are identified to produce 90–94% prediction accuracy globally for five types of precipitation (convective, stratiform, mixture, no precipitation, and other precipitation), which is much more robust than previous similar effort. One novelty of this work is to introduce data augmentation (subsampling and bootstrapping) to handle extremely unbalanced samples in each category. A careful evaluation of the impact matrices demonstrates that the polarization difference (PD), brightness temperature (Tc) and surface emissivity at high-frequency channels dominate the decision process, which is consistent with the physical understanding of polarized microwave radiative transfer over different surface types, as well as in snow and liquid clouds with different microphysical properties. Furthermore, the view-angle dependency artifact that the DPR’s precipitation flag bears with does not propagate into the conical-viewing GMI retrievals. This work provides a new and promising way for future physics-based ML retrieval algorithm development.

machine learning/artificial intelligence↗
NASA NTRS2022

GPM DPR Retrievals: Algorithm, Evaluation, and Validation

The primary goal of the Dual-frequency Precipitation Radar (DPR) aboard the Global Precipitation Measurement (GPM) core satellite is to infer precipitationrate and raindrop/particle size distributions (DSD/PSD). The focus of this paper is threefold: 1)description of the DPR retrieval algorithm that uses an adjustable relationship between rainrate (R) and the mass-weighted diameter (Dm) or an R-Dm relationship in solving for R and Dmsimultaneously; 2) evaluation of the DPR algorithm based on the physical simulationsthat employ measured DSD/PSD to understand the mechanism and errorcharacteristics of the retrieval method; 3) review of ground validation studies for theDPR product as well as analysis of the strengths and weaknesses of the ground radarand rain gauge/disdrometer validations. Overall, the DPR Version-6 algorithmprovides reasonably accurate estimates of R and Dm in rain. Non-uniformity in therain profile, however, tends to degrade the accuracy of the R and Dm estimates tosome extent as the range-independent assumption of the adjustable parameter () ofthe R-Dm relation is not able to fully account for natural variation of DSD in the verticalprofile. Underestimation of the DPR snow rate is found when compared with theindependent dual-frequency ratio (DFR) technique. This is possibly the result of theconstraint associated with the path integral attenuation (PIA)/differential PIA (dPIA)used in the DPR algorithm to find the best  and range-independent  assumption. Arange-variable  model, proposed in the DPR Version-7 algorithm, is expected toimprove rain and snow retrieval.

Liang Liao↗
NASA NTRS2020

Path Attenuation Estimates for the GPM Dual-Frequency Precipitation Radar (DPR)

Estimation of path attenuation is a critical part of retrieving precipitation parameters using measurements from the Dual-Frequency Precipitation Radar (DPR) on board the Global Precipitation Measurement Mission (GPM) satellite. In this paper, we describe the latest implementation of the Surface Reference Technique that uses surface scattering properties to infer path attenuation through the precipitation. Both single- and dual-frequency versions of this method are available and while the dual-frequency version appears to be more accurate at moderate rain rates, the single-frequency approach at Ku-band is needed when the Ka-band data are not available. Despite improvements afforded by the dual-frequency version of the method, other methods such as the Hitschfeld-Bordan and standard dual-frequency approaches offer advantages particularly at lighter rain rates and at near-nadir incidence angles over land. Weighted averages of the results from these methods appear to offer the best estimate of path attenuation presently available.

attenuation↗
NASA NTRS2020

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↗
NASA NTRS2020

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↗
NASA NTRS2020

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↗
NASA NTRS2019

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↗
NASA NTRS2019

Physical Evaluation of GPM DPR Single- and Dual-Wavelength Algorithms

A physical evaluation of the rain profiling retrieval algorithms for the Dual-frequency Precipitation Radar (DPR) aboard the Global Precipitation Measurement (GPM) core satellite is carried out by applying them to the hydrometeor profiles generated from measured raindrop size distributions (DSD). The DSD-simulated radar profiles are used as input to the algorithms, and their estimates of hydrometeors' parameters are compared with the same quantities derived directly from the DSD data (or truth). The retrieval accuracy is assessed by the degree to which the estimates agree with the truth. To check the validity and robustness of the retrievals, the profiles are constructed for cases ranging from fully correlated (or uniform) to totally uncorrelated DSDs along the columns. Investigation into the sensitivity of the retrieval results to the model assumptions is made to characterize retrieval uncertainties and identify error sources. Comparisons between the single- and dual-wavelength algorithm performance are carried out with either a single- or dual-wavelength constraint of the path integral or differential path integral attenuation. The results suggest that the DPR dual-wavelength algorithm generally provides accurate range-profiled estimates of rainfall rate and mass-weighted diameter with the dual-wavelength estimates superior in accuracy to those from the single-wavelength retrievals.

Liang Liao↗