## 1. Introduction

With the upgrade of the National Weather Service’s Weather Surveillance Radar-1988 Doppler (WSR-88D) network to dual-polarization capabilities comes additional information that operational meteorologists can use to better diagnose and predict the local weather (e.g., Zrnić and Ryzhkov 1999; Ryzhkov et al. 2005b). Because polarimetric radars are quite sensitive to certain phenomena, including hydrometeor phase transitions and size sorting, it is important to understand how such microphysical processes are detectable with measured polarimetric variables. For example, the important polarimetric “bright band” signature associated with melting graupel and snow is well documented (e.g., Brandes and Ikeda 2004; Giangrande et al. 2005, 2008). Additionally, polarimetric observations of size sorting can reveal important information about storm kinematics, including the low-level storm-relative helicity in supercell environments (Kumjian and Ryzhkov 2009).

In contrast, the impacts of warm rain precipitation processes on the polarimetric variables have received comparatively little attention. As rain falls from a precipitating cloud, the distribution of mass among different-sized drops is governed by several microphysical processes. The evolution of this drop size distribution (DSD) is a complex problem that has received considerable attention in the literature over the past few decades. Subcloud microphysical processes affecting the shape of the DSD at the ground include differential sedimentation, spontaneous breakup, collisional breakup, coalescence, and vapor diffusion (i.e., evaporation). Early works by Young (1975), Srivastava (1978), and Johnson (1982) determined that spontaneous breakup is relatively unimportant in comparison with collisional breakup, especially for higher rainwater contents. Laboratory studies by Low and List (1982a) and Beard and Ochs (1995) have facilitated our understanding of collisional processes such as coalescence and breakup. Based on their experimental data, Low and List (1982a,b) developed coalescence efficiencies and a parameterization for the fragment size distribution of particles resulting from filament, sheet, and disk breakup following the collision of two raindrops. These and similar parameterizations have been widely analyzed, improved, and utilized in zero-dimensional box models and one-dimensional rainshaft models (e.g., Gillespie and List 1976; Low and List 1982b; Brown 1986, 1987; List and McFarquhar 1990). Evaporation beneath the cloud base was included in the models by List et al. (1987), Tzivion et al. (1989), Brown (1993, 1994), Hu and Srivastava (1995, hereinafter HS), and Seifert (2008).

These studies and others have found that the collisional processes of coalescence and breakup tend to dominate the evolution of the DSD shape. These processes tend to drive an arbitrary initial DSD toward a family of equilibrium shapes that are related through simple multiplicative factors (List et al. 1987; Brown 1987, 1993; HS). The latter authors found that DSD evolution can be categorized into two phases: 1) collisional processes dominating the evolution of the spectrum as it approaches equilibrium, followed by 2) evaporation dominating the spectral evolution, smoothing maxima, and decreasing the overall mass (while only slowly changing the shape of the distribution). If the initial drop size distribution is already close to its equilibrium shape, HS found that the first stage does not occur. Note that some of the equilibrium DSDs in the aforementioned works are the result of artifacts in the original Low and List (1982a,b) parameterization; the parameterization by McFarquhar (2004) based on the Low and List data alleviates some of these shortcomings and produces equilibrium distributions that are quite different than those found in the previous studies.

The collisional processes and differential sedimentation by themselves do not contribute to or deplete the liquid water content in the subcloud layer: they simply redistribute mass to different sizes. However, net evaporation (the diffusion of water vapor away from drops) depletes the total rainwater content, which has significant implications for quantitative precipitation estimation (QPE). In addition, the generation of negative buoyancy via evaporational cooling plays an important role in storm evolution, including the production of severe downdrafts (e.g., Srivastava 1985, 1987) and even possibly affecting a supercell storm’s likelihood of producing a tornado (e.g., Markowski et al. 2002, 2003; Grzych et al. 2007). Despite the many studies that quantify the impacts of environmental conditions on the rate of evaporation and how evaporation affects the rainfall rate, drop size distribution, downdrafts, and radar reflectivity factor (e.g., Srivastava 1985, 1987; Rosenfeld and Mintz 1988; HS), there is a paucity of studies investigating how varying evaporation rates are manifest in the polarimetric variables. A notable exception is Li and Srivastava (2001), who quantified the impacts of evaporation on differential reflectivity, *Z*_{DR}. In contrast, the purpose of this paper is to quantify the sensitivity of all polarimetric variables to various environmental thermodynamic conditions, and variations in DSDs for light rainfall rates, as well as to formulate recommendations to aid in hydrometeorological rainfall estimation.

Most frequently, such rainfall estimates are made by using remote sensing techniques, especially radar (and, soon, dual-polarization radar). Hence, this paper will focus on evaporation. As such, this work should not be viewed as an attempt to investigate DSD evolution; rather, our goal is to quantify the impacts of evaporation on the polarimetric variables as a preliminary step in understanding how warm rain microphysics are manifested in the radar measurements. The inherent limitations of such an approach reduce the general applicability of such a study. Nonetheless, the results of this work should help with the aforementioned questions and may contribute to improving QPE.

The following section will outline the polarimetric variables and their physical interpretation as well as the physics of the evaporation of raindrops. The qualitative impacts of evaporation on the polarimetric variables will be conceptualized to provide a framework for the quantitative analysis conducted in the remainder of the paper. This quantification will be done using an explicit microphysics model that is described in section 3. Section 4 details the sensitivity tests, outlining the model sensitivity and exploring the parameter space. The model is then used to perform simulations in different evaporation scenarios in section 5. The impacts of evaporation on rainfall estimation are explored in the discussion in section 6, followed by a brief summary of the important conclusions in section 7.

## 2. Background

### a. Polarimetric variables

Following the nationwide upgrade of the WSR-88D radar network, meteorologists will have for their use a full slate of polarimetric variables in addition to the conventional radar reflectivity factor at horizontal polarization (*Z _{H}*), Doppler velocity (

*υ*), and spectral width (

_{r}*σ*). These additional variables are the differential reflectivity (

_{υ}*Z*

_{DR}), differential phase (Φ

_{DP}) and specific differential phase (

*K*

_{DP}), and the copolar cross-correlation coefficient at zero lag (

*ρ*

_{HV}). Here, we will briefly describe the meaning of these measurands, but the reader is directed to more thorough descriptions in the literature (e.g., Herzegh and Jameson 1992; Doviak and Zrnić 1993; Zrnić and Ryzhkov 1999; Straka et al. 2000; Ryzhkov et al. 2005b).

Seliga and Bringi (1976) first proposed the use of the ratio of backscattered power at orthogonal polarizations for observations of rainfall. The so-called differential reflectivity *Z*_{DR} is the difference between reflectivity factors at horizontal and vertical polarizations. Differential reflectivity *Z*_{DR} is dependent on particle size, shape, orientation, and dielectric constant but is independent of concentration. Since raindrops become increasingly oblate with increasing diameter (e.g., Pruppacher and Beard 1970; Pruppacher and Pitter 1971), *Z*_{DR} is considered a good measure of the median drop size of a distribution. For spherical or randomly oriented hydrometeors, *Z*_{DR} is 0 dB.

The use of a differential propagation phase Φ_{DP} for rainfall measurements was introduced by Sachidananda and Zrnić (1986, 1987). Due to the oblateness of raindrops, the horizontally polarized wave lags progressively behind the vertically polarized wave in rain, resulting in a positive differential phase shift. The radar measures the cumulative phase shift along each radial. Because Φ_{DP} is immune to radar calibration, attenuation, and beam blockage (Zrnić and Ryzhkov 1999), and is not biased by noise, it is an attractive variable for quantitative precipitation estimation and attenuation correction.

A more convenient way of interpreting the path-integrated differential phase shift is its range derivative, the specific differential phase *K*_{DP}, which provides the phase shift per unit distance in the radial direction. Here, *K*_{DP} increases as the oblateness and dielectric constant of hydrometeors increase and is not affected by spherical or isotropic hydrometeors such as tumbling dry hail (Balakrishnan and Zrnić 1990). Sachidananda and Zrnić (1987) found that *K*_{DP} is almost linearly related to rainfall rate. In addition, *K*_{DP} is less sensitive to the particle diameter than radar reflectivity factor (*Z _{H}* ∼

*D*

^{6}whereas

*K*

_{DP}∼

*D*

^{4.24}); thus,

*K*

_{DP}is more sensitive to the smaller drop sizes than

*Z*and

_{H}*Z*

_{DR}, which are dominated by larger particles within the radar sampling volume.

The correlation coefficient between the horizontally and vertically polarized signals at zero lag *ρ*_{HV} is useful for discriminating between meteorological and nonmeteorological echoes such as biological scatterers (Zrnić and Ryzhkov 1999), chaff (Zrnić and Ryzhkov 2004), tornadic debris (Ryzhkov et al. 2005c), and smoke plumes (Melnikov et al. 2008). In pure rain at S band, *ρ*_{HV} generally does not fall below 0.98. Further decreases in the measured *ρ*_{HV} are due to the diversity of the particle sizes, orientations, shapes, and irregularities, and to the phase compositions within the radar sampling volume. For example, observations of rain mixed with melting hail will produce lower *ρ*_{HV} than pure rain or pure dry hail. Hydrometeors of the size at which resonance scattering effects occur will also contribute to substantially lower values of *ρ*_{HV}. Resonance scattering, or Mie scattering, occurs when the combined microwave radiation reflected from the near side of the particle and the radiation that penetrates the particle and is reflected from its back side interfere in such a way that the resulting backscattered wave is augmented or diminished.

### b. Evaporation

*r*is the radius of the drop,

*f*is the ventilation coefficient for the vapor diffusivity,

_{υ}*D*is the molecular diffusion coefficient, and

_{υ}*ρ*and

_{υ}*ρ*are the vapor densities of the ambient environment and at the surface of the drop, respectively. Following Pruppacher and Klett (1978) and Rogers and Yau (1989), one can derive an approximate expression describing the rate of change of the drop radius:

_{υr}*S*is the saturation ratio, and

*F*and

_{K}*F*are terms related to heat conduction and vapor diffusion, respectively, that include ventilation effects and are detailed in the appendix. For now, we simply focus on a qualitative interpretation of (2). For subsaturated environments (where

_{D}*S*< 1), the drop radius will decrease, indicating decay through evaporation. Also note that

*dr*/

*dt*is inversely proportional to the drop radius, meaning that the radii of smaller drops will decrease more rapidly under evaporation than larger drops. This well-known fact is important when considering the impacts of evaporation on polarimetric variables.

Theoretically, the preferential depletion of smaller drops will result in a decrease in the observed *Z _{H}* and

*K*

_{DP}with an

*increase*in the observed

*Z*

_{DR}. The observed decrease in

*Z*has been well documented and is intuitive. Decreasing drop diameters across the spectrum and a decrease in the concentration of smaller drops (those that are totally evaporated) will result in a decreased magnitude of the backscattered signal. Recall that

_{H}*K*

_{DP}is less sensitive than

*Z*to large drops, and as a corollary, more sensitive to changes in the lower end of the drop size spectrum. Therefore, one can expect evaporation to affect

_{H}*K*

_{DP}more substantially than

*Z*

_{H}.The expected increase in *Z*_{DR} is less intuitive because all drops are losing mass (size). However, since *Z*_{DR} is a measure of the median drop size in a distribution, a preferential depletion of smaller drops (which generally have a large concentration) causes an *increase* in the median drop size of a given DSD. This effect is shown schematically in Fig. 1. At S band, the change in *ρ*_{HV} due to evaporation in pure rain is not expected to be significant for reasons discussed in a later section. The magnitude of the changes in all polarimetric variables should be dependent on the relative contributions of small drops and large drops and thus is strongly dependent on the DSD.

## 3. Evaporation model

### a. General description

In an effort to quantify the impacts of evaporation on the polarimetric variables under different conditions and assumptions, a simple numerical model is constructed. The idealized one-dimensional model explicitly computes the changes in size of raindrops falling through subsaturated air. The model domain is the subcloud layer, with 100-m vertical resolution. Eighty initial drop sizes are considered, ranging from 0.05 to 7.95 mm in 0.1-mm increments. Each drop size “bin” is tracked independently in order to isolate the effects of evaporation. Hence, no drop interactions such as collisions, coalescence, or breakup are taken into account. Note that coalescence and breakup significantly contribute to the evolution of the drop size spectrum, as found in the numerous theoretical and modeling studies mentioned in the introduction. These collisional processes become increasingly important in heavier rainfall. On the other hand, evaporation tends to only change the slope of the DSD slowly, instead mainly affecting the total water content (e.g., Srivastava 1978; HS). For rainfall estimation, evaporation is important since it is the only subcloud process that directly affects the total mass of rainwater reaching the ground. The relative contributions to the depletion of total water content from combinations of coalescence, breakup, and evaporation have been investigated previously (e.g., HS; Seifert 2008). HS found that the effects of coalescence and breakup tend to approximately balance. In their model, the total rainwater mass depletion in simulations that employed full microphysics (coalescence, breakup, and evaporation) was similar to that in which only evaporation was considered. Thus, only including evaporation in this model, although inherently limiting its applicability toward simulating the evolution of the DSD, is justifiable to improve the computational efficiency given our focus on the radar measurements and associated rainfall estimation.

At the top of the domain, or in the “cloud,” any DSD model can be prescribed. In the subcloud domain, any vertical profile of temperature and relative humidity can be administered. In our model, the feedback on the environmental thermodynamic profiles due to evaporation may be turned on or off. The feedback adjustment scheme is described below. Since the model is one-dimensional, no size sorting due to vertical shear (e.g., Kumjian and Ryzhkov 2009) is considered. At the initial time, drops begin to fall into the top of the domain, as in many rainshaft models. After numerous tests, the time step was selected to be 0.25 s, maximizing the computational efficiency while maintaining stable solutions.

### b. Equations

*α*= 3.78 m s

^{−1}mm

^{−0.67},

*β*= 0.67, and the equivalent spherical diameter

*D*is given in millimeters. The multiplicative factor (

*ρ*

_{0}/

*ρ*)

_{a}^{0.4}is a density correction, where

*ρ*

_{0}is the surface reference density value following Foote and duToit (1969). Though (3) is less accurate for large drops than more recent empirical models (e.g., Brandes et al. 2002), it well represents the fall speeds of smaller drops (<5 mm) for which the evaporation effects are most significant. Thus, performing the change of variables in (2) and using the velocity relation (3) yields an analytic expression for the change in diameter of a raindrop with initial size

*D*

_{0}as a function of height:

*ξ*as the right-hand side of (2). To numerically integrate this expression,

*ξ*and

*ρ*are assumed to be constant over each height step, Δ

_{a}*h*=

*h*−

_{i}*h*

_{i+1}(where

*h*

_{i+1}<

*h*), which should be relatively small (100 m is used in this study), resulting in

_{i}*n*(

*m*) in many of the aforementioned modeling studies is governed by the following equation:

*m*due to growth or decay by vapor diffusion,

*ṁ*(i.e.,

*dm*/

*dt*, which is the time rate of change of the drop mass due to condensation or evaporation). On the right-hand side, the terms

*B*and

*C*account for the collisional processes of breakup and coalescence, respectively; these will be ignored in this study. For the calculation of radar variables, it is more convenient to work in terms of diameter

*D*. Thus, we seek an expression of (6) in terms of

*N*(

*D*), where

*N*(

*D*)

*dD*is the number of raindrops of size

*D*to

*D*+

*dD*. The concentrations

*n*(

*m*) and

*N*(

*D*) are related by

*m*= (

*π*/6)

*ρ*

_{l}D^{3}, and thus

*ρ*is the density of liquid water. When computing the polarimetric radar variables, drop shape is important and must be dealt with appropriately. However, assuming a spherical shape of drops does not significantly impact the computation of evaporation (Straka and Gilmore 2006). Thus, the drop shape inconsistency will not affect the results of this study. The mass diffusion term must be written in terms of

_{l}*D*:

*D*:

*N*=

*N*(

*D*,

*z*,

*t*). For brevity, the functional dependence of

*N*will herein be omitted. Thus, expanding the derivatives in (10) results in

*υ*,

_{t}*Ḋ*= 4

*ξD*

^{−1}, we can write (12) in discrete form,

*i*, height level

*j*, and drop size bin

*k*.

As raindrops begin to evaporate in subsaturated air, the thermodynamic properties of the air in the rainshaft are affected. Specifically, evaporation causes the shaft to gradually cool and moisten in time, thereby decreasing the subsaturations and limiting the amount of evaporation that occurs. The decrease in temperature and precipitation loading contributes to negative buoyancy relative to the surrounding precipitation-free air, which forces a localized downdraft. The downdraft also serves to limit the evaporation as it carries drops downward faster than their terminal velocities, decreasing the amount of time available for evaporation. Srivastava (1985, 1987) has showed how intense downdrafts can occur due to a combination of evaporation and melting of ice particles (both of which have a cooling effect on the local environment).

*f*is the heat ventilation coefficient,

_{h}*K*is the thermal conductivity of the air, and

*T*is the temperature of the drop’s surface. For simplicity, we assume the raindrop temperatures to be in equilibrium at each level, neglecting the brief relaxation period toward their equilibrium temperature. Recently, Tardif and Rasmussen (2010) have investigated the validity of the equilibrium assumption and found that it may not hold if drops experience thermodynamic changes on time scales comparable to the equilibrium relaxation time. In our study, the height step (100 m) ensures that the drops remain within one grid box (in which thermodynamic conditions are constant) over the relaxation time scale, which tends to be on the order of 4–5 s for most drop sizes. Using the total amount of mass evaporated at each level, and assuming that all latent enthalpy from the environment contributes directly to cooling the drop’s surface (which immediately cools the surrounding air), the incremental change in the temperature of the air is calculated at each time step. Similarly, the entire mass of liquid water lost due to evaporation goes directly into vapor, increasing the vapor density in the air. Using the new temperature and vapor density, and the ideal gas law for water vapor, the new relative humidity can be calculated.

_{r}This environmental feedback scheme may be turned on or off in the model. Since the model is one-dimensional and is not tied to any dynamics, we have omitted the generation of a downdraft. The velocity relation (3) overestimates the fall speeds of the moderate and large drops already, which may partially account for the errors introduced when neglecting a downdraft. Also, to be realistic and dynamically consistent, the generation of a downdraft would require some parameterization of entrainment, a process that still has some uncertainty. Without this entrainment, the 1D column gradually moistens and cools until it becomes saturated, halting the evaporation.

*D*is in millimeters. See Brandes et al. (2002) for a discussion of the differences between this relation and others that are used frequently in the literature. Note that this relation produces slightly more spherical drops, especially in the range 1–4 mm. The polarimetric variables (especially

*Z*

_{DR}and

*ρ*

_{HV}) are more strongly affected by the larger drops, the shapes of which are similar to other relations. Thus, we expect only slight differences if other drop shape relations are used instead of (15). The raindrops are assumed to have a distribution of canting angles with mean 0° and standard deviation 20° with respect to the vertical. In general, 20° is too high for typical rainfall; a standard deviation of 10° is more appropriate (e.g., Ryzhkov 2001; Ryzhkov et al. 2002). However, the higher value was chosen to increase the effects of evaporation on

*ρ*

_{HV}. The impact of this wider distribution of canting angles is to dampen the effects of evaporation on the other variables very slightly. The wavelength of the radar will be discussed using convention, according to which “S band” refers to the 10.9-cm wavelength (as in the WSR-88D network) and “C band” refers to a wavelength of 5.3 cm.

The evaporation model calculations will be separated into two groups: sensitivity studies, which will focus on the exploration of parameter space, and simulations, which will model the vertical profiles of polarimetric variables using observed thermodynamic soundings from different climate regions. The sensitivity studies are designed to test the model sensitivity as well as the sensitivity of the polarimetric variables to different thermodynamic conditions, drop size distribution models, and rainfall rates. The simulations are intended to illustrate potential variations in vertical profiles of polarimetric variables in different regions of the United States that may be observed following the nationwide weather radar upgrade to polarization diversity. For the experiments, the gradual moistening and cooling feedback mechanism has been turned off in order to create “snapshots” of the profiles of polarimetric variables. For applied simulations, the mechanism should be turned on to avoid overaggressive evaporation.

## 4. Sensitivity studies

### a. Design

*N*

_{00}= 8000 m

^{−3}mm

^{−1}, is calculated for use at sea level; since we are using fixed rainfall rates aloft, we must use

*N*

_{0}=

*N*

_{00}[

*ρ*(

*h*

_{0})/

*ρ*

_{0}]

^{0.4}. For the MP distribution, with the velocity relation (3), the slope parameter Λ is calculated using

*μ*and the slope parameter Λ:

*μ*varying between −1 and 5 are used in the calculations. These values encompass the bulk of the observations in Oklahoma rain (e.g., Cao et al. 2008). A slightly broader range of values of

*μ*was found in rain in Florida (Zhang et al. 2001), so there is some dependence on climate region. Since Λ is determined from the preselected

*μ*values using (19) and we have a fixed rainfall rate

*R*,

*N*

_{0}using (17) for

*N*(

*D*):

*μ*= 1 distribution, which nearly overlap for all drop sizes except at the small-drop end of the spectrum.

For each set of model runs the polarimetric variables are calculated for both S and C bands. In the first set of experiments, idealized isothermal layers 2 km deep are considered. The temperature is fixed at 20°C, and the relative humidity (RH) is set constant in height but varies from 10% to 95% in each experiment. Though not necessarily realistic, these experiments simply explore the parameter space of the model. In the second set of sensitivity runs, a dry-adiabatic lapse rate is prescribed with a surface temperature of 30°C. The RH profile is assumed to increase linearly from the given surface value to 100% at the cloud base (top of the domain, taken to be 3 km). The surface values are varied from 55% to 95% in 5% increments. These idealized well-mixed boundary layer cases apply to warm season precipitation events. Simulations using observed thermodynamic profiles, including both warm and cool season events, are explored in section 5. In these two sets of sensitivity tests, each of the DSD models is used. The final set of sensitivity tests varies *R* from 0.1 to 20 mm h^{−1} using the MP DSD for a dry-adiabatic environment with a surface RH of 75% and a surface temperature of 30°C.

### b. Transient differential sedimentation

At the onset of precipitation, a transient effect known as differential sedimentation occurs, where the largest raindrops fall to the ground before the smaller raindrops. The reason for this is simply the difference in the terminal velocity of the large drops versus the smaller drops; the large drops fall faster and thus reach the ground before the smaller drops. Once enough time has elapsed, all drops that do not evaporate reach the ground, and the visible effects of differential sedimentation disappear.

Recall that size sorting is an extremely important effect as observed by polarimetric radars (e.g., see a review by Kumjian and Ryzhkov 2008). With the largest, most oblate drops initially falling farther from the cloud than smaller drops, the polarimetric variables (especially *Z*_{DR}) are significantly enhanced beneath the cloud, with *Z*_{DR} increasing toward the ground. At S band, *ρ*_{HV} tends to change very little, only slightly decreasing toward the ground. This is because the resonance effects are still minimal even for the largest drops at S band. At C band, however, the large drops are within the band of sizes where resonance effects are quite strong. Thus, in addition to increased *Z*_{DR} toward the ground, C-band measurements will show a significant decrease in *ρ*_{HV}. Figure 3 shows an example of the impacts of differential sedimentation on the polarimetric variables. Note that this effect dominates any signal from evaporation, which produces more subtle changes in the vertical profiles of the polarimetric variables. Thus, for the remainder of the sensitivity studies and experiments where we quantify the effects of evaporation, the calculations are carried out for long times (equivalently, waiting for the profiles to attain steady-state solutions where the environmental feedback mechanism has been turned off).

### c. Results of sensitivity tests

^{1}The evaporative change will be denoted as a Δ before the variable. As discussed above, evaporation will produce negative Δ

*Z*and Δ

_{H}*K*

_{DP}(indicative of a decrease) and a positive Δ

*Z*

_{DR}(indicative of an increase). To make a more meaningful comparison between Δ

*Z*and Δ

_{H}*K*

_{DP}, the evaporative change in

*K*

_{DP}will be converted into logarithmic units:

*Z*,

_{H}*Z*

_{DR}, and

*K*

_{DP}can be converted into relative changes (e.g., a 3-dB decrease corresponds to about a 50% reduction from its original value).

The results from the first set of experiments (considering isothermal layers) are summarized in Figs. 4 and 5 for S and C bands, respectively. It is clear that all of the variables are sensitive to the relative humidity in the layer, but perhaps more importantly the results are sensitive to the initial DSD model selected. The Γ, *μ* = 5 model exhibits much greater evaporative changes in *Z _{H}* and

*K*

_{DP}than do the other models (Figs. 4a,c and 5a,c). This is explained by two factors. First, the distribution contains a large concentration of small drops, which are preferentially evaporated, resulting in a substantial decrease in mass. Second, the Γ,

*μ*= 5 has fewer large drops than any other distribution. These large drops, which do not evaporate as efficiently as smaller drops, tend to overwhelm the contribution to the observed

*Z*(and, to a lesser extent,

_{H}*K*

_{DP}) for the other DSD models, which all exhibit lower magnitudes of Δ

*Z*and Δ

_{H}*K*

_{DP}(dB).

Since *Z*_{DR} is more sensitive to drop size than *Z _{H}*, it follows that the large Δ

*Z*

_{DR}values occur for the DSDs with relatively large concentrations of big drops. However, also playing an important role in producing the significant Δ

*Z*

_{DR}values is a large concentration of small drops. The preferential evaporation of a significant portion of the spectrum will substantially increase the median drop size of the spectrum. This is why the MP model has the highest Δ

*Z*

_{DR}(Figs. 4b and 5b): it has the highest concentration of small drops and a comparatively large concentration of big drops. Also note that the Δ

*Z*

_{DR}at C band (Fig. 5b) is somewhat higher for the DSD models due to resonance scattering effects.

At S band, the Δ*ρ*_{HV} magnitudes are quite small (<0.002 even for extreme evaporation) for all DSD shapes, at least for the electromagnetic scattering model employed in this study (Fig. 4d). Such small changes are insignificant and likely within the uncertainty of the WSR-88D measurements, or about ±0.005–0.01. At C band, the evaporative changes in *ρ*_{HV} are slightly larger in magnitude and negative for all models (Fig. 5d). For modest evaporation rates, changes in *ρ*_{HV} at both radar wavelengths are insignificant and probably difficult to detect operationally.

Next, idealized 3-km-deep mixed layers are considered for the calculations. The results for S and C bands are summarized in Fig. 6. The resulting Δ*Z _{H}* and Δ

*K*

_{DP}are dependent on the surface RH and the DSD prescribed aloft, as expected. However, a feature of note is that Δ

*Z*and Δ

_{H}*K*

_{DP}reverse sign for three of the DSD models (MP, Γ,

*μ*= −1, and Γ,

*μ*= 1) for high values of surface RH. This is because the relatively cool and moist conditions beneath the cloud base do cause enough evaporation to counteract the effects of “raindrop convergence,” a result of the drops encountering increasing air density with fall distance. In the absence of any evaporation, the concentration of drops increases slightly with decreasing height, which affects

*Z*and

_{H}*K*

_{DP}, but not

*Z*

_{DR}or

*ρ*

_{HV}. Also of note is that the difference between the S- and C-band values of Δ

*Z*

_{DR}and Δ

*ρ*

_{HV}increases with higher concentrations of large drops, a result of the resonance scattering effects prevalent in big drops at C band.

The results of varying the rainfall rate *R* are provided in Fig. 7. For the higher rainfall rates considered, one should expect a greater change in *K*_{DP} for a given amount of evaporation than is observed in *Z _{H}*. At both S and C bands, Δ

*ρ*

_{HV}is insignificant for all rainfall rates with the MP model. For the S-band calculations, Δ

*Z*

_{DR}decreases with increasing rainfall rate whereas the Δ

*Z*

_{DR}at C band increases slightly. The difference in behavior may be attributable to the enhanced resonance scattering effects at C band.

The sensitivity experiments have shown, for given thermodynamic conditions, that the largest Δ*Z _{H}* and Δ

*K*

_{DP}occur for the Γ,

*μ*= 5 DSD model, with the Γ,

*μ*= −1 model having the smallest evaporative changes. Again, these results are directly related to the shape of the DSD: the smaller the concentration of large (>4 mm) drops, the greater the impacts of evaporation will be on

*Z*and

_{H}*K*

_{DP}. The logarithmic changes in

*K*

_{DP}are slightly greater in magnitude than those of

*Z*since

_{H}*K*

_{DP}is less dependent on drop diameter than

*Z*; the preferential depletion of smaller drops has a greater relative impact on

_{H}*K*

_{DP}. Note that since

*Z*typically is displayed in logarithmic units (dB

_{H}*Z*), whereas

*K*

_{DP}generally is shown in degrees per kilometer, Δ

*K*

_{DP}may be more obvious to operational meteorologists. For example, an equivalent evaporative change of 3 dB for

*Z*(e.g., 33–30 dB

_{H}*Z*) “appears” to be less significant than a 3-dB relative change in

*K*

_{DP}(e.g., 1.0°–0.5° km

^{−1}). However, it should be noted that for the relatively light rainfall rates considered in this study and at S band,

*K*

_{DP}can be noisy and difficult to estimate. At shorter radar wavelengths, the

*K*

_{DP}estimates will be less noisy as the

*K*

_{DP}values will be larger (since

*K*

_{DP}is inversely proportional to radar wavelength), possibly making such evaporative changes easier to detect operationally with C- and X-band polarimetric radars.

For given conditions, the largest values of Δ*Z*_{DR} occur for the MP model, with Γ, *μ* = −1 having the smallest values. The changes depend on the relative contributions from both the small- and large-drop ends of the spectrum. While having the largest concentration of big drops of any of the models used in this study, the Γ, *μ* = −1 model resulted in the smallest Δ*Z*_{DR}. This is because of the comparatively low concentration of small drops (<2 mm). With fewer small drops being depleted, there is a lesser shift in the median drop size of the spectrum. Coupled with this is the fact that the substantial concentrations of big drops, which are not as affected by evaporation, dominate the contribution to the backscattered signal.

## 5. Simulations with observed soundings

### a. Design

Once the model sensitivity is explored, it is informative to look at real cases in an attempt to see the types of variations in vertical profiles of the polarimetric variables that may be observed in different evaporation scenarios. To accomplish this, we have selected three soundings from different regions that display different thermodynamic characteristics. These observed temperature and relative humidity profiles are used to initialize the model.

The three soundings are shown in Fig. 8. The first sounding comes from Albuquerque, New Mexico, at 0000 UTC 16 August 2007 (Fig. 8a) and represents the summer season in the western United States with a deep, dry well-mixed layer. Light precipitation echoes were observed by the Albuquerque WSR-88D about the time of the sounding. Next is a spring sounding (Fig. 8b), taken from Norman, Oklahoma, at 1800 UTC 24 May 2008 and may be said to represent an atmosphere primed for severe weather: a cyclic tornadic supercell developed north of Oklahoma City, Oklahoma, in the afternoon, producing at least 10 tornadoes. The third sounding (Fig. 8c) comes from Wallops Island, Virginia, at 0000 UTC 7 October 2009 and represents a cool, moist boundary layer in which light rain was observed throughout the day. Note that the environmental feedback mechanism has been turned off for these simulations.

### b. Results of simulations

The MP distribution with a 5 mm h^{−1} rainfall rate aloft is used in the simulations. The MP model is employed because it is perhaps the most widely known DSD model and provides intermediate magnitudes of evaporative changes. As previously demonstrated, the choice of DSD model can have a significant impact on the evaporative changes of the polarimetric variables. This should be considered when interpreting the outcome of the simulations and the computed vertical profiles of the polarimetric variables presented herein. Nonetheless, the purpose of these simulations is to illustrate the potential variations in vertical profiles of polarimetric observables that may occur in different climate regions and seasons throughout the United States.

Figure 9 displays the vertical profiles of polarimetric variables for the Albuquerque sounding (Fig. 8a). For both radar wavelengths, Δ*Z _{H}* is less than −6 dB. The Δ

*Z*

_{DR}is nearly 0.3 dB at S band, which may be detectable by the WSR-88D, and greater than 0.3 dB at C band. The magnitude of the change in

*ρ*

_{HV}is negligible at S band. At C band, the Δ

*ρ*

_{HV}is about −0.003, which may be observable with very high quality radar systems. The evaporative changes in all polarimetric variables for each of the simulated soundings are summarized in Table 1.

Because of the uncertainty involved with the choice of the initial DSD, the same simulations are performed with the Γ, *μ* = −1 distribution as well as the Γ, *μ* = 5 model. These models were selected for comparison since they encompass the extremes for evaporative changes in the different polarimetric variables (Table 2). The Γ, *μ* = −1 distribution exhibits much smaller magnitudes of Δ*Z _{H}* and Δ

*K*

_{DP}than do the MP simulations. Also note that Δ

*ρ*

_{HV}changes sign with the Γ,

*μ*= −1 model depending on the radar wavelength. In contrast, the Γ,

*μ*= 5 model produces the most significant Δ

*Z*and Δ

_{H}*K*

_{DP}values and very small changes in

*Z*

_{DR}and

*ρ*

_{HV}. These simulations can be viewed as the approximate quantitative bounds on the possible changes in polarimetric characteristics of light rain due to evaporation in differing environments.

## 6. Rainfall-rate estimation

*R*is given in millimeters per hour and the lowercase subscripts for

*Z*and

_{h}*Z*

_{dr}indicate values in linear units. Understanding how radar-estimated rainfall rates may change in height due to evaporation is important for QPE, especially in regions far from the radar or in regions of complex terrain where low levels may not be adequately sampled (e.g., Maddox et al. 2002). Computing rainfall rates based on accepted algorithms [(23)–(25)] for observed soundings provides an illustration of the types of variations possible across the country. Additionally, comparing these estimated rainfall rates with the actual rainfall rate demonstrates the performance of each of the algorithms in situations where the DSD is changing with height through evaporation.

Vertical profiles of the rainfall rate are plotted in Fig. 10. An MP distribution is assumed at the top of the domain with a rainfall rate of 5 mm h^{−1}. In all cases, the true rainfall rate changes more rapidly than do any of the estimates. Note that the initial large decrease in *R* at Wallops Island is due to the dry layer near 700 mb (see Fig. 8c). At the top of the domain, the *R*(*Z _{H}*,

*Z*

_{DR}) estimate is closest to the actual

*R*, but as the raindrops evaporate, the

*R*(

*K*

_{DP}) better matches the true

*R*. It should be noted, however, that actual measurements of

*K*

_{DP}can be noisy in light rain, so

*R*(

*K*

_{DP}) is not necessarily the best estimate for operational use in the case of evaporation and light rain. The evaporative changes at the ground using different rainfall rate estimates are summarized in Table 3. Overall, the

*R*(

*Z*,

_{H}*Z*

_{DR}) relation best captures the change in rainfall rate with height (i.e., it is least sensitive to the changes in the DSD at low rainfall rates, which was frequently claimed as an advantage of the polarimetric radar estimates), though all relations underestimate the change as calculated from the actual DSD. It is clear that substantial decreases in rainfall rate are possible in dry environments, and even in the moist Wallops Island environment the true rainfall rate decreased by almost 3 mm h

^{−1}. If the rainfall rate at low levels is not adequately estimated based on radar measurements, this error can accumulate with time.

*R*(

*Z*,

_{H}*Z*

_{DR}) relations suggested in the literature (see Ryzhkov et al. 2005a for a review). For comparison, the calculations are reproduced comparing the

*R*(

*Z*,

_{H}*Z*

_{DR}) relation from (25) to the one recommended by Giangrande and Ryzhkov (2008),

*R*on

*Z*

_{DR}in (25) and (27): the subtle changes in shape of the DSD are not captured by

*Z*, whereas

_{H}*Z*

_{DR}does provide some information about the shape of the DSD.

Reproducing the calculations with the Γ, *μ* = −1 model, the relative errors (between the estimates of evaporative change in the rainfall rate and the true evaporative change in the rainfall rate) depend on the sounding but are all within 2 mm h^{−1} of the true ΔR. In contrast, these errors are minimized with the *R*(*Z _{H}*,

*Z*

_{DR}) relation for the Γ,

*μ*= 5 and MP models. The difficulty for operational forecasters is in determining the initial DSD aloft. There are rigorous methods for DSD retrieval based on polarimetric radar measurements, such as the procedures described in Zhang et al. (2001, 2006) and Brandes et al. (2002, 2004). These methods mainly rely on

*Z*and

_{H}*Z*

_{DR}measurements.

Using the idealized well-mixed layer environments, and initializing the model with a variety of constrained gamma DSD models with varying rainfall rates (0.1−20 mm h^{−1}) and shape parameters (−1 ≤ *μ* ≤ 5), the relation between the initial *Z _{H}* and

*Z*

_{DR}aloft and the evaporative change in rainfall rate over the domain (in this case, 3 km) is explored. With the constrained gamma DSD models, the

*relative change in rainfall rate due to evaporation*(Δ

*R*/

*R*)

*is independent of the initial Z*. Thus, Δ

_{H}aloft*R*/

*R*is a function of the environment (including temperature, RH, and depth) and initial

*Z*

_{DR}aloft. Observations of

*Z*and

_{H}*Z*

_{DR}at cloud base can be used to identify the initial rainfall rate aloft. Obviously, these parameters may be a function of time in evolving storms, but for QPE one must assume a steady state for at least the duration of a radar volume scan. Next, the RH (or another moisture variable) and temperature profiles in the layer beneath the radar horizon can be determined from observations or numerical weather prediction model output [e.g., the Rapid Update Cycle (RUC); e.g., Benjamin et al. (1991), (1994), (2004)], assuming it is sufficiently accurate. In this manner, environmental conditions and the observed

*Z*

_{DR}(which does not change much with height) can be utilized efficiently to estimate Δ

*R*/

*R*(Fig. 11). Such a process may aid in QPE, especially in regions far from the radar or in poor low-level radar coverage.

## 7. Summary of conclusions

In this paper we have investigated and quantified the impacts of evaporation on the polarimetric variables that will be available following the national upgrade of the WSR-88D radar network, which include *Z _{H}*,

*Z*

_{DR},

*K*

_{DP}, and

*ρ*

_{HV}. Evaporation results in significant decreases in

*Z*and

_{H}*K*

_{DP}and subtle increases in

*Z*

_{DR}, with no significant change in

*ρ*

_{HV}at S band. At C band, changes in

*K*

_{DP},

*Z*

_{DR}, and

*ρ*

_{HV}are amplified for a given amount of evaporation. The change in the true rainfall rate due to evaporation is dependent on the initial shape of the DSD aloft, which illuminates the importance of polarimetric radar observations in reducing this uncertainty. To test the sensitivity of the polarimetric variables to the thermodynamic conditions in the subcloud environment, the initial DSD aloft, and the rainfall rate, a simple explicit bin microphysics model was constructed and numerous experiments were conducted. Because it neglects raindrop collisional processes, the model is limited to relatively light rainfall rates and should not be used to simulate DSD evolution. We summarize the conclusions reached from the model runs as follows:

As expected, warmer and/or drier environments produce more evaporation and thus more substantial changes in the polarimetric variables, as well as the rainfall rates inferred from these variables.

The resulting vertical profiles of polarimetric variables from the simplistic model are sensitive to the initial DSD aloft, which is of considerable uncertainty in various precipitation regimes. Thus, polarimetric radar measurements (

*Z*,_{H}*Z*_{DR}) aloft are needed to diminish the uncertainty due to DSD variability. Drop size distributions with large concentrations of smaller drops (<2 mm) and relatively low number concentrations of big drops (>4 mm) exhibit the largest evaporative changes in*Z*and_{H}*K*_{DP}. This is because the small drops are preferentially evaporated, while the backscattered signals are not overwhelmed by contributions from a large number of big drops. Distributions with a comparatively large concentration of big drops and a sufficiently high concentration of smaller drops (such as the MP model) will exhibit the largest increase in*Z*_{DR}. Evaporative changes in*ρ*_{HV}are generally quite small in magnitude and likely unobservable at both S and C bands.When the evaporative changes are converted into logarithmic units, the difference in

*K*_{DP}is larger in magnitude than is that of*Z*. This is because_{H}*K*_{DP}is less sensitive to large drops than*Z*. Since_{H}*K*_{DP}is generally displayed in linear units (° km^{−1}), evaporative changes may be more evident in*K*_{DP}than in*Z*, especially at shorter radar wavelengths. It is important to note that estimates of_{H}*K*_{DP}at S band can be noisy, especially in light rain as considered in this study. The relative changes of*K*_{DP}are about the same at S and C bands.Except for evaporative changes in

*Z*_{DR}at C band, the changes in the polarimetric variables become less significant with increasing rainfall rate. For rainfall rates larger than about 10 mm h^{−1}, the effects of evaporation on the evolution of the drop size spectrum as well as the vertical profiles of polarimetric variables are less significant than for other microphysical processes such as coalescence and breakup, which should be considered the dominant processes and should be included in more general applications.

Additionally, the use of observed soundings to initialize the model produces realistic vertical profiles of the polarimetric variables. Different environments produce different degrees of evaporation that are in principle measurable by operational S- and C-band polarimetric radars. Quantitative bounds to the evaporative changes in the polarimetric variables are suggested using the extremes in the variations in DSD shape. However, we should again emphasize the limitations of the model. Since drop coalescence and breakup are ignored, the model is applicable only to relatively light rain. To more accurately model subcloud microphysics for more general applicability, these effects must be considered. Also, the moistening and cooling of the shaft must be turned “on” for simulation studies, especially when coupled with dynamical models.

In situations where low levels are unobserved by radar, modeling based on thermodynamic information from soundings or numerical weather prediction models such as the RUC may provide guidance for precipitation estimates, allowing for adjustments to be made to observed rainfall rates aloft. Assuming the DSD can be described by a constrained gamma model, the relative change in rainfall rate due to evaporation Δ*R*/*R* is independent of *Z _{H}* aloft. Thus,

*Z*

_{DR}(together with the environmental thermodynamic profiles) can be used in a simple evaporation model to estimate Δ

*R*/

*R*. This method differs from conventional statistical or climatological techniques because the simple modeling approach in this study is physically based. The use of polarimetric measurements provides crucial information that is not available with single-polarization radar measurements of

*Z*. Such adjustments to conventional techniques may improve the estimates of rainfall rates at the surface, benefiting hydrology models and forecasts.

_{H}## Acknowledgments

This work is partially funded under NSF Grant ATM-0532107. Additional funding comes from the NOAA/Office of Oceanic and Atmospheric Research under NOAA/University of Oklahoma Cooperative Agreement NA17RJ1227, U.S. Department of Commerce. We thank Dr. Scott Giangrande for useful comments on the topic. Drs. Heather Reeves and J. J. Gourley provided useful reviews of the manuscript that helped the clarity and focus of the paper. Additionally, three anonymous reviewers provided numerous comments and suggestions that significantly improved the manuscript.

## REFERENCES

Atlas, D., and C. W. Ulbrich, 1977: Path- and area-integrated rainfall measurement by microwave attenuation in the 1–3 cm band.

,*J. Appl. Meteor.***16****,**1322–1331.Balakrishnan, N., and D. S. Zrnić, 1990: Use of polarization to characterize precipitation and discriminate large hail.

,*J. Atmos. Sci.***47****,**1525–1540.Beard, K. V., and H. T. Ochs III, 1995: Collisions between small precipitation drops. Part II: Formulas for coalescence, temporary coalescence, and satellites.

,*J. Atmos. Sci.***52****,**3977–3996.Benjamin, S. G., K. A. Brewster, R. L. Brümmer, B. F. Jewett, T. W. Schlatter, T. L. Smith, and P. A. Stamus, 1991: An isentropic three-hourly data assimilation system using ACARS aircraft observations.

,*Mon. Wea. Rev.***119****,**888–906.Benjamin, S. G., K. J. Brundage, P. A. Miller, T. L. Smith, G. A. Grell, D. Kim, J. M. Brown, and T. W. Schlatter, 1994: The Rapid Update Cycle at NMC. Preprints,

*10th Conf. on Numerical Weather Prediction,*Portland, OR, Amer. Meteor. Soc., 566–568.Benjamin, S. G., and Coauthors, 2004: An hourly assimilation–forecast cycle: The RUC.

,*Mon. Wea. Rev.***132****,**495–518.Brandes, E. A., and K. Ikeda, 2004: Freezing-level estimation with polarimetric radar.

,*J. Appl. Meteor.***43****,**1541–1553.Brandes, E. A., G. Zhang, and J. Vivekanandan, 2002: Experiments in rainfall estimation with a polarimetric radar in a subtropical environment.

,*J. Appl. Meteor.***41****,**674–685.Brandes, E. A., G. Zhang, and J. Vivekanandan, 2004: Drop size distribution retrieval with polarimetric radar: Model and application.

,*J. Appl. Meteor.***43****,**461–475.Bringi, V. N., and V. Chandrasekar, 2001:

*Polarimetric Doppler Weather Radar: Principles and Applications*. Cambridge University Press, 636 pp.Brown, P. S., 1986: Analysis of the Low and List drop-breakup formulation.

,*J. Climate Appl. Meteor.***25****,**313–321.Brown, P. S., 1987: Parameterization of drop-spectrum evolution due to coalescence and breakup.

,*J. Atmos. Sci.***44****,**242–249.Brown, P. S., 1993: Analysis and parameterization of the combined coalescence, breakup, and evaporation processes.

,*J. Atmos. Sci.***50****,**2940–2951.Brown, P. S., 1994: Vertical variation of the steady-state drop spectrum in a one-dimensional rain shaft.

,*J. Atmos. Sci.***51****,**2075–2085.Cao, Q., G. Zhang, E. Brandes, T. Schuur, A. V. Ryzhkov, and K. Ikeda, 2008: Analysis of video disdrometer and polarimetric radar data to characterize rain microphysics in Oklahoma.

,*J. Appl. Meteor. Climatol.***47****,**2238–2255.Doviak, R. J., and D. S. Zrnić, 1993:

*Doppler Radar and Weather Observations*. Academic Press, 562 pp.Foote, G. B., and P. S. duToit, 1969: Terminal velocity of raindrops aloft.

,*J. Appl. Meteor.***8****,**249–253.Giangrande, S. E., and A. V. Ryzhkov, 2008: Estimation of rainfall based on the results of polarimetric echo classification.

,*J. Appl. Meteor. Climatol.***47****,**2445–2462.Giangrande, S. E., A. V. Ryzhkov, and J. Krause, 2005: Automatic detection of the melting layer with a polarimetric prototype of the WSR-88D radar. Preprints,

*32nd Int. Conf. on Radar Meteorology,*Albuquerque, NM, Amer. Meteor. Soc., 11R.2. [Available online at http://ams.confex.com/ams/pdfpapers/95894.pdf].Giangrande, S. E., J. M. Krause, and A. V. Ryzhkov, 2008: Automatic designation of the melting layer with a polarimetric prototype of the WSR-88D radar.

,*J. Appl. Meteor. Climatol.***47****,**1354–1364.Gillespie, J. R., and R. List, 1976: Evolution of raindrop size distributions in steady-state rainshafts. Preprints,

*Int. Conf. Cloud Physics,*Boulder, CO, Amer. Meteor. Soc., 472–477.Grzych, M. L., B. D. Lee, and C. A. Finley, 2007: Thermodynamic analysis of supercell rear-flank downdrafts from project ANSWERS.

,*Mon. Wea. Rev.***135****,**240–246.Herzegh, P. H., and A. R. Jameson, 1992: Observing precipitation through dual-polarization radar measurements.

,*Bull. Amer. Meteor. Soc.***73****,**1365–1374.Hu, Z., and R. Srivastava, 1995: Evolution of raindrop size distribution by coalescence, breakup, and evaporation: Theory and observations.

,*J. Atmos. Sci.***52****,**1761–1783.Johnson, D. B., 1982: Raindrop multiplication by drop breakup.

,*J. Appl. Meteor.***21****,**1048–1050.Kumjian, M. R., and A. V. Ryzhkov, 2008: Microphysical size sorting revealed by dual-polarization Doppler radar. Preprints,

*Symp. on Recent Developments in Atmospheric Applications of Radar and Lidar,*New Orleans, LA, Amer. Meteor. Soc., P2.13. [Available online at http://ams.confex.com/ams/pdfpapers/133439.pdf].Kumjian, M. R., and A. V. Ryzhkov, 2009: Storm-relative helicity revealed from polarimetric radar measurements.

,*J. Atmos. Sci.***66****,**667–685.Li, X., and R. C. Srivastava, 2001: An analytical solution for raindrop evaporation and its application to radar rainfall measurements.

,*J. Appl. Meteor.***40****,**1607–1616.List, R., and G. M. McFarquhar, 1990: The evolution of three-peak raindrop size distributions in one-dimensional shaft models. Part I: Single-pulse rain.

,*J. Atmos. Sci.***47****,**2996–3006.List, R., N. R. Donaldson, and R. E. Stewart, 1987: Temporal evolution of drop spectra to collisional equilibrium in steady and pulsating rain.

,*J. Atmos. Sci.***44****,**362–372.Low, T. B., and R. List, 1982a: Collision, coalescence, and breakup of raindrops. Part I: Experimentally established coalescence efficiencies and fragment size distributions in breakup.

,*J. Atmos. Sci.***39****,**1591–1606.Low, T. B., and R. List, 1982b: Collision, coalescence, and breakup of raindrops. Part II: Parameterization of fragment size distributions.

,*J. Atmos. Sci.***39****,**1607–1618.Maddox, R. A., J. Zhang, J. J. Gourley, and K. Howard, 2002: Weather radar coverage over the contiguous United States.

,*Wea. Forecasting***17****,**927–934.Markowski, P. M., J. Straka, and E. Rasmussen, 2002: Direct surface thermodynamic observations within the rear-flank downdrafts of nontornadic and tornadic supercells.

,*Mon. Wea. Rev.***130****,**1692–1721.Markowski, P. M., J. M. Straka, and E. N. Rasmussen, 2003: Tornadogenesis resulting from the transport of circulation by a downdraft: Idealized numerical simulations.

,*J. Atmos. Sci.***60****,**795–823.McFarquhar, G. M., 2004: A new representation of collision-induced breakup of raindrops and its implications for the shapes of raindrop size distributions.

,*J. Atmos. Sci.***61****,**777–794.Melnikov, V. M., D. S. Zrnić, R. M. Rabin, and P. Zhang, 2008: Radar polarimetric signatures of fire plumes in Oklahoma.

,*Geophys. Res. Lett.***35****,**L14815. doi:10.1029/2008GL034311.Mishchenko, M. I., 2000: Calculation of the amplitude matrix for a nonspherical particle in a fixed orientation.

,*Appl. Opt.***39****,**1026–1031.Pruppacher, H. R., and K. V. Beard, 1970: A wind tunnel investigation of the internal circulation and shape of water drops falling at terminal velocity in air.

,*Quart. J. Roy. Meteor. Soc.***96****,**247–256.Pruppacher, H. R., and R. L. Pitter, 1971: A semi-empirical determination of the shape of cloud and raindrops.

,*J. Atmos. Sci.***28****,**86–94.Pruppacher, H. R., and J. D. Klett, 1978:

*Microphysics of Clouds and Precipitation*. 2nd ed. Oxford University Press, 953 pp.Rasmussen, R. M., and A. J. Heymsfield, 1987: Melting and shedding of graupel and hail. Part I: Model physics.

,*J. Atmos. Sci.***44****,**2754–2763.Rogers, R. R., and M. K. Yau, 1989:

*A Short Course in Cloud Physics*. 3rd ed. Elsevier Press, 290 pp.Rosenfeld, D., and Y. Mintz, 1988: Evaporation of rain falling from convective clouds as derived from radar measurements.

,*J. Appl. Meteor.***27****,**209–215.Ryzhkov, A. V., 2001: Interpretation of polarimetric radar covariance matrix for meteorological scatterers: Theoretical analysis.

,*J. Atmos. Oceanic Technol.***18****,**315–328.Ryzhkov, A. V., D. S. Zrnić, J. C. Hubbert, V. N. Bringi, J. Vivekanandan, and E. A. Brandes, 2002: Polarimetric radar observations and interpretation of co-cross-polar correlation coefficients.

,*J. Atmos. Oceanic Technol.***19****,**340–354.Ryzhkov, A. V., S. E. Giangrande, and T. J. Schuur, 2005a: Rainfall estimation with a polarimetric prototype of WSR-88D.

,*J. Appl. Meteor.***44****,**502–515.Ryzhkov, A. V., T. J. Schuur, D. W. Burgess, P. L. Heinselman, S. E. Giangrande, and D. S. Zrnić, 2005b: The Joint Polarization Experiment.

,*Bull. Amer. Meteor. Soc.***86****,**809–824.Ryzhkov, A. V., T. J. Schuur, D. W. Burgess, and D. S. Zrnić, 2005c: Polarimetric tornado detection.

,*J. Appl. Meteor.***44****,**557–570.Sachidananda, M., and D. S. Zrnić, 1986: Differential propagation phase shift and rainfall rate estimation.

,*Radio Sci.***21****,**235–247.Sachidananda, M., and D. S. Zrnić, 1987: Rain rate estimated from differential polarization measurements.

,*J. Atmos. Oceanic Technol.***4****,**588–598.Seifert, A., 2008: On the parameterization of evaporation of raindrops as simulated by a one-dimensional rainshaft model.

,*J. Atmos. Sci.***65****,**3608–3619.Seliga, T. A., and V. N. Bringi, 1976: Potential use of radar differential reflectivity measurements at orthogonal polarizations for measuring precipitation.

,*J. Appl. Meteor.***15****,**69–76.Srivastava, R. C., 1978: Parameterization of raindrop size distributions.

,*J. Atmos. Sci.***35****,**108–117.Srivastava, R. C., 1985: A simple model of evaporatively driven downdraft: Application to microburst downdraft.

,*J. Atmos. Sci.***42****,**1004–1023.Srivastava, R. C., 1987: A model of intense downdrafts driven by the melting and evaporation of precipitation.

,*J. Atmos. Sci.***44****,**1752–1774.Straka, J. M., and M. S. Gilmore, 2006: Does the influence of oblate-like distortions in larger raindrops make a difference in collection and evaporation parameterizations?

,*J. Appl. Meteor. Climatol.***45****,**1582–1591.Straka, J. M., D. S. Zrnić, and A. V. Ryzhkov, 2000: Bulk hydrometeor classification and quantification using polarimetric radar data: Synthesis of relations.

,*J. Appl. Meteor.***39****,**1341–1372.Tardif, R., and R. M. Rasmussen, 2010: Evaporation of nonequilibrium raindrops as a fog formation mechanism.

,*J. Atmos. Sci.***67****,**345–364.Tzivion, S., G. Feingold, and Z. Levin, 1989: The evolution of raindrop spectra. Part II: Collisional collection/breakup and evaporation in a rainshaft.

,*J. Atmos. Sci.***46****,**3312–3327.Young, K. C., 1975: The evolution of drop spectra due to condensation, coalescence, and breakup.

,*J. Atmos. Sci.***32****,**965–973.Zhang, G., J. Vivekanandan, and E. Brandes, 2001: A method for estimating rain rate and drop size distribution from polarimetric radar measurements.

,*IEEE Trans. Geosci. Remote Sens.***39****,**830–841.Zhang, G., J. Sun, and E. A. Brandes, 2006: Improving parameterization of rain microphysics with disdrometer and radar observations.

,*J. Atmos. Sci.***63****,**1273–1290.Zrnić, D. S., and A. V. Ryzhkov, 1999: Polarimetry for weather surveillance radars.

,*Bull. Amer. Meteor. Soc.***80****,**389–406.Zrnić, D. S., and A. V. Ryzhkov, 2004: Polarimetric properties of chaff.

,*J. Atmos. Oceanic Technol.***21****,**1017–1024.

## APPENDIX

### Equations for the Thermodynamic Terms

*L*is the latent heat of vaporization;

_{υ}*R*is the gas constant for water vapor;

_{υ}*T*is the air temperature;

*ρ*is the density of liquid water;

_{l}*K*and

*D*are the thermal conductivity of air and diffusivity of water vapor, respectively;

_{υ}*e*(

_{s}*T*) is the saturation vapor pressure as a function of

*T*; and

*f*and

_{υ}*f*are the ventilation coefficients for vapor and heat,

_{h}*N*

_{Sc}), Reynolds (

*N*

_{Re}), and Prandtl (

*N*

_{Pr}) numbers. The functional dependence of

*L*,

_{υ}*K*,

*D*, and

_{υ}*e*(

_{s}*T*) on temperature follows Rasmussen and Heymsfield (1987) and is provided in SI units here for convenience. The latent enthalpy of vaporization (J kg

^{−1}) as a function of temperature (K) is given as

*γ*is

^{−1}s

^{−1}K

^{−1}) is expressed as

^{2}s

^{−1}) is expressed as

*p*

_{0}is the reference level pressure, taken as 1000 hPa in this study. The saturation vapor pressure as a function of temperature is approximated following Rogers and Yau (1989):

*A*= 2.53 × 10

^{9}hPa. To calculate the Schmidt, Prandtl, and Reynolds numbers, the kinematic viscosity of the air (

*υ*) is required. The kinematic viscosity is dependent on the dynamic viscosity of the air (

_{a}*η*), as well as the air density,

_{a}*ρ*:

_{a}^{−1}s

^{−1}) for

*T*> 273 K is assumed to be

DSD models used for the simulations: the MP exponential distribution (solid black line); the *μ* = −1 gamma model (Γ, *μ* = −1, solid dark gray line), the *μ* = 1 gamma model (Γ, *μ* = 1, solid light gray line), the gamma model with *μ* = 3 (Γ, *μ* = 3, dashed black line), and the model with *μ* = 5 (Γ, *μ* = 5, dashed dark gray line). All models are for rainfall rates of 5 mm h^{−1}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

DSD models used for the simulations: the MP exponential distribution (solid black line); the *μ* = −1 gamma model (Γ, *μ* = −1, solid dark gray line), the *μ* = 1 gamma model (Γ, *μ* = 1, solid light gray line), the gamma model with *μ* = 3 (Γ, *μ* = 3, dashed black line), and the model with *μ* = 5 (Γ, *μ* = 5, dashed dark gray line). All models are for rainfall rates of 5 mm h^{−1}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

DSD models used for the simulations: the MP exponential distribution (solid black line); the *μ* = −1 gamma model (Γ, *μ* = −1, solid dark gray line), the *μ* = 1 gamma model (Γ, *μ* = 1, solid light gray line), the gamma model with *μ* = 3 (Γ, *μ* = 3, dashed black line), and the model with *μ* = 5 (Γ, *μ* = 5, dashed dark gray line). All models are for rainfall rates of 5 mm h^{−1}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Example simulation in a well-mixed boundary layer of 3-km depth, with surface temperature of 30°C and surface relative humidity of 70%. The simulation uses the MP DSD aloft with a rainfall rate of 5 mm h^{−1}. The left column is after 60 s, showing the transient effect of differential sedimentation in vertical profiles of each of the polarimetric variables. The right column shows the steady-state profiles for the same environment. Solid lines are for S band; dashed lines indicate C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Example simulation in a well-mixed boundary layer of 3-km depth, with surface temperature of 30°C and surface relative humidity of 70%. The simulation uses the MP DSD aloft with a rainfall rate of 5 mm h^{−1}. The left column is after 60 s, showing the transient effect of differential sedimentation in vertical profiles of each of the polarimetric variables. The right column shows the steady-state profiles for the same environment. Solid lines are for S band; dashed lines indicate C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Example simulation in a well-mixed boundary layer of 3-km depth, with surface temperature of 30°C and surface relative humidity of 70%. The simulation uses the MP DSD aloft with a rainfall rate of 5 mm h^{−1}. The left column is after 60 s, showing the transient effect of differential sedimentation in vertical profiles of each of the polarimetric variables. The right column shows the steady-state profiles for the same environment. Solid lines are for S band; dashed lines indicate C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Sensitivity of the evaporative change in the S-band dual-polarization variables to various constant RH profiles for the five DSD models. The domain depth is 2 km and is isothermal at 20°C. The plotting convention for the DSD models is the same as in Fig. 2. Shown are evaporative changes in (a) *Z _{H}*, (b)

*Z*

_{DR}, (c)

*K*

_{DP}, and (d)

*ρ*

_{HV}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Sensitivity of the evaporative change in the S-band dual-polarization variables to various constant RH profiles for the five DSD models. The domain depth is 2 km and is isothermal at 20°C. The plotting convention for the DSD models is the same as in Fig. 2. Shown are evaporative changes in (a) *Z _{H}*, (b)

*Z*

_{DR}, (c)

*K*

_{DP}, and (d)

*ρ*

_{HV}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Sensitivity of the evaporative change in the S-band dual-polarization variables to various constant RH profiles for the five DSD models. The domain depth is 2 km and is isothermal at 20°C. The plotting convention for the DSD models is the same as in Fig. 2. Shown are evaporative changes in (a) *Z _{H}*, (b)

*Z*

_{DR}, (c)

*K*

_{DP}, and (d)

*ρ*

_{HV}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

As in Fig. 4, but for C band.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

As in Fig. 4, but for C band.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

As in Fig. 4, but for C band.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Calculations of the evaporative changes in each of the polarimetric variables (each row is a variable) based on a 3-km-deep idealized well-mixed layer, with a surface temperature of 30°C. The surface RH values are shown along the abscissa. Each column represents a DSD model, as indicated. The solid lines represent S-band values, and the dashed lines represent C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Calculations of the evaporative changes in each of the polarimetric variables (each row is a variable) based on a 3-km-deep idealized well-mixed layer, with a surface temperature of 30°C. The surface RH values are shown along the abscissa. Each column represents a DSD model, as indicated. The solid lines represent S-band values, and the dashed lines represent C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Calculations of the evaporative changes in each of the polarimetric variables (each row is a variable) based on a 3-km-deep idealized well-mixed layer, with a surface temperature of 30°C. The surface RH values are shown along the abscissa. Each column represents a DSD model, as indicated. The solid lines represent S-band values, and the dashed lines represent C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Sensitivity of the evaporative change in the polarimetric variables to rainfall rate for the exponential DSD. The 3-km-deep domain has a dry-adiabatic lapse rate (surface temperature of 30°C) with a surface relative humidity of 75%. The S-band values are shown as solid lines and C-band values are shown as dashed lines. (left) The Δ*Z _{H}* (black) and Δ

*K*

_{DP}(gray) are shown, both in decibels. (right) Black represents Δ

*Z*

_{DR}and gray represents Δ

*ρ*

_{HV}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Sensitivity of the evaporative change in the polarimetric variables to rainfall rate for the exponential DSD. The 3-km-deep domain has a dry-adiabatic lapse rate (surface temperature of 30°C) with a surface relative humidity of 75%. The S-band values are shown as solid lines and C-band values are shown as dashed lines. (left) The Δ*Z _{H}* (black) and Δ

*K*

_{DP}(gray) are shown, both in decibels. (right) Black represents Δ

*Z*

_{DR}and gray represents Δ

*ρ*

_{HV}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Sensitivity of the evaporative change in the polarimetric variables to rainfall rate for the exponential DSD. The 3-km-deep domain has a dry-adiabatic lapse rate (surface temperature of 30°C) with a surface relative humidity of 75%. The S-band values are shown as solid lines and C-band values are shown as dashed lines. (left) The Δ*Z _{H}* (black) and Δ

*K*

_{DP}(gray) are shown, both in decibels. (right) Black represents Δ

*Z*

_{DR}and gray represents Δ

*ρ*

_{HV}.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Observed low-level (surface to approximately 3 km AGL) soundings used for the evaporation model simulations, representing three different evaporation scenarios: (a) Albuquerque, NM, at 0000 UTC 16 Aug 2007; (b) Norman, OK, at 1800 UTC 24 May 2008; and (c) Wallops Island, VA, at 0000 UTC 7 Oct 2009. Thick black solid (dashed) line indicates the temperature (dewpoint temperature) profile. The horizontal bar at the top of each sounding indicates the 3-km AGL height.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Observed low-level (surface to approximately 3 km AGL) soundings used for the evaporation model simulations, representing three different evaporation scenarios: (a) Albuquerque, NM, at 0000 UTC 16 Aug 2007; (b) Norman, OK, at 1800 UTC 24 May 2008; and (c) Wallops Island, VA, at 0000 UTC 7 Oct 2009. Thick black solid (dashed) line indicates the temperature (dewpoint temperature) profile. The horizontal bar at the top of each sounding indicates the 3-km AGL height.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Observed low-level (surface to approximately 3 km AGL) soundings used for the evaporation model simulations, representing three different evaporation scenarios: (a) Albuquerque, NM, at 0000 UTC 16 Aug 2007; (b) Norman, OK, at 1800 UTC 24 May 2008; and (c) Wallops Island, VA, at 0000 UTC 7 Oct 2009. Thick black solid (dashed) line indicates the temperature (dewpoint temperature) profile. The horizontal bar at the top of each sounding indicates the 3-km AGL height.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Simulated vertical profiles of polarimetric variables for the Albuquerque sounding (Fig. 8a). The solid line indicates S band and the dashed line indicates C band. Variables shown are (top) *Z _{H}*, (top middle)

*Z*

_{DR}, (bottom middle)

*K*

_{DP}, and (bottom)

*ρ*

_{HV}. The MP DSD model is used at the cloud base, with a rainfall rate of 5 mm h

^{−1}aloft.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Simulated vertical profiles of polarimetric variables for the Albuquerque sounding (Fig. 8a). The solid line indicates S band and the dashed line indicates C band. Variables shown are (top) *Z _{H}*, (top middle)

*Z*

_{DR}, (bottom middle)

*K*

_{DP}, and (bottom)

*ρ*

_{HV}. The MP DSD model is used at the cloud base, with a rainfall rate of 5 mm h

^{−1}aloft.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Simulated vertical profiles of polarimetric variables for the Albuquerque sounding (Fig. 8a). The solid line indicates S band and the dashed line indicates C band. Variables shown are (top) *Z _{H}*, (top middle)

*Z*

_{DR}, (bottom middle)

*K*

_{DP}, and (bottom)

*ρ*

_{HV}. The MP DSD model is used at the cloud base, with a rainfall rate of 5 mm h

^{−1}aloft.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Vertical profiles of rainfall rates for the three soundings. The thick solid line is the actual *R* determined from the DSD at each level. The dash–dot line is *R*(*Z _{H}*) [(23)], the dashed line is

*R*(

*K*

_{DP}) [(24)], and the dotted line is

*R*(

*Z*,

_{H}*Z*

_{DR}) [(25)]. At cloud base, the MP model is used for these simulations, with a rainfall rate of 5 mm h

^{−1}aloft.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Vertical profiles of rainfall rates for the three soundings. The thick solid line is the actual *R* determined from the DSD at each level. The dash–dot line is *R*(*Z _{H}*) [(23)], the dashed line is

*R*(

*K*

_{DP}) [(24)], and the dotted line is

*R*(

*Z*,

_{H}*Z*

_{DR}) [(25)]. At cloud base, the MP model is used for these simulations, with a rainfall rate of 5 mm h

^{−1}aloft.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Vertical profiles of rainfall rates for the three soundings. The thick solid line is the actual *R* determined from the DSD at each level. The dash–dot line is *R*(*Z _{H}*) [(23)], the dashed line is

*R*(

*K*

_{DP}) [(24)], and the dotted line is

*R*(

*Z*,

_{H}*Z*

_{DR}) [(25)]. At cloud base, the MP model is used for these simulations, with a rainfall rate of 5 mm h

^{−1}aloft.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

The dependence of the relative change in rainfall rate Δ*R*/*R* on *Z*_{DR} aloft and surface RH. Results are from calculations using a 3-km-deep idealized well-mixed layer with a surface temperature of 30°C. Each curve corresponds to a different value of the surface RH (solid black line for 55%, solid gray line for 75%, and dashed gray line for 95%). Only constrained gamma models (with 5 mm h^{−1} rainfall rate) are used in the calculations, where Δ*R*/*R* is independent of *Z _{H}* aloft. The gray numbers along the abscissa correspond to the shape parameter (

*μ*) values associated with the given

*Z*

_{DR}. (left) The S-band values, and (right) the C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

The dependence of the relative change in rainfall rate Δ*R*/*R* on *Z*_{DR} aloft and surface RH. Results are from calculations using a 3-km-deep idealized well-mixed layer with a surface temperature of 30°C. Each curve corresponds to a different value of the surface RH (solid black line for 55%, solid gray line for 75%, and dashed gray line for 95%). Only constrained gamma models (with 5 mm h^{−1} rainfall rate) are used in the calculations, where Δ*R*/*R* is independent of *Z _{H}* aloft. The gray numbers along the abscissa correspond to the shape parameter (

*μ*) values associated with the given

*Z*

_{DR}. (left) The S-band values, and (right) the C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

The dependence of the relative change in rainfall rate Δ*R*/*R* on *Z*_{DR} aloft and surface RH. Results are from calculations using a 3-km-deep idealized well-mixed layer with a surface temperature of 30°C. Each curve corresponds to a different value of the surface RH (solid black line for 55%, solid gray line for 75%, and dashed gray line for 95%). Only constrained gamma models (with 5 mm h^{−1} rainfall rate) are used in the calculations, where Δ*R*/*R* is independent of *Z _{H}* aloft. The gray numbers along the abscissa correspond to the shape parameter (

*μ*) values associated with the given

*Z*

_{DR}. (left) The S-band values, and (right) the C-band values.

Citation: Journal of Applied Meteorology and Climatology 49, 6; 10.1175/2010JAMC2243.1

Results of simulated evaporative changes in the polarimetric variables at S and C bands from the three soundings in Fig. 8. Simulations used the MP distribution with a rainfall rate aloft of 5 mm h^{−1}.

As in Table 1 but using the Γ, *μ* = −1 and Γ, *μ* = 5 distributions, respectively, with the values resulting from each DSD model separated by a comma in the table. Magnitudes of Δ*ρ*_{HV} less than 10^{−4} are indicated as ∼0.

^{1}

In all cases, the changes in the polarimetric variables due to evaporation are not constant in height and depend on the environmental thermodynamic profiles as expressed by the variable *ξ*(*z*).