A Field Study of Footprint-Scale Variability of Raindrop Size Distribution
Explore the source record for details and available documents.
Ali Tokay↗
SEARCH · Search NASA
Results for “raindrop size distribution”
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.
Rain Fallspeeds and Rates Derived from Airborne Nadir-Pointing Doppler Radar Measurements
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
Heymsfield, Gerald M.↗
Airborne Differential Doppler Weather Radar
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.
Meneghini, R.↗
Radar Reflectivity Simulated by a 2-D Spectra Bin Model: Sensitivity of Cloud-aerosol Interaction
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.
Li, Kiaowen↗
Simulation of radar reflectivity and surface measurements of rainfall
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.
Chandrasekar, V.↗
Evaluation of the New Version of the Laser-Optical Disdrometer, OTT Parsivel2
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.
disdrometers↗
Retrieval of Hydrometeor Drop Size Distributions from TRMM Field Campaign Profiler Doppler Velocity Spectra Observations
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.
Williams, Christopher R.↗
Utilizing the Vertical Variability of Precipitation to Improve Radar QPE
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.
Precipitation↗
Cumulative slant path rain attenuation associated with COMSTAR beacon at 28.56 GHz for Wallops Island, Virginia
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.
Goldhirsh, J.↗
Comparison of Instantaneous Frequency Scaling from Rain Attenuation and Optical Disdrometer Measurements at K/Q bands
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.
millimeter wave↗
Sensitivity of a Cloud-Resolving Model to the Bulk and Explicit Bin Microphysical Schemes
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.
Li, Xiao-Wen↗
Radar Observation of Large Attenuation in Convective Storms: Implications for the Dropsize Distribution
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).
Tian, Lin↗
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↗
To What Extent Can Raindrop Size Be Determined by a Multiple-Frequency Radar?
In this paper, an analytical treatment of the atmospheric remote sensing problem of determining the raindrop size distribution (DSD) with a spaceborne multifrequency microwave nadir-looking radar system is presented. It is typically assumed that with two radar measurements at different frequencies one ought to be able to calculate two state variables of the DSD: a bulk quantity, such as the rain rate, and a distribution shape parameter. To determine if this nonlinear problem can indeed be solved, the DSD is modeled as a Gamma distribution and quadratic approximations to the corresponding radar-rain relations are used to examine the invertibility of the resulting system of equations in the case of two as well as three radar frequencies. From the investigation, it is found that for regions of DSD state space multiple solutions exist for two or even three different frequency radar measurements. This should not be surprising given the nonlinear coupled nature of the problem.
rain characteristics↗
A standard verification for rainfall estimation from remote platforms
Measurements of reflectivity at horizontal (Zh) and vertical (Zy) polarizations provide adequate information necessary to infer the two parameters of an exponential raindrop size distribution (No, Do) where the distribution is given by N(D) = Noexp(-3.67 D/Do). This distribution enables computation of water content or still air rainfall rates. The physical basis of the radar technique is outlined and illustrated theoretically, and experimental results, comparing radar derived rainfall rates with raingage and disdrometer measurements, are reviewed. The technique is useful for many meteorological and hydrological purposes, including ground truth measurements of rainfall rate over the ocean for comparison with satellite related observations.
Meitin, J. G., Jr.↗
On the measurement of precipitation frequencies by passive microwave radiometery
Over most of the microwave spectrum, raindrops both absorb and scatter radiation producing large changes in brightness temperatures relative to clear or cloudy conditions. Since the structure of rain varies substantially for different rain rates and climatological backgrounds, the raindrop size distribution, the rain layer thickness and the ice clouds above the rain layer are all important inputs to the model computations. The subsequent modeling involves applying the Mie theory to derive the absorption and scattering effects and the radiative transfer calculation is based upon a variational iterative approach which takes account of the multiple scattering effect of the rain layer. Results over both ocean and land backgrounds are demonstrated. It is also demonstrated that by using discrimination tests of the radiometric data, the rain/no rain decision can be made and the rainfall rate can be retrieved from a statistical inversion technique.
Kidder, S. Q.↗
Passive and active remote sensing of atmospheric precipitation
Both passive and active remote sensing of atmospheric precipitation are studied with the vector radiative transfer equations by making use of the Mie scattering phase functions and incorporating the raindrop-size distributions. For passive remote sensing, the Gaussian quadrature method is employed to solve for the brightness temperatures. For active remote sensing, an iterative approach carrying out to the second order in albedo is used to calculate for the bistatic coefficients, the backscattering cross sections/unit volume, and the interchannel cross talks. The calculated results are plotted as a function of rainfall rates and compared to various available experimental data. The theoretical model is easily applied to the remote sensing of aerosol particles, smoke, fog, and haze at infrared and visible frequencies.
Jin, Y. Q.↗
Extinction of visible and infrared radiation in rain Comparison of theory and experiment
A critical review is given of the experimental and theoretical results concerning the measurement of rainfall using optical extinction, i.e., the attenuation of radiation with wavelength less than or equal to that of the infrared band. It is shown that rainfall rates found from an empirical relation involving optical extinction generally display average deviations without regard for sign of only 25 percent when compared with those measured by raingages directly beneath the optical beam. It is also shown that the differences between experimental and theoretical results can be explained in terms of variations of the shape of the raindrop size distribution, i.e., deviations from exponentiality.
Ulbrich, C. W.↗