An improved statistical approach to merge satellite rainfall estimates and raingauge data.

Summary Deriving high quality daily rainfall estimates are required not only for successful hydrological modelling but also for its application in ungauged basins. At present, there are two resources to estimate rainfall fields: raingauge stations and remote sensing devices (such as satellites and radar). The observations obtained from raingauges are accurate and reliable, but are usually limited to interpolate daily rainfall as input into hydrological models by spatial coverage. Hence, the areal precipitation estimates constructed solely by raingauges exhibit a great deal of uncertainty where the raingauge network is sparse, especially for local convective events. Satellite-based gridded precipitation estimates are produced at high temporal (3 h) and reasonably high spatial resolution (0.25° × 0.25°). However, the accuracy of these surfaces is lower than using ground based observations. The main aim of this paper is to present a novel method of generating a synthesis of a discrete set of point raingauge observations and a satellite derived rainfall product. This method stems from some ideas of data assimilation and focuses on the difference between gridded and point data. A nonparametric kernel smoothing method is employed in this merging strategy with emphasis on discontinuity correction and spatial interpolation adapting for sparse design. A cross-validation study was undertaken to blend observations from the Australian raingauge network and satellite derived TRMM Multisatellite Precipitation Analysis (TMPA) 3B42. In order to test the proposed method, several methods including ordinary kriging and cokriging were trialled. The comparison showed that the bias results produced by the proposed method were far superior to the other methods trialled and other error statistics such as root mean square error were comparable or better (although the kriging methods tended to perform better in densely gauged areas). In particular, the bias of the proposed method was reduced consistently over the whole country, whilst the kriging methods trialled produced relatively large biases across the country leading to near overall zero bias. In addition to improving quantitative efficiency, the estimated rainfall field showed good visual performance in the presence of sparse gauge stations.

[1]  T. N. Krishnamurti,et al.  The status of the tropical rainfall measuring mission (TRMM) after two years in orbit , 2000 .

[2]  Ezio Todini,et al.  A Bayesian technique for conditioning radar precipitation estimates to rain-gauge measurements , 2001 .

[3]  Frédéric Fabry,et al.  The effect of gauge sampling density on the accuracy of streamflow prediction for rural catchments , 1993 .

[4]  W. Tobler Smooth pycnophylactic interpolation for geographical regions. , 1979, Journal of the American Statistical Association.

[5]  Theo Gasser,et al.  Finite-Sample Variance of Local Polynomials: Analysis and Solutions , 1996 .

[6]  H. Müller,et al.  Kernels for Nonparametric Curve Estimation , 1985 .

[7]  Witold F. Krajewski,et al.  Cokriging radar‐rainfall and rain gage data , 1987 .

[8]  Terry A. Slocum Thematic Cartography and Geographic Visualization , 2004 .

[9]  G. P. Cressman AN OPERATIONAL OBJECTIVE ANALYSIS SYSTEM , 1959 .

[10]  Berwin A. Turlach,et al.  Interpolation Methods for Adapting to Sparse Design in Nonparametric Regression , 1997 .

[11]  C. K. Chu,et al.  An interpolation method for adapting to sparse design in multivariate nonparametric regression , 2003 .

[12]  Dong-Jun Seo,et al.  Real-time estimation of rainfall fields using rain gage data under fractional coverage conditions , 1998 .

[13]  Interpolation methods for adapting to sparse design in nonparametric regression. Comment. Rejoinder , 1997 .

[14]  Yang Hong,et al.  Evaluation of TRMM Multisatellite Precipitation Analysis (TMPA) and Its Utility in Hydrologic Prediction in the La Plata Basin , 2008 .

[15]  K. Crawford,et al.  Improving WSR-88D Hourly Rainfall Estimates , 1998 .

[16]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[17]  Edzer J. Pebesma,et al.  Multivariable geostatistics in S: the gstat package , 2004, Comput. Geosci..

[18]  J. Janowiak,et al.  COMPARISON OF NEAR-REAL-TIME PRECIPITATION ESTIMATES FROM SATELLITE OBSERVATIONS AND NUMERICAL MODELS , 2007 .

[19]  A. E. Hoerl,et al.  Ridge regression: biased estimation for nonorthogonal problems , 2000 .

[20]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[21]  S. Barnes,et al.  A Technique for Maximizing Details in Numerical Weather Map Analysis , 1964 .

[22]  John O. Carter,et al.  Using spatial interpolation to construct a comprehensive archive of Australian climate data , 2001, Environ. Model. Softw..

[23]  Dong-Jun Seo,et al.  Real-time estimation of rainfall fields using radar rainfall and rain gage data , 1998 .

[24]  Y. Hong,et al.  The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-Global, Multiyear, Combined-Sensor Precipitation Estimates at Fine Scales , 2007 .

[25]  David R. Brillinger,et al.  Examples of Scientific Problems and Data Analyses in Demography, Neurophysiology, and Seismology , 1994 .

[26]  W. Härdle,et al.  Statistical Theory and Computational Aspects of Smoothing , 1996 .

[27]  Amvrossios C. Bagtzoglou,et al.  Satellite-based Flood Modeling Using TRMM-based Rainfall Products , 2007, Sensors.

[28]  G. Pegram,et al.  Combining radar and rain gauge rainfall estimates using conditional merging , 2005 .

[29]  Walter Collischonn,et al.  Daily hydrological modeling in the Amazon basin using TRMM rainfall estimates , 2008 .

[30]  David L. Toll,et al.  Statistical Evaluation of Combined Daily Gauge Observations and Rainfall Satellite Estimates over Continental South America , 2009 .

[31]  C. Gotway,et al.  Combining Incompatible Spatial Data , 2002 .

[32]  Dawn Harrison,et al.  Improving precipitation estimates from weather radar using quality control and correction techniques , 2000 .

[33]  P. Kyriakidis A Geostatistical Framework for Area-to-Point Spatial Interpolation , 2004 .

[34]  R. Hwang,et al.  DOUBLE SMOOTHING ESTIMATION OF THE MULTIVARIATE REGRESSION FUNCTION IN NONPARAMETRIC REGRESSION , 2002 .

[35]  Laurent Delobbe,et al.  Evaluation of radar-gauge merging methods for quantitative precipitation estimates , 2009 .

[36]  J. Zajaczkowski A comparison of the BAWAP and SILO spatially interpolated daily rainfall datasets , 2009 .

[37]  S. Barnes,et al.  Mesoscale objective map analysis using weighted time-series observations , 1973 .

[38]  J. Janowiak,et al.  CMORPH: A Method that Produces Global Precipitation Estimates from Passive Microwave and Infrared Data at High Spatial and Temporal Resolution , 2004 .

[39]  Phillip Jordan,et al.  Rainfall-runoff modelling across the Murray-Darling Basin , 2008 .