Rain rate modeling of 1-min from various integration times in South Korea

Rain plays a major impairment factor for propagation of electromagnetic waves in atmosphere for systems operating at frequencies above 10 GHz. Several effects are noted such as depolarization, scintillation, interference due to scattering and extra attenuation which seems to increase with frequency. To mitigate its effect in satellite communication, knowledge of local rainfall statistics is necessary which act as milestone for design of radio link. Rain attenuation is best visualize by the 1-min rainfall rate statistic but the measurement of this rain rate distribution is rare on a worldwide basis and observation of rain rate are done with longer integration times typically 30 min or more. In this paper, efforts have been made to develop model that can convert rain rate complementary cumulative distribution function to shorter integration times. The average relative error margin of about 5, 14, 43, 71 and 115 % are noted for 5 to 1-, 10 to 1-, 20 to 1, 30 to 1- and 60 to 1-min respectively from ITU-R P.837-6 method which have been analyzed in further section of this article. The empirical natures of conversion methods as such Segal method, Burgueno’s method, Chebil and Rahman method and Logarithmic model are studied along with the proposed new model that seems to be applicable in derivation of 1-min rain rate of the South Korea rain rate statistics. International Telecommunication Union-Radio communication Sector (ITU-R) has developed a recommendation ITU-R P.837-6 that enables the user to estimate the local 1-min rainfall rate statistical distribution which is compared with calculated 1-min rain rate distribution from experimental 1-min rainfall accumulation. Unfortunately, ITU-R P.837-6 estimated 1-min values show greater error percentages. In order to get better approximation of local 1-min rain rate estimation, a novel method is proposed and it’s efficiency have been compared with rainfall rate statistics obtained from nine different locations in the South Korea.

The approach used in this study model the long term behavior of rain by analyzing the whole set of data without attempting to classify it according to event types or rain intensity. For the prediction of rain attenuation, accurate knowledge of 1-min integrated CCDF of rain rate is required as recommended by Series (2015). But difficulty to obtain such short interval rain rate have forced researcher worldwide to find prediction method that shall be able to predict local rain rate characteristics (Aris et al. 2013). Paper proposed a prediction model based on the rainfall data collected by Korea Meteorological Administration (KMA) which has developed and operated a digital system for collecting and storing rainfall data at 1-min interval since 2004. Performance of proposed model is compared with ITU-R P.837-6 (2012) and found to give better prediction performance as it makes reduction in error analysis. The material of the paper is covered as follows: "Background" section briefly introduces background of established relationships for derivation of 1-min rain rate. "Methodology and analyses of experimental data" section deals with experimental data collection procedure. "Results and discussion" section presents the result and discussion of statistical analyses. Finally, "Conclusion" section draws conclusion of highly reliable statistical results.

Background
Rainfall data of longer integration time relatively on hourly basis is readily available, 5or 10-min accumulation data are recorded by several weather services, but 1-min accumulation interval data is available from only special observations at a few locations. Under this scenario, derivation of model for statistical distribution of short time interval requires the limited number of available observation to be combined in order to provide statistically valid empirical distribution function (Crane 1996). As the earlier approach is presented by ITU-R with latest global model (2012), its usefulness is depicted by Emiliani et al. (2010), Capsoni and Lorenzo (2009), Capsoni et al. (2009) which highlight the benefits of using EXCELL RSC physical approach to the conversion of rainfall statistics. Furthermore, several methods based on physical approach and empirical ways are studied but when physical approach is used, all the input parameters needed for analysis is unavailable. Similarly, mathematical theory with is based on first principles for deintegrating T-min experimental probability distribution (PD) into corresponding 1-min PD are studied (Matricciani 2011). However, the contribution present the need for more efficient propagation planning, based on the use of number and duration of rain events along with the fraction of rainy time. The experimental system carried out by KMA provides the record for only experimental 1-min rainfall amount as discussed in further section of this paper. Due to the ease of simplicity and easier analysis purpose, empirical nature of rainfall rate method is chosen (Emiliani et al. 2009). The rain rate characteristic of the South Korea was studied to predict 1-min rain rate statistic in Jung et al. (2008), which was based on the 2 years of rain events. Unfortunately, it was found to be less effective and generated higher error percentage. Similarly, with three years of rainfall data, a conversion method for rainfall rate with various integration time was proposed (Lee et al. 1994) base on linear and logarithmic approach. In addition, the global empirical approach was also analyzed in Jung et al. (2008), which was compared with older version of ITU-R P.837-5 rain rate model. A model for rain drop size distribution (Park et al. 2002), was introduced which describe extended gamma distribution function.
In this scenario, this paper presents a novel work for prediction of 1-min rainfall rate. Empirical methods is equally been studied in other countries as noted in Ong and Zhu (1997), Singh et al. (2006Singh et al. ( , 2007, Mandeep and Hassan (2008), Segal (1986), Burgueno et al. (1988) from long term rainfall database. In addition, the global coefficients values are listed in Emiliani et al. (2008) which extend its application to rain rate conversion methods in temperate, tropical and cold climates.
The performance of proposed ITU-R P.837-6 (2012) method for derivation of 1-min rain rate is compared with globally applicable empirical methods namely Segal (1986), Burgueno et al. (1988), Chebil and Rahman (1999), logarithmic (Lee et al. 1994) along with global coefficients approach and polynomial fit analysis of the rainfall rate data. The use of polynomial relationship was found to be better for derivation of 1-min integration time as shown in Khairolanuar et al. (2014), Owolawi (2011). The brief overviews of applied models are presented below:

Segal method
This method was developed based on database of high resolution rainfall records prepared at the Communications Research Centre. The rainfall records were taken from ten years of daily tipping bucket rain gauge charts for each of the 47 stations in Canada. The conversion method is expressed as (Segal 1986): where, conversion factor, ρ τ (P) = aP b , R 1 (P) represents the rainfall rate in a 1-min integration time with the possibility of occurrence P, R τ (P) is the rainfall rate in τ-minutes integration time, and parameters a and b are regression coefficients that are derived from statistical analysis of rainfall data.

Burgueno et al. method
Based on 49 years of rainfall data measured at Barcelona, Spain, Burgueno et al used direct power law fit as (Burgueno et al. 1988): where R 1 (P) and R τ (P) are the precipitation rates with a sampling interval of 1-and τ-min respectively with equal probability of time percentage P, a and b represent the conversion variables.

Chebil and Rahman method
Chebil and Rahman introduced an experimental technique for estimating the precipitation rate conversion element by using the conversion process from 60-to 1-min integration time as (Chebil and Rahman 1999): where R 60 (P) is the precipitation rate in 60-min integration time. ρ 60 (P) is expressed as a mixed power-exponential law, ρ 60 (P) = aP b + ce (dP) with regression variables represented by a, b, c and d analyzed from statistical analysis of rainfall data. Suitability of this method has been further tested for other lower integration time intervals. (1) (3) ρ 60 (P) = R 1 (P)/R 60 (P)

Logarithmic model
The expression for this model is given as (Lee et al. 1994): where a is the regression variable derived from statistical analysis of rainfall rate.

ITU-R recommended model
The most acknowledge model of the latter kind included in Study Group 3 Report for the update of Rec. ITU-R P.837-5, Annex 3 (2012), which allows global 1-min rain rate prediction from the knowledge of the local rainfall rate. EXCELL Rainfall Statistics Conversion (ERSC) (henceforth EXCELL RSC) method is used for conversion of rainfall rate statistics from long to 1-min integration time. This method is based on the simulated movement of rain cells over a virtual rain gauge, with given integration time T, whose translation velocity depends both on the type of precipitation and on the observation period. The conversion of rainfall was obtained using a virtual rain gauge according to the local mean yearly wind velocity, which as extracted from the ERA-40 database. The model goes through an iterative inversion procedure that aims at identifying the local P(R) 1 (a set of P 0 , n, R a ) which when used as input to the rain gauge simulator, provides the best possible estimate of the measured P(R) T . The detail of this approach can be obtained from description of method adopted for the update of Rec. ITU-R P.837-5, Annex 3.

Proposed model
The new model is based on the curve fitting technique analyze from Matlab. Polynomial in one variable is expressed as a n x n + a n−1 x n−1 + · · · + a 1 x 1 + a 0 where x is a variable and exponents are non negative integers with real number coefficients and a n � = 0. A function in the form f(x) = a n x n + a n−1 x n−1 + · · · + a 1 x 1 + a 0 is a polynomial function. Rain cells characteristics have been found to be better shown by an exponential profile which is able to represent real single-peaked rain structure (Luini and Capsoni 2011). The rainfall pattern is observed to be visible through exponential coefficients. This paper present a comparison with several polynomial functions and proposed Model 1 of modified fourth order polynomial function as a suitable approach for South Korea's own numerical prediction model for 1-min rainfall rate derivation. These models are represented as: ] + c, represent second degree polynomial function (9) R 1 (P) = a[R τ (P)] 3 +b[R τ (P)] 2 +c[R τ (P)]+d, represent third degree polynomial function where R 1 (P) represents the rainfall rate in a 1-min integration time with the possibility of occurrence P, R τ (P) is the rainfall rate in τ-minutes integration time, and coefficients a, b, c and d are regression coefficients that are obtained through statistical analysis of rainfall data with the use of curve fitting technique derived from Matlab programming. The effectiveness of proposed Model 1 are measured through various error analyses which highlight better agreement with experimental 1-min rainfall data obtained from KMA for listed nine different regions of South Korea. The importance of constant values and modification on coefficients are judged through regression analysis and curve fitting approach which give minimum error values against the calculated 1-min rain rate values from experimental rainfall amount collection under 1-min duration.

Methodology and analyses of experimental data
The Republic of Korea lies in temperate zone with four distinct seasons. Geographically, country is located in the middle latitudes of the Northern Hemisphere, on the east coast of the Eurasian Continent and also adjacent to the Western Pacific. Thus, complex climate characteristics are observed in this belt which reveals both continental and oceanic features. The seasonal climate characteristic of the South Korea is shown in Table 1.
Rainfall data is effectively recorded by KMA, a central governmental organization of the Republic of Korea under the Ministry of Environment (MOE), which has developed a digital system for accurate measurement of 1-min rainfall amount since 2004 through the use of Tipping bucket rain gauge over several sites some of which are shown in Fig. 1 with the intensity to develop the South Korea's own numerical prediction model. Furthermore, Table 2 indicates the selected sites with their geographical co-ordinate values.
Tipping bucket rain gauge, which is used throughout the measurement sites, automatically record rainfall and facilitate the digitization of telemetric observation signals. The heater is installed inside the sewer for measurement under snowfall. Figure 2 shows the internal structure of rain gauge. KMA uses conducting vessel size of 0.5 mm to improve the shortcomings of the gutter. Once the collected water is more than 0.5 mm it eventually, fills the bucket. Bucket is mounted in particular axis of rotation to shift the center like a seesaw. This bucket is in contact with the Reed Switch with the rotation axis which is operated by electrical pulse occurred due to tipping phenomenon. Finally, signal generated through Reed Switch is recorded on recording device which provide measurement of 1-min rainfall amount. Table 3 shows specification of rain gauge used for recording of 1-min rainfall amount. The tipping bucket has unstably balanced twin-bucket with sensitivity of 0.1 mm per tips which trigger an electronic impulse and is stored in the data logger which scans the data at an interval over 1-min. The availability of the gauge is about 99.2 %. The 0.8 %unavailability is due to system maintenance. Figure 3 shows the overall operation of experimental system used for rainfall amount data logging parts where the accumulated rainfall amount data is first collected in a data logger which record the number of tips for every 1-min interval which is then converted to rainfall amount and finally stored in data storage devices. ITU rain attenuation prediction method, Series (2015), is based on 0.01 % of a time for rain rate parameter. Similarly, Recommendations ITUR (2015) emphasizes on use of rain rate exceeded for 0.01 % of the time with an integration time of 1-min for long term statistics of rain attenuation. The importance of 1-min rain rate has been further studied for satellite and terrestrial rain attenuation predictions (Abdulrahman et al.  and Chuncheon sites respectively. Unfortunately, ITU-R 837-6 (2012) model prediction overestimate the rain rate of these mentioned sites and does not satisfies the local statistical data. The received decade 1-min rainfall amounts are arranged for other time integrations period as 5-, 10-, 20-, 30-and 60-min. These data are sorted and required rainfall amount are extracted for several times percentages as mentioned by ITU-R P. 311-15 (2015). For example, for Gwangju site, at 0.01 % of time, 1-min rainfall amount for 526 ((10 × 365 × 24 × 60 × 0.01)/100) instance was taken which is converted to rain rate expression as mentioned in Kestwal et al. (2014). The result of 1-, 5-, 10-, 20-, 30-and 60-min rainfall rate of selected nine major cities as Gwangju, Daegu, Daejeon, Busan, Seogwipo, Seoul, Ulsan, Incheon and Chuncheon, are calculated from experimental 1-min rainfall amount as obtained from KMA which are summarized in Table 4. These data are used as a basis for comparison with estimated 1-min rain rate from ITU-R P.837-6 (2012) and to propose local prediction model. In order to generate estimated 1-min rain rate from the software as recommended by ITU-R P.837-6 Annex 3, we have selected Mode A operational mode in which rain rate values for several source integration times as stated in Table 4 along with sites latitude and longitude information as mentioned in Table 2

Results and discussion
In order to better visualize the 1-min rain rate distribution data against ITU-R P.837-6 (2012) predicted values for several time percentage over the nine sites in the South Korea at equiprobable exceedance probability (0.001 ≤ P ≤ 0.1 %), CCDFs of rain rate are plotted which are shown from Figs. 4a-c, 5a, b. The plots include the value of Daejeon (also called as Taejon) site, which was chosen by ITU-R P.837-6 (2012) model in its ERA-40 database (Uppala et al. 2005). For the experimental purpose, Taejon site is included within Daejeon site by KMA. ITU-R P.837-6 (2012) model does not accurately predict the 1-min rain rate distribution for nine sites even though this model shows fair well statistics at lower time conversion. In addition, this model dramatically overestimate 1-min rain rate at higher time conversion. For instance, as depicted in Fig. 4a for 5to 1-min integration time, ITU-R P.837-6 (2012) models seems to give fairy satisfactory result but chances of error is still remain, particularly for Seogwipo site. As integration times increase to 10-, 20-, 30-, and 60-min the probability of overestimating 1-min rain rate is increased as highlighted from Figs. 4b, c, 5a, b. The reason behind this difference could be the fact that the matrices used to obtain the parameters might have low spatial resolution. This indicates that ITU-R P.837-6 (2012) model performance statistics does not shows good pattern with calculated rainfall rate from experimental 1-min rainfall amount of the South Korea regions. In this concern, there is the immediate need for 1-min rain rate prediction model that can show greater efficiency against the local 1-min Fig. 4 a Cumulative distributions of 1-min rainfall rate compared with 5-min integration rainfall rate data. b Cumulative distributions of 1-min rainfall rate compared with 10-min integration rainfall rate data. c Cumulative distributions of 1-min rainfall rate compared with 20-min integration rainfall rate data rain rate distribution. Under this scenario, the paper presents new model that shall be applicable in analyzing the 1-min rain fall rate distribution pattern.
As an initial step, regression analysis is performed to match the data to known distribution. Regression analysis is a statistical method to estimate the values of dependent variables that correspond to certain values of new independent variables once the magnitude of the influence of independent variables on dependent variables is measured, thereby determining the regression plane or line with regard to the independent variables. This model summarizes the large amount of data with minimum modeling error (Crane 1996). The regression coefficients applicable for mentioned nine sites are generated through curve fitting approach using Matlab programming whose generated values are listed in Tables 5, 6, 7 and 8 along with the average coefficients generated out of nine   Table 5 continued sites. These regression coefficients are considered useful for obtaining 1-min rain rate when rain rates are available at different integration times, especially when long-term precipitation data from meteorological stations are utilized for obtaining short integration time rain rates for attenuation prediction. Thus using ten rain conversion methods and coefficients from Tables 5, 6, 7 and 8, rain rates at different integration times are converted to 1-min rain rate distribution.
The effectiveness of proposed model is observed from the coefficient of determination, R 2 , values as listed in Table 9. This statistical property of regression concerns the relationship between the PD of the parameter estimates and the true values of those parameters. The coefficient of determination, R 2 , describes the proportion of variance in measured data explained by the models. It is the portion of total variation in dependent variable that is explained by variation in independent variable (Steel and Torrie 1960). R 2 ranges from 0 to 1, with higher values indicating less error variance whose values are summarized in Table 9. Out of ten mentioned empirical methods only Burgueno et al., second order polynomial fit, third order polynomial fit, Model 1, Model 2 and Model 3 values are listed because of dependability on statistical analyses for regression values.
As noted from Table 9, the average regression values obtained while applying Model 1 are 0.9819, 0.9772, 0.9733, 0.9621 and 0.9495 for 5-, 10-, 20-, 30-and 60-to 1-min conversion times respectively. These values are closer to unity as observed against the other applied models. Hence, the proposed model gives less chances of error variance.

Evaluation of proposed method
In order to measure the goodness of fit of proposed model, paper present several error analyses. Mean, standard deviation (SD) and root mean square (RMS) values of error probability, ε(P), are gathered, where they are compared to the performance of the ITU-R P.   Table 7 Regression coefficients at different integration times for three empirical methods     Table 7 continued    Absolute percentage relative error figure is given as,

Table 8 Regression coefficients at different integration times for two empirical methods
where R e (P) T and R m (P) T are the rain rate values of the estimated and the measured T-min integrated rainfall CDF, respectively, at the same probability level P, in the percentage interval 10 −3 < P < 10 −1 %. RMS as defined by Owolawi and Afullo (2007), where N is the total number of available probability values, X est and X mea are the estimated and measured quantities respectively. Similarly, SD, is calculated as, where N is the total number of available probability values, ε(P) i and μ are each error value and arithmetic mean of error quantities respectively. The calculated error probabilities are presented in tabular form to accurately identify the obtained error values. The average error values thus obtained over all integration times using regional coefficient sets for each of the methods are listed in Table 10,    which indicates that proposed model gives less relative error percentages which are <1 % for all conversion times. This is justified from lower values of SD and RMS calculation which is <7 % and 6 % in aggregate for all integration times respectively. In contrast, ITU-R P.837-6 (2012) produces higher error percentage of 5. 19, 13.73, 43.23, 71.19 and 115.36 % for 5-, 10-, 20-, 30-and 60-to 1-min conversion time respectively. Hence, proposed model provides a better accuracy for all integration times. In addition, polynomial fits of third and second orders can be considered as a second and third preferred method because these models also show lower chances of error probabilities. Furthermore, Table 11 presents the results of evaluation for average error probability using regional coefficient sets at 0.01 %, which is considered to be suitable time percentage for calculation of rain induced attenuation (ITU-R 2015) and is very crucial for system designers to obtain preliminary design of the satellite microwave link, satellite payload design and to have broad idea of rain attenuation for microwave engineers.
Considering the variability of the rain rate predictions at 0.01 % time exceedance, the third order polynomial fit does better result presenting the relative error values of −1.06, 0.1, −0.49, −0.37 and −1.81 % for 5-, 10-, 20-, 30-and 60-1-min conversion times respectively. This is supported by less value of RMS errors as presented in Table 11. Interesting, fair variability in error chances are observed from proposed model against the third order polynomial fit. Under this condition, proposed model can be consider as a second choice which is followed by second order polynomial fit as an third preferred model at 0.01 % of time. In contrary, ITU-R P.837-6 and global coefficients method provide high error values for all time percentages. Even though, these models shows low values of error probability at lower time percentage especially, 5-to 1-min conversion time, but still this error is higher than other prediction models. In order to verify the prediction performance of the models, relative error percentages along with SD and RMS values are calculated using the average coefficients sets. Table 12 shows the average values of relative error percentage, SD and RMS results as obtained by using the average coefficient sets for all the measurement sites.
As noted from Table 12, the proposed model, have good result for 5-to 1-min time conversion and error probability dramatically increases for other higher times conversion especially at 60-to 1-min conversion time where negative values of rain rate are obtained. This might be due to the exponential function that it includes. Similarly, ITU-R P.837-6, Burgueno et al., global coefficients and model 3 show higher error chances and are impracticable to use for 1-min rain rate derivation. In addition, Chebil and Rahman method and Third order polynomial fit produced increased error chances as per the increasing conversion times. Segal method and second order polynomial fit result in similar nature of error probabilities. Interestingly, logarithmic model is found to be best while considering the average coefficient sets because the error probabilities are lower as 1. 01, 5.42, 9.15, 7.06 and 6.08 % for 5-, 10-, 20-, 30-, and 60-to 1-min conversion times respectively. This is verified from lower SD and RMS values of 9. 08, 28.75, 28.53, 23.66 and 21.48 % along with 10.47, 36.63, 36.07, 31.11 and 29.96 % for 5-, 10-, 20-, 30-and 60to 1-min respectively. Additionally, performances of models are graded at 0.01 % of time while using average coefficient sets whose results are shown in Table 13. As noted, third order polynomial fit produces less error chances as compared to other models. Table 11 Mean error obtained after testing over the 0.01 % of time      Table 14, third order polynomial model produces relatively less error probabilities 0.65,0.7,4.43 and 0.96 % along with RMS values of 4.49,4.55,8.1,-to 1-min respectively which indicates the model suitability. The second most suitable model shall be Logarithmic model with similar error probabilities values. In other hand, ITU-R P.837-6, Burgueno et al., Chebil and Rahman, global coefficients and model 3 result in greater error chances so these models are not preferable at this time percentage. In addition, Segal method gives better estimation from 5-to 1-min conversion time but the error chances increase as conversion times get increased. Proposed model 1 produces higher error chances at higher conversion time, especially 60-to 1-min and relative error values are <6 % for lower time conversion. Second order polynomial fit and model 2 result higher error values as compared to third order polynomial fit.

As indicated in
Furthermore, models performance is justified with the application of regional and average coefficient sets in Daejeon site. As previously mentioned, this site is also considered in ERA-40 data base. Daejeon is the fifth largest metropolis with elevation of about 77 m above sea level and average precipitation is above 300 mm during the month of July and August. The calculated rain rate from experimental 1-min rainfall amount in this site at 0.01 % of time is 79.8 mm/h. Performance of error analysis in this site is done by using regional coefficient sets and average coefficient sets which are listed in Tables 14  and 15 respectively while considering all integration times.
As noted from Table 14, proposed model 1 and third order polynomial fit exhibit relatively less error chances which is <1 % for all integration times. In contrast, ITU-R P.837-6 and global coefficients produce higher error probabilities which indicate that it is not suitable for 1-min conversion process.
Moreover, the regional and average coefficients are further tested at 0.01 % of time whose error values are depicted in Tables 16 and 17 respectively. As noted from Table 16, Third order polynomial fit and Segal method produced lower relative error values among which earlier method generate less error chances. This is indicated by Table 16 where Third order polynomial fit give values of −3.13, −2.3, 0.04, 3.58 and −1.27 % for 5-, 10-, 20-, 30-and 60-to 1-min conversion times respectively. This is supported by lower RMS values of 2.5, 1.84, 0.03, 2.86 and 1.01 % for 5-, 10-, 20-, 30-and 60-to 1-min conversion times respectively. In contrast, ITU-R P.837-6 and global coefficients generate larger error values for all integration times as indicted in Table 16. Table 14 Error obtained after testing over the interval [0.001-0.1 %] using regional coefficient sets      Table 16 Error obtained after testing at 0.01 % of time using regional coefficient sets   Table 17 Error obtained after testing at 0.01 % of time using average coefficient sets   . 6 a 1-min rainfall rate compared with 5-min integration time rainfall rate data. b 1-min rainfall rate compared with 10-min integration time rainfall rate data. c 1-min rainfall rate compared with 20-min integration time rainfall rate data. d 1-min rainfall rate compared with 30-min integration time rainfall rate data. e 1-min rainfall rate compared with 60-min integration time rainfall rate data As indicated in Table 17, logarithmic method and Second order polynomial fit give less error values among which former model shows less error variance for higher time integration. Although second and third order polynomial show better result at lower integration times especially 5-, 10-, 20-and 30-to 1-min but for higher integration time their error chances are increased. Hence, Logarithmic model is more preferable at 0.01 % of time. In contrary, ITU-R P.837-6, Burgueno et al., Chebil and Rahman methods and model 3 generate higher error values as depicted in Table 17 which signifies their unsuitability for prediction of 1-min rain rate at 0.01 % of time.
Estimated 1-min rainfall rate curves for Daejeon Sites from 5-, 10-, 20-, 30-and 60-to 1-min conversion times are highlighted in Fig. 6a-e respectively. These plots indicate that empirical nature of models markedly follows the calculated 1-min rain rate pattern. Figure 6a, b indicate that ITU-R P.837-6 method under estimate calculated 1-min rain rate for lower times conversion especially 5-and 10-to 1-min. Additionally, Fig. 6c-e highlight the overestimation shown by ITU-R P.837-6 method and global coefficient approach against the calculated 1-min rain rate.

Conclusion
The results obtained will serve as effective tools for communication system designer to understand the effect of rain in propagation medium. In this regard, this paper estimates the suitable empirical conversion model based on a decade long rainfall distribution data from KMA (2004KMA ( -2013 over nine regions of the South Korea. Several rainfall rate conversion processes have been carried out for various integration times using these cumulative rainfall distribution. Rain rates exceeded for 0.001 < P < 0.1 % of the time was compared with prominent rain rate models analyzed in nine sites rainfall data. In addition, specific to 0.01 % time percentage value, several error metrics are evaluated. The performance criteria are based on the estimated statistics for Error percentage, SD and root mean square error over the nine different region's rainfall data. The calculated 1-min rainfall rate from experimental 1-min rainfall amount is compared with ITU-R P.837-6 method and five existing empirical rainfall rate models along with five different polynomial fits. ITU-R P.837-6 method underestimate 1-min rain rate at lower integration times especially 5-, 10-to 1-min conversion time and overestimate at higher integration times mostly 20-, 30-and 60-to 1-min. Under regional coefficients set, proposed model 1, show better estimation of 1-min rain rate and at 0.01 % of time, third order polynomial fit hold satisfactory result with less error probabilities. Interestingly, proposed model 1 and third order polynomial fit shows same condition in Daejeon site. Similarly, while using average coefficients set, Logarithmic model hold better estimation of 1-min rain rate and at 0.01 % of time, Logarithmic model, second and third order polynomial fits give satisfactory result. Same condition is noted in Daejeon site too under average coefficients set approach. In overall, the paper emphasizes that ITU-R P.837-6 global 1-min rain rate estimation performance did not significantly reflect the South Korea's local rainfall characteristics. On the basis of overall result, it can be concluded that the proposed model seems to provide a better and more reliable alternative to the ITU-R P.837-6 method for better estimation of 1-min rainfall rate. We hope this work will be milestone