# Wavelet Filter Approach andZ-R Relationship in Meteorological Forecasting

Authors: Wallop Jiwlong; Anongrit Kangrang
Journal: Advances in Meteorology
 (2015)

---

## Abstract

The purpose of this study is to investigate theZ-R relationship for computing rainfall using conventional and wavelet filters technique. Wavelet filter technique was applied to data filtration process. The proposed model was applied to determine the rainfall of five rain gauge meteorological stations in Thailand. The three-hourly rainfall and radar reflectivity data were used in this study. The results indicated that the accumulative rainfall of wavelet filters technique was close to the observed rainfall data more than the results of conventional practice for both calibration and validation processes. Consequently, we are confident that a wavelet filters technique is a useful tool for estimating the rainfall.


---

## Body

## 1. Introduction

One of the most important attributes to flood mitigation and water resources management is accurate computation of precipitation over the given areas. Rainfall data is the most important meteorological parameter in hydrological studies and the crop-water requirement, because it constitutes the base of all the computations in water resources management and decision making problems. The obtained rainfall data from rain gauge is the point of rainfall that has been recorded in a particular time scale.Therefore, the lack of necessary climatological data made it difficult to obtain the reliable rainfall. In addition, the meteorological radar is a very powerful and useful tool for hydrological applications, since it can provide estimates of rainfall fields with a high spatial and temporal resolution over large areas. In the recent past, a few research studies have attempted to forecast the occurrence of precipitation. In measuring rainfall by radar,Z-R relationships are widely used to convert radar measured reflectivity to rainfall intensity; hence the accuracy of the estimation of Z-R relationship is important [1, 2]. In [3] applied a regression analysis technique to determine the relationship of synchronous datasets between measured rainfall intensity by rain gauge and measured or effective reflectivity by weather surveillance radar (Ze) at the pixel over the rain. However, in reality perfect synchronization between Ze and R is unachievable, except at the closest range and nearest to the ground [4].Wavelet analysis has been developed in order to provide a performing analysing tool for this kind of signal. It effectively renders possible a time-scale localisation of the process thanks to a projection on a class of functions which in turn makes it possible to extract data within a local neighbourhood. Being a signal’s time-scale analysis method, wavelet transform is better than Fourier transform due to its high-resolution features, good localization both in time and scale domains, and its capacity of analyzing signals at multi-time scales [5, 6]. In addition, wavelet analysis also has the capability of diagnosing the jump features of signals. Therefore, wavelet analysis has been applied to the studies of climate multi-time scales analysis. In the previous studies, wavelet transform is mainly used to analyze the multi-time scales characteristics of annual temperature and rainfall, and many interesting results have been found [7–9]. Basically, the wavelet decomposition uses a pair of filters to decompose iteratively the original time series such as rainfall and reflectivity. It results in a hierarchy of new time series that are easier to model and predict. These filters must satisfy some constraints such as causality and information lossless [10].In this study, the rainfall events of five meteorological stations in Thailand were analyzed forZ-R relationship. The Buffer Probability Technique (BPT) was developed based on GIS basic function and Probability Matching Method. The Z-R pairs can be analyzed for empirical power relationship. The main objective of this study is the comparison between conventional and wavelet filter technique in estimating rainfall using Z-R relationship; the focus was on two different models compared with observed data.


## 2. Materials and Methods

The selected three-hourly rainfall data of five rain gauge stations from Thai Meteorological Department and the radar reflectivity data from the Department of Royal Rainmaking and Agricultural Aviation at radar Pimai station, Nakhon Ratchasima province, were used in the study. The meteorological stations were sta.387401 Mahasarakham, sta.403201 Chaiyaphum, sta.405201 Roiet, sta.431201 Nakhon Ratchasima, and sta.436201 Buriram, most of which are located over Chi-Mun basin in the Northeastern region of Thailand. These locations are presented in Table1 and Figure 1. As the radar reflectivity, Z (mm6/mm3), commonly varies across many orders of magnitude, therefore Z used in this study is the reflectivity expressed in terms of dBz:(1)ZdBz=10log⁡Zmm6/mm3.Table 1
Locations of the five meteorological stations.
Stations
Latitude
Longitude
387401 Mahasarakham
16.247°N
103.068°E
403201 Chaiyaphum
15.800°N
102.033°E
405201 Roiet
16.050°N
103.683°E
431201 Nakhon Ratchasima
14.963°N
102.077°E
436201 Buriram
15.217°N
103.233°EFigure 1
Locations of the five meteorological stations.Radar reflectivity was obtained by an image analysis process from rainfall events that occurred in Northeastern region of Thailand, from June to August 2009 for long rainfall. The reflectivity was recorded from the Department of Royal Rainmaking and Agricultural Aviation at Pimai station, Nakhon Ratchasima Province, which corresponds to 2.5 km of CAPPI radar products at Pimai site from S-band polarmetric radar that transmits radiation with a wavelength of 10.7 cm and produces a beam width of 1.2 degrees with a maximum range of 480 km, as illustrated in Table2. The CAPPI products have six minutes and 1 sq.km temporal and spatial resolutions, respectively. This study assumes that there is no bias caused by the bright band effect and different observation altitudes at 1.5 km [11]. The CAPPI data lies within 240 km from the radar and the reflectivity values less than 10 dBz and greater than 50 dBz were excluded from the analysis to avoid the effect of noise in the measured radar reflectivity [12].Table 2
Characteristics of radar at Pimai station, Nakhon Ratchasima Province.
Details of radar
Characteristics
(i) Type of  radar
Doppler weather surveillanceRadar model DWSR-8500 S, S band
(ii) Wave length: cm
10.7
(iii) Beam width: degree
1.2
(iv) Pulse length: microsecond
0.8
(v) Resolution of record data
1 degree × 1 degree × 1 km
(vi) Maximum transmission power: Kw
850
(vii) Maximum range: km
480
(viii) Sequence of elevation angles
Operation A: 0.8, 1.7, 2.5Operation B: 3.4, 4.2, 5.1, 6.0, 7.4, 9.2, 11.6, 14.8, 18.4, 22.0


### 2.1.Z-R Relationship

The widely accepted relationship of the radar reflectivity factorZ and rainfall rate R is given by the empirical power law relationship [13] which can be defined by Z=aRb, where a and b are coefficients that depend on location and difference in climatology such as season and type of rain. The parameters a and b are related to the intercept and slope of the best-fit line through a plot of R versus Z on a log-log plot. Thus, if Z value can be measured and thus known, then R value can be found. Many researchers suggested that parameter b does not need to be varied as much as the parameter a [14, 15]. The most common values for a and b are 200 and 1.6, respectively. Reference [11] found that the existing Z-R relationship that has been used for the Pimai radar site is Z=294R1.33.


### 2.2.Z-R Matching Technique

In this study,Z-R matching analysis of Buffer Probability Technique (BPT) was developed based on GIS buffer technique and the concept of Probability Matching Method (PMM) to estimate Z-R relationship; however, this Z-R matching technique assumes that the raindrops fall absolutely vertical from the atmosphere to the rain gauge and radar reflectivity (Z) at time t is related to ground rainfall intensity (R) at time t too. The principal of relationship between radar reflectivity and ground surface rainfall diagram is shown in Figure 2(a). And the probability of cumulative density function between radar reflectivity and ground surface rainfall is matched as Figure 2(b); this means that the radar reflectivity has the same probability of occurrence as the gauge measured rain intensity [2, 16]. From past researches, several matching processes between reflectivity (Z) and rainfall (R) had been proposed such as Traditional Matching Method (TMM), Probability Matching Method (PMM), Window Probability Matching Method (WPMM), and Window Correlation Matching Method (WCMM). In all processes except TMM, PMM is the process that the sampling volume timing and location problem are not taken into account [1]. The matching is done between the Cumulative Distribution Functions (CDFs) of Z and R, described by the following formula:(2)∫Ri∞PRdR=∫Zi∞PZdZ,where P(R) is the probability density function of rainfall from gauge measurement and P(Z) is the probability density function of radar reflectivity. The matching of Z and R is done at the same probability level. According to the fact that the analysis is done on frequency domain, the timing errors are then eliminated in this method.Figure 2
Relationship between (a) radar reflectivity and ground surface rainfall (b) cumulative density function ofZ-R matching method.
(a)
(b)


### 2.3. Wavelet Transform for Signal Decomposition

Generally, a signal or functionf(t) can be expressed in linear decomposition by [17](3)ft=∑lalψlt,where l is an integer index for the finite or infinite sum, al are the real valued expansion coefficients, and ψl(t) are the set of real valued functions of t called the expansion set. For the wavelet transform, the two-parameter system is constructed and (3) becomes(4)ft=∑k∑jaj,kψj,kt,where both j and k are integer indices and ψj,k(t) are the wavelet expansion functions that usually form an orthogonal basis. The set of expansion coefficients aj,k are called the discrete wavelet transform (DWT) of f(t) and it is the inverse transform. The wavelet system is a two-dimensional expansion set (basis) for some class of one-dimensional signal. The wavelet expansion provides a time-frequency localization of the signal. First generation wavelet systems are generated from a single scaling function (wavelet) by simple scaling and translation. The two-dimensional parameterization is achieved from the function (mother wave) ψ(t) by(5)ψj,kt=2j/2ψ2jt-kj,k∈Z,where Z is the set of all integers and the factor 2j/2 maintains a constant norm independent of scale j. This parameterization of the time or space location by k and the frequency or the logarithm of scale by j turns out to be extraordinarily effective.


### 2.4. Haar Scaling Function and Wavelets

The multiresolution formulation needs two closely related basic functions. In addition to the waveletψ(t) that has been discussed and another basic function called the scaling function ϕ(t). The simplest possible orthogonal wavelet system is generated from the Haar scaling function and wavelet. These are shown in Figure 3.Figure 3
Haar scaling function and wavelet.
(a)
ϕ
(
t
)
(b)
φ
(
t
)In order to generate a set of expansion functions, the signal can be represented by the series in(6)ft=∑j,kaj,k2j/kψ2jt-kor using (5), as (7)ft=∑j,kaj,kψj,kt,where the two-dimensional set of coefficients aj,k is called the discrete wavelet transform (DWT) of f(t).An efficient way to implement DWT is to use a filter process. Normally, the low frequency content of the signal (approximation, A) is the most important part. It demonstrates the signal identity. The high frequency component (detail, D) is nuance. The original signal, S, passes through two complementary filters and emerges as two signals of A and D.The concept of the wavelet filters technique can be easily applied for signal analysis; this study used the technique to decompose the details (D) from the approximations (A) of rainfall records (y). In wavelet analysis, the approximations are the high-scale, low frequency components of the signal, and the details are low-scale, high frequency components. The decomposition process can be iterated with successive approximation being decomposed in turn, so that one signal is broken down into many lower resolution components. The wavelet decomposition can yield valuable information about signal. The approximations components are processed by theZ-R relationship for rainfall prediction.In this study the 1D Haar wavelet filters were applied to a rainfall and radar reflectivity series in order to denoise the original signal into the low frequency approximations (A) and the high frequency details (D). The approximations which are the smoothed signal represent the trend of the rainfall and radar reflectivity series whereas the details represent the nuance of them.


## 2.1.Z-R Relationship

The widely accepted relationship of the radar reflectivity factorZ and rainfall rate R is given by the empirical power law relationship [13] which can be defined by Z=aRb, where a and b are coefficients that depend on location and difference in climatology such as season and type of rain. The parameters a and b are related to the intercept and slope of the best-fit line through a plot of R versus Z on a log-log plot. Thus, if Z value can be measured and thus known, then R value can be found. Many researchers suggested that parameter b does not need to be varied as much as the parameter a [14, 15]. The most common values for a and b are 200 and 1.6, respectively. Reference [11] found that the existing Z-R relationship that has been used for the Pimai radar site is Z=294R1.33.


## 2.2.Z-R Matching Technique

In this study,Z-R matching analysis of Buffer Probability Technique (BPT) was developed based on GIS buffer technique and the concept of Probability Matching Method (PMM) to estimate Z-R relationship; however, this Z-R matching technique assumes that the raindrops fall absolutely vertical from the atmosphere to the rain gauge and radar reflectivity (Z) at time t is related to ground rainfall intensity (R) at time t too. The principal of relationship between radar reflectivity and ground surface rainfall diagram is shown in Figure 2(a). And the probability of cumulative density function between radar reflectivity and ground surface rainfall is matched as Figure 2(b); this means that the radar reflectivity has the same probability of occurrence as the gauge measured rain intensity [2, 16]. From past researches, several matching processes between reflectivity (Z) and rainfall (R) had been proposed such as Traditional Matching Method (TMM), Probability Matching Method (PMM), Window Probability Matching Method (WPMM), and Window Correlation Matching Method (WCMM). In all processes except TMM, PMM is the process that the sampling volume timing and location problem are not taken into account [1]. The matching is done between the Cumulative Distribution Functions (CDFs) of Z and R, described by the following formula:(2)∫Ri∞PRdR=∫Zi∞PZdZ,where P(R) is the probability density function of rainfall from gauge measurement and P(Z) is the probability density function of radar reflectivity. The matching of Z and R is done at the same probability level. According to the fact that the analysis is done on frequency domain, the timing errors are then eliminated in this method.Figure 2
Relationship between (a) radar reflectivity and ground surface rainfall (b) cumulative density function ofZ-R matching method.
(a)
(b)


## 2.3. Wavelet Transform for Signal Decomposition

Generally, a signal or functionf(t) can be expressed in linear decomposition by [17](3)ft=∑lalψlt,where l is an integer index for the finite or infinite sum, al are the real valued expansion coefficients, and ψl(t) are the set of real valued functions of t called the expansion set. For the wavelet transform, the two-parameter system is constructed and (3) becomes(4)ft=∑k∑jaj,kψj,kt,where both j and k are integer indices and ψj,k(t) are the wavelet expansion functions that usually form an orthogonal basis. The set of expansion coefficients aj,k are called the discrete wavelet transform (DWT) of f(t) and it is the inverse transform. The wavelet system is a two-dimensional expansion set (basis) for some class of one-dimensional signal. The wavelet expansion provides a time-frequency localization of the signal. First generation wavelet systems are generated from a single scaling function (wavelet) by simple scaling and translation. The two-dimensional parameterization is achieved from the function (mother wave) ψ(t) by(5)ψj,kt=2j/2ψ2jt-kj,k∈Z,where Z is the set of all integers and the factor 2j/2 maintains a constant norm independent of scale j. This parameterization of the time or space location by k and the frequency or the logarithm of scale by j turns out to be extraordinarily effective.


## 2.4. Haar Scaling Function and Wavelets

The multiresolution formulation needs two closely related basic functions. In addition to the waveletψ(t) that has been discussed and another basic function called the scaling function ϕ(t). The simplest possible orthogonal wavelet system is generated from the Haar scaling function and wavelet. These are shown in Figure 3.Figure 3
Haar scaling function and wavelet.
(a)
ϕ
(
t
)
(b)
φ
(
t
)In order to generate a set of expansion functions, the signal can be represented by the series in(6)ft=∑j,kaj,k2j/kψ2jt-kor using (5), as (7)ft=∑j,kaj,kψj,kt,where the two-dimensional set of coefficients aj,k is called the discrete wavelet transform (DWT) of f(t).An efficient way to implement DWT is to use a filter process. Normally, the low frequency content of the signal (approximation, A) is the most important part. It demonstrates the signal identity. The high frequency component (detail, D) is nuance. The original signal, S, passes through two complementary filters and emerges as two signals of A and D.The concept of the wavelet filters technique can be easily applied for signal analysis; this study used the technique to decompose the details (D) from the approximations (A) of rainfall records (y). In wavelet analysis, the approximations are the high-scale, low frequency components of the signal, and the details are low-scale, high frequency components. The decomposition process can be iterated with successive approximation being decomposed in turn, so that one signal is broken down into many lower resolution components. The wavelet decomposition can yield valuable information about signal. The approximations components are processed by theZ-R relationship for rainfall prediction.In this study the 1D Haar wavelet filters were applied to a rainfall and radar reflectivity series in order to denoise the original signal into the low frequency approximations (A) and the high frequency details (D). The approximations which are the smoothed signal represent the trend of the rainfall and radar reflectivity series whereas the details represent the nuance of them.


## 3. Methodology

The methodologies for this research were as follows:(1)
The quality control of radar image was performed; the errors due to radial anomality, ANAPROP, and other signals that did not represent reflectivity and errors caused by electronic problems investigation that had been done for both conventional and wavelet filter technique were removed from the measured reflectivity data.(2)
Buffer over rain gauge by Buffer Probability Technique transform process based on the buffer basic function of GIS is constructed to estimate the radar reflectivity over the rain gauge station at the same time.(3)
AverageZ value in Buffer zone and defined nonzero Z-R pairs are determined.(4)
The radar reflectivity was normalized into three-hourly data with the same time interval as rainfall measurement.(5)
From Cumulative Distribution Functions (CDFs), the radar reflectivity and rain rate at the same percentile can be obtained. These matching pairs were then used forZ-R relationship evaluation by conventional and wavelet filter technique.(6)
Calibration model by using rainfall event during 1 June–31 August 2009 and verification model by using rainfall event during 1 June–31 October 2011 were obtained.TheZ-R relationship with the studied process is shown in Figure 4.Figure 4
Studies process for rainfall estimating byZ-R relationship comparison between the conventional and wavelet filter technique.


## 4. Results

Three-hourly rainfall data and radar reflectivity at five stations in the upper Chi-Mun basin of Thailand during the period from June to August 2009 are filtered by Haar wavelet function as shown in Figures5 and 6.Figure 5
Rainfall data filtration by Haar wavelet filter type for each station.
(a)
Sta.387401 Mahasarakham
(b)
Sta.403201 Chaiyaphum
(c)
Sta.405201 Roiet
(d)
Sta.431201 Nakhon Ratchasima
(e)
Sta.436201 BuriramFigure 6
Radar reflectivity data filtration by Haar wavelet filter type for each station.
(a)
Sta.387401 Mahasarakham
(b)
Sta.403201 Chaiyaphum
(c)
Sta.405201 Roiet
(d)
Sta.431201 Nakhon Ratchasima
(e)
Sta.436201 Buriram


### 4.1. Model Calibration

The results of model calibration showed in Table3 and Figure 7 that the coefficients of Z-R relationship from the conventional technique were in the range from 218.50 to 376.00, being larger at sta.436201 Buriram and sta.403201 Chaiyaphum. The b coefficient did not vary significantly, remaining in the range from 1.025 to 1.491. The highest values occurred at sta.405201 Roiet and sta.431201 Nakhon Ratchasima. The coefficients of determination were all above 0.834. For the data filtration by wavelet technique, the values of the coefficient varied widely from 132.82 to 156.02 whereas the b coefficient values were high, ranging from 1.183 to 1.558. The coefficients of determination were all above 0.951 suggesting a good Z-R relationship. However, this accepted model is necessary for verification with other data for evaluating the performance of the model.Table 3
Result ofZ-R relationship for each station during 1 June–31 August 2009 period calibration.
Station
Conventional technique
Wavelet filter technique
Z-R relationship
R
2
Z
W
F
-
R
W
F relationship
R
2
387401 Mahasarakham
Z
=
246.46
R
1.2160
0.928
Z
W
F
=
139.46
R
W
F
1.3484
0.987
403201 Chaiyaphum
Z
=
254.10
R
1.0253
0.834
Z
W
F
=
135.77
R
W
F
1.1827
0.951
405201 Roiet
Z
=
218.53
R
1.4909
0.985
Z
W
F
=
132.82
R
W
F
1.5580
0.982
431201 Nakhon Ratchasima
Z
=
219.11
R
1.2491
0.972
Z
W
F
=
137.19
R
W
F
1.2546
0.987
436201 Buriram
Z
=
376.00
R
1.0660
0.947
Z
W
F
=
156.02
R
W
F
1.2768
0.974Figure 7
Scatter plot ofZ-R pairs for conventional practice and ZWF-RWF pairs for wavelet filters technique for each station.


### 4.2. Model Verification

The calibrated model parameter from the previous section was further validated using another climatological data set that occurred from 1 June 2011 to 31 October 2011, which were not considered in calibrating process. These data of five stations were used to calculate the rain rate; ifZ value can be measured and thus known, then rain rate or R value can be found. The verification results and performance statistics for Z-R relationship at each station were presented in Table 4. The correlation coefficients for conventional and wavelet filters technique were all above 0.892 and 0.926, respectively, whereas the standard error of estimate (SEE) values were ranging from 1.723 to 4.464 for conventional technique and 1.393 to 3.745 for wavelet filters technique. However, the cumulative observed rainfall data are more closed to wavelet filters technique than the conventional technique as show in Figures 8 to 17.Table 4
Result of statisticalZ-R relationship for each station during June–October 2011 period verification.
Station
Conventional technique
Wavelet filter technique
Z-R relationship
r
SEE
Z-R relationship
r
SEE
387401 Mahasarakham
Z
=
246.46
R
1.2160
0.959
2.926
Z
W
F
=
139.46
R
W
F
1.3484
0.966
2.671
403201 Chaiyaphum
Z
=
254.10
R
1.0253
0.892
4.464
Z
W
F
=
135.77
R
W
F
1.1827
0.926
3.745
405201 Roiet
Z
=
218.53
R
1.4909
0.973
2.054
Z
W
F
=
132.82
R
W
F
1.5580
0.987
1.444
431201 Nakhon Ratchasima
Z
=
219.11
R
1.2491
0.978
1.944
Z
W
F
=
137.19
R
W
F
1.2546
0.989
1.393
436201 Buriram
Z
=
376.00
R
1.0660
0.977
1.723
Z
W
F
=
156.02
R
W
F
1.2768
0.984
1.436Figure 8
Time series plot of mean observed rainfall at sta.387401 Mahasarakham compared with conventional and wavelet filters technique.Figure 9
Cumulative observed rainfall data at sta.387401 Mahasarakham during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.Figure 10
Time series plot of mean observed rainfall at sta.403201 Chaiyaphum compared with conventional and wavelet filters technique.Figure 11
Cumulative observed rainfall data at sta.403201 Chaiyaphum during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.Figure 12
Time series plot of mean observed rainfall at sta.405201 Roiet compared with conventional and wavelet filters technique.Figure 13
Cumulative observed rainfall data at sta.405201 Roiet during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.Figure 14
Time series plot of mean observed rainfall at sta.431201 Nakhon Ratchasima compared with conventional and wavelet filters technique.Figure 15
Cumulative observed rainfall data at sta.431201 Nakhon Ratchasima during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.Figure 16
Time series plot of mean observed rainfall at sta.436201 Buriram compared with conventional and wavelet filters technique.Figure 17
Cumulative observed rainfall data at sta.436201 Buriram during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.This clearly indicates that theZWF-RWF relationship by wavelet filter technique developed in the present study can be used with confidence in estimating rainfall for this study area. The results show that the overall correlation coefficients between the estimated and observed rainfall for all station were high range, which are acceptable.


## 4.1. Model Calibration

The results of model calibration showed in Table3 and Figure 7 that the coefficients of Z-R relationship from the conventional technique were in the range from 218.50 to 376.00, being larger at sta.436201 Buriram and sta.403201 Chaiyaphum. The b coefficient did not vary significantly, remaining in the range from 1.025 to 1.491. The highest values occurred at sta.405201 Roiet and sta.431201 Nakhon Ratchasima. The coefficients of determination were all above 0.834. For the data filtration by wavelet technique, the values of the coefficient varied widely from 132.82 to 156.02 whereas the b coefficient values were high, ranging from 1.183 to 1.558. The coefficients of determination were all above 0.951 suggesting a good Z-R relationship. However, this accepted model is necessary for verification with other data for evaluating the performance of the model.Table 3
Result ofZ-R relationship for each station during 1 June–31 August 2009 period calibration.
Station
Conventional technique
Wavelet filter technique
Z-R relationship
R
2
Z
W
F
-
R
W
F relationship
R
2
387401 Mahasarakham
Z
=
246.46
R
1.2160
0.928
Z
W
F
=
139.46
R
W
F
1.3484
0.987
403201 Chaiyaphum
Z
=
254.10
R
1.0253
0.834
Z
W
F
=
135.77
R
W
F
1.1827
0.951
405201 Roiet
Z
=
218.53
R
1.4909
0.985
Z
W
F
=
132.82
R
W
F
1.5580
0.982
431201 Nakhon Ratchasima
Z
=
219.11
R
1.2491
0.972
Z
W
F
=
137.19
R
W
F
1.2546
0.987
436201 Buriram
Z
=
376.00
R
1.0660
0.947
Z
W
F
=
156.02
R
W
F
1.2768
0.974Figure 7
Scatter plot ofZ-R pairs for conventional practice and ZWF-RWF pairs for wavelet filters technique for each station.


## 4.2. Model Verification

The calibrated model parameter from the previous section was further validated using another climatological data set that occurred from 1 June 2011 to 31 October 2011, which were not considered in calibrating process. These data of five stations were used to calculate the rain rate; ifZ value can be measured and thus known, then rain rate or R value can be found. The verification results and performance statistics for Z-R relationship at each station were presented in Table 4. The correlation coefficients for conventional and wavelet filters technique were all above 0.892 and 0.926, respectively, whereas the standard error of estimate (SEE) values were ranging from 1.723 to 4.464 for conventional technique and 1.393 to 3.745 for wavelet filters technique. However, the cumulative observed rainfall data are more closed to wavelet filters technique than the conventional technique as show in Figures 8 to 17.Table 4
Result of statisticalZ-R relationship for each station during June–October 2011 period verification.
Station
Conventional technique
Wavelet filter technique
Z-R relationship
r
SEE
Z-R relationship
r
SEE
387401 Mahasarakham
Z
=
246.46
R
1.2160
0.959
2.926
Z
W
F
=
139.46
R
W
F
1.3484
0.966
2.671
403201 Chaiyaphum
Z
=
254.10
R
1.0253
0.892
4.464
Z
W
F
=
135.77
R
W
F
1.1827
0.926
3.745
405201 Roiet
Z
=
218.53
R
1.4909
0.973
2.054
Z
W
F
=
132.82
R
W
F
1.5580
0.987
1.444
431201 Nakhon Ratchasima
Z
=
219.11
R
1.2491
0.978
1.944
Z
W
F
=
137.19
R
W
F
1.2546
0.989
1.393
436201 Buriram
Z
=
376.00
R
1.0660
0.977
1.723
Z
W
F
=
156.02
R
W
F
1.2768
0.984
1.436Figure 8
Time series plot of mean observed rainfall at sta.387401 Mahasarakham compared with conventional and wavelet filters technique.Figure 9
Cumulative observed rainfall data at sta.387401 Mahasarakham during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.Figure 10
Time series plot of mean observed rainfall at sta.403201 Chaiyaphum compared with conventional and wavelet filters technique.Figure 11
Cumulative observed rainfall data at sta.403201 Chaiyaphum during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.Figure 12
Time series plot of mean observed rainfall at sta.405201 Roiet compared with conventional and wavelet filters technique.Figure 13
Cumulative observed rainfall data at sta.405201 Roiet during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.Figure 14
Time series plot of mean observed rainfall at sta.431201 Nakhon Ratchasima compared with conventional and wavelet filters technique.Figure 15
Cumulative observed rainfall data at sta.431201 Nakhon Ratchasima during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.Figure 16
Time series plot of mean observed rainfall at sta.436201 Buriram compared with conventional and wavelet filters technique.Figure 17
Cumulative observed rainfall data at sta.436201 Buriram during 1 June–31 October 2011 for verification compared with conventional and wavelet filters technique.This clearly indicates that theZWF-RWF relationship by wavelet filter technique developed in the present study can be used with confidence in estimating rainfall for this study area. The results show that the overall correlation coefficients between the estimated and observed rainfall for all station were high range, which are acceptable.


## 5. Discussion

This paper applied a wavelet filter technique to signify the data preprocessing as an important step in the construction ofZ-R relationship for rainfall estimation. The wavelet filters can extract the chaotic signal components from the original rainfall and dBz data for the calibration. Therefore, this study demonstrates that the importance of wavelet filters technique is a practical tool for the prediction of three-hourly rainfall in Upper Chi-Mun basin of Thailand.



---

*Source: 858717-2015-08-31.xml*