By combining all Argo profiles for the period 2002 to present, a cumulative density function is constructed on a 3D grid of the global ocean. This function quantifies the statistics of isopycnals: time-averaged density, root- mean square (rms) of isopycnal displacement, and Eddy Available Potential Energy (EAPE). EAPE is the analogue of the eddy kinetic energy, but for the potential energy reservoir. Because it is essentially tied to the spatial struc- ture and magnitude of mesoscale activity, EAPE is an important quantity that should be useful to evaluate eddy resolving/permitting model turbu- lence and circulation. Among other striking features are the turbulent be- havior of Pacific and southern Atlantic Tsuchiya jets and subsurface EAPE maxima in some parts of the ocean, particularly in the Southern Ocean.