Weather radar

From Wikipedia, the free encyclopedia

Weather radar in Norman, Oklahoma with rainshaft (Source: NOAA)

A weather radar is a type of radar used to locate precipitation, calculate its motion, estimate its type (rain, snow, hail, etc.), and forecast its future position and intensity. Modern weather radars are mostly doppler radars, capable of detecting the motion of rain droplets in addition to intensity of the precipitation. Both types of data can be analyzed to determine the structure of storms and their potential to cause severe weather.

Environment Canada King City (CWKR) weather radar station.

1 History
2 How weather radar works
3 Data types
4 Main types of radar outputs
5 Limitations and artifacts
6 Solutions for now and the future
7 Specialized Applications
- 7.1 Avionics Weather Radar
8 References
9 Bibliography
10 See also
- 10.1 Related articles
- 10.2 External links

[edit] History

During World War II, military radar operators noticed noise in returned echoes due to weather elements like rain, snow, and sleet. Just after the war, military scientists returned to civilian life or continued in the Armed Forces and pursued their work in developing a use for those echoes. In the United States, David Atlas,^[1] for the Air Force group at first, and later for MIT, developed the first operational weather radars. In Canada, J.S. Marshall and R.H. Douglas formed the "Stormy Weather Group^[2]" in Montreal. Marshall and his doctoral student Walter Palmer are well known for their work on the drop size distribution in mid-latitude rain that led to understanding of the Z-R relation, which correlates a given radar reflectivity with the rate at which water is falling on the ground. In the United Kingdom, research continued to study the radar echo patterns and weather elements such as stratiform rain and convective clouds, and experiments were done to evaluate the potential of different wavelengths from 1 to 10 centimetres.

In 1953, Donald Staggs, an electrical engineer working for the Illinois State Water Survey, made the first recorded radar observation of a "hook echo" associated with a tornadic thunderstorm.

Between 1950 and 1980, reflectivity radars, which measure position and intensity of precipitation, were built by weather services around the world. The early meteorologists had to watch a cathode ray tube. During the 1970s, radars began to be standardized and organized into networks. The first devices to capture radar images were developed. The number of scanned angles was increased to get a three-dimensional view of the precipitation, so that horizontal cross-sections (CAPPI) and vertical ones could be performed. Studies of the organization of thunderstorms were then possible with the Alberta Hail Project in Canada and NSSL in the US in particular. NSSL, created in 1964, began experimentation on dual polarization signals and on Doppler effect uses.

Between 1980 and 2000, weather radar networks became the norm in North America, Europe, Japan and other developed countries. Conventional radars were replaced by Doppler radars, which in addition to position and intensity of could track the relative velocity of the particles in the air. In the United States, a network consisting of 10 cm wavelength radars, called NEXRAD or WSR-88D (Weather Service Radar 1988 Doppler), was started in 1988. In Canada, Environment Canada constructed the King City station,^[3] with a five centimeter research Doppler radar, by 1985;McGill University dopplerized its radar (CWMN) in 1993. This led to a complete Canadian Doppler network^[4] between 1998 and 2004. France and other European countries switched to Doppler network by the end of the 1990s to early 2000s. Meanwhile, rapid advances in computer technology led to algorithms to detect signs of severe weather and a plethora of "products" for media outlets and researchers.

After 2000, research on dual polarization technology has moved into operational use, increasing the amount of information available on precipitation type (e.g. rain vs. snow). "Dual polarization" means that microwave radiation which is polarized both horizontally and vertically (with respect to the ground) is emitted. Wide-scale deployment is expected by the end of the decade in some countries such as the United States, France,^[5] and Canada.

Since 2003, the U.S. National Oceanic and Atmospheric Administration has been experimenting with phased-array radar as a replacement for conventional parabolic antenna in order to provide more time resolution in atmospheric sounding. This would be very important in severe thunderstorms as their evolution can be better evaluated with more timely data.

[edit] How weather radar works

[edit] Sending radar pulses

A radar beam spreads out as it moves away from the radar station, covering an increasingly large volume.

Weather radars send directional pulses of microwave radiation, on the order of a microsecond long, using a cavity magnetron or klystron tube connected to a waveguide to a parabolic antenna. The wavelengths of 1 to 10 cm are approximately ten times the diameter of the droplets or ice particles of interest, because Rayleigh scattering occurs at these frequencies. This means that part of the energy of each pulse will bounce off these small particles, back in the direction of the radar station.

Shorter wavelengths are useful for smaller particles, but the signal is more quickly attenuated. Thus 10 cm (S-band) radar is preferred to but is more expensive than a 5 cm C-band system. 3 cm X-band radar is used only for very short distance purposes, and 1 cm Ka-band weather radar is used only for research on small-particle phenomena such as drizzle and fog.

Radar pulses spread out as they move away from the radar station. This means that the region of air any given pulse is moving through is larger for areas farther away from the station, and smaller for nearby areas, decreasing accuracy at high distances. At the end of a 150-200 km sounding range, the volume of air scanned by a single pulse might be on the order of a cubic kilometer.

The volume of air that a given pulse takes up at any point in time may be calculated by the formula $\, {v = h r^2 \theta^2}$ , where v is the volume enclosed by the pulse, h is pulse width (in e.g. meters, calculated from the duration in seconds of the pulse times the speed of light), r is the distance from the station that the pulse has already travelled (in e.g. meters), and $\,\theta$ is the beam width in radians).

[edit] Listening for return signals

Between each pulse, the radar station serves as a receiver and listens for return signals from particles in the air. The duration of the "listen" cycle is on the order of a millisecond, which is a thousand times longer than the pulse duration. The length of this phase is determined by the need for the microwave radiation (which travels at the speed of light) to propagate from the detector, to the weather target, and back again, a distance which could be several hundred kilometers. The horizontal distance from station to target is calculated simply from the amount of time that lapses from the initiation of the pulse to the detection of the return signal. (The time is converted into distance by multiplying by the speed of light). If pulses are emitted too frequently, the returns from one pulse will be confused with the returns from previous pulses, resulting in incorrect distance calculations.

[edit] Determining height

Echos height above ground

Typical angles scanned in Canada. The zigzag represent data angles used to make CAPPIs at 1.5 et 4 km of altitude

Assuming the Earth is round, with knowledge of the variation of the index of refraction through air and the distance to the target, one can calculate the height above ground of the target.

After each scanning rotation, the antenna elevation is changed for the next sounding. This scenario will be repeated on many angles in order to scan all the volume of air around the radar within the maximum range. Usually, this scanning strategy is completed within 5 to 10 minutes in order to have data within 15 km above ground and 250 km distance of the radar.

Due to the Earth curvature and change of index of refraction with height, the radar cannot “see” below the height above ground of the minimal angle or closer to the radar than the maximal one. This image shows the height of a series of typical angles done by a 5 cm weather radar in Canada. They range from 0.3 to 25 degrees.

[edit] Calibrating intensity of return

Because the targets are not unique in each volume, the radar equation has to be developed beyond the basic one:

$P_r = \left [P_t{{ G^2 \lambda^2 \sigma^0}\over{{(4\pi)}^3 R^4}} \right] \propto \frac {\sigma^0} {R^4}$

where $\,P_r$ is received power, $\,P_t$ is transmitted power, $\,G_t$ is the gain of the transmitting antenna, $\,\lambda$ is radar wavelength, $\,\sigma$ is the radar cross section of the target and $\,R$ is the distance from transmitter to target.

In this case, we have to add the cross sections of all the targets:

$\sigma^0 = \bar \sigma^0 = V \sum \sigma^0_j = V \eta$

$\begin{cases} V\quad= scanned\ volume \\ \qquad= pulse\ length\ \ X\ beam\ width \\ \qquad= \left[\frac {c\tau}{2} \right] \left[\frac {\pi R^2 \theta^2}{4} \right] \end{cases}$

where $\,c$ is the light speed, $\,\tau$ is temporal duration of a pulse and $\,\theta$ is the beam width in radians.

In combining the two equations :

$P_r = \left [P_t{{ G^2 \lambda^2 }\over{{(4\pi)}^3 R^4}} \right] \left[\frac {c\tau}{2} \right] \left[\frac {\pi R^2 \theta^2}{4} \right] \eta = \left [P_t \tau G^2 \lambda^2 \theta^2 \right] \left[\frac {c}{512(\pi^2)} \right] \frac {\eta} {R^2}$

Which leads to:

$P_r \propto \frac {\eta} {R^2}$

Notice that the return now varies inversely to $\, R^2$ instead of $\,R^4$ . In order to compare the data coming from different distances from the radar, one has to normalize them with this ratio.

[edit] Data types

[edit] Reflectivity (in decibel or dBZ)

Return echoes from targets, reflectivity, are analyzed for their intensities in order to establish the precipitation rate in the scanned volume. The wavelengths used (1 to 10 cm) ensure that this return is proportional to the rate because they are within the validity of Rayleigh scattering which states that the targets must be much smaller than the wavelength of the scanning wave (by a factor of 10).

Reflectivity perceived by the radar (Z_e) varies by the 6^th power of the rain droplets' diameter (D), the square of the dielectric constant (K) of the targets and the drop size distribution (e.g. N[D] of Marshall-Palmer) of the drops. This gives a truncated Gamma function,^[6] of the form:

$Z_e = \int_{0}^{Dmax}$ |K|² N₀e ^{$- Λ$ D} D⁶dD

Precipitation rate (R), on the other hand, is equal to the number of particles, their volume and their fall speed (v[D]) as:

R = $\int_{0}^{Dmax}$ N₀e^{$- Λ$ D} (

π

D³/6) v(D)dD

So Z_e and R have similar functions that can be resolved giving a relation between the two of the form:

Z = aR^b

Where a and b depend on the type of precipitations (snow, rain, convective or stratiform) which have different $Λ$ , K, N₀ and v.

As the antenna scan the atmosphere, on every angle of azimuth it obtains a certain number of return from each targets encountered. Reflectivity is then averaged for that target in order to have a better data set.

Since variation in diameter and dielectric constant of the targets can leads to large variability in power return to the radar, reflectivity is expressed in dBZ (10 times the logarithm of the ratio of the echo to a standard 1 mm diameter drop filling the same scanned volume).

[edit] How to read reflectivity on a radar display

Radar returns are usually described by colour or level. The colours in a radar image normally range from blue or green for weak returns, to red or magenta for very strong returns. The numbers in a verbal report increase with the severity of the returns.

For example, the U.S. National Doppler Radar sites use the following scale for different levels of reflectivity:[1]

magenta: 65 dBZ (extremely heavy precipitation)
red: 52 dBZ
yellow: 36 dBZ
green: 20 dBZ (light precipitation)

Strong returns (red or magenta) may indicate not only heavy rain but also thunderstorms, hail, strong winds, or tornadoes, but they need to be interpreted carefully, for reasons described later in this article.

[edit] Aviation conventions

When describing weather radar returns, pilots, dispatchers, and air traffic controllers will typically refer to three return levels:[2]

level 1 corresponds to a green radar return, indicating usually light precipitation and little to no turbulence, leading to a possibility of reduced visibility.
level 2 corresponds to a yellow radar return, indicating moderate precipitation, leading to the possibility of very low visibility, moderate turbulence and an uncomfortable ride for aircraft passengers.
level 3 corresponds to a red radar return, indicating heavy precipitation, leading to the possibility of thunderstorms and severe turbulence and serious structural damage to the aircraft.

Aircraft will try to avoid level 2 returns when possible, and will always avoid level 3 unless they are specially-designed research aircraft.

[edit] Velocity

Idealized example of Doppler output. Approaching velocities are in blue and receding velocities are in red. Notice the sinuosidal variation of speed when going around the display at a particular range. (Source: Environment Canada).

See also: Pulse-doppler radar and Doppler radar

[edit] Pulse pair

Any rain drops or snow flakes in motion affect the frequency of the returned radar beam according to the Doppler effect. With velocities of less than 70 m/s (150 miles/h) for weather echos and radar wavelength of 10 cm, it amounts to only 10^-5%. This difference is too small to be noted by electronic instruments. However, as the targets move slightly between each pulse, the returned wave has a noticeable phase difference or phase shift from pulse to pulse.

Doppler weather radars are using this phase difference (pulse pair difference) to calculate the precipitations motion. The intensity of the successively returning pulse from the same scanned volume where targets have slightly moved is :

$I = I_0 Sin \left(\frac{4\pi (x_0 + v \Delta t)}{\lambda}\right) = I_0 Sin \left(\Theta_0 + \Delta\Theta\right)\begin{cases} x = \mathrm{distance\ radar\ to\ target} \\ \lambda = \mathrm{radar\ wavelength} \\ \Delta t = \mathrm{time\ between\ two\ pulses} \end{cases}$

So $\Delta\Theta = \left(\frac{4\pi v \Delta t}{\lambda}\right)$

v = target speed = $\frac{\lambda\Delta\Theta}{4\pi \Delta t}$

This speed is called the radial Doppler velocity because it gives only the radial variation of distance versus time between the radar and the target. The real speed and direction of motion has to be extracted by the process described below.

[edit] Doppler dilemma

Maximum range from reflectivity (red) and unambiguous Doppler velocity range (blue) with pulse repetition rate.

The phase between pulse pairs can vary from - $π$ and + $π$ , so the unambiguous Doppler velocity range is

V_max = $\pm$ $\frac{\lambda}{4\Delta t}$

This is called the Nyquist velocity. This is inversely dependent on the time between successive pulses: the smaller the interval, the larger is the unambiguous velocity range. However, we know that the maximum range from reflectivity is directly proportional to $Δ t$ :

x = $\frac{c\Delta t}{2}$

The choice becomes increasing the range from reflectivity at the expense of velocity range, or increasing the latter at the expense of range from reflectivity. In general, the useful range compromise is 100 to 150 km for reflectivity. This means for a wavelength of 5 cm, like on the image, an unambiguous velocity range of 8.3 to 12.5 m/s but the double for a 10 cm radar like the NEXRAD.

Some techniques using two alternating pulse repetition rates (PRF) permit to extented the Doppler range. The velocities noted with the first pulse rate could be equal or different with the second. For instance, if the maximum velocity with a certain rate is 10 m/s and the one with the other rate is 15 m/s. The data coming from both will be the same up to 10 m/s and differ afterward. It is then possible to find a mathematical relation between the two returns and calculate the real velocity beyond the limitation of the two PRF.

[edit] Doppler interpretation

Radial component of real winds when scanning on 360 degrees.

If one thinks of an autumn rain uniformly filling the radar area coverage and moving from west to east, The radar beam pointing toward the west will “see” the rain drops moving toward it while looking east, it will notice them going away. On the other hand, looking north or south, since there is no motion toward the radar in those directions, the radial velocity is null.

As the beam is scanning 360 degrees around the radar, data will comes from all those angles and be the radial projection of the actual wind on the individual angle. The intensity pattern formed by this scan can be represented by a cosine curve, as seen on the right. One can then calculate the direction and the strength of the motion of particles as long as there is enough coverage on the radar screen.

However, the rain drops are falling. As the radar only sees the radial component and has a certain elevation from ground, the radial velocities are contaminated by some fraction of the falling speed. This component is negligible in small elevation angles, but must be taken into account for higher scanning angles.

[edit] Polarization

Targeting with dual-polarization will reveal the form of the droplet

Most liquid hydrometeors have a larger horizontal axis due to the drag coefficient of air while falling (water droplets). This causes the water molecule dipole to be oriented in that direction so radar beams are generally polarized horizontally to receive the maximal return.

If we decide to send simultaneously two pulses with orthogonal polarization: vertical and horizontal, we receive two sets of data proportional to the two axis of the droplets that are independent:^[7]

The difference between the intensities is called Zdr and gives information on the form of the target.
Electromagnetic waves change phase while passing through denser material, the phase differential with distance or specific phase differential ( $K d p$ ) can be used to estimate the amount of precipitation in the scanned volume of atmosphere: the rain rate. This measurement is not affected by attenuation.
The Zdr relation should be stable with drops of the same form, the return from a group of droplets of different forms or a mix of drops, snow flakes, hail, etc… continually changing position will have a Zdr that change with time. This variation ( $ρ$ _hv) will thus give an idea of the variety of forms in the scanned volume.

With this new knowledge added to the reflectivity and the Doppler data, researchers have been working on developing algorithms to differentiate precipitation types, non-meteorological targets, better accumulation estimates, etc. NCAR has been one of the world leaders in this field with Dusan S. Zrnic et Alexandre V. Ryzhkov.

NOAA has set up a test bed for operational radar since 2000 and plan to equip all its 10 cm wavelength NEXRAD with polarization by the end of the decade. McGill University J.S. Marshall Radar Observatory in Montreal, Canada has converted their instrument by 1999 and the data are used operationally by Environment Canada in Montreal. Another EC radar in King City (North of Toronto) has been polarized in 2005, this one work on a 5 cm wavelength which gives new challenges. EC hope to generalize this conversion to all its network eventually. Finally, Météo-France is working too on the subject and hope to set up their first polarized radars in 2008.

For more details:

US polarization development

McGill University operational output

[edit] Main types of radar outputs

All data from radar scans are displayed according to the need of the users. Different outputs have been developed through time to reach this. Here is a list of common and specialized outputs available.

[edit] Plan position indicator

Main article: Plan position indicator

Thunderstorm line viewed in reflectivity (dBZ) on a PPI (NOAA)

Since data are obtained one angle at a time, the first way of displaying them as been the Plan Position Indicator (PPI) which is only the layout of radar return on a two dimensional image. One has to remember that the data coming from different distances to the radar are at different heights above ground.

This is very important as a high rain rate seen near the radar is relatively close to what reach the ground but what is seen from 160 km (100 miles) away is about 1.5 km above ground and could be far different from the amount reaching the surface. It is thus difficult to compare weather echoes at different distance from the radar.

PPIs are afflicted with ground echoes near the radar as a supplemental problem. These can be misinterpreted as real echoes. So other products and further treatments of data have been developed to supplement its shortcomings.

USAGE: Reflectivity, Doppler and polarimetric data can use PPI.

N.B.: In the case of Doppler data, two points of view are possible: relative to the surface or the storm. When looking at the general motion of the rain to extract wind at different altitudes, it is better to use data relative to the radar. But when looking for rotation or wind shear under a thunderstorm, it is better to use the storm relative images that subtract the general motion of precipitation leaving the user to view the air motion as if he would be sitting on the cloud. Here are real time example: NWS Burlington radar, one can compare the BASE and STORM Doppler products

[edit] Constant Altitude Plan Position Indicator

Main article: Constant Altitude Plan Position Indicator

To avoid some of the problems on PPIs, the CAPPI or Constant Altitude Plan Position Indicator has been developed by researchers in Canada. It is basically a horizontal cross-section through radar data. This way, one can compare precipitation on an equal footing at difference distance from the radar and avoid ground echoes. Although data are taken at a certain height above ground, a relation can be inferred between ground stations reports and the radar data.

CAPPIs call for a large number of angles from near the horizontal to near the vertical of the radar in order to have a cut that is as close as possible at all distance to the height needed. But even then, after a certain distance, there isn’t any angle available and the CAPPI becomes the PPI of the lowest angle. The zigzag line on the angles diagram above shows the data used to produce a 1.5 and 4 km height CAPPIs. Notice that the section after 120 km is using the same data.

USAGE: Mostly for reflectivity data. McGill University is producing Doppler CAPPIs but the nature of velocity make the output a bit noisy as velocities can change rapidly in direction with height contrary to a relatively smooth pattern in reflectivity.

Real time examples:

[edit] Vertical Composite

Another solution to the PPI problems is to produce images of the maximum reflectivity in a layer above ground. This solution is usually taken when the number of angles available is small or variable. The American National Weather Service is using such Composite as their scanning scheme can vary from 4 to 14 angles, according to their need, which would make very coarse CAPPIs. The Composite make sure that no strong echo is missed in the layer and a treatment using Doppler velocities eliminate the ground echoes.

Real time example: NWS Burlington radar, one can compare the BASE and COMPOSITE products

[edit] Accumulations

24 hours rain accumulation on the Val d'Irène radar in Eastern Canada. Notice the zones without data in the East and Southwest caused by radar beam blocking from mountains. (Source: Environment Canada)

One of the main use of radar is to be able to assess the amount of precipitations fallen over large basins for hydrological purpose. For instance, river flood control, sewer management and dam construction are all areas where planners want accumulation data. It ideally completes surface stations data which they can use for calibration.

To produce radar accumulations, we have to estimate the rain rate over a point by the average value over that point between one PPI, or CAPPI, and the next; then multiply by the time between those images. If one wants for a longer period of time, one has to add up all the accumulations from images during that time.

[edit] Echotops

Aviation is a heavy user of radar data. One map particularly important in this field is the Echotops for flight planning and avoidance of dangerous weather. Most country weather radars are scanning enough angles to have a 3D set of data over the area of coverage. It is relatively easy to estimate the maximum altitude at which precipitation is found within the volume. However, those are not the tops of clouds as they extended to higher altitudes than the precipitation.

[edit] Vertical cross sections

To know the vertical structure of clouds, in particular thunderstorms or the level of the melting layer, a vertical cross sections product of the radar data is available to meteorologist.

[edit] Radar networks

Over the past few decades, radar networks have been extended to allow the production of composite views covering large areas. For instance, all major countries (e.g., the United States, Canada, much of Europe) produce images that include all of their radars. This is not as trivial a task.

In fact, such a network can consist of different types of radar with different characteristics like beam width, wavelength and calibration. These differences have to be taken into account when matching data across the network, particularly to decide what data to use when two radars cover the same point. If one uses the stronger echo but it comes from the most distant radar, one uses returns that are from higher altitude coming from rain or snow that might evaporate before reaching the ground(virga). If one uses data from the closest radar, it might be attenuated passing through a thunderstorm. Composite images of precipitations using a network of radars are made with all those limitations in mind.

Here are some national radar networks :

[edit] Automatic algorithms

The square in this Doppler image has been automatically placed by the radar program to spot the position of a mesocyclone. Notice the inbound/outbound doublet (blue/yellow) with the zero velocity line (gray) parallel to the radial to the radar (up right). It is noteworthy to mention that the change in wind direction here occurs over less than 10 km (6 miles) (Source: Environment Canada).

To help meteorologist to spot dangerous weather, mathematical algorithms have been introduced in the weather radar treatment programs. These are particularly important in the analyzing the Doppler velocity data as they are more complex. The polarization data will even need more algorithms.

Main algorithms for reflectivity:

Vertically Integrated Liquid (VIL) is an estimate of the total mass of precipitation in the clouds.
Potential wind gust, which can estimate the winds under a cloud (a downdraft) using the VIL and the height of the echotops (radar estimated top of the cloud) for a given storm cell.
Hail algorithm that estimate the presence and potential size.

Main algorithms for Doppler velocities:

Mesocyclone detection: it is triggered by a velocity change over a small circular area. The algorithm is searching for a "doublet" of inbound/outbound velocities with the zero line of velocities, between the two, along a radial line from the radar. Usually the mesocyclone detection must be found on two or more stacked progressive tilts of the beam to be significative of rotation into a thunderstom cloud.
TVS or Tornado Vortex Signature algorithm is essentially a mesocyclone with a large velocity threshold found through many scanning angles. This algorithm is used in NEXRAD to indicate the possibility of a tornado formation.
Wind shear in low levels. This algorithm detects variation of wind velocities from point to point in the data and looking for a doublet of inbound/outbound velocities with the zero line perpendicular to the radar beam. The wind shear is associated with downdraft, (downburst and microburst), gust fronts and turbulence under thunderstorms.
VAD Wind Profile (VWP) is a display that estimates the direction and speed of the horizontal wind at various upper levels of the atmoshpere, using the technique explained in the Doppler section.

[edit] Animations

PPI reflectivity loop(in dBZ) showing the evolution of a hurricane (NOAA)

The animation of radar products can show the evolution of reflectivity and velocity patterns. The user can extract information on the dynamics of the meteorological phenomena, including the ability to extrapolate the motion and observe development or dissipation. This can also reveal non-meteorological artifacts (false echoes) that will be discussed later.

[edit] Limitations and artifacts

Radar data interpretation depends on many hypotheses about the atmosphere and the weather targets. They are:

International Standard Atmosphere.
Target small enough they obey the Rayleigh scattering so the return is proportional to the precipitation rate.
The volume scanned by the beam is full of meteorological targets (rain, snow, etc..), all of the same variety and in a uniform concentration.
No attenuation
No amplification
Return from side lobes of the beam are negligible.
The beam is close to a Gaussian function curve with power decreasing to half at half the width.
The outgoing waves and returning one are both polarized similarly.
There is no return from multiple reflections.

One has to keep in mind that those hypotheses are not necessarily met in many circumstances and be able to recognize when the truth from the false echoes.

[edit] Anomalous propagation (non-standard atmosphere)

The first assumption is that the radar beam is moving through air that cools down at a certain rate with height. The position of the echoes depend heavily on this hypothesis. However, the real atmosphere can vary greatly from the norm.

[edit] Super refraction

It is very common to have temperature inversions forming near the ground, for instance air cooling at night while remaining warm aloft. This is not what is expected as the index of refraction of air increase and the radar beam bend toward the Earth instead of going up. Eventually, it will hit the ground and be reflected back toward the radar. The processing program will then wrongly place the return echoes at the height and distance it would have been in normal conditions.

This type of false return is relatively easy to spot on a time loop if it is due to night cooling or marine inversion as one sees very strong echoes developing over an area, spreading in size laterally but not moving and varying greatly in intensity. However, inversion of temperature exist ahead of warm fronts and the abnormal propagation echoes are then mixed with real rain.

The extreme of this problem is when the inversion is very strong and shallow and the radar beam reflects many time on the ground as it has to follow a waveguide path. This will create multiple bands of strong echoes on the radar images.

[edit] Under refraction

On the other hand, if the air is unstable and cool faster that the standard atmosphere with height, the beam ends up higher than expected. This places the precipitation at a much higher altitude they really are. This situation is very difficult to spot.

[edit] Non-Rayleigh targets

If we want to reliably estimate the precipitation rate, the targets have to be 10 times smaller than the radar wave according to Rayleigh scattering. This is due to the fact that the water molecule has to be excited by the radar wave in order to give a return. This is relatively true for rain or snow as 5 or 10 cm radars are used.

However, for very large hydrometeors, since the wavelength is of the order of stone, the return level off according to the Mie scattering. A return of more than 55 dBZ is likely to come from hail but won’t vary proportionally to the size. On the other hand, very small targets like cloud droplets are too small to be excited and don’t give a recordable return on usual weather radars.

[edit] Partially filled scanned volume

Profiler high resolution view of a thunderstorm (top) and by a weather radar (bottom).

As demonstrated at the start of the article, radar beams have a physical dimension and data are sampled every degree, not continuously, along each angle of elevation. This results in an averaging of the values of the returns for reflectivity, velocities and polarization data on the resolution volume scanned.

In the figure to the left, at the top is a view of a thunderstorm taken by a wind profiler when is passed overhead. This is like a vertical cross section through the cloud with 150 m vertical and 30 m horizontal resolution. We can see that the reflectivity has large variations in a short distance. Now compare this with a simulated view of what a regular weather radar would see at 60 km (40 miles) at the bottom. Everything has been smoothed out.

This shows how the output of weather radar is only an approximation of the reality. Naturally, resolution can be improve by newer equipment but some things cannot. As mentioned previously, the volume scanned increase with distance so the possibly that the beam is only partially filled increase too. This leads to underestimating of the precipitation rate at larger distance and fool the user into thinking that rain is lighter as it moves away.

[edit] Beam geometry

Idealized energy distribution of a radar bam(Central lobe at 0 and secondary lobes on each side)

$Diffraction by a circular slit simulating the energy viewed by weather targets$

Diffraction by a circular slit simulating the energy viewed by weather targets

The radar beam is not like a laser but has a distribution of energy similar to the diffraction pattern of a light passing through a slit. This due to the fact that the wave is transmitted to the parabolic antenna through a slit in the wave-guide at the focal point. Most of the energy is at the center of the beam and decrease along a curve close to a Gaussian function on each side as mentioned before. However, there are secondary peaks of emission that will sample the targets at off angles from the center. All is done to minimized the power sent by those lobes but they are never zero.

When a secondary lobe hits a very reflective target, like a mountain or a strong thunderstorm, some of the energy is sent back to the radar. This energy is relatively weak but arrives at the same time the central peak is illuminating a different azimuth. The echo is thus misplaced by the processing program. This has the effect of actually broadening the real weather echo making a smearing of weaker values on each side of it. This causes the user to overestimate the extent of the real echoes.

The strong echoes are returns of the central peak of the radar from a series of small hills (yellow and reds pixels). The weaker echoes on each sides of them are from secondary lobes (blue and green)

[edit] Non weather targets

In the sky there is more than rain and snow. Other objects can be misinterpreted as rain by a weather radar. The main one are:

Birds, especially in period of migration.
Insects at low altitude.
Thin metal strips dropped by military aircraft to fool enemies.
Solids obstacles as mountains, buildings, aircraft.
Ground and sea clutter.

Each of them has their own characteristics that make possible to distinguish them to the trained eye but they may fool a layman. It is possible to eliminate some of them with post-treatment of data using reflectivity, Doppler and polarization data.

[edit] Attenuation

Example of strong attenuation when a line of thunderstorms move over (from left to right images) a 5 cm wavelength weather radar (red arrow). Source: Environment Canada

Micro-waves used in weather radars can be absorbed by rain, depending on the wavelength used. For the 10 centimeter radars, this attenuation is negligible. That is the reason why countries with high water content storms are using 10 centimeter wavelength like in the United Sates with NEXRAD. The cost of a larger antenna, klystron and other related equipments is offset by this benefit.

For a 5 centimeter radar, absorption becomes important in very heavy rain and this attenuation leads to underestimation of echoes in and beyond a strong thunderstorms line. Canada and other northern countries use this less costly kind of radars as their precipitations are usually less intense. However, users have to remember this effect when interpreting data. The images above show how a strong line of echoes seems to vanish as it moves over the radar. To compensate for this behaviour, radar sites are often chosen to somewhat overlap in coverage in order to give different point of view to the same storms.

Shorter wavelength are even more attenuated and are only useful on short range. Many television stations in United States have 3 centimeters radars to cover their listening audience. Knowing their limitations and using them with the local NEXRAD can add information to a meteorologist.

[edit] Bright band

1.5 km altitude CAPPI at the top with strong contamination from the brightband (yellows). The vertical cut at the bottom show that this strong return is only above ground (Source: Environment Canada).

As we have seen previously, the reflectivity depends on the diameter of the target and its capacity to reflect. Snow flakes are large but weakly reflective while rain drops are small but highly reflective.

When snow falls through a layer above freezing temperature, it melts and eventually becomes rain. Using the reflectivity equation, one can demonstrate that the returns from the snow before melting and the rain after, are not too different as the change in dielectric constant compensate for the change in size. However, during the melting process, the radar wave “sees” something akin to very large droplets as snow flakes become coated with water.

This gives enhanced returns that can be mistaken for stronger precipitations. On a PPI, this will show up as an intense ring of precipitations at the altitude where the beam crosses the melting level while on a series of CAPPIs, only the ones near that level will have stronger echoes. A good way to confirm a bright band is to make a vertical cross section through the data like in the picture above.

[edit] Multiples reflections

It is assumed that the beam hits the weather targets and returns directly to the radar. If fact, there is energy reemitted in all directions. Most of it is weak and multiple reflections diminish it even further so what can eventually return to the radar from such an event is negligible. In some case though, this could not be.

For instance, when the beam hits hail, the energy spread toward the wet ground will be reflected back to the hail and then to the radar. The resulting echo is weak but noticeable. Due to the extra path it has to go, it arrives later at the antenna and is placed further than its source. This gives a kind of triangle of false weaker reflectivities radialy behind the hail.

[edit] Solutions for now and the future

Radar image of reflectivity with many non-weather echoes (Source: Environment Canada).

The same image but cleaned using the Doppler velocities. (Source: Environment Canada).

These two images show what can be achieved already to clean up radar data. The output on the left is made with the raw returns and it is difficult to spot the real weather. Since usually rain and snow clouds are moving, one can use the Doppler velocities to eliminate a good part of the clutter. The image on the right has been filtered using this property in a somewhat complex technique.

However, not all non-meteorological targets remain still; one can think of birds for instance. Others, like the bright band, depend on the structure of the precipitations. Polarization offers a direct typing of the echoes which could be used to filter more false data or produce separate images for specialized purposes. This recent development in this field is bound to improve the quality of radar outputs.

Another question is the resolution. As mentioned previously, radar data are an average of the scanned volume by the beam. Resolution can be improved by larger antenna or denser networks. A program by the Center for Collaborative Adaptive Sensing of the Atmosphere (CASA):^[8] aims to supplement regular NEXRAD using many low cost X band (3 cm) weather radar mounted on cellular telephone towers. These radars will subdivide the large area of the NEXRAD into smaller domains to look at altitudes below its lowest angle. These will give details not currently available.

Timeliness is also a point needing improvement. With 5 to 10 minutes time between complete scans of weather radar, a lot of things can be missed in the development of a thunderstorm. A Phased-array radar is been tested at the National Severe Storms Lab in Norman, Oklahoma, to speed up the gathering of data.

Meteorological instrumentation and equipment (Earth based)

[edit] Specialized Applications

[edit] Avionics Weather Radar

Aircraft application of radar systems include weather radar, collision avoidance, target tracking, ground proximity, and other systems. For commercial Weather Radar Systems, ARINC 708 is the primary weather radar system using an airborne pulse Doppler radar.

[edit] References

^ David Atlas, "Radar in Meteorology", published by American Meteorological Society
^ Stormy Weather Group. McGill University ((2000)). Retrieved on 2006-05-21.
^ The King City Operational Doppler Radar: Development, All-Season Applications and Forecasting (PDF). Canadian Meteorological and Oceanographic Society (1990). Retrieved on 2006-05-24.
^ Information about Canadian radar network. The National Radar Program. Environment Canada (2002). Retrieved on 2006-06-14.
^ The PANTHERE project and the evolution of the French operational radar network and products: Rain estimation, Doppler winds, and dual polarization. Météo-France. 32nd Radar Conference of the AMS, Albuquerque, NM (2005). Retrieved on 2006-06-24.
^ M K Yau and R. R. Rogers, "Short Course in Cloud Physics, Third Edition", published by Butterworth-Heinemann
^ Carey, Larry (2003). Lecture on Polarimetric Radar. Texas A&M University. Retrieved on 2006-05-21.
^ List of lectures on CASA. American Meteorological Society ((2005)). Retrieved on 2006-05-21.

[edit] Bibliography

David Atlas, Radar in Meteorology: Battan Memorial and 40th Anniversary Radar Meteorology Conference, published by American Meteorological Society, Boston, 1990, 806 pages, ISBN 0-933876-86-6, AMS Code RADMET.

Yves Blanchard, Le radar, 1904-2004: histoire d'un siècle d'innovations techniques et opérationnelles , published by Ellipses, Paris, France, 2004 ISBN 2-7298-1802-2

R. J. Doviak et D. S. Zrnic, Doppler Radar and Weather Observations, Academic Press. Seconde Edition, San Diego Cal., 1993 p. 562.

M K Yau and R.R. Rogers, Short Course in Cloud Physics, Third Edition, published by Butterworth-Heinemann, January 1, 1989, 304 pages. EAN 9780750632157 ISBN 0-7506-3215-1

Roger M. Wakimoto and Ramesh Srivastava, Radar and Atmospheric Science: A Collection of Essays in Honor of David Atlas, publié par l'American Meteorological Society, Boston, August 2003. Series: Meteorological Monograph, Volume 30, number 52, 270 pages, ISBN 1-878220-57-8; AMS Code MM52.

[edit] See also

[edit] Related articles

Wikimedia Commons has media related to:

Weather radar

[edit] External links

European OPERA project -links to all weather radar maps in the World
History of Operational Use of Weather Radar by U. S. Weather Services. Part I: The Pre-NEXRAD Era.
Weather Radar development (for the USA), highlight of National Severe Storms Laboratory’s first 40 years
Weather radars by Environment Canada
Commons errors in interpreting radar by Environment Canada
Short weather radar course for teachers by Environment Canada
Canadian weather radar FAQ
McGill radar homepage
NEXRAD Doppler radar network information:

Retrieved from "http://en.wikipedia.org../../../w/e/a/Weather_radar.html"

Categories: Radar | Weather radars | Meteorological instrumentation and equipment