The simultaneous occurrence of extreme events, such as simultaneous storms and floods at different locations, has a serious impact on risk assessment and mitigation strategies. The joint occurrence of extreme events can be measured by the so-called upper tail dependence (UTD) coefficient λ _U. In this study, we reconsider the properties of the most popular λ _U estimators and show that their strong bias and uncertainty make most of the empirical results reported in the hydrological literature questionable. In order to overcome the limits of λ _U analysis, we test several alternative tools such as a pool of formal statistical tests devised for recognizing upper tail independence and graphical diagnostics based on binary correlation and binary entropy. The reliability of all the methods is preliminarily checked by Monte Carlo experiments. Statistical tests and graphical diagnostics are therefore applied to three different rainfall data sets that allow us to explore the properties of the spatial dependence structure of rainfall extremes over a wide range of spatio-temporal scales ranging from 30 min and 1 km to 30 days and ≈ 3000 km. Results highlight that (1) classical estimators provide non zero tail dependence even for cases where it should be zero; (2) formal tests and binary correlation highlight that the pairwise spatial dependence structure can be weaker than Gaussian, thus excluding UTD calculated in a pairwise manner; (3) the binary entropy computed on triples of locations shows that the pairwise UTD is not enough to explain the spatial dependence structure of extreme rainfall, whose complexity becomes evident only after resorting to higher order correlation measures. The results concerning the bias and uncertainty of λ _U estimators are fully general and suggest avoiding their use especially for the short time series usually available in hydrology.