^{1}

^{1}

^{*}

^{1}

^{*}

Considering the rainfall’s importance in hydrological modeling, the objective of this study was the performance comparison, in convergence terms, of techniques often used to estimate the average rainfall over an area: Thiessen Polygon (TP) Method; Reciprocal Distance Squared (RDS) Method; Kriging Method (KM) and Multiquadric Equations (ME) Method. The comparison was done indirectly, using GORE and BALANCE index to assess the convergence results from each method by increasing the rain gauges density in a region, through six scenarios. The Coremas/Mae D’água Watershed employed as study area, with an area of 8385 km^{2}, is situated on Brazilian semi-arid. The results showed the TP, as RDS and ME techniques to be employed successfully to obtain the average rainfall over an area, highlighting the MEM. On the other hand, KM, using two variograms models, had an unstable behavior, pointing the prior study of data and variogram’s choice as a need to practical applying.

The average rainfall over an area may be considered as the main input on watershed modeling process, especially of those which deal with surface runoff, partly because, in general, the rain is the only climatic variable that can explain fast increasing flow [1-3]. Still, several studies show that spatial variability of rainfall over the basin and their distribution pattern, as well as its interaction with the basin, have considerable effect on runoff response generated [4-7].

In this sense, despite the development of radar technology has been observed, due to the limitations and characteristics of theses specialized measurements, point measurements made by rain gauge are still required for better modeling [1,8-10]. Besides, the point measures of rain are, in many places, the only available time-series source with enough spatial density for hydrological studies, so to runoff modeling, as well as to water resources planning. Thus, the study of techniques to estimate the average rainfall, using point measures, and the distribution pattern analysis, remains indispensable.

Considering the methods used to determine the average rainfall, which vary from the simple linear combination to geostatistics techniques, it’s possible to emphasize in many studies and applications: the Thiessen Polygon (TP) Method [11,12], the Reciprocal Distance Squared (RDS) Method [13-16] the Kriging Method (KM) [5, 17-19] and the Multiquadric Equations (ME) Method [16,20-25].

Besides the mean rainfall, RDS, KM and ME also estimate a continuous surface that adjusts itself to the known rainfall values, being useful on point rainfall values determination within the basin. This particular characteristic is very useful on rainfall spatial distribution evaluation [

However, the real value of average rainfall and its distribution is still an unknown variable [

Given that, the aim of this study is the comparison of techniques to estimate the average annual rainfall over a given watershed, through the analysis of results in various scenarios of data availability, in view of the water resources planning, for long and medium term, in regions with point measures predomination of. For this, it has chosen the Coremas Mãe D’Água Watershed, located in Brazilian semi-arid (

The semi-arid climatic aspect adds more peculiarities to water resources planning because of the high variability, in time and space, of rainfall [

It is not possible to compare values from any method with the real one, because it is unknown. It is possible to evaluate the variability (or convergence) of some model caused for a change in availability data. Thus, even that methods remain incomparable directly, find which of them would need less spatial data to achieve results supposedly more reliable (obtained in a better data scenario), it is feasible, and so comparing them indirectly, following this line: in a favorable situation which the rain gauge density was great enough. Probably all techniques would give great results, very close to real values of mean rainfall; so, it is reasonable to believe that methods with a good convergence behavior, at first, are more reliable in a less favorable situation. [

from the rainfall-runoff modeling perspective [1,6,7,24, 25,27]. In some cases, it was found that reducing the number of rain gauges used could improve the models accuracy, inclusive. However, in the water resources management case for great scales of time and space, it’s thought that a greater dataset will reflect in a better mean rainfall estimate and its behavior which, in some cases, is the objective of the specific study.

Thus, the previously cited methods were compared with the two indices, using five scenarios with a different spatial data density, comparing their results with those obtained from a sixth scenario with the major density available, the reference set.

For this task, 56 rain gauge stations were selected within the study area with daily records from ten consecutive years, from 1965 to 1974, being admitting just one missing year. These registers were obtained from ANA database. The gap filling was done just with the simple mean and after this the annual series for each gauge were constructed.

The Watershed of Coremas Mãe D’Água Dam was chosen as study area for this work (

From 56 selected stations, six scenarios were made, increasing the number of gauges until the maximum, as showed in

The GORE index [_{i} be the real average rainfall in a given interval of time (here considering being equal to results obtained in Scenario 6, because it’s unknown in practice) and let P_{i}^{E} be the estimated value of rainfall for the same interval i in a given scenario, so:

n is the time intervals number and P is the mean of P_{i} values for all time considered.

Developed by [

[

The application of quadric surfaces for points data interpolation was initially developed by [

Based on regionalized variables concept, developed by [

The average rainfall on the watershed for each year and scenario (1 to 6), estimated by each method tested, is given in

by each method, comparing Scenarios 1 to 5 with the reference 6. From

It can observe that all methods, at Scenario 1, gave not so good results, with negative values. In general, indices had a trend of improvement when the number of rain gauges was increased, but it should be noted the irregular behavior presented by both KG and KCS. Despite the KCS having shown the best results in Scenario 1, it just returned to give good results in Scenario 4, but not so good as TP, RDS nor ME. On the other hand, KG showed good results for Scenario 2, but with the worse estimated results for Scenarios 4 and 5, demonstrating a significant instability. Thus, differently of other techniques, the KM was the only method that not presented the improvement behavior in results expected when increasing the data availability. So, results demonstrates that applying of geostatistical techniques on rainfall data needs preliminary studies of data employed and variogram model applied.

It can also verify that TP got an excellent performance from Scenario 2 onwards, for both indices. At this point, it should be noted that care should be taken with this technique. TP ponders the rainfall measures based on the area of influence of each station within the basin, which implies that a good homogeneity combined with an enough density may result reliable values for average rainfall, as seen here. However, in tiny scales of time or bad spatial distribution of data, the rainfall variability, or even the existence of error in the records, can contribute greatly to a discrepant with reality because of method formulation. Moreover, the use of the TP is not suitable for the estimation of rain in a certain region or point in the basin.

ME and RDS also presented great results, very similar, with some emphasis for ME at Scenarios 2, 4 and 5. Both showed the expected behavior improves, as TP did, when increasing the stations numbers. The ME has the advantage, given its theoretical base, that weights estimation, for all stations, are determined simultaneously, which allows possible registry errors may be diluted, resulting in a rainfall surface more reliable than that given by RDS. In other words, ME assesses the spatial structure of the events. As to results, in specific for BALANCE index, RDS obtained better results just at scenarios 1 and 3, pointing out that in the first one all methods were flawed in determining the average rainfall, underestimating its value significantly.

The direct comparison among techniques of average rainfall estimation is not possible because the real values are generally unknown. However, this research brought another approach in order to compare indirectly some methods, using not just the results obtained, sometimes very similar, but analyzing the data requiring each one to reach better results given in the best spatial data scenario, in other words, comparing their convergence behavior using GORE and BALANCE indices. The expected behavior from each technique was the continuous improvement in estimated results when increasing spatial data density. Thus, it is reasonable to infer that some methods are appropriate, with a good performance, when their behavior shows that, even with a small amount of data in space, their results approach those that would be given in better conditions.

The results for GORE and BALANCE indices, indicate that TP, RDS and ME as methods are to be applied with satisfaction to obtain the average rainfall value over an area. On the other hand, KM, tested with two variograms models, had a not expected unstable behavior.

Reflecting the need of preliminary studies about data and variograms to be applied, it means a disadvantage given by an increased complexity, especially from the point of view of procedures automation and management tools.

Returning to the methods with a good performance, emphasis must be given to ME, by the great results obtained and its formulation, which helps to mitigate casual data errors allowing estimating a continuous rainfall surface. Still, the study was done in a semi-arid region of Brazil, where pluvial behavior presents a high variability, even in larger time scales, reinforcing the results reached.

It is suggested that more studies be done in this way of indirect comparison of techniques, despite the technological advances, there are many regions where the rainfall monitoring is still scarce and there is a need for reliable water resources planning tools. In the same vein presented here, larger scenarios combinations may be done, using different watersheds and spatial and temporal discretization, so that the methods could be evaluated under various conditions.

The authors are deeply grateful to CNPq for Masters Scholarship to Fernando D. Barbalho (Proc. 556625/ 2009-9) and a Productivity Research Grant to Klebber T. M. Formiga (Proc. 310389/2012-7); as well as to the FINEP/CT-HIDRO/PROCESSOS HIDRÁULICOS/2007 to support this research and Federal University of Goiás (UFG) for providing infrastructural conditions for the development of this research.