Polarimetric Radar Verification of GPM Satellite-Based Retrievals of the Raindrop Size Distribution
No abstract available
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
No abstract available
No abstract available
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The use of vertical-incidence Doppler velocity in addition to radar reflectivity may yield information on drop size distribution and therefore result in better rainrate estimates. Doppler velocity can provide useful information on the raindrop size distribution. Doppler velocities from a zenith-pointing radar represent the sum of the mean reflectivity-weighted hydrometeor fallspeed and the vertical air motion. Dual-parameter rain estimation methods using the Doppler velocity, require that the latter can be removed, or is negligible. Atlas et al. (1972) derived relations between Doppler velocity, reflectivity, and rain rate assuming an exponential size distribution for rain. Ulbrich (1994) expanded on this work by deriving the relation between the Doppler velocity and the reflectivity assuming a Gamma size distribution. This distribution provides a more realistic representation of the small rain drops. To get accurate information on raindrop size distributions with the above method, the air motions must be removed from the observed Doppler velocities
The Precipitation Radar aboard the Tropical Rain Measuring Mission (TRMM) Satellite has shown the potential for spaceborne sensing of snow and rain by means of an incoherent pulsed radar operating at 13.8 GHz. The primary advantage of radar relative to passive instruments arises from the fact that the radar can image the 3-dimensional structure of storms. As a consequence, the radar data can be used to determine the vertical rain structure, rain type (convective/stratiform) effective storm height, and location of the melting layer. The radar, moreover, can be used to detect snow and improve the estimation of rain rate over land. To move toward spaceborne weather radars that can be deployed routinely as part of an instrument set consisting of passive and active sensors will require the development of less expensive, lighter-weight radars that consume less power. At the same time, the addition of a second frequency and an upgrade to Doppler capability are features that are needed to retrieve information on the characteristics of the drop size distribution, vertical air motion and storm dynamics. One approach to the problem is to use a single broad-band transmitter-receiver and antenna where two narrow-band frequencies are spaced apart by 5% to 10% of the center frequency. Use of Ka-band frequencies (26.5 GHz - 40 GHz) affords two advantages: adequate spatial resolution can be attained with a relatively small antenna and the differential reflectivity and mean Doppler signals are directly related to the median mass diameter of the snow and raindrop size distributions. The differential mean Doppler signal has the additional property that this quantity depends only on that part of the radial speed of the hydrometeors that is drop-size dependent. In principle, the mean and differential mean Doppler from a near-nadir viewing radar can be used to retrieve vertical air motion as well as the total mean radial velocity. In the paper, we present theoretical calculations for the differential reflectivity and Doppler as functions of the center frequency, frequency difference, and median mass diameter. For a fixed pair of frequencies, the detectability of the differential signals can be expressed as the number of independent samples required to detect rain or snow with a particular median mass diameter. Because sampling numbers on the order of 1000 are needed to detect the differential signal over a range of size distributions, the instrument must be confined to a near-nadir, narrow swath. Radar measurements from a zenith directed radar operated at 9.1 GHz and 10 GHz are used to investigate the qualitative characteristics of the differential signals. Disdrometer and rain gauge data taken at the surface, just below the radar, are used to test whether the differential signals can be used to estimate characteristics of the raindrop size distribution.
The Goddard Cumulus Ensemble (GCE) model with bin spectra microphysics is used to simulate mesoscale convective systems.The model uses explicit bins to represent size spectra of cloud nuclei, water drops, ice crystals, snow and graupel. Each hydrometeorite category is described by 33 mass bins. The simulations provide a unique data set of simulated raindrop size distribution in a realistic dynamic frame. Calculations of radar parameters using simulated drop size distribution serve as an evaluation of numerical model performance. In addition, the GCE bin spectra modes is a very useful tool to study uncertainties related to radar observations; all the environmental parameters are precisely known. In this presentation, we concentrate on the discussion of Z-R (ZDR-R) relation in the simulated systems. Due to computational limitations, the spectra bin model has been run in two dimensions with 31 stretched vertical layers and 1026 horizontal grid points (1 km resolution). Two different cases, one in midlatitude continent, the other in tropical ocean, have been simulated. The continental case is a strong convection which lasted for two hours. The oceanic case is a persistent system with more than 10 hours' life span. It is shown that the simulated Z-R (ZDR-R) relations generally agree with observations using radar and rain gauge data. The spatial and temporal variations of Z-R relation in different locations are also analyzed. Impact of aerosols on cloud formation and raindrop size distribution was studied. Both clean (low CCN) and dirty (high CCN) cases are simulated. The Z-R relation is shown to vary considerable in the initial CCN concentrations.
Raindrop size distributions (RSDs) are often estimated using surface raindrop sampling devices (e.g., disdrometers) or optical array (2D-PMS) probes. A number of authors have used these measured distributions to compute certain higher-order RSD moments that correspond to radar reflectivity, attenuation, optical extinction, etc. Scatter plots of these RSD moments versus disdrometer-measured rainrates are then used to deduce physical relationships between radar reflectivity, attenuation, etc., which are measured by independent instruments (e.g., radar), and rainrate. In this paper RSDs of the gamma form as well as radar reflectivity (via time series simulation) are simulated to study the correlation structure of radar estimates versus rainrate as opposed to RSD moment estimates versus rainrate. The parameters N0, D0 and m of a gamma distribution are varied over the range normally found in rainfall, as well as varying the device sampling volume. The simulations are used to explain some possible features related to discrepancies which can arise when radar rainfall measurements are compared with surface or aircraft-based sampling devices.
A comparative study of raindrop size distribution measurements has been conducted at NASA's Goddard Space Flight Center where the focus was to evaluate the performance of the upgraded laser-optical OTT Particle Size Velocity (Parsivel2; P2) disdrometer. The experimental setup included a collocated pair of tipping-bucket rain gauges, OTT Parsivel (P1) and P2 disdrometers, and Joss-Waldvogel (JW) disdrometers. Excellent agreement between the two collocated rain gauges enabled their use as a relative reference for event rain totals. A comparison of event total showed that the P2 had a 6%absolute bias with respect to the reference gauges, considerably lower than the P1 and JW disdrometers. Good agreement was also evident between the JW and P2 in hourly raindrop spectra for drop diameters between 0.5 and 4 mm. The P2 drop concentrations mostly increased toward small sizes, and the peak concentrations were mostly observed in the first three measurable size bins. The P1, on the other hand, underestimated small drops and overestimated the large drops, particularly in heavy rain rates. From the analysis performed, it appears that the P2 is an improvement over the P1 model for both drop size and rainfall measurements. P2 mean fall velocities follow accepted terminal fall speed relationships at drop sizes less than 1 mm. As a caveat, the P2 had approximately 1ms21 slower mean fall speed with respect to the terminal fall speed near 1 mm, and the difference between the mean measured and terminal fall speeds reduced with increasing drop size. This caveat was recognized as a software bug by the manufacturer and is currently being investigated.
Consistent with the original proposal and work plan, this project focused on estimating the raindrop size distributions (DSDs) retrieved from vertically pointing Doppler radar profilers and analyzing the relationship of the retrieved DSDs with the dynamics of the precipitation processes. The first phase of this project focused on developing the model to retrieve the DSD from the observed Doppler velocity spectra. The second phase used this model to perform DSD retrievals from the profiler observations made during the TRMM Ground Validation Field Campaigns of TEFLUN-B, TRMM-LBA, and KWAJEX. The third phase of this project established collaborations with scientists involved with each field campaign in order to validate the profiler DSD estimates and to enable the profiler retrievals to be used in their research. Through these collaborations, the retrieved DSDs were placed into context with the dynamical processes of the observed precipitating cloud systems.
Disdrometers are useful instruments that measure the distribution of raindrop sizes (drop size distribution: DSD) and the associated rainfall rates/accumulation as those raindrops fall to the ground. These DSDs are a quantity of interest to members of the modeler and observational communities. However, use of disdrometer data for model evaluation, radar monitoring, or other activities requires careful quality control and processing for key DSD properties of interest (e.g., the number concentration of drops) to ensure appropriate physical (scattering, fall speed) assumptions. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Laser Disdrometer Quantities and Video Disdrometer Quantities Value-Added Product (LDQUANTS/VDISQUANTS VAP) uses standard methods from Tokay et al. (e.g., 2013, 2014) to filter drops with unrealistic fall speeds. Further, it estimates several microphysical/geophysical quantities of parameterized DSDs (gamma or exponential assumption type fitting methods) as in previous disdrometer studies and ARM long-term efforts (e.g., Testud et al. 2001, Giangrande et al. 2014, Thompson et al. 2015). Disdrometers are also beneficial for cross-checks with other instrumentation, including rain gauges and radars. To support research interests and related radar monitoring activities, this product calculates radar equivalent quantities, including dual polarization radar quantities (e.g., Z, Differential Reflectivity ZDR, etc.), using T-Matrix scattering and additional wavelength, temperature, and drop shape assumptions (e.g., Thurai et al. 2007). All of these efforts allow disdrometer data sets to become more useful and easily handled by modeling and observational studies, or routine instrumentation checks.
Characteristics of the melting layer and raindrop size distribution can be exploited to further improve radar quantitative precipitation estimation (QPE). Using dual-polarimetric radar and disdrometers, we found that the characteristic size of raindrops reaching the ground in stratiform precipitation often varies linearly with the depth of the melting layer. As a result, a radar rainfall estimator was formulated using D(sub m) that can be employed by polarimetric as well as dual-frequency radars (e.g., space-based radars such as the GPM DPR), to lower the bias and uncertainty of conventional single radar parameter rainfall estimates by as much as 20%. Polarimetric radar also suffers from issues associated with sampling the vertical distribution of precipitation. Hence, we characterized the vertical profile of polarimetric parameters (VP3)-a radar manifestation of the evolving size and shape of hydrometeors as they fall to the ground-on dual-polarimetric rainfall estimation. The VP3 revealed that the profile of ZDR in stratiform rainfall can bias dual-polarimetric rainfall estimators by as much as 50%, even after correction for the vertical profile of reflectivity (VPR). The VP3 correction technique that we developed can improve operational dual-polarimetric rainfall estimates by 13% beyond that offered by a VPR correction alone.
Yearly, monthly, and time of day fade statistics are presented and characterized. A 19.04 GHz yearly fade distribution, corresponding to a second COMSTAR beacon frequency, is predicted using the concept of effective path length, disdrometer, and rain rate results. The yearly attenuation and rain rate distributions follow with good approximation log normal variations for most fade and rain rate levels. Attenuations were exceeded for the longest and shortest periods of times for all fades in August and February, respectively. The eight hour time period showing the maximum and minimum number of minutes over the year for which fades exceeded 12 db were approximately between 1600 to 2400, and 0400 to 1200 hours, respectively. In employing the predictive method for obtaining the 19.04 GHz fade distribution, it is demonstrated theoretically that the ratio of attenuations at two frequencies is minimally dependent of raindrop size distribution providing these frequencies are not widely separated.
In this work, an intercomparison of raindrop mean diameter frequency distribution (RDFD) is performed for numerical simulations of precipitating cloud systems using an array of models and microphysics schemes. This includes results from the Regional Atmospheric Modeling System (RAMS) double-moment microphysics, the Hebrew University Cloud Model bin microphysics (HUCM) interfaced to the RAMS parent model, and the Weather Research and Forecasting (WRF) Model with the Thompson, Morrison, double-moment 6-class (WDM6), and National Severe Storms Laboratory (NSSL) double-moment schemes. Simulations are examined with respect to the raindrop size distribution (DSD) volume-number mean diameter (D m ) and intercept parameter (N w ). When compared to a suite of disdrometer observations, the RDFD resulting from each microphysics scheme exhibits varying degrees of mean drop size constraints and peaks in the frequency distribution of D m . A more detailed investigation of the peaked RDFD from the RAMS simulations suggests that the parameterization of raindrop collisional breakup can impose strong limitations on the evolution of simulated drop growth. As such, a summary and comparison of the drop breakup parameterizations among the aforementioned microphysics schemes is presented. While some drop breakup parameterizations are adjusted toward the observations by modifying the threshold diameter for the onset of breakup, this study explores the use of a modified maximum breakup efficiency. This method permits the parameterization to retain its threshold breakup diameter, while limiting the strength of drop breakup and permitting a broader range of drop sizes. As a result, the simulated mean drop sizes are in better agreement with observations.
Rain attenuation is strongly dependent on the rain rate, but also on the rain drop size distribution (DSD). Typically, models utilize an average drop size distribution, such as those developed by Laws and Parsons, or Marshall and Palmer. However, individual rain events may possess drop size distributions which could be significantly different from the average and will impact, for example, fade mitigation techniques which utilize channel performance estimates from a signal at a different frequency. Therefore, a good understanding of the characteristics and variability of the raindrop size distribution is extremely important in predicting rain attenuation and instantaneous frequency scaling parameters on an event-toevent basis. Since June 2014, NASA Glenn Research Center (GRC) and the Politecnico di Milano (POLIMI) have measured the attenuation due to rain in Milan, Italy, on the 20/40 GHz beacon signal broadcast from the Alphasat TDP#5 Aldo Paraboni Q/V-band Payload. Concomitant with these measurements are the measurements of drop size distribution and rain rate utilizing a Thies Clima laser precipitation monitor (disdrometer). In this paper, we discuss the comparison of the predicted rain attenuation at 20 and 40 GHz derived from the drop size distribution data with the measured rain attenuation. The results are compared on statistical and real-time bases. We will investigate the performance of the rain attenuation model, instantaneous frequency scaling, and the distribution of the scaling factor. Further, seasonal rain characteristics will be analysed.
A cloud-resolving model is used to study sensitivities of two different microphysical schemes, one is the traditional bulk type, and the other is an explicit bin scheme, in simulating a mid-latitude squall line case (PRE-STORM, June 10-1 1,1985). Simulations using different microphysical schemes are compared with each other and also with the observations. Both the bulk and bin models reproduce the general features during the developing and mature stage of the system. Furthermore, the observations and the well-proven bulk scheme simulation serve as validations for the newly incorporated bin scheme. However, it is also shown that the bulk and bin simulations have distinct differences, most notably in the stratiform region of the squall line system. Weak convective cells exist in the stratiform region in the bulk simulation, but not in the bin simulation. These weak convective cells in the stratiform region simulated in the bulk scheme model are remnants of the stronger convections previously at the leading edge of the system, sustained by horizontal vorticity generated by its own cool pool near the surface. The bin simulation, on the other hand, has a horizontally homogeneous stratiform cloud structure, which agrees better with the observations. Examinations of the downdraft core strength, the potential temperature perturbation, and the evaporative cooling rate show that the differences between the bulk and bin models are due mainly to the stronger low-level evaporative cooling in the convective zone simulated in the bulk microphysical scheme, which is unrealistic because of the assumptions made in raindrop size distribution. Further sensitivity tests that reduce the evaporation rate in bulk scheme artificially produce more upright convective core and less weak cores in stratiform region. However, they produce weaker upper level outflow and consequently less stratiform rain area. The addition of a more realistic raindrop breakup scheme in the bin scheme results more realistic radar reflectivity and stronger surface rainfall. Despite the increase of the rain evaporation and strengthening of the near surface cool pool, bin scheme with rain breakup shows homogeneous stratiform rain. These sensitivity tests prove the robustness of the bin microphysical scheme and the difficulty of tuning the limited parameters in the bulk microphysical scheme to realistically reproduce detail structures in a mid-latitude squall line case study.
Airborne meteorological radars typically operate at attenuating wavelengths. The path integrated attenuation (PIA) can be estimated using the surface reference technique (SRT). In this method, an initial value is determined for the radar cross section of the earth surface in a rain-free area in relatively close proximity to the rain cloud. During subsequent observations of precipitation any decrease 'in the observed surface cross section from the reference value s assumed to be a result of the two-way attenuation along the propagation path. In this paper we present selected instances of high PIA observed over land by an airborne radar. The observations were taken in Brazil and Florida during TRMM (Tropical Rainfall Measurement Mission) field campaigns. We compared these observations with collocated and nearly simultaneous ground-based radar observations by an S-band radar that is not subject to significant attenuation. In this preliminary evaluation, a systematic difference in the attenuation in the two storms is attributed to a difference in the raindrop size distributions; this is supported by observations of ZDR (differential reflectivity).
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.