Summary: Knowledge of global patterns of biodiversity and regulating variables is indispensable to develop predictive models. The present study used predictive modelling approaches to investigate hypotheses that explain the variation in fish species richness between estuaries over a worldwide spatial extent. Ultimately, such models will allow assessment of future changes in ecosystem structure and function as a result of environmental changes. A comprehensive worldwide data base was compiled of the fish assemblage composition and environmental characteristics of estuaries. Generalized Linear Models were used to quantify how variation in species richness among estuaries is related to historical events, energy dynamics and ecosystem characteristics, while controlling for sampling effects. At the global extent, species richness differed among marine biogeographic realms and continents and increased with mean sea surface temperature, terrestrial net primary productivity and the stability of connectivity with a marine ecosystem (open vs. temporarily open estuaries). At a smaller extent (within a marine biogeographic realm or continent), other characteristics were also important in predicting variation in species richness, with species richness increasing with estuary area and continental shelf width. The results suggest that species richness in an estuary is defined by predictors that are spatially hierarchical. Over the largest spatial extents, species richness is influenced by the broader distributions and habitat use patterns of marine and freshwater species that can colonize estuaries, which are in turn governed by history contingency, energy dynamics and productivity variables. Species richness is also influenced by more regional and local parameters that can further affect the process of community colonization in an estuary including the connectivity of the estuary with the adjacent marine habitat, and, over smaller spatial extents, the size of these habitats. In summary, patterns of species richness in estuaries across large spatial extents seem to reflect from global to local processes acting on community colonization. The importance of considering spatial extent, sampling effects and of combining history and contemporary environmental characteristics when exploring biodiversity is highlighted.