Weather radar, also called weather surveillance radar (WSR) and Doppler weather radar, is a type of radar used to locate precipitation, calculate its motion, and estimate its type (rain, snow, hail etc.). Modern weather radars are mostly pulse-Doppler radars, capable of detecting the motion of rain droplets in addition to the intensity of the precipitation. Both types of data can be analyzed to determine the structure of storms and their potential to cause severe weather.
During World War II, radar operators discovered that weather was causing echoes on their screen, masking potential enemy targets. Techniques were developed to filter them, but scientists began to study the phenomenon. Soon after the war, surplus radars were used to detect precipitation. Since then, weather radar has evolved on its own and is now used by national weather services, research departments in universities, and in television newscasts. Raw images are routinely used and specialized software can take radar data to make short term forecasts of future positions and intensities of rain, snow, hail, and other weather phenomena. Radar output is even incorporated into numerical weather prediction models to improve analyses and forecasts.
During World War II, military radar operators noticed noise in returned echoes due to rain, snow, and sleet. 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, at first working for the Air Force and later for MIT, developed the first operational weather radars. In Canada, J.S. Marshall and R.H. Douglas formed the "Stormy Weather Group" 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 rainwater is falling. 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 centimeters. By 1950 the UK company EKCO was demonstrating its airborne 'cloud and collision warning search radar equipment'.
Between 1950 and 1980, reflectivity radars, which measure position and intensity of precipitation, were incorporated 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 cross-sections could be performed. Studies of the organization of thunderstorms were then possible for the Alberta Hail Project in Canada and National Severe Storms Laboratory (NSSL) in the US in particular.
The NSSL, created in 1964, began experimentation on dual polarization signals and on Doppler effect uses. In May 1973, a tornado devastated Union City, Oklahoma, just west of Oklahoma City. For the first time, a Dopplerized 10 cm wavelength radar from NSSL documented the entire life cycle of the tornado. The researchers discovered a mesoscale rotation in the cloud aloft before the tornado touched the ground – the tornadic vortex signature. NSSL's research helped convince the National Weather Service that Doppler radar was a crucial forecasting tool. The Super Outbreak of tornadoes on 3–4 April 1974 and their devastating destruction might have helped to get funding for further developments.
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 could track the relative velocity of the particles in the air. In the United States, the construction of a network consisting of 10 cm radars, called NEXRAD or WSR-88D (Weather Surveillance Radar 1988 Doppler), was started in 1988 following NSSL's research. In Canada, Environment Canada constructed the King City station, with a 5 cm research Doppler radar, by 1985; McGill University dopplerized its radar (J. S. Marshall Radar Observatory) in 1993. This led to a complete Canadian Doppler network between 1998 and 2004. France and other European countries had switched to Doppler networks by the early 2000s. Meanwhile, rapid advances in computer technology led to algorithms to detect signs of severe weather, and many applications for media outlets and researchers.
After 2000, research on dual polarization technology 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 was done by the end of the decade or the beginning of the next in some countries such as the United States, France, and Canada. In April 2013, all National Weather Service NEXRADs were completely dual-polarized.
Since 2003, the U.S. National Oceanic and Atmospheric Administration has been experimenting with phased-array radar as a replacement for conventional parabolic antenna to provide more time resolution in atmospheric sounding. This could be significant with severe thunderstorms, as their evolution can be better evaluated with more timely data.
Also in 2003, the National Science Foundation established the Engineering Research Center for Collaborative Adaptive Sensing of the Atmosphere (CASA), a multidisciplinary, multi-university collaboration of engineers, computer scientists, meteorologists, and sociologists to conduct fundamental research, develop enabling technology, and deploy prototype engineering systems designed to augment existing radar systems by sampling the generally undersampled lower troposphere with inexpensive, fast scanning, dual polarization, mechanically scanned and phased array radars.
How a weather radar works
Sending radar pulses
Weather radars send directional pulses of microwave radiation, on the order of a microsecond long, using a cavity magnetron or klystron tube connected by a waveguide to a parabolic antenna. The wavelengths of 1 – 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 but is more expensive than a 5 cm C-band system. 3 cm X-band radar is used only for short-range units, and 1 cm Ka-band weather radar is used only for research on small-particle phenomena such as drizzle and fog. W-band weather radar systems have seen limited university use, but due to quicker attenuation, most data are not operational.
Radar pulses spread out as they move away from the radar station. Thus the volume of air that a radar pulse is traversing is larger for areas farther away from the station, and smaller for nearby areas, decreasing resolution at far 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. This is called the pulse volume
Listening for return signals
Between each pulse, the radar station serves as a receiver as it 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 elapses 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 in air:
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.
Assuming the Earth is round, the radar beam in vacuum would rise according to the reverse curvature of the Earth. However, the atmosphere has a refractive index that diminishes with height, due to its diminishing density. This bends the radar beam slightly toward the ground and with a standard atmosphere this is equivalent to considering that the curvature of the beam is 4/3 the actual curvature of the Earth. Depending on the elevation angle of the antenna and other considerations, the following formula may be used to calculate the target's height above ground:
- r = distance radar–target,
- ke = 4/3,
- ae = Earth radius,
- θe = elevation angle above the radar horizon,
- ha = height of the feedhorn above ground.
A weather radar network uses a series of typical angles that will be set according to the needs. After each scanning rotation, the antenna elevation is changed for the next sounding. This scenario will be repeated on many angles 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 to have data within 15 km above ground and 250 km distance of the radar. For instance in Canada, the 5 cm weather radars use angles ranging from 0.3 to 25 degrees. The image to the right shows the volume scanned when multiple angles are used.
Due to the Earth's curvature and change of index of refraction with height, the radar cannot "see" below the height above ground of the minimal angle (shown in green) or closer to the radar than the maximal one (shown as a red cone in the center).
Calibrating intensity of return
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 equations:
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.
Reflectivity (in decibel or dBZ)
Return echoes from targets ("reflectivity") are analyzed for their intensities to establish the precipitation rate in the scanned volume. The wavelengths used (1–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 (Ze) varies by the sixth 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, of the form:
Precipitation rate (R), on the other hand, is equal to the number of particles, their volume and their fall speed (v[D]) as:
So Ze and R have similar functions that can be resolved giving a relation between the two of the form:
- Z = aRb
- As the antenna scans the atmosphere, on every angle of azimuth it obtains a certain strength of return from each type of target encountered. Reflectivity is then averaged for that target to have a better data set.
- Since variation in diameter and dielectric constant of the targets can lead 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).
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:
- magenta: 65 dBZ (extremely heavy precipitation, possible hail)
- 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 below.
When describing weather radar returns, pilots, dispatchers, and air traffic controllers will typically refer to three return levels:
- 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 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.
Some displays provided by commercial weather sites, like The Weather Channel, show precipitation types during the winter month : rain, snow, mixed precipitations (sleet and freezing rain). This is not an analysis of the radar data itself but a post-treatment done with other data sources, the primary being surface reports (METAR).
Over the area covered by radar echoes, a program assigns a precipitation type according to the surface temperature and dew point reported at the underlying weather stations. Precipitation types reported by human operated stations and certain automatic ones (AWOS) will have higher weight. Then the program does interpolations to produce an image with defined zones. These will include interpolation errors due to the calculation. Mesoscale variations of the precipitation zones will also be lost. More sophisticated programs use the numerical weather prediction output from models, such as NAM and WRF, for the precipitation types and apply it as a first guess to the radar echoes, then use the surface data for final output.
Until dual-polarization (section below) data are widely available, any precipitation types on radar images are only indirect information and must be taken with care.
Precipitation is found in and below clouds. Light precipitation such as drops and flakes is subject to the air currents, and scanning radar can pick up the horizontal component of this motion, thus giving the possibility to estimate the wind speed and direction where precipitation is present.
A target's motion relative to the radar station causes a change in the reflected frequency of the radar pulse, due to the Doppler effect. With velocities of less than 70-metre/second for weather echos and radar wavelength of 10 cm, this amounts to a change only 0.1 ppm. 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 use this phase difference (pulse pair difference) to calculate the precipitation's 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.
The phase between pulse pairs can vary from - and +, so the unambiguous Doppler velocity range is
- Vmax =
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 :
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–150 km for reflectivity. This means for a wavelength of 5 cm (as shown in the diagram), an unambiguous velocity range of 12.5 to 18.75 metre/second is produced (for 150 km and 100 km, respectively). For a 10 cm radar such as the NEXRAD, the unambiguous velocity range would be doubled.
Some techniques using two alternating pulse repetition frequencies (PRF) allow a greater 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 metre/second 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 will differ thereafter. It is then possible to find a mathematical relation between the two returns and calculate the real velocity beyond the limitation of the two PRFs.
In a uniform rainstorm moving eastward, a radar beam pointing west will "see" the raindrops moving toward itself, while a beam pointing east will "see" the drops moving away. When the beam scans to the north or to the south, no relative motion is noted.
In the synoptic scale interpretation, the user can extract the wind at different levels over the radar coverage region. As the beam is scanning 360 degrees around the radar, data will come 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 (maximum in the precipitation motion and zero in the perpendicular direction). 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.
In the velocity data, there could be smaller zones in the radar coverage where the wind varies from the one mentioned above. For example, a thunderstorm is a mesoscale phenomenon which often includes rotations and turbulence. These may only cover few square kilometers but are visible by variations in the radial speed. Users can recognize velocity patterns in the wind associated with rotations, such as mesocyclone, convergence (outflow boundary) and divergence (downburst).
Droplets of falling liquid water tend to 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 in order to receive the maximal signal reflection.
If two pulses are sent simultaneously with orthogonal polarization (vertical and horizontal, ZV and ZH respectively), two independent sets of data will be received. These signals can be compared in several useful ways:
- Differential Reflectivity (Zdr) – The differential reflectivity is the ratio of the reflected vertical and horizontal power returns as ZV/ZH. Among other things, it is a good indicator of drop shape and drop shape is a good estimate of average drop size.
- Correlation Coefficient (ρhv) – A statistical correlation between the reflected horizontal and vertical power returns. High values, near one, indicate homogeneous precipitation types, while lower values indicate regions of mixed precipitation types, such as rain and snow, or hail, or in extreme cases debris aloft, usually coinciding with a Tornado vortex signature.
- Linear Depolarization Ratio (LDR) – This is a ratio of a vertical power return from a horizontal pulse or a horizontal power return from a vertical pulse. It can also indicate regions where there is a mixture of precipitation types.
- Differential Phase (θdp) – The differential phase is a comparison of the returned phase difference between the horizontal and vertical pulses. This change in phase is caused by the difference in the number of wave cycles (or wavelengths) along the propagation path for horizontal and vertically polarized waves. It should not be confused with the Doppler frequency shift, which is caused by the motion of the cloud and precipitation particles. Unlike the differential reflectivity, correlation coefficient and linear depolarization ratio, which are all dependent on reflected power, the differential phase is a "propagation effect." It is a very good estimator of rain rate and is not affected by attenuation. The range derivative of differential phase (specific differential phase, Kdp) can be used to localize areas of strong precipitation/attenuation.
With this new knowledge added to the reflectivity, velocity, and spectrum width produced by Doppler weather radars, researchers have been working on developing algorithms to differentiate precipitation types, non-meteorological targets, and to produce better rainfall accumulation estimates. In the U.S., NCAR and NSSL have been world leaders in this field.
NOAA established a test deployment for dual-polametric radar at NSSL and equipped all its 10 cm NEXRAD radars with dual-polarization, which was completed in April 2013. In 2004, ARMOR Doppler Weather Radar in Huntsville, Alabama was equipped with a SIGMET Antenna Mounted Receiver, giving Dual-Polarmetric capabilities to the operator. McGill University J. S. Marshall Radar Observatory in Montreal, Canada has converted its instrument (1999) and the data are used operationally by Environment Canada in Montreal. Another Environment Canada radar, in King City (North of Toronto), was dual-polarized in 2005; it uses a 5 cm wavelength, which experiences greater attenuation. Environment Canada is working on converting all of its radars to dual-polarization. Météo-France is planning on incorporating dual-polarizing Doppler radar in its network coverage.
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.
Plan position indicator
Since data are obtained one angle at a time, the first way of displaying them has 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 reaches the ground but what is seen from 160 km 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 distances 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 such shortcomings.
Usage: Reflectivity, Doppler and polarimetric data can use PPI.
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.
Constant-altitude plan position indicator
To avoid some of the problems on PPIs, the constant-altitude plan position indicator (CAPPI) has been developed by Canadian researchers. 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 to have a cut that is as close as possible at all distance to the height needed. 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 1.5 km and 4 km height CAPPIs. Notice that the section after 120 km is using the same data.
Since the CAPPI uses the closest angle to the desired height at each point from the radar, the data can originate from slightly different altitudes, as seen on the image, in different points of the radar coverage. It is therefore crucial to have a large enough number of sounding angles to minimize this height change. Furthermore, the type of data must be changing relatively gradually with height to produce an image that is not noisy.
Reflectivity data being relatively smooth with height, CAPPIs are mostly used for displaying them. Velocity data, on the other hand, can change rapidly in direction with height and CAPPIs of them are not common. It seems that only McGill University is producing regularly Doppler CAPPIs with the 24 angles available on their radar. However, some researchers have published papers using velocity CAPPIs to study tropical cyclones and development of NEXRAD products. Finally, polarimetric data are recent and often noisy. There doesn't seem to have regular use of CAPPI for them although the SIGMET company offer a software capable to produce those types of images.
- 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 assures that no strong echo is missed in the layer and a treatment using Doppler velocities eliminates the ground echoes. Comparing base and composite products, one can locate virga and updrafts zones.
Real time example:
Another important use of radar data is the ability to assess the amount of precipitation that has fallen over large basins, to be used in hydrological calculations; such data is useful in flood control, sewer management and dam construction. The computed data from radar weather may be used in conjunction with data from ground stations.
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 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 always extend above the precipitation.
Vertical cross sections
To know the vertical structure of clouds, in particular thunderstorms or the level of the melting layer, a vertical cross-section product of the radar data is available. This is done by displaying only the data along a line, from coordinates A to B, taken from the different angles scanned.
Range Height Indicator
When a weather radar is scanning in only one direction vertically, it obtains high resolution data along a vertical cut of the atmosphere. The output of this sounding is called a Range Height Indicator (RHI) which is excellent for viewing the detailed vertical structure of a storm. This is different from the vertical cross section mentioned above by the fact that the radar is making a vertical cut along specific directions and does not scan over the entire 360 degrees around the site. This kind of sounding and product is only available on research radars.
Over the past few decades, radar networks have been extended to allow the production of composite views covering large areas. For instance, many countries, including the United States, Canada and much of Europe, produce images that include all of their radars. This is not a trivial task.
In fact, such a network can consist of different types of radar with different characteristics such as 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 more 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 closer 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.
To help meteorologists spot dangerous weather, mathematical algorithms have been introduced in the weather radar treatment programs. These are particularly important in 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.
- VIL Density is VIL divided by the height of the cloud top. It is a clue to the possibility of large hail in thunderstorms.
- 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 algorithms that estimate the presence of hail and its probable 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 thunderstorm 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 atmosphere, using the technique explained in the Doppler section.
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.
Radar Integrated Display with Geospatial Elements
A new popular presentation of weather radar data in United States is via Radar Integrated Display with Geospatial Elements (RIDGE) in which the radar data is projected on a map with geospatial elements such as topography maps, highways, state/county boundaries and weather warnings. The projection often is flexible giving the user a choice of various geographic elements. It is frequently used in conjunction with animations of radar data over a time period.
Limitations and artifacts
Radar data interpretation depends on many hypotheses about the atmosphere and the weather targets, including:
- International Standard Atmosphere.
- Targets small enough to obey the Rayleigh scattering, resulting in the return being 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 and returning waves are similarly polarized.
- There is no return from multiple reflections.
These assumptions are not always met; one must be able to differentiate between reliable and dubious echoes.
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.
Temperature inversions often form near the ground, for instance by air cooling at night while remaining warm aloft. As the index of refraction of air decreases faster than normal the radar beam bends toward the ground instead of continuing upward. 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 exists 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, the radar beam reflects many times toward the ground as it has to follow a waveguide path. This will create multiple bands of strong echoes on the radar images.
This situation can be found with inversions of temperature aloft or rapid decrease of moisture with height. In the former case, it could be difficult to notice.
On the other hand, if the air is unstable and cools faster than the standard atmosphere with height, the beam ends up higher than expected. This indicates that precipitation is occurring higher than the actual height. Such an error is difficult to detect without additional temperature lapse-rate data for the area.
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 because the water molecule has to be excited by the radar wave to give a return. This is relatively true for rain or snow as 5 or 10 cm wavelength radars are usually employed.
However, for very large hydrometeors, since the wavelength is on the order of stone, the return levels off according to Mie theory. 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 such as cloud droplets are too small to be excited and do not give a recordable return on common weather radars.
Resolution and partially filled scanned volume
As demonstrated at the start of the article, radar beams have a physical dimension and data are sampled at discrete angles, 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 as it was passing overhead. This is like a vertical cross section through the cloud with 150-metre vertical and 30-metre horizontal resolution. The reflectivity has large variations in a short distance. Compare this with a simulated view of what a regular weather radar would see at 60 km, in the bottom of the figure. Everything has been smoothed out. Not only the coarser resolution of the radar blur the image but the sounding incorporates area that are echo free, thus extending the thunderstorm beyond its real boundaries.
This shows how the output of weather radar is only an approximation of reality. The image to the right compares real data from two radars almost colocated. The TDWR has about half the beamwidth of the other and one can see twice more details than with the NEXRAD.
Resolution can be improved by newer equipment but some things cannot. As mentioned previously, the volume scanned increases with distance so the possibility that the beam is only partially filled also increases. This leads to underestimation of the precipitation rate at larger distances and fools the user into thinking that rain is lighter as it moves away.
The radar beam has a distribution of energy similar to the diffraction pattern of a light passing through a slit. This is because 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 decreases along a curve close to a Gaussian function on each side. However, there are secondary peaks of emission that will sample the targets at off-angles from the center. Designers attempt to minimize the power transmitted by such lobes, but they cannot be completely eliminated.
When a secondary lobe hits a reflective target such as a mountain or a strong thunderstorm, some of the energy is reflected to the radar. This energy is relatively weak but arrives at the same time that 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.
There is more than rain and snow in the sky. Other objects can be misinterpreted as rain or snow by weather radars. Insects and arthropods are swept along by the prevailing winds, while birds follow their own course. As such, fine line patterns within weather radar imagery, associated with converging winds, are dominated by insect returns. Bird migration, which tends to occur overnight within the lowest 2000 metres of the Earth's atmosphere, contaminates wind profiles gathered by weather radar, particularly the WSR-88D, by increasing the environmental wind returns by 30–60 km/hr. Other objects within radar imagery include:
- Thin metal strips (chaff) dropped by military aircraft to fool enemies.
- Solid obstacles such as mountains, buildings, and aircraft.
- Ground and sea clutter.
- Reflections from nearby buildings ("urban spikes").
Such extraneous objects have characteristics that allow a trained eye to distinguish them. It is also possible to eliminate some of them with post-treatment of data using reflectivity, Doppler, and polarization data.
The rotating blades of windmills on modern wind farms can return the radar beam to the radar if they are in its path. Since the blades are moving, the echoes will have a velocity and can be mistaken for real precipitation. The closer the wind farm, the stronger the return, and the combined signal from many towers is stronger. In some conditions, the radar can even see toward and away velocities that generate false positives for the tornado vortex signature algorithm on weather radar; such an event occurred in 2009 in Dodge City, Kansas.
As with other structures that stand in the beam, attenuation of radar returns from beyond windmills may also lead to underestimation.
Microwaves used in weather radars can be absorbed by rain, depending on the wavelength used. For 10 cm radars, this attenuation is negligible. That is the reason why countries with high water content storms are using 10 cm wavelength, for example the US NEXRAD. The cost of a larger antenna, klystron and other related equipment is offset by this benefit.
For a 5 cm radar, absorption becomes important in heavy rain and this attenuation leads to underestimation of echoes in and beyond a strong thunderstorm. Canada and other northern countries use this less costly kind of radar as the precipitation in such areas is usually less intense. However, users must consider this characteristic 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 to give different points of view of the same storms.
Shorter wavelengths are even more attenuated and are only useful on short range radar. Many television stations in the United States have 5 cm radars to cover their audience area. Knowing their limitations and using them with the local NEXRAD can supplement the data available to a meteorologist.
A radar beam's reflectivity depends on the diameter of the target and its capacity to reflect. Snowflakes are large but weakly reflective while rain drops are small but highly reflective.
When snow falls through a layer above freezing temperature, it melts into 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 compensates 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 precipitation 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, as illustrated in the picture above.
An opposite problem is that drizzle (precipitation with small water droplet diameter) tends not to show up on radar because radar returns are proportional to the sixth power of droplet diameter.
It is assumed that the beam hits the weather targets and returns directly to the radar. In fact, there is energy reflected 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. However, some situations allow a multiple-reflected radar beam to be received by the radar antenna. 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 length it has to go through, it arrives later at the antenna and is placed further than its source. This gives a kind of triangle of false weaker reflections placed radially behind the hail.
Solutions for now and the future
These two images show what can be presently achieved 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 rain and snow clouds are usually moving, one can use the Doppler velocities to eliminate a good part of the clutter (ground echoes, reflections from buildings seen as urban spikes, anomalous propagation). The image on the right has been filtered using this property.
However, not all non-meteorological targets remain still (birds, insects, dust). Others, like the bright band, depend on the structure of the precipitation. 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 is expected to improve the quality of radar products.
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) aims to supplement the regular NEXRAD (a network in the United States) 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.
Using 3 cm radars, the antenna of each radar is small (about 1 meter diameter) but the resolution is similar at short distance to that of NEXRAD. The attenuation is significant due to the wavelength used but each point in the coverage area is seen by many radars, each viewing from a different direction and compensating for data lost from others.
The number of elevation scanned and the time taken for a complete cycle depend on the weather situation. For instance, with little or no precipitation, the scheme may be limited the lowest angles and using longer impulses in order to detect wind shift near the surface. On the other hand, in violent thundertorm situations, it is better to scan on a large number of angles in order to have a 3 dimensions view of the precipitations as often as possible. To mitigate those different demands, scanning strategies have been developed according to the type of radar, the wavelength used and the most commons weather situations in the area considered.
One example of scanning strategies is given by the US NEXRAD radar network which has evolved with time. For instance, in 2008, it added extra resolution of data, and in 2014, additional intra-cycle scanning of lowest level elevation (MESO-SAILS).
Timeliness is also a point needing improvement. With 5 to 10 minutes time between complete scans of weather radar, much data is lost as a thunderstorm develops. A Phased-array radar is being tested at the National Severe Storms Lab in Norman, Oklahoma, to speed the data gathering. A team in Japan has deployed a phased-array radar for 3D NowCasting at the RIKEN Advanced Institute for Computational Science (AICS).
Avionics weather radar
Aircraft application of radar systems include weather radar, collision avoidance, target tracking, ground proximity, and other systems. For commercial weather radar, ARINC 708 is the primary specification for weather radar systems using an airborne pulse-Doppler radar.
Unlike ground weather radar, which is set at a fixed angle, airborne weather radar is being utilized from the nose or wing of an aircraft. Not only will the aircraft be moving up, down, left, and right, but it will be rolling as well. To compensate for this, the antenna is linked and calibrated to the vertical gyro located on the aircraft. By doing this, the pilot is able to set a pitch or angle to the antenna that will enable the stabilizer to keep the antenna pointed in the right direction under moderate maneuvers. The small servo motors will not be able to keep up with abrupt maneuvers, but it will try. In doing this the pilot is able to adjust the radar so that it will point towards the weather system of interest. If the airplane is at a low altitude, the pilot would want to set the radar above the horizon line so that ground clutter is minimized on the display. If the airplane is at a very high altitude, the pilot will set the radar at a low or negative angle, to point the radar towards the clouds wherever they may be relative to the aircraft. If the airplane changes attitude, the stabilizer will adjust itself accordingly so that the pilot doesn't have to fly with one hand and adjust the radar with the other.
There are two major systems when talking about the receiver/transmitter: the first is high-powered systems, and the second is low-powered systems; both of which operate in the X-band frequency range (8,000 – 12,500 MHz). High-powered systems operate at 10,000 – 60,000 watts. These systems consist of magnetrons that are fairly expensive (approximately $1,700) and allow for considerable noise due to irregularities with the system. Thus, these systems are highly dangerous for arcing and are not safe to be used around ground personnel. However, the alternative would be the low-powered systems. These systems operate 100 – 200 watts, and require a combination of high gain receivers, signal microprocessors, and transistors to operate as effectively as the high-powered systems. The complex microprocessors help to eliminate noise, providing a more accurate and detailed depiction of the sky. Also, since there are fewer irregularities throughout the system, the low-powered radars can be used to detect turbulence via the Doppler Effect. Since low-powered systems operate at considerable less wattage, they are safe from arcing and can be used at virtually all times.
Digital radar systems now have capabilities far beyond that of their predecessors. Digital systems now offer thunderstorm tracking surveillance. This provides users with the ability to acquire detailed information of each storm cloud being tracked. Thunderstorms are first identified by matching precipitation raw data received from the radar pulse to some sort of template preprogrammed into the system. In order for a thunderstorm to be identified, it has to meet strict definitions of intensity and shape that set it apart from any non-convective cloud. Usually, it must show signs of organization in the horizontal and continuity in the vertical: a core or a more intense center to be identified and tracked by digital radar trackers. Once the thunderstorm cell is identified, speed, distance covered, direction, and Estimated Time of Arrival (ETA) are all tracked and recorded to be utilized later.
Doppler radar and bird migration
Using the Doppler weather radar is not limited to determine the location and velocity of precipitation, but it can track bird migrations as well as seen in the section. The radio waves sent out by the radars bounce off rain and birds alike (or even insects like butterflies). The US National Weather Service, for instance, have reported having the flights of birds appear on their radars as clouds and then fade away when the birds land. The U.S. National Weather Service St. Louis has even reported monarch butterflies appearing on their radars.
Different programs in North America use regular weather radars and specialized radar data to determine the paths, height of flight, and timing of migrations. This is useful information in planning for windmill farms placement and operation, to reduce bird fatalities, aviation safety and other wildlife management. In Europe, there has been similar developments and even a comprehensive forecast program for aviation safety, based on radar detection.
Meteorite fall detection
According to the American Meteor Society, meteorite falls occur on a daily basis somewhere on Earth. However, the database of worldwide meteorite falls maintained by the Meteoritical Society typically records only about 10-15 new meteorite falls annually
Meteorite occur when a meteoroid falls into the Earth's atmosphere, generating an optically bright meteor by ionization and frictional heating. If the meteoroid is large enough and infall velocity is low enough, surviving meteorites will reach the ground. When the falling meteorites decelerate below about 2–4 km/s, usually at an altitude between 15 and 25 km, they no longer generate an optically bright meteor and enter "dark flight". Because of this, most meteorite falls occurring into the oceans, during the day, or otherwise go unnoticed.
It is in dark flight that falling meteorites typically fall through the interaction volumes of most types of radars. It as been demonstrated that it is possible to identify falling meteorites in weather radar imagery by different studies. This is especially useful for meteorite recovery, as weather radar are part of widespread networks and scan the atmosphere continuously. Furthermore, the meteorites are cause a perturbance of local winds by turbulence, which is noticeable on Doppler outputs, and are falling nearly vertically so their resting place on the ground is close to their radar signature.
An example is shown in the image on the right, showing the Park Forest, IL meteorite fall which occurred on 26 March 2003. The red/green feature at upper left is the motion of clouds near the radar itself, and a signature of falling meteorites is seen inside the yellow ellipse at image center. The intermixed red and green pixels indicate turbulence, in this case arising from the wakes of falling, high-velocity meteorites.
- David Atlas, "Radar in Meteorology", published by American Meteorological Society
- . McGill University. 2000. Retrieved 2006-05-21.
- . Colorado State University. 2008. Retrieved 2008-01-30.
- Cobb, Susan (29 October 2004). . NOAA Magazine. National Oceanic and Atmospheric Administration. Retrieved 2009-03-07.
- . NSSL. Retrieved 1 March 2014.
- Crozier, C.L.; Joe, P.I.; Scott, J.W.; Herscovitch, H.N.; Nichols, T.R. (1991). (PDF). Atmosphere-Ocean. Canadian Meteorological and Oceanographic Society (CMOS). 29 (3): 479–516. doi:. Archived from on 2 October 2006. Retrieved 10 May 2012.
- . The National Radar Program. Environment Canada. 2002. Archived from on 29 June 2004. Retrieved 2006-06-14.
- Golbon-Haghighi, M.H.; Zhang G.; Li Y.; Doviak R. (June 2016). . Atmosphere. 7 (6): 83. Bibcode:. doi:.
- [url=] The PANTHERE project and the evolution of the French operational radar network and products: Rain estimation, Doppler winds, and dual polarization, Parent du Châtelet, Jacques et al. Météo-France (2005) 32nd Radar Conference of the American Meteorological Society, Albuquerque NM
- National Weather Service (25 April 2013). . NOAA. Retrieved 26 April 2013.
- Doviak, R. J.; Zrnic, D. S. (1993). Doppler Radar and Weather Observations (2nd ed.). San Diego CA: Academic Press. ISBN 0-12-221420-X.
- (in English) . Glossary of Meteorology. American Meteorological Society. 2009. Retrieved 2009-09-27.
- de Podesta, M (2002). . CRC Press. p. 131. ISBN 0-415-25788-3.
- Doviak, R.J.; Zrnic, D. S. (1993). (PDF).
- Airbus (14 March 2007). (PDF). SKYbrary. p. 2. Retrieved 2009-11-19.
- Skolnik, Merrill I. (22 January 2008). "1.2". (pdf) (3rd ed.). McGraw-Hill. ISBN 9780071485470. Retrieved 2016-04-01.
- Skolnik, Merrill I. (22 January 2008). "19.2". (pdf) (3rd ed.). McGraw-Hill. ISBN 9780071485470. Retrieved 2016-04-01.
- Yau, M.K.; Rogers, R.R. Short Course in Cloud Physics (3rd ed.). Butterworth-Heinemann. ISBN 0-08-034864-5.
- National Weather Service. . WSR-88D Radar FAQs. National Oceanic and Atmospheric Administration. Retrieved 2009-12-15.
- Stoen, Hal (27 November 2001). . Aviation Tutorials Index. stoenworks.com. Retrieved 2009-12-15.
- Haby, Jeff. . Nowcasting winter precipitation on the Internet. theweatherprediction.com. Retrieved 2009-12-14.
- . Types of Maps. The Weather Network. Retrieved 2009-12-14.
- Carey, Larry (2003). (PDF). Texas A&M University. Retrieved 2006-05-21.
- Schuur, Terry. . CIMMS. National Severe Storms Laboratory. Retrieved 19 April 2013.
- (PDF). 2012-08-03. Retrieved 2013-05-09.
- National Weather Service. (pdf). NOAA. Retrieved 18 April 2013.
- Schuur, Terry. . CIMMS. National Severe Storms Laboratory. Retrieved 19 April 2013.
- Schurr, Terry; Heinselman, P.; Scharfenberg, K. (October 2003). (PDF). NSSL and CIMMS. Retrieved 2013-04-19.
- Fabry, Frédéric; J. S. Marshall Radar Observatory. . McGill University. Archived from on 10 June 2008. Retrieved 2013-04-18.
- J. S. Marshall Radar Observatory. . McGill University. Retrieved 2013-04-18.
- Ryzhkov; Giangrande; Krause; Park; Schuur; Melnikov. . Doppler Weather Radar Research and Development. CIMMS. Retrieved 2009-02-12.
- Doviak, R. J.; Zrnic, D. S. (1993). Doppler Radar and Weather Observations. San Diego Cal.: Academic Press. p. 562.
- Government of Canada (25 January 2012). . Environnement Canada. Retrieved 29 October 2012.
- Parent du Châtelet, Jacques; Météo-France; et al. (2005). (pdf). 32nd Conférence radar, Albuquerque, NM. American Meteorological Society.
- Fabry, Frédéric (August 2010). . Examples of remote-sensed data by instrument. J.S. Marshall Radar Observatory. Retrieved 2010-06-14.
- Harasti, Paul R.; McAdie, Colin J.; Dodge, Peter P.; Lee, Wen-Chau; Tuttle, John; Murillo, Shirley T.; Marks, Frank D., Jr. (April 2004). . Weather and Forecasting. American Meteorological Society. 19 (2): 219–239. Bibcode:. doi:. Retrieved 2009-06-09.
- (PDF). IRIS Product & Display Manual : Configuring IRIS Products. SIGMET. November 2004. Retrieved 2009-06-09.
- National Weather Service. . National Oceanic and Atmospheric Administration. Retrieved 2011-07-12.
- , National Weather Service (Texas Geographic Society – 2007)
- National Weather Service (31 January 2011). . Jetstream Online School for Weather. National Oceanic and Atmospheric Administration. Retrieved 2011-07-12.
- . Environment Canada. Retrieved 2007-06-23.
- Herbster, Chris (3 September 2008). . Introduction to NEXRAD Anomalies. Embry-Riddle Aeronautical University. Retrieved 2010-10-11.
- Diana Yates (2008). University of Illinois at Urbana – Champaign. Retrieved 2009-04-26
- Bart Geerts and Dave Leon (2003). University of Wyoming. Retrieved 2009-04-26
- Thomas A. Niziol (1998). Eastern Region WSR-88D Operations Note No. 12, August 1998. Retrieved 2009-04-26
- National Weather Service Office, Buffalo NY (8 June 2009). . National Oceanic and Atmospheric Administration. Retrieved 2009-09-01.
- Lammers, Dirk (29 August 2009). . Houston Chronicle. Associated Press. Retrieved 2009-09-01.
- Testud, J.; Le Bouar, E.; Obligis, E.; Ali-Mehenni, M. (2000). (pdf). J. Atmos. Oceanic Technol. 17 (3): 332–356. Bibcode:. doi:.
- Vulpiani, G.; Tabary, P.; Parent-du-Chatelet, J.; Marzano, F. S. (2008). (pdf). J. Atmos. Oceanic Technol. 25 (7): 1118–1135. Bibcode:. doi:.
- Carey, L. D.; Rutledge, S. A.; Ahijevych, D. A.; Keenan, T. D. (2000). (pdf). J. Appl. Meteor. 39 (9): 1405–1433. Bibcode:. doi:.
- Lemon, Leslie R. (June 1998). . Weather and Forecasting. 13 (2): 327–340. Bibcode:. doi:. ISSN . Retrieved 2011-05-25.
- David, McLaughlin; et al. (December 2009). . Bulletin of the American Meteorological Society. Boston MA: American Meteorological Society. 90 (12): 1797–1817. Bibcode:. doi:. ISSN . Retrieved 2010-08-31.
- . American Meteorological Society. 2005. Retrieved 2010-08-31.
- . Radar Operations Center. National Oceanic and Atmospheric Administration.
- WDT Support (July 7, 2015). . Radarscope.
- National Severe Storms Laboratory. (PDF). National Oceanic and Atmospheric Administration. Retrieved 2009-09-29.
- Otsuka, Shigenori; Tuerhong, Gulanbaier; Kikuchi, Ryota; Kitano, Yoshikazu; Taniguchi, Yusuke; Ruiz, Juan Jose; Satoh, Shinsuke; Ushio, Tomoo; Miyoshi, Takemasa (February 2016). (PDF). Weather and Forecasting. AMS: 329–340. Bibcode:. doi:. ISSN . Retrieved 2017-07-05.
- Bendix Corporation. Avionics Division. RDR-1200 Weather Radar System. Rev. Jul/73 ed. Fort Lauderdale: Bendix, Avionics Division, 1973.
- Barr, James C. Airborne Weather Radar. 1st ed. Ames: Iowa State UP, 1993.
- . IntelliWeather Inc. 2008. Retrieved 2011-11-26.
- . srh.noaa.gov. Retrieved November 9, 2015.
- Diana Yates (2008). . Urbana – Champaign, IL.: University of Illinois. Retrieved November 9, 2015.
- . Smithsonian.com. Retrieved November 9, 2015.
- . aba blog. Retrieved November 9, 2015.
- . Monarch-Butterfly.com. Retrieved November 9, 2015.
- Diehl, Robert H.; Larkin, Ronald P.; Black, John E. (April 2003). (pdf). The Auk. The American Ornithologists' Union. 120 (2): 278–290. doi:. ISSN . Retrieved November 9, 2015.
- Gagnon, François; Bélisle, Marc; Ibarzabal, Jacques; Vaillancourt, Pierre; Savard, Jean-Pierre L. (January 2010). (pdf). The Auk. The American Ornithologists' Union. 127 (1): 119–128. doi:. ISSN . Retrieved November 9, 2015.
- . /www.flysafe-birdtam.eu. Retrieved November 9, 2015..
- . American Meteor Society. Retrieved 2017-02-28.
- . www.lpi.usra.edu. Retrieved 2017-02-28.
- Fries, Marc; Fries, Jeffrey (2010-09-01). . Meteoritics & Planetary Science. 45 (9): 1476–1487. Bibcode:. doi:. ISSN .
- Brown, P.; McCAUSLAND, P. J. A.; Fries, M.; Silber, E.; Edwards, W. N.; Wong, D. K.; Weryk, R. J.; Fries, J.; Krzeminski, Z. (2011-03-01). . Meteoritics & Planetary Science. 46 (3): 339–363. Bibcode:. doi:. ISSN .
- Jenniskens, Peter; Fries, Marc D.; Yin, Qing-Zhu; Zolensky, Michael; Krot, Alexander N.; Sandford, Scott A.; Sears, Derek; Beauford, Robert; Ebel, Denton S. (2012-12-21). . Science. 338 (6114): 1583–1587. Bibcode:. doi:. ISSN . PMID .
- Fries, M. D.; Fries, J. A. (2010-09-01). . Meteoritics and Planetary Science Supplement. 73: 5365. Bibcode:.
- Fries, M.; Fries, J. (2010-03-01). . 41: 1179.
- Fries, M.; Fries, J.; Schaefer, J. (2011-03-01). . 42: 1130.