SLIDE 1 Near real-time radar reflectivity calibration and monitoring using ground clutter
Valentin Louf (Monash, BOM)
- A. Protat (BOM), C. Jakob (Monash)
Radar Workshop 2016, Melbourne
SLIDE 2 2 Ground clutter? How about no! Ground clutter is responsible for:
- Spurious signal.
- Bias reflectivity.
- Skew Doppler velocity
estimates. What we (traditionally) want:
ground clutters.
Radar Workshop 2016, Melbourne
SLIDE 3
Can a constructive use of ground clutter be made?
SLIDE 4 4 Can a constructive use of ground clutter be made?
(Rinehart and Frush, 1983).
- Estimation of the attenuation
caused by rainfall (Delrieu et al., 1997).
- Estimation of the refractive
index of air (Fabry, 2003).
- Estimation of the topographic
structure (Mesnard et al. 2003)
calibration, proposed by Rinehart in 1978.
Figure from Louf et al (2015)
Radar Workshop 2016, Melbourne
SLIDE 5 5 Monitoring of radar calibration.
1960: Atlas and Mossop 1978: Rinehart 2008: Silberstein et al. 2009: Marks et al. 2015: Wolff et al.
Not a new technique⦠but not a common one.
Radar Workshop 2016, Melbourne
SLIDE 6 6 How can we use ground clutters?
- The meteorological radar equation:
ππ ππΆπ = 10 log10 π
π + 20 log10 π β 10 log10 π·
β ππ: reflectivity, β π
π : echo power,
β π: the range. β π·: the weather radar constant.
- 10 log10 π· represents the radar sensitivity (in dB).
Radar Workshop 2016, Melbourne
SLIDE 7 7 How can we use ground clutters?
- Clutter is due almost entirely to human-made structure:
β Building microphysics do not change drastically with time:
π = 0
β It is unlikely that the range itself changes significantly over time:
- Ξ20 log10 π = 0
- Therefor:
- Ξππ = Ξ10 log10 π·
The variations of ground echoes reflectivity are directly linked to the variations in radar calibration.
Radar Workshop 2016, Melbourne
SLIDE 8 8 How can we use ground clutters?
β Select all PPIs for a day without precipitation (within the 5, or 10km of the radar). β Keep the echoes with a reflectivity above 50 dBZ and a frequency
- f occurrence superior to 50%.
- Apply the mask.
β Collects the clutter area reflectivity for a given moment/day/month. β Computes the probability density function (PDF) calculate the 95th percentile of reflectivity.
SLIDE 9 9 Computation of the Relative Calibration Adjustment (RCA) offset
- The RCA is defined as being:
ππ·π΅ ππΆ = π·πΈπΊ95π’β(ππππ‘πππππ) β π·πΈπΊ95π’β(ππππ£π’π’ππ )
- The RCA value is the adjustment (in dB) needed to obtain
agreement to the baseline.
SLIDE 10
10 Why using the cumulative distribution function?
SLIDE 11
The results for CPOL.
SLIDE 12
12 The RCA on CPOL: average seasons
SLIDE 13
13 The good
SLIDE 14
14 The bad
SLIDE 15
15 The ugly
A software update that was not suppose to change anything.
SLIDE 16
16 Back to normal
SLIDE 17
17 Long-term evolution of calibration
SLIDE 18 18 Β« Since it is unlikely that the terrain itself changes significantly over such periods of time [a day], and if other explanations such as variations in radar calibration or electromagnetic properties of the ground can be ruled out,
- ne must conclude that it is the medium in which the radar
wave travels that experiences these changes, in this case,
Fabry (2003) in Β« Meteorological Value of Ground Target Measurements by Radar Β»
SLIDE 19
19 Diurnal cycle of the RCA
SLIDE 20
What future for this technique?
SLIDE 21
21 Implementation on Indian radarβs network
Amar Jyothiβs work.
SLIDE 22
22 Test on operational radars (Melbourne).
SLIDE 23 23 Conclusions
- Conclusion on the radar technique:
β Calibration β Monitoring β Easy checking on large dataset.
- Conclusion on the atmospheric research:
β Clutters are not noise, they contain valuable informations!
SLIDE 24
The new CPOL dataset
SLIDE 25 25
β Raw data calibrated and filtered data. (ZH, ZC, ZDR, KDP, PHIDP, RHOHV, WIDTH)
- Polar. (Every 10 mins)
- Cartesian. (Daily mean)
- Level 2:
β 2D: Conv/strat, Z(2.5 km), rainfall, 0dBZ, cloud top height 17 dBZ, 20 dB, 40 dB
- 3D: 3D winds, hydrometeor
classification.
- Level 3:
- Climatology.
- Area mean rainfall, Peak
rainfall 99.5th percentile), Area fraction of convective and stratiform, Area fraction vs height, Stats on wind and convective mass fluxβ¦
The new CPOL dataset (delivery in January)