A weather radar is a type of radar used in meteorology to locate precipitation, calculate its motion and estimate its type (rain, snow, hail, etc..). Furthermore, the tridimentionnal data obtained can be analyzed to extract the structure of storms and their potential for damages. Finally, precipitations and clear air echos obtained by weather radar permit to estimate wind direction and speed in the lowest part of the atmosphere.
Weather radar is often used together with lightning detectors to detect major storm activity.
Weather radars are pulsed radars. The micro-wave generators are usually magnetron or klystron with wavelength of 1 to 10 cm. The wave is transmitted through wave-guide to a parabolic antenna toward targets.
Contrary to aircraft detection, weather targets are numerous within the volume of the radar beam: (h is pulse width, r distance to the radar and the beam width).
With a typical radar pulse and beam width, the scanned volume vary greatly from close to the radar to the end of the sounding range of 150 to 200 km. For instance, the returns from a far away distance will be an average of the echoes in a volume of the order of a cube kilometer of air.
Because the targets are not unique in each volume, the radar equation has to be developed beyond the basic one:
In this case, we have to add the cross sections of all the targets:
where is the light speed, is temporal duration of a pulse and is the beam width in radians.
In combining the two eqations :
Which leads to:
Notice that the return now varies inversely to instead of . In order to compare the data coming from different distances from the radar, one has to normalize them with this ratio.
Between each pulse, the antenna serve as receptor for the return from targets. The distance of the latter are:
So the maximum non ambiguous range depends on t between pulses. Any pulse returning AFTER a new pulse has been emitted will be misplaced as it is assumed that it comes from the second pulse.
Assuming the Earth is round, the variation of the index of refraction through air and the distance to the target, one can calculate the height above ground of it.
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.
Reflectivity (Z) varies by the 6th power of the rain droplets diameter D and the square of the dielectric constant (K) of the target. As the drop distribution (N*) is a truncated Gamma functionM K Yau and R. R. Rogers, "Short Course in Cloud Physics, Third Edition", published by Butterworth-Heinemann, its equation becomes:
Precipitation rate (R), on the other hand, is equal to the number of particles, their volume and their fall speed (v*) as:
So Z and R are related by:
Where a and b depend on the type of precipitations (snow, rain, convective or stratiform) which have different , K, N0 and v.
The frequency difference in the return from moving rain droplets or snow flakes are too small to be noted by actual electronic instruments. 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%. However, as they move slightly between each pulse, the returned wave has a noticeable phase difference from pulse to pulse.
Doppler 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 :
So
v = target speed =
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.
If we now look at the maximum velocity that can be deduced from pulse pairs, a sinus can vary from - and +, so one cannot resolve a greater velocity than:
This is called the Nyquist velocity. This is directly dependant on the time between successive pulses: the smaller it is, the larger is the non ambiguous range of speed. However, we know that the maximum range from reflectivity is inversely dependant on :
So we have a dilemma : increasing the range for reflectivity at the expense of velocity definition or increasing the latter at the expense of range. With the wavelengths used, the compromise has been the use a Pulse Repetition Rate that gives 100 to 150 km range.
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 will be a Cosinus. One can then deduce 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. Luckily, this component is negligible in small elevation angles but must be taken into account for higher scanning angles.
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 :
With this new knowledge, the reflectivity and the Doppler data, researcher 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:
McGill University operational output
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.
See : Plan Position Indicator
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
See: 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:
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
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.
Aviation is a heavy user of radar data. One map particularly important in this field is the Echo tops 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. So it is easy to produce the maximum height at which precipitation are found in this volume. However one has to remember that those are not the tops of clouds since it extended further up than the precipitations.
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.
For the past few decades, radar networks have been extented to the point that composite views covering large areas can be produces. For instance, all major countries like United States, Canada, European countries, etc... produce images including all their radars. This is not as trivial a task as it may seem.
The fact is that such a network can consist of different types of radar that have different characteristics like beam width, wavelength and calibrations. This has to be taken into account when matching data from one end to the other of the network. What data to use when two radars cover the same point with their PPI? If one use 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 done with all those limitations in mind.
Here are some national radar networks :
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 has they are more complex. The polarization data will even need more algorithms.
Main algorithms for reflectivity:
Main algorithms for Doppler velocities:
Radar data interpretation depends on many hypotheses about the atmosphere and the weather targets. They are:
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.
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.
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.
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.
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.
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.
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:
Each of them has their own characteristics that make possible to distinguish them to the trained eye but they may fool a layment. It is possible to eliminate some of them with post-treatment of data using reflectivity, Doppler and polarization data.
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.
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.
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.
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 propriety in a somewhat complex technique.
However, not all non meteorological targets are remaining still, one can think of birds for instance. Others, like the bright band, depend on the structure of the precipitations. Polarization offer a direct typing of the echoes which could be used to filter more false data or produce separate images for specialized purposes. This recent developments 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) :aim 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. They will then be giving details not available at this moment.
Timeliness is also a point to improve. With 5 to 10 minutes time between complete scans of weather radar, 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.
Radar | Weather_radars | Meteorological instrumentation and equipment
Niederschlagsradar | Radar météorologique | Radar Doppler | Weerradar | ドップラー・レーダー | Dopplertutka
This article is licensed under the GNU Free Documentation License.
It uses material from the
"Weather radar".
Home Page • arts • business • computers • games • health • hospitals • home • kids & teens • news • physicians • recreation• reference • regional • science • shopping • society • sports • world