Abstract This study focuses on the use of the bivariate normal distribution model to describe spatial distributions of tropical cyclone forecast errors. In this connection, it is found that forecast errors from the entire Atlantic tropical cyclone basin (Gulf of Mexico, Caribbean and North Atlantic) are multimodal and the fitting of these collective data to the usual unimodal bivariate normal distribution will be judged invalid by the usual statistical goodness-of-fit tests. While this is a recognized pitfall in classical statistics, it is often overlooked in meteorological application. The isolation of the clusters (components) and their statistical characteristics permit the issuance of forecast positions accompanied by more representative error ellipses. The study continues with a bivariate clustering analysis of a set of 979 tropical cyclone 24 h forecast errors for the Gulf of Mexico, Caribbean and North Atlantic. These errors were collected from the entire tropical cyclone basin without regard to se...