Radar rainfall estimation for flash flood forecasting in small urban watersheds

Radar rainfall estimation for flash flood forecasting in small, urban catchments is examined through analyses of radar, rain gage and discharge observations from the 14.3 km² Dead Run drainage basin in Baltimore County, Maryland. The flash flood forecasting problem pushes the envelope of rainfall estimation to time and space scales that are commensurate with the scales at which the fundamental governing laws of land surface processes are derived. Analyses of radar rainfall estimates are based on volume scan WSR-88D reflectivity observations for 36 storms during the period 2003–2005. Gage-radar analyses show large spatial variability of storm total rainfall over the 14.3 km² basin for flash flood producing storms. The ability to capture the detailed spatial variation of rainfall for flash flood producing storms by WSR-88D rainfall estimates varies markedly from event to event. As spatial scale decreases from the 14.3 km² scale of the Dead Run watershed to 1 km² (and the characteristic time scale of flash flood producing rainfall decreases from 1 h to 15 min) the predictability of flash flood response from WSR-88D rainfall estimates decreases sharply. Storm to storm variability of multiplicative bias in storm total rainfall estimates is a dominant element of the error structure of radar rainfall estimates, and it varies systematically over the warm season and with flood magnitude. Analyses of the 7 July 2004 and 28 June 2005 storms illustrate microphysical and dynamical controls on radar estimation error for extreme flash flood producing storms.     

Hydrological appraisal of rainfall estimates from radar,satellite, raingauge and satellite–gauge combination on the Qinhuai River Basin,China

ABSTRACT

Multisource rainfall products can be used to overcome the absence of gauged precipitation data for hydrological applications. This study aims to evaluate rainfall estimates from the Chinese S-band weather radar (CINRAD-SA), operational raingauges, multiple satellites (CMORPH, ERA-Interim, GPM, TRMM-3B42RT) and the merged satellite–gauge rainfall products, CMORPH-GC, as inputs to a calibrated probability distribution model (PDM) on the Qinhuai River Basin in Nanjing, China. The Qinhuai is a middle-sized catchment with an area of 799 km². All sources used in this study are capable of recording rainfall at high spatial and temporal resolution (3 h). The discrepancies between satellite and radar data are analysed by statistical comparison with raingauge data. The streamflow simulation results from three flood events suggest that rainfall estimates using CMORPH-GC, TRMM-3B42RT and S-band radar are more accurate than those using the other rainfall sources. These findings indicate the potential to use satellite and radar data as alternatives to raingauge data in hydrological applications for ungauged or poorly gauged basins.     

Microwave link rainfall estimation: Effects of link length and frequency,temporal sampling,power resolution,and wet antenna attenuation

Issues associated with microwave link rainfall estimation such as the effects of spatial and temporal variation in rain, the nonlinearity of R–k$R''–''k$ relations, temporal sampling, power resolution, and wet antenna attenuation are investigated using more than 1.5 years of data from a high-resolution X-band weather radar. Microwave link signals are generated for different link frequencies and lengths from these radar data, so that retrieved path-averaged rainfall intensities can be compared to true path-averaged values. Results of these simulations can be linked to the space–time structure of rain. A frequency-dependent relation between the rainfall intensity at an antenna and the attenuation caused by its wetting is derived using microwave link and rain gauge data. It is shown that if the correct temporal sampling strategy is chosen, the effects of the degradation of power resolution and of wet antenna attenuation (if a correction is applied) are minor (i.e., MBE and bias-corrected RMSE are >−20% and <20% of the mean rainfall intensity, respectively) for link frequencies and lengths above ∼20 GHz and ∼2 km, respectively.     

Influence of rainfall and soil properties spatial aggregation on extreme flash flood response modelling: An evaluation based on the Sesia river basin,North Western Italy

High resolution radar rainfall fields and a distributed hydrologic model are used to evaluate the sensitivity of flood and flash flood simulations to spatial aggregation of rainfall and soil properties at catchment scales ranging from 75 to 983 km². Hydrologic modeling is based on a Hortonian infiltration model and a network-based representation of hillslope and channel flow. The investigation focuses on three extreme flood and flash flood events occurred on the Sesia river basin, North Western Italy, which are analysed by using four aggregation lengths ranging from 1 to 16 km. The influence of rainfall spatial aggregation is examined by using the flow distance as a spatial coordinate, hence emphasising the role of river network in the averaging of space–time rainfall. The effects of reduced and distorted rainfall spatial variability on peak discharge have been found particularly severe for the flash flood events, with peak errors up to 35% for rainfall aggregation of 16 km and at 983 km² catchment size. Effects are particularly remarkable when significant structured rainfall variability combines with relatively important infiltration volumes due to dry initial conditions, as this emphasises the non-linear character of the rainfall–runoff relationship. In general, these results confirm that the correct estimate of rainfall volume is not enough for the accurate reproduction of flash flood events characterised by large and structured rainfall spatial variability, even at catchment scales around 250 km². However, accurate rainfall volume estimation may suffice for less spatially variable flood events. Increasing the soil properties aggregation length exerts similar effects on peak discharge errors as increasing the rainfall aggregation length, for the cases considered here and after rescaling to preserve the rainfall volume. Moreover, peak discharge errors are roughly proportional to runoff volume errors, which indicates that the shape of the flood wave is influenced in a limited way by modifying the detail of the soil property spatial representation. Conversely, rainfall aggregation may exert a pronounced influence on the discharge peak by reshaping the spatial organisation of the runoff volumes and without a comparable impact on the runoff volumes.     

Improving NWP through radar rainfall-driven land surface parameters: A case study on convective precipitation forecasting

In this study we investigate the effect of forcing the land surface scheme of an atmospheric mesoscale model with radar rainfall data instead of the model-generated rainfall fields. The goal is to provide improved surface conditions for the atmospheric model in order to achieve accurate simulations of the mesoscale circulations that can significantly affect the timing, distribution and intensity of convective precipitation. The performance of the approach is evaluated in a set of numerical experiments on the basis of a 2-day-long mesoscale convective system that occurred over the US Great Plains in July 2004. The experimental design includes multiple runs covering a variety of forcing periods. Continuous data integration was initially used to investigate the sensitivity of the model’s performance in varying soil state conditions, while shorter time windows prior to the storm event were utilized to assess the effectiveness of the procedure for improving convective precipitation forecasting. Results indicate that continuous integration of radar rainfall data brings the simulated precipitation fields closer to the observed ones, as compared to the control simulation. The precipitation forecasts (up to 48 h) appear improved also in the cases of shorter integration periods (24 and 36 h), making this technique potentially useful for operational settings of weather forecasting systems. A physical interpretation of the results is provided on the basis of surface moisture and energy exchange.     

Exploring the effect of data assimilation by WRF‐3DVar for numerical rainfall prediction with different types of storm events

The mesoscale Numerical Weather Prediction (NWP) model is gaining popularity among the hydrometeorological community in providing high‐resolution rainfall forecasts at the catchment scale. Although the performance of the model has been verified in capturing the physical processes of severe storm events, the modelling accuracy is negatively affected by significant errors in the initial conditions used to drive the model. Several meteorological investigations have shown that the assimilation of real‐time observations, especially the radar data can help improve the accuracy of the rainfall predictions given by mesoscale NWP models. The aim of this study is to investigate the effect of data assimilation for hydrological applications at the catchment scale. Radar reflectivity together with surface and upper‐air meteorological observations is assimilated into the Weather Research and Forecasting (WRF) model using the three‐dimensional variational data‐assimilation technique. Improvement of the rainfall accumulation and its temporal variation after data assimilation is examined for four storm events in the Brue catchment (135.2 km²) located in southwest England. The storm events are selected with different rainfall distributions in space and time. It is found that the rainfall improvement is most obvious for the events with one‐dimensional evenness in either space or time. The effect of data assimilation is even more significant in the innermost domain which has the finest spatial resolution. However, for the events with two‐dimensional unevenness of rainfall, i.e. the rainfall is concentrated in a small area and in a short time period, the effect of data assimilation is not ideal. WRF fails in capturing the whole process of the highly convective storm with densely concentrated rainfall in a small area and a short time period. A shortened assimilation time interval together with more efficient utilisation of the weather radar data might help improve the effectiveness of data assimilation in such cases. Copyright © 2012 John Wiley & Sons, Ltd.     

Accuracy evaluation of weather data generation and disaggregation methods at finer timescales

Availability of weather data at finer timescales such as hourly is vital in the application of dynamic physical and biological models. In this study, we have examined the suitability of various approaches (deterministic periodic versus stochastic) of disaggregating daily weather data into hourly data in the Cedar Creek watershed, TX, USA. We found the cosine function suitable to disaggregate daily maximum and minimum temperatures and wind speed data into respective hourly data. We also used a common logarithmic equation to compute vapor pressures from temperature data, and hence relative humidity (the ratio between actual and saturated vapor pressures multiplied by 100). Disaggregation following uniform distribution of daily rainfall over 24 h did not reproduce most statistical parameters computed from observed hourly rainfall data onsite. Conversely, both stochastic models formulated based on univariate (Hyetos) and multivariate (MuDRain) processes mimicked the measured hourly rainfall distributions very well. Overall, we found the MuDRain model superior, compared to other models to disaggregate daily rainfall data into hourly data.     

Radar rainfall uncertainty modelling influenced by wind

Radar‐based estimates of rainfall are affected by many sources of uncertainties, which would propagate through the hydrological model when radar rainfall estimates are used as input or initial conditions. An elegant solution to quantify these uncertainties is to model the empirical relationship between radar measurements and rain gauge observations (as the ‘ground reference’). However, most current studies only use a fixed and uniform model to represent the uncertainty of radar rainfall, without consideration of its variation under different synoptic regimes. Wind is such a typical weather factor, as it not only induces error in rain gauge measurements but also causes the raindrops observed by weather radar to drift when they reach the ground. For this reason, as a first attempt, this study introduces the wind field into the uncertainty model and designs the radar rainfall uncertainty model under different wind conditions. We separate the original dataset into three subsamples according to wind speed, which are named as WDI (0–2 m/s), WDII (2–4 m/s) and WDIII (>4 m/s). The multivariate distributed ensemble generator is introduced and established for each subsample. Thirty typical events (10 at each wind range) are selected to explore the behaviours of uncertainty under different wind ranges. In each time step, 500 ensemble members are generated, and the values of 5th to 95th percentile values are used to produce the uncertainty bands. Two basic features of uncertainty bands, namely dispersion and ensemble bias, increase significantly with the growth of wind speed, demonstrating that wind speed plays a considerable role in influencing the behaviour of the uncertainty band. On the basis of these pieces of evidence, we conclude that the radar rainfall uncertainty model established under different wind conditions should be more realistic in representing the radar rainfall uncertainty. This study is only a start in incorporating synoptic regimes into rainfall uncertainty analysis, and a great deal of more effort is still needed to build a realistic and comprehensive uncertainty model for radar rainfall data. Copyright © 2014 John Wiley & Sons, Ltd.     

Analysis and modeling of spatial correlation structure in small-scale rainfall in Central Oklahoma

Spatial correlation structure in small-scale rainfall is analyzed based on a dense cluster of raingauges in Central Oklahoma. This cluster, called the EVAC PicoNet, consists of 53 gauges installed in 25 measurement stations covering an area of about 3 km by 3 km. Two raingauges are placed in 24 stations and five in the central station. Three aspects of the estimated spatial correlation functions are discussed: dependence on time-scale ranging from 1 min to 24 h, inter-storm variability, and dependence on rainfall intensity. The results show a regular dependence of the correlogram parameters on the averaging time-scale, large differences of the correlograms in the individual storms, and the dominance of storms with high spatial variability on the average large sample characteristics. The authors also demonstrate and discuss the ambiguities in correlation estimates conditioned on rainfall intensities. The findings of this study have implications for raingauge network design, rainfall modeling, and conclusive evaluation of radar and satellite estimates of rainfall.     

Coupling between weather radar rainfall data and a distributed hydrological model for real-time flood forecasting / Couplage de données pluviométriques issues de radars météorologiques et d’un modèle hydrologique distribué pour la prévision de crue en temps réel

Abstract

Abstract After the destructive flood in 1998, the Chinese government planned to build national weather radar networks and to use radar data for real-time flood forecasting. Hence, coupling of weather radar rainfall data and a hydrological (Xinanjiang) model became an important issue. The present study reports on experience in such coupling at the Shiguanhe watershed. After having corrected the radar reflectivity and the attenuation data, the weather radar rainfall was estimated and then corrected in real time using a Kalman filter. In general, the precipitation estimated from weather radar is reasonably accurate in most of the catchment investigated, after corrections as above. Compared to the results simulated by raingauge data, the simulations based on the weather radar data are of similar accuracy. Present research results show that rainfall estimated from the weather radar, the radar data correction method, the method of coupling, and the Xinanjiang model lend themselves well to application in operational real-time flood forecasting.     

A non-parametric automatic blending methodology to estimate rainfall fields from rain gauge and radar data

Quantitative estimation of rainfall fields has been a crucial objective from early studies of the hydrological applications of weather radar. Previous studies have suggested that flow estimations are improved when radar and rain gauge data are combined to estimate input rainfall fields. This paper reports new research carried out in this field. Classical approaches for the selection and fitting of a theoretical correlogram (or semivariogram) model (needed to apply geostatistical estimators) are avoided in this study. Instead, a non-parametric technique based on FFT is used to obtain two-dimensional positive-definite correlograms directly from radar observations, dealing with both the natural anisotropy and the temporal variation of the spatial structure of the rainfall in the estimated fields. Because these correlation maps can be automatically obtained at each time step of a given rainfall event, this technique might easily be used in operational (real-time) applications. This paper describes the development of the non-parametric estimator exploiting the advantages of FFT for the automatic computation of correlograms and provides examples of its application on a case study using six rainfall events. This methodology is applied to three different alternatives to incorporate the radar information (as a secondary variable), and a comparison of performances is provided. In particular, their ability to reproduce in estimated rainfall fields (i) the rain gauge observations (in a cross-validation analysis) and (ii) the spatial patterns of radar fields are analyzed. Results seem to indicate that the methodology of kriging with external drift [KED], in combination with the technique of automatically computing 2-D spatial correlograms, provides merged rainfall fields with good agreement with rain gauges and with the most accurate approach to the spatial tendencies observed in the radar rainfall fields, when compared with other alternatives analyzed.     

Spatial rainfall estimation by linear and non-linear co-kriging of radar-rainfall and raingage data

The feasibility of linear and nonlinear geostatistical estimation techniques for optimal merging of rainfall data from raingage and radar observations is investigated in this study by use of controlled numerical experiments. Synthetic radar and raingage data are generated with their hypothetical error structures that explicitly account for sampling characteristics of the two sensors. Numerically simulated rainfall fields considered to be ground-truth fields on 4×4 km grids are used in the generation of radar and raingage observations. Ground-truth rainfall fields consist of generated rainfall fields with various climatic characteristics that preserve the space-time covariance function of rainfall events in extratropical cyclonic storms. Optimal mean areal precipitation estimates are obtained based on the minimum variance, unbiased property of kriging techniques under the second order homogeneity assumption of rainfall fields. The evaluation of estimated rainfall fields is done based on the refinement of spatial predictability over what would be provided from each sensor individually. Attention is mainly given to removal of measurement error and bias that are synthetically introduced to radar measurements. The influence of raingage network density on estimated rainfall fields is also examined.     

Probabilistic radar rainfall nowcasts using empirical and theoretical uncertainty models

Weather radar has a potential to provide accurate short‐term (0–3 h) forecasts of rainfall (i.e. radar nowcasts), which are of great importance in warnings and risk management for hydro‐meteorological events. However, radar nowcasts are affected by large uncertainties, which are not only linked to limitations in the forecast method but also because of errors in the radar rainfall measurement. The probabilistic quantitative precipitation nowcasting approach attempts to quantify these uncertainties by delivering the forecasts in a probabilistic form. This study implements two forms of probabilistic quantitative precipitation nowcasting for a hilly area in the south of Manchester, namely, the theoretically based scheme [ensemble rainfall forecasts (ERF)‐TN] and the empirically based scheme (ERF‐EM), and explores which one exhibits higher predictive skill. The ERF‐TN scheme generates ensemble forecasts of rainfall in which each ensemble member is determined by the stochastic realisation of a theoretical noise component. The so‐called ERF‐EM scheme proposed and applied for the first time in this study, aims to use an empirically based error model to measure and quantify the combined effect of all the error sources in the radar rainfall forecasts. The essence of the error model is formulated into an empirical relation between the radar rainfall forecasts and the corresponding ‘ground truth’ represented by the rainfall field from rain gauges measurements. The ensemble members generated by the two schemes have been compared with the rain gauge rainfall. The hit rate and the false alarm rate statistics have been computed and combined into relative operating characteristic curves. The comparison of the performance scores for the two schemes shows that the ERF‐EM achieves better performance than the ERF‐TN at 1‐h lead time. The predictive skills of both schemes are almost identical when the lead time increases to 2 h. In addition, the relation between uncertainty in the radar rainfall forecasts and lead time is also investigated by computing the dispersion of the generated ensemble members. Copyright © 2013 John Wiley & Sons, Ltd.     

An evaluation of numerical weather prediction based rainfall forecasts

ABSTRACT

Assessment of forecast precipitation is required before it can be used as input to hydrological models. Using radar observations in southeastern Australia, forecast rainfall from the Australian Community Climate Earth-System Simulator (ACCESS) was evaluated for 2010 and 2011. Radar rain intensities were first calibrated to gauge rainfall data from four research rainfall stations at hourly time steps. It is shown that the Australian ACCESS model (ACCESS-A) overestimated rainfall in low precipitation areas and underestimated elevated accumulations in high rainfall areas. The forecast errors were found to be dependent on the rainfall magnitude. Since the cumulative rainfall observations varied across the area and through the year, the relative error (RE) in the forecasts varied considerably with space and time, such that there was no consistent bias across the study area. Moreover, further analysis indicated that both location and magnitude errors were the main sources of forecast uncertainties on hourly accumulations, while magnitude was the dominant error on the daily time scale. Consequently, the precipitation output from ACCESS-A may not be useful for direct application in hydrological modelling, and pre-processing approaches such as bias correction or exceedance probability correction will likely be necessary for application of the numerical weather prediction (NWP) outputs.

EDITOR M.C. Acreman ASSOCIATE EDITOR A. Viglione     

Ability of a dual polarized X-band radar to estimate rainfall

The aim of this study is to assess rainfall estimates by a dual polarized X-band radar. This study was part of the European project FRAMEA (Flood forecasting using Radar in Alpine and Mediterranean Areas). Two radars were set up near the small town of Collobrières in South Eastern France. The first radar was a dual polarized X-band radar (Hydrix^®) associated with a ZPHI^® algorithm while the second one was an S-band radar (Météo France). We compared radar rainfall data with measurements obtained by two rain gauge networks (Météo France and Cemagref). During the experiments from February 2006 to June 2007, four significant rainfall events occurred. The accuracy of the rain rate obtained with both S-band and X-band radars decreased significantly beyond 60 km, in particular for the X-band radar. At closer ranges, such as 30–60 km from the radars, the X-band and the S-band radar retrievals showed similar performance with Nash criteria around 0.80 for the X-band radar and 0.75 for the S-band radar. Furthermore, the X-band radar did not require calibration on rainfall records, which tends to make it a useful method to assess rainfall in areas without a rain gauge network.     

Quantitative assessment of short‐term rainfall forecasts from radar nowcasts and MM5 forecasts

Short‐term Quantitative Precipitation Forecasts (QPFs) can be achieved from numerical weather prediction (NWP) models or radar nowcasting, that is the extrapolation of the precipitation at a future time from consecutive radar scans. Hybrid forecasts obtained by merging rainfall forecasts from radar nowcasting and NWP models are potentially more skilful than either radar nowcasts or NWP rainfall forecasts alone. This paper provides an assessment of deterministic and probabilistic high‐resolution QPFs achieved by implementing the Short‐term Ensemble Prediction System developed by the UK Met Office. Both radar nowcasts and hybrid forecasts have been performed. The results show that the performance of both deterministic nowcasts and deterministic hybrid forecasts decreases with increasing rainfall intensity and spatial resolution. The results also show that the blending with the NWP forecasts improves the performance of the forecasting system. Probabilistic hybrid forecasts have been obtained through the modelling of a stochastic noise component to produce a number of equally likely ensemble members, and the comparative assessment of deterministic and probabilistic hybrid forecasts shows that the probabilistic forecasting system is characterised by a higher discrimination accuracy than the deterministic one. Copyright © 2011 John Wiley & Sons, Ltd.     

Storage of water on vegetation under simulated rainfall of varying intensity

Little is understood about how storage of water on forest canopies varies during rainfall, even though storage changes intensity of throughfall and thus affects a variety of hydrological processes. In this study, laboratory rainfall simulation experiments using varying intensities yielded a better understanding of dynamics of rainfall storage on woody vegetation. Branches of eight species generally retained more water at higher rainfall intensities than at lower intensities, but incremental storage gains decreased as rainfall intensity increased. Leaf area was the best predictor of storage, especially for broadleaved species. Stored water ranged from 0.05 to 1.1 mm effective depth on leaves, depending on species and rainfall intensity. Storage was generally about 0.2 mm greater at rainfall intensity 420 mm h⁻¹ than at 20 mm h⁻¹. Needle-leaved branches generally retained more water per leaf area than did branches from broadleaved species, but branches that stored most at lower rainfall intensities tended to accumulate less additional storage at higher intensities. A simple nonlinear model was capable of predicting both magnitude (good model performance) and temporal scale (fair model performance) of storage responses to varying rainfall intensities. We hypothesize a conceptual mechanical model of canopy storage during rainfall that includes the concepts of static and dynamic storage to account for intensity-driven changes in storage. Scaling up observations to the canopy scale using LAI resulted in an estimate of canopy storage that generally agrees with estimates by traditional methods.     

Hydrological model sensitivity to parameter and radar rainfall estimation uncertainty

Radar estimates of rainfall are being increasingly applied to flood forecasting applications. Errors are inherent both in the process of estimating rainfall from radar and in the modelling of the rainfall–runoff transformation. The study aims at building a framework for the assessment of uncertainty that is consistent with the limitations of the model and data available and that allows a direct quantitative comparison between model predictions obtained by using radar and raingauge rainfall inputs. The study uses radar data from a mountainous region in northern Italy where complex topography amplifies radar errors due to radar beam occlusion and variability of precipitation with height. These errors, together with other error sources, are adjusted by applying a radar rainfall estimation algorithm. Radar rainfall estimates, adjusted and not, are used as an input to TOPMODEL for flood simulation over the Posina catchment (116 km²). Hydrological model parameter uncertainty is explicitly accounted for by use of the GLUE (Generalized Likelihood Uncertainty Estimation). Statistics are proposed to evaluate both the wideness of the uncertainty limits and the percentage of observations which fall within the uncertainty bounds. Results show the critical importance of proper adjustment of radar estimates and the use of radar estimates as close to ground as possible. Uncertainties affecting runoff predictions from adjusted radar data are close to those obtained by using a dense raingauge network, at least for the lowest radar observations available. Copyright © 2004 John Wiley & Sons, Ltd.     

A stochastic framework for downscaling processes of spatial averages based on the property of spectral multiscaling and its statistical diagnosis on spatio-temporal rainfall fields

Spectral multi-scaling postulates a power-law type of scaling of spectral distribution functions of stationary processes of spatial averages, over nested and geometrically similar sub-regions of the spatial parameter space of a given spatio-temporal random field. Presently a new framework is formulated for down-scaling processes of spatial averages, following naturally from the postulate of spectral multi-scaling, and key ingredients required for its implementation are described. Moreover, results from an extensive diagnostic study are presented, seeking statistical evidence supportive of spectral multi-scaling. Such evidence emerges from two sources of data. One is a 13 year long historical record of radar observations of rainfall in southeastern UK (Chenies radar), with high spatial (2 km) and temporal (5 min) resolution. The other is an ensemble of rain rate fields simulated by a spatio-temporal random pulse model fitted to the historical data. The results are consistent between historical and simulated rainfall data, indicating frequency-dependent scaling relationships interpreted as evidence of spectral multi-scaling across a range of spatial scales.     

Estimation of tail dependence coefficient in rainfall accumulation fields

Extreme rainfall events are of particular importance due to their severe impacts on the economy, the environment and the society. Characterization and quantification of extremes and their spatial dependence structure may lead to a better understanding of extreme events. An important concept in statistical modeling is the tail dependence coefficient (TDC) that describes the degree of association between concurrent rainfall extremes at different locations. Accurate knowledge of the spatial characteristics of the TDC can help improve on the existing models of the occurrence probability of extreme storms. In this study, efficient estimation of the TDC in rainfall is investigated using a dense network of rain gauges located in south Louisiana, USA. The inter-gauge distances in this network range from about 1 km to 9 km. Four different nonparametric TDC estimators are implemented on samples of the rain gauge data and their advantages and disadvantages are discussed. Three averaging time-scales are considered: 1 h, 2 h and 3 h. The results indicate that a significant tail dependency may exist that cannot be ignored for realistic modeling of multivariate rainfall fields. Presence of a strong dependence among extremes contradicts with the assumption of joint normality, commonly used in hydrologic applications.