<br />374
<br />
<br />PJ. BRf::MAUD AND Y,B. POINTIN
<br />
<br />?n thei: motion as detected in successive radar pictures is only the first step
<br />m makmg a precipitation forecast, as the intensity of the extrapolated radar
<br />echoes must then be converted into rainfaII at the ground (Einfalt, 1991), This
<br />conversion can lead, for the required high spatial and temporal resolution, to
<br />larger uncertainties than for hourly averages (Messaoud and POintin, 1990),
<br />mamly because of the variability of the raindrop size distribution, This paper is
<br />concerned only with this first step of producing the best extrapolation method,
<br />The approaches used until now can conveniently be considered as falling into
<br />one of three categories: cross-correlation techniques. echo centroid methods
<br />and more complex methods (Elvander, 1976; ColJier, ] 978, 1989),
<br />Cross-correlation techniques were the first to be studied (Wilson and
<br />Kessler, 1963) because they are the simplest way to make an estimate of the
<br />motion of the radar echoes, The principle consists of overlaying, with a spatial
<br />shift, the radar picture observed at I an the picture observed at t - [JI and
<br />searching for the optimum relative position of the two pictures corresponding
<br />to the best value of the cross-correlation coefficient (Austin and Bellon, 1974),
<br />This optimum shift indicates the displacement of the radar echoes, provided
<br />that there afe no significant size, shape and intensity changes during the time
<br />interval Ot, The forecast is then made by extrapolating this detected motion
<br />and applying it unifonnly over the whole radar picture to the time I + 0/,
<br />This technique has been operationally used, with some success for the forecast
<br />of widespread precipitation over reasonably flat terrain, by Bellon and Austin
<br />(1978), but it does not alIow for ditTerent motions of individual radar echoes
<br />which can occur in the case of developing convective activity or orographic
<br />rainfall. As these are very important types of precipitation for urban
<br />hydrology, cross-correlation techniques are not the most appropriate for this
<br />purpose.
<br />In the case of echo centroid methods, the motion of each radar echo is
<br />deduced from the successive positions of its centroid (Barclay and Wilk, 1970;
<br />Wilk and Gray, 1970; Bjerkaas and Forsyth, 1980), Therefore, these methods
<br />are more elaborate than cross-correlation methods but they still remain too
<br />simple, Echo centroid methods perform well in the case of isolated radar
<br />echoes (e.g. isolated thunderstorms), but in other cases, such as convective
<br />cens embedded within a large stratiform rain area or within a frontal band,
<br />these methods may have problems (Zitt"], 1976), Indeed, the successive
<br />positions of the echo centroid are not always representative of the real echo
<br />motion when the morphological changes of the echo from one picture to the
<br />next are important, especially if splits or mergers of echoes OCCur. Because of
<br />these problems, echo centroid methods are also not well suited for short-term
<br />hydrological purposes, ,
<br />The third kind of technique, usually named the 'complex method', may
<br />
<br />~
<br />.' "
<br />:;
<br />!
<br />
<br />~
<br />
<br />~
<br />
<br />FORECASTING HEAVY RAINFALL FROM RAIN CELL MOTION
<br />
<br />315
<br />
<br />solve some of the problems of the cross-correlation and echo centroid
<br />methods, Blackmer et al. (1973) have used a local cross-correlation technique
<br />to estimate the local motion of the echoes which are defined by a fixed
<br />reflectivity threshold (afterwards called a 'T echo'), This estimation of the
<br />local echo motion leads to more realistic trajectories than those deduced from
<br />the successive positions of the centroid, because it is less sensitive to mOr-
<br />phological changes, Going further, Einfalt e,t ,al. (1990) d~fined 'structured
<br />echoes' which are the union of two or more SImple echoes (T echoes), The
<br />matching of these structured echoes with simple echoes permits the possible
<br />split or merger of simple echoes to be taken into account. Thi~ technique,
<br />which uses complex mathematical tools such as pattern recogmtlon, offers
<br />distinct improvements over earlier techniques for detailed rainfall forecasts at
<br />short time and spatial scales. Nevertheless, the main characteristics of most
<br />existing methods, especially the definition of the tracked,ntities (T echoes),
<br />are not based upon the physics of the cloud processes associated with strong
<br />radar echoes in contrast to the PARAPLUIE method which is described
<br />,
<br />below,
<br />
<br />PARAPLUIE: A HEAVY RAINFALL FORECASTING METHOD
<br />
<br />'The PARAPLUIE method (Bremaud, 1991) belongs to the complex
<br />method category, and is divided into four steps: (a) tracked entity definition;
<br />(b) tracked entity characterization; (c) tracked entity matching for the
<br />detection of the ditTerent motions; (d) forecast by extrapolation, The method
<br />is fully automated, so it does not require operator intervention.
<br />
<br />Tracked entity definition
<br />
<br />In contrast to a large majority of methods, which use the T echo defined by
<br />a connectivity property above a fixed given reflectivity threshold, the
<br />PARAPLUIE algorithm defines the tracked entities, called 'CEL echoes', by
<br />the contour constructed at 6 dB below each enclosed reflectivity maximum,
<br />Thus the reflectivity threshold defining the CEL echo is not fixed but varies
<br />with the rain type, This definition of our tracked entities, originally introduced
<br />by Crane (1979) for aircraft hazard identification, is illustrated in Figs, I and
<br />2 which show for a frontal rain band observed on 11 October 1988 in the
<br />Covennes are; of France (Andrieu ,et aI., 1989), a PPI map and a horizontal,
<br />reflectivity profile along an east-west line 10 km north of the radar, respec-
<br />tively, In Fig, 2, the small-scale variation of the reflectivity profile between 5
<br />and 20 km is due to ground echoes in this mountainous area. These ground
<br />echoes are detected by the PARAPLUIE algorithm both by the small size of
<br />
<br />
|