Photo by Bill Wetzel, Jr.
The Introduction page contains a discussion of the motivation
for the ingredients-based approach and background information including: a history of
the ingredients-based methodology, our definition of an ingredient and an explanation of
the choice of ingredients, an examination of
traditional empirical techniques from an ingredients perspective, and a discussion of
the advantages of such an approach.
We begin on this page
with the list of the five key ingredients for winter season precipitation,
The remainder of this page discusses one set of the many possible diagnostics that could be used to
assess the presence and strength of each ingredient.
The IM presented here does not require employment of these diagnostics as the only
means of assessing the ingredients. In fact, at the heart of the IM is its flexibility to easily incorporatenew means of diagnosing the ingredients. Many different sets of diagnostics can be used, without compromising
the utility of the IM, as long as those diagnostics serve to assess the presence and strength of the
ingredients. We invite extensions and improvements to the diagnostics employed for each ingredient,
recognizing that any choice comes with limitations and that any one set of diagnostics will not be
suitable for all forecasters in all regions. To date, several NWS WFOs have taken our ingredients scripts and tailored them to their own needs. We encourage other users to
do the same.
- Forcing for Ascent, Where and how strong is the forcing?
- Instability (conditional and potential, upright and slantwise), How strong will the response be to a forcing?
- Moisture, Where and how much moisture is available?
- Efficiency, How much of the precipitation potential will be realized?
- Temperature, What form will the precipitation take and what snow to water
ratio can be expected?
1. Forcing for Ascent
As an initially unsaturated column of air is lifted by some forcing mechanism, its temperature decreases until saturation
occurs. Provided sufficient condensation or ice nuclei are present, any additional lifting after saturation is achieved results in
the condensation of water vapor to liquid or ice which can eventually
fall out as precipitation. This section addresses the forcing ingredient
and presents diagnostics for quantifying the degree of quasi-geostrophic (QG)
Note that reliance on QG diagnostics as the sole method of evaluating
the forcing for ascent ingredient is limiting. Use of this diagnostic
excludes explicit consideration of the effects of agesostrophy on the
vertical motion forcing. We chose to employ the Q-vector diagnostic
because it offers the distinct advantage of simultaneously capturing a
considerable amount of the synoptic-scale forcing as well as some frontal
forcing in a simple and familiar diagnostic.
The non-QG diagnostic, frontogenesis (calculated using the full horizontal
wind) has now been added to the ingredients maps. A discussion of this quantity has not been
added yet to this section.
A useful reference that explores a number of
vertical motion forcing mechanisms, including orographic forcing,
can be found at the
Centre Meteo Montreal's DIONYSOS.
The adiabatic and frictionless quasi-geostrophic
omega equation allows for the instantaneous diagnosis of vertical velocity
based on geopotential and temperature. This equation can be written
| Response (Term 1)
|| Forcing (Term 2)
|| Forcing (Term 3) |
| Laplacian of Vertical Motion
|| Differential Horizontal Temperature Advection
|| Rate of Change of Horizontal Vorticity Advection with Height |
- Term 1 is similar to a 3-D laplacian of omega, ,
except that this pseudo-Laplacian operator is modulated by a stability
, and the Coriolis parameter f.
By analogy to the behavior of a 3-D laplacian, the following relationships
describe the response () to a given forcing (Terms 2 and 3):
where UVM represents upward vertical motion and DVM downward vertical motion.
- Terms 2 and 3 represent the quasi-geostrophic forcing for vertical
motion. These terms have traditionally been considered as the separate
physical processes of the
Laplacian of horizontal temperature
advection and the rate of change of vorticity
with height (e.g., WAA, PVA, respectively).
For the case of warm air advection (WAA), . Thus (considering
only term 2) > 0 and < 0,
and term 2 contributes to forcing for UVM.
For the case of positive vorticity advection aloft (PVA), vorticity advection increases
with height . Thus (considering only term 3), > 0 and < 0, and term 3
contributes to forcing for
upward vertical motion.
When this form of the QG-omega equation is used to forecast forcing for vertical motions,
terms 2 and 3 cannot always be considered independently, as summarized below from
- Differential vorticity advection and temperature advection can act in
opposition to each other in the atmosphere, so a qualitative evaluation of the
two terms might result in one positive and one negative contribution, making an
estimation of the total forcing difficult.
- Vorticity advection and temperature advection are not independent of
each other because once the geopotential height field is specified, the 3-D vorticity
field and temperature advection are also specified.
- Terms 2 and 3 are not Galilean invariant (ie. they depend on the reference
frame in which they are being measured). If these forcing functions represent
fundamental physical processes, then observations of these processes should be
independent of the observer's reference frame.
These concerns can be remedied with the use of the Q-vector form of the quasi-geostrophic
omega equation (Hoskins et al., 1978),
where the Q-vector is defined as
This is an exact form of the QG omega equation (neglecting friction and diabatic effects)
and it includes contributions due to terms 2 and 3 from the form shown earlier.
A quantitative estimate of the vertical velocity would require a formal solution of
the QG-omega equation with boundary conditions on omega. However,
often a qualitative estimate of forcing for vertical motion is all that is required
in an operational
environment. Recognizing that the forcing () and the response () are
related by the pseudo-Laplacian operator, the Q-vector form of the QG-omega equation
provides a means of making a qualitative estimate:
Convergence of the Q-vector,
corresponds to a negative and thus forcing for upward vertical motion. Divergence
of the Q-vector corresponds to forcing for downward vertical motion.
Thus, isobaric contours of computed at a number of levels,
in conjunction with a cross-sectional analysis, provide a qualitative estimate of the direction and
magnitude of the vertical motion forcing throughout the atmosphere.
It is important to remember that the use of the Q-vector as a diagnostic for vertical
motion forcing does impose limits on the analysis, namely, that it is only
applied in those regions where the atmosphere is well described by the QG
equations (i.e., where the Rossby Number is low). No accommodation for the role of
non-QG motions in modulating the observed vertical motions is made in the
approach presented here. Such effects may be particularly
significant in the vicinity of frontal zones or jet streaks. A useful addition
to the present analysis would include a quantification of the degree to which
each case agrees with its QG description.
Furthermore, this approach to estimating the QG forcing neglects any modulation by the stability parameter,
because we approximate the pseudo-Laplacian operator
with the true Laplacian . Since appears in the full form of the
QG- equation, accurate estimates of the magnitude
of the upward vertical motions resulting from a given Q-vector convergence should
include consideration of the atmospheric stability. For this reason we also
include an instability ingredient, discussed in more detail in the next section.
For a discussion of the types of moist instability (CI, CSI, PI, PSI) and
techniques for diagnosing moist instability, including PVes, refer to pp. 17-23 of
Wetzel (2000). This is available in postscript or pdf form
through the main download page, or here
as gif images:
[good reference: David M. Schultz (NOAA/National Severe Storms Laboratory,
Norman, Oklahoma) and Philip N. Schumacher (NOAA/National Weather Service,
Grand Forks, North Dakota), The Use and Abuse of Conditional Symmetric
Instability (CSI) As a Diagnostic Tool.
In the Forcing Ingredient section, was discussed
as a technique for identifying regions of forcing for ascent.
However, the QG-omega equation reveals that an analysis of alone
will not account for the influence of atmospheric stability () on the
magnitude of the response () to a forcing. Thus, in the ingredients
approach, the instability ingredient is considered whenever the Q-vector diagnostic
is employed to account for the modulation of the response to a given forcing
The stability parameter considers only gravitational instabilities and
doesn't identify regions of symmetric instabilities. In order to allow for
consideration of both types, we use PVes to diagnose the instability ingredient,
where PVes is defined as,
By choosing to use PVes
rather than PVe, only conditional gravitational (CI)
and conditional symmetric instabilities (CSI) are considered, rather than
potential gravitational (PI) and potential symmetric instabilities (PSI).
Where PVes is negative, CI exists if
and CSI exists if
the atmosphere is characterized by two-dimensional flow and
CI and CSI are measures of the susceptibility of the atmosphere
to moist gravitational and moist symmetric convection, not indicators
of the existence of such convection (Schultz and Schumacher, 1999).
Thus, for either instability to be realized and convection to occur,
saturation must be present locally and a mechanism to force an
infinitesimal upward vertical motion must exist.
It is important to distinguish between gravitational and symmetric instabilities only insofar as the atmosphere
may respond differently, particularly with respect to the organization of precipitation bands
in each type of instability. However, this distinction is not emphasized in the forecast technique
because both instability mechanisms have similar implications, namely increased snowfall
amounts and the potential for lightning and thunder. Instead, the ingredients approach
will focus on identifying a region of
instability, then leave it to forecasters to analyze cross-sections and decide if a vertical
or slantwise response can be expected.
Given ample moisture, even very weak forcing can be sufficient to generate precipitation.
However precipitation will not be produced in a dry atmosphere even in the presence
of strong forcing and instability. Thus, an evaluation of available moisture should be
included in any forecast for precipitation.
There are a variety of ways to assess the moisture availability in a system.
The Magic Chart (Sangster, 1986; Chaston, 1989) leaves the evaluation of the
moisture availability up to the forecaster, recommending an examination of the
500-1000 hPa layer-averaged relative humidity.
The Garcia Method (Garcia, 1994) is based on the assumption that moisture is often the
limiting factor in a snow event. It predicts a direct relationship between the mixing
ratio on a mid-tropospheric isentropic surface
and snowfall amounts. For a 12-hour period, in a region where snowfall is expected based
on an evaluation of forcing mechanisms, the Garcia Method predicts
that the maximum amount of snow (in inches) will equal twice the maximum mixing ratio
(in g/kg) on an isentropic surface which intersects the 700-750 hPa layer over the
forecast area. The Cook Method (Cook, 1980) draws from the
theory of Jacobson et al. (1956) which holds that it isn't necessary to explicitly
worry about moisture in a developing storm because if the dynamics are strong moisture will
find its way into the storm.
The technique for evaluating moisture availability in the ingredient technique
presented here first involves the inspection of
relative humidity at a number of levels throughout the lower and mid-troposphere to
determine the degree of saturation. Relative humidity is defined as the ratio of the
actual mixing ratio (l) to the saturation mixing ratio with respect to water at
the same temperature and pressure (ls) and is typically expressed as a percent
(Wallace and Hobbs, 1977):
After determining the degree of saturation, the available moisture must be quantified
because relative humidity does not provide information about the absolute moisture
content and, thus,
cannot describe the amount or intensity of the resulting precipitation. For this purpose,
we examine the mixing ratio l, defined as the ratio of the mass of
water vapor to the mass of dry air in a parcel:
Mixing ratio was chosen for this technique because of our intent to integrate the ingredients
technique with the Garcia Method for estimating snowfall accumulation. The Garcia Method
is well suited for this purpose due to its comprehensive treatment of the moisture ingredient,
including an analysis of mixing ratios on isentropic surfaces.
However, isentropic analysis is not appropriate for neutral or unstable conditions, nor
regions where diabatic heating or cooling is occurring. This excludes any areas of instability and
regions characterized by a phase change. By integrating the Garcia Method
into the ingredients
technique, we can accommodate situations when isentropic analysis cannot be accurately or
prudently applied. The use of the Garcia Method in conjunction with the ingredients
approach is discussed more under Step 5: Estimating Snowfall Accumulation on the
Ingredients Analysis Tool page.
In this section the cloud microphysical
processes that influence the efficiency
of a winter precipitation event are discussed
(For a good reference on the cloud microphysics, see the
Wintertime Cloud Microphysics Review web site).
First, the processes involved in the formation of ice from supercooled liquid
droplets are examined. Then, the growth of ice particles following initiation,
and the conditions that lead to maximum growth are discussed. Finally, ways
to incorporate precipitation efficiency into the ingredients-based approach
To initiate the precipitation process in cold clouds, supercooled liquid must
freeze into an ice crystal. This can occur in three ways: heterogeneous nucleation,
homogeneous nucleation, or deposition onto a pre-existing ice crystal. Heterogeneous
nucleation is the process by which
supercooled liquid condenses onto an ice nuclei to form an ice crystal.
In homogeneous nucleation,
supercooled liquid can form an ice crystal without any nuclei present, however this requires a
temperature below -40 C. Homogeneous nucleation is not considered here because
it is not typically characteristic of winter storms at mid-latitudes.
Furthermore, it is assumed that there are no ice crystals
introduced into the layer of supercooled liquid, thus limiting our discussion to the initiation
of ice crystals by heterogeneous nucleation.
The initiation of an ice crystal by heterogeneous nucleation requires a small particle, or nuclei.
The ability of supercooled liquid to condense on the nuclei and form an ice crystal depends on
the chemical makeup of the particle, and the temperature and relative humidity of the ambient air
(Rogers and Yau, 1989). Over the upper Midwest, 80-90% of the ice nuclei (IN) are some form of clay,
primarily vermiculite. Vermiculite requires an air temperature of -15 C for ice nucleation to
occur. The table below shows some typical IN and their corresponding
| IN Composition || Activation Temperature |
| Silver Iodide || -4 C|
|Kaolinite || -9 C|
|Volcanic Ash ||-13 C|
|Vermiculite || -15 C|
Table 1: Examples of some typical ice nuclei (IN) and the temperature required for ice crystal
activation on their surfaces, from Rogers and Yau (1989)
Because of the presence of multiple types of IN in a cloud, and the variation of initiation
temperature with ambient relative humidity and surface characteristics of the
IN, a specific value cannot be assigned to the temperature at which ice
nucleation is always initiated. Instead,
guidelines determined by statistical studies must be used to identify a range of
temperatures for which ice nucleation can be expected.
Based on the results of a number of studies (e.g., Borovikov et al., 1963;
Mossop et al., 1970), Baumgardt (1999) prepared guidelines for
determining if ice nuclei might be initiated within in a cloud. These are shown in Table 2.
|Cloud Temperature || Chance of Ice in the Cloud |
|less than -4 C || 0% (No Ice) |
|-10 C || 60% |
|-12 C || 70% |
|-15 C || 90% |
|-20 C || 100% |
Table 2:Guidelines for determining if ice nucleation has
occurred in a cloud, adapted from Baumgardt (1999).
Furthermore, Baumgardt (1999) concluded that -12C to -14C should be considered the
minimum range for a good chance of ice being in a cloud, and offered -10C as the operational
cutoff point for no ice.
Ice Particle Growth
After initiation, ice particles can grow in the cloud through the following mechanisms:
deposition, accretion, or aggregation. Growth by deposition involves the change from water
vapor to ice when water vapor is deposited on ice or ice nuclei.
This is governed by the Bergeron-Findeisen
process in which ice crystals grow at the expense of liquid
droplets in environments for which the relative humidity is 100%.
Since the saturation vapor pressure over ice is less
than that over water, the vapor is compelled toward the ice, depositing on the
ice when it makes contact. The growth rate by deposition is greater at colder
temperatures, with a maximum growth rate around -15 C (Rogers and Yau, 1989, see
Figure 1: Normalized ice crystal growth rate as a function of temperature, for
50 kPa and 100kPa. Adapted from Rogers and Yau (1989).
Growth by accretion occurs when an ice particle overtakes or captures
supercooled liquid droplets. This occurs mainly after the ice particle has grown
to a sufficient size to begin to fall and collect the supercooled droplets.
Hence, ice must initially grow by deposition at the top or mid-level of the cloud
and only later by accretion as it falls through the cloud.
Aggregation is the joining of multiple ice particles to form one large snowflake. This
process generally occurs when temperatures are warmer than -10 C, with the largest
snowflakes occurring when the temperature is near 0 C.
In the ingredients-based approach to forecasting winter season precipitation
presented here, temperature is included at each level in an analysis of the ingredients in order
to assess the likelihood of ice nucleation and to anticipate enhanced precipitation
rates associated with rapid snowflake growth. The minimum cloud-top temperature component
of the efficiency analysis is best performed using satellite-derived cloud-top temperature data
and inspection of soundings during the event; however, forecast soundings and temperature
on isobaric surfaces can provide an estimate of the likelihood that ice will be
initiated in an event. Snowflake growth can be assessed in the ingredients approach using
model-predicted temperature by evaluating whether an area of strong forcing for vertical motion
coincides with a region of maximum depositional growth (T approximately -15 C).
In this situation, heavy snowfall is likely.
An analysis of the temperature ingredient in a winter season mid-latitude
cyclone addresses the precipitation type, or whether the precipitation
will fall as rain, snow, ice pellets, or freezing rain.
The temperature ingredient is best diagnosed through the careful monitoring of observed temperature profiles as the storm develops; however, isobaric analyses of model-predicted temperature and moisture fields can help forecasters to anticipate precipitation type well in advance of the event.
As with other forecasting rules of thumb, the conventional guidelines for predicting precipitation type are best applied in the context of a physical understanding of the thermodynamics governing precipitation type.
In a winter season weather system, precipitation is usually formed initially as
snow and the thermal structure of the layer through which it falls determines
the type of precipitation that will reach the ground. The initiation and growth
of the original ice crystals requires special consideration and was the subject
of the efficiency ingredient section. The discussion begins here with ice crystals that have grown large enough to fall from the cloud.
The amount of heat transferred to a melting ice crystal is proportional to (1) the difference between the temperature of the crystal (generally 0 C) and the surrounding air, and (2) the amount of time the crystal spends at that temperature. Therefore, the ability of an ice crystal to melt depends on its fall velocity (related to its size and shape), crystal structure, and the temperature of the air through which it falls. If the air temperature remains below freezing between the level of ice crystal initiation and the ground, the precipitation will
fall to the ground as snow. However, if a layer of warm (above freezing) air exists, either elevated or extending to the surface, there is the potential for one or more phase changes to occur before the precipitation
reaches the ground.
Monotonically Decreasing Temperature Profile
For an atmospheric temperature profile which decreases monotonically with height, with an above-freezing layer extending to the surface, the precipitation will
remain frozen and fall as snow if the wet bulb temperature is below 0 C throughout the column. This is because ice (liquid) begins to sublimate (evaporate) if it falls through an unsaturated layer. The environmental temperature will cool due to this sublimation (evaporation) and the frost point (dew point) will rise, so that both approach the wet-bulb temperature. For the
case of a surface
warm layer wet-bulb temperature greater than 0 C, there is the potential for partial or
complete melting of the ice crystal. The probability that the precipitation will be snow
has previously been expressed simply in terms of the depth of the warm layer (or, the height of the
freezing level above ground). Table 3 presents the work of McNulty (1988) describing
the likelihood as a function of warm layer depth that the precipitation type will be snow. Notice that the temperature of the warm layer is not considered here.
|Freezing Level|| Chance of Snow |
|35mb or 290 m || 50 % |
|25mb or 201 m || 70 % |
|12mb or 96 m || 90 % |
Table 3: Relationship between the freezing level for a temperature profile decreasing monotonically with height and the likelihood that the precipitation
type will be snow, adapted from McNulty (1988).
McNulty (1988) also notes, On the average, the freezing level must be at least 1200
feet above the surface to insure that the snow will melt to rain. This
statement can be extended by recognizing that in the Western Plains regions of the United
States, the 850 hPa surface is approximately 1200 feet above ground. Therefore, the 0 C 850 hPa
isotherm marks approximately the rain edge of the rain-snow boundary in these locations. Areas with an 850 hPa
temperature greater than 0 C are nearly certain to experience rain, however
a close inspection
of McNulty's statement reveals that it does not comment on the likelihood of precipitation falling
as snow for temperatures less than 0 C. Moreover, Table 3
even for a freezing
level at 96 m (or approximately 12 hPa) above ground, there is still a 10% chance that the precipitation will not be snow.
Therefore, according to this interpretation of McNulty's work, the conventional
use of the
850 hPa 0 C contour as the rain-snow line must be
applied with caution, particularly in the snow regime.
Temperature Profile with an Elevated Warm Layer
In the event of an elevated layer of above-freezing temperature with below-freezing temperatures underneath it (Figure 2)
the potential exists for freezing rain or sleet/ice pellets.
Figure 2: Schematic of an Elevated Warm Layer. Thick solid line is
ambient air temperature, thin solid lines are -10 C, 0 C, and 10 C isotherms.
Adapted from Czys et al. (1996).
The degree of melting that occurs in the warm layer is the key difference between conditions
characteristic of freezing rain and those of sleet (Penn, 1957).
If the snow completely melts to a liquid in the warm layer, then freezing rain or rain is likely.
However, if the snow is only partially melted when it reaches the cold air beneath the warm layer,
it begins to re-freeze immediately due to the presence of the remaining ice nuclei and falls to the
ground as sleet or snow.
Czys et al. (1996) present a means of determining the degree of melting in an elevated warm
layer based on the temperature of the elevated warm layer and the depth of that
The authors introduce a nondimensional parameter to assist in the determination of whether
complete melting can be expected in an elevated warm layer and thus, in the diagnosis of
freezing rain and ice pellets. In their study, is defined for a particle of a given radius
as the ratio between the residence time of the particle in the warm
layer and the
time required for complete melting . is based on the vertical
distance the particle must fall (i.e., depth of the warm layer) and the velocity at which it falls.
The relationship describing is derived through a balance between the
rate of energy required
to transform the solid to a liquid and the rate that the environment can supply
the needed energy.
The resulting nondimensional parameter , can be used to determine the degree of melting
in an elevated warm layer for a given maximum particle radius.
For >= 1, conditions are
characteristic of complete melting and freezing rain is favored. For < 1, only partial
melting of the largest particles is expected and precipitation will likely fall
as ice pellets or snow.
Figure 3 shows isopleths of the threshold value, = 1, for a range of ice particle
radii, warm layer depths, and mean warm layer temperatures.
Figure 3: Isonomogram of = 1 for different ice particle
radii (microns) computed
over a range of warm layer depths and mean temperatures. Adapted from Czys et al. (1996).
Based on 17 case studies in the Midwestern winter season of 1995-1996, Czys et al. (1996) found that a critical
ice particle radius of 400 microns provides a reasonable level of accuracy.
Thus, if the mean warm layer temperature and warm layer depth intersect in Figure 3 above
the 400 micron = 1 isopleth, complete melting and freezing rain can be expected provided the
surface temperature is below 0 C. If these parameters intersect below the 400 micron = 1 isopleth,
only partial melting is predicted.
An alternative to the Czys et al. (1996) technique considers only the
warm layer maximum temperature. Because the depth of the warm layer is generally proportional to the
maximum warm layer temperature, this simplification has proven to be a good diagnostic to
determine the degree of melting (Baumgardt, 1999). Stewart and King (1987) modeled the melting of
snow as it falls through an elevated warm layer based on the empirical rates of
snowflake melting described
by Stewart (1985). The authors ran this model for various snowflake sizes, densities, maximum warm layer
temperatures, and lapse rates. Evaporation and condensation were not considered, as the layer was
assumed just saturated. Figure 4 presents these results
and Table 4 summarizes the degree of melting as a function of maximum warm layer
Figure 4: The degree of snowflake melting as a function of maximum
warm layer temperature. The curves correspond to particle size expressed in terms of
the raindrop-equivalent diameter, and for two different initial snowflake
densities. Adapted from Stewart and King (1987).
Max. T of Warm Layer || Degree of Melting |
| > 3-4 C || Melts Completely |
| 1-3 C || Partial Melting |
| < 1 C || Very Little Melting with Re-freezing to Snow |
Table 4: Summary of results from Stewart and King (1987),
relating the degree of snowflake melting in an elevated warm layer to the
maximum temperature of the warm layer.
For the case of complete melting (corresponding to a warm layer
temperature > 3-4 C),
it is also important to consider the temperature of the cold layer beneath the warm layer (McNulty, 1988).
If the minimum temperature of the cold layer is greater than -10 C, it is unlikely that
the freezing process can be initiated so the precipitation will remain as a liquid until it
reaches the ground. Provided the ground temperature is below freezing, this will result in
freezing rain. When the minimum temperature of the cold layer is less than
-10 C, it is more likely that freezing will be initiated and snow or sleet will result
provided there are enough ice nuclei. The depth of the cold air will determine
degree of re-freezing and, thus, whether the precipitation will be snow or sleet.
The precipitation type analyses described above are best performed with the use
of modeled or observed
soundings. Because observed sounding data is usually only available at 12-hour
intervals, it may be useful to infer changes in the lower tropospheric thermal structure between
balloon launches. Penn (1957) discussed the influences on the evolution of a temperature
profile: thermal advection, adiabatic temperature change associated with vertical displacement,
and non-adiabatic heating or cooling. Thermal advection by the horizontal wind
is generally the
dominant process, however this can act in opposition to an adiabatic temperature change.
Because the horizontal advection of warm air is generally associated with upward vertical motion,
warming associated with the advection can be reduced by adiabatic cooling associated with
the rising air parcels. The net effect can be a reduction in the rate of warming by advection
alone by roughly 50%. Additionally, diabatic processes can influence the temperature profile. Evaporational
cooling when precipitation falls through a dry layer can reduce the temperature
by 5-10 C/hour
until the air is saturated and the temperature reaches the wet bulb temperature. In very heavy
snow events, the freezing level may also be altered by the absorption of latent
heat when snow
is melted to rain in a warm layer.
Figure 5 presents a summary of the steps involved in the diagnosis of precipitation
Figure 5: Summary of the steps involved in a diagnosis of precipitation type.
This figure illustrates the complexity of the physical processes involved in determining precipitation type.
Given forecasted and observed soundings, the ideas presented above will allow
forecasters to perform a reasonable assessment of the
likelihood that precipitation will be snow, ice pellets, freezing rain, or rain. However, in order to
include a determination of precipitation type on isobaric surfaces with
the other ingredients, the analysis is limited to diagnostics that can be computed
at one level or layer in the atmosphere. Thus,
in the ingredients maps, the 850 hPa 0 C rain-snow line will be
used as a rough guideline for a
determining precipitation type. This guideline is only reliable when the forecast region is more than a few hundred kilometers to either side of the 0 C isotherm and more than a few degrees cooler than 0 C, and will be used
in conjunction with an analysis of forecasted and observed soundings when the forecast area is any closer or any warmer.
If ice or liquid falls through a dry (unsaturated) layer, it begins to evaporate.
The environmental temperature will cool due to this evaporation and the
dew point will rise, both approaching the wet-bulb temperature. This effect
should be considered when predicting the type of precipitation.
- If the wet-bulb temperature is below zero (Tw < 0) and the precipitation enters
in the form of ICE, expect no melting to occur and the precipitation
to fall as snow.
- If the wet-bulb temperature is below zero (Tw < 0) and the
precipitation enters as a LIQUID, it may freeze to form ice pellets.
However, warm advection below the clouds may cause a precipitation
change back to rain if the temperature is above 0C.
- In both of these situations, the issue of complete evaporation may also be addressed.
The values given below are derived in McIlveen (1989):
In an environment with 95% Relative Humidity,
Droplet Radius (micrometers)
Fall Distance Before Complete Evaporation|
Or, cloud droplets with radius less than 30 micrometers cannot survive
a fall of more than a few centimeters when the air is even mildly
subsaturated. Drizzle can fall a few hundred meters, and
rain could fall 240 km -- more than 10 times the depth
of the troposphere. For dryer conditions these distances will be
Some Rules of thumb for determining snow type:
- Cook (1980):
Precipitation will be snow or change to snow whenever
the wet-bulb temperature curve goes below freezing within
2500 feet of the surface and does not go above freezing at
a higher level.
50% probability that precipitation
will fall as snow if the 200mb heights are 11760m, with a higher
probability for lower heights. Note that this is for low elevations
with average surface pressures (1005-1015mb); for lower surface
pressures or in the immediate coastal sections the 50% probability
height is 11640m.
Rain/snow lines for low elevations above sea level
5400 m |
2560 m |
4100 m |
2840 m |
1300 m |
1540 m |
Return to Winter Forecast Ingredients Page