the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Citizen rain gauges improve hourly radar rainfall bias correction using a twostep Kalman filter
Monton Methaprayun
Thom Bogaard
Gerrit Schoups
MarieClaire Ten Veldhuis
The low density of conventional rain gauge networks is often a limiting factor for radar rainfall bias correction. Citizen rain gauges offer a promising opportunity to collect rainfall data at a higher spatial density. In this paper, hourly radar rainfall bias adjustment was applied using two different rain gauge networks: tipping buckets, measured by Thai Meteorological Department (TMD), and daily citizen rain gauges. The radar rainfall bias correction factor was sequentially updated based on TMD and citizen rain gauge data using a twostep Kalman filter to incorporate the two gauge datasets of contrasting quality. Radar reflectivity data from the Sattahip radar station, gauge rainfall data from the TMD, and data from citizen rain gauges located in the Tubma Basin, Thailand, were used in the analysis. Daily data from the citizen rain gauge network were downscaled to an hourly resolution based on temporal distribution patterns obtained from radar rainfall time series and the TMD gauge network. Results show that an improvement in radar rainfall estimates was achieved by including the downscaled citizen observations compared with bias correction based on the conventional rain gauge network alone. These outcomes emphasize the value of citizen rainfall observations for radar bias correction, in particular in regions where conventional rain gauge networks are sparse.
 Article
(7059 KB)  Fulltext XML
 BibTeX
 EndNote
Hydrometeorological hazards, like flash floods and landslides, cause severe damage to economies, properties, and human lives worldwide. In this context, flood forecasting and warning systems are a valuable nonstructural measure to mitigate damage. However, such systems require the input of rainfall data at a high spatial and temporal resolution. In most regions of the world, automatic rain gauge networks are insufficient for this purpose. Weather radar, which can better capture the variation in rainfall fields at fine spatial and temporal resolutions, could be used as an alternative rainfall product to improve the accuracy of flash flood estimates and warnings. (Collinge and Kirby, 1987; Sun et al., 2000; Uijlenhoet, 2001; Bedient et al., 2003; Creutin and Borga, 2003; Mapiam et al., 2009a, 2014; Mapiam and Chautsuk, 2018; Corral et al., 2019). However, weather radar provides an indirect measurement of backscattered electromagnetic waves called radar reflectivity data (Z), and quantitative estimation of radar rainfall data (R) is acknowledged as a complex process. Various sources of error affect radar rainfall estimates, mainly errors in reflectivity measurements and reflectivity–rainfall conversion (Jordan et al., 2000). Correction of these two sources of error is a crucial procedure to increase the accuracy of radar rainfall estimates. “Ground truthing” by rain gauge data is required to calibrate the Z–R relationship (Z=AR^{b}). The calibrated parameter A in the Z–R relationship will include any errors in the radar system caused by the electrical calibration of the radar (Seed et al., 2002).
The calibrated Z–R equation is used to convert the measured instantaneous reflectivity data to rainfall intensity and, subsequently, to accumulate them into the required temporal resolution. However, the parameters A and b vary significantly, even within a single storm event, depending on the rainfall characteristics, which can exhibit a highly dynamic raindrop size distribution (DSD) (Ulbrich, 1983: Smith et al., 2009). Additionally, past studies have found that the Z–R parameters are sensitive to the temporal resolution of the rain gauge rainfall data that are used for the Z–R calibration (Hitchfeld and Bordan, 1954; Smith et al., 1975; Wilson and Brandes, 1979; Klazura, 1981; Steiner et al., 1995; Mapiam and Sriwongsitanon, 2008; Mapiam et al., 2009b). Consequently, an important source of error remains associated with the Z–R conversion process (Jordan et al., 2000; Berne and Krajewski, 2013). Many researchers have attempted to correct this kind of error by classifying the measured reflectivity data into different storm types and, thereafter, constructing the Z–R equation corresponding to the classified storm characteristics (Joss and Waldvogel, 1970; Rogers, 1971; Battan, 1973; Klazura, 1981; Austin, 1987; Rosenfeld et al., 1992, 1993; Tokay and Short, 1996; Amitai, 2000; Arai et al., 2005; Fang et al., 2018). To account for the effect of using rain gauge data with different temporal resolutions on the Z–R relationships, Mapiam et al. (2009b) developed a universal scaling transformation function for converting the reference parameters A (obtained from using daily gauge rainfall data in the calibration) to the parameter A for subdaily resolutions. This improved the accuracy of the estimated subdaily radar rainfall, especially in locations with limited shortduration rain gauge measurements.
After Z–R conversion, bias is expected to remain between the assessed radar rainfall and the true rainfall amount at the rain gauge locations if a fixed Z–R relationship is used to estimate radar rainfall over the entire radar domain (Chumchean et al., 2006; Wang et al., 2015). An effective bias correction technique is key for enhancing the quality of radar rainfall estimates (Steiner et al., 1999) and for removing the residual errors between radar rainfall obtained from the Z–R relationship and from the rain gauge data. Mean field bias (MFB) adjustment is the conventional method used to obtain a static bias factor which assumes that the Z–R relationship is homogeneous in space but varies in time (Smith et al., 2007; Vieux and Bedient, 2004; Wilson, 1970). In this method, a multiplicative correction factor is applied uniformly across the radar coverage. As the MFB approach does not consider the noise and uncertainty of the rain gauge observations nor the spatial variability in observation bias, this can lead to large errors in radar rainfall estimates, particularly in areas where the density of rain gauge networks is limited. The Kalman filter (KF) is an efficient algorithm that has been applied to correct the spatially uniform mean field bias, especially in real time by accounting for the temporal variation in the mean bias as well as uncertainties in the ground rainfall measurements (Anhert et al., 1986; Smith and Krajewski, 1991; Anagnostou et al., 1998; Seo et al., 1999; Dinku et al., 2002; Chumchean et al., 2006).
Previous studies have used the KF to predict and correct the mean field bias in order to mitigate the observation error variances affecting the mean field bias estimate. Chumchean et al. (2006) found that the density of the rain gauge network also plays an important role in the radar rainfall bias adjustment. They found that lowering the density of rain gauge observations in the KF process reduced accuracy of radar rainfall estimates. Additionally, the KF approach outperforms the use of MFB if the rain gauge density is less than one gauge per 90 km^{2}, and both KF and MFB produce identical performance when the rain gauge density is greater than one gauge per 70 km^{2}.
In basins where a dense rainfall network is not available, citizen science (CS) offers a promising opportunity to enhance the density of rainfall observations (Davids et al., 2019). With the popularization of smartphones and the availability of (relatively) simple and cheap equipment, abundant mobile applications and projects have been initiated in water resources management to measure hydrometeorological variables like rainfall, water level height, or water quality as well as to ground truth remotely sensed information on, for example, land use (Srivastra et al., 2018; Davids et al., 2019; See, 2019; Seibert et al., 2019). In the current study, we focus on rainfall measured by local citizens using a network of cheap rain gauges and a specially designed mobile application. As citizen rainfall observations are typically provided at a daily scale, a temporal downscaling technique is needed for subdaily applications. A variety of temporal rainfall downscaling methods have been developed since the 1970s. The simplest approach is to distribute daily rainfall data to subdaily resolutions by assuming uniform distributions. Stochastically generating subperiod data or spatially transferring finerresolution rainfall from a nearby rain gauge station to the study area based on spatial correlations are alternative approaches (Koutsoyiannis, 2003; Debele et al., 2007). However, these methods are usually not designed for realtime data disaggregation over large areas. Instead, a common approach for such scenarios is to downscale daily rainfall based on a simple fraction technique by considering the distribution patterns of highresolution gridded rainfall products from radar or satellite sensors (Paulat et al., 2008; Wüest et al., 2009; Vormoor and Skaugen, 2013; Sideris et al., 2014; Barton et al., 2020). This study aimed to modify the KF logic by integrating hourly rain gauge data with daily citizen rain gauge data that are downscaled to an hourly timescale using a simple fraction method. The radar rainfall bias correction factor was sequentially updated using a twostep Kalman filter accounting for the contrasting quality of the hourly rain gauge data and downscaled citizen rain gauge data. The question that we set out to answer is as follows: to what extent do the downscaled citizen rainfall observations improve the accuracy of hourly radar rainfall estimates? Several scenarios of hourly rainfall distribution patterns were applied for downscaling in order to investigate the most suitable technique for hourly radar rainfall assessment. Tubma Basin, located in Rayong Province, eastern Thailand, was used as a case study area to test the approach.
2.1 Study area
The study area is the Tubma Basin (12.6789–12.8775^{∘} N, 101.0881–101.2975^{∘} E) located in Rayong Province, eastern Thailand (Fig. 1). It covers a catchment area of 197 km^{2} with a basin elevation ranging from 4 to 416 m a.m.s.l. (above mean sea level). The main river, the Klong, is 42 km in length; originates in the Khao Chom Hae, Khao Ket, and Khao Krabok mountains; and flows downstream to the northwest before meeting the Gulf of Thailand in Pak Nam, Mueang Rayong District. The Tubma watershed is susceptible to flooding, in particular Rayong. In Fig. 1, we show the climatological variation across the study area and its surroundings, based on the 30year (1987–2017) annual mean rainfall from a network of 311 daily rain gauges owned by the TMD and situated within 200 km of the Tubma Basin. Spatial rainfall patterns were generated using the inverse distance squared (IDS) method between the gauge locations. The map shows that while there is a small gradient in mean annual rainfall (1100–1700 mm mean annual rainfall) across the Rayong and Chonburi provinces (within 90 km of the study area), changes are more pronounced when the distance exceeds the 90 km boundary, especially to the east of the study area. This is because these areas are affected differently by the southwest monsoon. Consequently, the evaluation of the effectiveness of the bias correction techniques was carried out within a 90 km radius of the study area to ensure similar climatology.
2.2 Radar data
2.2.1 Reflectivity data collection
The Tubma Basin is located within the coverage of the Sattahip radar station. The Sattahip radar, which belongs to the Department of Royal Rainmaking and Agricultural Aviation (DRRAA), is an Sband Doppler radar that transmits radiation with a frequency of 2.9 GHz and a halfpower beamwidth of 1.0^{∘}. The radar reflectivity product is provided on a Cartesian grid covering a 240×240 km extent with a 0.6×0.6 km spatial and 6 min temporal resolution. The Sattahip radar provides the reflectivity data derived from the 2.5 km constant altitude plan position indicator (CAPPI). The CAPPI reflectivity data originate from below the climatological freezing level, so the effects of measurement errors caused by the bright band were considered negligible. The effects of ground clutter were removed from the reflectivity data by finding the clutter locations and discarding the radar measurements in these areas. Additionally, the noise power caused by various sources, such as emission from space (cosmic noise), radiation from electrical sources near the radar antenna, and internally generated noise, were eliminated by setting reflectivity values below 15 dBZ to zero (Doviak and Zrnic, 1984). While Fulton et al. (1998) suggested that measured reflectivity values that are greater than 53 dBZ be limited to 53 dBZ in radar rainfall estimation to mitigate false interpretation caused by hail, the hail cap can be seen as an adaptable threshold representing the maximum expected instantaneous rain rate, which is unfortunately quite difficult to determine for a particular storm. Note that slightly higher hail threshold values have also been reported in tropical environments. After data quality control, we separated the data into three datasets: the first dataset from May to October 2013 and from May to September 2014 was used for the climatological Z–R calibration; the second dataset from October 2014 was used for the Z–R verification; and the third dataset from August to October 2019 was used in the bias correction processes.
2.2.2 The Z–R calibration and radar rainfall aggregation
The Z–R conversion error is a crucial source of error in radar rainfall estimates. The Z–R relationship was used to convert the measured reflectivity data (Z, mm^{6} m^{−3}) into rainfall rates (R, mm h^{−1}). The Z–R calibration and verification are essential procedures to ascertain the A and b parameters in the relationship. Firstly, the instantaneous 6 min radar reflectivity was converted to rainfall intensity using the climatological relationship Z=200R^{1.6} proposed by Marshall and Palmer (1948). Secondly, the estimated 6 min initial instantaneous radar rainfall data were aggregated to an hourly rainfall resolution using the accumulation algorithm proposed by Fabry et al. (1994). Thirdly, gauge rainfall was aggregated to an hourly resolution. Finally, the optimal value of the parameter A was established by minimizing the mean absolute error (MAE) between the gauge and the radar rainfall estimates, whereas the exponent b was 1.5 in our study. This is because radar rainfall estimates are relatively insensitive to b with typical values between 1.2 and 1.8 (Battan, 1973; Ulbrich, 1983). Thus, 1.5 was generally suitable to represent the exponent b in the Z–R relation (Doelling et al., 1998; Steiner and Smith, 2000; Hagen and Yuter; 2003; Germann et al., 2006; Chantraket et al., 2016). The MAE is illustrated in Eq. (1):
where G_{i,t} is the gauge rainfall (mm h^{−1}) at gauge i for hour t, R_{i,t} is the hourly radar rainfall accumulation (mm) at the pixel corresponding to the ith rain gauge for hour t, N_{G} is the total number of rain gauges, and T is the total period used in the calculation. Results found that the optimal climatological Z–R relationship for the Sattahip radar used in this study is Z=251R^{1.5}. This relation is appropriate for both the calibration and validation datasets with an MAE of 1.36 and 1.47 mm, respectively.
2.3 Rain gauge data
2.3.1 Rainfall data collection
Data from the network of 297 continuous tippingbucket gauge stations located within the Sattahip radar radius were collected (Fig. 2). These rain gauges, which produce data with a temporal resolution of 15 min, are owned and operated by the Thai Meteorological Department (TMD). All of the continuous rain gauges used in this study have a resolution of 0.5 mm. The data quality screening was first carried out using the doublemass curve method. To avoid norainfall events and systematical underreporting by the tippingbucket rain gauge, hourly data above the tippingbucket resolution of 0.5 mm were selected in the next step. Rain gauges with more than 80 % of their recorded rainfall amounts below the 0.5 mm threshold at the daily scale were excluded from the analysis. It turns out that many of these faulty gauges recorded zero rainfall throughout most of the study period. Following the implementation of these criteria, rainfall data obtained from 134 rain gauges corresponding to the collected reflectivity datasets were used for the Z–R calibration and validation processes. For the bias adjustment computation, the selection of rain gauge networks with rainfall behavior similar to the study area is necessary. We selected 14 TMD rain gauges in the region surrounding the Tubma Basin (Rayong and Chonburi provinces), based a on spatial decorrelation analysis, for this the process.
2.3.2 Citizen rain observations
Out of the total TMD rain gauge network, only one rain gauge is located in the Tubma Basin. To increase the density of the rain gauge network in the basin, lowcost citizen rain gauges were implemented in this study to capture the spatial heterogeneity of rainfall in the region. A total of 16 citizen rain gauges were installed (Fig. 2) with local residents taking daily measurements. The additional 16 citizen rain gauges (with one station colocated with an existing TMD gauge) increased the density of rain gauges in the Tubma Basin to one gauge per ∼12 km^{2}. The citizen observations were made by installing a coneshaped transparent plastic rain gauge in an open area around a school, monastery, bridge or other building. This type of rain gauge is in standard use in South Africa (Fig. A1) and has a diameter of approximately 12.7 cm (5 in.) and a maximum capacity of 100 mm of rainfall. A mobile application developed by Mobile Water Management (MWM) (Mobile Water Management, 2020), the Netherlands, was used to record rainfall data for each rain gauge on a daily basis. The application has a an easily accessible and userfriendly interface where participants simply fill in the observed rainfall amount, take a photo of the rain gauge, and upload the information to the application. The photo and the rainfall data, along with the measurement location and time, are automatically stored in the database. Photos are used for visual validation of the recorded rainfall depth in order to eliminate errors.
In this study, the participants recruited were government officers, teachers, and local residents living close to the stations; these participants were trained to take measurements at around 07:00 LT (local time) daily according to the TMD standards. The quality of the collected data was assured by doublechecking the observations against the abovementioned highresolution photos as well as by the strict requirements with respect to measurement times, which ensured consistency with the TMD daily rainfall recording standard. Validation of the coneshaped citizen gauges was conducted based on a citizen gauge colocated with an automatic TMD gauge in the Tubma Basin, from August to October 2019. The citizen gauge installed at the same location as a TMD gauge (location R.3) showed good similarity with an RMSE and bias of 5.5 mm and 1.04, respectively (Fig. A2).
Quality control consisted of screening all citizen rain gauge data for errors and inconsistencies using doublemass curves. If citizen rain gauges reported >100 mm d^{−1} rainfall (maximum capacity of the citizen rain gauge), these data were excluded from the analysis. If days with norainfall data were found from all citizen rain gauges, the bias correction of that day was discarded from the dataset. Considering these data selection criteria for rainfall data recorded from August to October 2019, more than 80 % of the whole period for the bias adjustment process was then used for further analysis.
The methodology for radar rainfall bias correction using tippingbucket and citizen gauges consists of the following steps. Daily citizen rain gauge data were downscaled to an hourly timescale to be used as input for bias correction. The downscaling methods used in this paper are discussed in Sect. 3.1. Next, an hourly radar bias correction model was developed that combined rain gauge and downscaled citizen rain gauge data using a KF approach, as presented in Sect. 3.2.
3.1 Downscaling daily to hourly rainfall
To downscale the daily citizen rain gauge data to an hourly timescale, information on the hourly storm distribution pattern is needed. Methodologies to obtain the temporal rainfall distribution patterns are outlined in Table 1.
3.2 Hourly radar bias model
This section describes how radar bias is modeled (Sect. 3.2.1), how observations are assimilated in the model to correct the bias (Sect. 3.2.2), and how model parameters are estimated (Sect. 3.2.3). Our approach extends a previously used Kalman filter radar bias model by including two different types of rainfall observations (data from traditional and from citizen rain gauges) and using a likelihoodbased method for parameter estimation.
3.2.1 Kalman filter with two observations: model assumptions
Mean field bias (MFB) adjustment is a common technique for bias correction in radar rainfall relative to ground stations. It can be computed as the ratio of the mean hourly radar rainfall estimate to the rain gauge measurement (Anagnostou and Krajewski, 1999; Yoo and Yoon, 2010; Hanchoowong et al., 2013; Shi et al., 2018). However, direct application of the MFB does not account for uncertainty in the bias associated with each radargauge measurement. Alternatively, a KF has previously been used to estimate the spatially uniform MFB in real time in several studies, including Anhert et al. (1986), Smith and Krajewski (1991), Anagnostou et al. (1998), Seo et al. (1999), Chumchean et al. (2006), Kim and Yoo (2014), and Shi et al. (2018). The KF has the benefit of accounting for uncertainties in the observations by weighting the contribution of measurements by their respective error variances (Kalman, 1960). This is also the approach adopted here, but we extend it to include two types of observations with different error characteristics, i.e., hourly rain gauge data from TMD (y_{t}) and hourly downscaled citizen rain gauge data (z_{t}).
First, we define the logarithmic mean field radar rainfall bias β_{t} at hour t as follows (Smith and Krajewski, 1991; Anagostou et al., 1998):
where G_{i,t} and R_{i,t} are as defined in Eq. (1). Following previous studies (e.g., Smith and Krajewski, 1991; Chumchean et al., 2006), we model the logarithmic mean field radar rainfall bias as a firstorder autoregressive (AR1) process with stationary variance. Radar bias at time t is expressed in terms of the bias at previous time (β_{t−1}) and the process noise (W_{t}):
Here, r_{1} is the lagone correlation coefficient of the timevarying bias β_{t}, W_{t} is a zeromean random error with variance ${\mathit{\sigma}}_{\mathrm{W}}^{\mathrm{2}}$, and ${\mathit{\sigma}}_{\mathit{\beta}}^{\mathrm{2}}$ is stationary variance of the process. We consider two types of observations, y_{t} and z_{t}, of the unknown bias at time t, derived from the TMD and citizen rain gauges. Each observation is modeled as a random sample from a normal distribution conditioned on the underlying unknown bias with distinct measurement error variances ${\mathit{\sigma}}_{{\mathrm{M}}_{y,t}}^{\mathrm{2}}$ and ${\mathit{\sigma}}_{{\mathrm{M}}_{z,t}}^{\mathrm{2}}$:
Here, ${\mathit{\sigma}}_{{\mathrm{M}}_{y,t}}^{\mathrm{2}}$ and ${\mathit{\sigma}}_{{\mathrm{M}}_{z,t}}^{\mathrm{2}}$ are timevarying measurement error variances for the TMD and citizen rain gauges, respectively.
A factor graph representation of the radar bias and observation models is illustrated in Fig. 3a, with circles denoting random variables, and black squares denoting “factors” or relations between variables in the model.
The model contains two parameters, i.e., r_{1} and ${\mathit{\sigma}}_{\mathit{\beta}}^{\mathrm{2}}$, that need to be specified, together with values for the measurement error variances, ${\mathit{\sigma}}_{{\mathrm{M}}_{y,t}}^{\mathrm{2}}$ and ${\mathit{\sigma}}_{{\mathrm{M}}_{z,t}}^{\mathrm{2}}$. Parameter estimation will be discussed in Sect. 3.2.3. Section 3.2.2 first describes how observations are assimilated assuming that the values of the parameters and the measurement error variances have been specified.
3.2.2 Kalman filter with two observations: data assimilation
Having defined the model, we describe how observations are assimilated into the model, resulting in adjustments of the estimated bias. While the regular Kalman filter has two steps (prediction and update), assimilating two observations at each time step involves three steps, i.e., a prediction followed by two updates. Figure 3b shows a visual depiction of the prediction and update steps as probabilistic information propagating along the edges of the factor graph. These three steps are as follows:

Time update step (prediction). The first step for each hour t computes the respective prior mean ${\widehat{\mathit{\beta}}}_{t}^{}$ and variance ${P}_{t}^{}$ of β_{t} as
$$\begin{array}{}\text{(7)}& {\displaystyle}& {\displaystyle}{\widehat{\mathit{\beta}}}_{t}^{}={r}_{\mathrm{1}}{\widehat{\mathit{\beta}}}_{t\mathrm{1}},\text{(8)}& {\displaystyle}& {\displaystyle}{P}_{t}^{}={r}_{\mathrm{1}}^{\mathrm{2}}{P}_{t\mathrm{1}}+{\mathit{\sigma}}_{\mathrm{W}}^{\mathrm{2}}.\end{array}$$Here, P_{t−1} is the posterior variance of β_{t−1}. For the first time step 0 (t=0), we assume β_{0}=0 (climatological logarithmic mean field bias) and ${P}_{\mathrm{0}}=(\mathrm{1}{r}_{\mathrm{1}}^{\mathrm{2}}){\mathit{\sigma}}_{\mathit{\beta}}^{\mathrm{2}}$ (represents stationary process variance) (Smith and Krajewski, 1991; Chumchean et al., 2006).

The first measurement update step (first correction). The first update merges the bias prediction from step 1 with noisy observation y_{t} (with measurement error variance ${\mathit{\sigma}}_{{\mathrm{M}}_{y,t}}^{\mathrm{2}}$), resulting in posterior mean ${\widehat{\mathit{\beta}}}_{y,t}$ and variance P_{y,t} of β_{t} given by the Kalman update equations (Eqs. 9–11).
$$\begin{array}{}\text{(9)}& {\displaystyle}& {\displaystyle}{K}_{y,t}={P}_{t}^{}{\left({P}_{t}^{}+{\mathit{\sigma}}_{{\mathrm{M}}_{y,t}}^{\mathrm{2}}\right)}^{\mathrm{1}}\text{(10)}& {\displaystyle}& {\displaystyle}{\widehat{\mathit{\beta}}}_{y,t}={\widehat{\mathit{\beta}}}_{t}^{}+{K}_{y,t}\left({y}_{t}{\widehat{\mathit{\beta}}}_{t}^{}\right)\text{(11)}& {\displaystyle}& {\displaystyle}{P}_{y,t}=\left(\mathrm{1}{K}_{y,t}\right){P}_{t}^{}\end{array}$$Here, K_{y,t} is the Kalman gain for assimilating observation y_{t}. If there is no observation available at time t, K_{y,t} is zero (mathematically, a missing observation is equivalent to an observation with infinite variance ${\mathit{\sigma}}_{{\mathrm{M}}_{y,t}}^{\mathrm{2}}$ in Eq. 9), and the previous equations reduce to not performing any update, i.e., the posterior mean and variance are equal to the prior mean and variance.

The second measurement update step (second correction). The second update is done using the posterior values from the first correction (${\widehat{\mathit{\beta}}}_{y,t}$ and P_{y,t}) as the prior values. The resulting respective Kalman gain and posterior mean and variance are given by
$$\begin{array}{}\text{(12)}& {\displaystyle}& {\displaystyle}{K}_{z,t}={P}_{y,t}{\left({P}_{y,t}+{\mathit{\sigma}}_{{\mathrm{M}}_{z,t}}^{\mathrm{2}}\right)}^{\mathrm{1}},\text{(13)}& {\displaystyle}& {\displaystyle}{\widehat{\mathit{\beta}}}_{z,t}={\widehat{\mathit{\beta}}}_{y,t}+{K}_{z,t}\left({z}_{t}{\widehat{\mathit{\beta}}}_{y,t}\right),\text{(14)}& {\displaystyle}& {\displaystyle}{P}_{z,t}=\left(\mathrm{1}{K}_{z,t}\right){P}_{y,t}.\end{array}$$If there is no observation available at time t, Kalman gain K_{z,t} is zero, and these equations result in no update being applied, i.e., the posterior mean and variance are the same as after the first update.
As the logarithmic mean field radar rainfall bias β_{t} is assumed to be normally distributed, the mean field bias (B_{t}) is lognormally distributed with posterior mean at time t obtained from the posterior mean and variance of β_{t} (Smith and Krajewski, 1991):
$$\begin{array}{}\text{(15)}& {\widehat{B}}_{t}={\mathrm{10}}^{\left({\widehat{\mathit{\beta}}}_{z,t}+\mathrm{0.5}{P}_{z,t}\right)}.\end{array}$$The overall procedure of sequentially assimilating and downscaling the data is referred to as CKF (Kalman filter combined with the citizen rain gauge data), and it is summarized by the flowchart in Fig. 4. Operationally, it is implemented as follows. First, as the citizen rain gauge data were received in the last hour of day i, before receiving the downscaled hourly citizen rain gauge data of day i, the ordinary KF and observed hourly data of TMD were used to predict and correct the hourly bias adjustment factor. Second, if the citizen rain gauge data were available at the end of the day i, they were downscaled to an hourly timescale, as explained in Sect. 3.1. Third, the TMD and downscaled citizen rain gauge data were used together to conduct two steps of measurement update in the CKF process for all hourly time steps of day i. Finally, bias adjustment factors were applied at every hourly time step to obtain the final product of hourly radar rainfall estimation of day i. The bias factor for the last hour of the day i was used afterward as the initial value for calculating the ordinary KF of day i+1.
3.2.3 Kalman filter with two observations: parameter estimation
The model requires specification of the two parameters, r_{1} and $\mathit{\sigma}{}_{\mathit{\beta}}^{\mathrm{2}}$, and the measurement error variances, ${\mathit{\sigma}}_{{\mathrm{M}}_{y,t}}^{\mathrm{2}}$ and ${\mathit{\sigma}}_{{\mathrm{M}}_{z,t}}^{\mathrm{2}}$. We estimate the latter using formulas for the variance of the mean bias across individual TMD and citizen rain gauges:
Here, ${\mathit{\sigma}}_{{y}_{t}}^{\mathrm{2}}$ and ${\mathit{\sigma}}_{{z}_{t}}^{\mathrm{2}}$ are variances that quantify spatial variability of radar bias at time t at TMD and citizen rain gauge locations, respectively, and n_{y,t} and n_{z,t} are the corresponding number of observations at hour tfrom the two networks.
Once measurement error variances are specified, values for parameters r_{1} and ${\mathit{\sigma}}_{\mathit{\beta}}^{\mathrm{2}}$ are obtained by maximizing the marginal likelihood function (Bock and Aitkin, 1981; Harvey, 1990; Proietti and Luati, 2013; Pulido et al., 2018). As mentioned earlier, we have two sources of observed log mean field bias at hour t, from TMD (y_{t}) and citizen rain gauges (z_{t}). When only TMD data are available, the marginal likelihood is computed as follows:
where D is the data vector that contains all observed values, β_{t} is the true hidden state of log mean field bias at hour t, and T is the number of hourly time steps. When both TMD and citizen rain gauge data are available, the marginal likelihood is obtained as follows:
where $\mathit{\tau}=\frac{\mathrm{1}}{{\mathit{\sigma}}_{{\mathrm{M}}_{y,t}}^{\mathrm{2}}}+\frac{\mathrm{1}}{{\mathit{\sigma}}_{{\mathrm{M}}_{z,t}}^{\mathrm{2}}}$ and prior mean ${\widehat{\mathit{\beta}}}_{t}^{}$ and variance ${P}_{t}^{}$ are from the prediction step in Sect. 3.2.2.
To obtain optimal values of the two parameters r_{1} and ${\mathit{\sigma}}_{\mathit{\beta}}^{\mathrm{2}}$ that maximize the (logarithm of the) marginal likelihood, the Nelder–Mead simplex was used, which is an algorithm for searching a local optimum of a function (Lagarias et al., 1998; Luersen and Le Riche, 2004; Gao and Han, 2012).
3.3 Verification of the proposed bias correction approaches
To investigate which bias adjustment technique (MFB, KF, and CKF) gives the most suitable radar rainfall estimates for the Tubma Basin, the adjusted radar rainfall estimates were validated against measured rainfall data. There was only one automatic TMD rain gauge available in the basin, which was insufficient for validation purposes. Consequently, to test the performance of the hourly rainfall bias correction, data from 13 TMD stations located within a 90 km radius of the center of the Tubma Basin were used in addition to the one TMD station in the basin. Furthermore, daily timescale validation was conducted, using the daily rainfall data from 16 citizen rain gauges located in the Tubma Basin. A leaveoneout crossvalidation (LOOCV) algorithm was implemented to avoid bias in selecting the validation rain gauges. For each round of crossvalidation, one rain gauge was left out for validation, and the remaining rain gauges were used as the calibration rain gauges to calculate the bias adjustment factor using the three different techniques. This was repeated for all combinations, and the error of radar rainfall estimates after correction with the estimated bias factor at each radar pixel corresponding to the heldout gauge was then computed for all trials. In this study, the rootmeansquare error (RMSE) and mean bias error (MBE) were applied as statistical measures to evaluate the effectiveness of the different bias correction methods at each validation rain gauge. The number of possible combinations is equal to the total number of validated gauges (N_{G}). Data for the period from August to October 2019 were used in the evaluation. Four scenarios combining the three bias adjustment techniques were evaluated, as summarized in Table 2. These four scenarios are as follows:

KFTMDH – 1 TMD rain gauge from the 14 total gauges was left out for validation, and the remaining 13 gauges were separated to calculate the bias adjustment factors using the MFB and KF. This process was iterated 10 times until all 14 TMD rain gauges were left out for crossvalidation. Aggregated hourly rainfall values between the adjusted radar and the gauge rainfall data were compared to obtain the RMSE and MBE.

FTMDD – to identify whether the MFB or KF is more accurate for daily rainfall simulation in the Tubma Basin if there are only 14 TMD rain gauges available, 14 TMD and 16 citizen rain gauges were used for the analysis. All TMD gauges were used for assessing the MFB and KF, and estimated bias factors were applied at the daily timescale. An assessment of the RMSE and MBE of daily rainfall at all 16 citizen rain stations was used for validation.

CKFD – to evaluate the added value of using citizen rain gauges in the basin for bias correction, 15 citizen rain gauges (leave one citizen rain gauge out for validation) were used in addition to the TMD gauges following the CKF procedure explained in Sect. 3.2.2. Estimation of the daily RMSE and MBE was carried out at the heldout citizen rain gauge.

CKFH^{*} – to test whether the CKF with the most suitable storm pattern could benefit radar rainfall estimates in the area further away from the Tubma Basin, 14 TMD gauges were used to generate four cases of hourly rainfall distribution patterns, as described in Table 1, for downscaling the selected 16 daily citizen rain gauge data to an hourly timescale. The synthesized hourly citizen rain gauge data were later used as a second observation for the second correction of the CKF. A total of 13 TMD gauges (leave one TMD out) were used to produce estimates with the MFB and KF, and all 16 citizen rain gauges were merged for CKF computation.
All bias adjustment techniques evaluated the effectiveness at the heldout gauge for all possible combinations of the LOOCV procedure.
4.1 Simulation of the bias adjustment factor
4.1.1 Parameter estimation for the KF and CKF
Five scenarios were investigated for radar bias correction using the KF, based on TMD and citizen rain gauge observations, including four scenarios comparing different hourly downscaling approaches for the citizen rain gauge data (Table 1). Parameter estimates of the KF are shown in Table 3. These results indicate that the parameter r_{1}, the lagone correlation coefficient of the logarithmic mean field bias, ranges from 0.15 to 0.53, depending on the hourly downscaling approach, whereas ${\mathit{\sigma}}_{\mathit{\beta}}^{\mathrm{2}}$, representing the stationary variance of the logarithmic mean field bias, remains relatively invariant (ranging from 0.24 to 0.28) over the same period of simulation.
4.1.2 Bias adjustment factor comparison
To test the performance of the bias adjustment techniques among KFTMD, CKFR_{P}, CKFR_{MP}, CKFG_{MP}, and CKFG_{Tubma}, all approaches were used to assess the mean field bias for each hour using the data period from August to October 2019. The results were compared to the MFB calculated using the 14 TMD rain gauges (MFBTMD) in the Tubma Basin and within a 90 km radius of the basin. The results summarized in Fig. 5 show the following:

The daily observed bias is somewhat higher and shows larger variability for the citizen gauges compared with the TMD gauges. The hourly observed bias values based on downscaled citizen gauge data are in the same range as hourly bias values based on TMD gauges, with somewhat higher median values and spread (25–75 percentile range) for the R_{P} and G_{Tubma} downscaling scenarios.

The hourly observation error variance is smallest for the CKFR_{P} downscaling approach and somewhat larger for the other CKF approaches compared with the observation error variance for the TMD gauges.

Estimated hourly bias values based on KFTMD show a slightly higher mean and a smaller variability range compared with observations. The bias produced by the KFTMD is close to the MFBTMD if the observation error variance is small. In the case that no measured data are available for the bias update, the computed bias factor (B_{t}) progressively converges to 1.3, to meet the climatological logarithmic mean field bias.

Estimated bias values based on the CKF approaches are able to reproduce bias variability as observed by TMD gauges, with median values deviating by 0.2–0.4 and the value range being slightly larger for CKFR_{P} and smaller for CKFR_{MP}.

CKF gives different bias values according to the storm distribution pattern (see Appendix B) and the availability of the daily citizen rain gauge data used in combination with the KF. In the case that no citizen rain gauge data are available for updating, the bias generated by the CKF for every combination is close to the ordinary KF with small differences depending on their respective r_{1} and ${\mathit{\sigma}}_{\mathit{\beta}}^{\mathrm{2}}$ parameters.
4.2 Effectiveness evaluation of bias correction approaches
4.2.1 Hourly rainfall validation for the larger region (90 km radius) surrounding Tubma Basin using the MFB and KF approaches (case 1)
Figure 6a shows crossvalidation results based on the RMSE and MBE between TMD rain gauges and adjusted radar rainfall using the MFB and KF for hourly bias adjustment. Bias adjustment reduces the RMSE and, especially, the MBE, with KFTMD performing somewhat better than the MFBTMD (especially in terms of the RMSE). This confirms the ability of the KF approach that considers the error variance of observed hourly data as the weight for correcting the predicted mean bias instead of using only the calculated mean field bias (Smith and Krajewski, 1991; Chumchean et al., 2006).
4.2.2 Daily rainfall validation in the Tubma Basin using the MFB and KF approaches with citizen gauges for validation (case 2)
Results associated with validating the bias correction performance within the Tubma Basin are presented in Fig. 6b. This shows bias correction performance within the Tubma Basin for the MFB and KFbased daily bias adjustment. The two approaches show similar performance at the daily scale and improve the RMSE by 20 %–30 % and the MBE by 50 %–60 % (for median and upper 75percentile, respectively). The added value of a KFbased approach is limited for this case, as 14 TMD rain gauges in the region were used to compute observation variance which cannot represent the mean field bias behavior in the Tubma Basin.
4.2.3 Daily rainfall validation in the Tubma Basin using the CKF approach (case 3)
Figure 6c shows crossvalidation results at the daily scale for the Tubma Basin, comparing bias correction approaches using TMD only and TMD combined with citizen gauges. Following the CKF steps, citizen rain gauge data are downscaled to an hourly timescale using four different approaches, resulting in variation in hourly observed bias and error variances, as shown in Fig. 5b and c, respectively. Crossvalidation results after accumulation to the daily scale show that CKFR_{P} outperforms the other approaches (CKFR_{MP}, CKFG_{Tubma}, MFBTMD, KFTMD, and CKFG_{MP}) in terms of both the RMSE and MBE. The performance of the CKF techniques for radar rainfall simulation in the Tubma Basin relates to the reliability of the downscaled hourly observations. This is reflected in the variation in the estimated observation error variances for CKFR_{P}, as shown in Fig. 5b and c. The better performance of CKFR_{P} is explained by the smallest range in observation error variance, which is indicative of better consistency observation bias. In comparison with the “Nobias” scenario, CKFR_{P} can improve the RMSE by 32 %–25 % and can improve the MBE by 90 %–80 % for the median and upper 75 percentile, respectively. In contrast, CKFG_{MP} exhibits the worst performance compared with the other CKF approaches, with a 13 %–16 % improvement in the RMSE and an 57 %–56 % improvement in the MBE for the median and upper 75 percentile, respectively. This apparent decrease in the efficiency of the CKF can be confirmed by the highest median value of the estimated observation error variances of CKFG_{MP} (see Fig. 5c), which is 33 % higher than that of CKFR_{P}.
4.2.4 Hourly rainfall validation using the MFB, KF, and CKF approaches (case 4)
Results for this section are presented in Figs. 7 and 8. The performance for radar rainfall estimates using CKFR_{P} with different distances from the Tubma Basin is reported in Appendix C. Crossvalidation results at an hourly timescale show a strong improvement achieved by bias adjustment using citizen gauges, in particular close to the Tubma Basin where the citizen gauges are located. Figure 7b and c show validation results based on TMD gauges for gauges close to (0–40 km radius) and further away from (40–90 km radius) the center of Tubma Basin (see Fig. 7a); both ranges cover a similar number of TMD gauges. Figure 7b and c show that CKFR_{P} bias adjustment significantly improves radar rainfall estimates at an hourly timescale, compared with bias adjustment approaches based on TMD gauges only in the 0–40 km range closest to the Tubma Basin. While there is a modest improvement in the mean RMSE (see the black line connecting the mean values of the box plots from the MFBTMD to CKFR_{P}), the upper 75 percentile RMSE is reduced from about 6 to 3.5 mm h^{−1}. The mean MBE is changed from 0.1 to −0.15 mm h^{−1} (see the red dotted line connecting the mean values from the MFBTMD to CKFR_{P}). For the 40–90 km range, CKFR_{P} performs similarly to the MFBTMD and KFTMD.
Figure 8 illustrates that hourly rainfall distribution patterns of TMD rain gauges in the 0–40 km range, which are mainly influenced by the southwest monsoon, appear to be more similar to the mean citizen rain gauge data than the range beyond 40 km. Consequently, the application of CKFR_{P} based on combining the citizen rain gauge network with the TMD rain gauge network (with similarity of rainfall characteristic) is key for improving radar rainfall estimates.
The results in Fig. 7 also show that the MBE values in the 0–40 km range are explicitly lower than those in the 40–90 km range. Apparently, at shorter range, positive and negative errors represented in the MBE cancel out more frequently than they do for the gauges at larger distance. In other words, gauges more or less randomly over or underestimate rainfall values: we can see similar rainfall distribution patterns among all gauges, although there is high variation in the rainfall amount during the storm period in Fig. 8b. Conversely, in the 40–90 km range, bias correction at gauge locations consistently leads to over or underestimation of rainfall. This can be explained by gauges at larger distances from the center of the basin being affected by different rainfall generation patterns, associated with their location closer to the coast or mountains (see Fig. 7a). The influence of the southwest monsoon strongly affects all gauges located in the coastal region on the windward side of a mountain, whereas rain gauge locations on the leeward side experience less rainfall. Figure 8c shows that TMD gauges located on the leeward side (e.g., 4590009 and 4590011) obviously appear as a steady gradient of the mass curve, reflecting light rainfall accumulation, whereas gauges on the windward side (e.g., 4590002 and 4590003) show mass curves with a sharper gradient. Further details on the sensitivity bias adjustment techniques' accuracy with respect to rainfall characteristics are provided in Appendix D.
In this study, we introduced a modified KF approach in radar bias correction in the Tubma Basin, eastern Thailand. The twostep KF integrates daily data from a dense citizen rain gauge network with hourly data from a much sparser network of conventional higherquality tippingbucket rain gauges. Daily citizen rain gauge observations were downscaled to an hourly timescale using four different approaches. The question that we aimed to answer was as follows: to what extent do the downscaled citizen rainfall observations improve the accuracy of hourly radar rainfall estimates? Results showed that citizen rain gauges significantly improve the performance of radar rainfall bias adjustment, up to a range of about 40 km from the center of the Tubma Basin (197 km^{2}) where the citizen rain gauge network is located. While a modest improvement in the mean RMSE was obtained, the upper 75 percentile RMSE was reduced from 6 to 3.5 mm h^{−1}. The MBE was changed from 0.1 to −0.15 mm h^{−1} across the validation period (August–October 2019). In the Tubma Basin, beyond the 40 km range, no significant improvement from the inclusion of the citizen gauges was found. The rainfall distribution pattern is key for downscaling the daily measured citizen rain gauge observations to an hourly temporal resolution. In the study region, we found that downscaling based on the rainfall patterns derived from hourly radar rainfall at overlying radar pixels corresponding to the citizen gauge location was the most suitable technique, resulting in the smallest variation in observation error variances of the mean field bias. In the case of a sparse rain gauge network, the mean field bias and the Kalman filter approach both show improvement, and the degree of improvement was similar between the two approaches. In other words, in a sparse gauge network, the added value of error information represented in the Kalman filter is limited. Note that citizen rain gauge data are available only at the end of the day; consequently, the modified twostep Kalman filter, as used in this study, has limitations with respect to realtime applications. However, the method proposed here has clear potential when creating highquality historical radarrainfall time series for climatological studies and in postevent analysis. Moreover, nearrealtime assessment could be achievable if the citizen rain gauge data were collected at a subdaily timescale.
An example of installing a coneshaped transparent plastic citizen rain gauge (in standard use in South Africa) at location R.22 (Map Tong school, Rayong Province) is illustrated in Fig. A1. Validation of the coneshaped citizen gauges was done at location R.3 by comparing the measured daily citizen gauge data with a colocated automatic TMD gauge, as presented the result in Fig. A2. The rootmeansquare error (RMSE) and bias were applied as statistical measures for the validation as shown in Eqs. (A1) and (A2), respectively:
Here, G_{CR,t} is daily citizen gauge rainfall at location R.3 for day t, G_{TMD,t} is daily TMD gauge rainfall at the same location of R.3 for day t, and T is the total period used in the calculation.
Four hourly rainfall distribution patterns were obtained, as outlined in Table 1. Figure B1 illustrates the cumulative fraction of daily rainfall at the hourly scale during the simulation period from August to October 2019. It can be seen that most rainfall was concentrated in the afternoon hours, with very little rainfall falling at night. R_{P} and R_{MP} showed relatively more rainfall concentrated in the afternoon, while R_{P} and G_{MP} showed larger variability in the downscaled hourly data with substantial outliers in the box plots, which were associated with the variability in the locations underlying the rainfall distributions (multiple radar pixels within the Tubma Basin for R_{P} versus multiple TMD gauges surrounding the Tubma Basin for G_{MP}). G_{MP} showed the flattest distribution with the longest rainy period of around 11 h compared with the others which had a period of heavy rainfall of around 4–5h d^{−1}. This is explained by the larger spatial variability in the gauges covered by G_{MP}.
We chose the 40 km separation boundary to achieve an equal number of gauges in the “near” and “far” groups. Following the referee's comment, we investigated the RMSE between gauge rainfall and radar rainfall without bias adjustment (RMSE_{Nobias}) and with the CKFR_{P} (RMSE_{CKFRP}) for individual stations (located at distances of 5–80 km from the Tubma Basin). We computed the percentage improvement in radar rainfall estimates using CKFR_{P} compared with Nobias at each rain gauge, indicating that the relative errors change with distance from the Tubma Basin (Fig. C1). Figure C1 shows that the improvement percentage of using CKFR_{P} tends to decrease with increasing distance from the Tubma Basin, where the citizen gauges are located. The percentage reduction gradually decreases beyond a distance of about 40 km.
Note that gauge 4780010, situated nearest to the basin, is expected to provide the best improvement; however, the Sattahip radar temporarily stopped measuring for 3 h (from 16:00 to 19:00 LT on 24 September 2019) due to the presence of a heavy storm's center that was only observed at this gauge. This leads to significant degradation of the radar rainfall performance. Furthermore, the lower percentage improvement for station 4780005 is associated with a localized heavy rainfall event that was only recorded at this location and negatively affected its performance as a representative station for bias correction.
Figure 7b shows that the upper 75 percentile RMSE of the shorter range is remarkably high while using only TMD gauges for the bias adjustment. These errors occurred for 3 h at three different gauge locations when heavy rainfall data were only measured at the validated gauge location while there was relatively uniform light rainfall at all available surrounding TMD gauges used for the bias adjustment calculation. Consequently, the calculated bias factors from the available gauges could not represent the heavy rainfall at the tested location, leading to the significant RMSE. An analysis of hourly rainfall hyetographs obtained from the TMD rain gauge network compared with the validated rain gauge occurring on 3 different days is illustrated in Fig. D1. It shows a considerable RMSE for 3 h on 3 d – 15 September 2019, 12:00 LT; 21 September 2019, 15:00 LT; and 22 September 2019, 14:00 LT – associated with the validated gauges, 4780001, 4780005, and 4780003, respectively. However, these RMSE values decrease considerably if the CKFR_{P} was only implemented in the shorter range.
Indeed, Fig. 8a and b show that one of the citizen gauges (R.18) collected a cumulative monthly rainfall amount that was higher than the range of the other citizen and TMD rain gauges. This is associated with a storm event that occurred in September 2019 (the storm center was over R.18, while the surroundings received appreciably less rainfall). Figure D2 shows the radar reflectivity field at 13.00 LT on 22 September 2019, during the peak of the storm, confirming the heavy rainfall affecting gauge location R.18. This shows that the citizen gauge network is able to capture local storm features thanks to the high density of the network. The multiple reporting gaps visible in Fig. 8a are caused by time errors in the observations submitted by local residents which were removed from the analysis (as explained in Sect. 2.3.2).
All raw data can be provided by the corresponding authors upon request.
The study was designed by PPM, TAB, and MCTV. MM, PPM, and TAB were responsible for data collection and validation. MM, PPM, and MCTV developed and conducted the radar bias adjustment procedure. GS, MM, and PPM developed and implemented the twostep Kalman filter. The paper was written by PPM and edited by MCTV, TAB, and GS.
At least one of the (co)authors is a member of the editorial board of Hydrology and Earth System Sciences. The peerreview process was guided by an independent editor, and the authors also have no other competing interests to declare.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This work was undertaken within the framework of a research exchange between Kasetsart University (KU), Thailand, and Delft University of Technology (TU Delft), the Netherlands. Punpim Puttaraksa Mapiam, Monton Methaprayun, and Thom Bogaard gratefully acknowledge the Faculty of Engineering at KU for financially supporting this research, including Punpim Puttaraksa Mapiam's 3week guest visit to TU Delft, Monton Methaprayun's 3month guest visit to TU Delft, and Thom Adrianus Bogaard's 2month guest visit to KU. The authors wish to thank the Department of Royal Rainmaking and Agricultural Aviation (DRRAA) and the Thai Meteorological Department (TMD) for providing the radar reflectivity and rain gauge data used in this study. We also appreciate Mobile Water Management (MWM), the Netherlands, for providing the mobile application used to measure daily citizen rain gauge rainfall in the study area.
This paper was edited by Nadav Peleg and reviewed by five anonymous referees.
Amitai, E.: Systematic Variation of Observed Radar Reflectivity–Rainfall Rate Relations in the Tropics, J. Appl. Meteorol., 39, 2198–2208, https://doi.org/10.1175/15200450(2001)040<2198:SVOORR>2.0.CO;2, 2000.
Anagnostou, E. N. and Krajewski, W. F.: RealTime Radar Rainfall Estimation. Part I: Algorithm Formulation, J. Atmos. Ocean. Tech., 16, 189–197, https://doi.org/10.1175/15200426(1999)016<0189:RTRREP>2.0.CO;2, 1999.
Anagnostou, E. N., Krajewski, W. F., Seo, D.J., and Johnson, E. R.: MeanField Rainfall Bias Studies for WSR88D, J. Hydrol. Eng., 3, 149–159, https://doi.org/10.1061/(ASCE)10840699(1998)3:3(149), 1998.
Anhert, P., Krajewski, W. F., and Johnson, E. R.: Kalman Filter estimation of radarrainfall mean field bias, in: 23rd Radar Meteorology Conf., Amer. Meteor. Soc., 33–37, https://www.jstor.org/stable/26225332 (last access: 2 February 2022), 1986.
Arai, K., Liang, X., and Liu, Q.: Method for estimation of rain rate with Rayleigh and Mie scattering assumptions on the ZR relationship for different rainfall types, Adv. Space Res., 36, 813–817, https://doi.org/10.1016/j.asr.2005.04.102, 2005.
Austin, P. M.: Relation between Measured Radar Reflectivity and Surface Rainfall, Mon. Weather Rev., 115, 1053–1070, https://doi.org/10.1175/15200493(1987)115<1053:RBMRRA>2.0.CO;2, 1987.
Barton, Y., Sideris, I. V., Germann, U., and Martius, O.: A method for realtime temporal disaggregation of blended radar–rain gauge precipitation fields, Meteorol. Appl., 27, e1843, https://doi.org/10.1002/met.1843, 2020.
Battan, L. J.: Radar observation of the atmosphere, Rev. ed. Edn., University of Chicago press, Chicago, https://doi.org/10.1002/qj.49709942229, 1973.
Bedient, P. B., Holder, A., Beneavides, J. A., and Vieux, B.: Radarbased flood warning system applied to tropical storm Allison, J. Hydrol. Eng., 8, 308–318, https://doi.org/10.1061/(ASCE)10840699(2003)8:6(308), 2003.
Berne, A. and Krajewski, W. F.: Radar for hydrology: Unfulfilled promise or unrecognized potential?, Adv. Water Resour., 51, 357366, https://doi.org/10.1016/j.advwatres.2012.05.005, 2013.
Bock, R. D. and Aitkin, M.: Marginal maximum likelihood estimation of item parameters: Application of an EM algorithm, Psychometrika, 46, 443–459, https://doi.org/10.1007/BF02293801, 1981.
Chantraket, P., Detyothin, C., Pankaew, S., and Kirtsaeng, S.: An Operational Weather RadarBased Calibration of Z–R Relationship over Central Region of Thailand, Int. J. Eng. Issues, 2, 92–100, 2016.
Chumchean, S., Seed, A., and Sharma, A.: Correcting of realtime radar rainfall bias using a Kalman filtering approach, J. Hydrol., 317, 123–137, https://doi.org/10.1016/j.jhydrol.2005.05.013, 2006.
Collinge, V. K. and Kirby, C.: Weather Radar and Flood Forecasting, Wiley, Chichester, England, ISBN 0471912964, 2003.
Corral, C., Berenguer, M., SempereTorres, D., Poletti, L., Silvestro, F., and Rebora, N.: Comparison of two early warning systems for regional flash flood hazard forecasting, J. Hydrol., 572, 603–619, https://doi.org/10.1016/j.jhydrol.2019.03.026, 2019.
Creutin, J.D. and Borga, M.: Radar hydrology modifies the monitoring of flashflood hazard, Hydrol. Process., 17, 1453–1456, https://doi.org/10.1002/hyp.5122, 2003.
Davids, J. C., Devkota, N., Pandey, A., Prajapati, R., Ertis, B. A., Rutten, M. M., Lyon, S. W., Bogaard, T. A., and van de Giesen, N.: Soda Bottle ScienceCitizen Science Monsoon Precipitation Monitoring in Nepal, Front. Earth Sci., 7, 46, https://doi.org/10.3389/feart.2019.00046, 2019.
Debele, B., Srinivasan, R., and Yves Parlange, J.: Accuracy evaluation of weather data generation and disaggregation methods at finer timescales, Adv. Water Resour., 30, 1286–1300, https://doi.org/10.1016/j.advwatres.2006.11.009, 2007.
Dinku, T., Anagnostou, E. N., and Borga, M.: Improving RadarBased Estimation of Rainfall over Complex Terrain, J. Appl. Meteorol., 41, 1163–1178, https://doi.org/10.1175/15200450(2002)041<1163:IRBEOR>2.0.CO;2, 2002.
Doelling, I. G., Joss, J., and Riedl, J.: Systematic variations of Z–Rrelationships from drop size distributions measured in northern Germany during seven years, Atmos. Res., 47–48, 635–649, https://doi.org/10.1016/S01698095(98)00043X, 1998.
Doviak, R. J. and Zmic, D. S.: Doppler Radar and Weather Observations, 1st Edn., Acadamic Press Inc., San Diego, California, USA, 575 pp., ISBN 9780323149167, 1984.
Fabry, F., Bellon, A., Duncan, M. R., and Austin, G. L.: High resolution rainfall measurements by radar for very small basins: the sampling problem reexamined, J. Hydrol., 161, 415–428, https://doi.org/10.1016/00221694(94)901384, 1994.
Fang, X., Shao, A., Yue, X., and Liu, W.: Statistics of the Z–R Relationship for Strong Convective Weather over the YangtzeHuaihe River Basin and Its Application to Radar Reflectivity Data Assimilation for a Heavy Rain Event, J. Meteorol. Res., 32, 598–611, https://doi.org/0.1007/s1335101871631, 2018.
Fulton, R. A., Breidenbach, J. P., Seo, D.J., Miller, D. A., and O'Bannon, T.: The WSR88D Rainfall Algorithm, Weather Forecast., 13, 377–395, https://doi.org/10.1175/15200434(1998)013<0377:TWRA>2.0.CO;2, 1998.
Gao, F. and Han, L.: Implementing the NelderMead simplex algorithm with adaptive parameters, Comput. Optimiz. Appl., 51.1, 259–277, 2012.
Germann, U., Galli, G., Boscacci, M., and Bolliger, M.: Radar precipitation measurement in a mountainous region, Q. J. Roy. Meteorol. Soc., 132, 1669–1692, https://doi.org/10.1256/qj.05.190, 2006.
Hagen, M. and Yuter, S.: Relations between Radar Reflectivity, LiquidWater Content, and Rainfall Rate during the MAP SOP, Q. J. Roy. Meteorol. Soc., 129, 477–493, https://doi.org/10.1256/qj.02.23, 2003.
Hanchoowong, R., Weesakul, U., and Chumchean, S.: Bias correction of radar rainfall estimates based on a geostatistical technique, Science Asia, 38, 373–385, https://doi.org/10.2306/scienceasia15131874.2012.38.373, 2012.
Harvey, A. C.: Forecasting, structural time series models and the Kalman filter, Cambridge University Press, ISBN 9781107049994, https://doi.org/10.1017/CBO9781107049994 1990.
Hitschfeld, W. and Bordan, J.: Errors inherent in the radar measurement of rainfall at attenuating wavelengths, J. Meteorol., 11, 58–67, https://doi.org/10.1175/15200469(1954)011<0058:EIITRM>2.0.CO;2, 1954.
Jordan, P., Seed, A., and Austin, G.: Sampling errors in radar estimates of rainfall, J. Geophys. Res., 105, 2247–2250, https://doi.org/10.1029/1999JD900130, 2000.
Joss, J. and Waldvogel, A.: Raindrop size distribution and Doppler velocities, in: Proceedings 14th Radar Meteorology Conference, Boston, https://www.jstor.org/stable/26253228 (last access: 10 February 2020), 1970.
Kalman, R. E.: A New Approach to Linear Filtering and Prediction Problems, J. Basic Eng., 82, 35–45, https://doi.org/10.1115/1.3662552, 1960.
Kim, J. and Yoo, C.: Using Extended Kalman Filter for Realtime Decision of Parameters of Z–R Relationship, J. Korea Water Resour. Assoc., 47, 119–133, https://doi.org/10.3741/JKWRA.2014.47.2.119, 2014.
Klazura, G. E.: Differences between Some RadarRainfall Estimation Procedures in a High Rain Rate Gradient Storm, J. Appl. Meteorol. Clim., 20, 1376, https://doi.org/10.1175/15200450(1981)020<1376:Dbsrre>2.0.Co;2, 1981.
Koutsoyiannis, D.: Rainfall disaggregation methods: Theory and applications, Università di Roma “La Sapienza”, Rome, 1–23, https://doi.org/10.13140/RG.2.1.2840.8564, 2003.
Lagarias, J., Reeds, J., Wright, M., and Wright, P.: Convergence Properties of the Nelder–Mead Simplex Method in Low Dimensions, SIAM J. Optimiz., 9, 112–147, https://doi.org/10.1137/S1052623496303470, 1998.
Luersen, M. A. and Le Riche, R.: Globalized Nelder–Mead method engineering optimization, Comput. Struct., 82, 2251–2260, https://doi.org/10.1016/j.compstruc.2004.03.072, 2004.
Mapiam, P. P. and Chautsuk, S.: Improving Runoff Estimates by Increasing Catchment Subdivision Complexity and Resolution of Rainfall Data in the Upper Ping River Basin, Thailand, Chiang Mai Univers. J. Nat. Sci., 17, 127–143, https://doi.org/10.12982/CMUJNS.2018.0010, 2018.
Mapiam, P. P. and Sriwongsitanon, N.: Climatological Z–R relationship for radar rainfall estimation in the upper Ping river basin, Science Asia, 34, 215–222, https://doi.org/10.2306/scienceasia15131874.2008.34.215, 2008.
Mapiam, P. P., Sharma, A., Chumchean, S., and Sriwongsitanon, N.: Runoff estimation using radar and rain gage data, in: The 18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation, 13–17 July 2009, Cairns, Australia, 2009a.
Mapiam, P. P., Sriwongsitanon, N., Chumchean, S., and Sharma, A.: Effects of Rain Gauge Temporal Resolution on the Specification of a ZR Relationship, J. Atmos. Ocean. Tech., 26, 1302–1314, https://doi.org/10.1175/2009JTECHA1161.1, 2009b.
Mapiam, P. P., Sharma, A., and Sriwongsitanon, N.: Defining the Z–R Relationship Using Gauge Rainfall with Coarse Temporal Resolution: Implications for Flood Forecasting, J. Hydrol. Eng., 19, 04014003, https://doi.org/10.1061/(ASCE)HE.19435584.0000616, 2014.
Marshall, J. S. and Palmer, W. M. K.: The Distribution Of Raindrops With Size, J. Meteorol., 5, 165–166, https://doi.org/10.1175/15200469(1948)005<0165:TDORWS>2.0.CO;2, 1948.
Mobile Water Management: https://mobilewatermanagement.nl/, last access: 25 October 2020.
Paulat, M., Frei, C., Hagen, M., and Wernli, H.: A gridded dataset of hourly precipitation in Germany: Its construction, climatology and application, Meteorol. Z., 17, 719–732, https://doi.org/10.1127/09412948/2008/0332, 2008.
Proietti, T. and Luati, A.: Maximum likelihood estimation of time series models: the Kalman filter and beyond, in: Handbook of Research Methods and Applications in Empirical Macroeconomics, Edward Elgar Publishing, ISBN 9781782545071, 2013.
Pulido, M., Tandeo, P., Bocquet, M., Carrassi, A., and Lucini, M.: Stochastic parameterization identification using ensemble Kalman filtering combined with maximum likelihood methods, Tellus A, 70, 1–17, https://doi.org/10.1080/16000870.2018.1442099, 2018.
Rogers, R. R.: The effect of variable target reflectivity on weather radar measurements, Q. J. Roy. Meteorol. Soc., 97, 154–167, https://doi.org/10.1002/qj.49709741203, 1971.
Rosenfeld, D., Atlas, D., Wolff, D., and Amitai, E.: Beamwidth Effects on Z–R Relations and Areaintegrated Rainfall, J. Appl. Meteorol. Clim., 31, 454–464, https://doi.org/10.1175/15200450(1992)031<0454:BEOZRR>2.0.CO;2, 1992.
Rosenfeld, D., Wolff, D., and Atlas, D.: General Probabilitymatched Relations between Radar Reflectivity and Rain Rate, J. Appl. Meteorol. Clim., 32, 50–72, https://doi.org/10.1175/15200450(1993)032<0050:GPMRBR>2.0.CO;2, 1993.
See, L.: A Review of Citizen Science and Crowdsourcing in Applications of Pluvial Flooding, Front. Earth Sci., 7, 44, https://doi.org/10.3389/feart.2019.00044, 2019.
Seed, A., Siriwardena, L., Sun, X., Jordan, P., and Elliot, J.: On the calibration of Australian weather radars, Cooperative Research Centre for Catchment Hydrology, Australia, 49 pp., https://ewater.org.au/archive/crcch/archive/pubs/pdfs/technical200207.pdf (last access date: 2 February 2022), 2002.
Seibert, J., Strobl, B., Etter, S., Hummer, P., and van Meerveld, H. J.: Virtual Staff Gauges for CrowdBased Stream Level Observations, Front Earth Sci., 7, 70, https://doi.org/10.3389/feart.2019.00070, 2019.
Seo, D. J., Breidenbach, J. P., and Johnson, E. R.: Realtime estimation of mean field bias in radar rainfall data, J. Hydrol., 223, 131–147, https://doi.org/10.1016/S00221694(99)001067, 1999.
Shi, Z., Wei, F., and Chandrasekar, V.: Radarbased quantitative precipitation estimation for the identification of debris flow occurrence over earthquakeaffected regions in Sichuan, China, Nat. Hazards Earth Syst. Sci., 18, 765–780, https://doi.org/10.5194/nhess187652018, 2018.
Sideris, I. V., Gabella, M., Erdin, R., and Germann, U.: Realtime radar–raingauge merging using spatiotemporal cokriging with external drift in the alpine terrain of Switzerland, Q. J. Roy. Meteorol. Soc., 140, 1097–1111, https://doi.org/10.1002/qj.2188, 2014.
Smith, J., Baeck, M., Meierdiercks, K., Miller, A., and Krajewski, W.: Radar Rainfall Estimation for Flash Flood Forecasting in Small Urban Watersheds, Adv. Water Resour., 30, 2087–2097, https://doi.org/10.1016/j.advwatres.2006.09.007, 2007.
Smith, J., Hui, E., Steiner, M., Baeck, M., Krajewski, W., and Ntelekos, A.: Variability of rainfall rate and raindrop size distributions in heavy rain, Water Resour. Res., 45, W04430, https://doi.org/10.1029/2008WR006840, 2009.
Smith, J. A. and Krajewski, W. F.: Estimation of the Mean Field Bias of Radar Rainfall Estimates, J. Appl. Meteorol. Clim., 30, 397–412, https://doi.org/10.1175/15200450(1991)030<0397:Eotmfb>2.0.Co;2, 1991.
Smith, P. L., Myers, C. G., and Orville, H. D.: Radar Reflectivity Factor Calculations in Numerical Cloud Models Using Bulk Parameterization of Precipitation, J. Appl. Meteorol., 14, 1156–1165, 1975.
Srivastava, S., Vaddadi, S., and Sadistap, S.: Smartphonebased System for water quality analysis, Appl. Water Sci., 8, 130, https://doi.org/10.1007/s1320101807800, 2018.
Steiner, M. and Smith, J. A.: Reflectivity, Rain Rate, and Kinetic Energy Flux Relationships Based on Raindrop Spectra, J. Appl. Meteorol., 39, 1923–1940, https://doi.org/10.1175/15200450(2000)039<1923:RRRAKE>2.0.CO;2, 2000.
Steiner, M., Houze, R., and Yuter, S.: Climatological characterization of threedimensional storm structure from operational radar and rain gauge data, J. Appl. Meteorol., 34, 1978–2007, https://doi.org/10.1175/15200450(1995)034<1978:CCOTDS>2.0.CO;2, 1995.
Steiner, M., Smith, J., Kessinger, C., and Ferrier, B.: Evaluation of algorithm parameters for radar data quality control, https://collaborate.princeton.edu/en/publications/evaluationofalgorithmparametersforradardataqualitycontrol (last access: 10 February 2022), 1999.
Sun, X., Mein, R. G., Keenan, T. D., and Elliott, J. F.: Flood estimation using radar and raingauge data, J. Hydrol., 239, 4–18, https://doi.org/10.1016/S00221694(00)003504, 2000.
Tokay, A. and Short, D. A.: Evidence from Tropical Raindrop Spectra of the Origin of Rain from Stratiform versus Convective Clouds, J. Appl. Meteorol., 35, 355–371, 1996.
Uijlenhoet, R.: Raindrop size distributions and radar reflectivity–rain rate relationships for radar hydrology, Hydrol. Earth Syst. Sci., 5, 615–628, https://doi.org/10.5194/hess56152001, 2001.
Ulbrich, C. W.: Natural Variations in the Analytical Form of the Raindrop Size Distribution, J. Clim. Appl. Meteorol., 22, 1764–1775, 1983.
Vieux, B. E. and Bedient, P. B.: Assessing urban hydrologic prediction accuracy through event reconstruction, J. Hydrol., 299, 217–236, https://doi.org/10.1016/j.jhydrol.2004.08.005, 2004.
Vormoor, K. and Skaugen, T.: Temporal Disaggregation of Daily Temperature and Precipitation Grid Data for Norway, J. Hydrometeorol., 14, 989–999, https://doi.org/10.1175/JHMD120139.1, 2013.
Wang, L.P., OchoaRodríguez, S., Van Assel, J., Pina, R. D., Pessemier, M., Kroll, S., Willems, P., and Onof, C.: Enhancement of radar rainfall estimates for urban hydrology through optical flow temporal interpolation and Bayesian gaugebased adjustment, J. Hydrol., 531, 408–426, https://doi.org/10.1016/j.jhydrol.2015.05.049, 2015.
Wilson, J. W. and Brandes, E. A.: Radar Measurement of Rainfall – A Summary, B. Am. Meteorol. Soc., 60, 1048, https://doi.org/10.1175/15200477(1979)060<1048:Rmors>2.0.Co;2, 1979.
Wilson, J. W. J.: Integration of radar and raingage data for improved rainfall measurement, J. Appl. Meteorol. Clim., 9, 489–497, 1970.
Wüest, M., Frei, C., Altenhoff, A., Hagen, M., Litschi, M., and Schär, C.: A gridded hourly precipitation dataset for Switzerland using raingauge analysis and radarbased disaggregation, In. J. Climatol., 30, 1764–1775, https://doi.org/10.1002/joc.2025, 2009.
Yoo, C. and Yoon, J.H.: A proposal of quality evaluation methodology for radar data, J. Korean Soc. Civ. Eng., 30, 429–435, 2010.
 Abstract
 Introduction
 Study area and data
 Methods
 Results and discussion
 Conclusion
 Appendix A: Citizen rain observation and validation
 Appendix B: Hourly rainfall distribution patterns
 Appendix C: The performance of radar rainfall estimates using CKFR_{P} with different distances from the Tubma Basin
 Appendix D: Investigation of rainfall characteristics affecting the accuracy of bias adjustment techniques
 Code and data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Review statement
 References
 Abstract
 Introduction
 Study area and data
 Methods
 Results and discussion
 Conclusion
 Appendix A: Citizen rain observation and validation
 Appendix B: Hourly rainfall distribution patterns
 Appendix C: The performance of radar rainfall estimates using CKFR_{P} with different distances from the Tubma Basin
 Appendix D: Investigation of rainfall characteristics affecting the accuracy of bias adjustment techniques
 Code and data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Review statement
 References