“…The distribution of extremes arising from mixing different populations of ordinary events is described here using a modified Metastatistical Extreme Value Distribution (MEVD) (Marani & Ignaccolo, 2015), which provides flexibility in incorporating the joint effect of different statistical populations and leverages the added value of incorporating physical mechanisms into statistical analysis (Klemeš, 1974). The MEVD relaxes some of the restrictive assumptions of the traditional Extreme Value Theory (EVT) and has been shown to outperform it in a wide range of applications, from daily and hourly rainfall, to remotely sensed precipitation, to hurricane intensities in the Atlantic Ocean, to peak flood flows (Zorzetto et al, 2016; Marra et al, 2018; Zorzetto & Marani, 2019; Zorzetto & Marani, 2020; Schellander et al, 2019; Hosseini, Scaioni, & Marani, 2020; Miniussi et al, 2020). Here we apply the mixed and original formulations of the MEVD to long series of daily rainfall in several American metropolitan areas, which have a high likelihood of being struck by a TC.…”