Abstract: Global maps of the mesoscale eddy available potential energy (EAPE) field at a depth of 500 m are created using potential density anomalies in a high‐resolution 1/12.5° global ocean model. Maps made from both a free‐running simulation and a data‐assimilative reanalysis of the HYbrid Coordinate Ocean Model (HYCOM) are compared with maps made by other researchers from density anomalies in Argo profiles. The HYCOM and Argo maps display similar features, especially in the dominance of western boundary currents. The reanalysis maps match the Argo maps more closely, demonstrating the added value of data assimilation. Global averages of the simulation, reanalysis, and Argo EAPE all agree to within about 10%. The model and Argo EAPE fields are compared to EAPE computed from temperature anomalies in a data set of “moored historical observations” (MHO) in conjunction with buoyancy frequencies computed from a global climatology. The MHO data set allows for an estimate of the EAPE in high‐frequency motions that is aliased into the Argo EAPE values. At MHO locations, 15–32% of the EAPE in the Argo estimates is due to aliased motions having periods of 10 days or less. Spatial averages of EAPE in HYCOM, Argo, and MHO data agree to within 50% at MHO locations, with both model estimates lying within error bars observations. Analysis of the EAPE field in an idealized model, in conjunction with published theory, suggests that much of the scatter seen in comparisons of different EAPE estimates is to be expected given the chaotic, unpredictable nature of mesoscale eddies.