Bridge collapse risk can be evaluated more rigorously if the hydrologic characteristics of bridge collapse sites are demystified, particularly for peak flows. In this study, forty-two bridge collapse sites were analyzed to find any trend in the peak flows. Flood frequency and other statistical analyses were used to derive peak flow distribution parameters, identify trends linked to flood magnitude and flood behavior (how extreme), quantify the return periods of peak flows, and compare different approaches of flood frequency in deriving the return periods. The results indicate that most of the bridge collapse sites exhibit heavy tail distribution and flood magnitudes that are well consistent when regressed over the drainage area. A comparison of different flood frequency analyses reveals that there is no single approach that is best generally for the dataset studied. These results indicate a commonality in flood behavior (outliers are expected, not random; heavy-tail property) for the collapse dataset studied and provides some basis for extending the findings obtained for the 42 collapsed bridges to other sites to assess the risk of future collapses.Water 2020, 12, 52 2 of 25 are used to directly define an extreme value distribution; and peaks over-threshold (POT), in which distributions are fit to both the frequency of floods above a threshold and their magnitudes. The U.S. Water Resources Council requires the use of log-Pearson type 3 (LP3) distributions for the annual maxima approach, although the generalized extreme value distribution (GEV) seems to be the most widely used model for extreme events [7]. The major benefit of the GEV model is its ability to fit highly skewed data [7]. Therefore, recent hydrologic research has focused on explaining the parameters of the GEV distribution of streamflow extremes [8][9][10][11][12][13][14].If annual maximum exceedances are assumed to be GEV-distributed, the POT exceedances are assumed to be generalized Pareto (GP)-distributed [15], following recommendations in the field of statistics [16,17]. Fitting the GP distribution to exceedances over a high threshold and also estimating the frequency of exceeding the threshold by fitting a Poisson distribution allows for the simultaneous fitting of parameters concerning both the frequency and intensity of extreme events. Compared to the annual maxima approach, therefore, the main advantage of POT modeling is that it allows for a more rational selection of events to be considered as "floods" and is not confined to only one event per year. The POT approach considers a wide range of events and provides the possibility of controlling the number of flood occurrences to be included in the analysis by appropriate selection of the threshold. However, the POT approach remains under-employed mainly because of the complexities associated with the choice of threshold and the selection of criteria for retaining flood peaks (Lang et al. (1999)). Nonetheless, threshold selection is tightly linked to the choice of the process distribution, to t...